Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Attributes | List of all members
G4ParticleHPInelastic Class Reference

#include <G4ParticleHPInelastic.hh>

Inheritance diagram for G4ParticleHPInelastic:
Inheritance graph
[legend]
Collaboration diagram for G4ParticleHPInelastic:
Collaboration graph
[legend]

Public Member Functions

 G4ParticleHPInelastic (G4ParticleDefinition *projectile=G4Neutron::Neutron(), const char *name="NeutronHPInelastic")
 
 ~G4ParticleHPInelastic ()
 
G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
virtual const std::pair
< G4double, G4double > 		GetFatalEnergyCheckLevels () const
 
G4int GetVerboseLevel () const
 
void SetVerboseLevel (G4int)
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void ModelDescription (std::ostream &outFile) const
 
- Public Member Functions inherited from G4HadronicInteraction
 G4HadronicInteraction (const G4String &modelName="HadronicModel")
 
virtual ~G4HadronicInteraction ()
 
virtual G4double SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
 
virtual G4bool IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
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)
 
G4int GetVerboseLevel () const
 
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
 
virtual std::pair< G4double,
G4double > 		GetEnergyMomentumCheckLevels () const
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
virtual void InitialiseModel ()
 

Protected Attributes

std::vector
< G4ParticleHPChannelList * > * 		theInelastic
 
G4String dataDirVariable
 
G4String dirName
 
G4int numEle
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Additional Inherited Members

- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 

Detailed Description

Definition at line 89 of file G4ParticleHPInelastic.hh.

Constructor & Destructor Documentation

G4ParticleHPInelastic::G4ParticleHPInelastic ( G4ParticleDefinition *  projectile = G4Neutron::Neutron(),
const char *  name = "NeutronHPInelastic" 
)

Definition at line 47 of file G4ParticleHPInelastic.cc.

