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 48 of file G4DecayTable.hh.

Constructor & Destructor Documentation

G4DecayTable() [1/2]

G4DecayTable::G4DecayTable

(

)

Definition at line 36 of file G4DecayTable.cc.

{
  channels = new G4VDecayChannelVector;
}

~G4DecayTable()

G4DecayTable::~G4DecayTable

(

)

Definition at line 41 of file G4DecayTable.cc.

{
  // remove and delete all contents
  for (const auto channel : *channels) {
    delete channel;
  }
  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 111 of file G4DecayTable.cc.

{
  G4cout << "G4DecayTable:  " << parent->GetParticleName() << G4endl;
  G4int index = 0;
  for (const auto channel : *channels) {
    G4cout << index << ": ";
    channel->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 95 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 105 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(), and G4RadioactiveDecay::LoadDecayTable().

operator!=()

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

inline

Definition at line 90 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 85 of file G4DecayTable.hh.

{
  return (this == &right);
}

operator[]()

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

inline

Definition at line 100 of file G4DecayTable.hh.

{
  return (*channels)[index];
}

SelectADecayChannel()

G4VDecayChannel * G4DecayTable::SelectADecayChannel

(

G4double

parentMass = -1.

)

Definition at line 77 of file G4DecayTable.cc.

{
  // check if contents exist
  if (channels->empty()) return nullptr;
 
  if (parentMass < 0.) parentMass = parent->GetPDGMass();
 
  G4double sumBR = 0.;
  for (const auto channel : *channels) {
    if (!(channel->IsOKWithParentMass(parentMass))) continue;
    sumBR += channel->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 (const auto channel : *channels) {
      sum += channel->GetBR();
      if (!(channel->IsOKWithParentMass(parentMass))) continue;
      if (br < sum) return channel;
    }
  }
  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