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

#include <G4BetaPlusDecay.hh>

Inheritance diagram for G4BetaPlusDecay:
Collaboration diagram for G4BetaPlusDecay:

Public Member Functions

 G4BetaPlusDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &endpointE, const G4double &ex, const G4Ions::G4FloatLevelBase &flb, const G4BetaDecayType &type)
 
virtual ~G4BetaPlusDecay ()
 
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 G4BetaPlusDecay.hh.

Constructor & Destructor Documentation

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

Definition at line 45 of file G4BetaPlusDecay.cc.

50  : G4NuclearDecay("beta+ decay", BetaPlus, excitationE, flb),
51  endpointEnergy(e0 - 2.*CLHEP::electron_mass_c2)
52 {
53  SetParent(theParentNucleus); // Store name of parent nucleus, delete G4MT_parent
54  SetBR(branch);
55 
57  G4IonTable* theIonTable =
59  G4int daughterZ = theParentNucleus->GetAtomicNumber() - 1;
60  G4int daughterA = theParentNucleus->GetAtomicMass();
61  SetDaughter(0, theIonTable->GetIon(daughterZ, daughterA, excitationE, flb) );
62  SetUpBetaSpectrumSampler(daughterZ, daughterA, betaType);
63  SetDaughter(1, "e+");
64  SetDaughter(2, "nu_e");
65 }
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
static constexpr double electron_mass_c2
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:

G4BetaPlusDecay::~G4BetaPlusDecay ( )
virtual

Definition at line 68 of file G4BetaPlusDecay.cc.

69 {
70  delete spectrumSampler;
71 }

Member Function Documentation

G4DecayProducts * G4BetaPlusDecay::DecayIt ( G4double  )
virtual

Implements G4VDecayChannel.

Definition at line 74 of file G4BetaPlusDecay.cc.

75 {
76  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)
78 
79  // Fill G4MT_daughters with e-, nu and residual nucleus (stored by SetDaughter)
81 
82  G4double parentMass = G4MT_parent->GetPDGMass();
83  G4double eMass = G4MT_daughters[1]->GetPDGMass();
84  G4double nucleusMass = G4MT_daughters[0]->GetPDGMass();
85  // Set up final state
86  // parentParticle is set at rest here because boost with correct momentum
87  // is done later
88  G4DynamicParticle parentParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
89  G4DecayProducts* products = new G4DecayProducts(parentParticle);
90 
91  if (spectrumSampler) {
92  // Generate positron isotropic in angle, with energy from stored spectrum
93  G4double eKE = endpointEnergy*spectrumSampler->shoot(G4Random::getTheEngine() );
94  G4double eMomentum = std::sqrt(eKE*(eKE + 2.*eMass) );
95 
96  G4double cosTheta = 2.*G4UniformRand() - 1.0;
97  G4double sinTheta = std::sqrt(1.0 - cosTheta*cosTheta);
99  G4double sinPhi = std::sin(phi);
100  G4double cosPhi = std::cos(phi);
101 
102  G4ParticleMomentum eDirection(sinTheta*cosPhi, sinTheta*sinPhi, cosTheta);
103  G4DynamicParticle* dynamicPositron
104  = new G4DynamicParticle(G4MT_daughters[1], eDirection*eMomentum);
105  products->PushProducts(dynamicPositron);
106 
107  // Generate neutrino with angle relative to positron, and energy from
108  // energy-momentum conservation using endpoint energy of reaction
109  G4double cosThetaENu = 2.*G4UniformRand() - 1.;
110  G4double eTE = eMass + eKE;
111  G4double nuEnergy = ((endpointEnergy - eKE)*(parentMass + nucleusMass - eTE)
112  - eMomentum*eMomentum)/(parentMass - eTE + eMomentum*cosThetaENu)/2.;
113 
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  // Generate daughter nucleus from sum of positron and neutrino 4-vectors:
131  // p_D = - p_e - p_nu
132  G4DynamicParticle* dynamicDaughter =
134  -eDirection*eMomentum - nuDirection*nuEnergy);
135  products->PushProducts(dynamicDaughter);
136 
137  } else {
138  // positron 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 endpoint 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  Esum += temp->GetKineticEnergy();
152  }
153  G4double eCons = (endpointEnergy - Esum)/keV;
154  if (eCons > 0.001) G4cout << " Beta+ check: eCons (keV) = " << eCons << G4endl;
155  */
156  return products;
157 }
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 G4BetaPlusDecay::DumpNuclearInfo ( )
virtual

Implements G4NuclearDecay.

Definition at line 197 of file G4BetaPlusDecay.cc.

198 {
199  G4cout << " G4BetaPlusDecay for parent nucleus " << GetParentName() << G4endl;
200  G4cout << " decays to " << GetDaughterName(0) << " , " << GetDaughterName(1)
201  << " and " << GetDaughterName(2) << " with branching ratio " << GetBR()
202  << "% and endpoint energy " << endpointEnergy/keV << " keV " << G4endl;
203 }
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: