Inheritance diagram for HcalForwardAnalysis:

Classes
struct	Photon

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

	HcalForwardAnalysis (const HcalForwardAnalysis &)=delete

const HcalForwardAnalysis &	operator= (const HcalForwardAnalysis &)=delete

void	produce (edm::Event &, const edm::EventSetup &) override

	~HcalForwardAnalysis () override

Public Member Functions inherited from SimProducer
const SimProducer &	operator= (const SimProducer &)=delete

void	registerProducts (edm::ProducesCollector producesCollector)

	SimProducer ()

	SimProducer (const SimProducer &)=delete

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 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 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	clear ()

void	fillEvent ()

void	init ()

void	parseDetId (int id, int &tower, int &cell, int &fiber)

void	setPhotons (const EndOfEvent *evt)

void	update (const BeginOfRun *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 G4Step *step) 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...

Private Attributes
int	count

int	evNum

int	fiberId [10000]

bool	fillt

float	lambda [10000]

int	nphot

float	primMomX

float	primMomY

float	primMomZ

float	primT

float	primX

float	primY

float	primZ

float	t [10000]

int	theEventCounter

std::vector< std::string >	theNames

std::vector< Photon >	thePhotons

TTree *	theTree

float	x [10000]

float	y [10000]

float	z [10000]

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()

template<class T >
void	produces (const std::string &instanceName)

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

Detailed Description

Definition at line 39 of file HcalForwardAnalysis.cc.

Constructor & Destructor Documentation

◆ HcalForwardAnalysis() [1/2]

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

Definition at line 94 of file HcalForwardAnalysis.cc.

References fiberId, edm::ParameterSet::getParameter(), mps_fire::i, init(), lambda, nphot, AlCaHLTBitMon_ParallelJobs::p, primMomX, primMomY, primMomZ, primT, primX, primY, primZ, t, theEventCounter, theNames, x, y, and z.

                                                                  {
   edm::ParameterSet m_SLP = p.getParameter<edm::ParameterSet>("HFShowerLibraryProducer");
   theNames = m_SLP.getParameter<std::vector<std::string> >("Names");
   //LibVer = m_HS.getParameter<std::string> ("LibVer");
   //produces<HFShowerPhotonCollection> ();
   init();
   theEventCounter = 0;
   nphot = 0;
   for (int i = 0; i < 10000; ++i) {
     x[i] = 0.;
     y[i] = 0.;
     z[i] = 0.;
     t[i] = 0.;
     lambda[i] = 0.;
     fiberId[i] = 0;
   }
   primX = primY = primZ = primT = 0.;
   primMomX = primMomY = primMomZ = 0.;
 }

◆ HcalForwardAnalysis() [2/2]

HcalForwardAnalysis::HcalForwardAnalysis ( const HcalForwardAnalysis & )

delete

◆ ~HcalForwardAnalysis()

HcalForwardAnalysis::~HcalForwardAnalysis ( )

override

Definition at line 114 of file HcalForwardAnalysis.cc.

114 {}

Member Function Documentation

◆ clear()

void HcalForwardAnalysis::clear ( void )

private

Definition at line 362 of file HcalForwardAnalysis.cc.

References fiberId, mps_fire::i, lambda, nphot, primMomX, primMomY, primMomZ, primT, primX, primY, primZ, t, thePhotons, x, y, and z.

Referenced by init(), and update().

                                 {
   nphot = 0;
   for (int i = 0; i < 10000; ++i) {
     x[i] = 0.;
     y[i] = 0.;
     z[i] = 0.;
     t[i] = 0.;
     lambda[i] = 0.;
     fiberId[i] = 0;
   }
   primX = primY = primZ = primT = 0.;
   primMomX = primMomY = primMomZ = 0.;
 
   thePhotons.clear();
 }

◆ fillEvent()

void HcalForwardAnalysis::fillEvent ( )

private

Definition at line 340 of file HcalForwardAnalysis.cc.

References fiberId, mps_fire::i, createfilelist::int, lambda, nphot, t, thePhotons, theTree, x, y, and z.

Referenced by produce().

                                     {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis: =====> filledEvent";
 #endif
   nphot = int(thePhotons.size());
   for (int i = 0; i < nphot; ++i) {
     x[i] = thePhotons[i].x;
     y[i] = thePhotons[i].y;
     z[i] = thePhotons[i].z;
     t[i] = thePhotons[i].t;
     lambda[i] = thePhotons[i].lambda;
     fiberId[i] = thePhotons[i].fiberId;
   }
   theTree->Fill();
 }

◆ init()

void HcalForwardAnalysis::init ( void )

private

Definition at line 125 of file HcalForwardAnalysis.cc.

References clear(), count, evNum, fiberId, lambda, nphot, primMomX, primMomY, primMomZ, primT, primX, primY, primZ, t, interactiveExample::theFile, theTree, x, y, and z.

Referenced by HcalForwardAnalysis().

                                {
   edm::Service<TFileService> theFile;
   theTree = theFile->make<TTree>("CherenkovPhotons", "Cherenkov Photons");
   theTree->Branch("nphot", &nphot, "nphot/I");
   theTree->Branch("x", &x, "x[nphot]/F");
   theTree->Branch("y", &y, "y[nphot]/F");
   theTree->Branch("z", &z, "z[nphot]/F");
   theTree->Branch("t", &t, "t[nphot]/F");
   theTree->Branch("lambda", &lambda, "lambda[nphot]/F");
   theTree->Branch("fiberId", &fiberId, "fiberId[nphot]/I");
   theTree->Branch("primX", &primX, "primX/F");
   theTree->Branch("primY", &primY, "primY/F");
   theTree->Branch("primZ", &primZ, "primZ/F");
   theTree->Branch("primMomX", &primMomX, "primMomX/F");
   theTree->Branch("primMomY", &primMomY, "primMomY/F");
   theTree->Branch("primMomZ", &primMomZ, "primMomZ/F");
   theTree->Branch("primT", &primT, "primT/F");
 
   // counter
   count = 0;
   evNum = 0;
   clear();
 }

◆ operator=()

const HcalForwardAnalysis& HcalForwardAnalysis::operator= ( const HcalForwardAnalysis & )

delete

◆ parseDetId()

void HcalForwardAnalysis::parseDetId	(	int	id,
		int &	tower,
		int &	cell,
		int &	fiber
	)

private

Definition at line 356 of file HcalForwardAnalysis.cc.

References hgcalTowerProducer_cfi::tower.

Referenced by setPhotons().

                                                                               {
   tower = id / 10000;
   cell = id / 10 - tower * 10;
   fiber = id - tower * 10000 - cell * 10;
 }

◆ produce()

void HcalForwardAnalysis::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &
	)

overridevirtual

Implements SimProducer.

Definition at line 120 of file HcalForwardAnalysis.cc.

References fillEvent(), and fillt.

                                                                         {
   if (fillt)
     fillEvent();
 }

◆ setPhotons()

void HcalForwardAnalysis::setPhotons ( const EndOfEvent * evt )

private

Definition at line 184 of file HcalForwardAnalysis.cc.

References funct::cos(), FiberG4Hit::depth(), fillt, HFShowerG4Hit::globalPosition(), HFShowerG4Hit::hitId(), dqmiolumiharvest::j, dqmdumpme::k, HFShowerPhoton::lambda(), Skims_PA_cff::name, parseDetId(), phi, FiberG4Hit::photon(), nanoDQM_cff::Photon, HFShowerG4Hit::primaryMomDir(), primMomX, primMomY, primMomZ, primT, primX, primY, primZ, heavyFlavorDQMFirstStep_cff::sdName, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, HFShowerPhoton::t(), theNames, thePhotons, theta(), HFShowerG4Hit::time(), FiberG4Hit::towerId(), HFShowerG4Hit::trackId(), HFShowerPhoton::x(), geometryCSVtoXML::xx, HFShowerPhoton::y(), geometryCSVtoXML::yy, HFShowerPhoton::z(), gpuVertexFinder::zv, and geometryCSVtoXML::zz.

Referenced by update().

                                                           {
   fillt = true;
   int idHC, j;
   FiberG4HitsCollection* theHC;
   // Look for the Hit Collection of HCal
   G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis:: Has " << allHC->GetNumberOfCollections()
                                      << " collections";
   for (int k = 0; k < allHC->GetNumberOfCollections(); ++k) {
     G4String name = (allHC->GetHC(k) == nullptr) ? "Unknown" : allHC->GetHC(k)->GetName();
     G4String nameSD = (allHC->GetHC(k) == nullptr) ? "Unknown" : allHC->GetHC(k)->GetSDname();
     edm::LogVerbatim("HcalForwardLib") << "Collecttion[" << k << "] " << allHC->GetHC(k) << "  " << name << ":"
                                        << nameSD;
   }
 #endif
   std::string sdName = theNames[0];  //name for fiber hits
   idHC = G4SDManager::GetSDMpointer()->GetCollectionID(sdName);
   theHC = (FiberG4HitsCollection*)allHC->GetHC(idHC);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() Hit Collection for " << sdName << " of ID "
                                      << idHC << " is obtained at " << theHC;
 #endif
   std::vector<HFShowerPhoton> ShortFiberPhotons;
   std::vector<HFShowerPhoton> LongFiberPhotons;
   LongFiberPhotons.clear();
   ShortFiberPhotons.clear();
   if (idHC >= 0 && theHC != nullptr) {
     int thehc_entries = theHC->entries();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalForwardLib") << "FiberhitSize " << thehc_entries;
 #endif
     for (j = 0; j < thehc_entries; j++) {
       FiberG4Hit* aHit = (*theHC)[j];
       std::vector<HFShowerPhoton> thePhotonsFromHit = aHit->photon();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalForwardLib") << "Fiberhit " << j << " has " << thePhotonsFromHit.size() << " photons.";
 #endif
       int fTowerId = -1;
       int fCellId = -1;
       int fFiberId = -1;
       parseDetId(aHit->towerId(), fTowerId, fCellId, fFiberId);
       for (unsigned int iph = 0; iph < thePhotonsFromHit.size(); ++iph) {
         if (aHit->depth() == 1)
           LongFiberPhotons.push_back(thePhotonsFromHit[iph]);
         if (aHit->depth() == 2)
           ShortFiberPhotons.push_back(thePhotonsFromHit[iph]);
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() NbPhotons " << thePhotonsFromHit.size()
                                          << " towerId " << fTowerId << " cellId " << fCellId << " fiberId " << fFiberId
                                          << " depth " << aHit->depth();
 #endif
     }
   } else {
     fillt = false;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons(): No Photons!";
 #endif
     return;
   }
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() LongFibPhotons: " << LongFiberPhotons.size()
                                      << " ShortFibPhotons: " << ShortFiberPhotons.size();
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() LongFibPhotons: " << LongFiberPhotons.size()
                                      << " ShortFibPhotons: " << ShortFiberPhotons.size();
 
   //Chamber hits to find information about primary particle on surface
   HFShowerG4HitsCollection* theChamberHC;
   G4HCofThisEvent* allChamberHC = (*evt)()->GetHCofThisEvent();
   sdName = theNames[1];
   idHC = G4SDManager::GetSDMpointer()->GetCollectionID(sdName);
   theChamberHC = (HFShowerG4HitsCollection*)allChamberHC->GetHC(idHC);
   math::XYZPoint primPosOnSurf(0, 0, 0);
   math::XYZPoint primMomDirOnSurf(0, 0, 0);
   float primTimeOnSurf = 0;
   //    the chamber hit is for primary particle, but step size can be small
   //    (in newer Geant4 versions) and as a result primary particle may have
   //    multiple hits. We want to take last one which is close the HF absorber
   //  if (idHC >= 0 && theChamberHC != nullptr && theChamberHC->entries()>0) {
   if (idHC >= 0 && theChamberHC != nullptr) {
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() Chamber Hits size: "
                                        << theChamberHC->entries();
     int thec_hc_entries = theChamberHC->entries();
     for (j = 0; j < thec_hc_entries; ++j) {
       HFShowerG4Hit* aHit = (*theChamberHC)[j];
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() Chamber Hit id " << aHit->hitId()
                                          << " track id " << aHit->trackId() << " prim. pos. " << aHit->globalPosition()
                                          << " prom mom. dir. " << aHit->primaryMomDir() << " time " << aHit->time();
 #endif
       primPosOnSurf.SetXYZ(aHit->globalPosition().x(), aHit->globalPosition().y(), aHit->globalPosition().z());
       primMomDirOnSurf.SetXYZ(aHit->primaryMomDir().x(), aHit->primaryMomDir().y(), aHit->primaryMomDir().z());
       primTimeOnSurf = aHit->time();
     }
   } else {
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons(): No Chamber hits are stored";
     fillt = false;
     return;
   }
   primX = primPosOnSurf.x();
   primY = primPosOnSurf.y();
   primZ = primPosOnSurf.z();
   if (primZ < 990) {  // there were interactions before HF
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::setPhotons(): First interaction before HF";
     fillt = false;
     return;
   }
   primT = primTimeOnSurf;
   primMomX = primMomDirOnSurf.x();
   primMomY = primMomDirOnSurf.y();
   primMomZ = primMomDirOnSurf.z();
   //angles for rotation matrices
   double theta = primMomDirOnSurf.theta();
   double phi = primMomDirOnSurf.phi();
 
   // my insert ----------------------------------------------------------------
   double sphi = sin(phi);
   double cphi = cos(phi);
   double ctheta = cos(theta);
   double stheta = sin(theta);
 
   double pex = 0, pey = 0, zv = 0;
   double xx, yy, zz;
 
   for (unsigned int k = 0; k < LongFiberPhotons.size(); ++k) {
     HFShowerPhoton aPhoton = LongFiberPhotons[k];
     // global coordinates
     xx = aPhoton.x();
     yy = aPhoton.y();
     zz = aPhoton.z();
 
     // local coordinates in rotated to shower axis system and vs shower origin
     pex = xx * ctheta * cphi + yy * ctheta * sphi - zz * stheta;
     pey = -xx * sphi + yy * cphi;
     zv = xx * stheta * cphi + yy * stheta * sphi + zz * ctheta - primZ / ctheta;
 
     double photonProdTime = aPhoton.t() - primTimeOnSurf;
     thePhotons.push_back(Photon(1, pex, pey, zv, photonProdTime, aPhoton.lambda()));
   }
   for (unsigned int k = 0; k < ShortFiberPhotons.size(); ++k) {
     HFShowerPhoton aPhoton = ShortFiberPhotons[k];
     // global coordinates
     xx = aPhoton.x();
     yy = aPhoton.y();
     zz = aPhoton.z();
 
     // local coordinates in rotated to shower axis system and vs shower origin
     pex = xx * ctheta * cphi + yy * ctheta * sphi - zz * stheta;
     pey = -xx * sphi + yy * cphi;
     zv = xx * stheta * cphi + yy * stheta * sphi + zz * ctheta - primZ / ctheta;
 
     double photonProdTime = aPhoton.t() - primTimeOnSurf;
     thePhotons.push_back(Photon(2, pex, pey, zv, photonProdTime, aPhoton.lambda()));
   }
 }

◆ update() [1/4]

void HcalForwardAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 149 of file HcalForwardAnalysis.cc.

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

                                                       {
   int irun = (*run)()->GetRunID();
   edm::LogVerbatim("HcalForwardLib") << " =====> Begin of Run = " << irun;
 }

◆ update() [2/4]

void HcalForwardAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 154 of file HcalForwardAnalysis.cc.

References clear(), and evNum.

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

                                                         {
   evNum = (*evt)()->GetEventID();
   clear();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis: =====> Begin of event = " << evNum;
 #endif
 }

◆ update() [3/4]

void HcalForwardAnalysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 162 of file HcalForwardAnalysis.cc.

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

162 {}

◆ update() [4/4]

void HcalForwardAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 164 of file HcalForwardAnalysis.cc.

References count, and setPhotons().

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

                                                       {
   count++;
 
   //fill the buffer
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis::Fill event " << (*evt)()->GetEventID();
 #endif
   setPhotons(evt);
 
   int iEvt = (*evt)()->GetEventID();
   if (iEvt < 10)
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis:: Event " << iEvt;
   else if ((iEvt < 100) && (iEvt % 10 == 0))
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis:: Event " << iEvt;
   else if ((iEvt < 1000) && (iEvt % 100 == 0))
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis:: Event " << iEvt;
   else if ((iEvt < 10000) && (iEvt % 1000 == 0))
     edm::LogVerbatim("HcalForwardLib") << "HcalForwardAnalysis:: Event " << iEvt;
 }