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

#include <G4ParticleHPElastic.hh>

+ Inheritance diagram for G4ParticleHPElastic:

Public Member Functions

 G4ParticleHPElastic ()
 
 ~G4ParticleHPElastic ()
 
G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus, G4bool isFromTSL)
 
virtual const std::pair< G4double, G4double > GetFatalEnergyCheckLevels () const
 
G4int GetNiso ()
 
void DoNotSuspend ()
 
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 G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
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 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
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void InitialiseModel ()
 
 G4HadronicInteraction (const G4HadronicInteraction &right)=delete
 
const G4HadronicInteraction & operator= (const G4HadronicInteraction &right)=delete
 
G4bool operator== (const G4HadronicInteraction &right) const =delete
 
G4bool operator!= (const G4HadronicInteraction &right) const =delete
 

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 48 of file G4ParticleHPElastic.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPElastic()

G4ParticleHPElastic::G4ParticleHPElastic ( )

Definition at line 43 of file G4ParticleHPElastic.cc.

44 : G4HadronicInteraction("NeutronHPElastic"), theElastic(nullptr), numEle(0)
45{
46 overrideSuspension = false;
47 SetMinEnergy(0.*eV);
48 SetMaxEnergy(20.*MeV);
49}
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition: G4HadronicInteraction.hh:89
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition: G4HadronicInteraction.hh:102

◆ ~G4ParticleHPElastic()

G4ParticleHPElastic::~G4ParticleHPElastic ( )

Definition at line 52 of file G4ParticleHPElastic.cc.

53{
54 //the vectror is shared among threads, only master deletes
55 if ( ! G4Threading::IsWorkerThread() ) {
56 if ( theElastic != nullptr ) {
57 for ( auto it=theElastic->cbegin(); it!=theElastic->cend(); ++it ) {
58 delete *it;
59 }
60 theElastic->clear();
61 }
62 }
63}
G4Threading::IsWorkerThread
G4bool IsWorkerThread()
Definition: G4Threading.cc:123

Member Function Documentation

◆ ApplyYourself() [1/2]

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

Reimplemented from G4HadronicInteraction.

Definition at line 66 of file G4ParticleHPElastic.cc.

67{
68 return this->ApplyYourself(aTrack, aNucleus, 0);
69}
G4ParticleHPElastic::ApplyYourself
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
Definition: G4ParticleHPElastic.cc:66

Referenced by ApplyYourself().

◆ ApplyYourself() [2/2]

G4HadFinalState * G4ParticleHPElastic::ApplyYourself ( const G4HadProjectile &  aTrack,
G4Nucleus &  aTargetNucleus,
G4bool  isFromTSL 
)

Definition at line 75 of file G4ParticleHPElastic.cc.

76{
77 G4ParticleHPManager::GetInstance()->OpenReactionWhiteBoard();
78 const G4Material * theMaterial = aTrack.GetMaterial();
79 G4int n = (G4int)theMaterial->GetNumberOfElements();
80 std::size_t index = theMaterial->GetElement(0)->GetIndex();
81
82 if ( ! isFromTSL ) {
83 if ( n != 1 ) {
84 G4int i;
85 G4double* xSec = new G4double[n];
86 G4double sum=0;
87 const G4double * NumAtomsPerVolume = theMaterial->GetVecNbOfAtomsPerVolume();
88 G4double rWeight;
89 G4ParticleHPThermalBoost aThermalE;
90 for ( i = 0; i < n; ++i ) {
91 index = theMaterial->GetElement(i)->GetIndex();
92 rWeight = NumAtomsPerVolume[i];
93 xSec[i] = ((*theElastic)[index])->GetXsec(aThermalE.GetThermalEnergy(aTrack,
94 theMaterial->GetElement(i),
95 theMaterial->GetTemperature()));
96 xSec[i] *= rWeight;
97 sum+=xSec[i];
98 }
99 G4double random = G4UniformRand();
100 G4double running = 0;
101 for ( i = 0; i < n; ++i ) {
102 running += xSec[i];
103 index = theMaterial->GetElement(i)->GetIndex();
104 if ( sum == 0 || random <= running/sum ) break;
105 }
106 delete [] xSec;
107 }
108 } else {
109 G4int i;
110 if ( n != 1 ) {
111 for ( i = 0; i < n; ++i ) {
112 if ( aNucleus.GetZ_asInt() == (G4int)(theMaterial->GetElement(i)->GetZ()) ) {
113 index = theMaterial->GetElement(i)->GetIndex();
114 }
115 }
116 }
117 }
118
119 G4HadFinalState* finalState = ((*theElastic)[index])->ApplyYourself(aTrack);
120 if (overrideSuspension) finalState->SetStatusChange(isAlive);
121
122 // Overwrite target parameters
123 aNucleus.SetParameters(G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA(),G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargZ());
124 const G4Element* target_element = (*G4Element::GetElementTable())[index];
125 const G4Isotope* target_isotope=nullptr;
126 G4int iele = (G4int)target_element->GetNumberOfIsotopes();
127 for ( G4int j = 0 ; j != iele ; ++j ) {
128 target_isotope=target_element->GetIsotope( j );
129 if ( target_isotope->GetN() == G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->GetTargA() ) break;
130 }
131 aNucleus.SetIsotope( target_isotope );
132
133 G4ParticleHPManager::GetInstance()->CloseReactionWhiteBoard();
134 return finalState;
135}
isAlive
@ isAlive
Definition: G4HadFinalState.hh:42
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:403
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4Element::GetIsotope
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:170
G4Element::GetIndex
size_t GetIndex() const
Definition: G4Element.hh:182
G4Element::GetNumberOfIsotopes
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:159
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.hh:67
G4HadProjectile::GetMaterial
const G4Material * GetMaterial() const
Definition: G4HadProjectile.hh:81
G4Isotope::GetN
G4int GetN() const
Definition: G4Isotope.hh:93
G4Material::GetTemperature
G4double GetTemperature() const
Definition: G4Material.hh:177
G4Material::GetElement
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:197
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:181
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:201
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.cc:73
G4ParticleHPManager::GetReactionWhiteBoard
G4ParticleHPReactionWhiteBoard * GetReactionWhiteBoard()
Definition: G4ParticleHPManager.cc:88
G4ParticleHPThermalBoost::GetThermalEnergy
G4double GetThermalEnergy(const G4HadProjectile &aP, const G4Element *anE, G4double aT)
Definition: G4ParticleHPThermalBoost.hh:45

