#include <HcalSimHitsValidation.h>

Inheritance diagram for HcalSimHitsValidation:

Public Member Functions
void	analyze (edm::Event const &ev, edm::EventSetup const &c) override

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

void	endJob () override

	HcalSimHitsValidation (edm::ParameterSet const &conf)

	~HcalSimHitsValidation () override

Public Member Functions inherited from DQMOneEDAnalyzer<>
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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 &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	dPhiWsign (double phi1, double phi2)

double	dR (double eta1, double phi1, double eta2, double phi2)

double	phi12 (double phi1, double en1, double phi2, double en2)

Private Attributes
bool	auxPlots_

std::string	ebHits_

std::string	eeHits_

std::vector< MonitorElement * >	emean_vs_ieta_HB

std::vector< MonitorElement * >	emean_vs_ieta_HE

std::vector< MonitorElement * >	emean_vs_ieta_HF

MonitorElement *	emean_vs_ieta_HO

std::string	g4Label_

const CaloGeometry *	geometry_

std::string	hcalHits_

const HcalDDDRecConstants *	hcons_

double	hf1_

double	hf2_

int	maxDepthHB_

int	maxDepthHE_

int	maxDepthHF_

int	maxDepthHO_

MonitorElement *	meEnConeEtaProfile

MonitorElement *	meEnConeEtaProfile_E

MonitorElement *	meEnConeEtaProfile_EH

std::vector< MonitorElement * >	meSimHitsEnergyHB

std::vector< MonitorElement * >	meSimHitsEnergyHE

std::vector< MonitorElement * >	meSimHitsEnergyHF

MonitorElement *	meSimHitsEnergyHO

int	nevtot

std::vector< MonitorElement * >	Nhb

std::vector< MonitorElement * >	Nhe

std::vector< MonitorElement * >	Nhf

MonitorElement *	Nho

std::vector< MonitorElement * >	occupancy_vs_ieta_HB

std::vector< MonitorElement * >	occupancy_vs_ieta_HE

std::vector< MonitorElement * >	occupancy_vs_ieta_HF

MonitorElement *	occupancy_vs_ieta_HO

std::string	outputFile_

bool	testNumber_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_ecalEB_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_ecalEE_

edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hcal_

edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_HRNDC_

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMOneEDAnalyzer<>
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

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

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 42 of file HcalSimHitsValidation.h.

Constructor & Destructor Documentation

◆ HcalSimHitsValidation()

HcalSimHitsValidation::HcalSimHitsValidation ( edm::ParameterSet const & conf )

Definition at line 4 of file HcalSimHitsValidation.cc.

                                                                         {
   // DQM ROOT output
   outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
   testNumber_ = conf.getUntrackedParameter<bool>("TestNumber", false);
   auxPlots_ = conf.getUntrackedParameter<bool>("auxiliaryPlots", false);
  
   // register for data access
   g4Label_ = conf.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits");
   hcalHits_ = conf.getUntrackedParameter<std::string>("HcalHitCollection", "HcalHits");
   ebHits_ = conf.getUntrackedParameter<std::string>("EBHitCollection", "EcalHitsEB");
   eeHits_ = conf.getUntrackedParameter<std::string>("EEHitCollection", "EcalHitsEE");
  
   // import sampling factors
   hf1_ = conf.getParameter<double>("hf1");
   hf2_ = conf.getParameter<double>("hf2");
  
   tok_evt_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
   tok_hcal_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hcalHits_));
   tok_ecalEB_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, ebHits_));
   tok_ecalEE_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, eeHits_));
   tok_HRNDC_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
  
   if (!outputFile_.empty()) {
     edm::LogVerbatim("OutputInfo") << " Hcal SimHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
   } else {
     edm::LogVerbatim("OutputInfo") << " Hcal SimHit Task histograms will NOT be saved";
   }
  
   nevtot = 0;
 }

References auxPlots_, ebHits_, eeHits_, g4Label_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalHits_, hf1_, hf2_, nevtot, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, testNumber_, tok_ecalEB_, tok_ecalEE_, tok_evt_, tok_geom_, tok_hcal_, and tok_HRNDC_.

◆ ~HcalSimHitsValidation()

HcalSimHitsValidation::~HcalSimHitsValidation ( )

override

Definition at line 36 of file HcalSimHitsValidation.cc.

36 {}

Member Function Documentation

◆ analyze()

void HcalSimHitsValidation::analyze	(	edm::Event const &	ev,
		edm::EventSetup const &	c
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 287 of file HcalSimHitsValidation.cc.

                                                                               {
   //===========================================================================
   // Getting SimHits
   //===========================================================================
  
   double phi_MC = -999.;  // phi of initial particle from HepMC
   double eta_MC = -999.;  // eta of initial particle from HepMC
  
   edm::Handle<edm::HepMCProduct> evtMC;
   ev.getByToken(tok_evt_, evtMC);  // generator in late 310_preX
   if (!evtMC.isValid()) {
     edm::LogVerbatim("OutputInfo") << "no HepMCProduct found";
   }
  
   // MC particle with highest pt is taken as a direction reference
   double maxPt = -99999.;
   int npart = 0;
  
   const HepMC::GenEvent *myGenEvent = evtMC->GetEvent();
   for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
        ++p) {
     double phip = (*p)->momentum().phi();
     double etap = (*p)->momentum().eta();
     double pt = (*p)->momentum().perp();
     if (pt > maxPt) {
       npart++;
       maxPt = pt;
       phi_MC = phip;
       eta_MC = etap;
     }
   }
  
   double partR = 0.3;
  
   // Hcal SimHits
  
   // Approximate calibration constants
   const float calib_HB = 120.;
   const float calib_HE = 190.;
   const float calib_HF1 = hf1_;  // 1.0/0.383;
   const float calib_HF2 = hf2_;  // 1.0/0.368;
  
   edm::Handle<edm::PCaloHitContainer> hcalHits;
   ev.getByToken(tok_hcal_, hcalHits);
   const auto SimHitResult = hcalHits.product();
  
   float eta_diff;
   float etaMax = 9999;
   int ietaMax = 0;
  
   double HcalCone = 0;
  
   geometry_ = &c.getData(tok_geom_);
  
   for (std::vector<PCaloHit>::const_iterator SimHits = SimHitResult->begin(); SimHits != SimHitResult->end();
        ++SimHits) {
     HcalDetId cell;
     if (testNumber_)
       cell = HcalHitRelabeller::relabel(SimHits->id(), hcons_);
     else
       cell = HcalDetId(SimHits->id());
  
     auto cellGeometry = geometry_->getSubdetectorGeometry(cell)->getGeometry(cell);
     double etaS = cellGeometry->getPosition().eta();
     double phiS = cellGeometry->getPosition().phi();
     double en = SimHits->energy();
  
     int sub = cell.subdet();
     int depth = cell.depth();
     double ieta = cell.ieta();
  
     // Energy in Cone
     double r = dR(eta_MC, phi_MC, etaS, phiS);
  
     if (r < partR) {
       eta_diff = std::abs(eta_MC - etaS);
       if (eta_diff < etaMax) {
         etaMax = eta_diff;
         ietaMax = cell.ieta();
       }
       // Approximation of calibration
       if (sub == 1)
         HcalCone += en * calib_HB;
       else if (sub == 2)
         HcalCone += en * calib_HE;
       else if (sub == 4 && (depth == 1 || depth == 3))
         HcalCone += en * calib_HF1;
       else if (sub == 4 && (depth == 2 || depth == 4))
         HcalCone += en * calib_HF2;
     }
  
     if (auxPlots_) {
       // HB
       if (sub == 1) {
         meSimHitsEnergyHB[0]->Fill(en);
         meSimHitsEnergyHB[depth]->Fill(en);
  
         emean_vs_ieta_HB[0]->Fill(double(ieta), en);
         emean_vs_ieta_HB[depth]->Fill(double(ieta), en);
  
         occupancy_vs_ieta_HB[0]->Fill(double(ieta));
         occupancy_vs_ieta_HB[depth]->Fill(double(ieta));
       }
       // HE
       if (sub == 2 && maxDepthHE_ > 0) {
         meSimHitsEnergyHE[0]->Fill(en);
         meSimHitsEnergyHE[depth]->Fill(en);
  
         emean_vs_ieta_HE[0]->Fill(double(ieta), en);
         emean_vs_ieta_HE[depth]->Fill(double(ieta), en);
  
         occupancy_vs_ieta_HE[0]->Fill(double(ieta));
         occupancy_vs_ieta_HE[depth]->Fill(double(ieta));
       }
       // HO
       if (sub == 3) {
         meSimHitsEnergyHO->Fill(en);
  
         emean_vs_ieta_HO->Fill(double(ieta), en);
  
         occupancy_vs_ieta_HO->Fill(double(ieta));
       }
       // HF
       if (sub == 4) {
         meSimHitsEnergyHF[0]->Fill(en);
         meSimHitsEnergyHF[depth]->Fill(en);
  
         emean_vs_ieta_HF[0]->Fill(double(ieta), en);
         emean_vs_ieta_HF[depth]->Fill(double(ieta), en);
  
         occupancy_vs_ieta_HF[0]->Fill(double(ieta));
         occupancy_vs_ieta_HF[depth]->Fill(double(ieta));
       }
  
     }  // auxPlots_
  
   }  // Loop over SimHits
  
   // Ecal EB SimHits
   double EcalCone = 0;
  
   if (!ebHits_.empty()) {
     edm::Handle<edm::PCaloHitContainer> ecalEBHits;
     ev.getByToken(tok_ecalEB_, ecalEBHits);
     const auto SimHitResultEB = ecalEBHits.product();
  
     for (std::vector<PCaloHit>::const_iterator SimHits = SimHitResultEB->begin(); SimHits != SimHitResultEB->end();
          ++SimHits) {
       EBDetId EBid = EBDetId(SimHits->id());
  
       auto cellGeometry = geometry_->getSubdetectorGeometry(EBid)->getGeometry(EBid);
       double etaS = cellGeometry->getPosition().eta();
       double phiS = cellGeometry->getPosition().phi();
       double en = SimHits->energy();
  
       double r = dR(eta_MC, phi_MC, etaS, phiS);
  
       if (r < partR)
         EcalCone += en;
     }
   }  // ebHits_
  
   // Ecal EE SimHits
   if (!eeHits_.empty()) {
     edm::Handle<edm::PCaloHitContainer> ecalEEHits;
     ev.getByToken(tok_ecalEE_, ecalEEHits);
     const auto SimHitResultEE = ecalEEHits.product();
  
     for (std::vector<PCaloHit>::const_iterator SimHits = SimHitResultEE->begin(); SimHits != SimHitResultEE->end();
          ++SimHits) {
       EEDetId EEid = EEDetId(SimHits->id());
  
       auto cellGeometry = geometry_->getSubdetectorGeometry(EEid)->getGeometry(EEid);
       double etaS = cellGeometry->getPosition().eta();
       double phiS = cellGeometry->getPosition().phi();
       double en = SimHits->energy();
  
       double r = dR(eta_MC, phi_MC, etaS, phiS);
  
       if (r < partR)
         EcalCone += en;
     }
   }  // eeHits_
  
   if (ietaMax != 0) {  // If ietaMax == 0, there were no good HCAL SimHits
     meEnConeEtaProfile->Fill(double(ietaMax), HcalCone);
     meEnConeEtaProfile_E->Fill(double(ietaMax), EcalCone);
     meEnConeEtaProfile_EH->Fill(double(ietaMax), HcalCone + EcalCone);
   }
  
   nevtot++;
 }

