Geant4-11
Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes
G4PreCompoundIon Class Referenceabstract

#include <G4PreCompoundIon.hh>

Inheritance diagram for G4PreCompoundIon:
G4PreCompoundFragment G4VPreCompoundFragment G4PreCompoundAlpha G4PreCompoundDeuteron G4PreCompoundHe3 G4PreCompoundTriton

Public Member Functions

G4double CalcEmissionProbability (const G4Fragment &aFragment) override
 
 G4PreCompoundIon (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
 
 G4PreCompoundIon (const G4PreCompoundIon &right)=delete
 
G4int GetA () const
 
G4double GetBindingEnergy () const
 
G4double GetEmissionProbability () const
 
G4double GetEnergyThreshold () const
 
const G4LorentzVector & GetMomentum () const
 
G4double GetNuclearMass () const
 
G4ReactionProduct * GetReactionProduct () const
 
G4int GetRestA () const
 
G4double GetRestNuclearMass () const
 
G4int GetRestZ () const
 
G4int GetZ () const
 
void Initialize (const G4Fragment &aFragment)
 
G4bool IsItPossible (const G4Fragment &aFragment) const
 
G4bool operator!= (const G4PreCompoundFragment &right) const =delete
 
G4bool operator!= (const G4PreCompoundIon &right) const =delete
 
G4bool operator!= (const G4VPreCompoundFragment &right) const =delete
 
const G4PreCompoundIon & operator= (const G4PreCompoundIon &right)=delete
 
G4bool operator== (const G4PreCompoundFragment &right) const =delete
 
G4bool operator== (const G4PreCompoundIon &right) const =delete
 
G4bool operator== (const G4VPreCompoundFragment &right) const =delete
 
G4double SampleKineticEnergy (const G4Fragment &aFragment) override
 
void SetMomentum (const G4LorentzVector &value)
 
void SetOPTxs (G4int)
 
void UseSICB (G4bool)
 
virtual ~G4PreCompoundIon ()
 

Protected Member Functions

virtual G4double CoalescenceFactor (G4int A) const =0
 
G4double CrossSection (G4double ekin) const
 
virtual G4double FactorialFactor (G4int N, G4int P) const =0
 
virtual G4double GetAlpha () const =0
 
G4double GetBeta () const override
 
virtual G4double GetRj (G4int NumberParticles, G4int NumberCharged) const =0
 
G4double ProbabilityDistributionFunction (G4double eKin, const G4Fragment &) override
 

Protected Attributes

G4NuclearLevelData * fNucData
 
G4Pow * g4calc
 
G4int OPTxs
 
G4int theA
 
G4double theBindingEnergy
 
G4double theCoulombBarrier
 
G4double theEmissionProbability
 
G4int theFragA
 
G4int theFragZ
 
G4double theMass
 
G4double theMaxKinEnergy
 
G4double theMinKinEnergy
 
G4DeexPrecoParameters * theParameters
 
G4double theReducedMass
 
G4int theResA
 
G4double theResA13
 
G4double theResMass
 
G4int theResZ
 
G4int theZ
 
G4bool useSICB
 

Private Member Functions

G4double GetOpt0 (G4double ekin) const
 
G4double IntegrateEmissionProbability (G4double low, G4double up, const G4Fragment &aFragment)
 

Private Attributes

G4double fact
 
G4int index
 
G4double muu
 
const G4ParticleDefinition * particle
 
G4double probmax
 
G4VCoulombBarrier * theCoulombBarrierPtr
 
G4LorentzVector theMomentum
 

Detailed Description

Definition at line 41 of file G4PreCompoundIon.hh.

Constructor & Destructor Documentation

◆ G4PreCompoundIon() [1/2]

G4PreCompoundIon::G4PreCompoundIon ( const G4ParticleDefinition *  part,
G4VCoulombBarrier *  aCoulombBarrier 
)

Definition at line 47 of file G4PreCompoundIon.cc.

50 : G4PreCompoundFragment(part,aCoulombBarrier)
51{
52 G4double r0 = theParameters->GetR0();
53 fact = 0.75*CLHEP::millibarn/(CLHEP::pi*r0*r0*r0);
54}
G4double
double G4double
Definition: G4Types.hh:83
G4PreCompoundFragment::G4PreCompoundFragment
G4PreCompoundFragment(const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
Definition: G4PreCompoundFragment.cc:42
G4VPreCompoundFragment::theParameters
G4DeexPrecoParameters * theParameters
Definition: G4VPreCompoundFragment.hh:139
CLHEP::millibarn
static constexpr double millibarn
Definition: SystemOfUnits.h:87
CLHEP::pi
static constexpr double pi
Definition: SystemOfUnits.h:55

References fact, G4DeexPrecoParameters::GetR0(), CLHEP::millibarn, CLHEP::pi, and G4VPreCompoundFragment::theParameters.

◆ ~G4PreCompoundIon()

G4PreCompoundIon::~G4PreCompoundIon ( )
virtual

Definition at line 56 of file G4PreCompoundIon.cc.

57{}

◆ G4PreCompoundIon() [2/2]

G4PreCompoundIon::G4PreCompoundIon ( const G4PreCompoundIon &  right)
delete

Member Function Documentation

◆ CalcEmissionProbability()

G4double G4PreCompoundFragment::CalcEmissionProbability ( const G4Fragment &  aFragment)
overridevirtualinherited

Implements G4VPreCompoundFragment.

Definition at line 55 of file G4PreCompoundFragment.cc.

57{
58 //G4cout << theCoulombBarrier << " " << GetMaximalKineticEnergy() << G4endl;
59 // If theCoulombBarrier effect is included in the emission probabilities
60 // Coulomb barrier is the lower limit of integration over kinetic energy
61
62 theEmissionProbability = 0.0;
63
64 if (theMaxKinEnergy <= theMinKinEnergy) { return 0.0; }
65
66 // compute power once
67 if(0 < index) {
68 muu = G4KalbachCrossSection::ComputePowerParameter(theResA, index);
69 }
70
71 theEmissionProbability =
72 IntegrateEmissionProbability(theMinKinEnergy,theMaxKinEnergy,aFragment);
73 /*
74 G4cout << "## G4PreCompoundFragment::CalcEmisProb "
75 << "Z= " << aFragment.GetZ_asInt()
76 << " A= " << aFragment.GetA_asInt()
77 << " Elow= " << LowerLimit/MeV
78 << " Eup= " << UpperLimit/MeV
79 << " prob= " << theEmissionProbability
80 << G4endl;
81 */
82 return theEmissionProbability;
83}
G4KalbachCrossSection::ComputePowerParameter
static G4double ComputePowerParameter(G4int resA, G4int idx)
Definition: G4KalbachCrossSection.cc:69
G4PreCompoundFragment::IntegrateEmissionProbability
G4double IntegrateEmissionProbability(G4double low, G4double up, const G4Fragment &aFragment)
Definition: G4PreCompoundFragment.cc:86

References G4KalbachCrossSection::ComputePowerParameter(), G4PreCompoundFragment::index, G4PreCompoundFragment::IntegrateEmissionProbability(), G4PreCompoundFragment::muu, G4VPreCompoundFragment::theEmissionProbability, G4VPreCompoundFragment::theMaxKinEnergy, G4VPreCompoundFragment::theMinKinEnergy, and G4VPreCompoundFragment::theResA.

◆ CoalescenceFactor()

virtual G4double G4PreCompoundIon::CoalescenceFactor ( G4int  A) const
protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

◆ CrossSection()

G4double G4PreCompoundFragment::CrossSection ( G4double  ekin) const
protectedinherited

Definition at line 113 of file G4PreCompoundFragment.cc.

114{
115 G4double res;
116 if(OPTxs == 0 || (OPTxs == 4 && theMaxKinEnergy < 10.)) {
117 res = GetOpt0(ekin);
118
119 } else if(OPTxs <= 2) {
120 res = G4ChatterjeeCrossSection::ComputeCrossSection(ekin,
121 theCoulombBarrier,
122 theResA13, muu,
123 index, theZ, theResA);
124
125 } else {
126 res = G4KalbachCrossSection::ComputeCrossSection(ekin, theCoulombBarrier,
127 theResA13, muu, index,
128 theZ, theA, theResA);
129 }
130 return res;
131}
G4ChatterjeeCrossSection::ComputeCrossSection
static G4double ComputeCrossSection(G4double K, G4double cb, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int resA)
Definition: G4ChatterjeeCrossSection.cc:56
G4KalbachCrossSection::ComputeCrossSection
static G4double ComputeCrossSection(G4double K, G4double cb, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int A, G4int resA)
Definition: G4KalbachCrossSection.cc:75
G4PreCompoundFragment::GetOpt0
G4double GetOpt0(G4double ekin) const
Definition: G4PreCompoundFragment.cc:133

References G4KalbachCrossSection::ComputeCrossSection(), G4ChatterjeeCrossSection::ComputeCrossSection(), G4PreCompoundFragment::GetOpt0(), G4PreCompoundFragment::index, G4PreCompoundFragment::muu, G4VPreCompoundFragment::OPTxs, G4VPreCompoundFragment::theA, G4VPreCompoundFragment::theCoulombBarrier, G4VPreCompoundFragment::theMaxKinEnergy, G4VPreCompoundFragment::theResA, G4VPreCompoundFragment::theResA13, and G4VPreCompoundFragment::theZ.

Referenced by ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

◆ FactorialFactor()

virtual G4double G4PreCompoundIon::FactorialFactor ( G4int  N,
G4int  P 
) const
protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

◆ GetA()

G4int G4VPreCompoundFragment::GetA ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::PerformEmission(), and ProbabilityDistributionFunction().

◆ GetAlpha()

virtual G4double G4VPreCompoundFragment::GetAlpha ( ) const
protectedpure virtualinherited

Implemented in G4HETCAlpha, G4HETCDeuteron, G4HETCHe3, G4HETCNeutron, G4HETCProton, G4HETCTriton, G4HETCFragment, G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, G4PreCompoundNeutron, G4PreCompoundProton, and G4PreCompoundTriton.

◆ GetBeta()

G4double G4PreCompoundIon::GetBeta ( ) const
overrideprotectedvirtual

Implements G4VPreCompoundFragment.

Definition at line 101 of file G4PreCompoundIon.cc.

102{
103 return -theCoulombBarrier;
104}

References G4VPreCompoundFragment::theCoulombBarrier.

◆ GetBindingEnergy()

G4double G4VPreCompoundFragment::GetBindingEnergy ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::AngularDistribution().

◆ GetEmissionProbability()

G4double G4VPreCompoundFragment::GetEmissionProbability ( ) const
inlineinherited

◆ GetEnergyThreshold()

G4double G4VPreCompoundFragment::GetEnergyThreshold ( ) const
inlineinherited

Referenced by G4HETCFragment::CalcEmissionProbability().

◆ GetMomentum()

const G4LorentzVector & G4VPreCompoundFragment::GetMomentum ( ) const
inlineinherited

◆ GetNuclearMass()

G4double G4VPreCompoundFragment::GetNuclearMass ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::AngularDistribution(), and G4PreCompoundEmission::PerformEmission().

◆ GetOpt0()

G4double G4PreCompoundFragment::GetOpt0 ( G4double  ekin) const
privateinherited

Definition at line 133 of file G4PreCompoundFragment.cc.

135{
136 G4double r0 = theParameters->GetR0()*theResA13;
137 // cross section is now given in mb (r0 is in mm) for the sake of consistency
138 //with the rest of the options
139 return 1.e+25*CLHEP::pi*r0*r0*theResA13*GetAlpha()*(1.0 + GetBeta()/ekin);
140}
G4VPreCompoundFragment::GetBeta
virtual G4double GetBeta() const =0
G4VPreCompoundFragment::GetAlpha
virtual G4double GetAlpha() const =0

References CLHEP::pi.

Referenced by G4PreCompoundFragment::CrossSection().

◆ GetReactionProduct()

G4ReactionProduct * G4VPreCompoundFragment::GetReactionProduct ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::PerformEmission().

◆ GetRestA()

G4int G4VPreCompoundFragment::GetRestA ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::PerformEmission().

◆ GetRestNuclearMass()

G4double G4VPreCompoundFragment::GetRestNuclearMass ( ) const
inlineinherited

◆ GetRestZ()

G4int G4VPreCompoundFragment::GetRestZ ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::PerformEmission().

◆ GetRj()

virtual G4double G4PreCompoundIon::GetRj ( G4int  NumberParticles,
G4int  NumberCharged 
) const
protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

◆ GetZ()

G4int G4VPreCompoundFragment::GetZ ( ) const
inlineinherited

Referenced by G4PreCompoundEmission::PerformEmission().

◆ Initialize()

void G4VPreCompoundFragment::Initialize ( const G4Fragment &  aFragment)
inherited

Definition at line 79 of file G4VPreCompoundFragment.cc.

80{
81 theFragA = aFragment.GetA_asInt();
82 theFragZ = aFragment.GetZ_asInt();
83 theResA = theFragA - theA;
84 theResZ = theFragZ - theZ;
85
86 theMinKinEnergy = theMaxKinEnergy = theCoulombBarrier = 0.0;
87 if ((theResA < theResZ) || (theResA < theA) || (theResZ < theZ)) {
88 return;
89 }
90
91 theResA13 = g4calc->Z13(theResA);
92 theCoulombBarrier = theCoulombBarrierPtr->
93 GetCoulombBarrier(theResA,theResZ,aFragment.GetExcitationEnergy());
94
95 G4double elim = (0 == OPTxs) ? theCoulombBarrier : theCoulombBarrier*0.5;
96
97 // Calculate masses
98 theResMass = G4NucleiProperties::GetNuclearMass(theResA, theResZ);
99 theReducedMass = theResMass*theMass/(theResMass + theMass);
100
101 // Compute Binding Energies for fragments
102 // needed to separate a fragment from the nucleus
103 theBindingEnergy = theResMass + theMass - aFragment.GetGroundStateMass();
104
105 // Compute Maximal Kinetic Energy which can be carried by fragments
106 // after separation - the true assimptotic value
107 G4double Ecm = aFragment.GetMomentum().m();
108 G4double twoEcm = Ecm + Ecm;
109 theMaxKinEnergy = std::max(((Ecm-theResMass)*(Ecm+theResMass) + theMass*theMass)
110 /twoEcm - theMass,0.0);
111 theMinKinEnergy = (elim == 0.0) ? 0.0 :
112 std::max(((theMass+elim)*(twoEcm-theMass-elim) +
113 theMass*theMass)/twoEcm - theMass,0.0);
114}
G4Fragment::GetGroundStateMass
G4double GetGroundStateMass() const
Definition: G4Fragment.hh:304
G4Fragment::GetExcitationEnergy
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:299
G4Fragment::GetMomentum
const G4LorentzVector & GetMomentum() const
Definition: G4Fragment.hh:323
G4Fragment::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Fragment.hh:276
G4Fragment::GetA_asInt
G4int GetA_asInt() const
Definition: G4Fragment.hh:271
G4NucleiProperties::GetNuclearMass
static G4double GetNuclearMass(const G4double A, const G4double Z)
Definition: G4NucleiProperties.cc:50
G4Pow::Z13
G4double Z13(G4int Z) const
Definition: G4Pow.hh:123
G4VPreCompoundFragment::theCoulombBarrierPtr
G4VCoulombBarrier * theCoulombBarrierPtr
Definition: G4VPreCompoundFragment.hh:132
G4INCL::Math::max
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Definition: G4INCLGlobals.hh:112

References G4VPreCompoundFragment::g4calc, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4Fragment::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4Fragment::GetZ_asInt(), CLHEP::HepLorentzVector::m(), G4INCL::Math::max(), G4VPreCompoundFragment::OPTxs, G4VPreCompoundFragment::theA, G4VPreCompoundFragment::theBindingEnergy, G4VPreCompoundFragment::theCoulombBarrier, G4VPreCompoundFragment::theCoulombBarrierPtr, G4VPreCompoundFragment::theFragA, G4VPreCompoundFragment::theFragZ, G4VPreCompoundFragment::theMass, G4VPreCompoundFragment::theMaxKinEnergy, G4VPreCompoundFragment::theMinKinEnergy, G4VPreCompoundFragment::theReducedMass, G4VPreCompoundFragment::theResA, G4VPreCompoundFragment::theResA13, G4VPreCompoundFragment::theResMass, G4VPreCompoundFragment::theResZ, G4VPreCompoundFragment::theZ, and G4Pow::Z13().

◆ IntegrateEmissionProbability()

G4double G4PreCompoundFragment::IntegrateEmissionProbability ( G4double  low,
G4double  up,
const G4Fragment &  aFragment 
)
privateinherited

Definition at line 85 of file G4PreCompoundFragment.cc.

88{
89 static const G4double den = 1.0/CLHEP::MeV;
90 G4double del = (up - low);
91 G4int nbins = del*den;
92 nbins = std::max(nbins, 4);
93 del /= static_cast<G4double>(nbins);
94 G4double e = low + 0.5*del;
95 probmax = ProbabilityDistributionFunction(e, aFragment);
96 //G4cout << " 0. e= " << e << " y= " << probmax << G4endl;
97
98 G4double sum = probmax;
99 for (G4int i=1; i<nbins; ++i) {
100 e += del;
101
102 G4double y = ProbabilityDistributionFunction(e, aFragment);
103 probmax = std::max(probmax, y);
104 sum += y;
105 if(y < sum*0.01) { break; }
106 //G4cout <<" "<<i<<". e= "<<e<<" y= "<<y<<" sum= "<<sum<< G4endl;
107 }
108 sum *= del;
109 //G4cout << "Evap prob: " << sum << " probmax= " << probmax << G4endl;
110 return sum;
111}
G4int
int G4int
Definition: G4Types.hh:85
G4PreCompoundFragment::ProbabilityDistributionFunction
virtual G4double ProbabilityDistributionFunction(G4double ekin, const G4Fragment &aFragment)=0
CLHEP::MeV
static constexpr double MeV
Definition: SystemOfUnits.h:181

References G4INCL::Math::max(), CLHEP::MeV, G4PreCompoundFragment::ProbabilityDistributionFunction(), and G4PreCompoundFragment::probmax.

Referenced by G4PreCompoundFragment::CalcEmissionProbability().

◆ IsItPossible()

G4bool G4VPreCompoundFragment::IsItPossible ( const G4Fragment &  aFragment) const
inlineinherited

◆ operator!=() [1/3]

G4bool G4PreCompoundFragment::operator!= ( const G4PreCompoundFragment &  right) const
deleteinherited

◆ operator!=() [2/3]

G4bool G4PreCompoundIon::operator!= ( const G4PreCompoundIon &  right) const
delete

◆ operator!=() [3/3]

G4bool G4VPreCompoundFragment::operator!= ( const G4VPreCompoundFragment &  right) const
deleteinherited

◆ operator=()

const G4PreCompoundIon & G4PreCompoundIon::operator= ( const G4PreCompoundIon &  right)
delete

◆ operator==() [1/3]

G4bool G4PreCompoundFragment::operator== ( const G4PreCompoundFragment &  right) const
deleteinherited

◆ operator==() [2/3]

G4bool G4PreCompoundIon::operator== ( const G4PreCompoundIon &  right) const
delete

◆ operator==() [3/3]

G4bool G4VPreCompoundFragment::operator== ( const G4VPreCompoundFragment &  right) const
deleteinherited

◆ ProbabilityDistributionFunction()

G4double G4PreCompoundIon::ProbabilityDistributionFunction ( G4double  eKin,
const G4Fragment &  aFragment 
)
overrideprotectedvirtual

Implements G4PreCompoundFragment.

Definition at line 59 of file G4PreCompoundIon.cc.

62{
63 G4double efinal = eKin + theBindingEnergy;
64 if(efinal <= 0.0 ) { return 0.0; }
65
66 G4double U = aFragment.GetExcitationEnergy();
67 G4int P = aFragment.GetNumberOfParticles();
68 G4int H = aFragment.GetNumberOfHoles();
69 G4int A = GetA();
70 G4int N = P + H;
71
72 static const G4double sixoverpi2 = 6.0/CLHEP::pi2;
73 G4double g0 = sixoverpi2*fNucData->GetLevelDensity(theFragZ, theFragA, U);
74 G4double g1 = sixoverpi2*fNucData->GetLevelDensity(theResZ, theResA, 0.0);
75
76 G4double gj = g1;
77
78 G4double A0 = (P*P+H*H+P-3*H)/(4.0*g0);
79 G4double A1 = std::max(0.0,(A0*g0 + A*(A-2*P-1)*0.25)/g1);
80
81 G4double E0 = U - A0;
82 if (E0 <= 0.0) { return 0.0; }
83
84 G4double E1 = std::max(0.0,theMaxKinEnergy - eKin - A1);
85
86 G4double Aj = A*(A+1)/(4.0*gj);
87 G4double Ej = std::max(0.0,efinal - Aj);
88
89 G4double rj = GetRj(P, aFragment.GetNumberOfCharged());
90 G4double xs = CrossSection(eKin);
91
92 G4double pA = fact*eKin*xs*rj
93 * CoalescenceFactor(theFragA) * FactorialFactor(N,P)
94 * std::sqrt(2.0/(theReducedMass*efinal))
95 * g4calc->powN(g1*E1/(g0*E0), N-A-1)
96 * g4calc->powN(gj*Ej/(g0*E0), A-1)*gj*g1/(g0*g0*E0*theResA);
97
98 return pA;
99}
A
const G4double A[17]
Definition: G4WaterStopping.cc:53
G4Fragment::GetNumberOfParticles
G4int GetNumberOfParticles() const
Definition: G4Fragment.hh:361
G4Fragment::GetNumberOfHoles
G4int GetNumberOfHoles() const
Definition: G4Fragment.hh:381
G4Fragment::GetNumberOfCharged
G4int GetNumberOfCharged() const
Definition: G4Fragment.hh:366
G4NuclearLevelData::GetLevelDensity
G4double GetLevelDensity(G4int Z, G4int A, G4double U)
Definition: G4NuclearLevelData.cc:634
G4Pow::powN
G4double powN(G4double x, G4int n) const
Definition: G4Pow.cc:166
G4PreCompoundFragment::CrossSection
G4double CrossSection(G4double ekin) const
Definition: G4PreCompoundFragment.cc:113
G4PreCompoundIon::FactorialFactor
virtual G4double FactorialFactor(G4int N, G4int P) const =0
G4PreCompoundIon::CoalescenceFactor
virtual G4double CoalescenceFactor(G4int A) const =0
G4PreCompoundIon::GetRj
virtual G4double GetRj(G4int NumberParticles, G4int NumberCharged) const =0
G4VPreCompoundFragment::fNucData
G4NuclearLevelData * fNucData
Definition: G4VPreCompoundFragment.hh:138
CLHEP::pi2
static constexpr double pi2
Definition: SystemOfUnits.h:58
P
static double P[]
Definition: nf_gammaFunctions.cc:89

References A, CoalescenceFactor(), G4PreCompoundFragment::CrossSection(), fact, FactorialFactor(), G4VPreCompoundFragment::fNucData, G4VPreCompoundFragment::g4calc, G4VPreCompoundFragment::GetA(), G4Fragment::GetExcitationEnergy(), G4NuclearLevelData::GetLevelDensity(), G4Fragment::GetNumberOfCharged(), G4Fragment::GetNumberOfHoles(), G4Fragment::GetNumberOfParticles(), GetRj(), G4INCL::Math::max(), P, CLHEP::pi2, G4Pow::powN(), G4VPreCompoundFragment::theBindingEnergy, G4VPreCompoundFragment::theFragA, G4VPreCompoundFragment::theFragZ, G4VPreCompoundFragment::theMaxKinEnergy, G4VPreCompoundFragment::theReducedMass, G4VPreCompoundFragment::theResA, and G4VPreCompoundFragment::theResZ.

◆ SampleKineticEnergy()

G4double G4PreCompoundFragment::SampleKineticEnergy ( const G4Fragment &  aFragment)
overridevirtualinherited

Implements G4VPreCompoundFragment.

Definition at line 142 of file G4PreCompoundFragment.cc.

143{
144 G4double delta = theMaxKinEnergy - theMinKinEnergy;
145 static const G4double toler = 1.25;
146 probmax *= toler;
147 G4double prob, T(0.0);
148 CLHEP::HepRandomEngine* rndm = G4Random::getTheEngine();
149 G4int i;
150 for(i=0; i<100; ++i) {
151 T = theMinKinEnergy + delta*rndm->flat();
152 prob = ProbabilityDistributionFunction(T, fragment);
153 /*
154 if(prob > probmax) {
155 G4cout << "G4PreCompoundFragment WARNING: prob= " << prob
156 << " probmax= " << probmax << G4endl;
157 G4cout << "i= " << i << " Z= " << theZ << " A= " << theA
158 << " resZ= " << theResZ << " resA= " << theResA << "\n"
159 << " T= " << T << " Tmax= " << theMaxKinEnergy
160 << " Tmin= " << limit
161 << G4endl;
162 for(G4int i=0; i<N; ++i) { G4cout << " " << probability[i]; }
163 G4cout << G4endl;
164 }
165 */
166 // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
167 if(probmax*rndm->flat() <= prob) { break; }
168 }
169 /*
170 G4cout << "G4PreCompoundFragment: i= " << i << " Z= " << theZ
171 << " A= " << theA <<" T(MeV)= " << T << " Emin(MeV)= "
172 << theMinKinEnergy << " Emax= " << theMaxKinEnergy << G4endl;
173 */
174 return T;
175}
CLHEP::HepRandomEngine::flat
virtual double flat()=0
anonymous_namespace{G4QuasiElRatios.cc}::toler
const G4double toler
Definition: G4QuasiElRatios.cc:69

References CLHEP::HepRandomEngine::flat(), G4PreCompoundFragment::ProbabilityDistributionFunction(), G4PreCompoundFragment::probmax, G4VPreCompoundFragment::theMaxKinEnergy, G4VPreCompoundFragment::theMinKinEnergy, and anonymous_namespace{G4QuasiElRatios.cc}::toler.

◆ SetMomentum()

void G4VPreCompoundFragment::SetMomentum ( const G4LorentzVector &  value)
inlineinherited

Referenced by G4PreCompoundEmission::PerformEmission().

◆ SetOPTxs()

void G4VPreCompoundFragment::SetOPTxs ( G4int  )
inlineinherited

◆ UseSICB()

void G4VPreCompoundFragment::UseSICB ( G4bool  )
inlineinherited

Field Documentation

◆ fact

G4double G4PreCompoundIon::fact
private

Definition at line 71 of file G4PreCompoundIon.hh.

Referenced by G4PreCompoundIon(), and ProbabilityDistributionFunction().

◆ fNucData

G4NuclearLevelData* G4VPreCompoundFragment::fNucData
protectedinherited

Definition at line 138 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::G4VPreCompoundFragment(), G4HETCFragment::IntegrateEmissionProbability(), ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

◆ g4calc

G4Pow* G4VPreCompoundFragment::g4calc
protectedinherited

Definition at line 140 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::G4VPreCompoundFragment(), G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

◆ index

G4int G4PreCompoundFragment::index
privateinherited

Definition at line 84 of file G4PreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CalcEmissionProbability(), G4PreCompoundFragment::CrossSection(), and G4PreCompoundFragment::G4PreCompoundFragment().

◆ muu

G4double G4PreCompoundFragment::muu
privateinherited

Definition at line 86 of file G4PreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CalcEmissionProbability(), G4PreCompoundFragment::CrossSection(), and G4PreCompoundFragment::G4PreCompoundFragment().

◆ OPTxs

G4int G4VPreCompoundFragment::OPTxs
protectedinherited

Definition at line 161 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), and G4VPreCompoundFragment::Initialize().

◆ particle

const G4ParticleDefinition* G4VPreCompoundFragment::particle
privateinherited

Definition at line 131 of file G4VPreCompoundFragment.hh.

Referenced by emcalc_gui.MainWindow::cb_select_particle(), emcalc_gui.MainWindow::cb_show_plot(), and G4VPreCompoundFragment::G4VPreCompoundFragment().

◆ probmax

G4double G4PreCompoundFragment::probmax
privateinherited

Definition at line 87 of file G4PreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::G4PreCompoundFragment(), G4PreCompoundFragment::IntegrateEmissionProbability(), and G4PreCompoundFragment::SampleKineticEnergy().

◆ theA

G4int G4VPreCompoundFragment::theA
protectedinherited

Definition at line 142 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), G4PreCompoundFragment::G4PreCompoundFragment(), G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), and G4HETCTriton::K().

