G4MoleculeCounter Class Reference

#include <G4MoleculeCounter.hh>

Public Types
typedef std::map< G4Molecule, NbMoleculeAgainstTime >	CounterMapType
typedef std::auto_ptr < std::vector< G4Molecule > >	RecordedMolecules

Public Member Functions
const NbMoleculeAgainstTime &	GetNbMoleculeAgainstTime (const G4Molecule &molecule)
RecordedMolecules	GetRecordedMolecules ()
virtual void	DontRegister (const G4MoleculeDefinition *)
virtual void	RegisterAll ()
void	Use (G4bool flag=true)
G4bool	InUse ()
void	SetVerbose (G4int)
G4int	GetVerbose ()
void	SetTimeSlice (double)
virtual void	ResetCounter ()

Static Public Member Functions
static void	DeleteInstance ()
static G4MoleculeCounter *	GetMoleculeCounter ()

Protected Member Functions
	G4MoleculeCounter ()
virtual	~G4MoleculeCounter ()
virtual void	AddAMoleculeAtTime (const G4Molecule &, G4double)
virtual void	RemoveAMoleculeAtTime (const G4Molecule &, G4double)

Protected Attributes
CounterMapType	fCounterMap
std::map< const G4MoleculeDefinition *, G4bool >	fDontRegister
G4int	fVerbose

Static Protected Attributes
static G4ThreadLocal G4MoleculeCounter *	fpInstance = 0
static G4bool	fUse = FALSE

Friends
class	G4Molecule

Detailed Description

Definition at line 52 of file G4MoleculeCounter.hh.

Member Typedef Documentation

typedef std::map<G4Molecule, NbMoleculeAgainstTime> G4MoleculeCounter::CounterMapType

Definition at line 55 of file G4MoleculeCounter.hh.

typedef std::auto_ptr<std::vector<G4Molecule> > G4MoleculeCounter::RecordedMolecules

Definition at line 60 of file G4MoleculeCounter.hh.

Constructor & Destructor Documentation

G4MoleculeCounter::G4MoleculeCounter ( )

protected

Definition at line 38 of file G4MoleculeCounter.cc.

{
    fVerbose = 0 ;
}

virtual G4MoleculeCounter::~G4MoleculeCounter ( )

inlineprotectedvirtual

Definition at line 65 of file G4MoleculeCounter.hh.

{;}

Member Function Documentation

void G4MoleculeCounter::AddAMoleculeAtTime ( const G4Molecule &  molecule, G4double  time  )

protectedvirtual

Definition at line 67 of file G4MoleculeCounter.cc.

References G4BestUnit, G4cout, G4endl, G4Molecule::GetDefinition(), and G4Molecule::GetName().

Referenced by G4Molecule::BuildTrack().

{
    if(fDontRegister[molecule.GetDefinition()]) return ;

    if(fVerbose)
    {
        G4cout<<"G4MoleculeCounter::AddAMoleculeAtTime : "<< molecule.GetName()
              << " at time : " << G4BestUnit(time, "Time") <<G4endl;
    }

    CounterMapType::iterator counterMap_i = fCounterMap.find(molecule) ;

    if(counterMap_i == fCounterMap.end())
    {
        // DEBUG
//        if(fVerbose)  G4cout << " !! ***** Map is empty " << G4endl;
        fCounterMap[molecule][time] = 1;
    }
    else if(counterMap_i->second.empty())
    {
        // DEBUG
//        if(fVerbose)  G4cout << " !! ***** Map is empty " << G4endl;
        counterMap_i->second[time] = 1;
    }
    else
    {
        NbMoleculeAgainstTime::iterator end = counterMap_i->second.end();
        end--;

        // DEBUG
//        if(fVerbose)
//            G4cout<<"!! End Time = "<< G4BestUnit(end->first, "Time") <<G4endl;

        if(end->first <= time)
        {
            counterMap_i->second[time]=end->second + 1;
        }
        else
        {
            NbMoleculeAgainstTime::iterator it = counterMap_i->second.lower_bound(time);

            while(it->first > time && it!=counterMap_i->second.begin())
            {
                // DEBUG
//                if(fVerbose)
//                    G4cout<<"!!  ********** Is going back!!!!"<<G4endl;
                it--;
            }

            if(it==counterMap_i->second.begin() && it->first > time)
            {
                // DEBUG
//                if(fVerbose)
//                    G4cout<<"!!  ********** Illegal !!!!"<<G4endl;
                return ;
            }

            // DEBUG
//            if(fVerbose)
//            {
//                G4cout<<"!! PREVIOUS NB = "<< it->second <<G4endl;
//                G4cout<<"!! PREVIOUS TIME = "<< G4BestUnit(it->first,"Time") <<G4endl;
//            }
            counterMap_i->second[time]=it->second + 1;
        }
    }

    // DEBUG
//    if(fVerbose)
//        G4cout<<"!! NB = "<< fCounterMap[molecule][time]<<G4endl;
}