48  :G4HadronicInteraction(name)
49  ,theInelastic(NULL)
50  ,numEle(0)
51  ,theProjectile(projectile)
52 {
53  G4String baseName;
54  if ( getenv("G4PARTICLEHPDATA") ) {
55  baseName = getenv( "G4PARTICLEHPDATA" );
56  }
57  //const char* dataDirVariable;
58  G4String particleName;
59  if( theProjectile == G4Neutron::Neutron() ) {
60  dataDirVariable = "G4NEUTRONHPDATA";
61  }else if( theProjectile == G4Proton::Proton() ) {
62  dataDirVariable = "G4PROTONHPDATA";
63  particleName = "Proton";
64  }else if( theProjectile == G4Deuteron::Deuteron() ) {
65  dataDirVariable = "G4DEUTERONHPDATA";
66  particleName = "Deuteron";
67  }else if( theProjectile == G4Triton::Triton() ) {
68  dataDirVariable = "G4TRITONHPDATA";
69  particleName = "Triton";
70  }else if( theProjectile == G4He3::He3() ) {
71  dataDirVariable = "G4HE3HPDATA";
72  particleName = "He3";
73  }else if( theProjectile == G4Alpha::Alpha() ) {
74  dataDirVariable = "G4ALPHAHPDATA";
75  particleName = "Alpha";
76  } else {
77  G4String message("G4ParticleHPInelastic may only be called for neutron, proton, deuteron, triton, He3 or alpha, while it is called for " + theProjectile->GetParticleName());
78  throw G4HadronicException(__FILE__, __LINE__,message.c_str());
79  }
80 
81  SetMinEnergy( 0.0 );
82  SetMaxEnergy( 20.*MeV );
83 
84 // G4cout << " entering G4ParticleHPInelastic constructor"<<G4endl;
85  if ( !getenv("G4PARTICLEHPDATA") && !getenv(dataDirVariable) ) {
86  G4String message( "Please set the environement variable " + G4String(dataDirVariable) + " to point to the " + theProjectile->GetParticleName() + " cross-section files." );
87  throw G4HadronicException(__FILE__, __LINE__,message.c_str());
88  }
89  if ( getenv(dataDirVariable) ) {
90  dirName = getenv(dataDirVariable);
91  } else {
92  dirName = baseName + "/" + particleName;
93  }
94 G4cout << dirName << G4endl;
95 
96  G4String tString = "/Inelastic";
97  dirName = dirName + tString;
98  //numEle = G4Element::GetNumberOfElements();
99 
100  G4cout << "@@@ G4ParticleHPInelastic instantiated for particle " << theProjectile->GetParticleName() << " data directory variable is " << dataDirVariable << " pointing to " << dirName << G4endl;
101 
102 /*
103  theInelastic = new G4ParticleHPChannelList[numEle];
104  for (G4int i=0; i<numEle; i++)
105  {
106  theInelastic[i].Init((*(G4Element::GetElementTable()))[i], dirName);
107  G4int itry = 0;
108  do
109  {
110  theInelastic[i].Register(&theNFS, "F01"); // has
111  theInelastic[i].Register(&theNXFS, "F02");
112  theInelastic[i].Register(&the2NDFS, "F03");
113  theInelastic[i].Register(&the2NFS, "F04"); // has, E Done
114  theInelastic[i].Register(&the3NFS, "F05"); // has, E Done
115  theInelastic[i].Register(&theNAFS, "F06");
116  theInelastic[i].Register(&theN3AFS, "F07");
117  theInelastic[i].Register(&the2NAFS, "F08");
118  theInelastic[i].Register(&the3NAFS, "F09");
119  theInelastic[i].Register(&theNPFS, "F10");
120  theInelastic[i].Register(&theN2AFS, "F11");
121  theInelastic[i].Register(&the2N2AFS, "F12");
122  theInelastic[i].Register(&theNDFS, "F13");
123  theInelastic[i].Register(&theNTFS, "F14");
124  theInelastic[i].Register(&theNHe3FS, "F15");
125  theInelastic[i].Register(&theND2AFS, "F16");
126  theInelastic[i].Register(&theNT2AFS, "F17");
127  theInelastic[i].Register(&the4NFS, "F18"); // has, E Done
128  theInelastic[i].Register(&the2NPFS, "F19");
129  theInelastic[i].Register(&the3NPFS, "F20");
130  theInelastic[i].Register(&theN2PFS, "F21");
131  theInelastic[i].Register(&theNPAFS, "F22");
132  theInelastic[i].Register(&thePFS, "F23");
133  theInelastic[i].Register(&theDFS, "F24");
134  theInelastic[i].Register(&theTFS, "F25");
135  theInelastic[i].Register(&theHe3FS, "F26");
136  theInelastic[i].Register(&theAFS, "F27");
137  theInelastic[i].Register(&the2AFS, "F28");
138  theInelastic[i].Register(&the3AFS, "F29");
139  theInelastic[i].Register(&the2PFS, "F30");
140  theInelastic[i].Register(&thePAFS, "F31");
141  theInelastic[i].Register(&theD2AFS, "F32");
142  theInelastic[i].Register(&theT2AFS, "F33");
143  theInelastic[i].Register(&thePDFS, "F34");
144  theInelastic[i].Register(&thePTFS, "F35");
145  theInelastic[i].Register(&theDAFS, "F36");
146  theInelastic[i].RestartRegistration();
147  itry++;
148  }
149  //while(!theInelastic[i].HasDataInAnyFinalState());
150  while( !theInelastic[i].HasDataInAnyFinalState() && itry < 6 );
151  // 6 is corresponding to the value(5) of G4ParticleHPChannel. TK
152 
153  if ( itry == 6 )
154  {
155  // No Final State at all.
156  G4bool exceptional = false;
157  if ( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() == 1 )
158  {
159  if ( (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() == 1 && (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() == 1 ) exceptional = true; //1H
160  }
161  if ( !exceptional ) throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this element");
162  }
163  }
164 */
165 /*
166  for (G4int i=0; i<numEle; i++)
167  {
168  theInelastic.push_back( new G4ParticleHPChannelList );
169  (*theInelastic[i]).Init((*(G4Element::GetElementTable()))[i], dirName, theProjectile);
170  G4int itry = 0;
171  do
172  {
173  (*theInelastic[i]).Register(&theNFS, "F01"); // has
174  (*theInelastic[i]).Register(&theNXFS, "F02");
175  (*theInelastic[i]).Register(&the2NDFS, "F03");
176  (*theInelastic[i]).Register(&the2NFS, "F04"); // has, E Done
177  (*theInelastic[i]).Register(&the3NFS, "F05"); // has, E Done
178  (*theInelastic[i]).Register(&theNAFS, "F06");
179  (*theInelastic[i]).Register(&theN3AFS, "F07");
180  (*theInelastic[i]).Register(&the2NAFS, "F08");
181  (*theInelastic[i]).Register(&the3NAFS, "F09");
182  (*theInelastic[i]).Register(&theNPFS, "F10");
183  (*theInelastic[i]).Register(&theN2AFS, "F11");
184  (*theInelastic[i]).Register(&the2N2AFS, "F12");
185  (*theInelastic[i]).Register(&theNDFS, "F13");
186  (*theInelastic[i]).Register(&theNTFS, "F14");
187  (*theInelastic[i]).Register(&theNHe3FS, "F15");
188  (*theInelastic[i]).Register(&theND2AFS, "F16");
189  (*theInelastic[i]).Register(&theNT2AFS, "F17");
190  (*theInelastic[i]).Register(&the4NFS, "F18"); // has, E Done
191  (*theInelastic[i]).Register(&the2NPFS, "F19");
192  (*theInelastic[i]).Register(&the3NPFS, "F20");
193  (*theInelastic[i]).Register(&theN2PFS, "F21");
194  (*theInelastic[i]).Register(&theNPAFS, "F22");
195  (*theInelastic[i]).Register(&thePFS, "F23");
196  (*theInelastic[i]).Register(&theDFS, "F24");
197  (*theInelastic[i]).Register(&theTFS, "F25");
198  (*theInelastic[i]).Register(&theHe3FS, "F26");
199  (*theInelastic[i]).Register(&theAFS, "F27");
200  (*theInelastic[i]).Register(&the2AFS, "F28");
201  (*theInelastic[i]).Register(&the3AFS, "F29");
202  (*theInelastic[i]).Register(&the2PFS, "F30");
203  (*theInelastic[i]).Register(&thePAFS, "F31");
204  (*theInelastic[i]).Register(&theD2AFS, "F32");
205  (*theInelastic[i]).Register(&theT2AFS, "F33");
206  (*theInelastic[i]).Register(&thePDFS, "F34");
207  (*theInelastic[i]).Register(&thePTFS, "F35");
208  (*theInelastic[i]).Register(&theDAFS, "F36");
209  (*theInelastic[i]).RestartRegistration();
210  itry++;
211  }
212  while( !(*theInelastic[i]).HasDataInAnyFinalState() && itry < 6 );
213  // 6 is corresponding to the value(5) of G4ParticleHPChannel. TK
214 
215  if ( itry == 6 )
216  {
217  // No Final State at all.
218  G4bool exceptional = false;
219  if ( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() == 1 )
220  {
221  if ( (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() == 1 && (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() == 1 ) exceptional = true; //1H
222  }
223  if ( !exceptional ) {
224  G4cerr << " ELEMENT Z " << (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() << " N " << (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() << G4endl; //1H
225 throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this element");
226  }
227  }
228 
229  }
230 */
231 
232  }
name
const XML_Char * name
Definition: expat.h:151
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:90
G4cout
G4GLOB_DLL std::ostream G4cout
G4ParticleHPInelastic::theInelastic
std::vector< G4ParticleHPChannelList * > * theInelastic
Definition: G4ParticleHPInelastic.hh:109
G4Triton::Triton
static G4Triton * Triton()
Definition: G4Triton.cc:95
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93
G4HadronicInteraction::G4HadronicInteraction
G4HadronicInteraction(const G4String &modelName="HadronicModel")
Definition: G4HadronicInteraction.cc:42
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Deuteron::Deuteron
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:94
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:103
G4endl
#define G4endl
Definition: G4ios.hh:61
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4Alpha::Alpha
static G4Alpha * Alpha()
Definition: G4Alpha.cc:89
G4He3::He3
static G4He3 * He3()
Definition: G4He3.cc:94

Here is the call graph for this function:

G4ParticleHPInelastic::~G4ParticleHPInelastic ( )

Definition at line 234 of file G4ParticleHPInelastic.cc.

235  {
236 // delete [] theInelastic;
237  //Vector is shared, only master deletes
238  if ( !G4Threading::IsWorkerThread() ) {
239  if ( theInelastic != NULL ) {
240  for ( std::vector<G4ParticleHPChannelList*>::iterator
241  it = theInelastic->begin() ; it != theInelastic->end() ; it++ ) {
242  delete *it;
243  }
244  theInelastic->clear();
245  }
246  }
247  }
G4ParticleHPInelastic::theInelastic
std::vector< G4ParticleHPChannelList * > * theInelastic
Definition: G4ParticleHPInelastic.hh:109
G4Threading::IsWorkerThread
G4bool IsWorkerThread()
Definition: G4Threading.cc:145

Here is the call graph for this function:

Member Function Documentation

G4HadFinalState * G4ParticleHPInelastic::ApplyYourself ( const G4HadProjectile &  aTrack,
G4Nucleus &  aTargetNucleus 
)
virtual

Implements G4HadronicInteraction.

Definition at line 251 of file G4ParticleHPInelastic.cc.

252  {
253  G4ParticleHPManager::GetInstance()->OpenReactionWhiteBoard();
254  const G4Material * theMaterial = aTrack.GetMaterial();
255  G4int n = theMaterial->GetNumberOfElements();
256  G4int index = theMaterial->GetElement(0)->GetIndex();
257  G4int it=0;
258  if(n!=1)
259  {
260  G4double* xSec = new G4double[n];
261  G4double sum=0;
262  G4int i;
263  const G4double * NumAtomsPerVolume = theMaterial->GetVecNbOfAtomsPerVolume();
264  G4double rWeight;
265  G4ParticleHPThermalBoost aThermalE;
266  for (i=0; i<n; i++)
267  {
268  index = theMaterial->GetElement(i)->GetIndex();
269  rWeight = NumAtomsPerVolume[i];
270  if ( aTrack.GetDefinition() == G4Neutron::Neutron() ) {
271  xSec[i] = ((*theInelastic)[index])->GetXsec(aThermalE.GetThermalEnergy(aTrack,
272  theMaterial->GetElement(i),
273  theMaterial->GetTemperature()));
274  } else {
275  xSec[i] = ((*theInelastic)[index])->GetXsec(aTrack.GetKineticEnergy());
276  }
277  xSec[i] *= rWeight;
278  sum+=xSec[i];
279 #ifdef G4PHPDEBUG
280  if( getenv("G4ParticleHPDebug") ) G4cout << " G4ParticleHPInelastic XSEC ELEM " << i << " = " << xSec[i] << G4endl;
281 #endif
282 
283  }
284  G4double random = G4UniformRand();
285  G4double running = 0;
286  for (i=0; i<n; i++)
287  {
288  running += xSec[i];
289  index = theMaterial->GetElement(i)->GetIndex();
290  it = i;
291  //if(random<=running/sum) break;
292  if( sum == 0 || random<=running/sum) break;
293  }
294  delete [] xSec;
295  }
296 
297 #ifdef G4PHPDEBUG
298  if( getenv("G4ParticleHPDebug") ) G4cout << " G4ParticleHPInelastic SELECTED ELEM " << it << " = " << theMaterial->GetElement(it)->GetName() << " FROM MATERIAL " << theMaterial->GetName() << G4endl;
299 #endif
300  //return theInelastic[index].ApplyYourself(theMaterial->GetElement(it), aTrack);
301  G4HadFinalState* result = ((*theInelastic)[index])->ApplyYourself(theMaterial->GetElement(it), aTrack);
302  //
303  aNucleus.SetParameters(G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA(),G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargZ());
304  const G4Element* target_element = (*G4Element::GetElementTable())[index];
305  const G4Isotope* target_isotope=NULL;
306  G4int iele = target_element->GetNumberOfIsotopes();
307  for ( G4int j = 0 ; j != iele ; j++ ) {
308  target_isotope=target_element->GetIsotope( j );
309  if ( target_isotope->GetN() == G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA() ) break;
310  }
311  //G4cout << "Target Material of this reaction is " << theMaterial->GetName() << G4endl;
312  //G4cout << "Target Element of this reaction is " << target_element->GetName() << G4endl;
313  //G4cout << "Target Isotope of this reaction is " << target_isotope->GetName() << G4endl;
314  aNucleus.SetIsotope( target_isotope );
315 
316  G4ParticleHPManager::GetInstance()->CloseReactionWhiteBoard();
317 
318  //GDEB
319  if( getenv("G4PHPTEST") ) {
320  G4HadSecondary* seco = result->GetSecondary(0);
321  if(seco) {
322  G4ThreeVector secoMom = seco->GetParticle()->GetMomentum();
323  G4cout << " G4ParticleHPinelastic COS THETA " << std::cos(secoMom.theta()) <<" " << secoMom << G4endl;
324  }
325  }
326 
327  return result;
328  }
result
G4double G4ParticleHPJENDLHEData::G4double result
Definition: G4ParticleHPJENDLHEData.cc:257
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:56
G4Element::GetNumberOfIsotopes
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:159
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:79
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:178
G4Material::GetElement
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:202
G4int
int G4int
Definition: G4Types.hh:78
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:206
G4Isotope::GetN
G4int GetN() const
Definition: G4Isotope.hh:94
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:97
G4cout
G4GLOB_DLL std::ostream G4cout
G4HadProjectile::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4HadProjectile.hh:81
G4HadProjectile::GetKineticEnergy
G4double GetKineticEnergy() const
Definition: G4HadProjectile.hh:106
G4Element::GetIndex
size_t GetIndex() const
Definition: G4Element.hh:182
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
n
const G4int n
Definition: G4UrQMD1_3Interface.hh:144
CLHEP::Hep3Vector::theta
double theta() const
G4ParticleHPThermalBoost::GetThermalEnergy
G4double GetThermalEnergy(const G4HadProjectile &aP, const G4Element *anE, G4double aT)
Definition: G4ParticleHPThermalBoost.hh:45
G4HadSecondary::GetParticle
G4DynamicParticle * GetParticle()
Definition: G4HadSecondary.hh:42
G4Element::GetIsotope
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:170
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
G4double
double G4double
Definition: G4Types.hh:76
G4Element::GetName
const G4String & GetName() const
Definition: G4Element.hh:127
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:398
G4ParticleHPManager::GetReactionWhiteBoard
G4ParticleHPReactionWhiteBoard * GetReactionWhiteBoard()
Definition: G4ParticleHPManager.cc:84
G4ParticleHPInelastic::ApplyYourself
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
Definition: G4ParticleHPInelastic.cc:251
G4DynamicParticle::GetMomentum
G4ThreeVector GetMomentum() const

