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
G4GeneratorPrecompoundInterface Class Reference

#include <G4GeneratorPrecompoundInterface.hh>

+ Inheritance diagram for G4GeneratorPrecompoundInterface:

Public Member Functions

 G4GeneratorPrecompoundInterface (G4VPreCompoundModel *p=0)
 
virtual ~G4GeneratorPrecompoundInterface ()
 
virtual G4HadFinalState * ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
 
virtual G4ReactionProductVector * Propagate (G4KineticTrackVector *theSecondaries, G4V3DNucleus *theNucleus)
 
void SetCaptureThreshold (G4double)
 
virtual void PropagateModelDescription (std::ostream &) const
 
- Public Member Functions inherited from G4VIntraNuclearTransportModel
 G4VIntraNuclearTransportModel (const G4String &modelName="CascadeModel", G4VPreCompoundModel *ptr=0)
 
virtual ~G4VIntraNuclearTransportModel ()
 
virtual G4ReactionProductVector * Propagate (G4KineticTrackVector *theSecondaries, G4V3DNucleus *theNucleus)=0
 
void SetDeExcitation (G4VPreCompoundModel *ptr)
 
void Set3DNucleus (G4V3DNucleus *const value)
 
void SetPrimaryProjectile (const G4HadProjectile &aPrimary)
 
const G4String & GetModelName () const
 
virtual void ModelDescription (std::ostream &outFile) const
 
virtual void PropagateModelDescription (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)=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 G4VIntraNuclearTransportModel
G4V3DNucleus * Get3DNucleus () const
 
G4VPreCompoundModel * GetDeExcitation () const
 
const G4HadProjectile * GetPrimaryProjectile () const
 
- Protected Member Functions inherited from G4HadronicInteraction
void SetModelName (const G4String &nam)
 
G4bool IsBlocked () const
 
void Block ()
 
- Protected Attributes inherited from G4VIntraNuclearTransportModel
G4String theTransportModelName
 
G4V3DNucleus * the3DNucleus
 
G4VPreCompoundModel * theDeExcitation
 
const G4HadProjectile * thePrimaryProjectile
 
- Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState theParticleChange
 
G4int verboseLevel
 
G4double theMinEnergy
 
G4double theMaxEnergy
 
G4bool isBlocked
 

Detailed Description

Definition at line 58 of file G4GeneratorPrecompoundInterface.hh.

Constructor & Destructor Documentation

◆ G4GeneratorPrecompoundInterface()

G4GeneratorPrecompoundInterface::G4GeneratorPrecompoundInterface ( G4VPreCompoundModel *  p = 0)

Definition at line 59 of file G4GeneratorPrecompoundInterface.cc.

60: CaptureThreshold(10*MeV)
61{
62 proton = G4Proton::Proton();
63 neutron = G4Neutron::Neutron();
64 if(preModel) { SetDeExcitation(preModel); }
65 else {
66 G4HadronicInteraction* hadi =
67 G4HadronicInteractionRegistry::Instance()->FindModel("PRECO");
68 G4VPreCompoundModel* pre = static_cast<G4VPreCompoundModel*>(hadi);
69 if(!pre) { pre = new G4PreCompoundModel(); }
70 SetDeExcitation(pre);
71 }
72}
G4HadronicInteractionRegistry::FindModel
G4HadronicInteraction * FindModel(const G4String &name)
Definition: G4HadronicInteractionRegistry.cc:104
G4HadronicInteractionRegistry::Instance
static G4HadronicInteractionRegistry * Instance()
Definition: G4HadronicInteractionRegistry.cc:36
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93
G4VIntraNuclearTransportModel::SetDeExcitation
void SetDeExcitation(G4VPreCompoundModel *ptr)
Definition: G4VIntraNuclearTransportModel.hh:131

◆ ~G4GeneratorPrecompoundInterface()

G4GeneratorPrecompoundInterface::~G4GeneratorPrecompoundInterface ( )
virtual

Definition at line 74 of file G4GeneratorPrecompoundInterface.cc.

75{
76}

Member Function Documentation

◆ ApplyYourself()

