#include <G4MolecularDecayTable.hh>

Public Member Functions
	G4MolecularDecayTable ()

	~G4MolecularDecayTable ()

	G4MolecularDecayTable (const G4MolecularDecayTable &)

G4MolecularDecayTable &	operator= (const G4MolecularDecayTable &right)

void	AddExcitedState (const G4String &)

void	AddeConfToExcitedState (const G4String &, const G4ElectronOccupancy &)

void	AddDecayChannel (const G4String &, const G4MolecularDecayChannel *)

void	CheckDataConsistency ()

const std::vector< const G4MolecularDecayChannel * > *	GetDecayChannels (const G4ElectronOccupancy *) const

const std::vector< const G4MolecularDecayChannel * > *	GetDecayChannels (const G4String &) const

const G4String &	GetExcitedState (const G4ElectronOccupancy *) const

const G4ElectronOccupancy &	GetElectronOccupancy (const G4String &) const

const statesMap &	GetExcitedStateMaps () const

const channelsMap &	GetDecayChannelsMap () const

Detailed Description

Class Description G4MolecularDecayTable operates as a container of deexcitation modes for excited or ionized molecules

Definition at line 88 of file G4MolecularDecayTable.hh.

Constructor & Destructor Documentation

◆ G4MolecularDecayTable() [1/2]

G4MolecularDecayTable::G4MolecularDecayTable ( )

Definition at line 43 of file G4MolecularDecayTable.cc.

44{;}

◆ ~G4MolecularDecayTable()

G4MolecularDecayTable::~G4MolecularDecayTable ( )

Definition at line 46 of file G4MolecularDecayTable.cc.

{
    channelsMap::iterator it_map = fDecayChannelsMap.begin();
 
    for(;it_map != fDecayChannelsMap.end() ; it_map++)
    {
        vector<const G4MolecularDecayChannel*>& decayChannels = it_map->second;
        if(!decayChannels.empty())
        {
            for(int i = 0 ; i < (int) decayChannels.size() ; i++)
            {
                if(decayChannels [i])
                {
                    delete decayChannels[i];
                    decayChannels[i] = 0;
                }
            }
            decayChannels.clear();
        }
    }
    fDecayChannelsMap.clear();
}

◆ G4MolecularDecayTable() [2/2]

G4MolecularDecayTable::G4MolecularDecayTable ( const G4MolecularDecayTable & right )

Definition at line 69 of file G4MolecularDecayTable.cc.

{
    *this = right;
}

Member Function Documentation

◆ AddDecayChannel()

void G4MolecularDecayTable::AddDecayChannel	(	const G4String &	label,
		const G4MolecularDecayChannel *	channel
	)

Definition at line 158 of file G4MolecularDecayTable.cc.

{
    fDecayChannelsMap[label].push_back(channel);
}

Referenced by G4MoleculeDefinition::AddDecayChannel().

◆ AddeConfToExcitedState()

void G4MolecularDecayTable::AddeConfToExcitedState	(	const G4String &	label,
		const G4ElectronOccupancy &	conf
	)

Definition at line 143 of file G4MolecularDecayTable.cc.

{
    statesMap::iterator statesIter = fExcitedStatesMap.find(conf);
 
    if (statesIter == fExcitedStatesMap.end())
    {
        fExcitedStatesMap[conf] = label;
    }
    else
    {
        G4Exception("G4MolecularDecayTable::AddExcitedState","G4MolecularDecayTable004",
                    FatalErrorInArgument,"Electronic configuration already registered in the decay table");
    }
}

Referenced by G4MoleculeDefinition::AddeConfToExcitedState().

◆ AddExcitedState()

void G4MolecularDecayTable::AddExcitedState ( const G4String & label )

Definition at line 130 of file G4MolecularDecayTable.cc.

{
    channelsMap::iterator channelsIter = fDecayChannelsMap.find(label);
    if(channelsIter != fDecayChannelsMap.end())
    {
        G4String errMsg = "Excited state" + label + " already registered in the decay table.";
        G4Exception("G4MolecularDecayTable::AddExcitedState",
                    "G4MolecularDecayTable003",FatalErrorInArgument, errMsg);
        return;
    }
    fDecayChannelsMap[label] ;
}

Referenced by G4MoleculeDefinition::AddExcitedState().

◆ CheckDataConsistency()

