Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4NeutronHPorLEInelastic Class Reference

#include <G4NeutronHPorLEInelastic.hh>

+ Inheritance diagram for G4NeutronHPorLEInelastic:

Public Member Functions

 G4NeutronHPorLEInelastic ()
 
 ~G4NeutronHPorLEInelastic ()
 
G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels () const
 
G4bool IsThisElementOK (G4String)
 
G4VCrossSectionDataSet * GiveXSectionDataSet ()
 
- 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
 
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
 
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
 

Additional Inherited Members

- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 95 of file G4NeutronHPorLEInelastic.hh.

Constructor & Destructor Documentation

◆ G4NeutronHPorLEInelastic()

G4NeutronHPorLEInelastic::G4NeutronHPorLEInelastic ( )

Definition at line 41 of file G4NeutronHPorLEInelastic.cc.

42 :G4HadronicInteraction("NeutronHPorLEInelastic")
43{
44 SetMinEnergy(0.*eV);
45 SetMaxEnergy(20.*MeV);
46
47// G4NeutronHPInelasticFS * theFS = new G4NeutronHPInelasticFS;
48 if(!getenv("G4NEUTRONHPDATA"))
49 throw G4HadronicException(__FILE__, __LINE__, "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
50 dirName = getenv("G4NEUTRONHPDATA");
51 G4String tString = "/Inelastic/";
52 dirName = dirName + tString;
53// G4cout <<"G4NeutronHPorLEInelastic::G4NeutronHPorLEInelastic testit "<<dirName<<G4endl;
54 numEle = G4Element::GetNumberOfElements();
55 theInelastic = new G4NeutronHPChannelList[numEle];
56 unavailable_elements.clear();
57 for (G4int i=0; i<numEle; i++)
58 {
59 theInelastic[i].Init( (*(G4Element::GetElementTable()))[i] , dirName );
60 do
61 {
62
63 try
64 {
65 theInelastic[i].Register(&theNFS, "F01"); // has
66 theInelastic[i].Register(&theNXFS, "F02");
67 theInelastic[i].Register(&the2NDFS, "F03");
68 theInelastic[i].Register(&the2NFS, "F04"); // has, E Done
69 theInelastic[i].Register(&the3NFS, "F05"); // has, E Done
70 theInelastic[i].Register(&theNAFS, "F06");
71 theInelastic[i].Register(&theN3AFS, "F07");
72 theInelastic[i].Register(&the2NAFS, "F08");
73 theInelastic[i].Register(&the3NAFS, "F09");
74 theInelastic[i].Register(&theNPFS, "F10");
75 theInelastic[i].Register(&theN2AFS, "F11");
76 theInelastic[i].Register(&the2N2AFS, "F12");
77 theInelastic[i].Register(&theNDFS, "F13");
78 theInelastic[i].Register(&theNTFS, "F14");
79 theInelastic[i].Register(&theNHe3FS, "F15");
80 theInelastic[i].Register(&theND2AFS, "F16");
81 theInelastic[i].Register(&theNT2AFS, "F17");
82 theInelastic[i].Register(&the4NFS, "F18"); // has, E Done
83 theInelastic[i].Register(&the2NPFS, "F19");
84 theInelastic[i].Register(&the3NPFS, "F20");
85 theInelastic[i].Register(&theN2PFS, "F21");
86 theInelastic[i].Register(&theNPAFS, "F22");
87 theInelastic[i].Register(&thePFS, "F23");
88 theInelastic[i].Register(&theDFS, "F24");
89 theInelastic[i].Register(&theTFS, "F25");
90 theInelastic[i].Register(&theHe3FS, "F26");
91 theInelastic[i].Register(&theAFS, "F27");
92 theInelastic[i].Register(&the2AFS, "F28");
93 theInelastic[i].Register(&the3AFS, "F29");
94 theInelastic[i].Register(&the2PFS, "F30");
95 theInelastic[i].Register(&thePAFS, "F31");
96 theInelastic[i].Register(&theD2AFS, "F32");
97 theInelastic[i].Register(&theT2AFS, "F33");
98 theInelastic[i].Register(&thePDFS, "F34");
99 theInelastic[i].Register(&thePTFS, "F35");
100 theInelastic[i].Register(&theDAFS, "F36");
101 }
102
103 catch ( G4HadronicException )
104 {
105 unavailable_elements.insert ( (*(G4Element::GetElementTable()))[i]->GetName() );
106 }
107 theInelastic[i].RestartRegistration();
108 }
109 while( !theInelastic[i].HasDataInAnyFinalState());
110
111 }
112
113// delete theFS;
114 if ( unavailable_elements.size() > 0 )
115 {
116 std::set< G4String>::iterator it;
117 G4cout << "HP Inelastic data are not available for thess elements "<< G4endl;
118 for ( it = unavailable_elements.begin() ; it != unavailable_elements.end() ; it++ )
119 G4cout << *it << G4endl;
120 G4cout << "Low Energy Parameterization Models will be used."<< G4endl;
121 }
122
123 createXSectionDataSet();
124}
G4int
int G4int
Definition: G4Types.hh:66
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
G4Element::GetNumberOfElements
static size_t GetNumberOfElements()
Definition: G4Element.cc:406
G4Element::GetElementTable
static const G4ElementTable * GetElementTable()
Definition: G4Element.cc:399
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:90
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:103
G4NeutronHPChannelList::Register
void Register(G4NeutronHPFinalState *theFS, const G4String &aName)
Definition: G4NeutronHPChannelList.cc:174

◆ ~G4NeutronHPorLEInelastic()

G4NeutronHPorLEInelastic::~G4NeutronHPorLEInelastic ( )

Definition at line 126 of file G4NeutronHPorLEInelastic.cc.

127{
128 delete [] theInelastic;
129 delete theDataSet;
130}

Member Function Documentation

◆ ApplyYourself()

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

Implements G4HadronicInteraction.

Definition at line 134 of file G4NeutronHPorLEInelastic.cc.

135{
136 G4int it=0;
137 const G4Material * theMaterial = aTrack.GetMaterial();
138 G4int n = theMaterial->GetNumberOfElements();
139 G4int index = theMaterial->GetElement(0)->GetIndex();
140 if(n!=1)
141 {
142 G4int i;
143 xSec = new G4double[n];
144 G4double sum=0;
145 const G4double * NumAtomsPerVolume = theMaterial->GetVecNbOfAtomsPerVolume();
146 G4double rWeight;
147 G4NeutronHPThermalBoost aThermalE;
148 for (i=0; i<n; i++)
149 {
150 index = theMaterial->GetElement(i)->GetIndex();
151 rWeight = NumAtomsPerVolume[i];
152 G4double x = aThermalE.GetThermalEnergy(aTrack, theMaterial->GetElement(i), theMaterial->GetTemperature());
153
154 //xSec[i] = theInelastic[index].GetXsec(aThermalE.GetThermalEnergy(aTrack,
155 // theMaterial->GetElement(i),
156 // theMaterial->GetTemperature()));
157 xSec[i] = theInelastic[index].GetXsec(x);
158
159 xSec[i] *= rWeight;
160 sum+=xSec[i];
161 }
162 G4double random = G4UniformRand();
163 G4double running = 0;
164 for (i=0; i<n; i++)
165 {
166 running += xSec[i];
167 index = theMaterial->GetElement(i)->GetIndex();
168 it = i;
169 if(random<=running/sum) break;
170 }
171 delete [] xSec;
172 // it is element-wise initialised.
173 }
174 //return theInelastic[index].ApplyYourself(aTrack);
175 return theInelastic[index].ApplyYourself( theMaterial->GetElement(it) , aTrack );
176}
G4double
double G4double
Definition: G4Types.hh:64
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53
G4Element::GetIndex
size_t GetIndex() const
Definition: G4Element.hh:182
G4HadProjectile::GetMaterial
const G4Material * GetMaterial() const
Definition: G4HadProjectile.hh:76
G4Material::GetTemperature
G4double GetTemperature() const
Definition: G4Material.hh:181
G4Material::GetElement
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:201
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:185
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:205
G4NeutronHPChannelList::GetXsec
G4double GetXsec(G4double anEnergy)
Definition: G4NeutronHPChannelList.hh:64
G4NeutronHPChannelList::ApplyYourself
G4HadFinalState * ApplyYourself(const G4Element *theElement, const G4HadProjectile &aTrack)
Definition: G4NeutronHPChannelList.cc:70
G4NeutronHPThermalBoost::GetThermalEnergy
G4double GetThermalEnergy(const G4HadProjectile &aP, const G4Element *anE, G4double aT)
Definition: G4NeutronHPThermalBoost.hh:42

Referenced by G4NeutronHPorLEInelasticModel::ApplyYourself().

◆ GetFatalEnergyCheckLevels()

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

Reimplemented from G4HadronicInteraction.

Definition at line 194 of file G4NeutronHPorLEInelastic.cc.

195{
196 //return std::pair<G4double, G4double>(10*perCent,10*GeV);
197 return std::pair<G4double, G4double>(10*perCent,DBL_MAX);
198}
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83

◆ GiveXSectionDataSet()

G4VCrossSectionDataSet * G4NeutronHPorLEInelastic::GiveXSectionDataSet ( )
inline

Definition at line 158 of file G4NeutronHPorLEInelastic.hh.

158{ return theDataSet; };

Referenced by G4NeutronHPorLEInelasticModel::GiveHPXSectionDataSet().

◆ IsThisElementOK()

G4bool G4NeutronHPorLEInelastic::IsThisElementOK ( G4String  name)

Definition at line 180 of file G4NeutronHPorLEInelastic.cc.

181{
182 if ( unavailable_elements.find( name ) == unavailable_elements.end() )
183 return TRUE;
184 else
185 return FALSE;
186}
TRUE
#define TRUE
Definition: globals.hh:55
FALSE
#define FALSE
Definition: globals.hh:52

Referenced by G4NeutronHPorLEInelasticModel::ApplyYourself().

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