G4HadFinalState * G4GeneratorPrecompoundInterface::ApplyYourself ( const G4HadProjectile &  aTrack,
G4Nucleus &  targetNucleus 
)
virtual

Implements G4HadronicInteraction.

Definition at line 234 of file G4GeneratorPrecompoundInterface.cc.

236{
237 G4cout << "G4GeneratorPrecompoundInterface: ApplyYourself interface called stand-allone."
238 << G4endl;
239 G4cout << "This class is only a mediator between generator and precompound"<<G4endl;
240 G4cout << "Please remove from your physics list."<<G4endl;
241 throw G4HadronicException(__FILE__, __LINE__, "SEVERE: G4GeneratorPrecompoundInterface model interface called stand-allone.");
242 return new G4HadFinalState;
243}
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout

◆ Propagate()

G4ReactionProductVector * G4GeneratorPrecompoundInterface::Propagate ( G4KineticTrackVector *  theSecondaries,
G4V3DNucleus *  theNucleus 
)
virtual

Implements G4VIntraNuclearTransportModel.

Definition at line 84 of file G4GeneratorPrecompoundInterface.cc.

86{
87 G4ReactionProductVector * theTotalResult = new G4ReactionProductVector;
88
89 // decay the strong resonances
90 G4DecayKineticTracks decay(theSecondaries);
91
92 // prepare the fragment
93 G4int anA=theNucleus->GetMassNumber();
94 G4int aZ=theNucleus->GetCharge();
95 G4int numberOfEx = 0;
96 G4int numberOfCh = 0;
97 G4int numberOfHoles = 0;
98 G4double exEnergy = 0.0;
99 G4double R = theNucleus->GetNuclearRadius();
100 G4ThreeVector exciton3Momentum(0.,0.,0.);
101 G4ThreeVector captured3Momentum(0.,0.,0.);
102 G4ThreeVector wounded3Momentum(0.,0.,0.);
103
104 // loop over secondaries
105#ifdef exactExcitationEnergy
106 G4LorentzVector secondary4Momemtum(0,0,0,0);
107#endif
108 G4KineticTrackVector::iterator iter;
109 for(iter=theSecondaries->begin(); iter !=theSecondaries->end(); ++iter)
110 {
111 G4ParticleDefinition* part = (*iter)->GetDefinition();
112 G4double e = (*iter)->Get4Momentum().e();
113 G4double mass = (*iter)->Get4Momentum().mag();
114 G4ThreeVector mom = (*iter)->Get4Momentum().vect();
115 if((part != proton && part != neutron) ||
116 (e > mass + CaptureThreshold) ||
117 ((*iter)->GetPosition().mag() > R)) {
118 G4ReactionProduct * theNew = new G4ReactionProduct(part);
119 theNew->SetMomentum(mom);
120 theNew->SetTotalEnergy(e);
121 theTotalResult->push_back(theNew);
122#ifdef exactExcitationEnergy
123 secondary4Momemtum += (*iter)->Get4Momentum();
124#endif
125 } else {
126 // within the nucleus, neutron or proton
127 // now calculate A, Z of the fragment, momentum, number of exciton states
128 ++anA;
129 ++numberOfEx;
130 G4int Z = G4int(part->GetPDGCharge()/eplus + 0.1);
131 aZ += Z;
132 numberOfCh += Z;
133 captured3Momentum += mom;
134 exEnergy += (e - mass);
135 }
136 delete (*iter);
137 }
138 delete theSecondaries;
139
140 // loop over wounded nucleus
141 G4Nucleon * theCurrentNucleon =
142 theNucleus->StartLoop() ? theNucleus->GetNextNucleon() : 0;
143 while(theCurrentNucleon) {
144 if(theCurrentNucleon->AreYouHit()) {
145 ++numberOfHoles;
146 ++numberOfEx;
147 --anA;
148 aZ -= G4int(theCurrentNucleon->GetDefinition()->GetPDGCharge()/eplus + 0.1);
149 wounded3Momentum += theCurrentNucleon->Get4Momentum().vect();
150 //G4cout << "hit nucleon " << theCurrentNucleon->Get4Momentum() << G4endl;
151 exEnergy += theCurrentNucleon->GetBindingEnergy();
152 }
153 theCurrentNucleon = theNucleus->GetNextNucleon();
154 }
155 exciton3Momentum = captured3Momentum - wounded3Momentum;
156
157 if(anA>0 && aZ>0) {
158 G4double fMass = G4NucleiProperties::GetNuclearMass(anA, aZ);
159#ifdef exactExcitationEnergy
160 // recalculate exEnergy from Energy balance....
161 const G4HadProjectile * primary = GetPrimaryProjectile();
162 G4double Einitial= primary->Get4Momentum().e()
163 + G4NucleiProperties::GetNuclearMass(theNucleus->GetMassNumber(),theNucleus->GetCharge());
164 G4double Efinal = fMass + secondary4Momemtum.e();
165 if ( (Einitial - Efinal) > 0 ) {
166 // G4cout << "G4GPI::Propagate() : positive exact excitation Energy "
167 // << (Einitial - Efinal)/MeV << " MeV, exciton estimate " << exEnergy/MeV << " MeV" << G4endl;
168 exEnergy=Einitial - Efinal;
169 }
170 else {
171 // G4cout << "G4GeneratorPrecompoundInterface::Propagate() : negative exact excitation Energy "
172 // << (Einitial - Efinal)/MeV << " MeV, setting excitation to 0 MeV" << G4endl;
173 exEnergy=0.;
174 }
175#endif
176 fMass += exEnergy;
177 G4ThreeVector balance=primary->Get4Momentum().vect() - secondary4Momemtum.vect() - exciton3Momentum;
178 #ifdef G4GPI_debug_excitation
179 G4cout << "momentum balance init/final " << balance << " value " << balance.mag() << G4endl
180 << "primary / secondaries "<< primary->Get4Momentum() << " / "
181 << secondary4Momemtum << " captured/wounded: " << captured3Momentum << " / " << wounded3Momentum
182 << " exciton " << exciton3Momentum << G4endl
183 << secondary4Momemtum.vect() + exciton3Momentum << G4endl;
184 #endif
185#ifdef exactExcitationEnergy
186 G4LorentzVector exciton4Momentum(exciton3Momentum, fMass);
187#else
188 G4LorentzVector exciton4Momentum(exciton3Momentum,
189 std::sqrt(exciton3Momentum.mag2() + fMass*fMass));
190#endif
191 if ( exEnergy > 0.0 ) { // Need to de-excite the remnant nucleus only if excitation energy > 0.
192 G4Fragment anInitialState(anA, aZ, exciton4Momentum);
193 anInitialState.SetNumberOfParticles(numberOfEx-numberOfHoles);
194 anInitialState.SetNumberOfCharged(numberOfCh);
195 anInitialState.SetNumberOfHoles(numberOfHoles);
196
197 G4ReactionProductVector * aPrecoResult = theDeExcitation->DeExcite(anInitialState);
198 // fill pre-compound part into the result, and return
199 theTotalResult->insert(theTotalResult->end(),aPrecoResult->begin(),aPrecoResult->end() );
200 delete aPrecoResult;
201
202 } else { // No/negative excitation energy, we only need to create the remnant nucleus
203 // energy is not conserved, ignore exciton momentum, i.e. remnant nucleus will be at rest
204 G4ParticleDefinition* theKindOfFragment = 0;
205 if (anA == 1 && aZ == 0) {
206 theKindOfFragment = G4Neutron::NeutronDefinition();
207 } else if (anA == 1 && aZ == 1) {
208 theKindOfFragment = G4Proton::ProtonDefinition();
209 } else if (anA == 2 && aZ == 1) {
210 theKindOfFragment = G4Deuteron::DeuteronDefinition();
211 } else if (anA == 3 && aZ == 1) {
212 theKindOfFragment = G4Triton::TritonDefinition();
213 } else if (anA == 3 && aZ == 2) {
214 theKindOfFragment = G4He3::He3Definition();
215 } else if (anA == 4 && aZ == 2) {
216 theKindOfFragment = G4Alpha::AlphaDefinition();;
217 } else {
218 theKindOfFragment =
219 G4ParticleTable::GetParticleTable()->GetIonTable()->GetIon(aZ,anA,0.0);
220 }
221 if (theKindOfFragment != 0) {
222 G4ReactionProduct * theNew = new G4ReactionProduct(theKindOfFragment);
223 theNew->SetMomentum(G4ThreeVector(0.,0.,0.));
224 theNew->SetTotalEnergy(fMass);
225 //theNew->SetFormationTime(??0.??);
226 theTotalResult->push_back(theNew);
227 }
228 }
229 }
230
231 return theTotalResult;
232}
G4ReactionProductVector
std::vector< G4ReactionProduct * > G4ReactionProductVector
Definition: G4ReactionProductVector.hh:43
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
CLHEP::Hep3Vector::mag
double mag() const
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
G4Alpha::AlphaDefinition
static G4Alpha * AlphaDefinition()
Definition: G4Alpha.cc:84
G4Deuteron::DeuteronDefinition
static G4Deuteron * DeuteronDefinition()
Definition: G4Deuteron.cc:89
G4HadProjectile::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4HadProjectile.hh:86
G4He3::He3Definition
static G4He3 * He3Definition()
Definition: G4He3.cc:89
G4IonTable::GetIon
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int J=0)
Definition: G4IonTable.cc:267
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:53
G4Nucleon::AreYouHit
G4bool AreYouHit() const
Definition: G4Nucleon.hh:97
G4Nucleon::Get4Momentum
virtual const G4LorentzVector & Get4Momentum() const
Definition: G4Nucleon.hh:72
G4Nucleon::GetBindingEnergy
G4double GetBindingEnergy() const
Definition: G4Nucleon.hh:75
G4Nucleon::GetDefinition
virtual G4ParticleDefinition * GetDefinition() const
Definition: G4Nucleon.hh:85
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:119
G4ParticleTable::GetParticleTable
static G4ParticleTable * GetParticleTable()
Definition: G4ParticleTable.cc:63
G4ParticleTable::GetIonTable
G4IonTable * GetIonTable()
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4ReactionProduct::SetMomentum
void SetMomentum(const G4double x, const G4double y, const G4double z)
Definition: G4ReactionProduct.cc:166
G4ReactionProduct::SetTotalEnergy
void SetTotalEnergy(const G4double en)
Definition: G4ReactionProduct.hh:138
G4Triton::TritonDefinition
static G4Triton * TritonDefinition()
Definition: G4Triton.cc:90
G4V3DNucleus::GetNextNucleon
virtual G4Nucleon * GetNextNucleon()=0
G4V3DNucleus::GetCharge
virtual G4int GetCharge()=0
G4V3DNucleus::StartLoop
virtual G4bool StartLoop()=0
G4V3DNucleus::GetNuclearRadius
virtual G4double GetNuclearRadius()=0
G4V3DNucleus::GetMassNumber
virtual G4int GetMassNumber()=0
G4VIntraNuclearTransportModel::GetPrimaryProjectile
const G4HadProjectile * GetPrimaryProjectile() const
Definition: G4VIntraNuclearTransportModel.hh:138
G4VPreCompoundModel::DeExcite
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)=0

◆ PropagateModelDescription()

void G4GeneratorPrecompoundInterface::PropagateModelDescription ( std::ostream &  outFile) const
virtual

Reimplemented from G4VIntraNuclearTransportModel.

Definition at line 244 of file G4GeneratorPrecompoundInterface.cc.

245{
246 outFile << "G4GeneratorPrecompoundInterface interfaces a high\n"
247 << "energy model through the wounded nucleus to precompound de-excition.\n"
248 << "Low energy protons and neutron present among secondaries produced by \n"
249 << "the high energy generator and within the nucleus are captured. The wounded\n"
250 << "nucleus and the captured particles form an excited nuclear fragment. This\n"
251 << "fragment is passed to the Geant4 pre-compound model for de-excitation.\n"
252 << "Nuclear de-excitation:\n";
253 // preco
254
255}

◆ SetCaptureThreshold()

void G4GeneratorPrecompoundInterface::SetCaptureThreshold ( G4double  value)
inline

Definition at line 88 of file G4GeneratorPrecompoundInterface.hh.

89{
90 CaptureThreshold=value;
91}
The documentation for this class was generated from the following files: