Inheritance diagram for DoCastorAnalysis:

Public Member Functions
	DoCastorAnalysis (const edm::ParameterSet &p)

	~DoCastorAnalysis () override

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 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 EndOfRun *>
	Observer ()

void	slotForUpdate (const EndOfRun * 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 EndOfEvent *>
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

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

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

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

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

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

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

Private Attributes
TFile *	CastorOutputEventFile

TTree *	CastorTree

int	eventIndex

std::vector< double > *	psimhit_energy

std::vector< double > *	psimhit_eta

std::vector< int > *	psimhit_module

std::vector< double > *	psimhit_phi

std::vector< int > *	psimhit_sector

std::vector< double > *	psimhit_x

std::vector< double > *	psimhit_y

std::vector< double > *	psimhit_z

std::vector< double >	simhit_energy

std::vector< double >	simhit_eta

double	simhit_etot

std::vector< int >	simhit_module

std::vector< double >	simhit_phi

std::vector< int >	simhit_sector

std::vector< double >	simhit_x

std::vector< double >	simhit_y

std::vector< double >	simhit_z

std::string	TreeFileName

int	verbosity

Additional Inherited Members
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 67 of file DoCastorAnalysis.cc.

Constructor & Destructor Documentation

◆ DoCastorAnalysis()

DoCastorAnalysis::DoCastorAnalysis ( const edm::ParameterSet & p )

Definition at line 108 of file DoCastorAnalysis.cc.

References CastorOutputEventFile, CastorTree, edm::ParameterSet::getParameter(), AlCaHLTBitMon_ParallelJobs::p, psimhit_energy, psimhit_eta, psimhit_module, psimhit_phi, psimhit_sector, psimhit_x, psimhit_y, psimhit_z, simhit_etot, AlCaHLTBitMon_QueryRunRegistry::string, TreeFileName, and verbosity.

                                                            {
   edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("DoCastorAnalysis");
   verbosity = m_Anal.getParameter<int>("Verbosity");
 
   TreeFileName = m_Anal.getParameter<std::string>("CastorTreeFileName");
 
   if (verbosity > 0) {
     edm::LogVerbatim("ForwardSim") << std::endl;
     edm::LogVerbatim("ForwardSim") << "============================================================================";
     edm::LogVerbatim("ForwardSim") << "DoCastorAnalysis:: Initialized as observer";
 
     edm::LogVerbatim("ForwardSim") << " Castor Tree will be created";
     edm::LogVerbatim("ForwardSim") << " Castor Tree will be in file: " << TreeFileName;
 #ifdef EDM_ML_DEBUG
     getchar();
 #endif
     edm::LogVerbatim("ForwardSim") << "============================================================================";
     edm::LogVerbatim("ForwardSim") << std::endl;
   }
 
   edm::LogVerbatim("ForwardSim") << "DoCastorAnalysis: output event root file created";
   TString treefilename = TreeFileName;
   CastorOutputEventFile = new TFile(treefilename, "RECREATE");
 
   CastorTree = new TTree("Sim", "Sim");
 
   CastorTree->Branch("simhit_x", "std::vector<double>", &psimhit_x);
   CastorTree->Branch("simhit_y", "std::vector<double>", &psimhit_y);
   CastorTree->Branch("simhit_z", "std::vector<double>", &psimhit_z);
 
   CastorTree->Branch("simhit_eta", "std::vector<double>", &psimhit_eta);
   CastorTree->Branch("simhit_phi", "std::vector<double>", &psimhit_phi);
   CastorTree->Branch("simhit_energy", "std::vector<double>", &psimhit_energy);
 
   CastorTree->Branch("simhit_sector", "std::vector<int>", &psimhit_sector);
   CastorTree->Branch("simhit_module", "std::vector<int>", &psimhit_module);
 
   CastorTree->Branch("simhit_etot", &simhit_etot, "simhit_etot/D");
 }

◆ ~DoCastorAnalysis()

DoCastorAnalysis::~DoCastorAnalysis ( )

override

Definition at line 148 of file DoCastorAnalysis.cc.

References CastorOutputEventFile, CastorTree, and verbosity.

                                     {
   //destructor of DoCastorAnalysis
 
   CastorOutputEventFile->cd();
   //-- CastorOutputEventFile->Write();
   CastorTree->Write("", TObject::kOverwrite);
   edm::LogVerbatim("ForwardSim") << "DoCastorAnalysis: Ntuple event written";
 #ifdef EDM_ML_DEBUG
   getchar();
 #endif
   CastorOutputEventFile->Close();
   edm::LogVerbatim("ForwardSim") << "DoCastorAnalysis: Event file closed";
 #ifdef EDM_ML_DEBUG
   getchar();
 #endif
 
   if (verbosity > 0) {
     edm::LogVerbatim("ForwardSim") << std::endl << "DoCastorAnalysis: end of process";
 #ifdef EDM_ML_DEBUG
     getchar();
 #endif
   }
 }

Member Function Documentation

◆ update() [1/6]

void DoCastorAnalysis::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 174 of file DoCastorAnalysis.cc.

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

174 { edm::LogVerbatim("ForwardSim") << " Starting new job "; }

edm::LogVerbatim

Log< level::Info, true > LogVerbatim

Definition: MessageLogger.h:134

◆ update() [2/6]

void DoCastorAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 178 of file DoCastorAnalysis.cc.

References eventIndex.

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

                                                    {
   edm::LogVerbatim("ForwardSim") << std::endl << "DoCastorAnalysis: Starting Run";
 
   // edm::LogVerbatim("ForwardSim") << "DoCastorAnalysis: output event root file created";
   // TString treefilename = TreeFileName;
   // CastorOutputEventFile = new TFile(treefilename,"RECREATE");
 
   eventIndex = 1;
 }

◆ update() [3/6]

void DoCastorAnalysis::update ( const EndOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun *>.

Definition at line 297 of file DoCastorAnalysis.cc.

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

297 { ; }

◆ update() [4/6]

void DoCastorAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 188 of file DoCastorAnalysis.cc.

References eventIndex.

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

                                                      {
   edm::LogVerbatim("ForwardSim") << "DoCastorAnalysis: Processing Event Number: " << eventIndex;
   eventIndex++;
 }

◆ update() [5/6]

void DoCastorAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 195 of file DoCastorAnalysis.cc.

References CastorTree, HBHEDarkening_cff::energy, PVValHelper::eta, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getUnitID(), CrabHelper::log, phi, psimhit_energy, psimhit_eta, psimhit_module, psimhit_phi, psimhit_sector, psimhit_x, psimhit_y, psimhit_z, nano_mu_digi_cff::sector, simhit_energy, simhit_eta, simhit_etot, simhit_module, simhit_phi, simhit_sector, simhit_x, simhit_y, simhit_z, funct::tan(), tauSpinnerTable_cfi::theta, CastorNumberingScheme::unpackIndex(), and ecaldqm::zside().

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

                                                    {
   // Look for the Hit Collection
 
   // access to the G4 hit collections
   G4HCofThisEvent *allHC = (*evt)()->GetHCofThisEvent();
 
   int CAFIid = G4SDManager::GetSDMpointer()->GetCollectionID("CastorFI");
   CaloG4HitCollection *theCAFI = (CaloG4HitCollection *)allHC->GetHC(CAFIid);
 
   CastorNumberingScheme *theCastorNumScheme = new CastorNumberingScheme();
 
   unsigned int volumeID = 0;
   // std::map<int,float,std::less<int> > themap;
 
   int nentries = theCAFI->entries();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "nentries in CAFI: " << nentries;
   getchar();
 #endif
 
   psimhit_x = &simhit_x;
   psimhit_x->clear();
   psimhit_x->reserve(nentries);
 
   psimhit_y = &simhit_y;
   psimhit_y->clear();
   psimhit_y->reserve(nentries);
 
   psimhit_z = &simhit_z;
   psimhit_z->clear();
   psimhit_z->reserve(nentries);
 
   psimhit_eta = &simhit_eta;
   psimhit_eta->clear();
   psimhit_eta->reserve(nentries);
 
   psimhit_phi = &simhit_phi;
   psimhit_phi->clear();
   psimhit_phi->reserve(nentries);
 
   psimhit_energy = &simhit_energy;
   psimhit_energy->clear();
   psimhit_energy->reserve(nentries);
 
   psimhit_sector = &simhit_sector;
   psimhit_sector->clear();
   psimhit_sector->reserve(nentries);
 
   psimhit_module = &simhit_module;
   psimhit_module->clear();
   psimhit_module->reserve(nentries);
 
   simhit_etot = 0;
 
   if (nentries > 0) {
     for (int ihit = 0; ihit < nentries; ihit++) {
       CaloG4Hit *aHit = (*theCAFI)[ihit];
       volumeID = aHit->getUnitID();
 
       //themap[volumeID] += aHit->getEnergyDeposit();
       int zside, sector, zmodule;
 
       theCastorNumScheme->unpackIndex(volumeID, zside, sector, zmodule);
 
       double energy = aHit->getEnergyDeposit() / CLHEP::GeV;
       //double time     = aHit->getTimeSlice();
 
       math::XYZPoint pos = aHit->getPosition();
       double theta = pos.theta();
       double eta = -log(tan(theta / 2.));
       double phi = pos.phi();
 
       psimhit_x->push_back(pos.x());
       psimhit_y->push_back(pos.y());
       psimhit_z->push_back(pos.z());
 
       psimhit_eta->push_back(eta);
       psimhit_phi->push_back(phi);
       psimhit_energy->push_back(energy);
 
       psimhit_sector->push_back(sector);
       psimhit_module->push_back(zmodule);
 
       simhit_etot += energy;
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("ForwardSim") << "hit " << ihit + 1 << " : x = " << (*psimhit_x)[ihit]
                                      << " , eta =  " << (*psimhit_eta)[ihit] << " , phi = " << (*psimhit_phi)[ihit]
                                      << " , energy = " << (*psimhit_energy)[ihit];
 #endif
     }
 
     //if (debug) edm::LogVerbatim("ForwardSim") << " total energy = " << simhit_etot;
 #ifdef EDM_ML_DEBUG
     getchar();
 #endif
     CastorTree->Fill();
 
   }  // nentries > 0
   delete theCastorNumScheme;
 }

◆ update() [6/6]

void DoCastorAnalysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 299 of file DoCastorAnalysis.cc.

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

299 { ; }

Member Data Documentation

◆ CastorOutputEventFile

TFile* DoCastorAnalysis::CastorOutputEventFile

private

Definition at line 92 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and ~DoCastorAnalysis().

◆ CastorTree

TTree* DoCastorAnalysis::CastorTree

private

Definition at line 93 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), update(), and ~DoCastorAnalysis().

◆ eventIndex

int DoCastorAnalysis::eventIndex

private

Definition at line 95 of file DoCastorAnalysis.cc.

Referenced by update().

◆ psimhit_energy

std::vector<double> * DoCastorAnalysis::psimhit_energy

private

Definition at line 102 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_eta

std::vector<double>* DoCastorAnalysis::psimhit_eta

private

Definition at line 102 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_module

std::vector<int> * DoCastorAnalysis::psimhit_module

private

Definition at line 103 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_phi

std::vector<double> * DoCastorAnalysis::psimhit_phi

private

Definition at line 102 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_sector

std::vector<int>* DoCastorAnalysis::psimhit_sector

private

Definition at line 103 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_x

std::vector<double>* DoCastorAnalysis::psimhit_x

private

Definition at line 101 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_y

std::vector<double> * DoCastorAnalysis::psimhit_y

private

Definition at line 101 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_z

std::vector<double> * DoCastorAnalysis::psimhit_z

private

Definition at line 101 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ simhit_energy

std::vector<double> DoCastorAnalysis::simhit_energy

private

Definition at line 98 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_eta

std::vector<double> DoCastorAnalysis::simhit_eta

private

Definition at line 98 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_etot

double DoCastorAnalysis::simhit_etot

private

Definition at line 105 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and update().

◆ simhit_module

std::vector<int> DoCastorAnalysis::simhit_module

private

Definition at line 99 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_phi

std::vector<double> DoCastorAnalysis::simhit_phi

private

Definition at line 98 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_sector

std::vector<int> DoCastorAnalysis::simhit_sector

private

Definition at line 99 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_x

std::vector<double> DoCastorAnalysis::simhit_x

private

Definition at line 97 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_y

std::vector<double> DoCastorAnalysis::simhit_y

private

Definition at line 97 of file DoCastorAnalysis.cc.

Referenced by update().

◆ simhit_z

std::vector<double> DoCastorAnalysis::simhit_z

private

Definition at line 97 of file DoCastorAnalysis.cc.

Referenced by update().

◆ TreeFileName

std::string DoCastorAnalysis::TreeFileName

private

Definition at line 90 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis().

◆ verbosity

int DoCastorAnalysis::verbosity

private

Definition at line 88 of file DoCastorAnalysis.cc.

Referenced by DoCastorAnalysis(), and ~DoCastorAnalysis().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DoCastorAnalysis()

◆ ~DoCastorAnalysis()

Member Function Documentation

◆ update() [1/6]

◆ update() [2/6]

◆ update() [3/6]

◆ update() [4/6]

◆ update() [5/6]

◆ update() [6/6]

Member Data Documentation

◆ CastorOutputEventFile

◆ CastorTree

◆ eventIndex

◆ psimhit_energy

◆ psimhit_eta

◆ psimhit_module

◆ psimhit_phi

◆ psimhit_sector

◆ psimhit_x

◆ psimhit_y

◆ psimhit_z

◆ simhit_energy

◆ simhit_eta

◆ simhit_etot

◆ simhit_module

◆ simhit_phi

◆ simhit_sector

◆ simhit_x

◆ simhit_y

◆ simhit_z

◆ TreeFileName

◆ verbosity