ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
ECCE @ EIC Software
Deprecated List
Modules
Namespaces
Classes
Files
File List
acts
analysis
coresoftware
Doxygen_Assist
ecce-detectors
fun4all_eicdetectors
geant4
tree
geant4-10.6-release
config
environments
examples
source
analysis
digits_hits
error_propagation
event
externals
g3tog4
geometry
global
graphics_reps
intercoms
interfaces
materials
parameterisations
particles
persistency
physics_lists
processes
biasing
cuts
decay
electromagnetic
hadronic
cross_sections
management
models
abla
abrasion
binary_cascade
cascade
coherent_elastic
de_excitation
em_dissociation
fission
gamma_nuclear
im_r_matrix
inclxx
lend
lepto_nuclear
management
particle_hp
include
src
G4ENDFTapeRead.cc
G4ENDFYieldDataContainer.cc
G4FFGDebuggingMacros.cc
G4FFGVerboseMacros.cc
G4FissionFragmentGenerator.cc
G4FissionProductYieldDist.cc
G4FPYBiasedLightFragmentDist.cc
G4FPYNormalFragmentDist.cc
G4FPYSamplingOps.cc
G4InterpolationManager.cc
G4NRESP71M03.cc
G4ParticleHP2AInelasticFS.cc
G4ParticleHP2N2AInelasticFS.cc
G4ParticleHP2NAInelasticFS.cc
G4ParticleHP2NDInelasticFS.cc
G4ParticleHP2NInelasticFS.cc
G4ParticleHP2NPInelasticFS.cc
G4ParticleHP2PInelasticFS.cc
G4ParticleHP3AInelasticFS.cc
G4ParticleHP3NAInelasticFS.cc
G4ParticleHP3NInelasticFS.cc
G4ParticleHP3NPInelasticFS.cc
G4ParticleHP4NInelasticFS.cc
G4ParticleHPAInelasticFS.cc
G4ParticleHPAngular.cc
G4ParticleHPArbitaryTab.cc
G4ParticleHPBGGNucleonInelasticXS.cc
G4ParticleHPCapture.cc
G4ParticleHPCaptureData.cc
G4ParticleHPCaptureFS.cc
G4ParticleHPChannel.cc
G4ParticleHPChannelList.cc
G4ParticleHPContAngularPar.cc
G4ParticleHPContEnergyAngular.cc
G4ParticleHPD2AInelasticFS.cc
G4ParticleHPDAInelasticFS.cc
G4ParticleHPData.cc
G4ParticleHPDeExGammas.cc
G4ParticleHPDInelasticFS.cc
G4ParticleHPDiscreteTwoBody.cc
G4ParticleHPElastic.cc
G4ParticleHPElasticData.cc
G4ParticleHPElasticFS.cc
G4ParticleHPElementData.cc
G4ParticleHPEnAngCorrelation.cc
G4ParticleHPFastLegendre.cc
G4ParticleHPFastLegendre_14.cc
G4ParticleHPFastLegendre_18.cc
G4ParticleHPFastLegendre_21.cc
G4ParticleHPFastLegendre_24.cc
G4ParticleHPFastLegendre_26.cc
G4ParticleHPFastLegendre_28.cc
G4ParticleHPFastLegendre_30.cc
G4ParticleHPFCFissionFS.cc
G4ParticleHPFFFissionFS.cc
G4ParticleHPField.cc
G4ParticleHPFieldPoint.cc
G4ParticleHPFinalState.cc
G4ParticleHPFission.cc
G4ParticleHPFissionBaseFS.cc
G4ParticleHPFissionData.cc
G4ParticleHPFissionFS.cc
G4ParticleHPFSFissionFS.cc
G4ParticleHPGamma.cc
G4ParticleHPHe3InelasticFS.cc
G4ParticleHPInelastic.cc
G4ParticleHPInelasticBaseFS.cc
G4ParticleHPInelasticCompFS.cc
G4ParticleHPInelasticData.cc
G4ParticleHPInterpolator.cc
G4ParticleHPIsoData.cc
G4ParticleHPIsotropic.cc
G4ParticleHPJENDLHEData.cc
G4ParticleHPJENDLHEElasticData.cc
G4ParticleHPJENDLHEInelasticData.cc
G4ParticleHPKallbachMannSyst.cc
G4ParticleHPLabAngularEnergy.cc
G4ParticleHPLCFissionFS.cc
G4ParticleHPLegendreStore.cc
G4ParticleHPLevel.cc
G4ParticleHPList.cc
G4ParticleHPMadlandNixSpectrum.cc
G4ParticleHPManager.cc
G4ParticleHPMessenger.cc
G4ParticleHPN2AInelasticFS.cc
G4ParticleHPN2PInelasticFS.cc
G4ParticleHPN3AInelasticFS.cc
G4ParticleHPNAInelasticFS.cc
G4ParticleHPNames.cc
G4ParticleHPNBodyPhaseSpace.cc
G4ParticleHPND2AInelasticFS.cc
G4ParticleHPNDInelasticFS.cc
G4ParticleHPNHe3InelasticFS.cc
G4ParticleHPNInelasticFS.cc
G4ParticleHPNPAInelasticFS.cc
G4ParticleHPNPInelasticFS.cc
G4ParticleHPNT2AInelasticFS.cc
G4ParticleHPNTInelasticFS.cc
G4ParticleHPNXInelasticFS.cc
G4ParticleHPPAInelasticFS.cc
G4ParticleHPPartial.cc
G4ParticleHPPDInelasticFS.cc
G4ParticleHPPhotonDist.cc
G4ParticleHPPInelasticFS.cc
G4ParticleHPProduct.cc
G4ParticleHPPTInelasticFS.cc
G4ParticleHPReactionWhiteBoard.cc
G4ParticleHPSCFissionFS.cc
G4ParticleHPT2AInelasticFS.cc
G4ParticleHPTCFissionFS.cc
G4ParticleHPThermalScattering.cc
G4ParticleHPThermalScatteringData.cc
G4ParticleHPThermalScatteringNames.cc
G4ParticleHPThreadLocalManager.cc
G4ParticleHPTInelasticFS.cc
G4ParticleHPVector.cc
G4ParticleHPWattSpectrum.cc
G4ShiftedGaussian.cc
G4WendtFissionFragmentGenerator.cc
parton_string
pre_equilibrium
qmd
quasi_elastic
radioactive_decay
rpg
theo_high_energy
util
processes
stopping
util
management
optical
parameterisation
scoring
solidstate
transportation
readout
run
track
tracking
visualization
macros
online_distribution
tutorials
doxygen_mainpage.h
File Members
External Links
•
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
G4ParticleHPNBodyPhaseSpace.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4ParticleHPNBodyPhaseSpace.cc
1
//
2
// ********************************************************************
3
// * License and Disclaimer *
4
// * *
5
// * The Geant4 software is copyright of the Copyright Holders of *
6
// * the Geant4 Collaboration. It is provided under the terms and *
7
// * conditions of the Geant4 Software License, included in the file *
8
// * LICENSE and available at http://cern.ch/geant4/license . These *
9
// * include a list of copyright holders. *
10
// * *
11
// * Neither the authors of this software system, nor their employing *
12
// * institutes,nor the agencies providing financial support for this *
13
// * work make any representation or warranty, express or implied, *
14
// * regarding this software system or assume any liability for its *
15
// * use. Please see the license in the file LICENSE and URL above *
16
// * for the full disclaimer and the limitation of liability. *
17
// * *
18
// * This code implementation is the result of the scientific and *
19
// * technical work of the GEANT4 collaboration. *
20
// * By using, copying, modifying or distributing the software (or *
21
// * any work based on the software) you agree to acknowledge its *
22
// * use in resulting scientific publications, and indicate your *
23
// * acceptance of all terms of the Geant4 Software license. *
24
// ********************************************************************
25
//
26
//
27
//
28
// P. Arce, June-2014 Conversion neutron_hp to particle_hp
29
//
30
#include "
G4ParticleHPNBodyPhaseSpace.hh
"
31
#include "
G4PhysicalConstants.hh
"
32
#include "
Randomize.hh
"
33
#include "
G4ThreeVector.hh
"
34
#include "
G4Gamma.hh
"
35
#include "
G4Electron.hh
"
36
#include "
G4Positron.hh
"
37
#include "
G4Neutron.hh
"
38
#include "
G4Proton.hh
"
39
#include "
G4Deuteron.hh
"
40
#include "
G4Triton.hh
"
41
#include "
G4He3.hh
"
42
#include "
G4Alpha.hh
"
43
44
G4ReactionProduct
*
G4ParticleHPNBodyPhaseSpace::Sample
(
G4double
anEnergy,
G4double
massCode,
G4double
)
45
{
46
G4ReactionProduct
* result =
new
G4ReactionProduct
;
47
G4int
Z
=
static_cast<
G4int
>
(massCode/1000);
48
G4int
A
=
static_cast<
G4int
>
(massCode-1000*
Z
);
49
50
if
(massCode==0)
51
{
52
result->
SetDefinition
(
G4Gamma::Gamma
());
53
}
54
else
if
(A==0)
55
{
56
result->
SetDefinition
(
G4Electron::Electron
());
57
if
(Z==1) result->
SetDefinition
(
G4Positron::Positron
());
58
}
59
else
if
(A==1)
60
{
61
result->
SetDefinition
(
G4Neutron::Neutron
());
62
if
(Z==1) result->
SetDefinition
(
G4Proton::Proton
());
63
}
64
else
if
(A==2)
65
{
66
result->
SetDefinition
(
G4Deuteron::Deuteron
());
67
}
68
else
if
(A==3)
69
{
70
result->
SetDefinition
(
G4Triton::Triton
());
71
if
(Z==2) result->
SetDefinition
(
G4He3::He3
());
72
}
73
else
if
(A==4)
74
{
75
result->
SetDefinition
(
G4Alpha::Alpha
());
76
if
(Z!=2)
throw
G4HadronicException
(__FILE__, __LINE__,
"Unknown ion case 1"
);
77
}
78
else
79
{
80
throw
G4HadronicException
(__FILE__, __LINE__,
"G4ParticleHPNBodyPhaseSpace: Unknown ion case 2"
);
81
}
82
83
// Get the energy from phase-space distribution
84
// in CMS
85
// P = Cn*std::sqrt(E')*(Emax-E')**(3*n/2-4)
86
G4double
maxE
=
GetEmax
(anEnergy, result->
GetMass
());
87
if
(maxE<=0){
88
maxE=1.*
CLHEP::eV
;
89
}
90
G4double
energy
;
91
G4double
max
(0);
92
if
(
theTotalCount
<=3)
93
{
94
max = maxE/2.;
95
}
96
else
if
(
theTotalCount
==4)
97
{
98
max = maxE/5.;
99
}
100
else
if
(
theTotalCount
==5)
101
{
102
max = maxE/8.;
103
}
104
else
105
{
106
throw
G4HadronicException
(__FILE__, __LINE__,
"NeutronHP Phase-space distribution cannot cope with this number of particles"
);
107
}
108
G4double
testit;
109
G4double
rand0 =
Prob
(max, maxE,
theTotalCount
);
110
G4double
rand;
111
112
G4int
icounter=0;
113
G4int
icounter_max=1024;
114
do
115
{
116
icounter++;
117
if
( icounter > icounter_max ) {
118
G4cout
<<
"Loop-counter exceeded the threshold value at "
<< __LINE__ <<
"th line of "
<< __FILE__ <<
"."
<<
G4endl
;
119
break
;
120
}
121
rand = rand0*
G4UniformRand
();
122
energy = maxE*
G4UniformRand
();
123
testit =
Prob
(energy, maxE,
theTotalCount
);
124
}
125
while
(rand > testit);
// Loop checking, 11.05.2015, T. Koi
126
result->
SetKineticEnergy
(energy);
127
128
// now do random direction
129
G4double
cosTh = 2.*
G4UniformRand
()-1.;
130
G4double
phi
=
twopi
*
G4UniformRand
();
131
G4double
theta
= std::acos(cosTh);
132
G4double
sinth = std::sin(theta);
133
G4double
mtot = result->
GetTotalMomentum
();
134
G4ThreeVector
tempVector(mtot*sinth*std::cos(phi), mtot*sinth*std::sin(phi), mtot*std::cos(theta) );
135
result->
SetMomentum
(tempVector);
136
G4ReactionProduct
aCMS = *
GetTarget
()+*
GetProjectileRP
();
137
result->
Lorentz
(*result, -1.*aCMS);
138
return
result;
139
}
geant4
tree
geant4-10.6-release
source
processes
hadronic
models
particle_hp
src
G4ParticleHPNBodyPhaseSpace.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:48
using
1.8.2 with
ECCE GitHub integration