#include <G4ProtonDecay.hh>

Inheritance diagram for G4ProtonDecay:

Public Member Functions
	G4ProtonDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation, const G4Ions::G4FloatLevelBase &flb)

virtual	~G4ProtonDecay ()

virtual G4DecayProducts *	DecayIt (G4double)

virtual void	DumpNuclearInfo ()

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 44 of file G4ProtonDecay.hh.

Constructor & Destructor Documentation

◆ G4ProtonDecay()

G4ProtonDecay::G4ProtonDecay	(	const G4ParticleDefinition *	theParentNucleus,
		const G4double &	theBR,
		const G4double &	Qvalue,
		const G4double &	excitation,
		const G4Ions::G4FloatLevelBase &	flb )

Definition at line 45 of file G4ProtonDecay.cc.

 : G4NuclearDecay("proton decay", Proton, excitationE, flb), transitionQ(Qvalue)
{
  SetParent(theParentNucleus);  // Store name of parent nucleus, delete G4MT_parent
  SetBR(branch);
 
  SetNumberOfDaughters(2);
  G4IonTable* theIonTable =
    (G4IonTable*)(G4ParticleTable::GetParticleTable()->GetIonTable());
  G4int daughterZ = theParentNucleus->GetAtomicNumber() - 1;
  G4int daughterA = theParentNucleus->GetAtomicMass() - 1;
  SetDaughter(0, theIonTable->GetIon(daughterZ, daughterA, excitationE, flb) );
  SetDaughter(1, "proton");
}

◆ ~G4ProtonDecay()

G4ProtonDecay::~G4ProtonDecay ( )

virtual

Definition at line 64 of file G4ProtonDecay.cc.

65{}

Member Function Documentation

◆ DecayIt()

G4DecayProducts * G4ProtonDecay::DecayIt ( G4double )

virtual

Implements G4VDecayChannel.

Definition at line 68 of file G4ProtonDecay.cc.

{
  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)  
  CheckAndFillParent();
 
  // Fill G4MT_daughters with proton and residual nucleus (stored by SetDaughter)  
  CheckAndFillDaughters();
 
  G4double protonMass = G4MT_daughters[1]->GetPDGMass();
  // Excitation energy included in PDG mass
  G4double nucleusMass = G4MT_daughters[0]->GetPDGMass();
 
  // Q value was calculated from atomic masses.
  // Use it to get correct proton energy.
  G4double cmMomentum = std::sqrt(transitionQ*(transitionQ + 2.*protonMass)*
                                (transitionQ + 2.*nucleusMass)*
                                (transitionQ + 2.*protonMass + 2.*nucleusMass) )/
                                (transitionQ + protonMass + nucleusMass)/2.; 
 
  // 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);
 
  G4double costheta = 2.*G4UniformRand()-1.0;
  G4double sintheta = std::sqrt(1.0 - costheta*costheta);
  G4double phi  = twopi*G4UniformRand()*rad;
  G4ThreeVector direction(sintheta*std::cos(phi),sintheta*std::sin(phi),
                          costheta);
 
  G4double KE = std::sqrt(cmMomentum*cmMomentum + protonMass*protonMass)
              - protonMass;
  G4DynamicParticle* daughterparticle =
    new G4DynamicParticle(G4MT_daughters[1], direction, KE, protonMass);
  products->PushProducts(daughterparticle);
 
  KE = std::sqrt(cmMomentum*cmMomentum + nucleusMass*nucleusMass) - nucleusMass;
  daughterparticle =
    new G4DynamicParticle(G4MT_daughters[0], -1.0*direction, KE, nucleusMass);
  products->PushProducts(daughterparticle);
 
  // Energy conservation check
  // For proton decays, do final energy check against reaction Q value
  // which is well-measured using atomic mass differences.  Nuclear masses
  // should not be used since they are not usually directly measured and we
  // always decay atoms and not fully stripped nuclei.
  /*
  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 = (transitionQ - Esum)/keV;
  if (eCons > 1.e-07) G4cout << " Proton decay check: Ediff (keV) = " << eCons << G4endl;
  */
  return products;
}

◆ DumpNuclearInfo()

void G4ProtonDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 130 of file G4ProtonDecay.cc.

{
  G4cout << " G4ProtonDecay for parent nucleus " << GetParentName() << G4endl;
  G4cout << " decays to " << GetDaughterName(0) << " + " << GetDaughterName(1)
         << " with branching ratio " << GetBR() << "% and Q value "
         << transitionQ << G4endl;
}

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4ProtonDecay()

◆ ~G4ProtonDecay()

Member Function Documentation

◆ DecayIt()

◆ DumpNuclearInfo()