#include <HcalForwardAnalysis.h>

Inheritance diagram for HcalForwardAnalysis:

Classes
struct	Photon

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

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

	~HcalForwardAnalysis () override

Public Member Functions inherited from SimProducer
void	registerProducts (edm::ProducesCollector producesCollector)

	SimProducer ()

Public Member Functions inherited from SimWatcher
	SimWatcher ()

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

	HcalForwardAnalysis (const HcalForwardAnalysis &)=delete

void	init ()

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

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]

float	lambda [10000]

int	nphot

float	primMomX

float	primMomY

float	primMomZ

float	primT

float	primX

float	primY

float	primZ

float	t [10000]

int	theEventCounter

edm::Service< TFileService >	theFile

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)

Detailed Description

Definition at line 33 of file HcalForwardAnalysis.h.

Constructor & Destructor Documentation

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

Definition at line 27 of file HcalForwardAnalysis.cc.

References fiberId, edm::ParameterSet::getParameter(), mps_fire::i, init(), lambda, nphot, 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::~HcalForwardAnalysis ( )

override

Definition at line 47 of file HcalForwardAnalysis.cc.

47 {}

HcalForwardAnalysis::HcalForwardAnalysis ( const HcalForwardAnalysis & )

privatedelete

Member Function Documentation

void HcalForwardAnalysis::clear ( void )

private

Definition at line 246 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();
 }

void HcalForwardAnalysis::fillEvent ( )

private

Definition at line 225 of file HcalForwardAnalysis.cc.

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

Referenced by produce().

                                     {
   /*
      edm::LogInfo("HcalForwardLib") << "HcalForwardAnalysis: =====> filledEvent";
      */
   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();
 }

void HcalForwardAnalysis::init ( void )

private

Definition at line 64 of file HcalForwardAnalysis.cc.

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

Referenced by HcalForwardAnalysis().

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

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

privatedelete

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

private

Definition at line 240 of file HcalForwardAnalysis.cc.

Referenced by setPhotons().

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

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

overridevirtual

Implements SimProducer.

Definition at line 53 of file HcalForwardAnalysis.cc.

References fillEvent().

                                                                         {
   //std::auto_ptr<HFShowerPhotonCollection> product(new HFShowerPhotonCollection);
   //edm::LogInfo("HcalForwardLib") << "HcalForwardAnalysis: =====> Filling  event";
   //fillEvent(*product);
   //iEvent.put(product);
   //std::auto_ptr<PHcalForwardLibInfo> product(new PHcalForwardLibInfo);
   //fillEvent(*product);
   fillEvent();
   //iEvent.put(product);
 }

void HcalForwardAnalysis::setPhotons ( const EndOfEvent * evt )

private

Definition at line 118 of file HcalForwardAnalysis.cc.

References funct::cos(), gather_cfg::cout, FiberG4Hit::depth(), HFShowerG4Hit::globalPosition(), HFShowerG4Hit::hitId(), dqmiolumiharvest::j, dqmdumpme::k, HFShowerPhoton::lambda(), LogDebug, parseDetId(), phi, FiberG4Hit::photon(), nanoDQM_cff::Photon, HFShowerG4Hit::primaryMomDir(), primMomX, primMomY, primMomZ, primT, primX, primY, primZ, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, HFShowerPhoton::t(), theNames, thePhotons, theta(), HFShowerG4Hit::time(), FiberG4Hit::towerId(), HFShowerG4Hit::trackId(), HFShowerPhoton::x(), HFShowerPhoton::y(), and HFShowerPhoton::z().

Referenced by update().

                                                           {
   int idHC, j;
   FiberG4HitsCollection* theHC;
   // Look for the Hit Collection of HCal
   G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
   std::string sdName = theNames[0];  //name for fiber hits
   idHC = G4SDManager::GetSDMpointer()->GetCollectionID(sdName);
   theHC = (FiberG4HitsCollection*)allHC->GetHC(idHC);
   LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() Hit Collection for " << sdName << " of ID " << idHC
                              << " is obtained at " << theHC;
   std::vector<HFShowerPhoton> ShortFiberPhotons;
   std::vector<HFShowerPhoton> LongFiberPhotons;
   LongFiberPhotons.clear();
   ShortFiberPhotons.clear();
   if (idHC >= 0 && theHC != nullptr) {
     std::cout << "FiberhitSize " << theHC->entries() << std::endl;
     for (j = 0; j < theHC->entries(); j++) {
       FiberG4Hit* aHit = (*theHC)[j];
       std::vector<HFShowerPhoton> thePhotonsFromHit = aHit->photon();
       std::cout << "Fiberhit " << j << " has " << thePhotonsFromHit.size() << " photons." << std::endl;
       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]);
       }
       LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() NbPhotons " << thePhotonsFromHit.size()
                                  << " towerId " << fTowerId << " cellId " << fCellId << " fiberId " << fFiberId
                                  << " depth " << aHit->depth();
     }
   } else {
     LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons(): No Photons!";
     return;
   }
   LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() LongFibPhotons: " << LongFiberPhotons.size()
                              << " ShortFibPhotons: " << ShortFiberPhotons.size();
   std::cout << "HcalForwardAnalysis::setPhotons() LongFibPhotons: " << LongFiberPhotons.size()
             << " ShortFibPhotons: " << ShortFiberPhotons.size() << std::endl;
 
   //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) {
     LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() Chamber Hits size: " << theChamberHC->entries();
     for (j = 0; j < theChamberHC->entries(); ++j) {
       HFShowerG4Hit* aHit = (*theChamberHC)[j];
       LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons() Chamber Hit id " << aHit->hitId() << " track id "
                                  << aHit->trackId() << " prim. pos. " << aHit->globalPosition() << " prom mom. dir. "
                                  << aHit->primaryMomDir() << " time " << aHit->time();
       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 {
     LogDebug("HcalForwardLib") << "HcalForwardAnalysis::setPhotons(): No Chamber hits. Something is wrong!";
     return;
   }
   primX = primPosOnSurf.x();
   primY = primPosOnSurf.y();
   primZ = primPosOnSurf.z();
   primT = primTimeOnSurf;
   primMomX = primMomDirOnSurf.x();
   primMomY = primMomDirOnSurf.y();
   primMomZ = primMomDirOnSurf.z();
   //production point of Cherenkov photons in long fibers w.r.t local coordinates
   float photonProdX = 0.;
   float photonProdY = 0.;
   float photonProdZ = 0.;
   float photonProdTime = 0.;
   //angles for rotation matrices
   double theta = primMomDirOnSurf.theta();
   double phi = primMomDirOnSurf.phi();
 
   for (unsigned int k = 0; k < LongFiberPhotons.size(); ++k) {
     HFShowerPhoton aPhoton = LongFiberPhotons[k];
     photonProdX = aPhoton.x() * cos(theta) * cos(phi) + aPhoton.y() * cos(theta) * sin(phi) - aPhoton.z() * sin(theta) -
                   primPosOnSurf.x();
     photonProdY = -aPhoton.x() * sin(phi) + aPhoton.y() * cos(phi) - primPosOnSurf.y();
     photonProdZ = aPhoton.x() * sin(theta) * cos(phi) + aPhoton.y() * sin(theta) * sin(phi) + aPhoton.z() * cos(theta) -
                   primPosOnSurf.z();
     photonProdTime = aPhoton.t() - primTimeOnSurf;
     thePhotons.push_back(Photon(1, photonProdX, photonProdY, photonProdZ, photonProdTime, aPhoton.lambda()));
   }
   for (unsigned int k = 0; k < ShortFiberPhotons.size(); ++k) {
     HFShowerPhoton aPhoton = ShortFiberPhotons[k];
     photonProdX = aPhoton.x() * cos(theta) * cos(phi) + aPhoton.y() * cos(theta) * sin(phi) - aPhoton.z() * sin(theta) -
                   primPosOnSurf.x();
     photonProdY = -aPhoton.x() * sin(phi) + aPhoton.y() * cos(phi) - primPosOnSurf.y();
     photonProdZ = aPhoton.x() * sin(theta) * cos(phi) + aPhoton.y() * sin(theta) * sin(phi) + aPhoton.z() * cos(theta) -
                   primPosOnSurf.z();
     photonProdTime = aPhoton.t() - primTimeOnSurf;
     thePhotons.push_back(Photon(2, photonProdX, photonProdY, photonProdZ, photonProdTime, aPhoton.lambda()));
   }
 }

void HcalForwardAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 87 of file HcalForwardAnalysis.cc.

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

void HcalForwardAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 92 of file HcalForwardAnalysis.cc.

References clear(), and evNum.

                                                         {
   evNum = (*evt)()->GetEventID();
   clear();
   edm::LogInfo("HcalForwardLib") << "HcalForwardAnalysis: =====> Begin of event = " << evNum;
 }

void HcalForwardAnalysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 98 of file HcalForwardAnalysis.cc.

98 {}

void HcalForwardAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 100 of file HcalForwardAnalysis.cc.

References count, LogDebug, and setPhotons().

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