#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 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 ()

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

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	(	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");
 }

Here is the call graph for this function:

G4ECDecay::~G4ECDecay ( )

virtual

Definition at line 68 of file G4ECDecay.cc.

69 {}

Member Function Documentation

G4DecayProducts * G4ECDecay::DecayIt ( G4double )

virtual

Implements G4VDecayChannel.

Definition at line 72 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:

void G4ECDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 198 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:

void G4ECDecay::SetARM ( G4bool onoff )

inline

Definition at line 56 of file G4ECDecay.hh.

56 {applyARM = onoff;}

Here is the caller graph for this function:

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

geant4.10.03.p01/source/processes/hadronic/models/radioactive_decay/include/G4ECDecay.hh
geant4.10.03.p01/source/processes/hadronic/models/radioactive_decay/src/G4ECDecay.cc

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation