G4ECDecay Class Reference

#include <G4ECDecay.hh>

Inheritance diagram for G4ECDecay:

Collaboration diagram for G4ECDecay:

Public Member Functions
	G4ECDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation, 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)
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 Attributes
const G4double	transitionQ
G4bool	applyARM

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 42 of file G4ECDecay.hh.

Constructor & Destructor Documentation

G4ECDecay()

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

Definition at line 47 of file G4ECDecay.cc.

 : G4NuclearDecay("electron capture", mode, excitationE), 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) );
  SetDaughter(1, "nu_e");
}

Here is the call graph for this function:

~G4ECDecay()

G4ECDecay::~G4ECDecay ( )

virtual

Definition at line 67 of file G4ECDecay.cc.

{}

Member Function Documentation

DecayIt()

G4DecayProducts * G4ECDecay::DecayIt ( G4double  )

virtual

Implements G4VDecayChannel.

Definition at line 71 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;
  switch (theMode)
    {
    case KshellEC:
      shellIndex = 0;
      break;
    case LshellEC:
      shellIndex = G4int(G4UniformRand()*3) + 1;
      break;
    case MshellEC:
      shellIndex = G4int(G4UniformRand()*3) + 4;
      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 (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 < 100) {
        // 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 (G4int i = 0; i < G4int(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; 
  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);

  G4int nArm = armProducts.size();
  if (nArm > 0) {
    G4ThreeVector bst = daughterParticle->Get4Momentum().boostVector();
    for (G4int 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;
}

Here is the call graph for this function:

DumpNuclearInfo()

void G4ECDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 197 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";
  }
  G4cout << G4endl;
  G4cout << " to " << GetDaughterName(0) << " + " << GetDaughterName(1)
         << " with branching ratio " << GetBR() << "% and Q value "
         << transitionQ << G4endl;
}

Here is the call graph for this function:

SetARM()

void G4ECDecay::SetARM ( G4bool  onoff )

inline

Definition at line 55 of file G4ECDecay.hh.

{applyARM = onoff;}

Here is the caller graph for this function:

Member Data Documentation

applyARM

G4bool G4ECDecay::applyARM

private

Definition at line 59 of file G4ECDecay.hh.

transitionQ

const G4double G4ECDecay::transitionQ

private

Definition at line 58 of file G4ECDecay.hh.

---

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

• G4ECDecay.hh
• G4ECDecay.cc