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
G4RPGPionInelastic Class Reference

#include <G4RPGPionInelastic.hh>

Inheritance diagram for G4RPGPionInelastic:
G4RPGInelastic G4HadronicInteraction G4RPGPiMinusInelastic G4RPGPiPlusInelastic

Public Member Functions

 G4RPGPionInelastic (const G4String &modelName="RPGPionInelastic")
 
 ~G4RPGPionInelastic ()
 
- 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 GetMultiplicityT12 (G4double KE) const
 
G4int GetMultiplicityT32 (G4double KE) const
 
std::vector< G4int > GetFSPartTypesForT32 (G4int mult, G4double KE, G4int tindex) const
 
std::vector< G4int > GetFSPartTypesForT12 (G4int mult, G4double KE, G4int tindex) const
 
std::vector< G4int > GetFSPartTypesForPipP (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForPimN (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForPipN (G4int mult, G4double KE) const
 
std::vector< G4int > GetFSPartTypesForPimP (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 pipPindex [8][2]
 
static const G4int pimPindex [8][2]
 
static const G4int T32_2bfs [2][2][2]
 
static const G4int T32_3bfs [2][7][3]
 
static const G4int T32_4bfs [2][15][4]
 
static const G4int T32_5bfs [2][24][5]
 
static const G4int T32_6bfs [2][5][6]
 
static const G4int T32_7bfs [2][6][7]
 
static const G4int T32_8bfs [2][7][8]
 
static const G4int T32_9bfs [2][8][9]
 
static const G4int T12_2bfs [2][5][2]
 
static const G4int T12_3bfs [2][13][3]
 
static const G4int T12_4bfs [2][22][4]
 
static const G4int T12_5bfs [2][31][5]
 
static const G4int T12_6bfs [2][6][6]
 
static const G4int T12_7bfs [2][7][7]
 
static const G4int T12_8bfs [2][8][8]
 
static const G4int T12_9bfs [2][9][9]
 
static G4ThreadLocal G4double pipPtot [30]
 
static G4ThreadLocal G4double pimPtot [30]
 
static G4ThreadLocal G4double t12_dSigma_dMult [8][30]
 
static G4ThreadLocal G4double t32_dSigma_dMult [8][30]
 
static const G4float pipPCrossSections [74][30]
 
static const G4float pimPCrossSections [101][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 G4RPGPionInelastic.hh.

Constructor & Destructor Documentation

G4RPGPionInelastic::G4RPGPionInelastic ( const G4String &  modelName = "RPGPionInelastic")

Definition at line 33 of file G4RPGPionInelastic.cc.

References python.hepunit::GeV, pimPCrossSections, pimPindex, pimPtot, pipPCrossSections, pipPindex, pipPtot, G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), t12_dSigma_dMult, and t32_dSigma_dMult.

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

Definition at line 48 of file G4RPGPionInelastic.hh.

48 { }

Member Function Documentation

std::vector<G4int> G4RPGPionInelastic::GetFSPartTypesForPimN ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 65 of file G4RPGPionInelastic.hh.

References GetFSPartTypesForT32().

66  {return GetFSPartTypesForT32(mult, KE, 1); }
G4RPGPionInelastic::GetFSPartTypesForT32
std::vector< G4int > GetFSPartTypesForT32(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGPionInelastic.cc:138
std::vector<G4int> G4RPGPionInelastic::GetFSPartTypesForPimP ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 71 of file G4RPGPionInelastic.hh.

References GetFSPartTypesForT12().

72  {return GetFSPartTypesForT12(mult, KE, 0); }
G4RPGPionInelastic::GetFSPartTypesForT12
std::vector< G4int > GetFSPartTypesForT12(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGPionInelastic.cc:186
std::vector<G4int> G4RPGPionInelastic::GetFSPartTypesForPipN ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 68 of file G4RPGPionInelastic.hh.

References GetFSPartTypesForT12().

69  {return GetFSPartTypesForT12(mult, KE, 1); }
G4RPGPionInelastic::GetFSPartTypesForT12
std::vector< G4int > GetFSPartTypesForT12(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGPionInelastic.cc:186
std::vector<G4int> G4RPGPionInelastic::GetFSPartTypesForPipP ( G4int  mult,
G4double  KE 
) const
inlineprotected

Definition at line 62 of file G4RPGPionInelastic.hh.

References GetFSPartTypesForT32().

63  {return GetFSPartTypesForT32(mult, KE, 0); }
G4RPGPionInelastic::GetFSPartTypesForT32
std::vector< G4int > GetFSPartTypesForT32(G4int mult, G4double KE, G4int tindex) const
Definition: G4RPGPionInelastic.cc:138
std::vector< G4int > G4RPGPionInelastic::GetFSPartTypesForT12 ( G4int  mult,
G4double  KE,
G4int  tindex 
) const
protected

Definition at line 186 of file G4RPGPionInelastic.cc.

References G4cout, G4endl, G4RPGInelastic::interpolateEnergy(), pimPCrossSections, pimPindex, G4RPGInelastic::sampleFlat(), T12_2bfs, T12_3bfs, T12_4bfs, T12_5bfs, T12_6bfs, T12_7bfs, T12_8bfs, and T12_9bfs.

Referenced by GetFSPartTypesForPimP(), and GetFSPartTypesForPipN().

187 {
188  G4int i;
189  G4double sigint(0.);
190  std::vector<G4double> sigma;
191 
192  std::pair<G4int, G4double> epair = interpolateEnergy(KE);
193  G4int k = epair.first;
194  G4double fraction = epair.second;
195 
196  G4int start = pimPindex[mult-2][0];
197  G4int stop = pimPindex[mult-2][1];
198 
199  for(i = start; i < stop; i++) {
200  sigint = pimPCrossSections[i][k]
201  + fraction*(pimPCrossSections[i][k+1] - pimPCrossSections[i][k]);
202  sigma.push_back(sigint);
203  }
204 
205  G4int channel = sampleFlat(sigma);
206 
207  std::vector<G4int> kinds;
208 
209  if (mult == 2) {
210  for(i = 0; i < mult; i++) kinds.push_back(T12_2bfs[tzindex][channel][i]);
211  } else if (mult == 3) {
212  for(i = 0; i < mult; i++) kinds.push_back(T12_3bfs[tzindex][channel][i]);
213  } else if (mult == 4) {
214  for(i = 0; i < mult; i++) kinds.push_back(T12_4bfs[tzindex][channel][i]);
215  } else if (mult == 5) {
216  for(i = 0; i < mult; i++) kinds.push_back(T12_5bfs[tzindex][channel][i]);
217  } else if (mult == 6) {
218  for(i = 0; i < mult; i++) kinds.push_back(T12_6bfs[tzindex][channel][i]);
219  } else if (mult == 7) {
220  for(i = 0; i < mult; i++) kinds.push_back(T12_7bfs[tzindex][channel][i]);
221  } else if (mult == 8) {
222  for(i = 0; i < mult; i++) kinds.push_back(T12_8bfs[tzindex][channel][i]);
223  } else if (mult == 9) {
224  for(i = 0; i < mult; i++) kinds.push_back(T12_9bfs[tzindex][channel][i]);
225  } else {
226  G4cout << " Illegal multiplicity " << G4endl;
227  }
228 
229  return kinds;
230 }
G4RPGPionInelastic::T12_2bfs
static const G4int T12_2bfs[2][5][2]
Definition: G4RPGPionInelastic.hh:86
G4RPGPionInelastic::T12_3bfs
static const G4int T12_3bfs[2][13][3]
Definition: G4RPGPionInelastic.hh:87
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4RPGPionInelastic::T12_8bfs
static const G4int T12_8bfs[2][8][8]
Definition: G4RPGPionInelastic.hh:92
G4int
int G4int
Definition: G4Types.hh:78
G4cout
G4GLOB_DLL std::ostream G4cout
G4RPGPionInelastic::pimPindex
static const G4int pimPindex[8][2]
Definition: G4RPGPionInelastic.hh:75
G4RPGPionInelastic::T12_4bfs
static const G4int T12_4bfs[2][22][4]
Definition: G4RPGPionInelastic.hh:88
G4RPGPionInelastic::T12_6bfs
static const G4int T12_6bfs[2][6][6]
Definition: G4RPGPionInelastic.hh:90
G4RPGPionInelastic::pimPCrossSections
static const G4float pimPCrossSections[101][30]
Definition: G4RPGPionInelastic.hh:101
G4endl
#define G4endl
Definition: G4ios.hh:61
G4RPGPionInelastic::T12_7bfs
static const G4int T12_7bfs[2][7][7]
Definition: G4RPGPionInelastic.hh:91
G4double
double G4double
Definition: G4Types.hh:76
G4RPGPionInelastic::T12_9bfs
static const G4int T12_9bfs[2][9][9]
Definition: G4RPGPionInelastic.hh:93
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGPionInelastic::T12_5bfs
static const G4int T12_5bfs[2][31][5]
Definition: G4RPGPionInelastic.hh:89
std::vector< G4int > G4RPGPionInelastic::GetFSPartTypesForT32 ( G4int  mult,
G4double  KE,
G4int  tindex 
) const
protected

Definition at line 138 of file G4RPGPionInelastic.cc.

References G4cout, G4endl, G4RPGInelastic::interpolateEnergy(), pipPCrossSections, pipPindex, G4RPGInelastic::sampleFlat(), T32_2bfs, T32_3bfs, T32_4bfs, T32_5bfs, T32_6bfs, T32_7bfs, T32_8bfs, and T32_9bfs.

Referenced by GetFSPartTypesForPimN(), and GetFSPartTypesForPipP().

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 = pipPindex[mult-2][0];
149  G4int stop = pipPindex[mult-2][1];
150 
151  for(i = start; i < stop; i++) {
152  sigint = pipPCrossSections[i][k]
153  + fraction*(pipPCrossSections[i][k+1] - pipPCrossSections[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(T32_2bfs[tzindex][channel][i]);
163  } else if (mult == 3) {
164  for(i = 0; i < mult; i++) kinds.push_back(T32_3bfs[tzindex][channel][i]);
165  } else if (mult == 4) {
166  for(i = 0; i < mult; i++) kinds.push_back(T32_4bfs[tzindex][channel][i]);
167  } else if (mult == 5) {
168  for(i = 0; i < mult; i++) kinds.push_back(T32_5bfs[tzindex][channel][i]);
169  } else if (mult == 6) {
170  for(i = 0; i < mult; i++) kinds.push_back(T32_6bfs[tzindex][channel][i]);
171  } else if (mult == 7) {
172  for(i = 0; i < mult; i++) kinds.push_back(T32_7bfs[tzindex][channel][i]);
173  } else if (mult == 8) {
174  for(i = 0; i < mult; i++) kinds.push_back(T32_8bfs[tzindex][channel][i]);
175  } else if (mult == 9) {
176  for(i = 0; i < mult; i++) kinds.push_back(T32_9bfs[tzindex][channel][i]);
177  } else {
178  G4cout << " Illegal multiplicity " << G4endl;
179  }
180 
181  return kinds;
182 }
G4RPGPionInelastic::T32_3bfs
static const G4int T32_3bfs[2][7][3]
Definition: G4RPGPionInelastic.hh:78
G4RPGPionInelastic::T32_4bfs
static const G4int T32_4bfs[2][15][4]
Definition: G4RPGPionInelastic.hh:79
G4RPGPionInelastic::pipPindex
static const G4int pipPindex[8][2]
Definition: G4RPGPionInelastic.hh:74
G4RPGInelastic::sampleFlat
G4int sampleFlat(std::vector< G4double > sigma) const
Definition: G4RPGInelastic.cc:523
G4int
int G4int
Definition: G4Types.hh:78
G4RPGPionInelastic::T32_5bfs
static const G4int T32_5bfs[2][24][5]
Definition: G4RPGPionInelastic.hh:80
G4cout
G4GLOB_DLL std::ostream G4cout
G4RPGPionInelastic::T32_6bfs
static const G4int T32_6bfs[2][5][6]
Definition: G4RPGPionInelastic.hh:81
G4RPGPionInelastic::pipPCrossSections
static const G4float pipPCrossSections[74][30]
Definition: G4RPGPionInelastic.hh:100
G4RPGPionInelastic::T32_8bfs
static const G4int T32_8bfs[2][7][8]
Definition: G4RPGPionInelastic.hh:83
G4RPGPionInelastic::T32_7bfs
static const G4int T32_7bfs[2][6][7]
Definition: G4RPGPionInelastic.hh:82
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76
G4RPGPionInelastic::T32_2bfs
static const G4int T32_2bfs[2][2][2]
Definition: G4RPGPionInelastic.hh:77
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4RPGPionInelastic::T32_9bfs
static const G4int T32_9bfs[2][8][9]
Definition: G4RPGPionInelastic.hh:84
G4int G4RPGPionInelastic::GetMultiplicityT12 ( G4double  KE) const
protected

Definition at line 99 of file G4RPGPionInelastic.cc.

References G4RPGInelastic::interpolateEnergy(), G4RPGInelastic::sampleFlat(), and t12_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 = t12_dSigma_dMult[j][k]
110  + fraction*(t12_dSigma_dMult[j][k+1] - t12_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
G4RPGPionInelastic::t12_dSigma_dMult
static G4ThreadLocal G4double t12_dSigma_dMult[8][30]
Definition: G4RPGPionInelastic.hh:97
G4double
double G4double
Definition: G4Types.hh:76
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506
G4int G4RPGPionInelastic::GetMultiplicityT32 ( G4double  KE) const
protected

Definition at line 118 of file G4RPGPionInelastic.cc.

References G4RPGInelastic::interpolateEnergy(), G4RPGInelastic::sampleFlat(), and t32_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 = t32_dSigma_dMult[j][k]
129  + fraction*(t32_dSigma_dMult[j][k+1] - t32_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
G4RPGPionInelastic::t32_dSigma_dMult
static G4ThreadLocal G4double t32_dSigma_dMult[8][30]
Definition: G4RPGPionInelastic.hh:98
G4double
double G4double
Definition: G4Types.hh:76
G4RPGInelastic::interpolateEnergy
std::pair< G4int, G4double > interpolateEnergy(G4double ke) const
Definition: G4RPGInelastic.cc:506

Field Documentation

const G4float G4RPGPionInelastic::pimPCrossSections
staticprotected

Definition at line 101 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic(), and GetFSPartTypesForT12().

const G4int G4RPGPionInelastic::pimPindex
staticprotected
Initial value:
=
{{0, 4}, {5,17}, {18,39}, {40,70}, {71,76}, {77,83}, {84,91}, {92,100}}

Definition at line 75 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic(), and GetFSPartTypesForT12().

G4ThreadLocal G4double G4RPGPionInelastic::pimPtot
staticprotected

Definition at line 96 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic().

const G4float G4RPGPionInelastic::pipPCrossSections
staticprotected

Definition at line 100 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic(), and GetFSPartTypesForT32().

const G4int G4RPGPionInelastic::pipPindex
staticprotected
Initial value:
=
{{0, 1}, {2, 8}, {9,23}, {24,47}, {48,52}, {53,58}, {59,65}, {66,73}}

Definition at line 74 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic(), and GetFSPartTypesForT32().

G4ThreadLocal G4double G4RPGPionInelastic::pipPtot
staticprotected

Definition at line 95 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic().

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

Definition at line 86 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

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

Definition at line 87 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

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

Definition at line 88 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

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

Definition at line 89 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

const G4int G4RPGPionInelastic::T12_6bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pim,pim,pim}, {pro,pip,pim,pim,pi0,pi0},
{pro,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pim,pim,pi0},
{neu,pip,pim,pi0,pi0,pi0}, {neu,pi0,pi0,pi0,pi0,pi0}},
{{neu,pip,pip,pip,pim,pim}, {neu,pip,pip,pim,pi0,pi0},
{neu,pip,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pi0},
{pro,pip,pim,pi0,pi0,pi0}, {pro,pi0,pi0,pi0,pi0,pi0}}}

Definition at line 90 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

const G4int G4RPGPionInelastic::T12_7bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pim,pim,pim,pi0}, {pro,pip,pim,pim,pi0,pi0,pi0},
{pro,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim},
{neu,pip,pip,pim,pim,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0},
{neu,pi0,pi0,pi0,pi0,pi0,pi0}},
{{neu,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pim,pi0,pi0,pi0},
{neu,pip,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim},
{pro,pip,pip,pim,pim,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0},
{pro,pi0,pi0,pi0,pi0,pi0,pi0}}}

Definition at line 91 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

const G4int G4RPGPionInelastic::T12_8bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pip,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pi0,pi0},
{pro,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0},
{neu,pip,pip,pim,pim,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0}},
{{neu,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pim,pim,pi0,pi0},
{neu,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0},
{pro,pip,pip,pim,pim,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0}}}