◆ theBindingEnergy

G4double G4VPreCompoundFragment::theBindingEnergy
protectedinherited

Definition at line 150 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::Initialize(), ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

◆ theCoulombBarrier

G4double G4VPreCompoundFragment::theCoulombBarrier
protectedinherited

Definition at line 158 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCFragment::CalcEmissionProbability(), G4PreCompoundFragment::CrossSection(), G4HETCAlpha::GetBeta(), G4HETCDeuteron::GetBeta(), G4HETCHe3::GetBeta(), G4HETCProton::GetBeta(), G4HETCTriton::GetBeta(), GetBeta(), G4PreCompoundProton::GetBeta(), G4VPreCompoundFragment::Initialize(), and G4HETCChargedFragment::SampleKineticEnergy().

◆ theCoulombBarrierPtr

G4VCoulombBarrier* G4VPreCompoundFragment::theCoulombBarrierPtr
privateinherited

Definition at line 132 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::Initialize().

◆ theEmissionProbability

G4double G4VPreCompoundFragment::theEmissionProbability
protectedinherited

Definition at line 157 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCFragment::CalcEmissionProbability(), and G4PreCompoundFragment::CalcEmissionProbability().

◆ theFragA

G4int G4VPreCompoundFragment::theFragA
protectedinherited

