G4MolecularDecayTable Class Reference

#include <G4MolecularDissociationTable.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 G4MolecularDissociationTable.hh.

Constructor & Destructor Documentation

G4MolecularDecayTable::G4MolecularDecayTable

(

)

Definition at line 43 of file G4MolecularDissociationTable.cc.

{;}

G4MolecularDecayTable::~G4MolecularDecayTable

(

)

Definition at line 46 of file G4MolecularDissociationTable.cc.

References int().

{
    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::G4MolecularDecayTable

(

const G4MolecularDecayTable &

right

)

Definition at line 69 of file G4MolecularDissociationTable.cc.

References right.

{
    *this = right;
}

Member Function Documentation

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

Definition at line 158 of file G4MolecularDissociationTable.cc.

Referenced by G4MoleculeDefinition::AddDecayChannel().

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

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

Definition at line 143 of file G4MolecularDissociationTable.cc.

References FatalErrorInArgument, and G4Exception().

Referenced by G4MoleculeDefinition::AddeConfToExcitedState().

{
    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");
    }
}

void G4MolecularDecayTable::AddExcitedState

(

const G4String &

label

)

Definition at line 130 of file G4MolecularDissociationTable.cc.

References FatalErrorInArgument, and G4Exception().

Referenced by G4MoleculeDefinition::AddExcitedState().

{
    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] ;
}

void G4MolecularDecayTable::CheckDataConsistency

(

)

Definition at line 163 of file G4MolecularDissociationTable.cc.

References G4INCL::ClusterDecay::decay(), FatalErrorInArgument, G4String::first(), G4Exception(), G4MolecularDecayChannel::GetProbability(), and G4INCL::Math::max().

{
    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);
        }
    }

}

const vector< const G4MolecularDecayChannel * > * G4MolecularDecayTable::GetDecayChannels

(

const G4ElectronOccupancy *

conf

)

const

Definition at line 81 of file G4MolecularDissociationTable.cc.

Referenced by G4MoleculeDefinition::GetDecayChannels().

{
    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);
}

const vector< const G4MolecularDecayChannel * > * G4MolecularDecayTable::GetDecayChannels

(

const G4String &

exState

)

const

Definition at line 90 of file G4MolecularDissociationTable.cc.

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

const channelsMap & G4MolecularDecayTable::GetDecayChannelsMap ( ) const

inline

Definition at line 134 of file G4MolecularDissociationTable.hh.

Referenced by operator=().

{
    return fDecayChannelsMap;
}

const G4ElectronOccupancy & G4MolecularDecayTable::GetElectronOccupancy

(

const G4String &

exState

)

const

Definition at line 112 of file G4MolecularDissociationTable.cc.

References FatalErrorInArgument, and G4Exception().

{
    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;
}

const G4String & G4MolecularDecayTable::GetExcitedState

(

const G4ElectronOccupancy *

conf

)

const

Definition at line 97 of file G4MolecularDissociationTable.cc.

References FatalErrorInArgument, and G4Exception().

Referenced by G4MoleculeDefinition::GetExcitedState().

{
    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;
}

const statesMap & G4MolecularDecayTable::GetExcitedStateMaps ( ) const

inline

Definition at line 129 of file G4MolecularDissociationTable.hh.

{
    return fExcitedStatesMap;
}

G4MolecularDecayTable & G4MolecularDecayTable::operator=

(

const G4MolecularDecayTable &

right

)

Definition at line 74 of file G4MolecularDissociationTable.cc.

References GetDecayChannelsMap().

{
    fExcitedStatesMap = aMolecularDecayTable.fExcitedStatesMap;
    fDecayChannelsMap = channelsMap(aMolecularDecayTable.GetDecayChannelsMap());
    return *this;
}

---

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

• G4MolecularDissociationTable.hh
• G4MolecularDissociationTable.cc