Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Member Functions | Static Protected Attributes
G4RPGNucleonInelastic Class Reference

#include <G4RPGNucleonInelastic.hh>

Inheritance diagram for G4RPGNucleonInelastic:
G4RPGInelastic G4HadronicInteraction G4RPGNeutronInelastic G4RPGProtonInelastic

Public Member Functions

 G4RPGNucleonInelastic (const G4String &modelName="RPGNucleonInelastic")
 
 ~G4RPGNucleonInelastic ()
 
- Public Member Functions inherited from G4RPGInelastic
 G4RPGInelastic (const G4String &modelName="RPGInelastic")
 
virtual ~G4RPGInelastic ()
 
- Public Member Functions inherited from G4HadronicInteraction
 G4HadronicInteraction (const G4String &modelName="HadronicModel")
 
virtual ~G4HadronicInteraction ()
 
virtual G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)=0
 
virtual G4double SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
 
virtual G4bool IsApplicable (const G4HadProjectile &, G4Nucleus &)
 
G4double GetMinEnergy () const
 
G4double GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMinEnergy (G4double anEnergy)
 
void SetMinEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
 
G4double GetMaxEnergy () const
 
G4double GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
 
void SetMaxEnergy (const G4double anEnergy)
 
void SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
 
void SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
 
const G4HadronicInteraction * GetMyPointer () const
 
virtual G4int GetVerboseLevel () const
 
virtual void SetVerboseLevel (G4int value)
 
const G4String & GetModelName () const
 
void DeActivateFor (const G4Material *aMaterial)
 
void ActivateFor (const G4Material *aMaterial)
 
void DeActivateFor (const G4Element *anElement)
 
void ActivateFor (const G4Element *anElement)
 
G4bool IsBlocked (const G4Material *aMaterial) const
 
G4bool IsBlocked (const G4Element *anElement) const
 
void SetRecoilEnergyThreshold (G4double val)
 
G4double GetRecoilEnergyThreshold () const
 
G4bool operator== (const G4HadronicInteraction &right) const
 
G4bool operator!= (const G4HadronicInteraction &right) const
 
virtual const std::pair
< G4double, G4double > 		GetFatalEnergyCheckLevels () const
 
virtual std::pair< G4double,
G4double > 		GetEnergyMomentumCheckLevels () const
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
virtual void ModelDescription (std::ostream &outFile) const
 

Protected Member Functions

G4int GetMultiplicityT1 (G4double KE) const
 
G4int GetMultiplicityT0 (G4double KE) const
 
std::vector< G4int > GetFSPartTypesForT1 (G4int mult, G4double KE, G4int tindex) const
 
std::vector< G4int > GetFSPartTypesForT0 (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForPP (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForNN (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForPN (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForNP (G4int mult, G4double KE) const
 
- Protected Member Functions inherited from G4RPGInelastic
G4double Pmltpc (G4int np, G4int nm, G4int nz, G4int n, G4double b, G4double c)
 
G4int Factorial (G4int n)
 
G4bool MarkLeadingStrangeParticle (const G4ReactionProduct &currentParticle, const G4ReactionProduct &targetParticle, G4ReactionProduct &leadParticle)
 
void SetUpPions (const G4int np, const G4int nm, const G4int nz, G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen)
 
void GetNormalizationConstant (const G4double availableEnergy, G4double &n, G4double &anpn)
 
void CalculateMomenta (G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, const G4HadProjectile *originalIncident, const G4DynamicParticle *originalTarget, G4ReactionProduct &modifiedOriginal, G4Nucleus &targetNucleus, G4ReactionProduct &currentParticle, G4ReactionProduct &targetParticle, G4bool &incidentHasChanged, G4bool &targetHasChanged, G4bool quasiElastic)
 
void SetUpChange (G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, G4ReactionProduct &currentParticle, G4ReactionProduct &targetParticle, G4bool &incidentHasChanged)
 
std::pair< G4int, G4double > interpolateEnergy (G4double ke) const
 
G4int sampleFlat (std::vector< G4double > sigma) const
 
void CheckQnums (G4FastVector< G4ReactionProduct, 256 > &vec, G4int &vecLen, G4ReactionProduct &currentParticle, G4ReactionProduct &targetParticle, G4double Q, G4double B, G4double S)
 
- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 

Static Protected Attributes

static const G4int pPindex [8][2]
 
static const G4int pNindex [8][2]
 
static const G4int T1_2bfs [2][1][2]
 
static const G4int T1_3bfs [2][6][3]
 
static const G4int T1_4bfs [2][18][4]
 
static const G4int T1_5bfs [2][32][5]
 
static const G4int T1_6bfs [2][7][6]
 
static const G4int T1_7bfs [2][8][7]
 
static const G4int T1_8bfs [2][10][8]
 
static const G4int T1_9bfs [2][11][9]
 
static const G4int T0_2bfs [1][2]
 
static const G4int T0_3bfs [9][3]
 
static const G4int T0_4bfs [22][4]
 
static const G4int T0_5bfs [38][5]
 
static const G4int T0_6bfs [7][6]
 
static const G4int T0_7bfs [9][7]
 
static const G4int T0_8bfs [10][8]
 
static const G4int T0_9bfs [12][9]
 
static G4ThreadLocal G4double pPtot [30]
 
static G4ThreadLocal G4double pNtot [30]
 
static G4ThreadLocal G4double t1_dSigma_dMult [8][30]
 
static G4ThreadLocal G4double t0_dSigma_dMult [8][30]
 
static const G4float pPCrossSections [93][30]
 
static const G4float pNCrossSections [108][30]
 

Additional Inherited Members

- Protected Types inherited from G4RPGInelastic
enum  {
  pi0, pip, pim, kp,
  km, k0, k0b, pro,
  neu, lam, sp, s0,
  sm, xi0, xim, om,
  ap, an
}
 
- Protected Attributes inherited from G4RPGInelastic
G4RPGFragmentation fragmentation
 
G4RPGTwoCluster twoCluster
 
G4RPGPionSuppression pionSuppression
 
G4RPGStrangeProduction strangeProduction
 
G4RPGTwoBody twoBody
 
G4ParticleDefinition * particleDef [18]
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 42 of file G4RPGNucleonInelastic.hh.

Constructor & Destructor Documentation

G4RPGNucleonInelastic::G4RPGNucleonInelastic ( const G4String &  modelName = "RPGNucleonInelastic")

Definition at line 33 of file G4RPGNucleonInelastic.cc.

References python.hepunit::GeV, pNCrossSections, pNindex, pNtot, pPCrossSections, pPindex, pPtot, G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), t0_dSigma_dMult, and t1_dSigma_dMult.

34  :G4RPGInelastic(modelName)
35 {
36  SetMinEnergy( 0.0 );
37  SetMaxEnergy( 30.*GeV );
38 
39  // Initialize t1_dSigma_dMult, t0_dSigma_dMult,
40  // nucleon-nucleon inelastic cross sections for a given multiplicity
41  // for |T_z| = 1 and 0, respectively
42 
43  G4int i, j, k;
44  G4int start, stop;
45 
46  for (j = 0; j < 8; j++) {
47  start = pPindex[j][0];
48  stop = pPindex[j][1] + 1;
49  for (k = 0; k < 30; k++) {
50  t1_dSigma_dMult[j][k] = 0.0;
51  for (i = start; i < stop; i++) t1_dSigma_dMult[j][k] += pPCrossSections[i][k];
52  }
53 
54  start = pNindex[j][0];
55  stop = pNindex[j][1] + 1;
56  for (k = 0; k < 30; k++) {
57  t0_dSigma_dMult[j][k] = 0.0;
58  for (i = start; i < stop; i++) t0_dSigma_dMult[j][k] += pNCrossSections[i][k];
59  }
60  }
61 
62  // Initialize total cross section array
63 
64  for (k = 0; k < 30; k++) {
65  pPtot[k] = 0.0;
66  pNtot[k] = 0.0;
67  for (j = 0; j < 8; j++) {
68  pPtot[k] += t1_dSigma_dMult[j][k];
69  pNtot[k] += t0_dSigma_dMult[j][k];
70  }
71  }
72 
73  // printCrossSections();
74 }
G4int
int G4int
Definition: G4Types.hh:78
G4RPGNucleonInelastic::pPCrossSections
static const G4float pPCrossSections[93][30]
Definition: G4RPGNucleonInelastic.hh:101
G4RPGNucleonInelastic::pNtot
static G4ThreadLocal G4double pNtot[30]
Definition: G4RPGNucleonInelastic.hh:97
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:90
G4RPGNucleonInelastic::pPtot
static G4ThreadLocal G4double pPtot[30]
Definition: G4RPGNucleonInelastic.hh:96
G4RPGInelastic::G4RPGInelastic
G4RPGInelastic(const G4String &modelName="RPGInelastic")
Definition: G4RPGInelastic.cc:38
G4RPGNucleonInelastic::pNindex
static const G4int pNindex[8][2]
Definition: G4RPGNucleonInelastic.hh:76
G4RPGNucleonInelastic::pPindex
static const G4int pPindex[8][2]
Definition: G4RPGNucleonInelastic.hh:75
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:103
G4RPGNucleonInelastic::pNCrossSections
static const G4float pNCrossSections[108][30]
Definition: G4RPGNucleonInelastic.hh:102
G4RPGNucleonInelastic::t1_dSigma_dMult
static G4ThreadLocal G4double t1_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:98
G4RPGNucleonInelastic::t0_dSigma_dMult
static G4ThreadLocal G4double t0_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:99
G4RPGNucleonInelastic::~G4RPGNucleonInelastic ( )
inline

Definition at line 48 of file G4RPGNucleonInelastic.hh.

48 { }

Member Function Documentation

std::vector<G4int> G4RPGNucleonInelastic::GetFSPartTypesForNN ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 65 of file G4RPGNucleonInelastic.hh.

References GetFSPartTypesForT1().

66  {return GetFSPartTypesForT1(mult, KE, 1); }
G4RPGNucleonInelastic::GetFSPartTypesForT1
std::vector< G4int > GetFSPartTypesForT1(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGNucleonInelastic.cc:186
std::vector<G4int> G4RPGNucleonInelastic::GetFSPartTypesForNP ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 71 of file G4RPGNucleonInelastic.hh.

References GetFSPartTypesForT0().

72  {return GetFSPartTypesForT0(mult, KE); }
G4RPGNucleonInelastic::GetFSPartTypesForT0
std::vector< G4int > GetFSPartTypesForT0(G4int mult, G4double KE) const
Definition: G4RPGNucleonInelastic.cc:138
std::vector<G4int> G4RPGNucleonInelastic::GetFSPartTypesForPN ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 68 of file G4RPGNucleonInelastic.hh.

References GetFSPartTypesForT0().

69  {return GetFSPartTypesForT0(mult, KE); }
G4RPGNucleonInelastic::GetFSPartTypesForT0
std::vector< G4int > GetFSPartTypesForT0(G4int mult, G4double KE) const
Definition: G4RPGNucleonInelastic.cc:138
std::vector<G4int> G4RPGNucleonInelastic::GetFSPartTypesForPP ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 62 of file G4RPGNucleonInelastic.hh.

References GetFSPartTypesForT1().

63  {return GetFSPartTypesForT1(mult, KE, 0); }
G4RPGNucleonInelastic::GetFSPartTypesForT1
std::vector< G4int > GetFSPartTypesForT1(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGNucleonInelastic.cc:186
std::vector< G4int > G4RPGNucleonInelastic::GetFSPartTypesForT0 ( G4int  mult,
G4double  KE 
) const
protected

Definition at line 138 of file G4RPGNucleonInelastic.cc.

References G4cout, G4endl, G4RPGInelastic::interpolateEnergy(), pNCrossSections, pNindex, G4RPGInelastic::sampleFlat(), T0_2bfs, T0_3bfs, T0_4bfs, T0_5bfs, T0_6bfs, T0_7bfs, T0_8bfs, and T0_9bfs.

Referenced by GetFSPartTypesForNP(), and GetFSPartTypesForPN().

139 {
140  G4int i;
141  G4double sigint(0.);
142  std::vector<G4double> sigma;
143 
144  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
145  G4int k = epair.first;
146  G4double fraction = epair.second;
147 
148  G4int start = pNindex[mult-2][0];
149  G4int stop = pNindex[mult-2][1];
150 
151  for(i = start; i < stop; i++) {
152  sigint = pNCrossSections[i][k]
153  + fraction*(pNCrossSections[i][k+1] - pNCrossSections[i][k]);
154  sigma.push_back(sigint);
155  }
156 
157  G4int channel = sampleFlat(sigma);
158 
159  std::vector<G4int> kinds;
160 
161  if (mult == 2) {
162  for(i = 0; i < mult; i++) kinds.push_back(T0_2bfs[channel][i]);
163  } else if (mult == 3) {
164  for(i = 0; i < mult; i++) kinds.push_back(T0_3bfs[channel][i]);
165  } else if (mult == 4) {
166  for(i = 0; i < mult; i++) kinds.push_back(T0_4bfs[channel][i]);
167  } else if (mult == 5) {
168  for(i = 0; i < mult; i++) kinds.push_back(T0_5bfs[channel][i]);
169  } else if (mult == 6) {
170  for(i = 0; i < mult; i++) kinds.push_back(T0_6bfs[channel][i]);
171  } else if (mult == 7) {
172  for(i = 0; i < mult; i++) kinds.push_back(T0_7bfs[channel][i]);
173  } else if (mult == 8) {
174  for(i = 0; i < mult; i++) kinds.push_back(T0_8bfs[channel][i]);
175  } else if (mult == 9) {
176  for(i = 0; i < mult; i++) kinds.push_back(T0_9bfs[channel][i]);
177  } else {
178  G4cout << " Illegal multiplicity " << G4endl;
179  }
180 
181  return kinds;
182 }
G4RPGNucleonInelastic::T0_8bfs
static const G4int T0_8bfs[10][8]
Definition: G4RPGNucleonInelastic.hh:93
G4RPGNucleonInelastic::T0_2bfs
static const G4int T0_2bfs[1][2]
Definition: G4RPGNucleonInelastic.hh:87
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4RPGNucleonInelastic::T0_7bfs
static const G4int T0_7bfs[9][7]
Definition: G4RPGNucleonInelastic.hh:92
G4int
int G4int
Definition: G4Types.hh:78
G4cout
G4GLOB_DLL std::ostream G4cout
G4RPGNucleonInelastic::T0_9bfs
static const G4int T0_9bfs[12][9]
Definition: G4RPGNucleonInelastic.hh:94
G4RPGNucleonInelastic::T0_4bfs
static const G4int T0_4bfs[22][4]
Definition: G4RPGNucleonInelastic.hh:89
G4RPGNucleonInelastic::T0_3bfs
static const G4int T0_3bfs[9][3]
Definition: G4RPGNucleonInelastic.hh:88
G4RPGNucleonInelastic::pNindex
static const G4int pNindex[8][2]
Definition: G4RPGNucleonInelastic.hh:76
G4endl
#define G4endl
Definition: G4ios.hh:61
G4RPGNucleonInelastic::pNCrossSections
static const G4float pNCrossSections[108][30]
Definition: G4RPGNucleonInelastic.hh:102
G4double
double G4double
Definition: G4Types.hh:76
G4RPGNucleonInelastic::T0_6bfs
static const G4int T0_6bfs[7][6]
Definition: G4RPGNucleonInelastic.hh:91
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGNucleonInelastic::T0_5bfs
static const G4int T0_5bfs[38][5]
Definition: G4RPGNucleonInelastic.hh:90
std::vector< G4int > G4RPGNucleonInelastic::GetFSPartTypesForT1 ( G4int  mult,
G4double  KE,
G4int  tindex 
) const
protected

Definition at line 186 of file G4RPGNucleonInelastic.cc.

References G4cout, G4endl, G4RPGInelastic::interpolateEnergy(), pPCrossSections, pPindex, G4RPGInelastic::sampleFlat(), T1_2bfs, T1_3bfs, T1_4bfs, T1_5bfs, T1_6bfs, T1_7bfs, T1_8bfs, and T1_9bfs.

Referenced by GetFSPartTypesForNN(), and GetFSPartTypesForPP().

188 {
189  G4int i;
190  G4double sigint(0.);
191  std::vector<G4double> sigma;
192 
193  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
194  G4int k = epair.first;
195  G4double fraction = epair.second;
196 
197  G4int start = pPindex[mult-2][0];
198  G4int stop = pPindex[mult-2][1];
199 
200  for(i = start; i < stop; i++) {
201  sigint = pPCrossSections[i][k]
202  + fraction*(pPCrossSections[i][k+1] - pPCrossSections[i][k]);
203  sigma.push_back(sigint);
204  }
205 
206  G4int channel = sampleFlat(sigma);
207 
208  std::vector<G4int> kinds;
209 
210  if (mult == 2) {
211  for(i = 0; i < mult; i++) kinds.push_back(T1_2bfs[tzindex][channel][i]);
212  } else if (mult == 3) {
213  for(i = 0; i < mult; i++) kinds.push_back(T1_3bfs[tzindex][channel][i]);
214  } else if (mult == 4) {
215  for(i = 0; i < mult; i++) kinds.push_back(T1_4bfs[tzindex][channel][i]);
216  } else if (mult == 5) {
217  for(i = 0; i < mult; i++) kinds.push_back(T1_5bfs[tzindex][channel][i]);
218  } else if (mult == 6) {
219  for(i = 0; i < mult; i++) kinds.push_back(T1_6bfs[tzindex][channel][i]);
220  } else if (mult == 7) {
221  for(i = 0; i < mult; i++) kinds.push_back(T1_7bfs[tzindex][channel][i]);
222  } else if (mult == 8) {
223  for(i = 0; i < mult; i++) kinds.push_back(T1_8bfs[tzindex][channel][i]);
224  } else if (mult == 9) {
225  for(i = 0; i < mult; i++) kinds.push_back(T1_9bfs[tzindex][channel][i]);
226  } else {
227  G4cout << " Illegal multiplicity " << G4endl;
228  }
229 
230  return kinds;
231 }
G4RPGNucleonInelastic::T1_8bfs
static const G4int T1_8bfs[2][10][8]
Definition: G4RPGNucleonInelastic.hh:84
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4int
int G4int
Definition: G4Types.hh:78
G4RPGNucleonInelastic::pPCrossSections
static const G4float pPCrossSections[93][30]
Definition: G4RPGNucleonInelastic.hh:101
G4cout
G4GLOB_DLL std::ostream G4cout
G4RPGNucleonInelastic::T1_4bfs
static const G4int T1_4bfs[2][18][4]
Definition: G4RPGNucleonInelastic.hh:80
G4RPGNucleonInelastic::T1_5bfs
static const G4int T1_5bfs[2][32][5]
Definition: G4RPGNucleonInelastic.hh:81
G4RPGNucleonInelastic::T1_2bfs
static const G4int T1_2bfs[2][1][2]
Definition: G4RPGNucleonInelastic.hh:78
G4RPGNucleonInelastic::T1_6bfs
static const G4int T1_6bfs[2][7][6]
Definition: G4RPGNucleonInelastic.hh:82
G4RPGNucleonInelastic::pPindex
static const G4int pPindex[8][2]
Definition: G4RPGNucleonInelastic.hh:75
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76
G4RPGNucleonInelastic::T1_7bfs
static const G4int T1_7bfs[2][8][7]
Definition: G4RPGNucleonInelastic.hh:83
G4RPGNucleonInelastic::T1_3bfs
static const G4int T1_3bfs[2][6][3]
Definition: G4RPGNucleonInelastic.hh:79
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGNucleonInelastic::T1_9bfs
static const G4int T1_9bfs[2][11][9]
Definition: G4RPGNucleonInelastic.hh:85
G4int G4RPGNucleonInelastic::GetMultiplicityT0 ( G4double  KE) const
protected

Definition at line 99 of file G4RPGNucleonInelastic.cc.

References G4RPGInelastic::interpolateEnergy(), G4RPGInelastic::sampleFlat(), and t0_dSigma_dMult.

100 {
101  G4double multint(0.);
102  std::vector<G4double> sigma;
103 
104  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
105  G4int k = epair.first;
106  G4double fraction = epair.second;
107 
108  for(G4int j = 0; j < 8; j++) {
109  multint = t0_dSigma_dMult[j][k]
110  + fraction*(t0_dSigma_dMult[j][k+1] - t0_dSigma_dMult[j][k]);
111  sigma.push_back(multint);
112  }
113 
114  return sampleFlat(sigma) + 2;
115 }
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4int
int G4int
Definition: G4Types.hh:78
G4double
double G4double
Definition: G4Types.hh:76
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGNucleonInelastic::t0_dSigma_dMult
static G4ThreadLocal G4double t0_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:99
G4int G4RPGNucleonInelastic::GetMultiplicityT1 ( G4double  KE) const
protected

Definition at line 118 of file G4RPGNucleonInelastic.cc.

References G4RPGInelastic::interpolateEnergy(), G4RPGInelastic::sampleFlat(), and t1_dSigma_dMult.

119 {
120  G4double multint(0.);
121  std::vector<G4double> sigma;
122 
123  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
124  G4int k = epair.first;
125  G4double fraction = epair.second;
126 
127  for(G4int j = 0; j < 8; j++) {
128  multint = t1_dSigma_dMult[j][k]
129  + fraction*(t1_dSigma_dMult[j][k+1] - t1_dSigma_dMult[j][k]);
130  sigma.push_back(multint);
131  }
132 
133  return sampleFlat(sigma) + 2;
134 }
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4int
int G4int
Definition: G4Types.hh:78
G4double
double G4double
Definition: G4Types.hh:76
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGNucleonInelastic::t1_dSigma_dMult
static G4ThreadLocal G4double t1_dSigma_dMult[8][30]
Definition: G4RPGNucleonInelastic.hh:98

Field Documentation

const G4float G4RPGNucleonInelastic::pNCrossSections
staticprotected

Definition at line 102 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic(), and GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::pNindex
staticprotected
Initial value:
=
{{0, 0}, {1,9}, {10,31}, {32,69}, {70,76}, {77,85}, {86,95}, {96,107}}

Definition at line 76 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic(), and GetFSPartTypesForT0().

G4ThreadLocal G4double G4RPGNucleonInelastic::pNtot
staticprotected

Definition at line 97 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic().

const G4float G4RPGNucleonInelastic::pPCrossSections
staticprotected

Definition at line 101 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic(), and GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::pPindex
staticprotected
Initial value:
=
{{0, 0}, {1, 6}, {7,24}, {25,56}, {57,63}, {64,71}, {72,81}, {82,92}}

Definition at line 75 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic(), and GetFSPartTypesForT1().

G4ThreadLocal G4double G4RPGNucleonInelastic::pPtot
staticprotected

Definition at line 96 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic().

const G4int G4RPGNucleonInelastic::T0_2bfs
staticprotected
Initial value:
=
{{pro,neu}}

Definition at line 87 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_3bfs
staticprotected
Initial value:
=
{{pro,pro,pim},{pro,neu,pi0},{neu,neu,pip},{pro,lam,k0},
{pro,s0,k0}, {pro,sm,kp}, {neu,lam,kp}, {neu,s0,kp},
{neu,sp,k0}}

Definition at line 88 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_4bfs
staticprotected
Initial value:
=
{{pro,neu,pip,pim},{pro,pro,pim,pi0},{pro,neu,pi0,pi0},
{neu,neu,pip,pi0},{pro,lam,kp,pim}, {pro,s0,kp,pim},
{pro,lam,k0,pi0}, {pro,s0,k0,pi0}, {pro,sp,k0,pim},
{pro,sm,kp,pi0}, {pro,sm,k0,pip}, {neu,lam,kp,pi0},
{neu,lam,k0,pip}, {neu,sp,kp,pim}, {neu,sp,k0,pi0},
{neu,s0,kp,pi0}, {neu,s0,k0,pip}, {neu,sm,kp,pip},
{pro,neu,kp,km}, {pro,neu,k0,k0b}, {pro,pro,k0,km},
{neu,neu,kp,k0b}}

Definition at line 89 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_5bfs
staticprotected
Initial value:
=
{{pro,neu,pip,pim,pi0},{pro,neu,pi0,pi0,pi0},{pro,pro,pip,pim,pim},
{pro,pro,pim,pi0,pi0},{neu,neu,pip,pip,pim},{neu,neu,pip,pi0,pi0},
{pro,lam,kp,pim,pi0}, {pro,lam,k0,pip,pim}, {pro,lam,k0,pi0,pi0},
{pro,s0,k0,pip,pim}, {pro,s0,k0,pi0,pi0}, {pro,s0,kp,pim,pi0},
{pro,sp,kp,pim,pim}, {pro,sp,k0,pim,pi0}, {pro,sm,k0,pip,pi0},
{pro,sm,kp,pip,pim}, {pro,sm,kp,pi0,pi0}, {neu,lam,kp,pip,pim},
{neu,lam,kp,pi0,pi0}, {neu,lam,k0,pip,pi0}, {neu,s0,kp,pip,pim},
{neu,s0,kp,pi0,pi0}, {neu,s0,k0,pip,pi0}, {neu,sp,k0,pip,pim},
{neu,sp,k0,pi0,pi0}, {neu,sp,kp,pim,pi0}, {neu,sm,kp,pip,pi0},
{neu,sm,k0,pip,pip}, {pro,neu,kp,km,pi0}, {pro,neu,k0,k0b,pi0},
{pro,neu,k0,km,pip}, {pro,neu,kp,k0b,pim}, {pro,pro,k0,k0b,pim},
{pro,pro,kp,km,pim}, {pro,pro,k0,km,pi0}, {neu,neu,kp,km,pip},
{neu,neu,k0,k0b,pip}, {neu,neu,kp,k0b,pi0}}

Definition at line 90 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_6bfs
staticprotected
Initial value:
=
{{pro,neu,pip,pip,pim,pim},{pro,neu,pip,pim,pi0,pi0},
{pro,neu,pi0,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pi0},
{pro,pro,pim,pi0,pi0,pi0},{neu,neu,pip,pip,pim,pi0},
{neu,neu,pip,pi0,pi0,pi0}}

Definition at line 91 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_7bfs
staticprotected
Initial value:
=
{{pro,neu,pip,pip,pim,pim,pi0},{pro,neu,pip,pim,pi0,pi0,pi0},
{pro,neu,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim},
{pro,pro,pip,pim,pim,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0},
{neu,neu,pip,pip,pip,pim,pim},{neu,neu,pip,pip,pim,pi0,pi0},
{neu,neu,pip,pi0,pi0,pi0,pi0}}

Definition at line 92 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_8bfs
staticprotected
Initial value:
=
{{pro,neu,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pi0,pi0},
{pro,neu,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pro,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pi0,pi0,pi0},
{pro,pro,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim,pim,pi0},
{neu,neu,pip,pip,pim,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0}}

Definition at line 93 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

const G4int G4RPGNucleonInelastic::T0_9bfs
staticprotected
Initial value:
=
{{pro,neu,pip,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pip,pim,pim,pi0,pi0,pi0},
{pro,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pro,pip,pip,pip,pim,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pim,pi0,pi0},
{pro,pro,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,neu,pip,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pip,pim,pim,pi0,pi0},
{neu,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0}}

Definition at line 94 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT0().

G4ThreadLocal G4double G4RPGNucleonInelastic::t0_dSigma_dMult
staticprotected

Definition at line 99 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic(), and GetMultiplicityT0().

const G4int G4RPGNucleonInelastic::T1_2bfs
staticprotected
Initial value:
=
{{{pro,pro}},
{{neu,neu}}}

Definition at line 78 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_3bfs
staticprotected
Initial value:
=
{{{pro,pro,pi0}, {pro,neu,pip}, {pro,lam,kp},
{pro,s0,kp}, {pro,sp,k0}, {neu,sp,kp}},
{{neu,neu,pi0}, {pro,neu,pim}, {neu,lam,k0},
{neu,s0,k0}, {neu,sm,kp}, {pro,sm,k0}}}

Definition at line 79 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_4bfs
staticprotected
Initial value:
=
{{{pro,pro,pip,pim},{pro,neu,pip,pi0},{pro,pro,pi0,pi0},
{neu,neu,pip,pip},{pro,lam,kp,pi0}, {pro,lam,k0,pip},
{neu,lam,kp,pip}, {neu,s0,kp,pip}, {pro,s0,kp,pi0},
{pro,s0,k0,pip}, {pro,sm,kp,pip}, {pro,sp,k0,pi0},
{neu,sp,k0,pip}, {pro,sp,kp,pim}, {neu,sp,kp,pi0},
{pro,pro,k0,k0b}, {pro,pro,kp,km}, {pro,neu,kp,k0b}},
{{neu,neu,pip,pim},{pro,neu,pim,pi0},{neu,neu,pi0,pi0},
{pro,pro,pim,pim},{neu,lam,k0,pi0}, {neu,lam,kp,pim},
{pro,lam,k0,pim}, {pro,s0,k0,pim}, {neu,s0,k0,pi0},
{neu,s0,kp,pim}, {neu,sp,k0,pim}, {neu,sm,kp,pi0},
{pro,sm,kp,pim}, {neu,sm,k0,pip}, {pro,sm,k0,pi0},
{neu,neu,kp,km}, {neu,neu,k0,k0b}, {pro,neu,k0,km}}}

Definition at line 80 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_5bfs
staticprotected
Initial value:
=
{{{pro,pro,pip,pim,pi0},{pro,pro,pi0,pi0,pi0},{pro,neu,pip,pip,pim},
{pro,neu,pip,pi0,pi0},{neu,neu,pip,pip,pi0},{pro,lam,kp,pip,pim},
{pro,lam,kp,pi0,pi0}, {pro,lam,k0,pip,pi0}, {pro,s0,kp,pip,pim},
{pro,s0,kp,pi0,pi0}, {pro,s0,k0,pip,pi0}, {pro,sp,k0,pip,pim},
{pro,sp,k0,pi0,pi0}, {pro,sp,kp,pim,pi0}, {pro,sm,kp,pip,pi0},
{pro,sm,k0,pip,pip}, {neu,lam,kp,pip,pi0}, {neu,lam,k0,pip,pip},
{neu,s0,kp,pip,pi0}, {neu,s0,k0,pip,pip}, {neu,sp,k0,pip,pi0},
{neu,sp,kp,pip,pim}, {neu,sp,kp,pi0,pi0}, {neu,sm,kp,pip,pip},
{pro,pro,pip,k0,km}, {pro,pro,pim,kp,k0b}, {pro,pro,pi0,k0,k0b},
{pro,pro,pi0,kp,km}, {pro,neu,pip,k0,k0b}, {pro,neu,pip,kp,km},
{pro,neu,pi0,kp,k0b}, {neu,neu,pip,kp,k0b}},
{{neu,neu,pip,pim,pi0},{neu,neu,pi0,pi0,pi0},{pro,neu,pip,pim,pim},
{pro,neu,pim,pi0,pi0},{pro,pro,pim,pim,pi0},{neu,lam,k0,pip,pim},
{neu,lam,k0,pi0,pi0}, {neu,lam,kp,pim,pi0}, {neu,s0,k0,pip,pim},
{neu,s0,k0,pi0,pi0}, {neu,s0,kp,pim,pi0}, {neu,sm,kp,pip,pim},
{neu,sm,kp,pi0,pi0}, {neu,sm,k0,pip,pi0}, {neu,sp,k0,pim,pi0},
{neu,sp,kp,pim,pim}, {pro,lam,k0,pim,pi0}, {pro,lam,kp,pim,pim},
{pro,s0,k0,pim,pi0}, {pro,s0,kp,pim,pim}, {pro,sm,kp,pim,pi0},
{pro,sm,k0,pip,pim}, {pro,sm,k0,pi0,pi0}, {pro,sp,k0,pim,pim},
{neu,neu,pim,kp,k0b}, {neu,neu,pip,k0,km}, {neu,neu,pi0,kp,km},
{neu,neu,pi0,k0,k0b}, {pro,neu,pim,kp,km}, {pro,neu,pim,k0,k0b},
{pro,neu,pi0,k0,km}, {pro,pro,pim,k0,km}}}

Definition at line 81 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_6bfs
staticprotected
Initial value:
=
{{{pro,pro,pip,pip,pim,pim},{pro,pro,pip,pim,pi0,pi0},
{pro,pro,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pi0},
{pro,neu,pip,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pim},
{neu,neu,pip,pip,pi0,pi0}},
{{neu,neu,pip,pip,pim,pim},{neu,neu,pip,pim,pi0,pi0},
{neu,neu,pi0,pi0,pi0,pi0},{pro,neu,pip,pim,pim,pi0},
{pro,neu,pim,pi0,pi0,pi0},{pro,pro,pip,pim,pim,pim},
{pro,pro,pim,pim,pi0,pi0}}}

Definition at line 82 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_7bfs
staticprotected
Initial value:
=
{{{pro,pro,pip,pip,pim,pim,pi0},{pro,pro,pip,pim,pi0,pi0,pi0},
{pro,pro,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pip,pim,pim},
{pro,neu,pip,pip,pim,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0},
{neu,neu,pip,pip,pip,pim,pi0},{neu,neu,pip,pip,pi0,pi0,pi0}},
{{neu,neu,pip,pip,pim,pim,pi0},{neu,neu,pip,pim,pi0,pi0,pi0},
{neu,neu,pi0,pi0,pi0,pi0,pi0},{pro,neu,pip,pip,pim,pim,pim},
{pro,neu,pip,pim,pim,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0},
{pro,pro,pip,pim,pim,pim,pi0},{pro,pro,pim,pim,pi0,pi0,pi0}}}

Definition at line 83 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_8bfs
staticprotected
Initial value:
=
{{{pro,pro,pip,pip,pip,pim,pim,pim},{pro,pro,pip,pip,pim,pim,pi0,pi0},
{pro,pro,pip,pim,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,neu,pip,pip,pip,pim,pim,pi0},{pro,neu,pip,pip,pim,pi0,pi0,pi0},
{pro,neu,pip,pi0,pi0,pi0,pi0,pi0},{neu,neu,pip,pip,pip,pip,pim,pim},
{neu,neu,pip,pip,pip,pim,pi0,pi0},{neu,neu,pip,pip,pi0,pi0,pi0,pi0}},
{{neu,neu,pip,pip,pip,pim,pim,pim},{neu,neu,pip,pip,pim,pim,pi0,pi0},
{neu,neu,pip,pim,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,neu,pip,pip,pim,pim,pim,pi0},{pro,neu,pip,pim,pim,pi0,pi0,pi0},
{pro,neu,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pip,pip,pim,pim,pim,pim},
{pro,pro,pip,pim,pim,pim,pi0,pi0},{pro,pro,pim,pim,pi0,pi0,pi0,pi0}}}

Definition at line 84 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

const G4int G4RPGNucleonInelastic::T1_9bfs
staticprotected
Initial value:
=
{{{pro,pro,pip,pip,pip,pim,pim,pim,pi0},{pro,pro,pip,pip,pim,pim,pi0,pi0,pi0},
{pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0},{pro,pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,neu,pip,pip,pip,pip,pim,pim,pim},{pro,neu,pip,pip,pip,pim,pim,pi0,pi0},
{pro,neu,pip,pip,pim,pi0,pi0,pi0,pi0},{pro,neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,neu,pip,pip,pip,pip,pim,pim,pi0},{neu,neu,pip,pip,pip,pim,pi0,pi0,pi0},
{neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0}},
{{neu,neu,pip,pip,pip,pim,pim,pim,pi0},{neu,neu,pip,pip,pim,pim,pi0,pi0,pi0},
{neu,neu,pip,pim,pi0,pi0,pi0,pi0,pi0},{neu,neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,neu,pip,pip,pip,pim,pim,pim,pim},{pro,neu,pip,pip,pim,pim,pim,pi0,pi0},
{pro,neu,pip,pim,pim,pi0,pi0,pi0,pi0},{pro,neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pro,pip,pip,pim,pim,pim,pim,pi0},{pro,pro,pip,pim,pim,pim,pi0,pi0,pi0},
{pro,pro,pip,pim,pi0,pi0,pi0,pi0,pi0}}}

Definition at line 85 of file G4RPGNucleonInelastic.hh.

Referenced by GetFSPartTypesForT1().

G4ThreadLocal G4double G4RPGNucleonInelastic::t1_dSigma_dMult
staticprotected

Definition at line 98 of file G4RPGNucleonInelastic.hh.

Referenced by G4RPGNucleonInelastic(), and GetMultiplicityT1().

The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:57:01 for Geant4.10 by   doxygen 1.8.7