Definition at line 146 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

◆ theFragZ

G4int G4VPreCompoundFragment::theFragZ
protectedinherited

Definition at line 147 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCAlpha::GetAlpha(), G4HETCDeuteron::GetAlpha(), G4HETCHe3::GetAlpha(), G4HETCTriton::GetAlpha(), G4PreCompoundAlpha::GetAlpha(), G4PreCompoundDeuteron::GetAlpha(), G4PreCompoundHe3::GetAlpha(), G4PreCompoundTriton::GetAlpha(), G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

◆ theMass

G4double G4VPreCompoundFragment::theMass
protectedinherited

Definition at line 155 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::G4VPreCompoundFragment(), and G4VPreCompoundFragment::Initialize().

◆ theMaxKinEnergy

G4double G4VPreCompoundFragment::theMaxKinEnergy
protectedinherited

Definition at line 152 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCFragment::CalcEmissionProbability(), G4PreCompoundFragment::CalcEmissionProbability(), G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment::Initialize(), ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), G4HETCNeutron::SampleKineticEnergy(), and G4PreCompoundFragment::SampleKineticEnergy().

◆ theMinKinEnergy

G4double G4VPreCompoundFragment::theMinKinEnergy
protectedinherited

Definition at line 151 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CalcEmissionProbability(), G4VPreCompoundFragment::Initialize(), and G4PreCompoundFragment::SampleKineticEnergy().

◆ theMomentum

G4LorentzVector G4VPreCompoundFragment::theMomentum
privateinherited

Definition at line 134 of file G4VPreCompoundFragment.hh.

◆ theParameters

G4DeexPrecoParameters* G4VPreCompoundFragment::theParameters
protectedinherited

Definition at line 139 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCFragment::G4HETCFragment(), G4PreCompoundIon(), and G4VPreCompoundFragment::G4VPreCompoundFragment().

◆ theReducedMass

G4double G4VPreCompoundFragment::theReducedMass
protectedinherited

Definition at line 154 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

◆ theResA

G4int G4VPreCompoundFragment::theResA
protectedinherited

Definition at line 144 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CalcEmissionProbability(), G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), G4HETCProton::K(), G4HETCTriton::K(), ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

◆ theResA13

G4double G4VPreCompoundFragment::theResA13
protectedinherited

Definition at line 149 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4HETCNeutron::GetAlpha(), G4PreCompoundNeutron::GetAlpha(), G4HETCNeutron::GetBeta(), G4PreCompoundNeutron::GetBeta(), and G4VPreCompoundFragment::Initialize().

◆ theResMass

G4double G4VPreCompoundFragment::theResMass
protectedinherited

Definition at line 153 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment::Initialize().

◆ theResZ

G4int G4VPreCompoundFragment::theResZ
protectedinherited

Definition at line 145 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCProton::GetAlpha(), G4PreCompoundProton::GetAlpha(), G4VPreCompoundFragment::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), G4HETCProton::K(), G4HETCTriton::K(), ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

◆ theZ

G4int G4VPreCompoundFragment::theZ
protectedinherited

Definition at line 143 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), G4PreCompoundFragment::G4PreCompoundFragment(), G4VPreCompoundFragment::Initialize(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), G4HETCProton::K(), and G4HETCTriton::K().

◆ useSICB

G4bool G4VPreCompoundFragment::useSICB
protectedinherited

Definition at line 163 of file G4VPreCompoundFragment.hh.

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3