#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::ProducerBase &iProd)

	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 35 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 49 of file HcalForwardAnalysis.cc.

49 {

50 }

HcalForwardAnalysis::HcalForwardAnalysis ( const HcalForwardAnalysis & )

privatedelete

Member Function Documentation

void HcalForwardAnalysis::clear ( void )

private

Definition at line 251 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 230 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 68 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 245 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 56 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 129 of file HcalForwardAnalysis.cc.

References funct::cos(), gather_cfg::cout, FiberG4Hit::depth(), HFShowerG4Hit::globalPosition(), HFShowerG4Hit::hitId(), gen::k, HFShowerPhoton::lambda(), LogDebug, parseDetId(), phi, FiberG4Hit::photon(), nanoDQM_cfi::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 92 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 99 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 106 of file HcalForwardAnalysis.cc.

106 {

107 }

void HcalForwardAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 109 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;
 }