Definition at line 92 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

const G4int G4RPGPionInelastic::T12_9bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pip,pim,pim,pim,pi0,pi0,pi0},
{pro,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pim,pim,pim,pi0,pi0},
{neu,pip,pip,pip,pip,pim,pim,pim,pim}},
{{neu,pip,pip,pip,pip,pim,pim,pim,pi0}, {neu,pip,pip,pip,pim,pim,pi0,pi0,pi0},
{neu,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pi0,pi0,pi0,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pim,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pip,pip,pim,pim,pi0,pi0,pi0,pi0}, {pro,pip,pip,pip,pim,pim,pim,pi0,pi0},
{pro,pip,pip,pip,pip,pim,pim,pim,pim}}}

Definition at line 93 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT12().

G4ThreadLocal G4double G4RPGPionInelastic::t12_dSigma_dMult
staticprotected

Definition at line 97 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic(), and GetMultiplicityT12().

const G4int G4RPGPionInelastic::T32_2bfs
staticprotected
Initial value:
=
{{{pro,pip}, {sp,kp}},
{{neu,pim}, {sm,k0}}}

Definition at line 77 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

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

Definition at line 78 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

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

Definition at line 79 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

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

Definition at line 80 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

const G4int G4RPGPionInelastic::T32_6bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pip,pim,pim}, {pro,pip,pip,pim,pi0,pi0},
{pro,pip,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0},
{neu,pip,pip,pip,pim,pi0}},
{{neu,pip,pip,pim,pim,pim}, {neu,pip,pim,pim,pi0,pi0},
{neu,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0},
{pro,pip,pim,pim,pim,pi0}}}

