#include <SimHitsValidationHcal.h>

Inheritance diagram for SimHitsValidationHcal:

Classes
struct	energysum

struct	etaRange

struct	idType

Public Member Functions
	SimHitsValidationHcal (const edm::ParameterSet &ps)

	~SimHitsValidationHcal ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)

void	analyzeHits (std::vector< PCaloHit > &)

virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Member Functions
std::vector< std::pair < std::string, std::string > >	getHistogramTypes ()

etaRange	getLimits (idType)

std::pair< int, int >	histId (int subdet, int eta, int depth, unsigned int dep)

Private Attributes
std::string	g4Label_

std::string	hcalHits_

const HcalDDDRecConstants *	hcons

bool	initialized

int	maxDepthHB_

int	maxDepthHE_

int	maxDepthHF_

int	maxDepthHO_

MonitorElement *	meEnergy_HB

MonitorElement *	meEnergy_HE

MonitorElement *	meEnergy_HF

MonitorElement *	meEnergy_HO

std::vector< MonitorElement * >	meHcalEnergyl100_

std::vector< MonitorElement * >	meHcalEnergyl250_

std::vector< MonitorElement * >	meHcalEnergyl25_

std::vector< MonitorElement * >	meHcalEnergyl50_

std::vector< MonitorElement * >	meHcalHitEta_

std::vector< MonitorElement * >	meHcalHitTimeEta_

MonitorElement *	metime_enweighted_HB

MonitorElement *	metime_enweighted_HE

MonitorElement *	metime_enweighted_HF

MonitorElement *	metime_enweighted_HO

MonitorElement *	metime_HB

MonitorElement *	metime_HE

MonitorElement *	metime_HF

MonitorElement *	metime_HO

bool	testNumber_

edm::EDGetTokenT < edm::PCaloHitContainer >	tok_hits_

std::vector< idType >	types

bool	verbose_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 29 of file SimHitsValidationHcal.h.

Constructor & Destructor Documentation

SimHitsValidationHcal::SimHitsValidationHcal ( const edm::ParameterSet & ps )

Definition at line 8 of file SimHitsValidationHcal.cc.

References g4Label_, edm::ParameterSet::getParameter(), hcalHits_, HLT_FULL_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, testNumber_, tok_hits_, and verbose_.

                                                                       {
 
   g4Label_  = ps.getParameter<std::string>("ModuleLabel");
   hcalHits_ = ps.getParameter<std::string>("HitCollection");
   verbose_  = ps.getParameter<bool>("Verbose");
   testNumber_= ps.getParameter<bool>("TestNumber");
 
   tok_hits_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_,hcalHits_));
 
   edm::LogInfo("HitsValidationHcal") << "Module Label: " << g4Label_ 
                      << "   Hits: " << hcalHits_
                      << " TestNumbering " << testNumber_;
 }

SimHitsValidationHcal::~SimHitsValidationHcal ( )

Definition at line 22 of file SimHitsValidationHcal.cc.

22 {}

Member Function Documentation

void SimHitsValidationHcal::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

protectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 149 of file SimHitsValidationHcal.cc.

References analyzeHits(), HLT_FULL_cff::depth, eta, edm::EventID::event(), edm::Event::getByToken(), HcalDDDRecConstants::getHCID(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, i, edm::EventBase::id(), edm::HandleBase::isValid(), phi, DetId::rawId(), edm::EventID::run(), jetcorrextractor::sign(), testNumber_, tok_hits_, HcalTestNumbering::unpackHcalIndex(), and z.

                                            {
   
   
 #ifdef DebugLog  
   edm::LogInfo("HitsValidationHcal") << "Run = " << e.id().run() << " Event = "
                      << e.id().event();
 #endif
   std::vector<PCaloHit>               caloHits;
   edm::Handle<edm::PCaloHitContainer> hitsHcal;
   
   bool getHits = false;
   e.getByToken(tok_hits_,hitsHcal); 
   if (hitsHcal.isValid()) getHits = true;
 #ifdef DebugLog  
   edm::LogInfo("HitsValidationHcal") << "HitsValidationHcal.: Input flags Hits "
                      << getHits;
 #endif
   if (getHits) {
     caloHits.insert(caloHits.end(),hitsHcal->begin(),hitsHcal->end());
 #ifdef DebugLog
     edm::LogInfo("HitsValidationHcal") << "testNumber_:" << testNumber_;
 #endif
     if (testNumber_) {
       for (unsigned int i=0; i<caloHits.size(); ++i) {
     unsigned int id_ = caloHits[i].id();
     int subdet, z, depth0, eta0, phi0, lay;
     HcalTestNumbering::unpackHcalIndex(id_, subdet, z, depth0, eta0, phi0, lay);
     int sign = (z==0) ? (-1):(1);
 #ifdef DebugLog
     edm::LogInfo("HitsValidationHcal") << "Hit[" << i 
                        << "] subdet|z|depth|eta|phi|lay " 
                        << subdet << "|" << z << "|" 
                        << depth0 << "|" << eta0 << "|" 
                        << phi0 << "|" << lay;
 #endif
     HcalDDDRecConstants::HcalID id = hcons->getHCID(subdet, eta0, phi0, lay, depth0);
     
     HcalDetId hid;
     if (subdet==int(HcalBarrel)) {
       hid = HcalDetId(HcalBarrel,sign*id.eta,id.phi,id.depth);        
     } else if (subdet==int(HcalEndcap)) {
       hid = HcalDetId(HcalEndcap,sign*id.eta,id.phi,id.depth);    
     } else if (subdet==int(HcalOuter)) {
       hid = HcalDetId(HcalOuter,sign*id.eta,id.phi,id.depth);    
     } else if (subdet==int(HcalForward)) {
       hid = HcalDetId(HcalForward,sign*id.eta,id.phi,id.depth);
     }
     caloHits[i].setID(hid.rawId());
 #ifdef DebugLog
     edm::LogInfo("HitsValidationHcal") << "Hit[" << i << "] " << hid;
 #endif
       }
     }
 #ifdef DebugLog 
     edm::LogInfo("HitsValidationHcal") << "HitsValidationHcal: Hit buffer " 
                        << caloHits.size(); 
 #endif
     analyzeHits (caloHits);
   }
 }

void SimHitsValidationHcal::analyzeHits ( std::vector< PCaloHit > & hits )

protected

Definition at line 211 of file SimHitsValidationHcal.cc.

References TauDecayModes::dec, HLT_FULL_cff::depth, SimHitsValidationHcal::energysum::e100, SimHitsValidationHcal::energysum::e25, SimHitsValidationHcal::energysum::e250, SimHitsValidationHcal::energysum::e50, relval_parameters_module::energy, eta, MonitorElement::Fill(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, histId(), i, meEnergy_HB, meEnergy_HE, meEnergy_HF, meEnergy_HO, meHcalEnergyl100_, meHcalEnergyl250_, meHcalEnergyl25_, meHcalEnergyl50_, meHcalHitEta_, meHcalHitTimeEta_, metime_enweighted_HB, metime_enweighted_HE, metime_enweighted_HF, metime_enweighted_HO, metime_HB, metime_HE, metime_HF, metime_HO, phi, and types.

Referenced by analyze().

                                                                   {
 
   int nHit = hits.size();
   double entotHB = 0, entotHE = 0, entotHF = 0, entotHO = 0; 
   double timetotHB = 0, timetotHE = 0, timetotHF = 0, timetotHO = 0; 
   int    nHB=0, nHE=0, nHO=0, nHF=0;
   
   std::map<std::pair<HcalDetId,unsigned int>,energysum> map_try;
   map_try.clear();
   std::map<std::pair<HcalDetId,unsigned int>,energysum>::iterator itr;
   
   for (int i=0; i<nHit; i++) {
     double energy    = hits[i].energy();
     double time      = hits[i].time();
     HcalDetId id     = HcalDetId(hits[i].id());
     int itime        = (int)(time);
     int subdet       = id.subdet();
     int depth        = id.depth();
     int eta          = id.ieta();
     unsigned int dep = hits[i].depth();
         
     std::pair<int,int> types = histId(subdet, eta, depth, dep);
     if (subdet == static_cast<int>(HcalBarrel)) {
       entotHB += energy;
       timetotHB += time;
       nHB++;
     } else if (subdet == static_cast<int>(HcalEndcap)) {
       entotHE += energy;
       timetotHE += time;
       nHE++;
     } else if (subdet == static_cast<int>(HcalOuter)) {
       entotHO += energy;
       timetotHO += time;
       nHO++;
     } else if (subdet == static_cast<int>(HcalForward)) {
       entotHF += energy;
       timetotHF += time;
       nHF++;
     }
 
     std::pair<HcalDetId,unsigned int> id0(id,dep);
     energysum ensum;
     if (map_try.count(id0) != 0) ensum =  map_try[id0];
     if (itime<250) {
       ensum.e250 += energy;
       if (itime<100) {
     ensum.e100 += energy;
     if (itime<50) {
       ensum.e50 += energy;
       if (itime<25) ensum.e25 += energy;
     }
       }
     }
     map_try[id0] = ensum;
     
 #ifdef DebugLog
     edm::LogInfo("HitsValidationHcal") << "Hit[" << i << "] ID " << std::dec 
                        << " " << id << std::dec << " Det " 
                        << id.det() << " Sub " << subdet 
                        << " depth " << depth << " depthX "
                        << dep << " Eta " << eta << " Phi " 
                        << id.iphi() << " E " << energy 
                        << " time " << time << " type " 
                        << types.first << " " << types.second;
 #endif
     
     double etax = eta - 0.5;
     if (eta < 0) etax += 1;
     if (types.first >= 0) {
       meHcalHitEta_[types.first]->Fill(etax,energy);
       meHcalHitTimeEta_[types.first]->Fill(etax,time);
     }
     if (types.second >= 0) {
       meHcalHitEta_[types.second]->Fill(etax,energy);
       meHcalHitTimeEta_[types.second]->Fill(etax,time);
     }
   }
   
   meEnergy_HB->Fill(entotHB);
   meEnergy_HE->Fill(entotHE);
   meEnergy_HF->Fill(entotHF);
   meEnergy_HO->Fill(entotHO);
 
   metime_HB->Fill(timetotHB);
   metime_HE->Fill(timetotHE);
   metime_HF->Fill(timetotHF);
   metime_HO->Fill(timetotHO);
   
   metime_enweighted_HB->Fill(timetotHB,entotHB);
   metime_enweighted_HE->Fill(timetotHE,entotHE);
   metime_enweighted_HF->Fill(timetotHF,entotHF);
   metime_enweighted_HO->Fill(timetotHO,entotHO);
   
   
   for (itr = map_try.begin(); itr != map_try.end(); ++itr) {
     HcalDetId id    = (*itr).first.first;
     energysum ensum = (*itr).second;
     std::pair<int,int> types = histId((int)(id.subdet()), id.ieta(), id.depth(), (*itr).first.second);
     int eta = id.ieta();
     int phi = id.iphi();
     double etax= eta-0.5;
     double phix= phi-0.5;
     if (types.first >= 0) {
       meHcalEnergyl25_[types.first]->Fill(etax,phix,ensum.e25);
       meHcalEnergyl50_[types.first]->Fill(etax,phix,ensum.e50);
       meHcalEnergyl100_[types.first]->Fill(etax,phix,ensum.e100);
       meHcalEnergyl250_[types.first]->Fill(etax,phix,ensum.e250);
     }
     if (types.second >= 0) {
       meHcalEnergyl25_[types.second]->Fill(etax,phix,ensum.e25);
       meHcalEnergyl50_[types.second]->Fill(etax,phix,ensum.e50);
       meHcalEnergyl100_[types.second]->Fill(etax,phix,ensum.e100);
       meHcalEnergyl250_[types.second]->Fill(etax,phix,ensum.e250);
     }
 
 #ifdef DebugLog
     edm::LogInfo("HitsValidationHcal") << " energy of tower ="   << (*itr).first.first 
                        << " in time 25ns is == " << (*itr).second.e25 
                        << " in time 25-50ns == " << (*itr).second.e50 
                        << " in time 50-100ns == " << (*itr).second.e100 
                        << " in time 100-250 ns == " << (*itr).second.e250;
 #endif
   }
   
 }

void SimHitsValidationHcal::bookHistograms	(	DQMStore::IBooker &	ib,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

protectedvirtual

Implements DQMEDAnalyzer.

Definition at line 24 of file SimHitsValidationHcal.cc.

References SimHitsValidationHcal::etaRange::bins, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), plotBeamSpotDB::first, edm::EventSetup::get(), HcalDDDRecConstants::getEtaRange(), getHistogramTypes(), getLimits(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNPhi(), hcons, SimHitsValidationHcal::etaRange::high, i, MessageLogger_cff::limit, SimHitsValidationHcal::etaRange::low, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meEnergy_HB, meEnergy_HE, meEnergy_HF, meEnergy_HO, meHcalEnergyl100_, meHcalEnergyl250_, meHcalEnergyl25_, meHcalEnergyl50_, meHcalHitEta_, meHcalHitTimeEta_, metime_enweighted_HB, metime_enweighted_HE, metime_enweighted_HF, metime_enweighted_HO, metime_HB, metime_HE, metime_HF, metime_HO, mergeVDriftHistosByStation::name, edm::second(), DQMStore::IBooker::setCurrentFolder(), SiPixelPhase1TrackClusters_cfi::title, and types.

                                                                                                           {
   
   edm::ESHandle<HcalDDDRecConstants> pHRNDC;
   es.get<HcalRecNumberingRecord>().get( pHRNDC );
   hcons = &(*pHRNDC);
   maxDepthHB_ = hcons->getMaxDepth(0);
   maxDepthHE_ = hcons->getMaxDepth(1);
   maxDepthHF_ = hcons->getMaxDepth(2);
   maxDepthHO_ = hcons->getMaxDepth(3);
 
   //Get Phi segmentation from geometry, use the max phi number so that all iphi values are included.
 
   int NphiMax = hcons->getNPhi(0);
 
   NphiMax = (hcons->getNPhi(1) > NphiMax ? hcons->getNPhi(1) : NphiMax);
   NphiMax = (hcons->getNPhi(2) > NphiMax ? hcons->getNPhi(2) : NphiMax);
   NphiMax = (hcons->getNPhi(3) > NphiMax ? hcons->getNPhi(3) : NphiMax);
 
   //Center the iphi bins on the integers
   float iphi_min = 0.5;
   float iphi_max = NphiMax + 0.5;
   int iphi_bins = (int) (iphi_max - iphi_min);
 
   int iEtaHBMax = hcons->getEtaRange(0).second;
   int iEtaHEMax = hcons->getEtaRange(1).second;
   int iEtaHFMax = hcons->getEtaRange(2).second;
   int iEtaHOMax = hcons->getEtaRange(3).second;
 
   //Retain classic behavior, all plots have same ieta range.
   //Comment out code to allow each subdetector to have its on range
 
   int iEtaMax = (iEtaHBMax > iEtaHEMax ? iEtaHBMax : iEtaHEMax);
   iEtaMax = (iEtaMax > iEtaHFMax ? iEtaMax : iEtaHFMax);
   iEtaMax = (iEtaMax > iEtaHOMax ? iEtaMax : iEtaHOMax);
 
   iEtaHBMax = iEtaMax;
   iEtaHEMax = iEtaMax;
   iEtaHFMax = iEtaMax;
   iEtaHOMax = iEtaMax;
 
   //Give an empty bin around the subdet ieta range to make it clear that all ieta rings have been included
   //float ieta_min_HB = -iEtaHBMax - 1.5;
   //float ieta_max_HB = iEtaHBMax + 1.5;
   //int ieta_bins_HB = (int) (ieta_max_HB - ieta_min_HB);
 
   //float ieta_min_HE = -iEtaHEMax - 1.5;
   //float ieta_max_HE = iEtaHEMax + 1.5;
   //int ieta_bins_HE = (int) (ieta_max_HE - ieta_min_HE);
 
   //float ieta_min_HF = -iEtaHFMax - 1.5;
   //float ieta_max_HF = iEtaHFMax + 1.5;
   //int ieta_bins_HF = (int) (ieta_max_HF - ieta_min_HF);
 
   //float ieta_min_HO = -iEtaHOMax - 1.5;
   //float ieta_max_HO = iEtaHOMax + 1.5;
   //int ieta_bins_HO = (int) (ieta_max_HO - ieta_min_HO);
 
 #ifdef DebugLog
   edm::LogInfo("HitsValidationHcal") << " Maximum Depths HB:"<< maxDepthHB_ 
                      << " HE:" << maxDepthHE_  << " HO:" 
                      << maxDepthHO_ << " HF:" << maxDepthHF_;
 #endif  
   std::vector<std::pair<std::string,std::string> > divisions = getHistogramTypes();
 
   edm::LogInfo("HitsValidationHcal") << "Booking the Histograms";
   ib.setCurrentFolder("HcalHitsV/SimHitsValidationHcal");
 
   //Histograms for Hits
 
   char name[100], title[200];
   for (unsigned int i=0; i<types.size(); ++i) {
     etaRange limit = getLimits(types[i]);
     sprintf (name, "HcalHitEta%s", divisions[i].first.c_str());
     sprintf (title, "Hit energy as a function of eta tower index in %s", divisions[i].second.c_str());
     meHcalHitEta_.push_back(ib.book1D(name, title, limit.bins, limit.low, limit.high));
       
     sprintf (name, "HcalHitTimeAEta%s", divisions[i].first.c_str());
     sprintf (title, "Hit time as a function of eta tower index in %s", divisions[i].second.c_str());
     meHcalHitTimeEta_.push_back(ib.book1D(name, title, limit.bins, limit.low, limit.high));
       
     sprintf (name, "HcalHitE25%s", divisions[i].first.c_str());
     sprintf (title, "Energy in time window 0 to 25 for a tower in %s", divisions[i].second.c_str());
     meHcalEnergyl25_.push_back(ib.book2D(name, title, limit.bins, limit.low, limit.high, iphi_bins, iphi_min, iphi_max));
       
     sprintf (name, "HcalHitE50%s", divisions[i].first.c_str());
     sprintf (title, "Energy in time window 0 to 50 for a tower in %s", divisions[i].second.c_str());
     meHcalEnergyl50_.push_back(ib.book2D(name, title, limit.bins, limit.low, limit.high, iphi_bins, iphi_min, iphi_max));
       
     sprintf (name, "HcalHitE100%s", divisions[i].first.c_str());
     sprintf (title, "Energy in time window 0 to 100 for a tower in %s", divisions[i].second.c_str());
     meHcalEnergyl100_.push_back(ib.book2D(name, title, limit.bins, limit.low, limit.high, iphi_bins, iphi_min, iphi_max));
       
     sprintf (name, "HcalHitE250%s", divisions[i].first.c_str());
     sprintf (title, "Energy in time window 0 to 250 for a tower in %s", divisions[i].second.c_str());
     meHcalEnergyl250_.push_back(ib.book2D(name, title, limit.bins, limit.low, limit.high, iphi_bins, iphi_min, iphi_max));
   }
 
   sprintf (name, "Energy_HB");
   meEnergy_HB = ib.book1D(name, name, 100,0,1);
   sprintf (name, "Energy_HE");
   meEnergy_HE = ib.book1D(name, name, 100,0,1);
   sprintf (name, "Energy_HO");
   meEnergy_HO = ib.book1D(name, name, 100,0,1);
   sprintf (name, "Energy_HF");
   meEnergy_HF = ib.book1D(name, name, 100,0,50);
       
   sprintf (name, "Time_HB");
   metime_HB = ib.book1D(name, name, 300,-150,150);
   sprintf (name, "Time_HE");
   metime_HE = ib.book1D(name, name, 300,-150,150);
   sprintf (name, "Time_HO");
   metime_HO = ib.book1D(name, name, 300,-150, 150);
   sprintf (name, "Time_HF");
   metime_HF = ib.book1D(name, name, 300,-150,150);
 
   sprintf (name, "Time_Enweighted_HB");
   metime_enweighted_HB = ib.book1D(name, name, 300,-150,150);
   sprintf (name, "Time_Enweighted_HE");
   metime_enweighted_HE = ib.book1D(name, name, 300,-150,150);
   sprintf (name, "Time_Enweighted_HO");
   metime_enweighted_HO = ib.book1D(name, name, 300,-150, 150);
   sprintf (name, "Time_Enweighted_HF");
   metime_enweighted_HF = ib.book1D(name, name, 300,-150,150);
 }

void SimHitsValidationHcal::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 495 of file SimHitsValidationHcal.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addDefault(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                        {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("ModuleLabel","g4SimHits");
   desc.add<std::string>("HitCollection","HcalHits");
   desc.add<bool>("Verbose",false);
   desc.add<bool>("TestNumber",false);
 
   descriptions.addDefault(desc);
 }

std::vector< std::pair< std::string, std::string > > SimHitsValidationHcal::getHistogramTypes ( )

private

Definition at line 418 of file SimHitsValidationHcal.cc.

References HLT_FULL_cff::depth, HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalOuter, relval_2017::k, bookConverter::max, HLT_FULL_cff::maxDepth, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, cscdqm::h::names, AlCaHLTBitMon_QueryRunRegistry::string, and types.

Referenced by bookHistograms().

                                                                                   {
 
   int maxDepth = std::max(maxDepthHB_,maxDepthHE_);
   maxDepth = std::max(maxDepth,maxDepthHF_);
   maxDepth = std::max(maxDepth,maxDepthHO_);
 
   std::vector<std::pair<std::string,std::string> > divisions;
   std::pair<std::string,std::string> names;
   char                               name1[40], name2[40];
   SimHitsValidationHcal::idType      type;
   //first overall Hcal
   for (int depth=0; depth<maxDepth; ++depth) {
     snprintf (name1, 40, "HC%d", depth);
     snprintf (name2, 40, "HCAL depth%d", depth+1);
     names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
     type  = SimHitsValidationHcal::idType(HcalEmpty,0,depth+1,depth+1);
     divisions.push_back(names);
     types.push_back(type);
   }
   //HB
   for (int depth=0; depth<maxDepthHB_; ++depth) {
     snprintf (name1, 40, "HB%d", depth);
     snprintf (name2, 40, "HB depth%d", depth+1);
     names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
     type  = SimHitsValidationHcal::idType(HcalBarrel,0,depth+1,depth+1);
     divisions.push_back(names);
     types.push_back(type);
   }
   //HE
   for (int depth=0; depth<maxDepthHE_; ++depth) {
     snprintf (name1, 40, "HE%d+z", depth);
     snprintf (name2, 40, "HE +z depth%d", depth+1);
     names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
     type  = SimHitsValidationHcal::idType(HcalEndcap,1,depth+1,depth+1);
     divisions.push_back(names);
     types.push_back(type);
     snprintf (name1, 40, "HE%d-z", depth);
     snprintf (name2, 40, "HE -z depth%d", depth+1);
     names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
     type  = SimHitsValidationHcal::idType(HcalEndcap,-1,depth+1,depth+1);
     divisions.push_back(names);
     types.push_back(type);
   }
   //HO
   {
     int depth = maxDepthHO_;
     snprintf (name1, 40, "HO%d", depth);
     snprintf (name2, 40, "HO depth%d", depth);
     names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
     type  = SimHitsValidationHcal::idType(HcalOuter,0,depth,depth);
     divisions.push_back(names);
     types.push_back(type);
   }
   //HF (first absorber, then different types of abnormal hits)
   std::string hfty1[4] = {"A","W","B","J"};
   std::string hfty2[4] = {"Absorber", "Window", "Bundle", "Jungle"};
   int         dept0[4] = {0, 1, 2, 3};
   for (int k=0; k<4; ++k) {
     for (int depth=0; depth<maxDepthHF_; ++depth) {
       snprintf (name1, 40, "HF%s%d+z", hfty1[k].c_str(), depth);
       snprintf (name2, 40, "HF (%s) +z depth%d", hfty2[k].c_str(), depth+1);
       names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
       type  = SimHitsValidationHcal::idType(HcalForward,1,depth+1,dept0[k]);
       divisions.push_back(names);
       types.push_back(type);
       snprintf (name1, 40, "HF%s%d-z", hfty1[k].c_str(), depth);
       snprintf (name2, 40, "HF (%s) -z depth%d", hfty2[k].c_str(), depth+1);
       names = std::pair<std::string,std::string>(std::string(name1),std::string(name2));
       type  = SimHitsValidationHcal::idType(HcalForward,-1,depth+1,dept0[k]);
       divisions.push_back(names);
       types.push_back(type);
     }
   }
 
   return divisions;
 }

SimHitsValidationHcal::etaRange SimHitsValidationHcal::getLimits ( idType type )

private

Definition at line 337 of file SimHitsValidationHcal.cc.

References MuonErrorMatrixValues_cff::etaRange, HcalDDDRecConstants::getEtaRange(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, SimHitsValidationHcal::idType::subdet, and SimHitsValidationHcal::idType::z.

Referenced by bookHistograms().

                                                                           {
 
   int    bins;
   std::pair<int,int> range;
   double low, high;
   
   if (type.subdet == HcalBarrel) {
     range = hcons->getEtaRange(0);
     low  =-range.second;
     high = range.second;
     bins = (high-low);
   } else if (type.subdet == HcalEndcap) {
     range = hcons->getEtaRange(1);
     bins = range.second- range.first;
     if (type.z == 1) {
       low  = range.first;
       high = range.second;
     } else {
       low  =-range.second;
       high =-range.first;
     }
   } else if (type.subdet == HcalOuter) {
     range = hcons->getEtaRange(3);
     low  =-range.second;
     high = range.second;
     bins = high-low;
   } else if (type.subdet == HcalForward) {
     range = hcons->getEtaRange(2);
     bins = range.second-range.first;
     if (type.z == 1) {
       low  = range.first;
       high = range.second;
     } else {
       low  =-range.second;
       high =-range.first;
     }
   } else {
     bins = 82;
     low  =-41;
     high = 41;
   }
 #ifdef DebugLog
   edm::LogInfo("HitsValidationHcal") << "Subdetector:" << type.subdet << " z:" 
                      << type.z << " range.first:" << range.first
                      << " and second:" << range.second;
   edm::LogInfo("HitsValidationHcal") << "bins: " << bins << " low:" << low 
                      << " high:" << high;
 #endif
   return SimHitsValidationHcal::etaRange(bins, low, high);
 }

std::pair< int, int > SimHitsValidationHcal::histId	(	int	subdet,
		int	eta,
		int	depth,
		unsigned int	dep
	)

private

Definition at line 388 of file SimHitsValidationHcal.cc.

References HcalEmpty, HcalEndcap, HcalForward, relval_2017::k, types, and z.

Referenced by analyzeHits().

                                                                                                {
 
   int id1(-1), id2(-1);
   for (unsigned int k=0; k<types.size(); ++k) {
     if (subdet == HcalForward) {
       if (subdet == (int)(types[k].subdet) && depth == types[k].depth1 &&
       eta*types[k].z > 0 && dep == (unsigned int)(types[k].depth2)) {
     id1 = k; break;
       }
     } else if (subdet == HcalEndcap) {
       if (subdet == (int)(types[k].subdet) && depth == types[k].depth1 &&
       eta*types[k].z > 0) {
     id1 = k; break;
       }
     } else {
       if (subdet == (int)(types[k].subdet) && depth == types[k].depth1) {
     id1 = k; break;
       }
     }
   }
   if (subdet == HcalForward) depth += 2*dep;
   for (unsigned int k=0; k<types.size(); ++k) {
     if (types[k].subdet == HcalEmpty && types[k].depth1 == depth) {
       id2 = k; 
       break;
     }
   }
   return std::pair<int,int>(id1,id2);
 }