Here is the call graph for this function:

void G4ParticleHPInelastic::BuildPhysicsTable ( const G4ParticleDefinition &  projectile)
virtual

Reimplemented from G4HadronicInteraction.

Definition at line 454 of file G4ParticleHPInelastic.cc.

454  {
455 
456  G4ParticleHPManager* hpmanager = G4ParticleHPManager::GetInstance();
457 
458  theInelastic = hpmanager->GetInelasticFinalStates( &projectile );
459 
460  if ( G4Threading::IsMasterThread() ) {
461 
462  if ( theInelastic == NULL ) theInelastic = new std::vector<G4ParticleHPChannelList*>;
463 
464  if ( numEle == (G4int)G4Element::GetNumberOfElements() ) return;
465 
466  if ( theInelastic->size() == G4Element::GetNumberOfElements() ) {
467  numEle = G4Element::GetNumberOfElements();
468  return;
469  }
470 /*
471  const char* dataDirVariable;
472  if( &projectile == G4Neutron::Neutron() ) {
473  dataDirVariable = "G4NEUTRONHPDATA";
474  }else if( &projectile == G4Proton::Proton() ) {
475  dataDirVariable = "G4PROTONHPDATA";
476  }else if( &projectile == G4Deuteron::Deuteron() ) {
477  dataDirVariable = "G4DEUTERONHPDATA";
478  }else if( &projectile == G4Triton::Triton() ) {
479  dataDirVariable = "G4TRITONHPDATA";
480  }else if( &projectile == G4He3::He3() ) {
481  dataDirVariable = "G4HE3HPDATA";
482  }else if( &projectile == G4Alpha::Alpha() ) {
483  dataDirVariable = "G4ALPHAHPDATA";
484  } else {
485  G4String message("G4ParticleHPInelastic may only be called for neutron, proton, deuteron, triton, He3 or alpha, while it is called for " + projectile.GetParticleName());
486  throw G4HadronicException(__FILE__, __LINE__,message.c_str());
487  }
488  if(!getenv(dataDirVariable)){
489  G4String message("Please set the environement variable " + G4String(dataDirVariable) + " to point to the " + projectile.GetParticleName() + " cross-section files.");
490  throw G4HadronicException(__FILE__, __LINE__,message.c_str());
491  }
492  dirName = getenv(dataDirVariable);
493  G4cout << dirName << G4endl;
494 
495  G4String tString = "/Inelastic";
496  dirName = dirName + tString;
497 
498 */
499  G4cout << "@@@ G4ParticleHPInelastic instantiated for particle " << projectile.GetParticleName() << " data directory variable is " << dataDirVariable << " pointing to " << dirName << G4endl;
500  for (G4int i = numEle ; i < (G4int)G4Element::GetNumberOfElements(); i++)
501  {
502  theInelastic->push_back( new G4ParticleHPChannelList );
503  ((*theInelastic)[i])->Init((*(G4Element::GetElementTable()))[i], dirName, const_cast<G4ParticleDefinition*>(&projectile));
504  G4int itry = 0;
505  do
506  {
507  ((*theInelastic)[i])->Register( new G4ParticleHPNInelasticFS , "F01"); // has
508  ((*theInelastic)[i])->Register( new G4ParticleHPNXInelasticFS , "F02");
509  ((*theInelastic)[i])->Register( new G4ParticleHP2NDInelasticFS , "F03");
510  ((*theInelastic)[i])->Register( new G4ParticleHP2NInelasticFS , "F04"); // has, E Done
511  ((*theInelastic)[i])->Register( new G4ParticleHP3NInelasticFS , "F05"); // has, E Done
512  ((*theInelastic)[i])->Register( new G4ParticleHPNAInelasticFS , "F06");
513  ((*theInelastic)[i])->Register( new G4ParticleHPN3AInelasticFS , "F07");
514  ((*theInelastic)[i])->Register( new G4ParticleHP2NAInelasticFS , "F08");
515  ((*theInelastic)[i])->Register( new G4ParticleHP3NAInelasticFS , "F09");
516  ((*theInelastic)[i])->Register( new G4ParticleHPNPInelasticFS , "F10");
517  ((*theInelastic)[i])->Register( new G4ParticleHPN2AInelasticFS , "F11");
518  ((*theInelastic)[i])->Register( new G4ParticleHP2N2AInelasticFS , "F12");
519  ((*theInelastic)[i])->Register( new G4ParticleHPNDInelasticFS , "F13");
520  ((*theInelastic)[i])->Register( new G4ParticleHPNTInelasticFS , "F14");
521  ((*theInelastic)[i])->Register( new G4ParticleHPNHe3InelasticFS , "F15");
522  ((*theInelastic)[i])->Register( new G4ParticleHPND2AInelasticFS , "F16");
523  ((*theInelastic)[i])->Register( new G4ParticleHPNT2AInelasticFS , "F17");
524  ((*theInelastic)[i])->Register( new G4ParticleHP4NInelasticFS , "F18"); // has, E Done
525  ((*theInelastic)[i])->Register( new G4ParticleHP2NPInelasticFS , "F19");
526  ((*theInelastic)[i])->Register( new G4ParticleHP3NPInelasticFS , "F20");
527  ((*theInelastic)[i])->Register( new G4ParticleHPN2PInelasticFS , "F21");
528  ((*theInelastic)[i])->Register( new G4ParticleHPNPAInelasticFS , "F22");
529  ((*theInelastic)[i])->Register( new G4ParticleHPPInelasticFS , "F23");
530  ((*theInelastic)[i])->Register( new G4ParticleHPDInelasticFS , "F24");
531  ((*theInelastic)[i])->Register( new G4ParticleHPTInelasticFS , "F25");
532  ((*theInelastic)[i])->Register( new G4ParticleHPHe3InelasticFS , "F26");
533  ((*theInelastic)[i])->Register( new G4ParticleHPAInelasticFS , "F27");
534  ((*theInelastic)[i])->Register( new G4ParticleHP2AInelasticFS , "F28");
535  ((*theInelastic)[i])->Register( new G4ParticleHP3AInelasticFS , "F29");
536  ((*theInelastic)[i])->Register( new G4ParticleHP2PInelasticFS , "F30");
537  ((*theInelastic)[i])->Register( new G4ParticleHPPAInelasticFS , "F31");
538  ((*theInelastic)[i])->Register( new G4ParticleHPD2AInelasticFS , "F32");
539  ((*theInelastic)[i])->Register( new G4ParticleHPT2AInelasticFS , "F33");
540  ((*theInelastic)[i])->Register( new G4ParticleHPPDInelasticFS , "F34");
541  ((*theInelastic)[i])->Register( new G4ParticleHPPTInelasticFS , "F35");
542  ((*theInelastic)[i])->Register( new G4ParticleHPDAInelasticFS , "F36");
543  ((*theInelastic)[i])->RestartRegistration();
544  itry++;
545  }
546  while( !((*theInelastic)[i])->HasDataInAnyFinalState() && itry < 6 ); // Loop checking, 11.05.2015, T. Koi
547  // 6 is corresponding to the value(5) of G4ParticleHPChannel. TK
548 
549  if ( itry == 6 ) {
550  // No Final State at all.
551  /*
552  G4bool exceptional = false;
553  if ( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() == 1 )
554  {
555  if ( (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() == 1 && (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() == 1 ) exceptional = true; //1H
556  }
557  if ( !exceptional ) {
558  G4cerr << " ELEMENT Z " << (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetZ() << " N " << (*(G4Element::GetElementTable()))[i]->GetIsotope( 0 )->GetN() << G4endl; //1H
559  throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this element");
560  }
561  */
562  if ( G4ParticleHPManager::GetInstance()->GetVerboseLevel() > 1 ) {
563  G4cout << "ParticleHP::Inelastic for " << projectile.GetParticleName() << ". Do not know what to do with element of \"" << (*(G4Element::GetElementTable()))[i]->GetName() << "\"." << G4endl;
564  G4cout << "The components of the element are" << G4endl;
565  //G4cout << "TKDB dataDirVariable = " << dataDirVariable << G4endl;
566  for ( G4int ii = 0 ; ii < (G4int)( (*(G4Element::GetElementTable()))[i]->GetNumberOfIsotopes() ) ; ii++ ) {
567  G4cout << " Z: " << (*(G4Element::GetElementTable()))[i]->GetIsotope( ii )->GetZ() << ", A: " << (*(G4Element::GetElementTable()))[i]->GetIsotope( ii )->GetN() << G4endl;
568  }
569  G4cout << "No possible final state data of the element is found in " << dataDirVariable << "." << G4endl;
570  }
571  }
572  }
573  hpmanager->RegisterInelasticFinalStates( &projectile , theInelastic );
574  }
575  numEle = G4Element::GetNumberOfElements();
576 }
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:56
lightions::Init
void Init()
Definition: G4IonTable.cc:90
G4int
int G4int
Definition: G4Types.hh:78
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4AutoDelete::Register
void Register(T *inst)
Definition: G4AutoDelete.hh:65
G4cout
G4GLOB_DLL std::ostream G4cout
G4Element::GetNumberOfElements
static size_t GetNumberOfElements()
Definition: G4Element.cc:405
G4ParticleHPInelastic::theInelastic
std::vector< G4ParticleHPChannelList * > * theInelastic
Definition: G4ParticleHPInelastic.hh:109
G4ParticleHPManager::GetInelasticFinalStates
std::vector< G4ParticleHPChannelList * > * GetInelasticFinalStates(const G4ParticleDefinition *)
Definition: G4ParticleHPManager.cc:224
G4endl
#define G4endl
Definition: G4ios.hh:61
G4Threading::IsMasterThread
G4bool IsMasterThread()
Definition: G4Threading.cc:146
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:398
G4ParticleHPManager::RegisterInelasticFinalStates
void RegisterInelasticFinalStates(const G4ParticleDefinition *, std::vector< G4ParticleHPChannelList * > *)
Definition: G4ParticleHPManager.cc:231

