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

#include <G4LENDInelastic.hh>

+ Inheritance diagram for G4LENDInelastic:

Public Member Functions

 G4LENDInelastic (G4ParticleDefinition *pd)
 
 ~G4LENDInelastic ()
 
G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
- Public Member Functions inherited from G4LENDModel
 G4LENDModel (G4String name="LENDModel")
 
 ~G4LENDModel ()
 
virtual G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
 
void ChangeDefaultEvaluation (G4String name)
 
void AllowNaturalAbundanceTarget ()
 
void AllowAnyCandidateTarget ()
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
void DumpLENDTargetInfo (G4bool force=false)
 
- 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 G4LENDModel
void create_used_target_map ()
 
void recreate_used_target_map ()
 
G4GIDI_target * get_target_from_map (G4int nuclear_code)
 
G4HadFinalState * returnUnchanged (const G4HadProjectile &aTrack, G4HadFinalState *theResult)
 
- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4LENDModel
G4ParticleDefinition * proj
 
G4LENDManager * lend_manager
 
std::map< G4int, G4LENDUsedTarget * > usedTarget_map
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 47 of file G4LENDInelastic.hh.

Constructor & Destructor Documentation

◆ G4LENDInelastic()

G4LENDInelastic::G4LENDInelastic ( G4ParticleDefinition *  pd)
inline

Definition at line 52 of file G4LENDInelastic.hh.

53 :G4LENDModel( "LENDInelastic" )
54 {
55 proj = pd;
56
57// theModelName = "LENDInelastic for ";
58// theModelName += proj->GetParticleName();
59 create_used_target_map();
60
61 G4HadronicInteraction* p =
62 G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
63 G4VPreCompoundModel* pre = static_cast<G4VPreCompoundModel*>(p);
64 if(!pre) { pre = new G4PreCompoundModel(); }
65 preco = pre;
66 };
G4HadronicInteractionRegistry::FindModel
G4HadronicInteraction * FindModel(const G4String &name)
Definition: G4HadronicInteractionRegistry.cc:110
G4HadronicInteractionRegistry::Instance
static G4HadronicInteractionRegistry * Instance()
Definition: G4HadronicInteractionRegistry.cc:36
G4LENDModel::create_used_target_map
void create_used_target_map()
Definition: G4LENDModel.cc:93
G4LENDModel::proj
G4ParticleDefinition * proj
Definition: G4LENDModel.hh:83

◆ ~G4LENDInelastic()

G4LENDInelastic::~G4LENDInelastic ( )
inline

Definition at line 68 of file G4LENDInelastic.hh.

68{;};

Member Function Documentation

◆ ApplyYourself()

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

Reimplemented from G4LENDModel.

Definition at line 66 of file G4LENDInelastic.cc.

