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

#include <G4LowEIonFragmentation.hh>

+ Inheritance diagram for G4LowEIonFragmentation:

Public Member Functions

 G4LowEIonFragmentation (G4ExcitationHandler *const value)
 
 G4LowEIonFragmentation ()
 
virtual ~G4LowEIonFragmentation ()
 
virtual G4HadFinalState * ApplyYourself (const G4HadProjectile &thePrimary, G4Nucleus &theNucleus)
 
- 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
 

Static Public Member Functions

static G4double GetCrossSection ()
 

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 50 of file G4LowEIonFragmentation.hh.

Constructor & Destructor Documentation

◆ G4LowEIonFragmentation() [1/2]

G4LowEIonFragmentation::G4LowEIonFragmentation ( G4ExcitationHandler *const  value)

Definition at line 53 of file G4LowEIonFragmentation.cc.

54{
55 theHandler = value;
56 theModel = new G4PreCompoundModel(theHandler);
57 proton = G4Proton::Proton();
58}
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93

◆ G4LowEIonFragmentation() [2/2]

G4LowEIonFragmentation::G4LowEIonFragmentation ( )

Definition at line 60 of file G4LowEIonFragmentation.cc.

61{
62 theHandler = new G4ExcitationHandler;
63 theModel = new G4PreCompoundModel(theHandler);
64 proton = G4Proton::Proton();
65}

◆ ~G4LowEIonFragmentation()

G4LowEIonFragmentation::~G4LowEIonFragmentation ( )
virtual

Definition at line 67 of file G4LowEIonFragmentation.cc.

68{
69 delete theModel;
70}

Member Function Documentation

◆ ApplyYourself()

G4HadFinalState * G4LowEIonFragmentation::ApplyYourself ( const G4HadProjectile &  thePrimary,
G4Nucleus &  theNucleus 
)
virtual

Implements G4HadronicInteraction.

Definition at line 72 of file G4LowEIonFragmentation.cc.

74{
75 area = 0;
76 // initialize the particle change
77 theResult.Clear();
78 theResult.SetStatusChange( stopAndKill );
79 theResult.SetEnergyChange( 0.0 );
80
81 // Get Target A, Z
82 G4int aTargetA = theNucleus.GetA_asInt();
83 G4int aTargetZ = theNucleus.GetZ_asInt();
84
85 // Get Projectile A, Z
86 G4int aProjectileA = thePrimary.GetDefinition()->GetBaryonNumber();
87 G4int aProjectileZ = G4lrint(thePrimary.GetDefinition()->GetPDGCharge()/eplus);
88
89 // Get Maximum radius of both
90
91 G4Fancy3DNucleus aPrim;
92 aPrim.Init(aProjectileA, aProjectileZ);
93 G4double projectileOuterRadius = aPrim.GetOuterRadius();
94
95 G4Fancy3DNucleus aTarg;
96 aTarg.Init(aTargetA, aTargetZ);
97 G4double targetOuterRadius = aTarg.GetOuterRadius();
98
99 // Get the Impact parameter
100 G4int particlesFromProjectile = 0;
101 G4int chargedFromProjectile = 0;
102 G4double impactParameter = 0;
103 G4double x,y;
104 G4Nucleon * pNucleon;
105 // need at lease one particle from the projectile model beyond the
106 // projectileHorizon.
107 while(0==particlesFromProjectile)
108 {
109 do
110 {
111 x = 2*G4UniformRand() - 1;
112 y = 2*G4UniformRand() - 1;
113 }
114 while(x*x + y*y > 1);
115 impactParameter = std::sqrt(x*x+y*y)*(targetOuterRadius+projectileOuterRadius);
116 ++totalTries;
117 area = pi*(targetOuterRadius+projectileOuterRadius)*
118 (targetOuterRadius+projectileOuterRadius);
119 G4double projectileHorizon = impactParameter-targetOuterRadius;
120
121 // Empirical boundary transparency.
122 G4double empirical = G4UniformRand();
123 if(projectileHorizon > empirical*projectileOuterRadius) { continue; }
124
125 // Calculate the number of nucleons involved in collision
126 // From projectile
127 aPrim.StartLoop();
128 while((pNucleon = aPrim.GetNextNucleon()))
129 {
130 if(pNucleon->GetPosition().y()>projectileHorizon)
131 {
132 // We have one
133 ++particlesFromProjectile;
134 if(pNucleon->GetParticleType() == proton)
135 {
136 ++chargedFromProjectile;
137 }
138 }
139 }
140 }
141 ++hits;
142
143 // From target:
144 G4double targetHorizon = impactParameter-projectileOuterRadius;
145 G4int chargedFromTarget = 0;
146 G4int particlesFromTarget = 0;
147 aTarg.StartLoop();
148 while((pNucleon = aTarg.GetNextNucleon()))
149 {
150 if(pNucleon->GetPosition().y()>targetHorizon)
151 {
152 // We have one
153 ++particlesFromTarget;
154 if(pNucleon->GetParticleType() == proton)
155 {
156 ++chargedFromTarget;
157 }
158 }
159 }
160
161 // Energy sharing between projectile and target.
162 // Note that this is a quite simplistic kinetically.
163 G4ThreeVector momentum = thePrimary.Get4Momentum().vect();
164 G4double w = (G4double)particlesFromProjectile/(G4double)aProjectileA;
165
166 G4double projTotEnergy = thePrimary.GetTotalEnergy();
167 G4double targetMass = G4NucleiProperties::GetNuclearMass(aTargetA, aTargetZ);
168 G4LorentzVector fragment4Momentum(momentum*w, projTotEnergy*w + targetMass);
169
170 // take the nucleons and fill the Fragments
171 G4Fragment anInitialState(aTargetA+particlesFromProjectile,
172 aTargetZ+chargedFromProjectile,
173 fragment4Momentum);
174 // M.A. Cortes fix
175 //anInitialState.SetNumberOfParticles(particlesFromProjectile);
176 anInitialState.SetNumberOfExcitedParticle(particlesFromProjectile+particlesFromTarget,
177 chargedFromProjectile + chargedFromTarget);
178 anInitialState.SetNumberOfHoles(particlesFromProjectile+particlesFromTarget,
179 chargedFromProjectile + chargedFromTarget);
180 G4double time = thePrimary.GetGlobalTime();
181 anInitialState.SetCreationTime(time);
182
183 // Fragment the Fragment using Pre-compound
184 G4ReactionProductVector* thePreCompoundResult = theModel->DeExcite(anInitialState);
185
186 // De-excite the projectile using ExcitationHandler
187 G4ReactionProductVector * theExcitationResult = 0;
188 if(particlesFromProjectile < aProjectileA)
189 {
190 G4LorentzVector residual4Momentum(momentum*(1.0-w), projTotEnergy*(1.0-w));
191
192 G4Fragment initialState2(aProjectileA-particlesFromProjectile,
193 aProjectileZ-chargedFromProjectile,
194 residual4Momentum );
195
196 // half of particles are excited (?!)
197 G4int pinit = (aProjectileA-particlesFromProjectile)/2;
198 G4int cinit = (aProjectileZ-chargedFromProjectile)/2;
199
200 initialState2.SetNumberOfExcitedParticle(pinit,cinit);
201 initialState2.SetNumberOfHoles(pinit,cinit);
202 initialState2.SetCreationTime(time);
203
204 theExcitationResult = theHandler->BreakItUp(initialState2);
205 }
206
207 // Fill the particle change and clear intermediate vectors
208 G4int nexc = 0;
209 G4int npre = 0;
210 if(theExcitationResult) { nexc = theExcitationResult->size(); }
211 if(thePreCompoundResult) { npre = thePreCompoundResult->size();}
212
213 if(nexc > 0) {
214 for(G4int k=0; k<nexc; ++k) {
215 G4ReactionProduct* p = (*theExcitationResult)[k];
216 theResult.AddSecondary(new G4DynamicParticle(p->GetDefinition(),p->GetMomentum()));
217 delete p;
218 }
219 }
220
221 if(npre > 0) {
222 for(G4int k=0; k<npre; ++k) {
223 G4ReactionProduct* p = (*thePreCompoundResult)[k];
224 theResult.AddSecondary(new G4DynamicParticle(p->GetDefinition(),p->GetMomentum()));
225 delete p;
226 }
227 }
228
229 delete thePreCompoundResult;
230 delete theExcitationResult;
231
232 // return the particle change
233 return &theResult;
234
235}
stopAndKill
@ stopAndKill
Definition: G4HadFinalState.hh:42
G4ReactionProductVector
std::vector< G4ReactionProduct * > G4ReactionProductVector
Definition: G4ReactionProductVector.hh:43
G4double
double G4double
Definition: G4Types.hh:64
G4int
int G4int
Definition: G4Types.hh:66
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53
CLHEP::Hep3Vector::y
double y() const
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
G4ExcitationHandler::BreakItUp
G4ReactionProductVector * BreakItUp(const G4Fragment &theInitialState) const
Definition: G4ExcitationHandler.cc:118
G4Fancy3DNucleus::GetNextNucleon
G4Nucleon * GetNextNucleon()
Definition: G4Fancy3DNucleus.cc:124
G4Fancy3DNucleus::GetOuterRadius
G4double GetOuterRadius()
Definition: G4Fancy3DNucleus.cc:175
G4Fancy3DNucleus::Init
void Init(G4int theA, G4int theZ)
Definition: G4Fancy3DNucleus.cc:76
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.cc:81
G4HadFinalState::AddSecondary
void AddSecondary(G4DynamicParticle *aP)
Definition: G4HadFinalState.cc:70
G4HadFinalState::SetEnergyChange
void SetEnergyChange(G4double anEnergy)
Definition: G4HadFinalState.cc:42
G4HadProjectile::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4HadProjectile.hh:81
G4HadProjectile::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4HadProjectile.hh:86
G4HadProjectile::GetGlobalTime
G4double GetGlobalTime() const
Definition: G4HadProjectile.hh:113
G4HadProjectile::GetTotalEnergy
G4double GetTotalEnergy() const
Definition: G4HadProjectile.hh:96
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:53
G4Nucleon::GetParticleType
G4ParticleDefinition * GetParticleType() const
Definition: G4Nucleon.hh:84
G4Nucleon::GetPosition
virtual const G4ThreeVector & GetPosition() const
Definition: G4Nucleon.hh:68
G4Nucleus::GetA_asInt
G4int GetA_asInt() const
Definition: G4Nucleus.hh:109
G4Nucleus::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Nucleus.hh:115
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:119
G4PreCompoundModel::DeExcite
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)
Definition: G4PreCompoundModel.cc:184
G4ReactionProduct::GetMomentum
G4ThreeVector GetMomentum() const
Definition: G4ReactionProduct.hh:120
G4ReactionProduct::GetDefinition
G4ParticleDefinition * GetDefinition() const
Definition: G4ReactionProduct.hh:104
G4INCL::Math::pi
const G4double pi
Definition: G4INCLGlobals.hh:70
G4lrint
int G4lrint(double ad)
Definition: templates.hh:163

◆ GetCrossSection()

static G4double G4LowEIonFragmentation::GetCrossSection ( )
inlinestatic

Definition at line 63 of file G4LowEIonFragmentation.hh.

64 {
65 // clog << "area/millibarn = "<<area/millibarn<<G4endl;
66 // clog << "hits = "<<hits<<G4endl;
67 // clog << "totalTries = "<<totalTries<<G4endl;
68 return area*hits/(static_cast<G4double>(totalTries)*CLHEP::millibarn);
69 }
The documentation for this class was generated from the following files: