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

#include <G4GEMProbability.hh>

+ Inheritance diagram for G4GEMProbability:

Public Member Functions

 G4GEMProbability (G4int anA, G4int aZ, G4double aSpin)
 
virtual ~G4GEMProbability ()
 
virtual G4double EmissionProbability (const G4Fragment &fragment, G4double maxKineticEnergy)
 
void Dump () const
 
G4double GetSpin (void) const
 
void SetCoulomBarrier (const G4VCoulombBarrier *aCoulombBarrierStrategy)
 
G4double GetCoulombBarrier (const G4Fragment &fragment) const
 
G4double CalcAlphaParam (const G4Fragment &) const
 
G4double CalcBetaParam (const G4Fragment &) const
 
- Public Member Functions inherited from G4VEmissionProbability
 G4VEmissionProbability (G4int Z, G4int A)
 
virtual ~G4VEmissionProbability ()=default
 
void Initialise ()
 
virtual G4double EmissionProbability (const G4Fragment &fragment, G4double anEnergy)
 
virtual G4double ComputeProbability (G4double anEnergy, G4double CB)
 
G4int GetZ (void) const
 
G4int GetA (void) const
 
void SetDecayKinematics (G4int rZ, G4int rA, G4double rmass, G4double fmass)
 
G4double GetRecoilExcitation () const
 
void SetEvapExcitation (G4double exc)
 
G4double GetProbability () const
 
void ResetProbability ()
 
G4double SampleEnergy ()
 
 G4VEmissionProbability (const G4VEmissionProbability &right)=delete
 
const G4VEmissionProbability & operator= (const G4VEmissionProbability &right)=delete
 
G4bool operator== (const G4VEmissionProbability &right) const =delete
 
G4bool operator!= (const G4VEmissionProbability &right) const =delete
 

Protected Attributes

G4double fPlanck
 
std::vector< G4double > ExcitEnergies
 
std::vector< G4double > ExcitSpins
 
std::vector< G4double > ExcitLifetimes
 
- Protected Attributes inherited from G4VEmissionProbability
G4NuclearLevelData * pNuclearLevelData
 
G4Pow * pG4pow
 
G4int OPTxs
 
G4int pVerbose
 
G4int theZ
 
G4int theA
 
G4int resZ = 0
 
G4int resA = 0
 
G4double pMass = 0.0
 
G4double pEvapMass = 0.0
 
G4double pResMass = 0.0
 
G4double pProbability = 0.0
 
G4double pTolerance = 0.0
 

Additional Inherited Members

- Protected Member Functions inherited from G4VEmissionProbability
void ResetIntegrator (size_t nbin, G4double de, G4double eps)
 
G4double IntegrateProbability (G4double elow, G4double ehigh, G4double CB)
 

Detailed Description

Definition at line 53 of file G4GEMProbability.hh.

Constructor & Destructor Documentation

◆ G4GEMProbability()

G4GEMProbability::G4GEMProbability ( G4int  anA,
G4int  aZ,
G4double  aSpin 
)

Definition at line 58 of file G4GEMProbability.cc.

59 : G4VEmissionProbability(aZ, anA), Spin(aSpin),
60 theCoulombBarrierPtr(nullptr)
61{
62 theEvapLDPptr = new G4EvaporationLevelDensityParameter;
63 fG4pow = G4Pow::GetInstance();
64 fPlanck= CLHEP::hbar_Planck*fG4pow->logZ(2);
65 fNucData = G4NuclearLevelData::GetInstance();
66}
G4NuclearLevelData::GetInstance
static G4NuclearLevelData * GetInstance()
Definition: G4NuclearLevelData.cc:423
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:41
G4Pow::logZ
G4double logZ(G4int Z) const
Definition: G4Pow.hh:137

◆ ~G4GEMProbability()

G4GEMProbability::~G4GEMProbability ( )
virtual

Definition at line 68 of file G4GEMProbability.cc.

69{
70 delete theEvapLDPptr;
71}

Member Function Documentation

◆ CalcAlphaParam()

G4double G4GEMProbability::CalcAlphaParam ( const G4Fragment &  fragment) const
inline

Definition at line 165 of file G4GEMProbability.hh.

166{
167 //JMQ 190709 values according to Furihata's paper (based on notes added
168 //on proof in Dostrovskii's paper)
169 G4double res;
170 if(theZ == 0) {
171 res = 0.76+1.93/fG4pow->Z13(fragment.GetA_asInt()-theA);
172 } else {
173 res = 1.0 + CCoeficient(fragment.GetZ_asInt()-theZ);
174 }
175 return res;
176}
G4double
double G4double
Definition: G4Types.hh:83
G4Fragment::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Fragment.hh:289
G4Fragment::GetA_asInt
G4int GetA_asInt() const
Definition: G4Fragment.hh:284
G4Pow::Z13
G4double Z13(G4int Z) const
Definition: G4Pow.hh:123

Referenced by CalcBetaParam().

◆ CalcBetaParam()

G4double G4GEMProbability::CalcBetaParam ( const G4Fragment &  fragment) const
inline

Definition at line 179 of file G4GEMProbability.hh.

180{
181 //JMQ 190709 values according to Furihata's paper (based on notes added
182 //on proof in Dostrovskii's paper)
183 G4double res;
184 if(theZ == 0) {
185 res = (1.66/fG4pow->Z23(fragment.GetA_asInt()-theA)-0.05)*CLHEP::MeV/
186 CalcAlphaParam(fragment);
187 } else {
188 res = -GetCoulombBarrier(fragment);
189 }
190 return res;
191}
G4GEMProbability::GetCoulombBarrier
G4double GetCoulombBarrier(const G4Fragment &fragment) const
Definition: G4GEMProbability.hh:136
G4GEMProbability::CalcAlphaParam
G4double CalcAlphaParam(const G4Fragment &) const
Definition: G4GEMProbability.hh:165
G4Pow::Z23
G4double Z23(G4int Z) const
Definition: G4Pow.hh:125

◆ Dump()

void G4GEMProbability::Dump ( ) const

Definition at line 261 of file G4GEMProbability.cc.

262{
263 G4double mass = G4NucleiProperties::GetNuclearMass(theA, theZ);
264 G4double efermi = 0.0;
265 if(theA > 1) {
266 efermi = G4NucleiProperties::GetNuclearMass(theA-1, theZ)
267 + neutron_mass_c2 - mass;
268 }
269 std::size_t nlev = ExcitEnergies.size();
270 G4cout << "GEM: List of Excited States for Isotope Z= "
271 << theZ << " A= " << theA << " Nlevels= " << nlev
272 << " Efermi(MeV)= " << efermi
273 << G4endl;
274 for(std::size_t i=0; i< nlev; ++i) {
275 G4cout << "Z= " << theZ << " A= " << theA
276 << " Mass(GeV)= " << mass/GeV
277 << " Eexc(MeV)= " << ExcitEnergies[i]
278 << " Time(ns)= " << ExcitLifetimes[i]/ns
279 << G4endl;
280 }
281 G4cout << G4endl;
282}
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4GEMProbability::ExcitEnergies
std::vector< G4double > ExcitEnergies
Definition: G4GEMProbability.hh:114
G4GEMProbability::ExcitLifetimes
std::vector< G4double > ExcitLifetimes
Definition: G4GEMProbability.hh:120
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:50
ns
#define ns(x)
Definition: xmltok.c:1649

Referenced by G4GEMChannel::Dump().

◆ EmissionProbability()

G4double G4GEMProbability::EmissionProbability ( const G4Fragment &  fragment,
G4double  maxKineticEnergy 
)
virtual

Reimplemented from G4VEmissionProbability.

Definition at line 73 of file G4GEMProbability.cc.

75{
76 G4double probability = 0.0;
77
78 if (MaximalKineticEnergy > 0.0 && fragment.GetExcitationEnergy() > 0.0) {
79 G4double CoulombBarrier = GetCoulombBarrier(fragment);
80
81 probability =
82 CalcProbability(fragment,MaximalKineticEnergy,CoulombBarrier);
83
84 // Next there is a loop over excited states for this channel
85 // summing probabilities
86 std::size_t nn = ExcitEnergies.size();
87 if (0 < nn) {
88 G4double SavedSpin = Spin;
89 for (std::size_t i = 0; i <nn; ++i) {
90 Spin = ExcitSpins[i];
91 // substract excitation energies
92 G4double Tmax = MaximalKineticEnergy - ExcitEnergies[i];
93 if (Tmax > 0.0) {
94 G4double width = CalcProbability(fragment,Tmax,CoulombBarrier);
95 //JMQ April 2010 added condition to prevent reported crash
96 // update probability
97 if (width > 0. && fPlanck < width*ExcitLifetimes[i]) {
98 probability += width;
99 }
100 }
101 }
102 // Restore Spin
103 Spin = SavedSpin;
104 }
105 }
106 // Normalization = probability;
107 return probability;
108}
G4Fragment::GetExcitationEnergy
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:312
G4GEMProbability::ExcitSpins
std::vector< G4double > ExcitSpins
Definition: G4GEMProbability.hh:117

◆ GetCoulombBarrier()

G4double G4GEMProbability::GetCoulombBarrier ( const G4Fragment &  fragment) const
inline

Definition at line 136 of file G4GEMProbability.hh.

137{
138 G4double res = 0.0;
139 if (theCoulombBarrierPtr) {
140 G4int Acomp = fragment.GetA_asInt();
141 G4int Zcomp = fragment.GetZ_asInt();
142 res = theCoulombBarrierPtr->GetCoulombBarrier(Acomp-theA, Zcomp-theZ,
143 fragment.GetExcitationEnergy() -
144 fNucData->GetPairingCorrection(Zcomp,Acomp));
145 }
146 return res;
147}
G4int
int G4int
Definition: G4Types.hh:85
G4NuclearLevelData::GetPairingCorrection
G4PairingCorrection * GetPairingCorrection()
Definition: G4NuclearLevelData.cc:610
G4VCoulombBarrier::GetCoulombBarrier
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U=0.0) const =0

Referenced by CalcBetaParam(), and EmissionProbability().

◆ GetSpin()

G4double G4GEMProbability::GetSpin ( void  ) const
inline

Definition at line 124 of file G4GEMProbability.hh.

125{
126 return Spin;
127}

◆ SetCoulomBarrier()

void G4GEMProbability::SetCoulomBarrier ( const G4VCoulombBarrier *  aCoulombBarrierStrategy)
inline

Definition at line 130 of file G4GEMProbability.hh.

131{
132 theCoulombBarrierPtr = aCoulombBarrierStrategy;
133}

Referenced by G4GEMChannel::G4GEMChannel().

Member Data Documentation

◆ ExcitEnergies

std::vector<G4double> G4GEMProbability::ExcitEnergies
protected

Definition at line 114 of file G4GEMProbability.hh.

Referenced by Dump(), EmissionProbability(), G4AlphaGEMProbability::G4AlphaGEMProbability(), G4B10GEMProbability::G4B10GEMProbability(), G4B11GEMProbability::G4B11GEMProbability(), G4B12GEMProbability::G4B12GEMProbability(), G4B13GEMProbability::G4B13GEMProbability(), G4B8GEMProbability::G4B8GEMProbability(), G4Be10GEMProbability::G4Be10GEMProbability(), G4Be11GEMProbability::G4Be11GEMProbability(), G4Be7GEMProbability::G4Be7GEMProbability(), G4Be9GEMProbability::G4Be9GEMProbability(), G4C10GEMProbability::G4C10GEMProbability(), G4C11GEMProbability::G4C11GEMProbability(), G4C12GEMProbability::G4C12GEMProbability(), G4C13GEMProbability::G4C13GEMProbability(), G4C14GEMProbability::G4C14GEMProbability(), G4C15GEMProbability::G4C15GEMProbability(), G4F17GEMProbability::G4F17GEMProbability(), G4F18GEMProbability::G4F18GEMProbability(), G4F19GEMProbability::G4F19GEMProbability(), G4F20GEMProbability::G4F20GEMProbability(), G4F21GEMProbability::G4F21GEMProbability(), G4He6GEMProbability::G4He6GEMProbability(), G4Li6GEMProbability::G4Li6GEMProbability(), G4Li7GEMProbability::G4Li7GEMProbability(), G4Li8GEMProbability::G4Li8GEMProbability(), G4Li9GEMProbability::G4Li9GEMProbability(), G4Mg22GEMProbability::G4Mg22GEMProbability(), G4Mg23GEMProbability::G4Mg23GEMProbability(), G4Mg24GEMProbability::G4Mg24GEMProbability(), G4Mg25GEMProbability::G4Mg25GEMProbability(), G4Mg26GEMProbability::G4Mg26GEMProbability(), G4Mg27GEMProbability::G4Mg27GEMProbability(), G4Mg28GEMProbability::G4Mg28GEMProbability(), G4N12GEMProbability::G4N12GEMProbability(), G4N13GEMProbability::G4N13GEMProbability(), G4N14GEMProbability::G4N14GEMProbability(), G4N15GEMProbability::G4N15GEMProbability(), G4N16GEMProbability::G4N16GEMProbability(), G4N17GEMProbability::G4N17GEMProbability(), G4Na21GEMProbability::G4Na21GEMProbability(), G4Na22GEMProbability::G4Na22GEMProbability(), G4Na23GEMProbability::G4Na23GEMProbability(), G4Na24GEMProbability::G4Na24GEMProbability(), G4Na25GEMProbability::G4Na25GEMProbability(), G4Ne18GEMProbability::G4Ne18GEMProbability(), G4Ne19GEMProbability::G4Ne19GEMProbability(), G4Ne20GEMProbability::G4Ne20GEMProbability(), G4Ne21GEMProbability::G4Ne21GEMProbability(), G4Ne22GEMProbability::G4Ne22GEMProbability(), G4Ne23GEMProbability::G4Ne23GEMProbability(), G4Ne24GEMProbability::G4Ne24GEMProbability(), G4O14GEMProbability::G4O14GEMProbability(), G4O15GEMProbability::G4O15GEMProbability(), G4O16GEMProbability::G4O16GEMProbability(), G4O17GEMProbability::G4O17GEMProbability(), G4O18GEMProbability::G4O18GEMProbability(), and G4O19GEMProbability::G4O19GEMProbability().

◆ ExcitLifetimes

std::vector<G4double> G4GEMProbability::ExcitLifetimes
protected

Definition at line 120 of file G4GEMProbability.hh.

Referenced by Dump(), EmissionProbability(), G4AlphaGEMProbability::G4AlphaGEMProbability(), G4B10GEMProbability::G4B10GEMProbability(), G4B11GEMProbability::G4B11GEMProbability(), G4B12GEMProbability::G4B12GEMProbability(), G4B13GEMProbability::G4B13GEMProbability(), G4B8GEMProbability::G4B8GEMProbability(), G4Be10GEMProbability::G4Be10GEMProbability(), G4Be11GEMProbability::G4Be11GEMProbability(), G4Be7GEMProbability::G4Be7GEMProbability(), G4Be9GEMProbability::G4Be9GEMProbability(), G4C10GEMProbability::G4C10GEMProbability(), G4C11GEMProbability::G4C11GEMProbability(), G4C12GEMProbability::G4C12GEMProbability(), G4C13GEMProbability::G4C13GEMProbability(), G4C14GEMProbability::G4C14GEMProbability(), G4C15GEMProbability::G4C15GEMProbability(), G4F17GEMProbability::G4F17GEMProbability(), G4F18GEMProbability::G4F18GEMProbability(), G4F19GEMProbability::G4F19GEMProbability(), G4F20GEMProbability::G4F20GEMProbability(), G4F21GEMProbability::G4F21GEMProbability(), G4He6GEMProbability::G4He6GEMProbability(), G4Li6GEMProbability::G4Li6GEMProbability(), G4Li7GEMProbability::G4Li7GEMProbability(), G4Li8GEMProbability::G4Li8GEMProbability(), G4Li9GEMProbability::G4Li9GEMProbability(), G4Mg22GEMProbability::G4Mg22GEMProbability(), G4Mg23GEMProbability::G4Mg23GEMProbability(), G4Mg24GEMProbability::G4Mg24GEMProbability(), G4Mg25GEMProbability::G4Mg25GEMProbability(), G4Mg26GEMProbability::G4Mg26GEMProbability(), G4Mg27GEMProbability::G4Mg27GEMProbability(), G4Mg28GEMProbability::G4Mg28GEMProbability(), G4N12GEMProbability::G4N12GEMProbability(), G4N13GEMProbability::G4N13GEMProbability(), G4N14GEMProbability::G4N14GEMProbability(), G4N15GEMProbability::G4N15GEMProbability(), G4N16GEMProbability::G4N16GEMProbability(), G4N17GEMProbability::G4N17GEMProbability(), G4Na21GEMProbability::G4Na21GEMProbability(), G4Na22GEMProbability::G4Na22GEMProbability(), G4Na23GEMProbability::G4Na23GEMProbability(), G4Na24GEMProbability::G4Na24GEMProbability(), G4Na25GEMProbability::G4Na25GEMProbability(), G4Ne18GEMProbability::G4Ne18GEMProbability(), G4Ne19GEMProbability::G4Ne19GEMProbability(), G4Ne20GEMProbability::G4Ne20GEMProbability(), G4Ne21GEMProbability::G4Ne21GEMProbability(), G4Ne22GEMProbability::G4Ne22GEMProbability(), G4Ne23GEMProbability::G4Ne23GEMProbability(), G4Ne24GEMProbability::G4Ne24GEMProbability(), G4O14GEMProbability::G4O14GEMProbability(), G4O15GEMProbability::G4O15GEMProbability(), G4O16GEMProbability::G4O16GEMProbability(), G4O17GEMProbability::G4O17GEMProbability(), G4O18GEMProbability::G4O18GEMProbability(), and G4O19GEMProbability::G4O19GEMProbability().

◆ ExcitSpins

std::vector<G4double> G4GEMProbability::ExcitSpins
protected

Definition at line 117 of file G4GEMProbability.hh.

Referenced by EmissionProbability(), G4AlphaGEMProbability::G4AlphaGEMProbability(), G4B10GEMProbability::G4B10GEMProbability(), G4B11GEMProbability::G4B11GEMProbability(), G4B12GEMProbability::G4B12GEMProbability(), G4B13GEMProbability::G4B13GEMProbability(), G4B8GEMProbability::G4B8GEMProbability(), G4Be10GEMProbability::G4Be10GEMProbability(), G4Be11GEMProbability::G4Be11GEMProbability(), G4Be7GEMProbability::G4Be7GEMProbability(), G4Be9GEMProbability::G4Be9GEMProbability(), G4C10GEMProbability::G4C10GEMProbability(), G4C11GEMProbability::G4C11GEMProbability(), G4C12GEMProbability::G4C12GEMProbability(), G4C13GEMProbability::G4C13GEMProbability(), G4C14GEMProbability::G4C14GEMProbability(), G4C15GEMProbability::G4C15GEMProbability(), G4F17GEMProbability::G4F17GEMProbability(), G4F18GEMProbability::G4F18GEMProbability(), G4F19GEMProbability::G4F19GEMProbability(), G4F20GEMProbability::G4F20GEMProbability(), G4F21GEMProbability::G4F21GEMProbability(), G4He6GEMProbability::G4He6GEMProbability(), G4Li6GEMProbability::G4Li6GEMProbability(), G4Li7GEMProbability::G4Li7GEMProbability(), G4Li8GEMProbability::G4Li8GEMProbability(), G4Li9GEMProbability::G4Li9GEMProbability(), G4Mg22GEMProbability::G4Mg22GEMProbability(), G4Mg23GEMProbability::G4Mg23GEMProbability(), G4Mg24GEMProbability::G4Mg24GEMProbability(), G4Mg25GEMProbability::G4Mg25GEMProbability(), G4Mg26GEMProbability::G4Mg26GEMProbability(), G4Mg27GEMProbability::G4Mg27GEMProbability(), G4Mg28GEMProbability::G4Mg28GEMProbability(), G4N12GEMProbability::G4N12GEMProbability(), G4N13GEMProbability::G4N13GEMProbability(), G4N14GEMProbability::G4N14GEMProbability(), G4N15GEMProbability::G4N15GEMProbability(), G4N16GEMProbability::G4N16GEMProbability(), G4N17GEMProbability::G4N17GEMProbability(), G4Na21GEMProbability::G4Na21GEMProbability(), G4Na22GEMProbability::G4Na22GEMProbability(), G4Na23GEMProbability::G4Na23GEMProbability(), G4Na24GEMProbability::G4Na24GEMProbability(), G4Na25GEMProbability::G4Na25GEMProbability(), G4Ne18GEMProbability::G4Ne18GEMProbability(), G4Ne19GEMProbability::G4Ne19GEMProbability(), G4Ne20GEMProbability::G4Ne20GEMProbability(), G4Ne21GEMProbability::G4Ne21GEMProbability(), G4Ne22GEMProbability::G4Ne22GEMProbability(), G4Ne23GEMProbability::G4Ne23GEMProbability(), G4Ne24GEMProbability::G4Ne24GEMProbability(), G4O14GEMProbability::G4O14GEMProbability(), G4O15GEMProbability::G4O15GEMProbability(), G4O16GEMProbability::G4O16GEMProbability(), G4O17GEMProbability::G4O17GEMProbability(), G4O18GEMProbability::G4O18GEMProbability(), and G4O19GEMProbability::G4O19GEMProbability().

◆ fPlanck

G4double G4GEMProbability::fPlanck
protected

Definition at line 111 of file G4GEMProbability.hh.

Referenced by EmissionProbability(), G4B10GEMProbability::G4B10GEMProbability(), G4B11GEMProbability::G4B11GEMProbability(), G4B12GEMProbability::G4B12GEMProbability(), G4B8GEMProbability::G4B8GEMProbability(), G4Be10GEMProbability::G4Be10GEMProbability(), G4Be11GEMProbability::G4Be11GEMProbability(), G4Be7GEMProbability::G4Be7GEMProbability(), G4Be9GEMProbability::G4Be9GEMProbability(), G4C10GEMProbability::G4C10GEMProbability(), G4C11GEMProbability::G4C11GEMProbability(), G4C12GEMProbability::G4C12GEMProbability(), G4C13GEMProbability::G4C13GEMProbability(), G4C14GEMProbability::G4C14GEMProbability(), G4C15GEMProbability::G4C15GEMProbability(), G4F17GEMProbability::G4F17GEMProbability(), G4F18GEMProbability::G4F18GEMProbability(), G4F19GEMProbability::G4F19GEMProbability(), G4F20GEMProbability::G4F20GEMProbability(), G4GEMProbability(), G4He6GEMProbability::G4He6GEMProbability(), G4Li6GEMProbability::G4Li6GEMProbability(), G4Li7GEMProbability::G4Li7GEMProbability(), G4Li8GEMProbability::G4Li8GEMProbability(), G4Li9GEMProbability::G4Li9GEMProbability(), G4Mg24GEMProbability::G4Mg24GEMProbability(), G4Mg25GEMProbability::G4Mg25GEMProbability(), G4N12GEMProbability::G4N12GEMProbability(), G4N13GEMProbability::G4N13GEMProbability(), G4N14GEMProbability::G4N14GEMProbability(), G4N15GEMProbability::G4N15GEMProbability(), G4N16GEMProbability::G4N16GEMProbability(), G4Na21GEMProbability::G4Na21GEMProbability(), G4Na22GEMProbability::G4Na22GEMProbability(), G4Na23GEMProbability::G4Na23GEMProbability(), G4Ne18GEMProbability::G4Ne18GEMProbability(), G4Ne19GEMProbability::G4Ne19GEMProbability(), G4Ne20GEMProbability::G4Ne20GEMProbability(), G4Ne21GEMProbability::G4Ne21GEMProbability(), G4O14GEMProbability::G4O14GEMProbability(), G4O15GEMProbability::G4O15GEMProbability(), G4O16GEMProbability::G4O16GEMProbability(), G4O17GEMProbability::G4O17GEMProbability(), G4O18GEMProbability::G4O18GEMProbability(), and G4O19GEMProbability::G4O19GEMProbability().

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