G4DecayTable Class Reference

#include <G4DecayTable.hh>

Public Types
typedef std::vector < G4VDecayChannel * >	G4VDecayChannelVector

Public Member Functions
	G4DecayTable ()
	~G4DecayTable ()
G4int	operator== (const G4DecayTable &right) const
G4int	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 50 of file G4DecayTable.hh.

Member Typedef Documentation

typedef std::vector<G4VDecayChannel*> G4DecayTable::G4VDecayChannelVector

Definition at line 58 of file G4DecayTable.hh.

Constructor & Destructor Documentation

G4DecayTable::G4DecayTable

(

)

Definition at line 43 of file G4DecayTable.cc.

                          :parent(0)
{
  channels =  new G4VDecayChannelVector;
}

G4DecayTable::~G4DecayTable

(

)

Definition at line 48 of file G4DecayTable.cc.

{
  // remove and delete all contents  
  G4VDecayChannelVector::iterator i;
  for (i = channels->begin(); i!= channels->end(); ++i) {
    delete (*i);
  }
  channels->clear();
  delete  channels;
  channels = 0;
}    

Member Function Documentation

void G4DecayTable::DumpInfo

(

)

const

Definition at line 102 of file G4DecayTable.cc.

{
  G4cout << "G4DecayTable:  " << parent->GetParticleName() << G4endl;
  G4int index =0;
  G4VDecayChannelVector::iterator i;
  for (i = channels->begin(); i!= channels->end(); ++i) {
    G4cout << index << ": ";
    (*i)->DumpInfo();
    index +=1;
  }
  G4cout << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4DecayTable::entries ( ) const

inline

Definition at line 100 of file G4DecayTable.hh.

{
  return channels->size();
}

Here is the caller graph for this function:

G4VDecayChannel * G4DecayTable::GetDecayChannel ( G4int  index ) const

inline

Definition at line 113 of file G4DecayTable.hh.

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

Here is the caller graph for this function:

void G4DecayTable::Insert

(

G4VDecayChannel *

aChannel

)

Definition at line 60 of file G4DecayTable.cc.

                                                    {
  if (parent == 0) { 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 r = aChannel->GetBR();
    G4VDecayChannelVector::iterator i;
    for (i = channels->begin(); i!= channels->end(); ++i) {
      if (r > (*i)->GetBR()) {
    channels->insert(i,aChannel);
    return;
      }
    }
    channels->push_back(aChannel);
  }
}

Here is the call graph for this function:

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

inline

Definition at line 74 of file G4DecayTable.hh.

{return (this != &right);};

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

inline

Definition at line 73 of file G4DecayTable.hh.

{return (this == &right);};

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

inline

Definition at line 106 of file G4DecayTable.hh.

{
  return (*channels)[index];
}

G4VDecayChannel * G4DecayTable::SelectADecayChannel

(

G4double

parentMass = -1.

)

Definition at line 81 of file G4DecayTable.cc.

{
  // check if contents exist
  if (channels->size()<1) return 0;

  if(parentMass<0.) parentMass=parent->GetPDGMass(); 
  const size_t MAX_LOOP = 10000;
  for (size_t loop_counter=0; loop_counter <MAX_LOOP; ++loop_counter){
    G4double sumBR = 0.0;
    G4double r= G4UniformRand();
    // select decay channel
    G4VDecayChannelVector::iterator i;
    for (i = channels->begin(); i!= channels->end(); ++i) {
      sumBR += (*i)->GetBR();
      if ( !((*i)->IsOKWithParentMass(parentMass)) ) continue;
      if (r < sumBR) return (*i);
    }
  }
  return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• geant4.10.03.p01/source/particles/management/include/G4DecayTable.hh
• geant4.10.03.p01/source/particles/management/src/G4DecayTable.cc