References funct::abs(), auxPlots_, c, LEDCalibrationChannels::depth, dR(), ebHits_, eeHits_, emean_vs_ieta_HB, emean_vs_ieta_HE, emean_vs_ieta_HF, emean_vs_ieta_HO, ALCARECOTkAlBeamHalo_cff::etaMax, ev, dqm::impl::MonitorElement::Fill(), geometry_, edm::HepMCProduct::GetEvent(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), hcons_, hf1_, hf2_, LEDCalibrationChannels::ieta, edm::HandleBase::isValid(), maxDepthHE_, MuonErrorMatrixAnalyzer_cfi::maxPt, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meSimHitsEnergyHB, meSimHitsEnergyHE, meSimHitsEnergyHF, meSimHitsEnergyHO, nevtot, npart, occupancy_vs_ieta_HB, occupancy_vs_ieta_HE, occupancy_vs_ieta_HF, occupancy_vs_ieta_HO, AlCaHLTBitMon_ParallelJobs::p, edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, alignCSCRings::r, HcalHitRelabeller::relabel(), HcalDetId::subdet(), testNumber_, tok_ecalEB_, tok_ecalEE_, tok_evt_, tok_geom_, and tok_hcal_.

◆ bookHistograms()

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

overridevirtual

Implements DQMOneEDAnalyzer<>.

Definition at line 38 of file HcalSimHitsValidation.cc.

                                                                                                           {
   hcons_ = &es.getData(tok_HRNDC_);
   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 = std::max(hcons_->getEtaRange(1).second, 1);
   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);
  
   Char_t histo[200];
  
   ib.setCurrentFolder("HcalHitsV/HcalSimHitTask");
  
   if (auxPlots_) {
     // General counters
     for (int depth = 0; depth <= maxDepthHB_; depth++) {
       if (depth == 0) {
         sprintf(histo, "N_HB");
       } else {
         sprintf(histo, "N_HB%d", depth);
       }
  
       Nhb.push_back(ib.book1D(histo, histo, 2600, 0., 2600.));
     }
     for (int depth = 0; depth <= maxDepthHE_; depth++) {
       if (depth == 0) {
         sprintf(histo, "N_HE");
       } else {
         sprintf(histo, "N_HE%d", depth);
       }
  
       Nhe.push_back(ib.book1D(histo, histo, 2600, 0., 2600.));
     }
  
     sprintf(histo, "N_HO");
     Nho = ib.book1D(histo, histo, 2200, 0., 2200.);
  
     for (int depth = 0; depth <= maxDepthHF_; depth++) {
       if (depth == 0) {
         sprintf(histo, "N_HF");
       } else {
         sprintf(histo, "N_HF%d", depth);
       }
  
       Nhf.push_back(ib.book1D(histo, histo, 1800, 0., 1800.));
     }
  
     // Mean energy vs iEta TProfiles
     for (int depth = 0; depth <= maxDepthHB_; depth++) {
       if (depth == 0) {
         sprintf(histo, "emean_vs_ieta_HB");
       } else {
         sprintf(histo, "emean_vs_ieta_HB%d", depth);
       }
  
       emean_vs_ieta_HB.push_back(
           ib.bookProfile(histo, histo, ieta_bins_HB, ieta_min_HB, ieta_max_HB, -10., 2000., " "));
     }
     for (int depth = 0; depth <= maxDepthHE_; depth++) {
       if (depth == 0) {
         sprintf(histo, "emean_vs_ieta_HE");
       } else {
         sprintf(histo, "emean_vs_ieta_HE%d", depth);
       }
  
       emean_vs_ieta_HE.push_back(
           ib.bookProfile(histo, histo, ieta_bins_HE, ieta_min_HE, ieta_max_HE, -10., 2000., " "));
     }
  
     sprintf(histo, "emean_vs_ieta_HO");
     emean_vs_ieta_HO = ib.bookProfile(histo, histo, ieta_bins_HO, ieta_min_HO, ieta_max_HO, -10., 2000., " ");
  
     for (int depth = 0; depth <= maxDepthHF_; depth++) {
       if (depth == 0) {
         sprintf(histo, "emean_vs_ieta_HF");
       } else {
         sprintf(histo, "emean_vs_ieta_HF%d", depth);
       }
  
       emean_vs_ieta_HF.push_back(
           ib.bookProfile(histo, histo, ieta_bins_HF, ieta_min_HF, ieta_max_HF, -10., 2000., " "));
     }
  
     // Occupancy vs. iEta TH1Fs
     for (int depth = 0; depth <= maxDepthHB_; depth++) {
       if (depth == 0) {
         sprintf(histo, "occupancy_vs_ieta_HB");
       } else {
         sprintf(histo, "occupancy_vs_ieta_HB%d", depth);
       }
  
       occupancy_vs_ieta_HB.push_back(ib.book1D(histo, histo, ieta_bins_HB, ieta_min_HB, ieta_max_HB));
     }
     for (int depth = 0; depth <= maxDepthHE_; depth++) {
       if (depth == 0) {
         sprintf(histo, "occupancy_vs_ieta_HE");
       } else {
         sprintf(histo, "occupancy_vs_ieta_HE%d", depth);
       }
  
       occupancy_vs_ieta_HE.push_back(ib.book1D(histo, histo, ieta_bins_HE, ieta_min_HE, ieta_max_HE));
     }
  
     sprintf(histo, "occupancy_vs_ieta_HO");
     occupancy_vs_ieta_HO = ib.book1D(histo, histo, ieta_bins_HO, ieta_min_HO, ieta_max_HO);
  
     for (int depth = 0; depth <= maxDepthHF_; depth++) {
       if (depth == 0) {
         sprintf(histo, "occupancy_vs_ieta_HF");
       } else {
         sprintf(histo, "occupancy_vs_ieta_HF%d", depth);
       }
  
       occupancy_vs_ieta_HF.push_back(ib.book1D(histo, histo, ieta_bins_HF, ieta_min_HF, ieta_max_HF));
     }
  
     // Energy spectra
     for (int depth = 0; depth <= maxDepthHB_; depth++) {
       if (depth == 0) {
         sprintf(histo, "HcalSimHitTask_energy_of_simhits_HB");
       } else {
         sprintf(histo, "HcalSimHitTask_energy_of_simhits_HB%d", depth);
       }
  
       meSimHitsEnergyHB.push_back(ib.book1D(histo, histo, 510, -0.1, 5.));
     }
     for (int depth = 0; depth <= maxDepthHE_; depth++) {
       if (depth == 0) {
         sprintf(histo, "HcalSimHitTask_energy_of_simhits_HE");
       } else {
         sprintf(histo, "HcalSimHitTask_energy_of_simhits_HE%d", depth);
       }
  
       meSimHitsEnergyHE.push_back(ib.book1D(histo, histo, 510, -0.1, 5.));
     }
  
     sprintf(histo, "HcalSimHitTask_energy_of_simhits_HO");
     meSimHitsEnergyHO = ib.book1D(histo, histo, 510, -0.1, 5.);
  
     for (int depth = 0; depth <= maxDepthHF_; depth++) {
       if (depth == 0) {
         sprintf(histo, "HcalSimHitTask_energy_of_simhits_HF");
       } else {
         sprintf(histo, "HcalSimHitTask_energy_of_simhits_HF%d", depth);
       }
  
       meSimHitsEnergyHF.push_back(ib.book1D(histo, histo, 1010, -5., 500.));
     }
  
   }  // auxPlots_
  
   // Energy in Cone
   sprintf(histo, "HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths");
   meEnConeEtaProfile = ib.bookProfile(histo, histo, ieta_bins_HF, ieta_min_HF, ieta_max_HF, -10., 200., " ");
  
   sprintf(histo, "HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths_E");
   meEnConeEtaProfile_E = ib.bookProfile(histo, histo, ieta_bins_HF, ieta_min_HF, ieta_max_HF, -10., 200., " ");
  
   sprintf(histo, "HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths_EH");
   meEnConeEtaProfile_EH = ib.bookProfile(histo, histo, ieta_bins_HF, ieta_min_HF, ieta_max_HF, -10., 200., " ");
 }

References auxPlots_, LEDCalibrationChannels::depth, emean_vs_ieta_HB, emean_vs_ieta_HE, emean_vs_ieta_HF, emean_vs_ieta_HO, edm::EventSetup::getData(), HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNPhi(), hcons_, timingPdfMaker::histo, cuy::ib, hcalTTPDigis_cfi::iEtaMax, createfilelist::int, SiStripPI::max, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meSimHitsEnergyHB, meSimHitsEnergyHE, meSimHitsEnergyHF, meSimHitsEnergyHO, Nhb, Nhe, Nhf, Nho, occupancy_vs_ieta_HB, occupancy_vs_ieta_HE, occupancy_vs_ieta_HF, occupancy_vs_ieta_HO, and tok_HRNDC_.

◆ dPhiWsign()

double HcalSimHitsValidation::dPhiWsign	(	double	phi1,
		double	phi2
	)

private

Definition at line 513 of file HcalSimHitsValidation.cc.

                                                                 {
   // clockwise      phi2 w.r.t phi1 means "+" phi distance
   // anti-clockwise phi2 w.r.t phi1 means "-" phi distance
  
   double PI = 3.1415926535898;
   double a1 = phi1;
   double a2 = phi2;
   double tmp = a2 - a1;
   if (a1 * a2 < 0.) {
     if (a1 > 0.5 * PI)
       tmp += 2. * PI;
     if (a2 > 0.5 * PI)
       tmp -= 2. * PI;
   }
   return tmp;
 }

References testProducerWithPsetDescEmpty_cfi::a2, PI, and createJobs::tmp.

◆ dR()

double HcalSimHitsValidation::dR	(	double	eta1,
		double	phi1,
		double	eta2,
		double	phi2
	)

private

Definition at line 480 of file HcalSimHitsValidation.cc.

                                                                                    {
   double PI = 3.1415926535898;
   double deltaphi = phi1 - phi2;
   if (phi2 > phi1) {
     deltaphi = phi2 - phi1;
   }
   if (deltaphi > PI) {
     deltaphi = 2. * PI - deltaphi;
   }
   double deltaeta = eta2 - eta1;
   double tmp = sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
   return tmp;
 }

References HLT_FULL_cff::eta1, HLT_FULL_cff::eta2, mathSSE::sqrt(), and createJobs::tmp.

Referenced by analyze().

◆ endJob()

void HcalSimHitsValidation::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 247 of file HcalSimHitsValidation.cc.

                                    {
   if (auxPlots_) {
     for (int i = 1; i <= occupancy_vs_ieta_HB[0]->getNbinsX(); i++) {
       int ieta = i - 43;  // -41 -1, 1 41
  
       float phi_factor;
  
       if (std::abs(ieta) <= 20)
         phi_factor = 72.;
       else if (std::abs(ieta) < 40)
         phi_factor = 36.;
       else
         phi_factor = 18.;
  
       float cnorm;
  
       // Occupancy vs. iEta TH1Fs
       for (int depth = 0; depth <= maxDepthHB_; depth++) {
         cnorm = occupancy_vs_ieta_HB[depth]->getBinContent(i) / (phi_factor * nevtot);
         occupancy_vs_ieta_HB[depth]->setBinContent(i, cnorm);
       }
       for (int depth = 0; depth <= maxDepthHE_; depth++) {
         cnorm = occupancy_vs_ieta_HE[depth]->getBinContent(i) / (phi_factor * nevtot);
         occupancy_vs_ieta_HE[depth]->setBinContent(i, cnorm);
       }
  
       cnorm = occupancy_vs_ieta_HO->getBinContent(i) / (phi_factor * nevtot);
       occupancy_vs_ieta_HO->setBinContent(i, cnorm);
  
       for (int depth = 0; depth <= maxDepthHF_; depth++) {
         cnorm = occupancy_vs_ieta_HF[depth]->getBinContent(i) / (phi_factor * nevtot);
         occupancy_vs_ieta_HF[depth]->setBinContent(i, cnorm);
       }
     }
   }
  
   // let's see if this breaks anything
   // if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
 }

References funct::abs(), auxPlots_, LEDCalibrationChannels::depth, dqm::impl::MonitorElement::getBinContent(), mps_fire::i, LEDCalibrationChannels::ieta, maxDepthHB_, maxDepthHE_, maxDepthHF_, nevtot, occupancy_vs_ieta_HB, occupancy_vs_ieta_HE, occupancy_vs_ieta_HF, occupancy_vs_ieta_HO, and dqm::impl::MonitorElement::setBinContent().

◆ phi12()

double HcalSimHitsValidation::phi12	(	double	phi1,
		double	en1,
		double	phi2,
		double	en2
	)

private

Definition at line 494 of file HcalSimHitsValidation.cc.

                                                                                     {
   // weighted mean value of phi1 and phi2
  
   double tmp;
   double PI = 3.1415926535898;
   double a1 = phi1;
   double a2 = phi2;
  
   if (a1 > 0.5 * PI && a2 < 0.)
     a2 += 2 * PI;
   if (a2 > 0.5 * PI && a1 < 0.)
     a1 += 2 * PI;
   tmp = (a1 * en1 + a2 * en2) / (en1 + en2);
   if (tmp > PI)
     tmp -= 2. * PI;
  
   return tmp;
 }