G4ECDecay Class Reference

#include <G4ECDecay.hh>

Inheritance diagram for G4ECDecay:

Public Member Functions
	G4ECDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation, const G4Ions::G4FloatLevelBase &flb, const G4RadioactiveDecayMode &mode)
virtual	~G4ECDecay ()
virtual G4DecayProducts *	DecayIt (G4double)
virtual void	DumpNuclearInfo ()
void	SetARM (G4bool onoff)
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 ()
G4ParticleDefinition *	GetDaughterNucleus ()
void	SetHLThreshold (G4double HLT)
G4double	GetHLThreshold ()
virtual void	DumpNuclearInfo ()=0
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
virtual G4DecayProducts *	DecayIt (G4double parentMass=-1.0)=0
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

Additional Inherited Members
Protected Member Functions inherited from G4VDecayChannel
void	ClearDaughtersName ()
void	CheckAndFillDaughters ()
void	CheckAndFillParent ()
G4double	DynamicalMass (G4double massPDG, G4double width, G4double maxDev=1.0) const
	G4VDecayChannel ()
	G4VDecayChannel (const G4VDecayChannel &)
G4VDecayChannel &	operator= (const G4VDecayChannel &)
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 42 of file G4ECDecay.hh.

Constructor & Destructor Documentation

G4ECDecay()

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

Definition at line 47 of file G4ECDecay.cc.

 : G4NuclearDecay("electron capture", mode, excitationE, flb), transitionQ(Qvalue),
   applyARM(true)
{
  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(); 
  SetDaughter(0, theIonTable->GetIon(daughterZ, daughterA, excitationE, flb) );
  SetDaughter(1, "nu_e");
  DefineSubshellProbabilities(daughterZ, daughterZ);
}

~G4ECDecay()

G4ECDecay::~G4ECDecay ( )

virtual

Definition at line 69 of file G4ECDecay.cc.

{}

Member Function Documentation

DecayIt()

G4DecayProducts * G4ECDecay::DecayIt ( G4double  )

virtual

Implements G4VDecayChannel.

Definition at line 73 of file G4ECDecay.cc.

