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

#include <G4VPreCompoundFragment.hh>

Inheritance diagram for G4VPreCompoundFragment:
G4HETCFragment G4PreCompoundFragment G4HETCChargedFragment G4HETCNeutron G4PreCompoundIon G4PreCompoundNucleon G4HETCAlpha G4HETCDeuteron G4HETCHe3 G4HETCProton G4HETCTriton G4PreCompoundAlpha G4PreCompoundDeuteron G4PreCompoundHe3 G4PreCompoundTriton G4PreCompoundNeutron G4PreCompoundProton

Public Member Functions

virtual G4double CalcEmissionProbability (const G4Fragment &aFragment)=0
 
 G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
 
 G4VPreCompoundFragment (const G4VPreCompoundFragment &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 G4VPreCompoundFragment &right) const =delete
 
const G4VPreCompoundFragment & operator= (const G4VPreCompoundFragment &right)=delete
 
G4bool operator== (const G4VPreCompoundFragment &right) const =delete
 
virtual G4double SampleKineticEnergy (const G4Fragment &aFragment)=0
 
void SetMomentum (const G4LorentzVector &value)
 
void SetOPTxs (G4int)
 
void UseSICB (G4bool)
 
virtual ~G4VPreCompoundFragment ()
 

Protected Member Functions

virtual G4double GetAlpha () const =0
 
virtual G4double GetBeta () const =0
 

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 Attributes

const G4ParticleDefinition * particle
 
G4VCoulombBarrier * theCoulombBarrierPtr
 
G4LorentzVector theMomentum
 

Friends

std::ostream & operator<< (std::ostream &, const G4VPreCompoundFragment &)
 
std::ostream & operator<< (std::ostream &, const G4VPreCompoundFragment *)
 

Detailed Description

Definition at line 55 of file G4VPreCompoundFragment.hh.

Constructor & Destructor Documentation

◆ G4VPreCompoundFragment() [1/2]

G4VPreCompoundFragment::G4VPreCompoundFragment ( const G4ParticleDefinition *  part,
G4VCoulombBarrier *  aCoulombBarrier 
)
explicit

Definition at line 39 of file G4VPreCompoundFragment.cc.

41 : particle(part), theCoulombBarrierPtr(aCoulombBarrier),
42 theMomentum(0.,0.,0.,0.),
43 theA(particle->GetBaryonNumber()),
44 theZ(G4lrint(particle->GetPDGCharge())),
45 theResA(0),theResZ(0),theFragA(0),theFragZ(0),theBindingEnergy(0.0),
46 theMinKinEnergy(0.0),theMaxKinEnergy(0.0),theResMass(0.0),
47 theReducedMass(0.0),
48 theEmissionProbability(0.0),theCoulombBarrier(0.0),
49 OPTxs(3),useSICB(true)
50{
51 theMass = particle->GetPDGMass();
52 fNucData = G4NuclearLevelData::GetInstance();
53 theParameters = fNucData->GetParameters();
54 g4calc = G4Pow::GetInstance();
55 theResA13 = 0.0;
56}
G4NuclearLevelData::GetParameters
G4DeexPrecoParameters * GetParameters()
Definition: G4NuclearLevelData.cc:619
G4NuclearLevelData::GetInstance
static G4NuclearLevelData * GetInstance()
Definition: G4NuclearLevelData.cc:421
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:113
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:41
G4VPreCompoundFragment::fNucData
G4NuclearLevelData * fNucData
Definition: G4VPreCompoundFragment.hh:138
G4VPreCompoundFragment::theCoulombBarrierPtr
G4VCoulombBarrier * theCoulombBarrierPtr
Definition: G4VPreCompoundFragment.hh:132
G4VPreCompoundFragment::theParameters
G4DeexPrecoParameters * theParameters
Definition: G4VPreCompoundFragment.hh:139
G4VPreCompoundFragment::particle
const G4ParticleDefinition * particle
Definition: G4VPreCompoundFragment.hh:131
G4lrint
int G4lrint(double ad)
Definition: templates.hh:134

References fNucData, g4calc, G4Pow::GetInstance(), G4NuclearLevelData::GetInstance(), G4NuclearLevelData::GetParameters(), G4ParticleDefinition::GetPDGMass(), particle, theMass, theParameters, and theResA13.

◆ ~G4VPreCompoundFragment()

G4VPreCompoundFragment::~G4VPreCompoundFragment ( )
virtual

Definition at line 58 of file G4VPreCompoundFragment.cc.

59{}

◆ G4VPreCompoundFragment() [2/2]

G4VPreCompoundFragment::G4VPreCompoundFragment ( const G4VPreCompoundFragment &  right)
delete

Member Function Documentation

◆ CalcEmissionProbability()

virtual G4double G4VPreCompoundFragment::CalcEmissionProbability ( const G4Fragment &  aFragment)
pure virtual

Implemented in G4HETCFragment, and G4PreCompoundFragment.

◆ GetA()

G4int G4VPreCompoundFragment::GetA ( ) const
inline

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

◆ GetAlpha()

virtual G4double G4VPreCompoundFragment::GetAlpha ( ) const
protectedpure virtual

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

◆ GetBeta()

virtual G4double G4VPreCompoundFragment::GetBeta ( ) const
protectedpure virtual

Implemented in G4HETCAlpha, G4HETCDeuteron, G4HETCHe3, G4HETCNeutron, G4HETCProton, G4HETCTriton, G4HETCFragment, G4PreCompoundIon, G4PreCompoundNeutron, and G4PreCompoundProton.

◆ GetBindingEnergy()

G4double G4VPreCompoundFragment::GetBindingEnergy ( ) const
inline

Referenced by G4PreCompoundEmission::AngularDistribution().

◆ GetEmissionProbability()

G4double G4VPreCompoundFragment::GetEmissionProbability ( ) const
inline

◆ GetEnergyThreshold()

G4double G4VPreCompoundFragment::GetEnergyThreshold ( ) const
inline

Referenced by G4HETCFragment::CalcEmissionProbability().

◆ GetMomentum()

const G4LorentzVector & G4VPreCompoundFragment::GetMomentum ( ) const
inline

◆ GetNuclearMass()

G4double G4VPreCompoundFragment::GetNuclearMass ( ) const
inline

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

◆ GetReactionProduct()

G4ReactionProduct * G4VPreCompoundFragment::GetReactionProduct ( ) const
inline

Referenced by G4PreCompoundEmission::PerformEmission().

◆ GetRestA()

G4int G4VPreCompoundFragment::GetRestA ( ) const
inline

Referenced by G4PreCompoundEmission::PerformEmission().

◆ GetRestNuclearMass()

G4double G4VPreCompoundFragment::GetRestNuclearMass ( ) const
inline

◆ GetRestZ()

G4int G4VPreCompoundFragment::GetRestZ ( ) const
inline

Referenced by G4PreCompoundEmission::PerformEmission().

◆ GetZ()

G4int G4VPreCompoundFragment::GetZ ( ) const
inline

Referenced by G4PreCompoundEmission::PerformEmission().

◆ Initialize()

void G4VPreCompoundFragment::Initialize ( const G4Fragment &  aFragment)

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}
G4double
double G4double
Definition: G4Types.hh:83
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
G4INCL::Math::max
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Definition: G4INCLGlobals.hh:112

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

◆ IsItPossible()

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

◆ operator!=()

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

◆ operator=()

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

◆ operator==()

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

◆ SampleKineticEnergy()

virtual G4double G4VPreCompoundFragment::SampleKineticEnergy ( const G4Fragment &  aFragment)
pure virtual

Implemented in G4HETCChargedFragment, G4HETCNeutron, and G4PreCompoundFragment.

Referenced by G4PreCompoundEmission::PerformEmission().

◆ SetMomentum()

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

Referenced by G4PreCompoundEmission::PerformEmission().

◆ SetOPTxs()

void G4VPreCompoundFragment::SetOPTxs ( G4int  )
inline

◆ UseSICB()

void G4VPreCompoundFragment::UseSICB ( G4bool  )
inline

Friends And Related Function Documentation

◆ operator<< [1/2]

std::ostream & operator<< ( std::ostream &  out,
const G4VPreCompoundFragment &  theFragment 
)
friend

Definition at line 61 of file G4VPreCompoundFragment.cc.

63{
64 out << &theFragment;
65 return out;
66}

◆ operator<< [2/2]

std::ostream & operator<< ( std::ostream &  out,
const G4VPreCompoundFragment *  theFragment 
)
friend

Definition at line 68 of file G4VPreCompoundFragment.cc.

70{
71 out
72 << "PreCompoundModel Emitted Fragment: Z= " << theFragment->GetZ()
73 << " A= " << theFragment->GetA()
74 << " Mass(GeV)= " << theFragment->GetNuclearMass()/CLHEP::GeV;
75 return out;
76}
G4VPreCompoundFragment::GetNuclearMass
G4double GetNuclearMass() const
CLHEP::GeV
static constexpr double GeV
Definition: SystemOfUnits.h:184

Field Documentation

◆ fNucData

G4NuclearLevelData* G4VPreCompoundFragment::fNucData
protected

Definition at line 138 of file G4VPreCompoundFragment.hh.

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

◆ g4calc

G4Pow* G4VPreCompoundFragment::g4calc
protected

Definition at line 140 of file G4VPreCompoundFragment.hh.

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

◆ OPTxs

G4int G4VPreCompoundFragment::OPTxs
protected

Definition at line 161 of file G4VPreCompoundFragment.hh.

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

◆ particle

const G4ParticleDefinition* G4VPreCompoundFragment::particle
private

Definition at line 131 of file G4VPreCompoundFragment.hh.

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

◆ theA

G4int G4VPreCompoundFragment::theA
protected

Definition at line 142 of file G4VPreCompoundFragment.hh.

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

◆ theBindingEnergy

G4double G4VPreCompoundFragment::theBindingEnergy
protected

Definition at line 150 of file G4VPreCompoundFragment.hh.

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

◆ theCoulombBarrier

G4double G4VPreCompoundFragment::theCoulombBarrier
protected

Definition at line 158 of file G4VPreCompoundFragment.hh.

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

◆ theCoulombBarrierPtr

G4VCoulombBarrier* G4VPreCompoundFragment::theCoulombBarrierPtr
private

Definition at line 132 of file G4VPreCompoundFragment.hh.

Referenced by Initialize().

◆ theEmissionProbability

G4double G4VPreCompoundFragment::theEmissionProbability
protected

Definition at line 157 of file G4VPreCompoundFragment.hh.

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

◆ theFragA

G4int G4VPreCompoundFragment::theFragA
protected

Definition at line 146 of file G4VPreCompoundFragment.hh.

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

◆ theFragZ

G4int G4VPreCompoundFragment::theFragZ
protected

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(), Initialize(), G4HETCFragment::IntegrateEmissionProbability(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

◆ theMass

G4double G4VPreCompoundFragment::theMass
protected

Definition at line 155 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment(), and Initialize().

◆ theMaxKinEnergy

G4double G4VPreCompoundFragment::theMaxKinEnergy
protected

Definition at line 152 of file G4VPreCompoundFragment.hh.

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

◆ theMinKinEnergy

G4double G4VPreCompoundFragment::theMinKinEnergy
protected

Definition at line 151 of file G4VPreCompoundFragment.hh.

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

◆ theMomentum

G4LorentzVector G4VPreCompoundFragment::theMomentum
private

Definition at line 134 of file G4VPreCompoundFragment.hh.

◆ theParameters

G4DeexPrecoParameters* G4VPreCompoundFragment::theParameters
protected

Definition at line 139 of file G4VPreCompoundFragment.hh.

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

◆ theReducedMass

G4double G4VPreCompoundFragment::theReducedMass
protected

Definition at line 154 of file G4VPreCompoundFragment.hh.

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

◆ theResA

G4int G4VPreCompoundFragment::theResA
protected

Definition at line 144 of file G4VPreCompoundFragment.hh.

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

◆ theResA13

G4double G4VPreCompoundFragment::theResA13
protected

Definition at line 149 of file G4VPreCompoundFragment.hh.

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

◆ theResMass

G4double G4VPreCompoundFragment::theResMass
protected

Definition at line 153 of file G4VPreCompoundFragment.hh.

Referenced by Initialize().

◆ theResZ

G4int G4VPreCompoundFragment::theResZ
protected

Definition at line 145 of file G4VPreCompoundFragment.hh.

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

◆ theZ

G4int G4VPreCompoundFragment::theZ
protected

Definition at line 143 of file G4VPreCompoundFragment.hh.

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

◆ useSICB

G4bool G4VPreCompoundFragment::useSICB
protected

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