◆ BuildPhysicsTable()

void G4ParticleHPElastic::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4HadronicInteraction.

Definition at line 157 of file G4ParticleHPElastic.cc.

158{
159
160 G4ParticleHPManager* hpmanager = G4ParticleHPManager::GetInstance();
161
162 theElastic = hpmanager->GetElasticFinalStates();
163
164 if ( G4Threading::IsMasterThread() ) {
165
166 if ( theElastic == nullptr ) theElastic = new std::vector<G4ParticleHPChannel*>;
167
168 if ( numEle == (G4int)G4Element::GetNumberOfElements() ) return;
169
170 if ( theElastic->size() == G4Element::GetNumberOfElements() ) {
171 numEle = (G4int)G4Element::GetNumberOfElements();
172 return;
173 }
174
175 G4ParticleHPElasticFS * theFS = new G4ParticleHPElasticFS;
176 if(!G4FindDataDir("G4NEUTRONHPDATA"))
177 throw G4HadronicException(__FILE__, __LINE__, "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
178 dirName = G4FindDataDir("G4NEUTRONHPDATA");
179 G4String tString = "/Elastic";
180 dirName = dirName + tString;
181 for ( G4int i = numEle; i < (G4int)G4Element::GetNumberOfElements(); ++i ) {
182 theElastic->push_back( new G4ParticleHPChannel );
183 ((*theElastic)[i])->Init((*(G4Element::GetElementTable()))[i], dirName);
184 //while(!((*theElastic)[i])->Register(theFS)) ;
185 ((*theElastic)[i])->Register(theFS) ;
186 }
187 delete theFS;
188 hpmanager->RegisterElasticFinalStates( theElastic );
189
190 }
191 numEle = (G4int)G4Element::GetNumberOfElements();
192}
G4FindDataDir
const char * G4FindDataDir(const char *)
Definition: G4FindDataDir.cc:93
G4Element::GetNumberOfElements
static size_t GetNumberOfElements()
Definition: G4Element.cc:410
G4ParticleHPManager::RegisterElasticFinalStates
void RegisterElasticFinalStates(std::vector< G4ParticleHPChannel * > *val)
Definition: G4ParticleHPManager.hh:109
G4ParticleHPManager::GetElasticFinalStates
std::vector< G4ParticleHPChannel * > * GetElasticFinalStates()
Definition: G4ParticleHPManager.hh:108
G4AutoDelete::Register
void Register(T *inst)
Definition: G4AutoDelete.hh:65
G4Threading::IsMasterThread
G4bool IsMasterThread()
Definition: G4Threading.cc:124
lightions::Init
void Init()
Definition: G4IonTable.cc:77

Referenced by G4ParticleHPThermalScattering::BuildPhysicsTable().

◆ DoNotSuspend()

void G4ParticleHPElastic::DoNotSuspend ( )
inline

Definition at line 64 of file G4ParticleHPElastic.hh.

64{overrideSuspension = true;}

◆ GetFatalEnergyCheckLevels()

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

Reimplemented from G4HadronicInteraction.

Definition at line 138 of file G4ParticleHPElastic.cc.

139{
140 // max energy non-conservation is mass of heavy nucleus
141 return std::pair<G4double, G4double>(10.0*perCent, 350.0*CLHEP::GeV);
142}

◆ GetNiso()

G4int G4ParticleHPElastic::GetNiso ( )
inline

Definition at line 62 of file G4ParticleHPElastic.hh.

62{return ((*theElastic)[0])->GetNiso();};

◆ GetVerboseLevel()

G4int G4ParticleHPElastic::GetVerboseLevel ( ) const

Definition at line 145 of file G4ParticleHPElastic.cc.

146{
147 return G4ParticleHPManager::GetInstance()->GetVerboseLevel();
148}

◆ ModelDescription()

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

Reimplemented from G4HadronicInteraction.

Definition at line 195 of file G4ParticleHPElastic.cc.

196{
197 outFile << "High Precision model based on Evaluated Nuclear Data Files (ENDF) for inelastic reaction of neutrons below 20MeV\n";
198}

◆ SetVerboseLevel()

void G4ParticleHPElastic::SetVerboseLevel ( G4int  newValue)

Definition at line 151 of file G4ParticleHPElastic.cc.

152{
153 G4ParticleHPManager::GetInstance()->SetVerboseLevel(newValue);
154}
The documentation for this class was generated from the following files: