#include <HcalTestAnalysis.h>

Inheritance diagram for HcalTestAnalysis:

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

	~HcalTestAnalysis () 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 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	fill (const EndOfEvent *ev)

void	layerAnalysis ()

std::vector< int >	layerGrouping (int)

void	qieAnalysis (CLHEP::HepRandomEngine *)

double	timeOfFlight (int det, int layer, double eta)

std::vector< int >	towersToAdd (int centre, int nadd)

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

Private Attributes
int	addTower

std::vector< CaloHit >	caloHitCache

int	centralTower

unsigned int	count

double	edepEB

double	edepEE

double	edepHB

double	edepHE

double	edepHO

double	edepl [20]

double	eta0

std::string	fileName

std::vector< int >	group_

HcalDDDSimConstants *	hcons

double	mudist [20]

HcalQie *	myqie

std::vector< std::string >	names

int	nGroup

int	nTower

HcalNumberingFromDDD *	numberingFromDDD

HcalTestNumberingScheme *	org

double	phi0

std::vector< int >	tower_

HcalTestHistoClass *	tuples

std::unique_ptr< HcalTestHistoManager >	tuplesManager

Additional Inherited Members

Detailed Description

Definition at line 35 of file HcalTestAnalysis.h.

Constructor & Destructor Documentation

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

Definition at line 32 of file HcalTestAnalysis.cc.

References addTower, centralTower, count, eta0, fileName, edm::ParameterSet::getParameter(), group_, mps_fire::i, layerGrouping(), myqie, names, nGroup, nTower, phi0, AlCaHLTBitMon_QueryRunRegistry::string, tower_, towersToAdd(), and tuplesManager.

                                                           : 
   addTower(3),tuples(nullptr),numberingFromDDD(nullptr),hcons(nullptr),org(nullptr) {
 
   edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("HcalTestAnalysis");
   eta0         = m_Anal.getParameter<double>("Eta0");
   phi0         = m_Anal.getParameter<double>("Phi0");
   int laygroup = m_Anal.getParameter<int>("LayerGrouping");
   centralTower = m_Anal.getParameter<int>("CentralTower");
   names        = m_Anal.getParameter<std::vector<std::string> >("Names");
   fileName     = m_Anal.getParameter<std::string>("FileName");
 
   tuplesManager.reset(nullptr); 
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Initialised as observer of "
               << "begin/end events and of G4step";
 
   count  = 0;
   group_ = layerGrouping(laygroup);
   nGroup = 0;
   for (unsigned int i=0; i<group_.size(); i++) 
     if (group_[i]>nGroup) nGroup = group_[i];
   tower_ = towersToAdd(centralTower, addTower);
   nTower = tower_.size()/2;
 
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: initialised for " << nGroup 
               << " Longitudinal groups and " << nTower 
               << " towers";
 
   // qie
   myqie  = new HcalQie(p);
 } 

HcalTestAnalysis::~HcalTestAnalysis ( )

override

Definition at line 63 of file HcalTestAnalysis.cc.

References count, myqie, numberingFromDDD, org, and tuples.

                                     {
   edm::LogInfo("HcalSim") << "HcalTestAnalysis: -------->  Total number of "
               << "selected entries : " << count;
   edm::LogInfo("HcalSim") << "HcalTestAnalysis: Pointers:: HcalQie " << myqie 
               << ", HistoClass " << tuples << ", Numbering Scheme "
               << org << " and FromDDD " << numberingFromDDD;
   if (myqie)  {
     edm::LogInfo("HcalSim") << "HcalTestAnalysis: Delete HcalQie";
     delete myqie;
   }
   if (numberingFromDDD) {
     edm::LogInfo("HcalSim") << "HcalTestAnalysis: Delete HcalNumberingFromDDD";
     delete numberingFromDDD;
   }
 }

Member Function Documentation

void HcalTestAnalysis::fill ( const EndOfEvent * ev )

private

Definition at line 344 of file HcalTestAnalysis.cc.

References caloHitCache, MillePedeFileConverter_cfg::e, PVValHelper::eta, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getTimeSlice(), HcalTestNumberingScheme::getUnitID(), CaloG4Hit::getUnitID(), GeV, HcalEndcap, HcalForward, cmsBatch::log, LogDebug, names, numberingFromDDD, org, HcalTestNumberingScheme::packHcalIndex(), phi, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), ntuplemaker::time, timeOfFlight(), HcalNumberingFromDDD::unitID(), HcalTestNumberingScheme::unpackHcalIndex(), and HcalDigiParam_cfi::zside.

Referenced by update().

                                                   {
 
   LogDebug("HcalSim") << "HcalTestAnalysis: Fill event " 
               << (*evt)()->GetEventID();
   
   // access to the G4 hit collections 
   G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
   
   int nhc = 0, neb = 0, nef = 0, j = 0;    
   caloHitCache.erase (caloHitCache.begin(), caloHitCache.end());
  
   // Hcal
   int HCHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[0]);
   CaloG4HitCollection* theHCHC = (CaloG4HitCollection*) allHC->GetHC(HCHCid);
   LogDebug("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names[0] 
               << " of ID " << HCHCid << " is obtained at " << theHCHC;
   if (HCHCid >= 0 && theHCHC != nullptr) {
     for (j = 0; j < theHCHC->entries(); j++) {
 
       CaloG4Hit* aHit = (*theHCHC)[j]; 
     
       double e        = aHit->getEnergyDeposit()/GeV;
       double time     = aHit->getTimeSlice();
       
       math::XYZPoint pos   = aHit->getPosition();
       double theta    = pos.theta();
       double   eta    = -log(tan(theta * 0.5));
       double   phi    = pos.phi();
     
       uint32_t unitID = aHit->getUnitID();
       int subdet, zside, layer, etaIndex, phiIndex, lay;
       org->unpackHcalIndex(unitID,subdet,zside,layer,etaIndex,phiIndex,lay);
       double jitter   = time-timeOfFlight(subdet,lay,eta);
       if (jitter<0) jitter = 0;
       CaloHit hit(subdet,lay,e,eta,phi,jitter,unitID);
       caloHitCache.push_back(hit);
       nhc++;
 
       std::string det =  "HB";
       if (subdet == static_cast<int>(HcalForward)) {
     det = "HF";
       } else if (subdet == static_cast<int>(HcalEndcap)) {
     if (etaIndex <= 20) {
       det  = "HES";
     } else {
       det = "HED";
     }
       }
 
       LogDebug("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) 
               << layer  << " time " << std::setw(6) << time 
               << " theta "  << std::setw(8) << theta << " eta " 
               << std::setw(8) << eta  << " phi " << std::setw(8) 
               << phi << " e " << std::setw(8) << e;
     }
   }
   LogDebug("HcalSim") << "HcalTestAnalysis::HCAL hits : " << nhc << "\n";
 
   // EB
   int EBHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[1]);
   CaloG4HitCollection* theEBHC = (CaloG4HitCollection*) allHC->GetHC(EBHCid);
   LogDebug("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names[1]
               << " of ID " << EBHCid << " is obtained at " << theEBHC;
   if (EBHCid >= 0 && theEBHC != nullptr) {
     for (j = 0; j < theEBHC->entries(); j++) {
 
       CaloG4Hit* aHit = (*theEBHC)[j]; 
     
       double e        = aHit->getEnergyDeposit()/GeV;
       double time     = aHit->getTimeSlice();
       std::string det =  "EB";
       
       math::XYZPoint pos  = aHit->getPosition();
       double theta    = pos.theta();
       double   eta    = -log(tan(theta/2.));
       double   phi    = pos.phi();
 
       HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta,phi,1,1);
       uint32_t unitID = org->getUnitID(id);
       int subdet, zside, layer, ieta, iphi, lay;
       org->unpackHcalIndex(unitID,subdet,zside,layer,ieta,iphi,lay);
       subdet = 10;
       layer  = 0;
       unitID = org->packHcalIndex(subdet,zside,layer,ieta,iphi,lay);
       CaloHit hit(subdet,lay,e,eta,phi,time,unitID);
       caloHitCache.push_back(hit);
       neb++;
       LogDebug("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) 
               << layer << " time " << std::setw(6) << time 
               << " theta " << std::setw(8) << theta << " eta " 
               << std::setw(8) << eta << " phi " << std::setw(8) 
               << phi << " e " << std::setw(8) << e;
     }
   }
   LogDebug("HcalSim") << "HcalTestAnalysis::EB hits : " << neb << "\n";
 
   // EE
   int EEHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[2]);
   CaloG4HitCollection* theEEHC = (CaloG4HitCollection*) allHC->GetHC(EEHCid);
   LogDebug("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names[2]
               << " of ID " << EEHCid << " is obtained at " << theEEHC;
   if (EEHCid >= 0 && theEEHC != nullptr) {
     for (j = 0; j < theEEHC->entries(); j++) {
 
       CaloG4Hit* aHit = (*theEEHC)[j]; 
     
       double e        = aHit->getEnergyDeposit()/GeV;
       double time     = aHit->getTimeSlice();
       std::string det = "EE";
       
       math::XYZPoint pos  = aHit->getPosition();
       double theta    = pos.theta();
       double   eta    = -log(tan(theta/2.));
       double   phi    = pos.phi();
 
       HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta,phi,1,1);
       uint32_t unitID = org->getUnitID(id);
       int subdet, zside, layer, ieta, iphi, lay;
       org->unpackHcalIndex(unitID,subdet,zside,layer,ieta,iphi,lay);
       subdet = 11;
       layer  = 0;
       unitID = org->packHcalIndex(subdet,zside,layer,ieta,iphi,lay);
       CaloHit hit(subdet,lay,e,eta,phi,time,unitID);
       caloHitCache.push_back(hit);
       nef++;
       LogDebug("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2)
               << layer << " time " << std::setw(6) << time 
               << " theta " << std::setw(8) << theta << " eta " 
               << std::setw(8) << eta  << " phi " << std::setw(8) 
               << phi << " e " << std::setw(8) << e;
     }
   }
   LogDebug("HcalSim") << "HcalTestAnalysis::EE hits : " << nef << "\n";
 }

void HcalTestAnalysis::layerAnalysis ( )

private

Definition at line 600 of file HcalTestAnalysis.cc.

References edepEB, edepEE, edepHB, edepHE, edepHO, edepl, HcalTestHistoClass::fillLayers(), mps_fire::i, LogDebug, MeV, mudist, and tuples.

Referenced by update().

                                     {
 
   int i = 0;
   LogDebug("HcalSim") << "\n ===>>> HcalTestAnalysis: Energy deposit in MeV " 
               << "\n at EB : " << std::setw(6) << edepEB/MeV 
               << "\n at EE : " << std::setw(6) << edepEE/MeV 
               << "\n at HB : " << std::setw(6) << edepHB/MeV
               << "\n at HE : " << std::setw(6) << edepHE/MeV
               << "\n at HO : " << std::setw(6) << edepHO/MeV
               << "\n ---- HcalTestAnalysis: Energy deposit in Layers"; 
   for (i = 0; i < 20; i++) 
     LogDebug("HcalSim") << " Layer " << std::setw(2) << i << " E " 
             << std::setw(8) << edepl[i]/MeV  << " MeV";
 
   tuples->fillLayers(edepl, edepHO, edepHB+edepHE, mudist);  
 }

std::vector< int > HcalTestAnalysis::layerGrouping ( int group )

private

Definition at line 79 of file HcalTestAnalysis.cc.

References mps_fire::i, and groupFilesInBlocks::temp.

Referenced by HcalTestAnalysis().

                                                         {
 
   std::vector<int> temp(19);
   if (group <= 1) {
     int grp[19] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2};
     for (int i=0; i<19; i++)
       temp[i] = grp[i];
   } else  if (group == 2) {
     int grp[19] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19};
     for (int i=0; i<19; i++)
       temp[i] = grp[i];
   } else if (group == 3) {
     int grp[19] = {1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7};
     for (int i=0; i<19; i++)
       temp[i] = grp[i];
   } else if (group == 4) {
     int grp[19] = {1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 7, 7};
     for (int i=0; i<19; i++)
       temp[i] = grp[i];
   } else {
     int grp[19] = {1, 1, 2, 3, 3, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 7, 7};
     for (int i=0; i<19; i++)
       temp[i] = grp[i];
   }
 
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Layer Grouping ";
   for (int i=0; i<19; i++)
     edm::LogInfo("HcalSim") << "HcalTestAnalysis: Group[" << i << "] = "
                   << temp[i];
   return temp;
 }

void HcalTestAnalysis::qieAnalysis ( CLHEP::HepRandomEngine * engine )

private

Definition at line 480 of file HcalTestAnalysis.cc.

References addTower, caloHitCache, centralTower, CaloHit::det(), CaloHit::e(), eta0, HcalTestHistoClass::fillHits(), HcalTestHistoClass::fillQie(), HcalQie::getCode(), HcalQie::getEnergy(), HcalTestNumberingScheme::getUnitID(), edmIntegrityCheck::group, group_, HcalBarrel, HcalEndcap, hfClusterShapes_cfi::hits, mps_fire::i, CaloHit::id(), training_settings::idx, CaloHit::layer(), LogDebug, myqie, nGroup, nt, nTower, numberingFromDDD, org, phi0, AlCaHLTBitMon_QueryRunRegistry::string, tower_, tuples, HcalNumberingFromDDD::unitID(), HcalTestNumberingScheme::unpackHcalIndex(), and HcalDigiParam_cfi::zside.

Referenced by update().

                                                                {
 
   //Fill tuple with hit information
   int hittot = caloHitCache.size();
   tuples->fillHits(caloHitCache);
 
   //Get the index of the central tower
   HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta0,phi0,1,1);
   uint32_t unitID = org->getUnitID(id);
   int subdet, zside, layer, ieta, iphi, lay;
   org->unpackHcalIndex(unitID,subdet,zside,layer,ieta,iphi,lay);
   int      laymax = 0;
   std::string det = "Unknown";
   if (subdet == static_cast<int>(HcalBarrel)) {
     laymax = 4; det = "HB";
   } else if (subdet == static_cast<int>(HcalEndcap)) {
     laymax = 2; det = "HES";
   }
   LogDebug("HcalSim") << "HcalTestAnalysis::Qie: " << det << " Eta " << ieta 
               << " Phi " << iphi << " Laymax " << laymax << " Hits "
               << hittot;
 
   if (laymax>0 && hittot>0) {
     std::vector<CaloHit>  hits(hittot);
     std::vector<double> eqielay(80,0.0),  esimlay(80,0.0),  esimtot(4,0.0);
     std::vector<double> eqietow(200,0.0), esimtow(200,0.0), eqietot(4,0.0);
     int etac = (centralTower/100)%100;
     int phic = (centralTower%100);
 
     for (int layr=0; layr<nGroup; layr++) {
       /*
       int layx, layy=20;
       for (int i=0; i<20; i++) 
     if (group_[i] == layr+1 && i < layy) layy = i+1;
       if (subdet == static_cast<int>(HcalBarrel)) {
     if      (layy < 2)    layx = 0;
     else if (layy < 17)   layx = 1;
     else if (layy == 17)  layx = 2;
     else                  layx = 3;
       } else {
     if      (layy < 2)    layx = 0;
     else                  layx = 1;
       }
       */
       for (int it=0; it<nTower; it++) {
     int    nhit = 0;
     double esim = 0;
     for (int k1 = 0; k1 < hittot; k1++) {
       CaloHit hit = caloHitCache[k1];
       int     subdetc = hit.det();
       int     layer   = hit.layer();
       int     group   = 0;
       if (layer > 0 && layer < 20) group = group_[layer];
       if (subdetc == subdet && group == layr+1) {
         int zsidec, ietac, iphic, idx;
         unitID = hit.id();
         org->unpackHcalIndex(unitID,subdetc,zsidec,layer,ietac,iphic,lay);
         if (etac > 0 && phic > 0) {
           idx = ietac*100 + iphic;
         } else if (etac > 0) {
           idx = ietac*100;
         } else if (phic > 0) {
           idx = iphic;
         } else {
           idx = 0;
         }
         if (zsidec==zside && idx==tower_[it]) {
           hits[nhit] = hit;
           LogDebug("HcalSim") << "HcalTest: Hit " << nhit << " " << hit;
           nhit++;
           esim += hit.e();
         }
       }
     }
 
     std::vector<int> cd = myqie->getCode(nhit, hits, engine);
     double         eqie = myqie->getEnergy(cd);
 
     LogDebug("HcalSim") << "HcalTestAnalysis::Qie: Energy in layer " 
                 << layr << " Sim " << esim << " After QIE " <<eqie;
     for (int i=0; i<4; i++) {
       if (tower_[nTower+it] <= i) {
         esimtot[i]         += esim;
         eqietot[i]         += eqie;
         esimlay[20*i+layr] += esim;
         eqielay[20*i+layr] += eqie;
         esimtow[50*i+it]   += esim;
         eqietow[50*i+it]   += eqie;
       }
     }
       }
     }
     LogDebug("HcalSim") << "HcalTestAnalysis::Qie: Total energy " << esimtot[3]
             << " (SimHit) " << eqietot[3] << " (After QIE)";
 
     std::vector<double> latphi(10);
     int nt = 2*addTower + 1;
     for (int it=0; it<nt; it++)
       latphi[it] = it-addTower;
     for (int i=0; i<4; i++) {
       double scals=1, scalq=1;
       std::vector<double> latfs(10,0.), latfq(10,0.), longs(20), longq(20);
       if (esimtot[i]>0) scals = 1./esimtot[i];
       if (eqietot[i]>0) scalq = 1./eqietot[i];
       for (int it=0; it<nTower; it++) {
     int phib = it%nt; 
     latfs[phib] += scals*esimtow[50*i+it];
     latfq[phib] += scalq*eqietow[50*i+it];
       }
       for (int layr=0; layr<=nGroup; layr++) {
     longs[layr] = scals*esimlay[20*i+layr];
     longq[layr] = scalq*eqielay[20*i+layr];
       }
       tuples->fillQie(i,esimtot[i],eqietot[i],nGroup,longs,longq,
               nt,latphi,latfs,latfq);
     }
   }
 }

double HcalTestAnalysis::timeOfFlight	(	int	det,
		int	layer,
		double	eta
	)

private

Definition at line 619 of file HcalTestAnalysis.cc.

References funct::cos(), JetChargeProducer_cfi::exp, HcalBarrel, LogDebug, funct::sin(), theta(), and tmp.

Referenced by fill().

                                                                     {
 
   double theta = 2.0*atan(exp(-eta));
   double dist  = 0.;
   if (det == static_cast<int>(HcalBarrel)) {
     const double rLay[19] = {
       1836.0, 1902.0, 1962.0, 2022.0, 2082.0, 2142.0, 2202.0, 2262.0, 2322.0, 
       2382.0, 2448.0, 2514.0, 2580.0, 2646.0, 2712.0, 2776.0, 2862.5, 3847.0,
       4052.0};
     if (layer>0 && layer<20) dist += rLay[layer-1]*mm/sin(theta);
   } else {
     const double zLay[19] = {
       4034.0, 4032.0, 4123.0, 4210.0, 4297.0, 4384.0, 4471.0, 4558.0, 4645.0, 
       4732.0, 4819.0, 4906.0, 4993.0, 5080.0, 5167.0, 5254.0, 5341.0, 5428.0,
       5515.0};
     if (layer>0 && layer<20) dist += zLay[layer-1]*mm/cos(theta);
   }
   double tmp = dist/c_light/ns;
   LogDebug("HcalSim") << "HcalTestAnalysis::timeOfFlight " << tmp 
               << " for det/lay " << det << " " << layer 
               << " eta/theta " << eta << " " << theta/deg << " dist " 
               << dist;
   return tmp;
 }