Here is the call graph for this function:

const std::pair< G4double, G4double > G4ParticleHPInelastic::GetFatalEnergyCheckLevels ( ) const
virtual

Reimplemented from G4HadronicInteraction.

Definition at line 330 of file G4ParticleHPInelastic.cc.

331 {
332  // max energy non-conservation is mass of heavy nucleus
333 // if ( getenv("G4PHP_DO_NOT_ADJUST_FINAL_STATE") ) return std::pair<G4double, G4double>(5*perCent,250*GeV);
334  // This should be same to the hadron default value
335 // return std::pair<G4double, G4double>(10*perCent,10*GeV);
336  return std::pair<G4double, G4double>(10*perCent,DBL_MAX);
337 }
perCent
static constexpr double perCent
Definition: G4SIunits.hh:332
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4int G4ParticleHPInelastic::GetVerboseLevel ( ) const

Definition at line 408 of file G4ParticleHPInelastic.cc.

409 {
410  return G4ParticleHPManager::GetInstance()->GetVerboseLevel();
411 }
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:56

Here is the call graph for this function:

void G4ParticleHPInelastic::ModelDescription ( std::ostream &  outFile) const
virtual

Reimplemented from G4HadronicInteraction.

Definition at line 578 of file G4ParticleHPInelastic.cc.

579 {
580  outFile << "Extension of High Precision model for inelastic reaction of proton, deuteron, triton, He3 and alpha below 20MeV\n";
581 }
void G4ParticleHPInelastic::SetVerboseLevel ( G4int  newValue)

Definition at line 412 of file G4ParticleHPInelastic.cc.

413 {
414  G4ParticleHPManager::GetInstance()->SetVerboseLevel(newValue);
415 }
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:56

Here is the call graph for this function:

Member Data Documentation

G4String G4ParticleHPInelastic::dataDirVariable
protected

Definition at line 110 of file G4ParticleHPInelastic.hh.

G4String G4ParticleHPInelastic::dirName
protected

Definition at line 111 of file G4ParticleHPInelastic.hh.

G4int G4ParticleHPInelastic::numEle
protected

Definition at line 112 of file G4ParticleHPInelastic.hh.

std::vector<G4ParticleHPChannelList*>* G4ParticleHPInelastic::theInelastic
protected

Definition at line 109 of file G4ParticleHPInelastic.hh.

The documentation for this class was generated from the following files: