#include <G4StepStatistics.h>

Inheritance diagram for G4StepStatistics:

Public Member Functions
	G4StepStatistics (const edm::ParameterSet &pSet)

void	update (const DDDWorld *)
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfJob *)
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfRun *)
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfEvent *)
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfTrack *)
	This routine will be called when the appropriate signal arrives. More...

void	update (const G4Step *iStep)
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfRun *)
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfEvent *iRun)
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfTrack *)
	This routine will be called when the appropriate signal arrives. More...

Public Member Functions inherited from SimWatcher
	SimWatcher ()

virtual	~SimWatcher ()

Public Member Functions inherited from Observer< const DDDWorld * >
	Observer ()

void	slotForUpdate (const DDDWorld *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()

void	slotForUpdate (const BeginOfJob *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()

void	slotForUpdate (const BeginOfRun *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()

void	slotForUpdate (const BeginOfEvent *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()

void	slotForUpdate (const BeginOfTrack *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const G4Step * >
	Observer ()

void	slotForUpdate (const G4Step *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfRun * >
	Observer ()

void	slotForUpdate (const EndOfRun *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent *iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfTrack * >
	Observer ()

void	slotForUpdate (const EndOfTrack *iT)

virtual	~Observer ()

Private Attributes
unsigned int	Event

edm::Service< TFileService >	fs

std::map< const StepID, unsigned int * >	G4StatsMap

Int_t	G4StepFreq [100000]

TTree *	G4StepTree

bool	m_verbose

Int_t	PDGID [100000]

TClonesArray *	Process

TClonesArray *	Region

Additional Inherited Members

Detailed Description

Definition at line 130 of file G4StepStatistics.h.

Constructor & Destructor Documentation

G4StepStatistics::G4StepStatistics ( const edm::ParameterSet & pSet )

inline

Definition at line 142 of file G4StepStatistics.h.

References fs, G4StepFreq, G4StepTree, TFileService::make(), PDGID, Process, and Region.

                                                  : 
      m_verbose(pSet.getUntrackedParameter<bool>("verbose",false)),
      Event(0)
      {
        //Adding TFile Service output
        G4StepTree = fs->make<TTree>("G4StepTree","G4Step Tree ");
        G4StepTree->Branch("Event",&Event,"Event/I");
        G4StepTree->Branch("PDGID",&PDGID,"PDGID[100000]/I");
        Region = new TClonesArray("TObjString",100000);
        G4StepTree->Branch("Region",&Region);
        Process = new TClonesArray("TObjString",100000);
        G4StepTree->Branch("Process",&Process);
        G4StepTree->Branch("G4StepFreq",&G4StepFreq,"G4StepFreq[100000]/I");
      }

Member Function Documentation

void G4StepStatistics::update ( const DDDWorld * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const DDDWorld * >.

Definition at line 156 of file G4StepStatistics.h.

164 {

void G4StepStatistics::update ( const BeginOfJob * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 157 of file G4StepStatistics.h.

164 {

void G4StepStatistics::update ( const BeginOfRun * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 158 of file G4StepStatistics.h.

164 {

void G4StepStatistics::update ( const BeginOfEvent * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 162 of file G4StepStatistics.h.

164 {

void G4StepStatistics::update ( const BeginOfTrack * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 163 of file G4StepStatistics.h.

164 {

void G4StepStatistics::update ( const G4Step * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 164 of file G4StepStatistics.h.

References gather_cfg::cout, G4StatsMap, StepID::GetParticlePDGID(), StepID::GetProcessName(), StepID::GetRegionName(), and m_verbose.

                                     { 
    std::cout <<"++ signal G4Step " ;
    //Dump the relevant information from the G4 step in the object mysteptest
    StepID mysteptest(iStep);
    //Add the StepID to the map, or increment if it already exists:
    if ( G4StatsMap.find(mysteptest) == G4StatsMap.end() )
      {
        //Allocating new memory for a pointer to associate with the key mysteptest
        //in our map. Initializing it to 1,will be incremented working on the value of the pointer.
        unsigned int* MyValue = new unsigned int(1);
        //Inserting the new key,value pair
        G4StatsMap.insert(std::make_pair(mysteptest, MyValue));
      }
    else
      {
        //Incrementing the value of the pointer by 1
        *G4StatsMap[mysteptest] = *G4StatsMap[mysteptest] + 1;
         }
    
    //If the verbose flag is set, then dump the information
    if (m_verbose) {
      std::cout << " StepID RegionName: "<< mysteptest.GetRegionName();
      std::cout << " StepID ProcessName: "<< mysteptest.GetProcessName();
      std::cout << " StepID ParticlePDGID: "<< mysteptest.GetParticlePDGID();
    }
    std::cout<<std::endl;
 }

void G4StepStatistics::update ( const EndOfRun * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun * >.

Definition at line 192 of file G4StepStatistics.h.

198 {

void G4StepStatistics::update ( const EndOfEvent * )

inlinevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 198 of file G4StepStatistics.h.

References gather_cfg::cout, G4StatsMap, G4StepFreq, G4StepTree, cmsHarvester::index, m_verbose, PDGID, and relval_parameters_module::step.

                                       {
   std::cout <<"++ signal EndOfEvent " <<std::endl;
   Event++;
   
   //Dumping the map in the log if verbose is chosen:
   if(m_verbose) {
     std::cout <<" G4StatsMap size is: "<<G4StatsMap.size()<<std::endl;
   }
   int index(0);
   for (std::map<const StepID,unsigned int*>::const_iterator step = G4StatsMap.begin(); step != G4StatsMap.end(); ++step, ++index){
     if(m_verbose) {
       std::cout <<" G4StatsMap step is: "<<step->first.GetRegionName()<<" "<<step->first.GetProcessName()<<" "<<step->first.GetParticlePDGID();//<<" "<<step->first.GetTrackID() ;
       std::cout <<" Number of such steps: "<< *step->second <<std::endl;
     }
     //Rolling the map into 5 "arrays", containing the StepID information and the G4Step statistics
     PDGID[index]=step->first.GetParticlePDGID();
     new ((*Region)[index]) TObjString (step->first.GetRegionName());
     new ((*Process)[index]) TObjString (step->first.GetProcessName());
     G4StepFreq[index]=*step->second;
   }
   
   G4StepTree->Fill();
 }