std::vector< int > HcalTestAnalysis::towersToAdd	(	int	centre,
		int	nadd
	)

private

Definition at line 111 of file HcalTestAnalysis.cc.

References funct::abs(), PVValHelper::eta, mps_fire::i, SiStripPI::max, phi, phimax, phimin, and groupFilesInBlocks::temp.

Referenced by HcalTestAnalysis().

                                                                  {
 
   int etac = (centre/100)%100;
   int phic = (centre%100);
   int etamin, etamax, phimin, phimax;
   if (etac>0) {
     etamin = etac-nadd;
     etamax = etac+nadd;
   } else {
     etamin = etac;
     etamax = etac;
   }
   if (phic>0) {
     phimin = phic-nadd;
     phimax = phic+nadd;
   } else {
     phimin = phic;
     phimax = phic;
   }
 
   int nbuf, kount=0;
   nbuf = (etamax-etamin+1)*(phimax-phimin+1);
   std::vector<int> temp(2*nbuf);
   for (int eta=etamin; eta<=etamax; eta++) {
     for (int phi=phimin; phi<=phimax; phi++) {
       temp[kount] = (eta*100 + phi);
       temp[kount+nbuf] = std::max(abs(eta-etac),abs(phi-phic));
       kount++;
     }
   }
 
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Towers to be considered for"
               << " Central " << centre << " and " << nadd 
               << " on either side";
   for (int i=0; i<nbuf; i++)
     edm::LogInfo("HcalSim") << "HcalTestAnalysis: Tower[" << std::setw(3) << i 
                 << "] " << temp[i] << " " << temp[nbuf+i];
   return temp;
 }

void HcalTestAnalysis::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 153 of file HcalTestAnalysis.cc.

References fileName, hcons, names, numberingFromDDD, org, and tuplesManager.

                                                     {
 
   // Numbering From DDD
   edm::ESHandle<HcalDDDSimConstants>    hdc;
   (*job)()->get<HcalSimNumberingRecord>().get(hdc);
   if(!hcons) { hcons = (HcalDDDSimConstants*)(&(*hdc)); }
   //if(!hcons) { hcons = &(*hdc); }
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Initialise "
               << "HcalNumberingFromDDD for " << names[0];
   numberingFromDDD = new HcalNumberingFromDDD(hcons);
 
   // Ntuples
   tuplesManager.reset(new HcalTestHistoManager(fileName));
 
   // Numbering scheme
   org    = new HcalTestNumberingScheme(false);
 
 }

void HcalTestAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 173 of file HcalTestAnalysis.cc.

References HcalDDDSimConstants::cell(), centralTower, HcalCellType::HcalCell::eta, PVValHelper::eta, eta0, HcalBarrel, HcalEndcap, HcalForward, hcons, heppy_loop::loop, names, HcalCellType::HcalCell::ok, org, phi, HcalCellType::HcalCell::phi, phi0, sd, HCalSD::setNumberingScheme(), AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

                                                     {
 
   int irun = (*run)()->GetRunID();
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Begin of Run = " << irun;
 
   bool loop = true, eta = true, phi = true;
   int  etac = (centralTower/100)%100;
   if (etac == 0) {
     etac = 1;
     eta  = false;
   }
   int  phic = (centralTower%100);
   if (phic == 0) {
     phic = 1;
     phi  = false;
   }
   int  idet = static_cast<int>(HcalBarrel);
   while (loop) {
     HcalCellType::HcalCell tmp = hcons->cell(idet,1,1,etac,phic);
     if (tmp.ok) {
       if (eta) eta0 = tmp.eta;
       if (phi) phi0 = tmp.phi;
       loop = false;
     } else if (idet == static_cast<int>(HcalBarrel)) {
       idet = static_cast<int>(HcalEndcap);
     } else if (idet == static_cast<int>(HcalEndcap)) {
       idet = static_cast<int>(HcalForward);
     } else {
       loop = false;
     }
   }
 
   edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Central Tower " 
               << centralTower << " corresponds to eta0 = " << eta0 
               << " phi0 = " << phi0;
  
   std::string sdname = names[0];
   G4SDManager* sd = G4SDManager::GetSDMpointerIfExist();
   if (sd != nullptr) {
     G4VSensitiveDetector* aSD = sd->FindSensitiveDetector(sdname);
     if (aSD==nullptr) {
       edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: No SD with "
                  << "name " << sdname << " in this Setup";
     } else {
       HCalSD* theCaloSD = dynamic_cast<HCalSD*>(aSD);
       edm::LogInfo("HcalSim") << "HcalTestAnalysis::beginOfRun: Finds SD with "
                   << "name " << theCaloSD->GetName() 
                   << " in this Setup";
       if (org) {
         theCaloSD->setNumberingScheme(org);
     edm::LogInfo("HcalSim") << "HcalTestAnalysis::beginOfRun: set a new "
                 << "numbering scheme";
       }
     }
   } else {
     edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: Could not get"
                    << " SD Manager!";
   }
 
 }

void HcalTestAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 235 of file HcalTestAnalysis.cc.

References edepEB, edepEE, edepHB, edepHE, edepHO, edepl, mps_fire::i, LogDebug, mudist, HcalTestHistoClass::setCounters(), and tuples.

                                                       {
  
   // create tuple object
   tuples = new HcalTestHistoClass();
   // Reset counters
   tuples->setCounters();
  
   int i = 0;
   edepEB = edepEE = edepHB = edepHE = edepHO = 0.;
   for (i = 0; i < 20; i++) edepl[i] = 0.;
   for (i = 0; i < 20; i++) mudist[i] = -1.;
 
   int iev = (*evt)()->GetEventID();
   LogDebug("HcalSim") <<"HcalTestAnalysis: Begin of event = " << iev;
 }

void HcalTestAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 320 of file HcalTestAnalysis.cc.

References count, fill(), layerAnalysis(), LogDebug, qieAnalysis(), tuples, and tuplesManager.

                                                     {
 
   count++;
   // Fill event input 
   fill(evt);
   LogDebug("HcalSim") << "HcalTestAnalysis:: ---  after Fill";
 
   // Qie analysis
   CLHEP::HepRandomEngine* engine = G4Random::getTheEngine();
   qieAnalysis(engine);
   LogDebug("HcalSim") << "HcalTestAnalysis:: ---  after QieAnalysis";
 
   // Layers tuples filling
   layerAnalysis();
   LogDebug("HcalSim") << "HcalTestAnalysis:: ---  after LayerAnalysis";
 
   // Writing the data to the Tree
   tuplesManager.get()->fillTree(tuples); // (no need to delete it...)
   tuples = nullptr; // but avoid to reuse it...
   LogDebug("HcalSim") << "HcalTestAnalysis:: --- after fillTree";
 
 }

void HcalTestAnalysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 252 of file HcalTestAnalysis.cc.

References edepEB, edepEE, edepHB, edepHE, edepHO, edepl, PVValHelper::eta, eta0, cmsBatch::log, LogDebug, MeV, mudist, dataset::name, phi, phi0, mathSSE::sqrt(), funct::tan(), and theta().

                                                   {
 
   if (aStep != nullptr) {
     G4VPhysicalVolume* curPV  = aStep->GetPreStepPoint()->GetPhysicalVolume();
     G4String name = curPV->GetName();
     name.assign(name,0,3);
     double edeposit = aStep->GetTotalEnergyDeposit();
     int    layer=-1;
     if (name == "EBR") {
       edepEB += edeposit;
     } else if (name == "EFR") {
       edepEE += edeposit;
     } else if (name == "HBS") {
       layer = (curPV->GetCopyNo()/10)%100;
       if (layer >= 0 && layer < 17) {
     edepHB += edeposit;
       } else {
     edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HB "
                    << curPV->GetName() << curPV->GetCopyNo();
     layer = -1;
       }
     } else if (name == "HES") {
       layer = (curPV->GetCopyNo()/10)%100;
       if (layer >= 0 && layer < 19) {
     edepHE += edeposit;
       } else {
     edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HE " 
                    << curPV->GetName() << curPV->GetCopyNo();
     layer = -1;
       }
     } else if (name == "HTS") {
       layer = (curPV->GetCopyNo()/10)%100;
       if (layer >= 17 && layer < 20) {
     edepHO += edeposit;
        } else {
      edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HO " 
                     << curPV->GetName() << curPV->GetCopyNo();
     layer = -1;
        }
     }
     if (layer >= 0 && layer < 20) {
       edepl[layer] += edeposit;
 
       // Calculate the distance if it is a muon
       G4String part = aStep->GetTrack()->GetDefinition()->GetParticleName();
       if ((part == "mu-" || part == "mu+") && mudist[layer] < 0) {
         math::XYZPoint pos(aStep->GetPreStepPoint()->GetPosition().x(),
                            aStep->GetPreStepPoint()->GetPosition().y(),
                            aStep->GetPreStepPoint()->GetPosition().z());
         double theta   = pos.theta();
         double   eta   = -log(tan(theta * 0.5));
         double   phi   = pos.phi();
         double  dist   = sqrt ((eta-eta0)*(eta-eta0) + (phi-phi0)*(phi-phi0));
         mudist[layer]  = dist*std::sqrt(pos.perp2());
       }
     }
 
     if (layer >= 0 && layer < 20) {
       LogDebug("HcalSim") << "HcalTestAnalysis:: G4Step: " << name << " Layer "
               << std::setw(3) << layer << " Edep " << std::setw(6) 
               << edeposit/MeV << " MeV";
     }
   } else {
     edm::LogInfo("HcalSim") << "HcalTestAnalysis:: G4Step: Null Step";
   }
 }