#include <G4StepStatistics.h>

Inheritance diagram for G4StepStatistics:

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

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

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

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

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

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

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

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

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

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

Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)

bool	isMT () const

const SimWatcher &	operator= (const SimWatcher &)=delete

virtual void	registerConsumes (edm::ConsumesCollector)

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

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
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 125 of file G4StepStatistics.h.

Constructor & Destructor Documentation

◆ G4StepStatistics()

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

inline

Definition at line 136 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

◆ update() [1/9]

void G4StepStatistics::update ( const DDDWorld * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const DDDWorld *>.

Definition at line 148 of file G4StepStatistics.h.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

156 {

◆ update() [2/9]

void G4StepStatistics::update ( const BeginOfJob * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 149 of file G4StepStatistics.h.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

156 {

◆ update() [3/9]

void G4StepStatistics::update ( const BeginOfRun * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 150 of file G4StepStatistics.h.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

156 {

◆ update() [4/9]

void G4StepStatistics::update ( const BeginOfEvent * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 154 of file G4StepStatistics.h.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

156 {

◆ update() [5/9]

void G4StepStatistics::update ( const BeginOfTrack * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 155 of file G4StepStatistics.h.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

156 {

◆ update() [6/9]

void G4StepStatistics::update ( const G4Step * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 156 of file G4StepStatistics.h.

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

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

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

◆ update() [7/9]

void G4StepStatistics::update ( const EndOfRun * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun *>.

Definition at line 182 of file G4StepStatistics.h.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

188 {

◆ update() [8/9]

void G4StepStatistics::update ( const EndOfEvent * )

inlineoverridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 188 of file G4StepStatistics.h.

References gather_cfg::cout, G4StatsMap, G4StepFreq, G4StepTree, m_verbose, and PDGID.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

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