G4DecayTable Class Reference

#include <G4DecayTable.hh>

Public Types
using	G4VDecayChannelVector = std::vector< G4VDecayChannel * >

Public Member Functions
	G4DecayTable ()
	~G4DecayTable ()
	G4DecayTable (const G4DecayTable &)=delete
G4DecayTable &	operator= (const G4DecayTable &)=delete
G4bool	operator== (const G4DecayTable &right) const
G4bool	operator!= (const G4DecayTable &right) const
void	Insert (G4VDecayChannel *aChannel)
G4int	entries () const
G4VDecayChannel *	SelectADecayChannel (G4double parentMass=-1.)
G4VDecayChannel *	GetDecayChannel (G4int index) const
G4VDecayChannel *	operator[] (G4int index)
void	DumpInfo () const

Detailed Description

Definition at line 45 of file G4DecayTable.hh.

Member Typedef Documentation

G4VDecayChannelVector

using G4DecayTable::G4VDecayChannelVector = std::vector<G4VDecayChannel*>

Definition at line 49 of file G4DecayTable.hh.

Constructor & Destructor Documentation

G4DecayTable() [1/2]

G4DecayTable::G4DecayTable

(

)

Definition at line 35 of file G4DecayTable.cc.

{
  channels =  new G4VDecayChannelVector;
}

~G4DecayTable()

G4DecayTable::~G4DecayTable

(

)

Definition at line 40 of file G4DecayTable.cc.

{
  // remove and delete all contents  
  for (auto iCh = channels->cbegin(); iCh!= channels->cend(); ++iCh)
  {
    delete (*iCh);
  }
  channels->clear();
  delete  channels;
  channels = nullptr;
  parent = nullptr;
}    

G4DecayTable() [2/2]

G4DecayTable::G4DecayTable ( const G4DecayTable &  )

delete

Member Function Documentation

DumpInfo()

void G4DecayTable::DumpInfo

(

)

const

Definition at line 120 of file G4DecayTable.cc.

{
  G4cout << "G4DecayTable:  " << parent->GetParticleName() << G4endl;
  G4int index =0;
  for (auto iCh = channels->cbegin(); iCh!= channels->cend(); ++iCh)
  {
    G4cout << index << ": ";
    (*iCh)->DumpInfo();
    index += 1;
  }
  G4cout << G4endl;
}

Referenced by G4Radioactivation::DecayIt(), G4IsotopeProperty::DumpInfo(), G4ParticleDefinition::DumpTable(), G4RadioactiveDecay::LoadDecayTable(), and G4DecayTableMessenger::SetNewValue().

entries()

G4int G4DecayTable::entries ( ) const

inline

Definition at line 100 of file G4DecayTable.hh.

{
  return G4int(channels->size());
}

Referenced by G4DecayWithSpin::AtRestDoIt(), G4Radioactivation::CalculateChainsFromParent(), G4IonTable::CreateIon(), G4Decay::DecayIt(), G4Radioactivation::DecayIt(), G4RadioactiveDecay::DecayIt(), G4KineticTrack::G4KineticTrack(), G4TextPPReporter::GeneratePropertyTable(), G4SampleResonance::GetMinimumMass(), G4RadioactiveDecay::LoadDecayTable(), G4TextPPRetriever::ModifyPropertyTable(), G4DecayWithSpin::PostStepDoIt(), G4MTRunManagerKernel::SetUpDecayChannels(), and G4TaskRunManagerKernel::SetUpDecayChannels().

GetDecayChannel()

G4VDecayChannel * G4DecayTable::GetDecayChannel ( G4int  index ) const

inline

Definition at line 112 of file G4DecayTable.hh.

{
  G4VDecayChannel* selectedChannel = nullptr;
  if ( (index>=0) && (index<G4int(channels->size())) )
  {
    selectedChannel = (*channels)[index];
  }
  return selectedChannel;
}

Referenced by G4DecayWithSpin::AtRestDoIt(), G4Radioactivation::CalculateChainsFromParent(), G4IonTable::CreateIon(), G4KineticTrack::Decay(), G4Decay::DecayIt(), G4Radioactivation::DecayIt(), G4KineticTrack::G4KineticTrack(), G4TextPPReporter::GeneratePropertyTable(), G4SampleResonance::GetMinimumMass(), G4RadioactiveDecay::LoadDecayTable(), G4TextPPRetriever::ModifyPropertyTable(), G4DecayWithSpin::PostStepDoIt(), G4DecayTableMessenger::SetNewValue(), G4MTRunManagerKernel::SetUpDecayChannels(), and G4TaskRunManagerKernel::SetUpDecayChannels().

Insert()

void G4DecayTable::Insert

(

G4VDecayChannel *

aChannel

)

Definition at line 53 of file G4DecayTable.cc.

{
  if (parent == nullptr)
  { 
    parent = (G4ParticleDefinition*)(aChannel->GetParent()); 
  }
  if (parent != aChannel->GetParent())
  {
#ifdef G4VERBOSE
    G4cout << " G4DecayTable::Insert :: bad G4VDecayChannel (mismatch parent) "
           << "       " << parent->GetParticleName()
           << " input:" << aChannel->GetParent()->GetParticleName() << G4endl;
#endif
  }
  else
  {
    G4double br = aChannel->GetBR();
    for (auto iCh = channels->cbegin(); iCh!= channels->cend(); ++iCh)
    {
      if (br > (*iCh)->GetBR())
      {
        channels->insert(iCh,aChannel);
        return;
      }
    }
    channels->push_back(aChannel);
  }
}

Referenced by G4ExcitedMesonConstructor::Add2EtaMode(), G4ExcitedMesonConstructor::Add2KMode(), G4ExcitedMesonConstructor::Add2KPiMode(), G4ExcitedMesonConstructor::Add2PiEtaMode(), G4ExcitedMesonConstructor::Add2PiMode(), G4ExcitedMesonConstructor::Add2PiOmegaMode(), G4ExcitedMesonConstructor::Add2PiRhoMode(), G4ExcitedMesonConstructor::Add3PiMode(), G4ExcitedMesonConstructor::Add4PiMode(), G4ExcitedMesonConstructor::AddKEtaMode(), G4ExcitedMesonConstructor::AddKKStarMode(), G4ExcitedMesonConstructor::AddKOmegaMode(), G4ExcitedMesonConstructor::AddKPiMode(), G4ExcitedMesonConstructor::AddKRhoMode(), G4ExcitedMesonConstructor::AddKStar2PiMode(), G4ExcitedMesonConstructor::AddKStarPiMode(), G4ExcitedMesonConstructor::AddKTwoPiMode(), G4ExcitedMesonConstructor::AddPiA2Mode(), G4ExcitedMesonConstructor::AddPiEtaMode(), G4ExcitedMesonConstructor::AddPiF0Mode(), G4ExcitedMesonConstructor::AddPiF2Mode(), G4ExcitedMesonConstructor::AddPiGammaMode(), G4ExcitedMesonConstructor::AddPiOmegaMode(), G4ExcitedMesonConstructor::AddPiRhoMode(), G4ExcitedMesonConstructor::AddRhoEtaMode(), G4ExcitedMesonConstructor::AddRhoGammaMode(), G4HadronicBuilder::BuildDecayTableForBCHadrons(), G4Radioactivation::CalculateChainsFromParent(), G4ShortLivedConstructor::ConstructBaryons(), G4ShortLivedConstructor::ConstructMesons(), G4MuonicAtomHelper::ConstructMuonicAtom(), G4AntiLambda::Definition(), G4AntiNeutron::Definition(), G4AntiOmegaMinus::Definition(), G4AntiSigmabMinus::Definition(), G4AntiSigmabPlus::Definition(), G4AntiSigmabZero::Definition(), G4AntiSigmacPlus::Definition(), G4AntiSigmacPlusPlus::Definition(), G4AntiSigmacZero::Definition(), G4AntiSigmaMinus::Definition(), G4AntiSigmaPlus::Definition(), G4AntiSigmaZero::Definition(), G4AntiXiMinus::Definition(), G4AntiXiZero::Definition(), G4Lambda::Definition(), G4Neutron::Definition(), G4OmegaMinus::Definition(), G4SigmabMinus::Definition(), G4SigmabPlus::Definition(), G4SigmabZero::Definition(), G4SigmacPlus::Definition(), G4SigmacPlusPlus::Definition(), G4SigmacZero::Definition(), G4SigmaMinus::Definition(), G4SigmaPlus::Definition(), G4SigmaZero::Definition(), G4XiMinus::Definition(), G4XiZero::Definition(), G4AntiDoubleHyperDoubleNeutron::Definition(), G4AntiDoubleHyperH4::Definition(), G4AntiHyperAlpha::Definition(), G4AntiHyperH4::Definition(), G4AntiHyperHe5::Definition(), G4AntiHyperTriton::Definition(), G4DoubleHyperDoubleNeutron::Definition(), G4DoubleHyperH4::Definition(), G4HyperAlpha::Definition(), G4HyperH4::Definition(), G4HyperHe5::Definition(), G4HyperTriton::Definition(), G4AntiKaonZero::Definition(), G4Eta::Definition(), G4EtaPrime::Definition(), G4KaonMinus::Definition(), G4KaonPlus::Definition(), G4KaonZero::Definition(), G4KaonZeroLong::Definition(), G4KaonZeroShort::Definition(), G4PionMinus::Definition(), G4PionPlus::Definition(), G4PionZero::Definition(), G4MuonMinus::Definition(), G4MuonPlus::Definition(), G4TauMinus::Definition(), G4TauPlus::Definition(), and G4RadioactiveDecay::LoadDecayTable().

