#include <G4ITDecay.hh>

Inheritance diagram for G4ITDecay:

Collaboration diagram for G4ITDecay:

Public Member Functions
	G4ITDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation)

virtual	~G4ITDecay ()

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

G4int	parentZ

G4int	parentA

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 46 of file G4ITDecay.hh.

Constructor & Destructor Documentation

◆ G4ITDecay()

G4ITDecay::G4ITDecay	(	const G4ParticleDefinition *	theParentNucleus,
		const G4double &	theBR,
		const G4double &	Qvalue,
		const G4double &	excitation
	)

Definition at line 50 of file G4ITDecay.cc.

  : G4NuclearDecay("IT decay", IT, excitationE), transitionQ(Qvalue), 
    applyARM(true)
 {
   SetParent(theParentNucleus);  // Store name of parent nucleus, delete G4MT_parent 
   SetBR(branch);
 
   parentZ = theParentNucleus->GetAtomicNumber();
   parentA = theParentNucleus->GetAtomicMass(); 
 
   SetNumberOfDaughters(1);
   G4IonTable* theIonTable =
     (G4IonTable*)(G4ParticleTable::GetParticleTable()->GetIonTable());
   SetDaughter(0, theIonTable->GetIon(parentZ, parentA, excitationE) );
 }

Here is the call graph for this function:

◆ ~G4ITDecay()

G4ITDecay::~G4ITDecay ( )

virtual

Definition at line 69 of file G4ITDecay.cc.

70 {}

Member Function Documentation

◆ DecayIt()

G4DecayProducts * G4ITDecay::DecayIt ( G4double )

virtual

Implements G4VDecayChannel.

Definition at line 73 of file G4ITDecay.cc.

 {
   // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)  
   CheckAndFillParent();
 
   // Set up final state
   // parentParticle is set at rest here because boost with correct momentum 
   // is done later
   G4LorentzVector atRest(G4MT_parent->GetPDGMass(),
                          G4ThreeVector(0.,0.,0.) );
   G4DynamicParticle parentParticle(G4MT_parent, atRest);
   G4DecayProducts* products = new G4DecayProducts(parentParticle);
 
   // Let G4PhotonEvaporation do the decay 
   G4PhotonEvaporation* photoEvap = new G4PhotonEvaporation;
   photoEvap->RDMForced(true);
   photoEvap->SetICM(true);
 
   G4Fragment* nucleus = new G4Fragment(parentA, parentZ, atRest);
   G4Fragment* eOrGamma = photoEvap->EmittedFragment(nucleus);
 
   // Daughter nuclide is returned in nucleus pointer
   G4double finalDaughterExcitation = nucleus->GetExcitationEnergy();
   if (finalDaughterExcitation < 1*keV) finalDaughterExcitation = 0.0;
   G4LorentzVector daughterMomentum = nucleus->GetMomentum();
   G4IonTable* theIonTable =
     (G4IonTable*)(G4ParticleTable::GetParticleTable()->GetIonTable() );
   G4ParticleDefinition* daughterIon =
     theIonTable->GetIon(parentZ, parentA, finalDaughterExcitation);
   G4DynamicParticle* dynDaughter = new G4DynamicParticle(daughterIon,
                                                          daughterMomentum); 
   delete nucleus;
 
   if (eOrGamma) {
     G4DynamicParticle* eOrGammaDyn =
       new G4DynamicParticle(eOrGamma->GetParticleDefinition(),
                             eOrGamma->GetMomentum() );
     eOrGammaDyn->SetProperTime(eOrGamma->GetCreationTime() );
     products->PushProducts(eOrGammaDyn);
     delete eOrGamma;
 
     // Now do atomic relaxation
     if (applyARM) {
       G4int shellIndex = photoEvap->GetVacantShellNumber();
       if (shellIndex > -1) {
         G4VAtomDeexcitation* atomDeex =
           G4LossTableManager::Instance()->AtomDeexcitation();
         if (atomDeex->IsFluoActive() && parentZ > 5 && parentZ < 100) {
           G4int nShells = G4AtomicShells::GetNumberOfShells(parentZ);
           if (shellIndex >= nShells) shellIndex = nShells;
           G4AtomicShellEnumerator as = G4AtomicShellEnumerator(shellIndex);
           const G4AtomicShell* shell = atomDeex->GetAtomicShell(parentZ, as);
           std::vector<G4DynamicParticle*> armProducts;
 
           // 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, parentZ, 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);
           } 
 
           G4int nArm = armProducts.size();
           if (nArm > 0) {
             G4ThreeVector bst = dynDaughter->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);
             }
           }
         }
       }
     } // if ARM on 
   } // eOrGamma
 
   products->PushProducts(dynDaughter);
 
   // 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 = G4MT_parent->GetPDGMass() - dynDaughter->GetMass() - KEsum;
   G4cout << " IT check: Ediff (keV) = " << eCons/keV << G4endl; 
   */
   delete photoEvap;
 
   return products;
 }

Here is the call graph for this function:

◆ DumpNuclearInfo()

void G4ITDecay::DumpNuclearInfo ( )

virtual

Implements G4NuclearDecay.

Definition at line 190 of file G4ITDecay.cc.

 {
   G4cout << " G4ITDecay for parent nucleus " << GetParentName() << G4endl;
   G4cout << " decays to " << GetDaughterName(0)
          << " + gammas (or electrons), with branching ratio " << GetBR()
          << "% and Q value " << transitionQ << G4endl;
 }

Here is the call graph for this function:

◆ SetARM()

void G4ITDecay::SetARM ( G4bool onoff )

inline

Definition at line 59 of file G4ITDecay.hh.

59 {applyARM = onoff;}

G4ITDecay::applyARM

G4bool applyARM

Definition: G4ITDecay.hh:65

Here is the caller graph for this function:

Member Data Documentation

◆ applyARM

G4bool G4ITDecay::applyARM

private

Definition at line 65 of file G4ITDecay.hh.

◆ parentA

G4int G4ITDecay::parentA

private

Definition at line 64 of file G4ITDecay.hh.

◆ parentZ

G4int G4ITDecay::parentZ

private

Definition at line 63 of file G4ITDecay.hh.

◆ transitionQ

const G4double G4ITDecay::transitionQ

private

Definition at line 62 of file G4ITDecay.hh.

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4ITDecay()

◆ ~G4ITDecay()

Member Function Documentation

◆ DecayIt()

◆ DumpNuclearInfo()

◆ SetARM()

Member Data Documentation

◆ applyARM

◆ parentA

◆ parentZ

◆ transitionQ