#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 G4BetaDecayType &type)

virtual	~G4BetaPlusDecay ()

virtual G4DecayProducts *	DecayIt (G4double)

virtual void	DumpNuclearInfo ()

Public Member Functions inherited from G4NuclearDecay
	G4NuclearDecay (const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation)

virtual	~G4NuclearDecay ()

G4RadioactiveDecayMode	GetDecayMode ()

G4double	GetDaughterExcitation ()

G4ParticleDefinition *	GetDaughterNucleus ()

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 G4String &	GetKinematicsName () const

G4double	GetBR () const

G4int	GetNumberOfDaughters () const

G4ParticleDefinition *	GetParent ()

G4ParticleDefinition *	GetDaughter (G4int anIndex)

G4int	GetAngularMomentum ()

const G4String &	GetParentName () const

const G4String &	GetDaughterName (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 G4ThreeVector &	GetPolarization () const

Private Member Functions
void	SetUpBetaSpectrumSampler (const G4int &parentZ, const G4int &parentA, const G4BetaDecayType &type)

Private Attributes
const G4double	endpointEnergy

G4RandGeneral *	spectrumSampler

Additional Inherited Members
Protected Member Functions inherited from G4VDecayChannel
void	ClearDaughtersName ()

void	CheckAndFillDaughters ()

void	CheckAndFillParent ()

	G4VDecayChannel ()

	G4VDecayChannel (const G4VDecayChannel &)

G4VDecayChannel &	operator= (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

G4String *	parent_name

G4String **	daughters_name

G4double	rangeMass

G4ThreeVector	parent_polarization

G4ParticleTable *	particletable

G4int	verboseLevel

G4ParticleDefinition *	G4MT_parent

G4ParticleDefinition **	G4MT_daughters

G4double	G4MT_parent_mass

G4double *	G4MT_daughters_mass

G4double *	G4MT_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::G4BetaPlusDecay	(	const G4ParticleDefinition *	theParentNucleus,
		const G4double &	theBR,
		const G4double &	endpointE,
		const G4double &	ex,
		const G4BetaDecayType &	type
	)

Definition at line 45 of file G4BetaPlusDecay.cc.

  : G4NuclearDecay("beta+ decay", BetaPlus, excitationE),
    endpointEnergy(e0 - 2.*CLHEP::electron_mass_c2)
 {
   SetParent(theParentNucleus);  // Store name of parent nucleus, delete G4MT_parent 
   SetBR(branch);
 
   SetNumberOfDaughters(3);
   G4IonTable* theIonTable =
     (G4IonTable*)(G4ParticleTable::GetParticleTable()->GetIonTable());
   G4int daughterZ = theParentNucleus->GetAtomicNumber() - 1;
   G4int daughterA = theParentNucleus->GetAtomicMass(); 
   SetDaughter(0, theIonTable->GetIon(daughterZ, daughterA, excitationE) );
   SetUpBetaSpectrumSampler(daughterZ, daughterA, betaType);
   SetDaughter(1, "e+");
   SetDaughter(2, "nu_e");
 }

Here is the call graph for this function:

◆ ~G4BetaPlusDecay()

G4BetaPlusDecay::~G4BetaPlusDecay ( )

virtual

Definition at line 67 of file G4BetaPlusDecay.cc.

 {
     delete spectrumSampler;
 }

Member Function Documentation

◆ DecayIt()

G4DecayProducts * G4BetaPlusDecay::DecayIt ( G4double )

virtual

Implements G4VDecayChannel.

Definition at line 73 of file G4BetaPlusDecay.cc.

 {
   // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)  
   CheckAndFillParent();
 
   // Fill G4MT_daughters with e-, nu and residual nucleus (stored by SetDaughter)  
   CheckAndFillDaughters();
 
   G4double parentMass = G4MT_parent->GetPDGMass();
   G4double eMass = G4MT_daughters[1]->GetPDGMass();
   G4double nucleusMass = G4MT_daughters[0]->GetPDGMass();
 
   // Set up final state
   // parentParticle is set at rest here because boost with correct momentum 
   // is done later
   G4DynamicParticle parentParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
   G4DecayProducts* products = new G4DecayProducts(parentParticle);
 
   if (spectrumSampler) {
     // Generate positron isotropic in angle, with energy from stored spectrum
     G4double eKE = endpointEnergy*spectrumSampler->shoot(G4Random::getTheEngine() );
     G4double eMomentum = std::sqrt(eKE*(eKE + 2.*eMass) );
 
     G4double cosTheta = 2.*G4UniformRand() - 1.0;
     G4double sinTheta = std::sqrt(1.0 - cosTheta*cosTheta);
     G4double phi = twopi*G4UniformRand()*rad;
     G4double sinPhi = std::sin(phi);
     G4double cosPhi = std::cos(phi);
 
     G4ParticleMomentum eDirection(sinTheta*cosPhi, sinTheta*sinPhi, cosTheta);
     G4DynamicParticle* dynamicPositron
       = new G4DynamicParticle(G4MT_daughters[1], eDirection*eMomentum);
     products->PushProducts(dynamicPositron);
 
     // Generate neutrino with angle relative to positron, and energy from
     // energy-momentum conservation using endpoint energy of reaction
     G4double cosThetaENu = 2.*G4UniformRand() - 1.;
     G4double eTE = eMass + eKE;
     G4double nuEnergy = ((endpointEnergy - eKE)*(parentMass + nucleusMass - eTE)
           - eMomentum*eMomentum)/(parentMass - eTE + eMomentum*cosThetaENu)/2.;
 
     G4double sinThetaENu = std::sqrt(1.0 - cosThetaENu*cosThetaENu);
     phi = twopi*G4UniformRand()*rad;
     G4double sinPhiNu = std::sin(phi);
     G4double cosPhiNu = std::cos(phi);
 
     G4ParticleMomentum nuDirection;
     nuDirection.setX(sinThetaENu*cosPhiNu*cosTheta*cosPhi -
                      sinThetaENu*sinPhiNu*sinPhi + cosThetaENu*sinTheta*cosPhi);
     nuDirection.setY(sinThetaENu*cosPhiNu*cosTheta*sinPhi +
                      sinThetaENu*sinPhiNu*cosPhi + cosThetaENu*sinTheta*sinPhi);
     nuDirection.setZ(-sinThetaENu*cosPhiNu*sinTheta + cosThetaENu*cosTheta);
 
     G4DynamicParticle* dynamicNeutrino
       = new G4DynamicParticle(G4MT_daughters[2], nuDirection*nuEnergy);
     products->PushProducts(dynamicNeutrino);
 
     // Generate daughter nucleus from sum of positron and neutrino 4-vectors:
     // p_D = - p_e - p_nu
     G4DynamicParticle* dynamicDaughter =
       new G4DynamicParticle(G4MT_daughters[0],
                             -eDirection*eMomentum - nuDirection*nuEnergy);
     products->PushProducts(dynamicDaughter);
 
   } else {
     // positron energy below threshold -> no decay
     G4DynamicParticle* noDecay =
       new G4DynamicParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
     products->PushProducts(noDecay);
   }
 
   // Check energy conservation against endpoint value, not nuclear masses
   /*
   G4int nProd = products->entries();
   G4DynamicParticle* temp = 0;
   G4double Esum = 0.0;
   for (G4int i = 0; i < nProd; i++) {
     temp = products->operator[](i);
     Esum += temp->GetKineticEnergy();
   }
   G4double eCons = (endpointEnergy - Esum)/keV;
   if (eCons > 0.001) G4cout << " Beta+ check: eCons (keV) = " << eCons << G4endl;
   */
   return products;
 }

Here is the call graph for this function:

◆ DumpNuclearInfo()

void G4BetaPlusDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 197 of file G4BetaPlusDecay.cc.

 {
   G4cout << " G4BetaPlusDecay for parent nucleus " << GetParentName() << G4endl;
   G4cout << " decays to " << GetDaughterName(0) << " , " << GetDaughterName(1) 
          << " and " << GetDaughterName(2) << " with branching ratio " << GetBR()
          << "% and endpoint energy " << endpointEnergy/keV << " keV " << G4endl;
 }

Here is the call graph for this function:

◆ SetUpBetaSpectrumSampler()

void G4BetaPlusDecay::SetUpBetaSpectrumSampler	(	const G4int &	parentZ,
		const G4int &	parentA,
		const G4BetaDecayType &	type
	)

private

Definition at line 161 of file G4BetaPlusDecay.cc.

 {
   G4double e0 = endpointEnergy/CLHEP::electron_mass_c2;
   G4BetaDecayCorrections corrections(daughterZ, daughterA);
   spectrumSampler = 0;
 
   // Check for cases in which Q < 2Me (e.g. z67.a162) 
   if (e0 > 0.) {
     // Array to store spectrum pdf
     G4int npti = 100;
     G4double* pdf = new G4double[npti];
 
     G4double e;  // Total positron energy in units of electron mass
     G4double p;  // Positron momentum in units of electron mass
     G4double f;  // Spectral shap function
     for (G4int ptn = 0; ptn < npti; ptn++) {
       // Calculate simple phase space
       e = 1. + e0*(ptn + 0.5)/G4double(npti);
       p = std::sqrt(e*e - 1.);
       f = p*e*(e0 - e + 1.)*(e0 - e + 1.);
 
       // Apply Fermi factor to get allowed shape
       f *= corrections.FermiFunction(e);
 
       // Apply shape factor for forbidden transitions
       f *= corrections.ShapeFactor(betaType, p, e0-e+1.);
       pdf[ptn] = f;
     }
     spectrumSampler = new G4RandGeneral(pdf, npti);
     delete[] pdf;
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

◆ endpointEnergy

const G4double G4BetaPlusDecay::endpointEnergy

private

Definition at line 61 of file G4BetaPlusDecay.hh.

◆ spectrumSampler

G4RandGeneral* G4BetaPlusDecay::spectrumSampler

private

Definition at line 62 of file G4BetaPlusDecay.hh.

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4BetaPlusDecay()

◆ ~G4BetaPlusDecay()

Member Function Documentation

◆ DecayIt()

◆ DumpNuclearInfo()

◆ SetUpBetaSpectrumSampler()

Member Data Documentation

◆ endpointEnergy

◆ spectrumSampler