void G4MolecularDecayTable::CheckDataConsistency ( )

Definition at line 163 of file G4MolecularDecayTable.cc.

{
    channelsMap::iterator channelsIter;
 
    //Let's check probabilities
 
    for (channelsIter=fDecayChannelsMap.begin(); channelsIter!=fDecayChannelsMap.end(); channelsIter++ )
    {
 
        vector<const G4MolecularDecayChannel*>& decayVect = channelsIter->second;
        G4double sum=0;
 
        G4double max = decayVect.size();
 
        for (size_t i=0; i<max; i++)
        {
            const G4MolecularDecayChannel* decay = decayVect[i];
            const G4double prob = decay->GetProbability();
            sum += prob;
        }
 
        if (sum != 1)
        {
            G4String errMsg = "Deexcitation Channels probabilities in " + channelsIter->first
                    + "excited state don't sum up to 1";
            G4Exception("G4MolecularDecayTable::CheckDataConsistency",
                        "G4MolecularDecayTable005",FatalErrorInArgument, errMsg);
        }
    }
 
}

◆ GetDecayChannels() [1/2]

const vector< const G4MolecularDecayChannel * > * G4MolecularDecayTable::GetDecayChannels ( const G4ElectronOccupancy * conf ) const

Definition at line 81 of file G4MolecularDecayTable.cc.

{
    statesMap::const_iterator it_exstates  = fExcitedStatesMap.find(*conf);
    if(it_exstates == fExcitedStatesMap.end()) return 0;
    channelsMap::const_iterator it_decchannel =  fDecayChannelsMap.find(it_exstates->second) ;
    if(it_decchannel == fDecayChannelsMap.end()) return 0;
    return &(it_decchannel->second);
}

Referenced by G4MoleculeDefinition::GetDecayChannels().

◆ GetDecayChannels() [2/2]

const vector< const G4MolecularDecayChannel * > * G4MolecularDecayTable::GetDecayChannels ( const G4String & exState ) const

Definition at line 90 of file G4MolecularDecayTable.cc.

{
    channelsMap::const_iterator it_decchannel = fDecayChannelsMap.find(exState);
    if(it_decchannel == fDecayChannelsMap.end()) return 0;
    return &(it_decchannel->second);
}

◆ GetDecayChannelsMap()

const channelsMap & G4MolecularDecayTable::GetDecayChannelsMap ( ) const

inline

Definition at line 134 of file G4MolecularDecayTable.hh.

{
    return fDecayChannelsMap;
}

Referenced by operator=().

◆ GetElectronOccupancy()

const G4ElectronOccupancy & G4MolecularDecayTable::GetElectronOccupancy ( const G4String & exState ) const

Definition at line 112 of file G4MolecularDecayTable.cc.

{
    statesMap::const_iterator statesIter;
    const G4ElectronOccupancy* conf(0);
    for (statesIter=fExcitedStatesMap.begin(); statesIter!=fExcitedStatesMap.end(); statesIter++ )
    {
        if(exState == statesIter->second) conf = &(statesIter->first);
    }
 
    if(statesIter == fExcitedStatesMap.end())
    {
        G4String errMsg = "Excited state" + exState + " not found";
        G4Exception("G4MolecularDecayTable::GetElectronOccupancy(const G4String&)",
                    "G4MolecularDecayTable002",FatalErrorInArgument, errMsg);
    }
    return *conf;
}

◆ GetExcitedState()

const G4String & G4MolecularDecayTable::GetExcitedState ( const G4ElectronOccupancy * conf ) const

Definition at line 97 of file G4MolecularDecayTable.cc.

{
    statesMap::const_iterator it_exstates  = fExcitedStatesMap.find(*conf);
 
    if(it_exstates == fExcitedStatesMap.end())
    {
        G4String errMsg = "Excited state not found";
        G4Exception("G4MolecularDecayTable::GetExcitedState(const G4ElectronOccupancy*)",
                    "G4MolecularDecayTable001",FatalErrorInArgument, errMsg);
//        return *(new G4String("IM FAKE"));  // fake return statement
    }
 
    return it_exstates->second;
}

Referenced by G4MoleculeDefinition::GetExcitedState().

◆ GetExcitedStateMaps()

const statesMap & G4MolecularDecayTable::GetExcitedStateMaps ( ) const