G4ITDecay Class Reference

#include <G4ITDecay.hh>

Inheritance diagram for G4ITDecay:

Public Member Functions
	G4ITDecay (G4PhotonEvaporation *aPhotonEvap)
	G4ITDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation)
	~G4ITDecay () override=default
void	SetupDecay (const G4ParticleDefinition *)
G4DecayProducts *	DecayIt (G4double) override
void	DumpNuclearInfo () override
void	SetARM (G4bool onoff)
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 45 of file G4ITDecay.hh.

Constructor & Destructor Documentation

G4ITDecay() [1/2]

G4ITDecay::G4ITDecay

(

G4PhotonEvaporation *

aPhotonEvap

)

Definition at line 50 of file G4ITDecay.cc.

  : G4NuclearDecay("IT decay", IT, 0.0, noFloat), photonEvaporation(ptr)
{}

G4ITDecay() [2/2]

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

Definition at line 54 of file G4ITDecay.cc.

  : G4NuclearDecay("IT decay", IT, excitationE, noFloat)
{
  SetParent(theParentNucleus);  // Store name of parent nucleus, delete G4MT_parent 
  SetBR(branch);
  SetNumberOfDaughters(1);
  SetDaughter(0, theParentNucleus);
 
  SetupDecay(theParentNucleus);
}

~G4ITDecay()

G4ITDecay::~G4ITDecay ( )

overridedefault

Member Function Documentation

DecayIt()

G4DecayProducts * G4ITDecay::DecayIt ( G4double )

overridevirtual

Implements G4VDecayChannel.

Definition at line 74 of file G4ITDecay.cc.

{
  // Set up final state
  // parentParticle is set at rest here because boost with correct momentum 
  // is done later
  G4LorentzVector atRest(theParent->GetPDGMass(), G4ThreeVector(0.,0.,0.) );
  G4DynamicParticle parentParticle(theParent, atRest);
  G4DecayProducts* products = new G4DecayProducts(parentParticle);
 
  // Let G4PhotonEvaporation do the decay
  G4Fragment parentNucleus(parentA, parentZ, atRest);
 
  // one emission, parent nucleaus become less excited
  G4Fragment* eOrGamma = photonEvaporation->EmittedFragment(&parentNucleus);
 
  // Modified nuclide is returned as dynDaughter
  auto theIonTable = G4ParticleTable::GetParticleTable()->GetIonTable();
  G4ParticleDefinition* daughterIon =
    theIonTable->GetIon(parentZ, parentA, parentNucleus.GetExcitationEnergy(), 
                        G4Ions::FloatLevelBase(parentNucleus.GetFloatingLevelNumber()));
  G4DynamicParticle* dynDaughter = new G4DynamicParticle(daughterIon,
                                                         parentNucleus.GetMomentum());
 
  if (nullptr != eOrGamma) {
    G4DynamicParticle* eOrGammaDyn =
      new G4DynamicParticle(eOrGamma->GetParticleDefinition(),
                            eOrGamma->GetMomentum() );
    eOrGammaDyn->SetProperTime(eOrGamma->GetCreationTime() );
    products->PushProducts(eOrGammaDyn);
    delete eOrGamma;
 
    // Now do atomic relaxation if e- is emitted
    if (applyARM) {
      G4int shellIndex = photonEvaporation->GetVacantShellNumber();
      if (shellIndex > -1) {
        G4VAtomDeexcitation* atomDeex =
          G4LossTableManager::Instance()->AtomDeexcitation();
        if (atomDeex->IsFluoActive() && parentZ > 5 && parentZ < 105) {
          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
          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);
          } 
 
          std::size_t nArm = armProducts.size();
          if (nArm > 0) {
            G4ThreeVector bst = dynDaughter->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);
            }
          }
        }
      }
    } // 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; 
  */
  return products;
}

Referenced by G4Radioactivation::AddDeexcitationSpectrumForBiasMode(), and G4RadioactiveDecay::DoDecay().

DumpNuclearInfo()

void G4ITDecay::DumpNuclearInfo ( )

overridevirtual

Implements G4NuclearDecay.

Definition at line 178 of file G4ITDecay.cc.

{
  if (theParent != nullptr) {
    G4cout << " G4ITDecay for parent nucleus " << theParent->GetParticleName() << G4endl;
  }
}

SetARM()

void G4ITDecay::SetARM ( G4bool onoff )

inline

Definition at line 65 of file G4ITDecay.hh.

{applyARM = onoff;}

Referenced by G4RadioactiveDecay::BuildPhysicsTable().

SetupDecay()

void G4ITDecay::SetupDecay

(

const G4ParticleDefinition *

theParentNucleus

)

Definition at line 67 of file G4ITDecay.cc.

{
  theParent = theParentNucleus;
  parentZ = theParentNucleus->GetAtomicNumber();
  parentA = theParentNucleus->GetAtomicMass();
}

Referenced by G4Radioactivation::AddDeexcitationSpectrumForBiasMode(), G4RadioactiveDecay::DoDecay(), and G4ITDecay().

---

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

• G4ITDecay.hh
• G4ITDecay.cc