Geant4  10.01.p01
G4ParticleHPInelastic.cc
Go to the documentation of this file.
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 // this code implementation is the intellectual property of
27 // neutron_hp -- source file
28 // J.P. Wellisch, Nov-1996
29 // A prototype of the low energy neutron transport model.
30 //
31 // By copying, distributing or modifying the Program (or any work
32 // based on the Program) you indicate your acceptance of this statement,
33 // and all its terms.
34 //
35 // $Id$
36 //
37 // 070523 bug fix for G4FPE_DEBUG on by A. Howard (and T. Koi)
38 // 081203 limit maximum trial for creating final states add protection for 1H isotope case by T. Koi
39 //
40 // P. Arce, June-2014 Conversion neutron_hp to particle_hp
41 //
42 #include "G4ParticleHPInelastic.hh"
43 #include "G4SystemOfUnits.hh"
44 
45 #include "G4ParticleHPManager.hh"
46 
47 G4ParticleHPInelastic::G4ParticleHPInelastic(G4ParticleDefinition* projectile, const char* name )
48  :G4HadronicInteraction(name), theProjectile(projectile)
49 {
50  const char* dataDirVariable;
51  if( projectile == G4Neutron::Neutron() ) {
52  dataDirVariable = "G4NEUTRONHPDATA";
53  }else if( projectile == G4Proton::Proton() ) {
54  dataDirVariable = "G4PROTONHPDATA";
55  }else if( projectile == G4Deuteron::Deuteron() ) {
56  dataDirVariable = "G4DEUTERONHPDATA";
57  }else if( projectile == G4Triton::Triton() ) {
58  dataDirVariable = "G4TRITONHPDATA";
59  }else if( projectile == G4He3::He3() ) {
60  dataDirVariable = "G4HE3HPDATA";
61  }else if( projectile == G4Alpha::Alpha() ) {
62  dataDirVariable = "G4ALPHAHPDATA";
63  } else {
64  G4String message("G4ParticleHPInelastic may only be called for neutron, proton, deuteron, triton, He3 or alpha, while it is called for " + projectile->GetParticleName());
65  throw G4HadronicException(__FILE__, __LINE__,message.c_str());
66  }
67 
68  SetMinEnergy( 0.0 );
69  SetMaxEnergy( 20.*MeV );
70 
71 // G4cout << " entering G4ParticleHPInelastic constructor"<<G4endl;
72  if(!getenv(dataDirVariable)){
73  G4String message("Please set the environement variable " + G4String(dataDirVariable) + " to point to the " + projectile->GetParticleName() + " cross-section files.");
74  throw G4HadronicException(__FILE__, __LINE__,message.c_str());
75  }
76  dirName = getenv(dataDirVariable);
77  G4cout << dirName << G4endl;
78 
79  G4String tString = "/Inelastic";
80  dirName = dirName + tString;
81  numEle = G4Element::GetNumberOfElements();
82 
83  G4cout << "@@@ G4ParticleHPInelastic instantiated for particle " << projectile->GetParticleName() << " data directory variable is " << dataDirVariable << " pointing to " << dirName << G4endl;
84 
85 /*
86  theInelastic = new G4ParticleHPChannelList[numEle];
87  for (G4int i=0; i<numEle; i++)
88  {
89  theInelastic[i].Init((*(G4Element::GetElementTable()))[i], dirName);
90  G4int itry = 0;
91  do
92  {
93  theInelastic[i].Register(&theNFS, "F01"); // has
94  theInelastic[i].Register(&theNXFS, "F02");
95  theInelastic[i].Register(&the2NDFS, "F03");
96  theInelastic[i].Register(&the2NFS, "F04"); // has, E Done
97  theInelastic[i].Register(&the3NFS, "F05"); // has, E Done
98  theInelastic[i].Register(&theNAFS, "F06");
99  theInelastic[i].Register(&theN3AFS, "F07");
100  theInelastic[i].Register(&the2NAFS, "F08");
101  theInelastic[i].Register(&the3NAFS, "F09");
102  theInelastic[i].Register(&theNPFS, "F10");
103  theInelastic[i].Register(&theN2AFS, "F11");
104  theInelastic[i].Register(&the2N2AFS, "F12");
105  theInelastic[i].Register(&theNDFS, "F13");
106  theInelastic[i].Register(&theNTFS, "F14");
107  theInelastic[i].Register(&theNHe3FS, "F15");
108  theInelastic[i].Register(&theND2AFS, "F16");
109  theInelastic[i].Register(&theNT2AFS, "F17");
110  theInelastic[i].Register(&the4NFS, "F18"); // has, E Done
111  theInelastic[i].Register(&the2NPFS, "F19");
112  theInelastic[i].Register(&the3NPFS, "F20");
113  theInelastic[i].Register(&theN2PFS, "F21");
114  theInelastic[i].Register(&theNPAFS, "F22");
115  theInelastic[i].Register(&thePFS, "F23");
116  theInelastic[i].Register(&theDFS, "F24");
117  theInelastic[i].Register(&theTFS, "F25");
118  theInelastic[i].Register(&theHe3FS, "F26");
119  theInelastic[i].Register(&theAFS, "F27");
120  theInelastic[i].Register(&the2AFS, "F28");
121  theInelastic[i].Register(&the3AFS, "F29");
122  theInelastic[i].Register(&the2PFS, "F30");
123  theInelastic[i].Register(&thePAFS, "F31");
124  theInelastic[i].Register(&theD2AFS, "F32");
125  theInelastic[i].Register(&theT2AFS, "F33");
126  theInelastic[i].Register(&thePDFS, "F34");
127  theInelastic[i].Register(&thePTFS, "F35");
128  theInelastic[i].Register(&theDAFS, "F36");
129  theInelastic[i].RestartRegistration();
130  itry++;
131  }
132  //while(!theInelastic[i].HasDataInAnyFinalState());
133  while( !theInelastic[i].HasDataInAnyFinalState() && itry < 6 );
134  // 6 is corresponding to the value(5) of G4ParticleHPChannel. TK
135 
136  if ( itry == 6 )
137  {
138  // No Final State at all.
139  G4bool exceptional = false;
140  if ( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() == 1 )
141  {
142  if ( (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() == 1 && (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() == 1 ) exceptional = true; //1H
143  }
144  if ( !exceptional ) throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this element");
145  }
146  }
147 */
148 
149  for (G4int i=0; i<numEle; i++)
150  {
151  theInelastic.push_back( new G4ParticleHPChannelList );
152  (*theInelastic[i]).Init((*(G4Element::GetElementTable()))[i], dirName, theProjectile);
153  G4int itry = 0;
154  do
155  {
156  (*theInelastic[i]).Register(&theNFS, "F01"); // has
157  (*theInelastic[i]).Register(&theNXFS, "F02");
158  (*theInelastic[i]).Register(&the2NDFS, "F03");
159  (*theInelastic[i]).Register(&the2NFS, "F04"); // has, E Done
160  (*theInelastic[i]).Register(&the3NFS, "F05"); // has, E Done
161  (*theInelastic[i]).Register(&theNAFS, "F06");
162  (*theInelastic[i]).Register(&theN3AFS, "F07");
163  (*theInelastic[i]).Register(&the2NAFS, "F08");
164  (*theInelastic[i]).Register(&the3NAFS, "F09");
165  (*theInelastic[i]).Register(&theNPFS, "F10");
166  (*theInelastic[i]).Register(&theN2AFS, "F11");
167  (*theInelastic[i]).Register(&the2N2AFS, "F12");
168  (*theInelastic[i]).Register(&theNDFS, "F13");
169  (*theInelastic[i]).Register(&theNTFS, "F14");
170  (*theInelastic[i]).Register(&theNHe3FS, "F15");
171  (*theInelastic[i]).Register(&theND2AFS, "F16");
172  (*theInelastic[i]).Register(&theNT2AFS, "F17");
173  (*theInelastic[i]).Register(&the4NFS, "F18"); // has, E Done
174  (*theInelastic[i]).Register(&the2NPFS, "F19");
175  (*theInelastic[i]).Register(&the3NPFS, "F20");
176  (*theInelastic[i]).Register(&theN2PFS, "F21");
177  (*theInelastic[i]).Register(&theNPAFS, "F22");
178  (*theInelastic[i]).Register(&thePFS, "F23");
179  (*theInelastic[i]).Register(&theDFS, "F24");
180  (*theInelastic[i]).Register(&theTFS, "F25");
181  (*theInelastic[i]).Register(&theHe3FS, "F26");
182  (*theInelastic[i]).Register(&theAFS, "F27");
183  (*theInelastic[i]).Register(&the2AFS, "F28");
184  (*theInelastic[i]).Register(&the3AFS, "F29");
185  (*theInelastic[i]).Register(&the2PFS, "F30");
186  (*theInelastic[i]).Register(&thePAFS, "F31");
187  (*theInelastic[i]).Register(&theD2AFS, "F32");
188  (*theInelastic[i]).Register(&theT2AFS, "F33");
189  (*theInelastic[i]).Register(&thePDFS, "F34");
190  (*theInelastic[i]).Register(&thePTFS, "F35");
191  (*theInelastic[i]).Register(&theDAFS, "F36");
192  (*theInelastic[i]).RestartRegistration();
193  itry++;
194  }
195  while( !(*theInelastic[i]).HasDataInAnyFinalState() && itry < 6 );
196  // 6 is corresponding to the value(5) of G4ParticleHPChannel. TK
197 
198  if ( itry == 6 )
199  {
200  // No Final State at all.
201  G4bool exceptional = false;
202  if ( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() == 1 )
203  {
204  if ( (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() == 1 && (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() == 1 ) exceptional = true; //1H
205  }
206  if ( !exceptional ) {
207  G4cerr << " ELEMENT Z " << (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() << " N " << (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() << G4endl; //1H
208 throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this element");
209  }
210  }
211 
212  }
213 
214  }
215 
216  G4ParticleHPInelastic::~G4ParticleHPInelastic()
217  {
218 // delete [] theInelastic;
219  for ( std::vector<G4ParticleHPChannelList*>::iterator
220  it = theInelastic.begin() ; it != theInelastic.end() ; it++ )
221  {
222  delete *it;
223  }
224  theInelastic.clear();
225  }
226 
227  #include "G4ParticleHPThermalBoost.hh"
228 
229  G4HadFinalState * G4ParticleHPInelastic::ApplyYourself(const G4HadProjectile& aTrack, G4Nucleus& aNucleus )
230  {
231  if ( numEle < (G4int)G4Element::GetNumberOfElements() ) addChannelForNewElement();
232  G4ParticleHPManager::GetInstance()->OpenReactionWhiteBoard();
233  const G4Material * theMaterial = aTrack.GetMaterial();
234  G4int n = theMaterial->GetNumberOfElements();
235  G4int index = theMaterial->GetElement(0)->GetIndex();
236  G4int it=0;
237  if(n!=1)
238  {
239  xSec = new G4double[n];
240  G4double sum=0;
241  G4int i;
242  const G4double * NumAtomsPerVolume = theMaterial->GetVecNbOfAtomsPerVolume();
243  G4double rWeight;
244  G4ParticleHPThermalBoost aThermalE;
245  for (i=0; i<n; i++)
246  {
247  index = theMaterial->GetElement(i)->GetIndex();
248  rWeight = NumAtomsPerVolume[i];
249  if(aTrack.GetDefinition() == G4Neutron::Neutron() ) {
250  xSec[i] = (*theInelastic[index]).GetXsec(aThermalE.GetThermalEnergy(aTrack,
251  theMaterial->GetElement(i),
252  theMaterial->GetTemperature()));
253  } else {
254  xSec[i] = (*theInelastic[index]).GetXsec(aTrack.GetKineticEnergy());
255  }
256  xSec[i] *= rWeight;
257  sum+=xSec[i];
258 #ifdef G4PHPDEBUG
259  if( getenv("G4ParticleHPDebug") ) G4cout << " G4ParticleHPInelastic XSEC ELEM " << i << " = " << xSec[i] << G4endl;
260 #endif
261 
262  }
263  G4double random = G4UniformRand();
264  G4double running = 0;
265  for (i=0; i<n; i++)
266  {
267  running += xSec[i];
268  index = theMaterial->GetElement(i)->GetIndex();
269  it = i;
270  //if(random<=running/sum) break;
271  if( sum == 0 || random<=running/sum) break;
272  }
273  delete [] xSec;
274  }
275 
276 #ifdef G4PHPDEBUG
277  if( getenv("G4ParticleHPDebug") ) G4cout << " G4ParticleHPInelastic SELECTED ELEM " << it << " = " << theMaterial->GetElement(it)->GetName() << " FROM MATERIAL " << theMaterial->GetName() << G4endl;
278 #endif
279  //return theInelastic[index].ApplyYourself(theMaterial->GetElement(it), aTrack);
280  G4HadFinalState* result = (*theInelastic[index]).ApplyYourself(theMaterial->GetElement(it), aTrack);
281  //
282  aNucleus.SetParameters(G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA(),G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargZ());
283  const G4Element* target_element = (*G4Element::GetElementTable())[index];
284  const G4Isotope* target_isotope=NULL;
285  G4int iele = target_element->GetNumberOfIsotopes();
286  for ( G4int j = 0 ; j != iele ; j++ ) {
287  target_isotope=target_element->GetIsotope( j );
288  if ( target_isotope->GetN() == G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA() ) break;
289  }
290  //G4cout << "Target Material of this reaction is " << theMaterial->GetName() << G4endl;
291  //G4cout << "Target Element of this reaction is " << target_element->GetName() << G4endl;
292  //G4cout << "Target Isotope of this reaction is " << target_isotope->GetName() << G4endl;
293  aNucleus.SetIsotope( target_isotope );
294 
295  G4ParticleHPManager::GetInstance()->CloseReactionWhiteBoard();
296 
297  //GDEB
298  if( getenv("G4PHPTEST") ) {
299  G4HadSecondary* seco = result->GetSecondary(0);
300  if(seco) {
301  G4ThreeVector secoMom = seco->GetParticle()->GetMomentum();
302  G4cout << " G4ParticleHPinelastic COS THETA " << std::cos(secoMom.theta()) <<" " << secoMom << G4endl;
303  }
304  }
305 
306  return result;
307  }
308 
309 const std::pair<G4double, G4double> G4ParticleHPInelastic::GetFatalEnergyCheckLevels() const
310 {
311  // max energy non-conservation is mass of heavy nucleus
312 // if ( getenv("G4PHP_DO_NOT_ADJUST_FINAL_STATE") ) return std::pair<G4double, G4double>(5*perCent,250*GeV);
313  // This should be same to the hadron default value
314 // return std::pair<G4double, G4double>(10*perCent,10*GeV);
315  return std::pair<G4double, G4double>(10*perCent,DBL_MAX);
316 }
317 
318 void G4ParticleHPInelastic::addChannelForNewElement()
319 {
320  for ( G4int i = numEle ; i < (G4int)G4Element::GetNumberOfElements() ; i++ )
321  {
322  G4cout << "G4ParticleHPInelastic Prepairing Data for the new element of " << (*(G4Element::GetElementTable()))[i]->GetName() << G4endl;
323 
324  theInelastic.push_back( new G4ParticleHPChannelList );
325  (*theInelastic[i]).Init((*(G4Element::GetElementTable()))[i], dirName, theProjectile);
326  G4int itry = 0;
327  do
328  {
329  (*theInelastic[i]).Register(&theNFS, "F01"); // has
330  (*theInelastic[i]).Register(&theNXFS, "F02");
331  (*theInelastic[i]).Register(&the2NDFS, "F03");
332  (*theInelastic[i]).Register(&the2NFS, "F04"); // has, E Done
333  (*theInelastic[i]).Register(&the3NFS, "F05"); // has, E Done
334  (*theInelastic[i]).Register(&theNAFS, "F06");
335  (*theInelastic[i]).Register(&theN3AFS, "F07");
336  (*theInelastic[i]).Register(&the2NAFS, "F08");
337  (*theInelastic[i]).Register(&the3NAFS, "F09");
338  (*theInelastic[i]).Register(&theNPFS, "F10");
339  (*theInelastic[i]).Register(&theN2AFS, "F11");
340  (*theInelastic[i]).Register(&the2N2AFS, "F12");
341  (*theInelastic[i]).Register(&theNDFS, "F13");
342  (*theInelastic[i]).Register(&theNTFS, "F14");
343  (*theInelastic[i]).Register(&theNHe3FS, "F15");
344  (*theInelastic[i]).Register(&theND2AFS, "F16");
345  (*theInelastic[i]).Register(&theNT2AFS, "F17");
346  (*theInelastic[i]).Register(&the4NFS, "F18"); // has, E Done
347  (*theInelastic[i]).Register(&the2NPFS, "F19");
348  (*theInelastic[i]).Register(&the3NPFS, "F20");
349  (*theInelastic[i]).Register(&theN2PFS, "F21");
350  (*theInelastic[i]).Register(&theNPAFS, "F22");
351  (*theInelastic[i]).Register(&thePFS, "F23");
352  (*theInelastic[i]).Register(&theDFS, "F24");
353  (*theInelastic[i]).Register(&theTFS, "F25");
354  (*theInelastic[i]).Register(&theHe3FS, "F26");
355  (*theInelastic[i]).Register(&theAFS, "F27");
356  (*theInelastic[i]).Register(&the2AFS, "F28");
357  (*theInelastic[i]).Register(&the3AFS, "F29");
358  (*theInelastic[i]).Register(&the2PFS, "F30");
359  (*theInelastic[i]).Register(&thePAFS, "F31");
360  (*theInelastic[i]).Register(&theD2AFS, "F32");
361  (*theInelastic[i]).Register(&theT2AFS, "F33");
362  (*theInelastic[i]).Register(&thePDFS, "F34");
363  (*theInelastic[i]).Register(&thePTFS, "F35");
364  (*theInelastic[i]).Register(&theDAFS, "F36");
365  (*theInelastic[i]).RestartRegistration();
366  itry++;
367  }
368  while( !(*theInelastic[i]).HasDataInAnyFinalState() && itry < 6 );
369  // 6 is corresponding to the value(5) of G4ParticleHPChannel. TK
370 
371  if ( itry == 6 )
372  {
373  // No Final State at all.
374  G4bool exceptional = false;
375  if ( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() == 1 )
376  {
377  if ( (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() == 1 && (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() == 1 ) exceptional = true; //1H
378  }
379  if ( !exceptional ) throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this element");
380  }
381  }
382 
383  numEle = (G4int)G4Element::GetNumberOfElements();
384 }
385 
386 G4int G4ParticleHPInelastic::GetVerboseLevel() const
387 {
388  return G4ParticleHPManager::GetInstance()->GetVerboseLevel();
389 }
390 void G4ParticleHPInelastic::SetVerboseLevel( G4int newValue )
391 {
392  G4ParticleHPManager::GetInstance()->SetVerboseLevel(newValue);
393 }
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:46
G4Element::GetNumberOfIsotopes
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:158
G4Nucleus::SetIsotope
void SetIsotope(const G4Isotope *iso)
Definition: G4Nucleus.hh:122
lightions::Init
void Init()
Definition: G4IonTable.cc:87
MeV
static const double MeV
Definition: G4SIunits.hh:193
G4ParticleHPInelastic.hh
G4ParticleHPInelastic::theDFS
G4ParticleHPDInelasticFS theDFS
Definition: G4ParticleHPInelastic.hh:128
G4ParticleHPInelastic::theND2AFS
G4ParticleHPND2AInelasticFS theND2AFS
Definition: G4ParticleHPInelastic.hh:134
G4ParticleHPInelastic::theProjectile
G4ParticleDefinition * theProjectile
Definition: G4ParticleHPInelastic.hh:150
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:79
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4ParticleHPInelastic::theDAFS
G4ParticleHPDAInelasticFS theDAFS
Definition: G4ParticleHPInelastic.hh:127
G4ParticleHPInelastic::theN2AFS
G4ParticleHPN2AInelasticFS theN2AFS
Definition: G4ParticleHPInelastic.hh:130
G4ParticleHPInelastic::the3AFS
G4ParticleHP3AInelasticFS the3AFS
Definition: G4ParticleHPInelastic.hh:120
G4ParticleHPInelastic::theNHe3FS
G4ParticleHPNHe3InelasticFS theNHe3FS
Definition: G4ParticleHPInelastic.hh:136
G4ParticleHPInelastic::theHe3FS
G4ParticleHPHe3InelasticFS theHe3FS
Definition: G4ParticleHPInelastic.hh:129
name
G4String name
Definition: TRTMaterials.hh:40
G4ParticleHPInelastic::theN3AFS
G4ParticleHPN3AInelasticFS theN3AFS
Definition: G4ParticleHPInelastic.hh:132
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:178
G4ParticleHPInelastic::the2N2AFS
G4ParticleHP2N2AInelasticFS the2N2AFS
Definition: G4ParticleHPInelastic.hh:114
G4ParticleHPInelastic::theNPAFS
G4ParticleHPNPAInelasticFS theNPAFS
Definition: G4ParticleHPInelastic.hh:138
G4ParticleHPInelastic::the3NAFS
G4ParticleHP3NAInelasticFS the3NAFS
Definition: G4ParticleHPInelastic.hh:121
G4ParticleHPInelastic::theNFS
G4ParticleHPNInelasticFS theNFS
Definition: G4ParticleHPInelastic.hh:137
G4ParticleHPInelastic::theNT2AFS
G4ParticleHPNT2AInelasticFS theNT2AFS
Definition: G4ParticleHPInelastic.hh:140
G4ParticleHPInelastic::theNTFS
G4ParticleHPNTInelasticFS theNTFS
Definition: G4ParticleHPInelastic.hh:141
G4ParticleHPInelastic::theInelastic
std::vector< G4ParticleHPChannelList * > theInelastic
Definition: G4ParticleHPInelastic.hh:106
G4Material::GetElement
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:202
G4ParticleHPInelastic::the3NFS
G4ParticleHP3NInelasticFS the3NFS
Definition: G4ParticleHPInelastic.hh:122
G4int
int G4int
Definition: G4Types.hh:78
G4ParticleHPInelastic::G4ParticleHPInelastic
G4ParticleHPInelastic(G4ParticleDefinition *projectile=G4Neutron::Neutron(), const char *name="NeutronHPInelastic")
Definition: G4ParticleHPInelastic.cc:47
G4ParticleHPInelastic::thePDFS
G4ParticleHPPDInelasticFS thePDFS
Definition: G4ParticleHPInelastic.hh:144
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:159
G4ParticleHPInelastic::thePFS
G4ParticleHPPInelasticFS thePFS
Definition: G4ParticleHPInelastic.hh:145
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:89
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:206
G4AutoDelete::Register
void Register(T *inst)
Definition: G4AutoDelete.hh:65
G4ParticleHPInelastic::theNPFS
G4ParticleHPNPInelasticFS theNPFS
Definition: G4ParticleHPInelastic.hh:139
G4Isotope::GetN
G4int GetN() const
Definition: G4Isotope.hh:94
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:95
G4cout
G4GLOB_DLL std::ostream G4cout
G4ParticleHPInelastic::the3NPFS
G4ParticleHP3NPInelasticFS the3NPFS
Definition: G4ParticleHPInelastic.hh:123
G4Element::GetNumberOfElements
static size_t GetNumberOfElements()
Definition: G4Element.cc:410
G4HadProjectile::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4HadProjectile.hh:81
G4ParticleHPInelastic::thePAFS
G4ParticleHPPAInelasticFS thePAFS
Definition: G4ParticleHPInelastic.hh:143
G4bool
bool G4bool
Definition: G4Types.hh:79
G4HadProjectile::GetKineticEnergy
G4double GetKineticEnergy() const
Definition: G4HadProjectile.hh:106
G4ParticleHPInelastic::the2AFS
G4ParticleHP2AInelasticFS the2AFS
Definition: G4ParticleHPInelastic.hh:113
G4ParticleHPInelastic::theNDFS
G4ParticleHPNDInelasticFS theNDFS
Definition: G4ParticleHPInelastic.hh:135
G4Triton::Triton
static G4Triton * Triton()
Definition: G4Triton.cc:95
G4ParticleHPInelastic::thePTFS
G4ParticleHPPTInelasticFS thePTFS
Definition: G4ParticleHPInelastic.hh:146
G4Element::GetIndex
size_t GetIndex() const
Definition: G4Element.hh:181
G4ParticleHPInelastic::theTFS
G4ParticleHPTInelasticFS theTFS
Definition: G4ParticleHPInelastic.hh:148
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93
perCent
static const double perCent
Definition: G4SIunits.hh:296
G4ParticleHPInelastic::the2NDFS
G4ParticleHP2NDInelasticFS the2NDFS
Definition: G4ParticleHPInelastic.hh:116
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
n
const G4int n
Definition: G4UrQMD1_3Interface.hh:144
G4Deuteron::Deuteron
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:94
G4ParticleHPInelastic::theT2AFS
G4ParticleHPT2AInelasticFS theT2AFS
Definition: G4ParticleHPInelastic.hh:147
G4ParticleHPThermalBoost::GetThermalEnergy
G4double GetThermalEnergy(const G4HadProjectile &aP, const G4Element *anE, G4double aT)
Definition: G4ParticleHPThermalBoost.hh:45
G4ParticleHPInelastic::the4NFS
G4ParticleHP4NInelasticFS the4NFS
Definition: G4ParticleHPInelastic.hh:124
G4ParticleHPInelastic::theNAFS
G4ParticleHPNAInelasticFS theNAFS
Definition: G4ParticleHPInelastic.hh:133
G4ParticleHPInelastic::the2NFS
G4ParticleHP2NInelasticFS the2NFS
Definition: G4ParticleHPInelastic.hh:117
G4ParticleHPInelastic::theN2PFS
G4ParticleHPN2PInelasticFS theN2PFS
Definition: G4ParticleHPInelastic.hh:131
G4ParticleHPInelastic::theNXFS
G4ParticleHPNXInelasticFS theNXFS
Definition: G4ParticleHPInelastic.hh:142
G4ParticleHPInelastic::theAFS
G4ParticleHPAInelasticFS theAFS
Definition: G4ParticleHPInelastic.hh:125
G4HadSecondary::GetParticle
G4DynamicParticle * GetParticle()
Definition: G4HadSecondary.hh:42
G4ParticleHPInelastic::the2PFS
G4ParticleHP2PInelasticFS the2PFS
Definition: G4ParticleHPInelastic.hh:119
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:102
G4Element::GetIsotope
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:169
G4Material::GetTemperature
G4double GetTemperature() const
Definition: G4Material.hh:182
G4endl
#define G4endl
Definition: G4ios.hh:61
G4HadProjectile::GetMaterial
const G4Material * GetMaterial() const
Definition: G4HadProjectile.hh:76
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
G4Alpha::Alpha
static G4Alpha * Alpha()
Definition: G4Alpha.cc:89
G4double
double G4double
Definition: G4Types.hh:76
G4Element::GetName
const G4String & GetName() const
Definition: G4Element.hh:127
G4SystemOfUnits.hh
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:403
G4ParticleHPManager::GetReactionWhiteBoard
G4ParticleHPReactionWhiteBoard * GetReactionWhiteBoard()
Definition: G4ParticleHPManager.cc:55
G4ParticleHPInelastic::the2NAFS
G4ParticleHP2NAInelasticFS the2NAFS
Definition: G4ParticleHPInelastic.hh:115
G4ParticleHPInelastic::ApplyYourself
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
Definition: G4ParticleHPInelastic.cc:229
G4He3::He3
static G4He3 * He3()
Definition: G4He3.cc:94
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4ParticleHPInelastic::the2NPFS
G4ParticleHP2NPInelasticFS the2NPFS
Definition: G4ParticleHPInelastic.hh:118
G4Nucleus::SetParameters
void SetParameters(const G4double A, const G4double Z)
Definition: G4Nucleus.cc:198
G4ParticleHPThermalBoost.hh
G4ParticleHPInelastic::theD2AFS
G4ParticleHPD2AInelasticFS theD2AFS
Definition: G4ParticleHPInelastic.hh:126
G4DynamicParticle::GetMomentum
G4ThreeVector GetMomentum() const
G4cerr
G4GLOB_DLL std::ostream G4cerr
G4ParticleHPInelastic::GetFatalEnergyCheckLevels
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels() const
Definition: G4ParticleHPInelastic.cc:309
G4ParticleHPManager.hh