void G4MoleculeCounter::DeleteInstance ( )

static

Definition at line 53 of file G4MoleculeCounter.cc.

{
    if(fpInstance)
    {
        delete fpInstance;
        fpInstance = 0;
    }
}

void G4MoleculeCounter::DontRegister ( const G4MoleculeDefinition *  molDef )

inlinevirtual

Definition at line 137 of file G4MoleculeCounter.hh.

References fDontRegister.

{
    fDontRegister[molDef] = true ;
}

G4MoleculeCounter * G4MoleculeCounter::GetMoleculeCounter ( )

static

Definition at line 45 of file G4MoleculeCounter.cc.

Referenced by G4Molecule::BuildTrack(), and G4Molecule::~G4Molecule().

{
    if(!fpInstance)
        fpInstance = new G4MoleculeCounter();

    return fpInstance;
}

const NbMoleculeAgainstTime & G4MoleculeCounter::GetNbMoleculeAgainstTime ( const G4Molecule &  molecule )

inline

Definition at line 122 of file G4MoleculeCounter.hh.

References fCounterMap.

{
    return fCounterMap[molecule];
}

G4MoleculeCounter::RecordedMolecules G4MoleculeCounter::GetRecordedMolecules

(

)

Definition at line 253 of file G4MoleculeCounter.cc.

References G4cout, and G4endl.

{
    if(fVerbose > 1)
    {
        G4cout<<"Entering in G4MoleculeCounter::RecordMolecules"<<G4endl;
    }

    CounterMapType::iterator it;
    RecordedMolecules output (new vector<G4Molecule>) ;

    for(it = fCounterMap.begin() ; it != fCounterMap.end() ; it++)
    {
        output->push_back(it->first);
    }
    return output;
}

G4int G4MoleculeCounter::GetVerbose ( )

inline

Definition at line 132 of file G4MoleculeCounter.hh.

References fVerbose.

{
    return fVerbose ;
}

G4bool G4MoleculeCounter::InUse ( )

inline

Definition at line 101 of file G4MoleculeCounter.hh.

References fUse.

Referenced by G4Molecule::BuildTrack().

    {
        return fUse;
    }

void G4MoleculeCounter::RegisterAll ( )

inlinevirtual

Definition at line 142 of file G4MoleculeCounter.hh.

References fDontRegister.

{
    fDontRegister.clear();
}

void G4MoleculeCounter::RemoveAMoleculeAtTime ( const G4Molecule &  molecule, G4double  time  )

protectedvirtual

Definition at line 139 of file G4MoleculeCounter.cc.

References FatalErrorInArgument, compDoubleWithPrecision::fPrecision, G4BestUnit, G4cout, G4endl, G4Exception(), G4Molecule::GetDefinition(), G4Molecule::GetName(), and G4Molecule::PrintState().

Referenced by G4Molecule::~G4Molecule().

{
    if(fDontRegister[molecule.GetDefinition()]) return ;

    if(fVerbose)
    {
        G4cout<<"G4MoleculeCounter::RemoveAMoleculeAtTime : "<< molecule.GetName()
             << " at time : " << G4BestUnit(time,"Time") <<G4endl;
    }

    NbMoleculeAgainstTime& nbMolPerTime = fCounterMap[molecule];

    if(nbMolPerTime.empty())
    {
        molecule.PrintState();
        G4String errMsg = "You are trying to remove molecule "
                          + molecule.GetName()
                          +" from the counter while this kind of molecules has not been registered yet";
        G4Exception("G4MoleculeCounter::RemoveAMoleculeAtTime","",FatalErrorInArgument, errMsg);

        return;
    }
    else
    {
        NbMoleculeAgainstTime::iterator it ;

        if(nbMolPerTime.size() == 1)
        {
            it = nbMolPerTime.begin() ;
            // DEBUG
//            if(fVerbose)
//                G4cout << "!! fCounterMap[molecule].size() == 1" << G4endl;
        }
        else
        {
            it = nbMolPerTime.lower_bound(time);
        }

        if(it==nbMolPerTime.end())
        {
            // DEBUG
//            if(fVerbose)
//                G4cout << " ********** NO ITERATOR !!!!!!!!! " << G4endl;
            it--;

            if(time<it->first)
            {
                G4String errMsg = "There was no "+ molecule.GetName()
                        + " record at the time or even before the time asked";
                G4Exception("G4MoleculeCounter::RemoveAMoleculeAtTime","",FatalErrorInArgument, errMsg);
            }
        }

        // DEBUG
//        if(fVerbose)
//        {
////            G4cout << "G4MoleculeCounter::RemoveAMoleculeAtTime " << G4endl;
//            G4cout<<"!! Molecule = " << molecule.GetName() << G4endl;
//            G4cout<<"!! At Time = "<< G4BestUnit(time,"Time") <<G4endl;
//            G4cout<<"!! PREVIOUS TIME = "<< G4BestUnit(it->first,"Time")<<G4endl;
//            G4cout<<"!! PREVIOUS Nb = "<< it->second <<G4endl;
//        }

        // If valgrind problem on the line below, it means that the pointer "it"
        // points nowhere
        if(nbMolPerTime.value_comp()(*it, *nbMolPerTime.begin()))
        {
            // DEBUG
//            if(fVerbose)
//                G4cout<<"!! ***** In value_comp ... " << G4endl;
            it++;
            if(time<it->first)
            {
                G4String errMsg = "There was no "+ molecule.GetName()
                        + " record at the time or even before the time asked";
                G4Exception("G4MoleculeCounter::RemoveAMoleculeAtTime","",FatalErrorInArgument, errMsg);
            }
        }

        while(it->first - time > compDoubleWithPrecision::fPrecision  && it!=nbMolPerTime.begin())
        {
            // DEBUG
//            if(fVerbose)
//            {
//                G4cout<<"!! ***** Is going back!!!!"<<G4endl;
//                G4cout<<"!! PREVIOUS TIME = "<< G4BestUnit(it-> first,"Time") <<G4endl;
//            }
            it--;
        }

        if(it==nbMolPerTime.begin() && it->first > time)
        {
            // DEBUG
//            if(fVerbose)
//                G4cout<<"!!  ********** Illegal !!!!"<<G4endl;
            return ;
        }

        // DEBUG
//        if(fVerbose)
//        {
//            G4cout<<"!! PREVIOUS NB = "<< (*it).second <<G4endl;
//            G4cout<<"!! PREVIOUS TIME = "<< G4BestUnit(it->first,"Time")<<G4endl;
//        }
        nbMolPerTime[time]=it->second - 1;
    }

    // DEBUG
//    if(fVerbose)
//    {
//        G4cout<<"!! NB = "<< nbMolPerTime[time]<<G4endl;
//    }
}

void G4MoleculeCounter::ResetCounter ( )

inlinevirtual

Definition at line 117 of file G4MoleculeCounter.hh.

References fCounterMap.

{
    fCounterMap.clear();
}

void G4MoleculeCounter::SetTimeSlice

(

double

timeSlice

)

Definition at line 62 of file G4MoleculeCounter.cc.

References compDoubleWithPrecision::fPrecision.

{
    compDoubleWithPrecision::fPrecision = timeSlice;
}

void G4MoleculeCounter::SetVerbose ( G4int  level )

inline

Definition at line 127 of file G4MoleculeCounter.hh.

References fVerbose.

{
    fVerbose = level;
}

void G4MoleculeCounter::Use ( G4bool  flag = true )

inline

Definition at line 97 of file G4MoleculeCounter.hh.

References fUse.

    {
        fUse=flag;
    }

Friends And Related Function Documentation

friend class G4Molecule

friend

Definition at line 74 of file G4MoleculeCounter.hh.

Field Documentation

CounterMapType G4MoleculeCounter::fCounterMap

protected

Definition at line 68 of file G4MoleculeCounter.hh.

Referenced by GetNbMoleculeAgainstTime(), and ResetCounter().

std::map<const G4MoleculeDefinition*, G4bool> G4MoleculeCounter::fDontRegister

protected

Definition at line 69 of file G4MoleculeCounter.hh.

Referenced by DontRegister(), and RegisterAll().

G4ThreadLocal G4MoleculeCounter * G4MoleculeCounter::fpInstance = 0

staticprotected

Definition at line 66 of file G4MoleculeCounter.hh.

G4bool G4MoleculeCounter::fUse = FALSE

staticprotected

Definition at line 70 of file G4MoleculeCounter.hh.

Referenced by InUse(), and Use().

G4int G4MoleculeCounter::fVerbose

protected

Definition at line 72 of file G4MoleculeCounter.hh.

Referenced by GetVerbose(), and SetVerbose().

---

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

• G4MoleculeCounter.hh
• G4MoleculeCounter.cc