Definition at line 81 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

const G4int G4RPGPionInelastic::T32_7bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pip,pim,pim,pi0}, {pro,pip,pip,pim,pi0,pi0,pi0},
{pro,pip,pi0,pi0,pi0,pi0,pi0}, {neu,pip,pip,pip,pip,pim,pim},
{neu,pip,pip,pip,pim,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0}},
{{neu,pip,pip,pim,pim,pim,pi0}, {neu,pip,pim,pim,pi0,pi0,pi0},
{neu,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pip,pim,pim,pim,pim},
{pro,pip,pim,pim,pim,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0}}}

Definition at line 82 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

const G4int G4RPGPionInelastic::T32_8bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pip,pip,pim,pim,pim}, {pro,pip,pip,pip,pim,pim,pi0,pi0},
{pro,pip,pip,pim,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,pip,pip,pip,pip,pim,pim,pi0}, {neu,pip,pip,pip,pim,pi0,pi0,pi0},
{neu,pip,pip,pi0,pi0,pi0,pi0,pi0}},
{{neu,pip,pip,pip,pim,pim,pim,pim}, {neu,pip,pip,pim,pim,pim,pi0,pi0},
{neu,pip,pim,pim,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pip,pip,pim,pim,pim,pim,pi0}, {pro,pip,pim,pim,pim,pi0,pi0,pi0},
{pro,pim,pim,pi0,pi0,pi0,pi0,pi0}}}

Definition at line 83 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

const G4int G4RPGPionInelastic::T32_9bfs
staticprotected
Initial value:
=
{{{pro,pip,pip,pip,pip,pim,pim,pim,pi0}, {pro,pip,pip,pip,pim,pim,pi0,pi0,pi0},
{pro,pip,pip,pim,pi0,pi0,pi0,pi0,pi0}, {pro,pip,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{neu,pip,pip,pip,pip,pip,pim,pim,pim}, {neu,pip,pip,pip,pip,pim,pim,pi0,pi0},
{neu,pip,pip,pip,pim,pi0,pi0,pi0,pi0}, {neu,pip,pip,pi0,pi0,pi0,pi0,pi0,pi0}},
{{neu,pip,pip,pip,pim,pim,pim,pim,pi0}, {neu,pip,pip,pim,pim,pim,pi0,pi0,pi0},
{neu,pip,pim,pim,pi0,pi0,pi0,pi0,pi0}, {neu,pim,pi0,pi0,pi0,pi0,pi0,pi0,pi0},
{pro,pip,pip,pip,pim,pim,pim,pim,pim}, {pro,pip,pip,pim,pim,pim,pim,pi0,pi0},
{pro,pip,pim,pim,pim,pi0,pi0,pi0,pi0}, {pro,pim,pim,pi0,pi0,pi0,pi0,pi0,pi0}}}

Definition at line 84 of file G4RPGPionInelastic.hh.

Referenced by GetFSPartTypesForT32().

G4ThreadLocal G4double G4RPGPionInelastic::t32_dSigma_dMult
staticprotected

Definition at line 98 of file G4RPGPionInelastic.hh.

Referenced by G4RPGPionInelastic(), and GetMultiplicityT32().

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