operator!=()

G4bool G4DecayTable::operator!= ( const G4DecayTable &  right ) const

inline

Definition at line 94 of file G4DecayTable.hh.

{
  return (this != &right);
}

operator=()

G4DecayTable & G4DecayTable::operator= ( const G4DecayTable &  )

delete

operator==()

G4bool G4DecayTable::operator== ( const G4DecayTable &  right ) const

inline

Definition at line 88 of file G4DecayTable.hh.

{
  return (this == &right);
}

operator[]()

G4VDecayChannel * G4DecayTable::operator[] ( G4int  index )

inline

Definition at line 106 of file G4DecayTable.hh.

{
  return (*channels)[index];
}

SelectADecayChannel()

G4VDecayChannel * G4DecayTable::SelectADecayChannel

(

G4double

parentMass = -1.

)

Definition at line 82 of file G4DecayTable.cc.

{
  // check if contents exist
  if (channels->size()<1) return nullptr;
 
  if(parentMass < 0.) parentMass=parent->GetPDGMass(); 
 
  G4double sumBR = 0.;
  for (auto iCh = channels->cbegin(); iCh!= channels->cend(); ++iCh)
  {
    if ( !((*iCh)->IsOKWithParentMass(parentMass)) ) continue;
    sumBR += (*iCh)->GetBR();
  }
  if (sumBR <= 0.0)
  {
#ifdef G4VERBOSE
    G4cout << " G4DecayTable::SelectADecayChannel :: no possible DecayChannel"
           << "       " << parent->GetParticleName() << G4endl;
#endif
    return nullptr;
  }
 
  const std::size_t MAX_LOOP = 10000;
  for (std::size_t loop_counter=0; loop_counter<MAX_LOOP; ++loop_counter)
  {
    G4double sum = 0.0;
    G4double br= sumBR*G4UniformRand();
    // select decay channel
    for (auto iCh = channels->cbegin(); iCh!= channels->cend(); ++iCh)
    {
      sum += (*iCh)->GetBR();
      if ( !((*iCh)->IsOKWithParentMass(parentMass)) ) continue;
      if (br < sum) return (*iCh);
    }
  }
  return nullptr;
}

Referenced by G4Decay::DecayIt(), G4IntraNucleiCascader::decayTrappedParticle(), and G4RadioactiveDecay::DoDecay().

---

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

• G4DecayTable.hh
• G4DecayTable.cc