68{
69 G4ThreeVector projMom = aTrack.Get4Momentum().vect();
70 G4double temp = aTrack.GetMaterial()->GetTemperature();
71
72 G4int iZ = aTarg.GetZ_asInt();
73 G4int iA = aTarg.GetA_asInt();
74 G4int iM = 0;
75 if (aTarg.GetIsotope() != nullptr) iM = aTarg.GetIsotope()->Getm();
76
77 G4double ke = aTrack.GetKineticEnergy();
78
79 G4HadFinalState* theResult = &theParticleChange;
80 theResult->Clear();
81
82 G4GIDI_target* aGIDITarget =
83 get_target_from_map(lend_manager->GetNucleusEncoding(iZ, iA, iM) );
84 if (aGIDITarget == nullptr) {
85// G4cout << " No target found " << G4endl;
86 theParticleChange.Clear();
87 theParticleChange.SetStatusChange(isAlive);
88 theParticleChange.SetEnergyChange(aTrack.GetKineticEnergy());
89 theParticleChange.SetMomentumChange(aTrack.Get4Momentum().vect().unit());
90 return &theParticleChange;
91 }
92// return returnUnchanged(aTrack, theResult);
93
94 // Get GIDI final state products for givern target and projectile
95 G4int loop(0);
96 G4int loopMax = 1000;
97 std::vector<G4GIDI_Product>* products;
98 do {
99 products = aGIDITarget->getOthersFinalState(ke*MeV, temp, MyRNG, NULL);
100 loop++;
101 } while (products == nullptr && loop < loopMax);
102
103 // G4LENDInelastic accepts all light fragments and gammas from GIDI (A < 5)
104 // and removes any heavy fragments which cause large baryon number violation.
105 // Charge and energy non-conservation still occur, but over a large number
106 // of events, this improves on average.
107
108 if (loop > loopMax - 1) {
109// G4cout << " too many loops, return initial state " << G4endl;
110
111 theParticleChange.Clear();
112 theParticleChange.SetStatusChange(isAlive);
113 theParticleChange.SetEnergyChange(aTrack.GetKineticEnergy());
114 theParticleChange.SetMomentumChange(aTrack.Get4Momentum().vect().unit());
115 return &theParticleChange;
116
117// if (aTrack.GetDefinition() == G4Proton::Proton() ||
118// aTrack.GetDefinition() == G4Neutron::Neutron() ) {
119// theResult = preco->ApplyYourself(aTrack, aTarg);
120// } else {
121// theResult = returnUnchanged(aTrack, theResult);
122// }
123
124 } else {
125 G4int iTotZ = iZ + aTrack.GetDefinition()->GetAtomicNumber();
126 G4int iTotA = iA + aTrack.GetDefinition()->GetAtomicMass();
127
128 // Loop over GIDI products and separate light from heavy fragments
129 G4int GZtot(0);
130 G4int GAtot(0);
131 G4int productA(0);
132 G4int productZ(0);
133 std::vector<G4int> lightProductIndex;
134 std::vector<G4int> heavyProductIndex;
135 for (G4int i = 0; i < int( products->size() ); i++ ) {
136 productA = (*products)[i].A;
137 if (productA < 5) {
138 lightProductIndex.push_back(i);
139 GZtot += (*products)[i].Z;
140 GAtot += productA;
141 } else {
142 heavyProductIndex.push_back(i);
143 }
144 }
145
146 // Randomize order of heavies to correct somewhat for sampling bias
147 // std::random_shuffle(heavyProductIndex.begin(), heavyProductIndex.end() );
148 // std::cout << " Heavy product index before shuffle : " ;
149 // for (G4int i = 0; i < int(heavyProductIndex.size() ); i++) std::cout << heavyProductIndex[i] << ", " ;
150 // std::cout << std::endl;
151
152 auto rng = std::default_random_engine {};
153 std::shuffle(heavyProductIndex.begin(), heavyProductIndex.end(), rng);
154
155 // std::cout << " Heavy product index after shuffle : " ;
156 // for (G4int i = 0; i < int(heavyProductIndex.size() ); i++) std::cout << heavyProductIndex[i] << ", " ;
157 // std::cout << std::endl;
158
159 std::vector<G4int> savedHeavyIndex;
160 G4int itest(0);
161 for (G4int i = 0; i < int(heavyProductIndex.size() ); i++) {
162 itest = heavyProductIndex[i];
163 productA = (*products)[itest].A;
164 productZ = (*products)[itest].Z;
165 if ((GAtot + productA <= iTotA) && (GZtot + productZ <= iTotZ) ) {
166 savedHeavyIndex.push_back(itest);
167 GZtot += productZ;
168 GAtot += productA;
169 }
170 }
171
172/*
173 G4cout << " saved light products = ";
174 for (G4int k = 0; k < int(lightProductIndex.size() ); k++ ) {
175 itest = lightProductIndex[k];
176 G4cout << "(" << (*products)[itest].Z << ", " << (*products)[itest].A << "), ";
177 }
178 G4cout << G4endl;
179
180 G4cout << " saved heavy products = ";
181 for (G4int k = 0; k < int(savedHeavyIndex.size() ); k++ ) {
182 itest = savedHeavyIndex[k];
183 G4cout << "(" << (*products)[itest].Z << ", " << (*products)[itest].A << "), ";
184 }
185 G4cout << G4endl;
186*/
187 // Now convert saved products to Geant4 particles
188 // Note that, at least for heavy fragments, GIDI masses and Geant4 masses
189 // have slightly different values.
190
191 G4DynamicParticle* theSec = nullptr;
192 G4ThreeVector Psum;
193 for (G4int i = 0; i < int(lightProductIndex.size() ); i++) {
194 itest = lightProductIndex[i];
195 productZ = (*products)[itest].Z;
196 productA = (*products)[itest].A;
197 theSec = new G4DynamicParticle();
198 if (productA == 1 && productZ == 0) {
199 theSec->SetDefinition(G4Neutron::Neutron() );
200 } else if (productA == 1 && productZ == 1) {
201 theSec->SetDefinition(G4Proton::Proton() );
202 } else if (productA == 2 && productZ == 1) {
203 theSec->SetDefinition(G4Deuteron::Deuteron() );
204 } else if (productA == 3 && productZ == 1) {
205 theSec->SetDefinition(G4Triton::Triton() );
206 } else if (productA == 4 && productZ == 2) {
207 theSec->SetDefinition(G4Alpha::Alpha() );
208 } else {
209 theSec->SetDefinition(G4Gamma::Gamma() );
210 }
211
212 G4ThreeVector momentum((*products)[itest].px*MeV,
213 (*products)[itest].py*MeV,
214 (*products)[itest].pz*MeV );
215 Psum += momentum;
216 theSec->SetMomentum(momentum);
217// theResult->AddSecondary(theSec);
218 theParticleChange.AddSecondary(theSec);
219 }
220
221 G4int productM(0);
222 for (G4int i = 0; i < int(savedHeavyIndex.size() ); i++) {
223 itest = savedHeavyIndex[i];
224 productZ = (*products)[itest].Z;
225 productA = (*products)[itest].A;
226 productM = (*products)[itest].m;
227 theSec = new G4DynamicParticle();
228 theSec->SetDefinition(G4IonTable::GetIonTable()->GetIon(productZ,
229 productA,
230 productM) );
231 G4ThreeVector momentum((*products)[itest].px*MeV,
232 (*products)[itest].py*MeV,
233 (*products)[itest].pz*MeV );
234 Psum += momentum;
235 theSec->SetMomentum(momentum);
236// theResult->AddSecondary(theSec);
237 theParticleChange.AddSecondary(theSec);
238 }
239
240 // Create heavy fragment if necessary to try to balance A, Z
241 // Note: this step is only required to prevent warnings at the process level
242 // where "catastrophic" non-conservation tolerances are set to 1 GeV.
243 // The residual generated will not necessarily be the one that would
244 // occur in the actual reaction.
245 if (iTotA - GAtot > 1) {
246 theSec = new G4DynamicParticle();
247 if (iTotZ == GZtot) {
248 // Special case when a nucleus of only neutrons is requested
249 // Violate charge conservation and set Z = 1
250 // G4cout << " Z = 1, A = "<< iTotA - GAtot << " created " << G4endl;
251 theSec->SetDefinition(G4IonTable::GetIonTable()->GetIon(1, iTotA-GAtot, 0) );
252 } else {
253 theSec->SetDefinition(G4IonTable::GetIonTable()->GetIon(iTotZ-GZtot, iTotA-GAtot, 0) );
254 }
255 theSec->SetMomentum(projMom - Psum);
256// theResult->AddSecondary(theSec);
257 theParticleChange.AddSecondary(theSec);
258 }
259 } // loop OK
260
261 delete products;
262// theResult->SetStatusChange( stopAndKill );
263 theParticleChange.SetStatusChange( stopAndKill );
264// return theResult;
265 return &theParticleChange;
266}
isAlive
@ isAlive
Definition: G4HadFinalState.hh:42
stopAndKill
@ stopAndKill
Definition: G4HadFinalState.hh:42
MyRNG
double MyRNG(void *)
Definition: G4LENDModel.cc:45
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
G4Alpha::Alpha
static G4Alpha * Alpha()
Definition: G4Alpha.cc:88
G4Deuteron::Deuteron
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:93
G4DynamicParticle::SetDefinition
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
Definition: G4DynamicParticle.cc:325
G4DynamicParticle::SetMomentum
void SetMomentum(const G4ThreeVector &momentum)
Definition: G4DynamicParticle.cc:393
G4GIDI_target
Definition: G4GIDI_target.hh:58
G4GIDI_target::getOthersFinalState
std::vector< G4GIDI_Product > * getOthersFinalState(double e_in, double temperature, double(*rng)(void *), void *rngState)
Definition: G4GIDI_target.cc:385
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:85
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.hh:67
G4HadFinalState::AddSecondary
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
Definition: G4HadFinalState.hh:61
G4HadFinalState::SetEnergyChange
void SetEnergyChange(G4double anEnergy)
Definition: G4HadFinalState.cc:39
G4HadFinalState::SetMomentumChange
void SetMomentumChange(const G4ThreeVector &aV)
Definition: G4HadFinalState.hh:56
G4HadProjectile::GetMaterial
const G4Material * GetMaterial() const
Definition: G4HadProjectile.hh:81
G4HadProjectile::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4HadProjectile.hh:86
G4HadProjectile::GetKineticEnergy
G4double GetKineticEnergy() const
Definition: G4HadProjectile.hh:116
G4HadProjectile::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4HadProjectile.hh:91
G4IonTable::GetIonTable
static G4IonTable * GetIonTable()
Definition: G4IonTable.cc:170
G4LENDManager::GetNucleusEncoding
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
Definition: G4LENDManager.cc:343
G4LENDModel::lend_manager
G4LENDManager * lend_manager
Definition: G4LENDModel.hh:84
G4LENDModel::get_target_from_map
G4GIDI_target * get_target_from_map(G4int nuclear_code)
Definition: G4LENDModel.cc:267
G4Material::GetTemperature
G4double GetTemperature() const
Definition: G4Material.hh:180
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:103
G4ParticleDefinition::GetAtomicNumber
G4int GetAtomicNumber() const
G4ParticleDefinition::GetAtomicMass
G4int GetAtomicMass() const
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:92
G4Triton::Triton
static G4Triton * Triton()
Definition: G4Triton.cc:94
The documentation for this class was generated from the following files: