#include <DoCastorAnalysis.h>

Inheritance diagram for DoCastorAnalysis:

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

	~DoCastorAnalysis () override

Public Member Functions inherited from SimWatcher
	SimWatcher ()

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 BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...

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 EndOfEvent *evt) 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 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

Detailed Description

Definition at line 62 of file DoCastorAnalysis.h.

Constructor & Destructor Documentation

◆ DoCastorAnalysis()

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

Definition at line 28 of file DoCastorAnalysis.cc.

                                                            {
   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) {
     std::cout << std::endl;
     std::cout << "============================================================================" << std::endl;
     std::cout << "DoCastorAnalysis:: Initialized as observer" << std::endl;
  
     std::cout << " Castor Tree will be created" << std::endl;
     std::cout << " Castor Tree will be in file: " << TreeFileName << std::endl;
     if (debug)
       getchar();
  
     std::cout << "============================================================================" << std::endl;
     std::cout << std::endl;
   }
  
   std::cout << "DoCastorAnalysis: output event root file created" << std::endl;
   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_time","std::vector<double>",&psimhit_time);
   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");
 }

References CastorOutputEventFile, CastorTree, gather_cfg::cout, debug, 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.

◆ ~DoCastorAnalysis()

DoCastorAnalysis::~DoCastorAnalysis ( )

override

Definition at line 69 of file DoCastorAnalysis.cc.

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

References CastorOutputEventFile, CastorTree, gather_cfg::cout, debug, and verbosity.

Member Function Documentation

◆ update() [1/6]

void DoCastorAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 106 of file DoCastorAnalysis.cc.

                                                      {
   std::cout << "DoCastorAnalysis: Processing Event Number: " << eventIndex << std::endl;
   eventIndex++;
 }

References gather_cfg::cout, and eventIndex.

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

◆ update() [2/6]

void DoCastorAnalysis::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 92 of file DoCastorAnalysis.cc.

92 { std::cout << " Starting new job " << std::endl; }

References gather_cfg::cout.

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

◆ update() [3/6]

void DoCastorAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 96 of file DoCastorAnalysis.cc.

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

References gather_cfg::cout, and eventIndex.

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

◆ update() [4/6]

void DoCastorAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 113 of file DoCastorAnalysis.cc.

                                                    {
   // 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();
   if (debug)
     std::cout << "nentries in CAFI: " << nentries << std::endl;
   if (debug)
     getchar();
  
   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_time=&simhit_time;
   //psimhit_time->clear();
   //psimhit_time->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() / 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_time->push_back(time);
       psimhit_sector->push_back(sector);
       psimhit_module->push_back(zmodule);
  
       simhit_etot += energy;
  
       if (debug)
         std::cout << "hit " << ihit + 1 << " : x = " << (*psimhit_x)[ihit] << " , eta =  " << (*psimhit_eta)[ihit]
                   << " , phi = " << (*psimhit_phi)[ihit] << " , energy = " << (*psimhit_energy)[ihit] << std::endl;
     }
  
     //if(debug) std::cout<<" total energy = "<<simhit_etot<<std::endl;
     if (debug)
       getchar();
     CastorTree->Fill();
  
   }  // nentries > 0
   delete theCastorNumScheme;
 }

References CastorTree, gather_cfg::cout, debug, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getUnitID(), GeV, dqm-mbProfile::log, phi, 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, funct::tan(), theta(), CastorNumberingScheme::unpackIndex(), and ecaldqm::zside().

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

◆ update() [5/6]

void DoCastorAnalysis::update ( const EndOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun * >.

Definition at line 217 of file DoCastorAnalysis.cc.

217 { ; }

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

◆ 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 219 of file DoCastorAnalysis.cc.

219 { ; }

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

Member Data Documentation

◆ CastorOutputEventFile

TFile* DoCastorAnalysis::CastorOutputEventFile

private

Definition at line 87 of file DoCastorAnalysis.h.

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

◆ CastorTree

TTree* DoCastorAnalysis::CastorTree

private

Definition at line 88 of file DoCastorAnalysis.h.

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

◆ eventIndex

int DoCastorAnalysis::eventIndex

private

Definition at line 90 of file DoCastorAnalysis.h.

Referenced by update().

◆ psimhit_energy

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

private

Definition at line 98 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_eta

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

private

Definition at line 98 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_module

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

private

Definition at line 99 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_phi

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

private

Definition at line 98 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_sector

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

private

Definition at line 99 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_x

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

private

Definition at line 97 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_y

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

private

Definition at line 97 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ psimhit_z

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

private

Definition at line 97 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ simhit_energy

std::vector<double> DoCastorAnalysis::simhit_energy

private

Definition at line 93 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_eta

std::vector<double> DoCastorAnalysis::simhit_eta

private

Definition at line 93 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_etot

double DoCastorAnalysis::simhit_etot

private

Definition at line 102 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis(), and update().

◆ simhit_module

std::vector<int> DoCastorAnalysis::simhit_module

private

Definition at line 94 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_phi

std::vector<double> DoCastorAnalysis::simhit_phi

private

Definition at line 93 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_sector

std::vector<int> DoCastorAnalysis::simhit_sector

private

Definition at line 94 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_x

std::vector<double> DoCastorAnalysis::simhit_x

private

Definition at line 92 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_y

std::vector<double> DoCastorAnalysis::simhit_y

private

Definition at line 92 of file DoCastorAnalysis.h.

Referenced by update().

◆ simhit_z

std::vector<double> DoCastorAnalysis::simhit_z

private

Definition at line 92 of file DoCastorAnalysis.h.

Referenced by update().

◆ TreeFileName

std::string DoCastorAnalysis::TreeFileName

private

Definition at line 85 of file DoCastorAnalysis.h.

Referenced by DoCastorAnalysis().

◆ verbosity

int DoCastorAnalysis::verbosity

private

Definition at line 83 of file DoCastorAnalysis.h.

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