#include <HcalForwardAnalysis.h>

Inheritance diagram for HcalForwardAnalysis:

Classes
struct	Photon

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

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

virtual	~HcalForwardAnalysis ()

Public Member Functions inherited from SimProducer
void	registerProducts (edm::EDProducer &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 &)

void	init ()

const HcalForwardAnalysis &	operator= (const HcalForwardAnalysis &)

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

void	setPhotons (const EndOfEvent *evt)

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

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

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

void	update (const EndOfEvent *evt)
	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(), 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 ( )

virtual

Definition at line 49 of file HcalForwardAnalysis.cc.

49 {

50 }

HcalForwardAnalysis::HcalForwardAnalysis ( const HcalForwardAnalysis & )

private

Member Function Documentation

void HcalForwardAnalysis::clear ( void )

private

Definition at line 251 of file HcalForwardAnalysis.cc.

References fiberId, 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, i, 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, TFileDirectory::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 & )

private

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

virtual

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(), j, gen::k, HFShowerPhoton::lambda(), LogDebug, parseDetId(), phi, FiberG4Hit::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> 0) {
         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> 0) {
         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 * )

privatevirtual

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

privatevirtual

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

privatevirtual

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

privatevirtual

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;
 }