{
  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)  
  CheckAndFillParent();
 
  // Fill G4MT_daughters with alpha and residual nucleus (stored by SetDaughter)  
  CheckAndFillDaughters();
 
  // Get shell number of captured electron
  G4int shellIndex = -1;
  G4double ran;
  switch (theMode)
    {
    case KshellEC:
      shellIndex = 0;
      break;
    case LshellEC: // PL1+PL2+PL3=1
      ran=G4UniformRand();
      if (ran <= PL1) shellIndex =1;
      else if (ran<= (PL1+PL2)) shellIndex =2;
      else shellIndex =3;
      break;
    case MshellEC:  // PM1+PM2+PM3=1
      ran=G4UniformRand();
      if (ran < PM1) shellIndex =4;
      else if (ran< (PM1+PM2)) shellIndex =5;
      else shellIndex = 6;
      break;
    case NshellEC: // PN1+PN2+PN3=1
      ran=G4UniformRand();
      if (ran < PN1) shellIndex = 9;
      else if (ran<= (PN1+PN2)) shellIndex =2;
      else shellIndex =10;
      break;
    default:
      G4Exception("G4ECDecay::DecayIt()", "HAD_RDM_009",
                  FatalException, "Invalid electron shell selected");
    }
 
  // Initialize decay products with parent nucleus at rest
  G4DynamicParticle parentParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
  G4DecayProducts* products = new G4DecayProducts(parentParticle);
  G4double eBind = 0.0;
 
  // G4LossTableManager must already be initialized with G4UAtomicDeexcitation
  // This is currently done in G4RadioactiveDecay::BuildPhysicsTable
  G4VAtomDeexcitation* atomDeex = 
          G4LossTableManager::Instance()->AtomDeexcitation();
  std::vector<G4DynamicParticle*> armProducts;
 
  if (applyARM) {
    if (nullptr != atomDeex) {
      G4int aZ = G4MT_daughters[0]->GetAtomicNumber();
      G4int nShells = G4AtomicShells::GetNumberOfShells(aZ);
      if (shellIndex >= nShells) shellIndex = nShells;
      G4AtomicShellEnumerator as = G4AtomicShellEnumerator(shellIndex);
      const G4AtomicShell* shell = atomDeex->GetAtomicShell(aZ, as);
      eBind = shell->BindingEnergy(); 
      if (atomDeex->IsFluoActive() && aZ > 5 && aZ < 105) {
        // Do atomic relaxation
    // VI, SI
    // Allows fixing of Bugzilla 1727
    //const G4double deexLimit = 0.1*keV;
    G4double deexLimit = 0.1*keV;
    if (G4EmParameters::Instance()->DeexcitationIgnoreCut()) deexLimit = 0.;
 
        atomDeex->GenerateParticles(&armProducts, shell, aZ, deexLimit, deexLimit);
      }
 
      G4double productEnergy = 0.;
      for (std::size_t i = 0; i < armProducts.size(); ++i) {
        productEnergy += armProducts[i]->GetKineticEnergy();
      }
      G4double deficit = shell->BindingEnergy() - productEnergy;
      if (deficit > 0.0) {
        // Add a dummy electron to make up extra energy
        G4double cosTh = 1.-2.*G4UniformRand();
        G4double sinTh = std::sqrt(1.- cosTh*cosTh);
        G4double phi = twopi*G4UniformRand();
 
        G4ThreeVector electronDirection(sinTh*std::sin(phi),
                                        sinTh*std::cos(phi), cosTh);
        G4DynamicParticle* extra =
          new G4DynamicParticle(G4Electron::Electron(), electronDirection,
                                deficit);
        armProducts.push_back(extra);
      }
    } // atomDeex
  }  // applyARM
 
  G4double daughterMass = G4MT_daughters[0]->GetPDGMass();
 
  // CM momentum using Q value corrected for binding energy of captured electron
  G4double Q = transitionQ - eBind;
 
  // Negative transitionQ values for some rare nuclides require this
  // Absolute values in these cases are small 
  if (Q < 0.0) Q = 0.0;
 
  G4double cmMomentum = Q*(Q + 2.*daughterMass)/(Q + daughterMass)/2.;
 
  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 = cmMomentum;
  G4DynamicParticle* daughterParticle =
    new G4DynamicParticle(G4MT_daughters[1], direction, KE, 0.0);
  products->PushProducts(daughterParticle);
 
  KE = std::sqrt(cmMomentum*cmMomentum + daughterMass*daughterMass) - daughterMass;
  daughterParticle =
    new G4DynamicParticle(G4MT_daughters[0], -1.0*direction, KE, daughterMass);
  products->PushProducts(daughterParticle);
 
  std::size_t nArm = armProducts.size();
  if (nArm > 0) {
    G4ThreeVector bst = daughterParticle->Get4Momentum().boostVector();
    for (std::size_t i = 0; i < nArm; ++i) {
      G4DynamicParticle* dp = armProducts[i];
      G4LorentzVector lv = dp->Get4Momentum().boost(bst);
      dp->Set4Momentum(lv);
      products->PushProducts(dp);
    }
  }
 
  // Energy conservation check
  /*
  G4int newSize = products->entries();
  G4DynamicParticle* temp = 0;
  G4double KEsum = 0.0;
  for (G4int i = 0; i < newSize; i++) {
    temp = products->operator[](i);
    KEsum += temp->GetKineticEnergy();
  }
 
  G4double eCons = (transitionQ - KEsum)/keV; 
  G4cout << " EC check: Ediff (keV) = " << eCons << G4endl; 
  */
  return products;
}

DumpNuclearInfo()

void G4ECDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 217 of file G4ECDecay.cc.

{
  G4cout << " G4ECDecay of parent nucleus " << GetParentName() << " from ";
  if (theMode == 3) {
    G4cout << "K shell";
  } else if (theMode == 4) {
    G4cout << "L shell";
  } else if (theMode == 5) {
    G4cout << "M shell";
  }
  else if (theMode == 6) {
    G4cout << "N shell";
  }
  G4cout << G4endl;
  G4cout << " to " << GetDaughterName(0) << " + " << GetDaughterName(1)
         << " with branching ratio " << GetBR() << "% and Q value "
         << transitionQ << G4endl;
}

SetARM()

void G4ECDecay::SetARM ( G4bool  onoff )

inline

Definition at line 56 of file G4ECDecay.hh.

{applyARM = onoff;}

Referenced by G4RadioactiveDecay::LoadDecayTable().

---

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

• G4ECDecay.hh
• G4ECDecay.cc