G4BetaPlusDecay Class Reference

#include <G4BetaPlusDecay.hh>

Inheritance 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)
	~G4BetaPlusDecay () override=default
G4DecayProducts *	DecayIt (G4double) override
void	DumpNuclearInfo () override
Public Member Functions inherited from G4NuclearDecay
	G4NuclearDecay (const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation, const G4Ions::G4FloatLevelBase &floatingLevel)
	~G4NuclearDecay () override=default
G4bool	IsOKWithParentMass (G4double parentMass) override
G4RadioactiveDecayMode	GetDecayMode () const
G4double	GetDaughterExcitation () const
G4Ions::G4FloatLevelBase	GetFloatingLevel () const
G4ParticleDefinition *	GetDaughterNucleus ()
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="", const G4String &theDaughterName5="")
virtual	~G4VDecayChannel ()
G4bool	operator== (const G4VDecayChannel &r) const
G4bool	operator!= (const G4VDecayChannel &r) const
G4bool	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)
void	SetPolarization (const G4ThreeVector &)
const G4ThreeVector &	GetPolarization () const

Additional Inherited Members
Protected Member Functions inherited from G4VDecayChannel
	G4VDecayChannel ()
	G4VDecayChannel (const G4VDecayChannel &)
G4VDecayChannel &	operator= (const G4VDecayChannel &)
void	ClearDaughtersName ()
void	CheckAndFillDaughters ()
void	CheckAndFillParent ()
G4double	DynamicalMass (G4double massPDG, G4double width, G4double maxDev=1.0) const
Protected Attributes inherited from G4NuclearDecay
const G4RadioactiveDecayMode	theMode
Protected Attributes inherited from G4VDecayChannel
G4String	kinematics_name = ""
G4double	rbranch = 0.0
G4String *	parent_name = nullptr
G4String **	daughters_name = nullptr
G4double	rangeMass = 2.5
G4ThreeVector	parent_polarization
G4ParticleTable *	particletable = nullptr
G4ParticleDefinition *	G4MT_parent = nullptr
G4ParticleDefinition **	G4MT_daughters = nullptr
G4double	G4MT_parent_mass = 0.0
G4double *	G4MT_daughters_mass = nullptr
G4double *	G4MT_daughters_width = nullptr
G4Mutex	daughtersMutex
G4Mutex	parentMutex
G4int	numberOfDaughters = 0
G4int	verboseLevel = 1
Static Protected Attributes inherited from G4VDecayChannel
static const G4String	noName = " "

Detailed Description

Definition at line 45 of file G4BetaPlusDecay.hh.

Constructor & Destructor Documentation

G4BetaPlusDecay()

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

Definition at line 56 of file G4BetaPlusDecay.cc.

 : G4NuclearDecay("beta+ decay", BetaPlus, excitationE, flb),
   maxEnergy(e0/eMass - 2.0),
   estep(maxEnergy/(G4double)(npti - 1))
{
  SetParent(theParentNucleus);  // Store name of parent nucleus, delete G4MT_parent 
  SetBR(branch);
  SetNumberOfDaughters(3);
 
  fPrimaryIon = theParentNucleus;
  fLepton = G4Positron::Positron();
  fNeutrino = G4NeutrinoE::NeutrinoE();
 
  G4IonTable* theIonTable = G4ParticleTable::GetParticleTable()->GetIonTable();
  G4int daughterZ = theParentNucleus->GetAtomicNumber() - 1;
  G4int daughterA = theParentNucleus->GetAtomicMass();
  fResIon = const_cast<const G4ParticleDefinition*>(theIonTable->GetIon(daughterZ, daughterA,
                                                                        excitationE, flb));
 
  parentMass = theParentNucleus->GetPDGMass();
  resMass = fResIon->GetPDGMass();
 
  SetUpBetaSpectrumSampler(daughterZ, daughterA, betaType);
 
  SetDaughter(0, fResIon);
  SetDaughter(1, fLepton);
  SetDaughter(2, fNeutrino);
 
  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)  
  CheckAndFillParent();
 
  // Fill G4MT_daughters with e+, nu and residual nucleus (stored by SetDaughter)  
  CheckAndFillDaughters();
}

~G4BetaPlusDecay()

G4BetaPlusDecay::~G4BetaPlusDecay ( )

overridedefault

Member Function Documentation

DecayIt()

G4DecayProducts * G4BetaPlusDecay::DecayIt ( G4double )

overridevirtual

Implements G4VDecayChannel.

Definition at line 95 of file G4BetaPlusDecay.cc.

{
  // Set up final state
  // parentParticle is set at rest here because boost with correct momentum 
  // is done later
  G4DynamicParticle prim(fPrimaryIon, G4ThreeVector(0,0,1), 0.0);
  G4DecayProducts* products = new G4DecayProducts(prim);
 
  // Generate positron isotropic in angle, with energy from stored spectrum
  const G4double eKE = eMass*G4BetaSpectrumSampler::shoot(npti, cdf, estep);
 
  G4double eMomentum = std::sqrt(eKE*(eKE + 2.*eMass));
  G4ThreeVector dir = G4RandomDirection();
  G4DynamicParticle* dp = new G4DynamicParticle(fLepton, dir, eKE);
  products->PushProducts(dp);
  /*
  G4cout << "G4BetaPlusDecay::DecayIt: " << fPrimaryIon->GetParticleName() 
     << " -> " << fResIon->GetParticleName() << " + " << fLepton->GetParticleName()
     << " + " << fNeutrino->GetParticleName() << " Ee(MeV)=" << eKE
         << G4endl;
  */
  // 4-momentum of residual ion and neutrino
  G4LorentzVector lv(-eMomentum*dir.x(), -eMomentum*dir.y(), -eMomentum*dir.z(),
                     parentMass - eKE - eMass);
 
  const G4double elim = CLHEP::eV;
  // centrum of mass system
  G4double M = lv.mag();
  G4double edel = M - resMass;
  // Free energy should be above limit
  if (edel >= elim) {
    // neutrino
    G4double eNu = 0.5*(M - resMass*resMass/M);
    G4LorentzVector lvnu(eNu*G4RandomDirection(), eNu);
    lvnu.boost(lv.boostVector());
    dir = lvnu.vect().unit();
    dp = new G4DynamicParticle(fNeutrino, dir, lvnu.e());
    products->PushProducts(dp);
 
    // residual
    lv -= lvnu;
    dir = lv.vect().unit();
    G4double ekin = std::max(lv.e() - resMass, 0.0);
    dp = new G4DynamicParticle(fResIon, dir, ekin);
    products->PushProducts(dp);
 
  } else {
    // neglecting relativistic kinematic and giving all energy to neutrino
    dp = new G4DynamicParticle(fNeutrino, G4RandomDirection(), elim);
    products->PushProducts(dp);
    dp = new G4DynamicParticle(fResIon, G4ThreeVector(0.0,0.0,1.0), 0.0);
    products->PushProducts(dp);
  }
 
  return products;
}

DumpNuclearInfo()

void G4BetaPlusDecay::DumpNuclearInfo ( )

overridevirtual

Implements G4NuclearDecay.

Definition at line 190 of file G4BetaPlusDecay.cc.

{
  G4cout << " G4BetaPlusDecay  " << fPrimaryIon->GetParticleName()
     << " -> " << fResIon->GetParticleName() << " + " << fLepton->GetParticleName()
     << " + " << fNeutrino->GetParticleName() << " Eemax(MeV)="
     << maxEnergy*eMass << " BR=" << GetBR() << "%" << G4endl;
}

---

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

• G4BetaPlusDecay.hh
• G4BetaPlusDecay.cc