Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4BetaMinusDecay Class Reference

#include <G4BetaMinusDecay.hh>

Inheritance diagram for G4BetaMinusDecay:
Collaboration diagram for G4BetaMinusDecay:

Public Member Functions

 G4BetaMinusDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &endpointE, const G4double &ex, const G4Ions::G4FloatLevelBase &flb, const G4BetaDecayType &type)
 
virtual ~G4BetaMinusDecay ()
 
virtual G4DecayProductsDecayIt (G4double)
 
virtual void DumpNuclearInfo ()
 
- Public Member Functions inherited from G4NuclearDecay
 G4NuclearDecay (const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation, const G4Ions::G4FloatLevelBase &floatingLevel)
 
virtual ~G4NuclearDecay ()
 
G4RadioactiveDecayMode GetDecayMode ()
 
G4double GetDaughterExcitation ()
 
G4Ions::G4FloatLevelBase GetFloatingLevel ()
 
G4ParticleDefinitionGetDaughterNucleus ()
 
void SetHLThreshold (G4double HLT)
 
G4double GetHLThreshold ()
 
- Public Member Functions inherited from G4VDecayChannel
 G4VDecayChannel (const G4String &aName, G4int Verbose=1)
 
 G4VDecayChannel (const G4String &aName, const G4String &theParentName, G4double theBR, G4int theNumberOfDaughters, const G4String &theDaughterName1, const G4String &theDaughterName2="", const G4String &theDaughterName3="", const G4String &theDaughterName4="")
 
virtual ~G4VDecayChannel ()
 
G4int operator== (const G4VDecayChannel &right) const
 
G4int operator!= (const G4VDecayChannel &right) const
 
G4int operator< (const G4VDecayChannel &right) const
 
const G4StringGetKinematicsName () const
 
G4double GetBR () const
 
G4int GetNumberOfDaughters () const
 
G4ParticleDefinitionGetParent ()
 
G4ParticleDefinitionGetDaughter (G4int anIndex)
 
G4int GetAngularMomentum ()
 
const G4StringGetParentName () const
 
const G4StringGetDaughterName (G4int anIndex) const
 
G4double GetParentMass () const
 
G4double GetDaughterMass (G4int anIndex) const
 
void SetParent (const G4ParticleDefinition *particle_type)
 
void SetParent (const G4String &particle_name)
 
void SetBR (G4double value)
 
void SetNumberOfDaughters (G4int value)
 
void SetDaughter (G4int anIndex, const G4ParticleDefinition *particle_type)
 
void SetDaughter (G4int anIndex, const G4String &particle_name)
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
void DumpInfo ()
 
G4double GetRangeMass () const
 
void SetRangeMass (G4double val)
 
virtual G4bool IsOKWithParentMass (G4double parentMass)
 
void SetPolarization (const G4ThreeVector &)
 
const G4ThreeVectorGetPolarization () const
 

Additional Inherited Members

- Protected Member Functions inherited from G4VDecayChannel
void ClearDaughtersName ()
 
void CheckAndFillDaughters ()
 
void CheckAndFillParent ()
 
 G4VDecayChannel ()
 
 G4VDecayChannel (const G4VDecayChannel &)
 
G4VDecayChanneloperator= (const G4VDecayChannel &)
 
G4double DynamicalMass (G4double massPDG, G4double width, G4double maxDev=+1.) const
 
- Protected Attributes inherited from G4NuclearDecay
const G4RadioactiveDecayMode theMode
 
- Protected Attributes inherited from G4VDecayChannel
G4String kinematics_name
 
G4double rbranch
 
G4int numberOfDaughters
 
G4Stringparent_name
 
G4String ** daughters_name
 
G4double rangeMass
 
G4ThreeVector parent_polarization
 
G4ParticleTableparticletable
 
G4int verboseLevel
 
G4ParticleDefinitionG4MT_parent
 
G4ParticleDefinition ** G4MT_daughters
 
G4double G4MT_parent_mass
 
G4doubleG4MT_daughters_mass
 
G4doubleG4MT_daughters_width
 
G4Mutex daughtersMutex
 
G4Mutex parentMutex
 
- Static Protected Attributes inherited from G4VDecayChannel
static const G4String noName = " "
 

Detailed Description

Definition at line 44 of file G4BetaMinusDecay.hh.

Constructor & Destructor Documentation

G4BetaMinusDecay::G4BetaMinusDecay ( const G4ParticleDefinition theParentNucleus,
const G4double theBR,
const G4double endpointE,
const G4double ex,
const G4Ions::G4FloatLevelBase flb,
const G4BetaDecayType type 
)

Definition at line 44 of file G4BetaMinusDecay.cc.

49  : G4NuclearDecay("beta- decay", BetaMinus, excitationE, flb), endpointEnergy(e0)
50 {
51  SetParent(theParentNucleus); // Store name of parent nucleus, delete G4MT_parent
52  SetBR(branch);
53 
55  G4IonTable* theIonTable =
57  G4int daughterZ = theParentNucleus->GetAtomicNumber() + 1;
58  G4int daughterA = theParentNucleus->GetAtomicMass();
59  SetDaughter(0, theIonTable->GetIon(daughterZ, daughterA, excitationE, flb) );
60  SetDaughter(1, "e-");
61  SetDaughter(2, "anti_nu_e");
62 
63  SetUpBetaSpectrumSampler(daughterZ, daughterA, betaType);
64 }
void SetBR(G4double value)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
Definition: G4IonTable.cc:503
int G4int
Definition: G4Types.hh:78
G4int GetAtomicNumber() const
G4IonTable * GetIonTable() const
void SetNumberOfDaughters(G4int value)
G4int GetAtomicMass() const
G4NuclearDecay(const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation, const G4Ions::G4FloatLevelBase &floatingLevel)
static G4ParticleTable * GetParticleTable()
void SetParent(const G4ParticleDefinition *particle_type)
void SetDaughter(G4int anIndex, const G4ParticleDefinition *particle_type)

Here is the call graph for this function:

G4BetaMinusDecay::~G4BetaMinusDecay ( )
virtual

Definition at line 67 of file G4BetaMinusDecay.cc.

68 {
69  delete spectrumSampler;
70 }

Member Function Documentation

G4DecayProducts * G4BetaMinusDecay::DecayIt ( G4double  )
virtual

Implements G4VDecayChannel.

Definition at line 73 of file G4BetaMinusDecay.cc.

74 {
75  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)
77 
78  // Fill G4MT_daughters with e-, nu and residual nucleus (stored by SetDaughter)
80 
81  G4double parentMass = G4MT_parent->GetPDGMass();
82  G4double eMass = G4MT_daughters[1]->GetPDGMass();
83  G4double nucleusMass = G4MT_daughters[0]->GetPDGMass();
84  // Set up final state
85  // parentParticle is set at rest here because boost with correct momentum
86  // is done later
87  G4DynamicParticle parentParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
88  G4DecayProducts* products = new G4DecayProducts(parentParticle);
89 
90  if (spectrumSampler) {
91  // Electron, neutrino and daughter nucleus energies
92  G4double eKE = endpointEnergy*spectrumSampler->shoot(G4Random::getTheEngine() );
93  G4double eMomentum = std::sqrt(eKE*(eKE + 2.*eMass) );
94 
95  G4double cosThetaENu = 2.*G4UniformRand() - 1.;
96  G4double eTE = eMass + eKE;
97  G4double nuEnergy = ((endpointEnergy - eKE)*(parentMass + nucleusMass - eTE)
98  - eMomentum*eMomentum)/(parentMass - eTE + eMomentum*cosThetaENu)/2.;
99 
100  // Electron 4-vector, isotropic angular distribution
101  G4double cosTheta = 2.*G4UniformRand() - 1.0;
102  G4double sinTheta = std::sqrt(1.0 - cosTheta*cosTheta);
103 
105  G4double sinPhi = std::sin(phi);
106  G4double cosPhi = std::cos(phi);
107 
108  G4ParticleMomentum eDirection(sinTheta*cosPhi, sinTheta*sinPhi, cosTheta);
109  G4DynamicParticle* dynamicElectron
110  = new G4DynamicParticle(G4MT_daughters[1], eDirection*eMomentum);
111  products->PushProducts(dynamicElectron);
112 
113  // Neutrino 4-vector
114  G4double sinThetaENu = std::sqrt(1.0 - cosThetaENu*cosThetaENu);
115  phi = twopi*G4UniformRand()*rad;
116  G4double sinPhiNu = std::sin(phi);
117  G4double cosPhiNu = std::cos(phi);
118 
119  G4ParticleMomentum nuDirection;
120  nuDirection.setX(sinThetaENu*cosPhiNu*cosTheta*cosPhi -
121  sinThetaENu*sinPhiNu*sinPhi + cosThetaENu*sinTheta*cosPhi);
122  nuDirection.setY(sinThetaENu*cosPhiNu*cosTheta*sinPhi +
123  sinThetaENu*sinPhiNu*cosPhi + cosThetaENu*sinTheta*sinPhi);
124  nuDirection.setZ(-sinThetaENu*cosPhiNu*sinTheta + cosThetaENu*cosTheta);
125 
126  G4DynamicParticle* dynamicNeutrino
127  = new G4DynamicParticle(G4MT_daughters[2], nuDirection*nuEnergy);
128  products->PushProducts(dynamicNeutrino);
129 
130  // Daughter nucleus 4-vector
131  // p_D = - p_e - p_nu
132  G4DynamicParticle* dynamicDaughter =
134  -eDirection*eMomentum - nuDirection*nuEnergy);
135  products->PushProducts(dynamicDaughter);
136 
137  } else {
138  // electron energy below threshold -> no decay
139  G4DynamicParticle* noDecay =
140  new G4DynamicParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
141  products->PushProducts(noDecay);
142  }
143 
144  // Check energy conservation against Q value, not nuclear masses
145  /*
146  G4int nProd = products->entries();
147  G4DynamicParticle* temp = 0;
148  G4double Esum = 0.0;
149  for (G4int i = 0; i < nProd; i++) {
150  temp = products->operator[](i);
151  // G4cout << temp->GetParticleDefinition()->GetParticleName() << " has "
152  // << temp->GetTotalEnergy()/keV << " keV " << G4endl;
153  Esum += temp->GetKineticEnergy();
154  }
155  G4double eCons = (endpointEnergy - Esum)/keV;
156  if (std::abs(eCons) > 0.001) G4cout << " Beta- check: eCons = " << eCons << G4endl;
157  */
158  return products;
159 }
void CheckAndFillDaughters()
CLHEP::Hep3Vector G4ThreeVector
G4int PushProducts(G4DynamicParticle *aParticle)
G4ParticleDefinition * G4MT_parent
G4ParticleDefinition ** G4MT_daughters
static constexpr double rad
Definition: G4SIunits.hh:149
void setY(double)
void setZ(double)
void setX(double)
static constexpr double twopi
Definition: G4SIunits.hh:76
#define G4UniformRand()
Definition: Randomize.hh:97
G4double GetPDGMass() const
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

void G4BetaMinusDecay::DumpNuclearInfo ( )
virtual

Implements G4NuclearDecay.

Definition at line 198 of file G4BetaMinusDecay.cc.

199 {
200  G4cout << " G4BetaMinusDecay for parent nucleus " << GetParentName() << G4endl;
201  G4cout << " decays to " << GetDaughterName(0) << " , " << GetDaughterName(1)
202  << " and " << GetDaughterName(2) << " with branching ratio " << GetBR()
203  << "% and endpoint energy " << endpointEnergy/keV << " keV " << G4endl;
204 }
G4double GetBR() const
const G4String & GetParentName() const
G4GLOB_DLL std::ostream G4cout
const G4String & GetDaughterName(G4int anIndex) const
#define G4endl
Definition: G4ios.hh:61
static constexpr double keV
Definition: G4SIunits.hh:216

Here is the call graph for this function:


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