HcalRecHitsAnalyzer Class Reference

#include <HcalRecHitsAnalyzer.h>

Inheritance diagram for HcalRecHitsAnalyzer:

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	dqmBeginRun (const edm::Run &run, const edm::EventSetup &c) override
	HcalRecHitsAnalyzer (edm::ParameterSet const &conf)
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) override
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override
	~DQMEDAnalyzer () override=default
Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () 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
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)
	~ProducerBase () 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 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 &	itemsToGetFrom (BranchType iType) 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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	dPhiWsign (double phi1, double phi2)
double	dR (double eta1, double phi1, double eta2, double phi2)
virtual void	fillRecHitsTmp (int subdet_, edm::Event const &ev)
int	hcalSevLvl (const CaloRecHit *hit)
double	phi12 (double phi1, double en1, double phi2, double en2)

Private Attributes
std::vector< uint32_t >	cauxstwd
std::vector< double >	cchi2
std::vector< int >	cdepth
std::vector< double >	cen
std::vector< double >	cenM0
std::vector< double >	cenM3
std::vector< double >	ceta
std::vector< int >	cieta
std::vector< int >	ciphi
std::vector< double >	cphi
std::vector< uint32_t >	cstwd
std::vector< int >	csub
std::vector< double >	ctime
std::vector< double >	cz
std::string	ecalselector_
std::vector< MonitorElement * >	emap
MonitorElement *	emap_HO
std::vector< MonitorElement * >	emean_vs_ieta_HB
std::vector< MonitorElement * >	emean_vs_ieta_HBM0
std::vector< MonitorElement * >	emean_vs_ieta_HBM3
std::vector< MonitorElement * >	emean_vs_ieta_HE
std::vector< MonitorElement * >	emean_vs_ieta_HEM0
std::vector< MonitorElement * >	emean_vs_ieta_HEM3
std::vector< MonitorElement * >	emean_vs_ieta_HEP17
std::vector< MonitorElement * >	emean_vs_ieta_HEP17M0
std::vector< MonitorElement * >	emean_vs_ieta_HEP17M3
std::vector< MonitorElement * >	emean_vs_ieta_HF
MonitorElement *	emean_vs_ieta_HO
int	etype_
std::string	eventype_
bool	famos_
edm::ESHandle< CaloGeometry >	geometry
std::vector< int >	hcalHBSevLvlVec
std::vector< int >	hcalHESevLvlVec
std::vector< int >	hcalHFSevLvlVec
std::vector< int >	hcalHOSevLvlVec
std::string	hcalselector_
const HcalDDDRecConstants *	hcons
bool	hep17_
int	ieta_bins_
float	ieta_max_
float	ieta_min_
int	imc
int	iphi_bins_
float	iphi_max_
float	iphi_min_
int	iz
MonitorElement *	map_ecal
int	maxDepthAll_
int	maxDepthHB_
int	maxDepthHE_
int	maxDepthHF_
int	maxDepthHO_
std::string	mc_
MonitorElement *	meDeltaEta
MonitorElement *	meDeltaPhi
MonitorElement *	meEcalHcalEnergyConeHB
MonitorElement *	meEcalHcalEnergyConeHE
MonitorElement *	meEcalHcalEnergyConeHF
MonitorElement *	meEcalHcalEnergyConeHO
MonitorElement *	meEcalHcalEnergyHB
MonitorElement *	meEcalHcalEnergyHE
MonitorElement *	meEnConeEtaProfile
MonitorElement *	meEnConeEtaProfile_E
MonitorElement *	meEnConeEtaProfile_EH
MonitorElement *	meEnergyHcalVsEcalHB
MonitorElement *	meEnergyHcalVsEcalHE
MonitorElement *	meLog10Chi2profileHB
MonitorElement *	meLog10Chi2profileHE
MonitorElement *	meNumEcalRecHitsConeHB
MonitorElement *	meNumEcalRecHitsConeHE
MonitorElement *	meRecHitsEnergyHB
MonitorElement *	meRecHitsEnergyHBM0
MonitorElement *	meRecHitsEnergyHBM3
MonitorElement *	meRecHitsEnergyHE
MonitorElement *	meRecHitsEnergyHEM0
MonitorElement *	meRecHitsEnergyHEM3
std::vector< MonitorElement * >	meRecHitsEnergyHEP17
std::vector< MonitorElement * >	meRecHitsEnergyHEP17M0
std::vector< MonitorElement * >	meRecHitsEnergyHEP17M3
MonitorElement *	meRecHitsEnergyHF
MonitorElement *	meRecHitsEnergyHO
MonitorElement *	meRecHitsEnergyM2vM0HB
MonitorElement *	meRecHitsEnergyM2vM0HE
MonitorElement *	meRecHitsEnergyM3vM0HB
MonitorElement *	meRecHitsEnergyM3vM0HE
MonitorElement *	meRecHitsEnergyM3vM2HB
MonitorElement *	meRecHitsEnergyM3vM2HE
MonitorElement *	meRecHitsM2Chi2HB
MonitorElement *	meRecHitsM2Chi2HE
MonitorElement *	meSumRecHitsEnergyConeHB
MonitorElement *	meSumRecHitsEnergyConeHE
MonitorElement *	meSumRecHitsEnergyConeHF
MonitorElement *	meSumRecHitsEnergyConeHFL
MonitorElement *	meSumRecHitsEnergyConeHFS
MonitorElement *	meSumRecHitsEnergyConeHO
MonitorElement *	meSumRecHitsEnergyHB
MonitorElement *	meSumRecHitsEnergyHE
MonitorElement *	meSumRecHitsEnergyHF
MonitorElement *	meSumRecHitsEnergyHO
MonitorElement *	meTE_HB
MonitorElement *	meTE_HE
MonitorElement *	meTE_HF
MonitorElement *	meTE_High_HB
MonitorElement *	meTE_High_HO
MonitorElement *	meTE_HO
MonitorElement *	meTE_Low_HB
MonitorElement *	meTE_Low_HE
MonitorElement *	meTE_Low_HF
MonitorElement *	meTEprofileHB
MonitorElement *	meTEprofileHB_High
MonitorElement *	meTEprofileHB_Low
MonitorElement *	meTEprofileHE
MonitorElement *	meTEprofileHE_Low
MonitorElement *	meTEprofileHF
MonitorElement *	meTEprofileHF_Low
MonitorElement *	meTEprofileHO
MonitorElement *	meTEprofileHO_High
MonitorElement *	meTimeHB
MonitorElement *	meTimeHE
MonitorElement *	meTimeHF
MonitorElement *	meTimeHO
int	nChannels_ [5]
int	nevtot
std::vector< MonitorElement * >	Nhb
std::vector< MonitorElement * >	Nhe
std::vector< MonitorElement * >	Nhf
std::vector< MonitorElement * >	Nho
std::vector< MonitorElement * >	nrechits_vs_iphi_HBM
std::vector< MonitorElement * >	nrechits_vs_iphi_HBP
std::vector< MonitorElement * >	nrechits_vs_iphi_HEM
std::vector< MonitorElement * >	nrechits_vs_iphi_HEP
std::vector< MonitorElement * >	nrechits_vs_iphi_HFM
std::vector< MonitorElement * >	nrechits_vs_iphi_HFP
MonitorElement *	nrechits_vs_iphi_HOM
MonitorElement *	nrechits_vs_iphi_HOP
std::vector< MonitorElement * >	occupancy_map_HB
std::vector< MonitorElement * >	occupancy_map_HE
std::vector< MonitorElement * >	occupancy_map_HF
MonitorElement *	occupancy_map_HO
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_
MonitorElement *	RecHit_Aux_StatusWord_HB
MonitorElement *	RecHit_Aux_StatusWord_HE
MonitorElement *	RecHit_Aux_StatusWord_HF
MonitorElement *	RecHit_Aux_StatusWord_HO
MonitorElement *	RecHit_StatusWord_HB
MonitorElement *	RecHit_StatusWord_HE
MonitorElement *	RecHit_StatusWord_HF
MonitorElement *	RecHit_StatusWord_HF67
MonitorElement *	RecHit_StatusWord_HO
MonitorElement *	RecHit_StatusWordCorr_HB
MonitorElement *	RecHit_StatusWordCorr_HE
MonitorElement *	sevLvl_HB
MonitorElement *	sevLvl_HE
MonitorElement *	sevLvl_HF
MonitorElement *	sevLvl_HO
std::string	sign_
int	subdet_
const HcalChannelQuality *	theHcalChStatus
const HcalSeverityLevelComputer *	theHcalSevLvlComputer
const HcalTopology *	theHcalTopology
edm::EDGetTokenT< EBRecHitCollection >	tok_EB_
edm::EDGetTokenT< EERecHitCollection >	tok_EE_
edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_
edm::EDGetTokenT< HFRecHitCollection >	tok_hf_
edm::EDGetTokenT< HORecHitCollection >	tok_ho_
std::string	topFolderName_

Additional Inherited Members
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 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)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 58 of file HcalRecHitsAnalyzer.h.

Constructor & Destructor Documentation

HcalRecHitsAnalyzer::HcalRecHitsAnalyzer

(

edm::ParameterSet const &

conf

)

Definition at line 6 of file HcalRecHitsAnalyzer.cc.

References isolatedTracks_cfi::EBRecHitCollectionLabel, ecalselector_, isolatedTracks_cfi::EERecHitCollectionLabel, etype_, eventype_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalselector_, hep17_, imc, iz, nevtot, outputFile_, sign_, AlCaHLTBitMon_QueryRunRegistry::string, subdet_, tok_EB_, tok_EE_, tok_hbhe_, tok_hf_, and tok_ho_.

  : topFolderName_           ( conf.getParameter<std::string>("TopFolderName") )
{

  // DQM ROOT output
  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
  
  if ( !outputFile_.empty() ) {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
  } else {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
  }
  
  nevtot = 0;
 
  hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");
  ecalselector_ = conf.getUntrackedParameter<std::string>("ecalselector", "yes");
  eventype_     = conf.getUntrackedParameter<std::string>("eventype", "single");
  sign_         = conf.getUntrackedParameter<std::string>("sign", "*");
  //useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);

  //HEP17 configuration
  hep17_        = conf.getUntrackedParameter<bool>("hep17");

  //Collections
  tok_hbhe_ = consumes<HBHERecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HBHERecHitCollectionLabel"));
  tok_hf_  = consumes<HFRecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HFRecHitCollectionLabel"));
  tok_ho_ = consumes<HORecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HORecHitCollectionLabel"));
  edm::InputTag EBRecHitCollectionLabel = conf.getParameter<edm::InputTag>("EBRecHitCollectionLabel");
  tok_EB_ = consumes<EBRecHitCollection>(EBRecHitCollectionLabel);
  edm::InputTag EERecHitCollectionLabel = conf.getParameter<edm::InputTag>("EERecHitCollectionLabel");
  tok_EE_ = consumes<EERecHitCollection>(EERecHitCollectionLabel);

  subdet_ = 5;
  if (hcalselector_ == "noise") subdet_ = 0;
  if (hcalselector_ == "HB"   ) subdet_ = 1;
  if (hcalselector_ == "HE"   ) subdet_ = 2;
  if (hcalselector_ == "HO"   ) subdet_ = 3;
  if (hcalselector_ == "HF"   ) subdet_ = 4;
  if (hcalselector_ == "all"  ) subdet_ = 5;
  if (hcalselector_ == "ZS"   ) subdet_ = 6;

  etype_ = 1;
  if (eventype_ == "multi") etype_ = 2;

  iz = 1;
  if(sign_ == "-") iz = -1;
  if(sign_ == "*") iz = 0;

  imc = 0;

  }

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 507 of file HcalRecHitsAnalyzer.cc.

References cauxstwd, cchi2, cdepth, cen, cenM0, cenM3, ceta, vertices_cff::chi2, cieta, ciphi, cphi, cstwd, csub, ctime, egammaForCoreTracking_cff::depth, dR(), ecalselector_, emap, emap_HO, emean_vs_ieta_HB, emean_vs_ieta_HBM0, emean_vs_ieta_HBM3, emean_vs_ieta_HE, emean_vs_ieta_HEM0, emean_vs_ieta_HEM3, emean_vs_ieta_HEP17, emean_vs_ieta_HEP17M0, emean_vs_ieta_HEP17M3, emean_vs_ieta_HF, emean_vs_ieta_HO, PVValHelper::eta, MonitorElement::Fill(), HcalObjRepresent::Fill(), fillRecHitsTmp(), edm::EventSetup::get(), edm::Event::getByToken(), hcalHBSevLvlVec, hcalHESevLvlVec, hcalHFSevLvlVec, hcalHOSevLvlVec, hep17_, mps_fire::i, imc, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meLog10Chi2profileHB, meLog10Chi2profileHE, meRecHitsEnergyHB, meRecHitsEnergyHBM0, meRecHitsEnergyHBM3, meRecHitsEnergyHE, meRecHitsEnergyHEM0, meRecHitsEnergyHEM3, meRecHitsEnergyHEP17, meRecHitsEnergyHEP17M0, meRecHitsEnergyHEP17M3, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitsEnergyM2vM0HB, meRecHitsEnergyM2vM0HE, meRecHitsEnergyM3vM0HB, meRecHitsEnergyM3vM0HE, meRecHitsEnergyM3vM2HB, meRecHitsEnergyM3vM2HE, meRecHitsM2Chi2HB, meRecHitsM2Chi2HE, meTE_HB, meTE_HE, meTE_HF, meTE_High_HB, meTE_High_HO, meTE_HO, meTE_Low_HB, meTE_Low_HE, meTE_Low_HF, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, meTimeHB, meTimeHE, meTimeHF, meTimeHO, nevtot, Nhb, Nhe, Nhf, Nho, nrechits_vs_iphi_HBM, nrechits_vs_iphi_HBP, nrechits_vs_iphi_HEM, nrechits_vs_iphi_HEP, nrechits_vs_iphi_HFM, nrechits_vs_iphi_HFP, nrechits_vs_iphi_HOM, nrechits_vs_iphi_HOP, occupancy_map_HB, occupancy_map_HE, occupancy_map_HF, occupancy_map_HO, phi, phi12(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), rpcPointValidation_cfi::recHit, RecHit_Aux_StatusWord_HB, RecHit_Aux_StatusWord_HE, RecHit_Aux_StatusWord_HF, RecHit_Aux_StatusWord_HO, RecHit_StatusWord_HB, RecHit_StatusWord_HE, RecHit_StatusWord_HF, RecHit_StatusWord_HO, RecHit_StatusWordCorr_HB, RecHit_StatusWordCorr_HE, sevLvl_HB, sevLvl_HE, sevLvl_HF, sevLvl_HO, subdet_, lumiQTWidget::t, theHcalChStatus, theHcalSevLvlComputer, theHcalTopology, tok_EB_, and tok_EE_.

                                                                            {

  using namespace edm;

  // cuts for each subdet_ector mimiking  "Scheme B"
  //  double cutHB = 0.9, cutHE = 1.4, cutHO = 1.1, cutHFL = 1.2, cutHFS = 1.8; 

  // energy in HCAL
  double eHcal        = 0.;
  // Total numbet of RecHits in HCAL, in the cone, above 1 GeV theshold
  int nrechits       = 0;
  int nrechitsThresh = 0;

  // energy in ECAL
  double eEcal       = 0.;
  double eEcalB      = 0.;
  double eEcalE      = 0.;
  double eEcalCone   = 0.;

  // HCAL energy around MC eta-phi at all depths;
  double partR = 0.3;

  // Single particle samples: actual eta-phi position of cluster around
  // hottest cell
  double etaHot  = 99999.; 
  double phiHot  = 99999.; 

  //   previously was:  c.get<IdealGeometryRecord>().get (geometry);
  c.get<CaloGeometryRecord>().get (geometry);

  // HCAL Topology **************************************************
  edm::ESHandle<HcalTopology> topo;
  c.get<HcalRecNumberingRecord>().get(topo);
  theHcalTopology = topo.product();

  // HCAL channel status map ****************************************
  edm::ESHandle<HcalChannelQuality> hcalChStatus;
  c.get<HcalChannelQualityRcd>().get( "withTopo", hcalChStatus );
  theHcalChStatus = hcalChStatus.product();

  // Assignment of severity levels **********************************
  edm::ESHandle<HcalSeverityLevelComputer> hcalSevLvlComputerHndl;
  c.get<HcalSeverityLevelComputerRcd>().get(hcalSevLvlComputerHndl);
  theHcalSevLvlComputer = hcalSevLvlComputerHndl.product(); 

  // Fill working vectors of HCAL RecHits quantities (all of these are drawn)
  fillRecHitsTmp(subdet_, ev); 

  // HB   
  if( subdet_ ==5 || subdet_ == 1 ){ 
     for(unsigned int iv=0; iv<hcalHBSevLvlVec.size(); iv++){
        sevLvl_HB->Fill(hcalHBSevLvlVec[iv]);
     }    
  }
  // HE   
  if( subdet_ ==5 || subdet_ == 2 ){
     for(unsigned int iv=0; iv<hcalHESevLvlVec.size(); iv++){
        sevLvl_HE->Fill(hcalHESevLvlVec[iv]);
     }
  }
  // HO 
  if( subdet_ ==5 || subdet_ == 3 ){
     for(unsigned int iv=0; iv<hcalHOSevLvlVec.size(); iv++){
        sevLvl_HO->Fill(hcalHOSevLvlVec[iv]);
     }
  }
  // HF 
  if( subdet_ ==5 || subdet_ == 4 ){
     for(unsigned int iv=0; iv<hcalHFSevLvlVec.size(); iv++){
        sevLvl_HF->Fill(hcalHFSevLvlVec[iv]);
     }
  } 

  //===========================================================================
  // IN ALL other CASES : ieta-iphi maps 
  //===========================================================================

  // ECAL 
  if(ecalselector_ == "yes" && (subdet_ == 1 || subdet_ == 2 || subdet_ == 5)) {

    Handle<EBRecHitCollection> rhitEB;
    if ( ev.getByToken(tok_EB_, rhitEB) ) {
      
      for ( const auto & recHit : *(rhitEB.product()) ) {
    
    double en  = recHit.energy();
    eEcal  += en;
    eEcalB += en;
      }      
    }
  
    
    Handle<EERecHitCollection> rhitEE;
    if ( ev.getByToken(tok_EE_, rhitEE) ) {
    
      for ( const auto & recHit : *(rhitEE.product()) ) {
    
    double en   = recHit.energy();
    eEcal  += en;
    eEcalE += en;
    
      }
    }
  }     // end of ECAL selection 

  // Counting, including ZS items
  // Filling HCAL maps  ----------------------------------------------------
  //   double maxE = -99999.;
  
  // element 0: any depth. element 1,2,..: depth 1,2
  std::vector<int> nhb_v(maxDepthHB_+1,0);
  std::vector<int> nhe_v(maxDepthHE_+1,0);
  std::vector<int> nho_v(maxDepthHO_+1,0);
  std::vector<int> nhf_v(maxDepthHF_+1,0);

  for (unsigned int i = 0; i < cen.size(); i++) {
    
    int sub       = csub[i];
    int depth     = cdepth[i];
    int ieta      = cieta[i]; 
    int iphi      = ciphi[i]; 
    double en     = cen[i]; 
    double enM0   = cenM0[i];
    double enM3   = cenM3[i];
    //     double eta    = ceta[i]; 
    //     double phi    = cphi[i]; 
    uint32_t stwd = cstwd[i];
    uint32_t auxstwd = cauxstwd[i];
    //    double z   = cz[i];

    //This will be true if hep17 == "yes" and the rechit is in the hep17 wedge
    bool isHEP17 = (sub == 2) && (iphi >= 63) && (iphi <= 66) && (ieta > 0) && (hep17_);

    //Make sure that an invalid depth won't cause an error. We should probably report the problem as well.
    if( depth < 1 ) continue;
    if( sub == 1 && depth > maxDepthHB_ ) continue;
    if( sub == 2 && depth > maxDepthHE_ ) continue;
    if( sub == 3 && depth > maxDepthHO_ ) continue;
    if( sub == 4 && depth > maxDepthHF_ ) continue;

    if( sub ==1 ){ nhb_v[depth]++; nhb_v[0]++;} // element 0: any depth, element 1,2,..: depth 1,2,...
    if( sub ==2 ){ nhe_v[depth]++; nhe_v[0]++;} //
    if( sub ==3 ){ nho_v[depth]++; nho_v[0]++;} //
    if( sub ==4 ){ nhf_v[depth]++; nhf_v[0]++;} //

    if( subdet_ == 6) {                                    // ZS specific
    }

    if( subdet_ != 6) {  
      int ieta2 = ieta;
      int depth2 = depth;
      if(sub == 4){
    if (ieta2 < 0) ieta2--;
        else ieta2++;
      }
      if(sub == 3) emap_HO->Fill(double(ieta2),double(iphi),en); //HO   
      else emap[depth2-1]->Fill(double(ieta2),double(iphi),en); // HB+HE+HF

      // to distinguish HE and HF
      if( depth == 1 || depth == 2 ) {
        int ieta1 =  ieta;
    if(sub == 4) { 
      if (ieta1 < 0) ieta1--;
          else  ieta1++;   
    }
      }

      if ( sub == 1){
     emean_vs_ieta_HB[depth-1]->Fill(double(ieta), en);
     emean_vs_ieta_HBM0[depth-1]->Fill(double(ieta), enM0);
     emean_vs_ieta_HBM3[depth-1]->Fill(double(ieta), enM3);
     occupancy_map_HB[depth-1]->Fill(double(ieta),double(iphi));
     if (ieta>0) nrechits_vs_iphi_HBP[depth-1]->Fill(double(iphi));
     else        nrechits_vs_iphi_HBM[depth-1]->Fill(double(iphi));
      }
      if ( sub == 2){
     if(!isHEP17){
        emean_vs_ieta_HE[depth-1]->Fill(double(ieta), en);
        emean_vs_ieta_HEM0[depth-1]->Fill(double(ieta), enM0);
        emean_vs_ieta_HEM3[depth-1]->Fill(double(ieta), enM3);
         }else{
        emean_vs_ieta_HEP17[depth-1]->Fill(double(ieta), en);
        emean_vs_ieta_HEP17M0[depth-1]->Fill(double(ieta), enM0);
        emean_vs_ieta_HEP17M3[depth-1]->Fill(double(ieta), enM3);
         }
     occupancy_map_HE[depth-1]->Fill(double(ieta),double(iphi));
     if (ieta>0) nrechits_vs_iphi_HEP[depth-1]->Fill(double(iphi));
     else        nrechits_vs_iphi_HEM[depth-1]->Fill(double(iphi));
      }
      if ( sub == 3){
     emean_vs_ieta_HO->Fill(double(ieta), en);
     occupancy_map_HO->Fill(double(ieta),double(iphi));
     if (ieta>0) nrechits_vs_iphi_HOP->Fill(double(iphi));
     else        nrechits_vs_iphi_HOM->Fill(double(iphi));
      }
      if ( sub == 4){
     emean_vs_ieta_HF[depth-1]->Fill(double(ieta), en);
     occupancy_map_HF[depth-1]->Fill(double(ieta),double(iphi));
     if (ieta>0) nrechits_vs_iphi_HFP[depth-1]->Fill(double(iphi));
     else        nrechits_vs_iphi_HFM[depth-1]->Fill(double(iphi));
      }
    }


    
    //32-bit status word  
    uint32_t statadd;
    unsigned int isw67 = 0;

    //Statusword correlation
    unsigned int sw27 = 27;
    unsigned int sw13 = 13;

    uint32_t statadd27 = 0x1<<sw27;
    uint32_t statadd13 = 0x1<<sw13;

    float status27 = 0;
    float status13 = 0;

    if(stwd & statadd27) status27 = 1;
    if(stwd & statadd13) status13 = 1;

    if        (sub == 1){
      RecHit_StatusWordCorr_HB->Fill(status13, status27);
    } else if (sub == 2){
      RecHit_StatusWordCorr_HE->Fill(status13, status27);
    }


    for (unsigned int isw = 0; isw < 32; isw++){
      statadd = 0x1<<(isw);
      if (stwd & statadd){
    if      (sub == 1) RecHit_StatusWord_HB->Fill(isw);
    else if (sub == 2) RecHit_StatusWord_HE->Fill(isw);
    else if (sub == 3) RecHit_StatusWord_HO->Fill(isw);
    else if (sub == 4){
      RecHit_StatusWord_HF->Fill(isw);
      if (isw == 6) isw67 += 1;
      if (isw == 7) isw67 += 2;
    }
      }
    }

    for (unsigned int isw =0; isw < 32; isw++){
      statadd = 0x1<<(isw);
      if( auxstwd & statadd ){
        if      (sub == 1) RecHit_Aux_StatusWord_HB->Fill(isw);
        else if (sub == 2) RecHit_Aux_StatusWord_HE->Fill(isw);
        else if (sub == 3) RecHit_Aux_StatusWord_HO->Fill(isw);
        else if (sub == 4) RecHit_Aux_StatusWord_HF->Fill(isw);
      }

    }

  } 

    for(int depth = 0; depth <= maxDepthHB_; depth++) Nhb[depth]->Fill(double(nhb_v[depth]));
    for(int depth = 0; depth <= maxDepthHE_; depth++) Nhe[depth]->Fill(double(nhe_v[depth]));
    for(int depth = 0; depth <= maxDepthHO_; depth++) Nho[depth]->Fill(double(nho_v[depth]));
    for(int depth = 0; depth <= maxDepthHF_; depth++) Nhf[depth]->Fill(double(nhf_v[depth]));

  //===========================================================================
  // SUBSYSTEMS,  
  //===========================================================================
  
  if ((subdet_ != 6) && (subdet_ != 0)) {

    double clusEta = 999.;
    double clusPhi = 999.; 
    double clusEn  = 0.;
    
    double HcalCone    = 0.;

    int ietaMax   =  9999;
    //     double enMax1 = -9999.;
    //     double enMax2 = -9999.;
    //     double enMax3 = -9999.;
    //     double enMax4 = -9999.;
    //     double enMax  = -9999.;
    //     double etaMax =  9999.;

    //   CYCLE over cells ====================================================

    for (unsigned int i = 0; i < cen.size(); i++) {
      int sub    = csub[i];
      double eta = ceta[i]; 
      double phi = cphi[i]; 
      double ieta = cieta[i]; 
      double iphi = ciphi[i]; 
      double en  = cen[i]; 
      double enM0  = cenM0[i]; 
      double enM3  = cenM3[i]; 
      double chi2  = cchi2[i];
      double chi2_log10=9.99; // initial value - stay with this value if chi2<0.
      if (chi2>0.) chi2_log10=log10(chi2);
      double t   = ctime[i];
      double depth = cdepth[i];

      bool isHEP17 = (sub == 2) && (iphi >= 63) && (iphi <= 66) && (ieta > 0) && (hep17_);

//       int   ieta = cieta[i];

      double rhot = dR(etaHot, phiHot, eta, phi); 
      if(rhot < partR && en > 1.) { 
    clusEta = (clusEta * clusEn + eta * en)/(clusEn + en);
        clusPhi = phi12(clusPhi, clusEn, phi, en); 
        clusEn += en;
      }

      nrechits++;       
      eHcal += en;

      if(en > 1. ) nrechitsThresh++;
      
      //The energy and overall timing histos are drawn while
      //the ones split by depth are not
      if(sub == 1 && (subdet_ == 1 || subdet_ == 5)) {  
    meTimeHB->Fill(t);
    meRecHitsEnergyHB->Fill(en);
    meRecHitsEnergyHBM0->Fill(enM0);
    meRecHitsEnergyHBM3->Fill(enM3);

        meRecHitsEnergyM2vM0HB->Fill(enM0,en);
        meRecHitsEnergyM3vM0HB->Fill(enM0,enM3);
        meRecHitsEnergyM3vM2HB->Fill(en,enM3);

        meRecHitsM2Chi2HB->Fill(chi2_log10);
        meLog10Chi2profileHB->Fill(en,chi2_log10);
    
    meTE_Low_HB->Fill( en, t);
    meTE_HB->Fill( en, t);
    meTE_High_HB->Fill( en, t);
    meTEprofileHB_Low->Fill(en, t);
    meTEprofileHB->Fill(en, t);
    meTEprofileHB_High->Fill(en, t);
      }     
      if(sub == 2 && (subdet_ == 2 || subdet_ == 5)) {  
    meTimeHE->Fill(t);
        if(!isHEP17){
       meRecHitsEnergyHE->Fill(en);
       meRecHitsEnergyHEM0->Fill(enM0);
       meRecHitsEnergyHEM3->Fill(enM3);
        }else{
       meRecHitsEnergyHEP17[0]->Fill(en);
       meRecHitsEnergyHEP17M0[0]->Fill(enM0);
       meRecHitsEnergyHEP17M3[0]->Fill(enM3);
       meRecHitsEnergyHEP17[depth]->Fill(en);
       meRecHitsEnergyHEP17M0[depth]->Fill(enM0);
       meRecHitsEnergyHEP17M3[depth]->Fill(enM3);
        }

        meRecHitsEnergyM2vM0HE->Fill(enM0,en);
        meRecHitsEnergyM3vM0HE->Fill(enM0,enM3);
        meRecHitsEnergyM3vM2HE->Fill(en,enM3);

        meRecHitsM2Chi2HE->Fill(chi2_log10);
        meLog10Chi2profileHE->Fill(en,chi2_log10);  

    meTE_Low_HE->Fill( en, t);
    meTE_HE->Fill( en, t);
    meTEprofileHE_Low->Fill(en, t);
    meTEprofileHE->Fill(en, t);
      }
      if(sub == 4 && (subdet_ == 4 || subdet_ == 5)) {  
    meTimeHF->Fill(t);
    meRecHitsEnergyHF->Fill(en);      

    meTE_Low_HF->Fill(en, t);
    meTE_HF->Fill(en, t);
    meTEprofileHF_Low->Fill(en, t);
    meTEprofileHF->Fill(en, t);

      }
      if(sub == 3 && (subdet_ == 3 || subdet_ == 5)) {  
    meTimeHO->Fill(t);
    meRecHitsEnergyHO->Fill(en);

    meTE_HO->Fill( en, t);
    meTE_High_HO->Fill( en, t);
    meTEprofileHO->Fill(en, t);
    meTEprofileHO_High->Fill(en, t);
      }
    }

    if(imc != 0) {
      //Cone by depth are not drawn, the others are used for pion scan
      meEnConeEtaProfile       ->Fill(double(ietaMax),  HcalCone);   // 
      meEnConeEtaProfile_E     ->Fill(double(ietaMax), eEcalCone);   
      meEnConeEtaProfile_EH    ->Fill(double(ietaMax),  HcalCone+eEcalCone); 
    }

    // Single particle samples ONLY !  ======================================
    // Fill up some histos for "integrated" subsustems. 
    // These are not drawn
  }

  nevtot++;
}

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 120 of file HcalRecHitsAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), egammaForCoreTracking_cff::depth, emap, emap_HO, emean_vs_ieta_HB, emean_vs_ieta_HBM0, emean_vs_ieta_HBM3, emean_vs_ieta_HE, emean_vs_ieta_HEM0, emean_vs_ieta_HEM3, emean_vs_ieta_HEP17, emean_vs_ieta_HEP17M0, emean_vs_ieta_HEP17M3, emean_vs_ieta_HF, emean_vs_ieta_HO, hep17_, trackerHits::histo, ieta_bins_, ieta_max_, ieta_min_, imc, createfilelist::int, iphi_bins_, iphi_max_, iphi_min_, maxDepthAll_, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meLog10Chi2profileHB, meLog10Chi2profileHE, meRecHitsEnergyHB, meRecHitsEnergyHBM0, meRecHitsEnergyHBM3, meRecHitsEnergyHE, meRecHitsEnergyHEM0, meRecHitsEnergyHEM3, meRecHitsEnergyHEP17, meRecHitsEnergyHEP17M0, meRecHitsEnergyHEP17M3, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitsEnergyM2vM0HB, meRecHitsEnergyM2vM0HE, meRecHitsEnergyM3vM0HB, meRecHitsEnergyM3vM0HE, meRecHitsEnergyM3vM2HB, meRecHitsEnergyM3vM2HE, meRecHitsM2Chi2HB, meRecHitsM2Chi2HE, meTE_HB, meTE_HE, meTE_HF, meTE_High_HB, meTE_High_HO, meTE_HO, meTE_Low_HB, meTE_Low_HE, meTE_Low_HF, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, meTimeHB, meTimeHE, meTimeHF, meTimeHO, nChannels_, Nhb, Nhe, Nhf, Nho, nrechits_vs_iphi_HBM, nrechits_vs_iphi_HBP, nrechits_vs_iphi_HEM, nrechits_vs_iphi_HEP, nrechits_vs_iphi_HFM, nrechits_vs_iphi_HFP, nrechits_vs_iphi_HOM, nrechits_vs_iphi_HOP, occupancy_map_HB, occupancy_map_HE, occupancy_map_HF, occupancy_map_HO, occupancy_vs_ieta_HB, occupancy_vs_ieta_HE, occupancy_vs_ieta_HF, occupancy_vs_ieta_HO, RecHit_Aux_StatusWord_HB, RecHit_Aux_StatusWord_HE, RecHit_Aux_StatusWord_HF, RecHit_Aux_StatusWord_HO, RecHit_StatusWord_HB, RecHit_StatusWord_HE, RecHit_StatusWord_HF, RecHit_StatusWord_HO, RecHit_StatusWordCorr_HB, RecHit_StatusWordCorr_HE, DQMStore::IBooker::setCurrentFolder(), sevLvl_HB, sevLvl_HE, sevLvl_HF, sevLvl_HO, subdet_, and topFolderName_.

{ 

    Char_t histo[200];

    ibooker.setCurrentFolder(topFolderName_);

    // General counters (drawn)

    //Produce both a total per subdetector, and number of rechits per subdetector depth
    // The bins are 1 unit wide, and the range is determined by the number of channels per subdetector

    for(int depth = 0; depth <= maxDepthHB_; depth++){
      if(depth == 0){ sprintf  (histo, "N_HB" );}
      else{           sprintf  (histo, "N_HB_depth%d",depth );}
      int NBins = (int) (nChannels_[1] * 1.1);
      Nhb.push_back( ibooker.book1D(histo, histo, NBins, 0., (float)NBins) );
    } 
    for(int depth = 0; depth <= maxDepthHE_; depth++){
      if(depth == 0){ sprintf  (histo, "N_HE" );}
      else{           sprintf  (histo, "N_HE_depth%d",depth );}
      int NBins = (int) (nChannels_[2] * 1.1);
      Nhe.push_back( ibooker.book1D(histo, histo, NBins,0., (float)NBins) );
    } 
    for(int depth = 0; depth <= maxDepthHO_; depth++){
      if(depth == 0){ sprintf  (histo, "N_HO" );}
      else{           sprintf  (histo, "N_HO_depth%d",depth );}
      int NBins = (int) (nChannels_[3] * 1.1);
      Nho.push_back( ibooker.book1D(histo, histo, NBins,0., (float)NBins) );
    } 
    for(int depth = 0; depth <= maxDepthHF_; depth++){
      if(depth == 0){ sprintf  (histo, "N_HF" );}
      else{           sprintf  (histo, "N_HF_depth%d",depth );}
      int NBins = (int) (nChannels_[4] * 1.1);
      Nhf.push_back( ibooker.book1D(histo, histo, NBins,0., (float)NBins) );
    } 

    // ZS
    if(subdet_ == 6) {

    }

    // ALL others, except ZS
    else {  
      for(int depth = 1; depth <= maxDepthAll_; depth++){
        sprintf  (histo, "emap_depth%d",depth );
        emap.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
      }
      sprintf(histo, "emap_HO");
      emap_HO = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_); 

      //The mean energy histos are drawn, but not the RMS or emean seq
      
      for (int depth = 1; depth <= maxDepthHB_; depth++) {
    sprintf  (histo, "emean_vs_ieta_HB%d",depth );
    emean_vs_ieta_HB.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );

    sprintf  (histo, "emean_vs_ieta_M0_HB%d",depth );
    emean_vs_ieta_HBM0.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );

    sprintf  (histo, "emean_vs_ieta_M3_HB%d",depth );
    emean_vs_ieta_HBM3.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
      }
      for (int depth = 1; depth <= maxDepthHE_; depth++) {
    sprintf  (histo, "emean_vs_ieta_HE%d",depth );
    emean_vs_ieta_HE.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );

    sprintf  (histo, "emean_vs_ieta_M0_HE%d",depth );
    emean_vs_ieta_HEM0.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );

    sprintf  (histo, "emean_vs_ieta_M3_HE%d",depth );
    emean_vs_ieta_HEM3.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
      }

      if(hep17_){
         for (int depth = 1; depth <= maxDepthHE_; depth++) {
            sprintf  (histo, "emean_vs_ieta_HEP17_depth%d",depth );
            emean_vs_ieta_HEP17.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );

        sprintf  (histo, "emean_vs_ieta_M0_HEP17_depth%d",depth );
        emean_vs_ieta_HEP17M0.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );

        sprintf  (histo, "emean_vs_ieta_M3_HEP17_depth%d",depth );
        emean_vs_ieta_HEP17M3.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
         }
      }

      for (int depth = 1; depth <= maxDepthHF_; depth++) {
    sprintf  (histo, "emean_vs_ieta_HF%d",depth );
    emean_vs_ieta_HF.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
      }
      sprintf  (histo, "emean_vs_ieta_HO" );
      emean_vs_ieta_HO = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " " );

      //The only occupancy histos drawn are occupancy vs. ieta
      //but the maps are needed because this is where the latter are filled from

      for (int depth = 1; depth <= maxDepthHB_; depth++) {
         sprintf  (histo, "occupancy_map_HB%d",depth );
         occupancy_map_HB.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
      }

      for (int depth = 1; depth <= maxDepthHE_; depth++) {
         sprintf  (histo, "occupancy_map_HE%d",depth );
         occupancy_map_HE.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
      }

      sprintf  (histo, "occupancy_map_HO" );
      occupancy_map_HO = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);      

      for (int depth = 1; depth <= maxDepthHF_; depth++) {
         sprintf  (histo, "occupancy_map_HF%d",depth );
         occupancy_map_HF.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
      }
      
      // nrechits vs iphi
      for (int depth = 1; depth <= maxDepthHB_; depth++) {
    
     sprintf  (histo, "occupancy_vs_ieta_HB%d",depth );
     occupancy_vs_ieta_HB.push_back( ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_) );
         sprintf  (histo, "nrechits_vs_iphi_HBP_d%d",depth );
         nrechits_vs_iphi_HBP.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
         sprintf  (histo, "nrechits_vs_iphi_HBM_d%d",depth );
         nrechits_vs_iphi_HBM.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
      }

      for (int depth = 1; depth <= maxDepthHE_; depth++) {
         sprintf  (histo, "occupancy_vs_ieta_HE%d",depth );
         occupancy_vs_ieta_HE.push_back( ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_) );
         sprintf  (histo, "nrechits_vs_iphi_HEP_d%d",depth );
         nrechits_vs_iphi_HEP.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
         sprintf  (histo, "nrechits_vs_iphi_HEM_d%d",depth );
         nrechits_vs_iphi_HEM.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
      }

      sprintf  (histo, "occupancy_vs_ieta_HO" );
      occupancy_vs_ieta_HO = ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_);
      sprintf  (histo, "nrechits_vs_iphi_HOP" );
      nrechits_vs_iphi_HOP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "nrechits_vs_iphi_HOM" );
      nrechits_vs_iphi_HOM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);

      for (int depth = 1; depth <= maxDepthHF_; depth++) {
         sprintf  (histo, "occupancy_vs_ieta_HF%d",depth );
         occupancy_vs_ieta_HF.push_back( ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_) );
         sprintf  (histo, "nrechits_vs_iphi_HFP_d%d",depth );
         nrechits_vs_iphi_HFP.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
         sprintf  (histo, "nrechits_vs_iphi_HFM_d%d",depth );
         nrechits_vs_iphi_HFM.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
      }

      //All status word histos except HF67 are drawn
      sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HB" ) ;
      RecHit_StatusWord_HB = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
      
      sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HE" ) ;
      RecHit_StatusWord_HE = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 

      sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HF" ) ;
      RecHit_StatusWord_HF = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 

      sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HO" ) ;
      RecHit_StatusWord_HO = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 

      //Aux status word histos
      sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HB" ) ;
      RecHit_Aux_StatusWord_HB = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
      
      sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HE" ) ;
      RecHit_Aux_StatusWord_HE = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 

      sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HF" ) ;
      RecHit_Aux_StatusWord_HF = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 

      sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HO" ) ;
      RecHit_Aux_StatusWord_HO = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 

    }  // end-of (subdet_ =! 6)

      //Status word correlations
      sprintf (histo, "HcalRecHitTask_RecHit_StatusWordCorr_HB");
      RecHit_StatusWordCorr_HB = ibooker.book2D(histo, histo, 2, -0.5, 1.5, 2, -0.5, 1.5);

      sprintf (histo, "HcalRecHitTask_RecHit_StatusWordCorr_HE");
      RecHit_StatusWordCorr_HE = ibooker.book2D(histo, histo, 2, -0.5, 1.5, 2, -0.5, 1.5);


    //======================= Now various cases one by one ===================

    //Histograms drawn for single pion scan
    if(subdet_ != 0 && imc != 0) { // just not for noise  
      sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths");
      meEnConeEtaProfile = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -100., 2000., " ");  
      
      sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
      meEnConeEtaProfile_E = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -100., 2000., " ");  
      
      sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
      meEnConeEtaProfile_EH = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -100., 2000., " ");  
    }

    // ************** HB **********************************
    if (subdet_ == 1 || subdet_ == 5 ){

      //Only severity level, energy of rechits and overall HB timing histos are drawn  

      sprintf(histo, "HcalRecHitTask_severityLevel_HB");
      sevLvl_HB = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HB" ) ;
      meRecHitsEnergyHB = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HB" ) ;
      meRecHitsEnergyHBM0 = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HB" ) ;
      meRecHitsEnergyHBM3 = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M2vM0_HB" ) ;
      meRecHitsEnergyM2vM0HB = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM0_HB" ) ;
      meRecHitsEnergyM3vM0HB = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM2_HB" ) ;
      meRecHitsEnergyM3vM2HB = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
      
      sprintf (histo, "HcalRecHitTask_M2Log10Chi2_of_rechits_HB" ) ;
      meRecHitsM2Chi2HB = ibooker.book1D(histo, histo, 120 , -2. , 10.); 
      
      sprintf (histo, "HcalRecHitTask_timing_HB" ) ;
      meTimeHB = ibooker.book1D(histo, histo, 70, -48., 92.); 

      //High, medium and low histograms to reduce RAM usage
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HB" ) ;
      meTE_Low_HB = ibooker.book2D(histo, histo, 50, -5., 45.,  70, -48., 92.);

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_HB" ) ;
      meTE_HB = ibooker.book2D(histo, histo, 150, -5., 295.,  70, -48., 92.);

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_High_HB" ) ;
      meTE_High_HB = ibooker.book2D(histo, histo, 150, -5., 2995.,  70, -48., 92.);
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HB" ) ;
      meTEprofileHB_Low = ibooker.bookProfile(histo, histo, 50, -5., 45., -48., 92., " "); 

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HB" ) ;
      meTEprofileHB = ibooker.bookProfile(histo, histo, 150, -5., 295., -48., 92., " "); 

      sprintf (histo, "HcalRecHitTask_Log10Chi2_vs_energy_profile_HB" ) ;
      meLog10Chi2profileHB = ibooker.bookProfile(histo, histo, 150, -5., 295., -2., 9.9, " "); 

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HB" ) ;
      meTEprofileHB_High = ibooker.bookProfile(histo, histo, 150, -5., 2995., -48., 92., " "); 

    }
    
    // ********************** HE ************************************
    if ( subdet_ == 2 || subdet_ == 5 ){


      //Only severity level, energy of rechits and overall HB timing histos are drawn  
      sprintf(histo, "HcalRecHitTask_severityLevel_HE");
      sevLvl_HE = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HE" ) ;
      meRecHitsEnergyHE = ibooker.book1D(histo, histo, 2010, -10., 2000.);
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HE" ) ;
      meRecHitsEnergyHEM0 = ibooker.book1D(histo, histo, 2010, -10., 2000.);
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HE" ) ;
      meRecHitsEnergyHEM3 = ibooker.book1D(histo, histo, 2010, -10., 2000.);
      
      if(hep17_){
         sprintf (histo, "HcalRecHitTask_energy_of_rechits_HEP17" ) ;
         meRecHitsEnergyHEP17.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
      
         sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HEP17" ) ;
         meRecHitsEnergyHEP17M0.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
      
         sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HEP17" ) ;
         meRecHitsEnergyHEP17M3.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
         for (int depth = 1; depth <= maxDepthHE_; depth++) {
            sprintf (histo, "HcalRecHitTask_energy_of_rechits_HEP17_depth%d",depth ) ;
            meRecHitsEnergyHEP17.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
      
            sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HEP17_depth%d",depth ) ;
            meRecHitsEnergyHEP17M0.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
      
            sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HEP17_depth%d",depth ) ;
            meRecHitsEnergyHEP17M3.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
         }
      }

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M2vM0_HE" ) ;
      meRecHitsEnergyM2vM0HE = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM0_HE" ) ;
      meRecHitsEnergyM3vM0HE = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM2_HE" ) ;
      meRecHitsEnergyM3vM2HE = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
      
      sprintf (histo, "HcalRecHitTask_M2Log10Chi2_of_rechits_HE" ) ;
      meRecHitsM2Chi2HE = ibooker.book1D(histo, histo, 120 , -2. , 10.); 
      
      sprintf (histo, "HcalRecHitTask_timing_HE" ) ;
      meTimeHE = ibooker.book1D(histo, histo, 70, -48., 92.); 
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HE" ) ;
      meTE_Low_HE = ibooker.book2D(histo, histo, 80, -5., 75.,  70, -48., 92.);

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_HE" ) ;
      meTE_HE = ibooker.book2D(histo, histo, 200, -5., 395.,  70, -48., 92.);
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HE" ) ;
      meTEprofileHE_Low = ibooker.bookProfile(histo, histo, 80, -5., 75., -48., 92., " "); 

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HE" ) ;
      meTEprofileHE = ibooker.bookProfile(histo, histo, 200, -5., 395., -48., 92., " "); 
      
      sprintf (histo, "HcalRecHitTask_Log10Chi2_vs_energy_profile_HE" ) ;
      meLog10Chi2profileHE = ibooker.bookProfile(histo, histo, 200, -5., 395., -2., 9.9, " ");       

    }

    // ************** HO ****************************************
    if ( subdet_ == 3 || subdet_ == 5  ){
      
      //Only severity level, energy of rechits and overall HB timing histos are drawn  

      sprintf(histo, "HcalRecHitTask_severityLevel_HO");
      sevLvl_HO = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HO" ) ;
      meRecHitsEnergyHO = ibooker.book1D(histo, histo, 2010 , -10. , 2000.);
      
      sprintf (histo, "HcalRecHitTask_timing_HO" ) ;
      meTimeHO = ibooker.book1D(histo, histo, 70, -48., 92.); 
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_HO" ) ;
      meTE_HO= ibooker.book2D(histo, histo, 60, -5., 55., 70, -48., 92.);

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_High_HO" ) ;
      meTE_High_HO= ibooker.book2D(histo, histo, 100, -5., 995., 70, -48., 92.);
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HO" ) ;
      meTEprofileHO = ibooker.bookProfile(histo, histo, 60, -5., 55., -48., 92., " "); 

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HO" ) ;
      meTEprofileHO_High = ibooker.bookProfile(histo, histo, 100, -5., 995., -48., 92., " "); 
      
    }   
  
    // ********************** HF ************************************
    if ( subdet_ == 4 || subdet_ == 5 ){

      //Only severity level, energy of rechits and overall HB timing histos are drawn  
      
      sprintf(histo, "HcalRecHitTask_severityLevel_HF");
      sevLvl_HF = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HF" ) ;
      meRecHitsEnergyHF = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 

      sprintf (histo, "HcalRecHitTask_timing_HF" ) ;
      meTimeHF = ibooker.book1D(histo, histo, 70, -48., 92.); 
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HF" ) ;
      meTE_Low_HF = ibooker.book2D(histo, histo, 100, -5., 195., 70, -48., 92.);

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_HF" ) ;
      meTE_HF = ibooker.book2D(histo, histo, 200, -5., 995., 70, -48., 92.);
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HF" ) ;
      meTEprofileHF_Low = ibooker.bookProfile(histo, histo, 100, -5., 195., -48., 92., " "); 

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HF" ) ;
      meTEprofileHF = ibooker.bookProfile(histo, histo, 200, -5., 995., -48., 92., " "); 

    }

}

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

private

Definition at line 1113 of file HcalRecHitsAnalyzer.cc.

References PI, and tmp.

                                                              {
  // 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;

}

void HcalRecHitsAnalyzer::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  c  )

override

Definition at line 59 of file HcalRecHitsAnalyzer.cc.

References geometry, edm::EventSetup::get(), HcalDDDRecConstants::getEtaRange(), HcalGeometry::getHxSize(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNPhi(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, ieta_bins_, ieta_max_, ieta_min_, hcalTTPDigis_cfi::iEtaMax, createfilelist::int, iphi_bins_, iphi_max_, iphi_min_, maxDepthAll_, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, nChannels_, and edm::ESHandle< T >::product().

                                                                                 {
  
    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);

    edm::ESHandle<CaloGeometry> geometry;

    es.get<CaloGeometryRecord > ().get(geometry);

    const CaloGeometry* geo = geometry.product();
    const HcalGeometry* gHB = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel));
    const HcalGeometry* gHE = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap));
    const HcalGeometry* gHO = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalOuter));
    const HcalGeometry* gHF = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalForward));

    nChannels_[1] = gHB->getHxSize(1); 
    nChannels_[2] = gHE->getHxSize(2); 
    nChannels_[3] = gHO->getHxSize(3); 
    nChannels_[4] = gHF->getHxSize(4); 

    nChannels_[0] = nChannels_[1] + nChannels_[2] + nChannels_[3] + nChannels_[4];

    //std::cout << "Channels HB:" << nChannels_[1] << " HE:" << nChannels_[2] << " HO:" << nChannels_[3] << " HF:" << nChannels_[4] << std::endl;

    //We hardcode the HF depths because in the dual readout configuration, rechits are not defined for depths 3&4
    maxDepthHF_ = (maxDepthHF_ > 2 ? 2 : maxDepthHF_); //We retain the dynamic possibility that HF might have 0 or 1 depths

    maxDepthAll_ = ( maxDepthHB_ + maxDepthHO_ > maxDepthHE_ ? maxDepthHB_ + maxDepthHO_ : maxDepthHE_ );
    maxDepthAll_ = ( maxDepthAll_ > maxDepthHF_ ? maxDepthAll_ : maxDepthHF_ );

    //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
    iphi_min_ = 0.5;
    iphi_max_ = NphiMax + 0.5;
    iphi_bins_ = (int) (iphi_max_ - iphi_min_);

    //Retain classic behavior, all plots have same ieta range.

    int iEtaMax = (hcons->getEtaRange(0).second > hcons->getEtaRange(1).second ? hcons->getEtaRange(0).second : hcons->getEtaRange(1).second);
    iEtaMax = (iEtaMax > hcons->getEtaRange(2).second ? iEtaMax : hcons->getEtaRange(2).second);
    iEtaMax = (iEtaMax > hcons->getEtaRange(3).second ? iEtaMax : hcons->getEtaRange(3).second);

    //Give an empty bin around the subdet ieta range to make it clear that all ieta rings have been included
    ieta_min_ = -iEtaMax - 1.5;
    ieta_max_ = iEtaMax + 1.5;
    ieta_bins_ = (int) (ieta_max_ - ieta_min_);

  }

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

private

Definition at line 1087 of file HcalRecHitsAnalyzer.cc.

References PI, mathSSE::sqrt(), and tmp.

Referenced by analyze().

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

void HcalRecHitsAnalyzer::fillRecHitsTmp ( int  subdet_, edm::Event const &  ev  )

privatevirtual

Definition at line 908 of file HcalRecHitsAnalyzer.cc.

References edm::SortedCollection< T, SORT >::begin(), cauxstwd, cchi2, cdepth, cen, cenM0, cenM3, ceta, vertices_cff::chi2, cieta, ciphi, cphi, cstwd, csub, ctime, cz, egammaForCoreTracking_cff::depth, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, edm::Event::getByToken(), ecaldqm::getGeometry(), HcalBarrel, HcalEndcap, HcalForward, hcalHBSevLvlVec, hcalHESevLvlVec, hcalHFSevLvlVec, hcalHOSevLvlVec, HcalOuter, hcalSevLvl(), iz, phi, interestingDetIdCollectionProducer_cfi::severityLevel, lumiQTWidget::t, tok_hbhe_, tok_hf_, and tok_ho_.

Referenced by analyze().

                                                                       {
  
  using namespace edm;
  
  
  // initialize data vectors
  csub.clear();
  cen.clear();
  cenM0.clear();
  cenM3.clear();
  cchi2.clear();
  ceta.clear();
  cphi.clear();
  ctime.clear();
  cieta.clear();
  ciphi.clear();
  cdepth.clear();
  cz.clear();
  cstwd.clear();
  cauxstwd.clear();
  hcalHBSevLvlVec.clear();
  hcalHESevLvlVec.clear();
  hcalHFSevLvlVec.clear();
  hcalHOSevLvlVec.clear(); 

  if( subdet_ == 1 || subdet_ == 2  || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    
    //HBHE
    edm::Handle<HBHERecHitCollection> hbhecoll;
    if ( ev.getByToken(tok_hbhe_, hbhecoll) ) {
    
      for (HBHERecHitCollection::const_iterator j=hbhecoll->begin(); j != hbhecoll->end(); j++) {
    HcalDetId cell(j->id());
    const HcalGeometry* cellGeometry = 
      dynamic_cast<const HcalGeometry*>(geometry->getSubdetectorGeometry(cell));
    double eta  = cellGeometry->getPosition(cell).eta () ;
    double phi  = cellGeometry->getPosition(cell).phi () ;
    double zc   = cellGeometry->getPosition(cell).z ();
    int sub     = cell.subdet();
    int depth   = cell.depth();
    int inteta  = cell.ieta();
    int intphi  = cell.iphi();
    double en   = j->energy();
    double enM0 = j->eraw();
    double enM3 = j->eaux();
    double chi2 = j->chi2();
    double t    = j->time();
    int stwd    = j->flags();
    int auxstwd = j->aux();
    
    int severityLevel = hcalSevLvl( (CaloRecHit*) &*j );
    if( cell.subdet()==HcalBarrel ){
      hcalHBSevLvlVec.push_back(severityLevel);
    }else if (cell.subdet()==HcalEndcap ){
      hcalHESevLvlVec.push_back(severityLevel);
    } 
    
    if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
      
      csub.push_back(sub);
      cen.push_back(en);
      cenM0.push_back(enM0);
      cenM3.push_back(enM3);
      cchi2.push_back(chi2);
      ceta.push_back(eta);
      cphi.push_back(phi);
      ctime.push_back(t);
      cieta.push_back(inteta);
      ciphi.push_back(intphi);
      cdepth.push_back(depth);
      cz.push_back(zc);
      cstwd.push_back(stwd);
      cauxstwd.push_back(auxstwd);
    }
      }
      
    }
  }

  if( subdet_ == 4 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {

    //HF
    edm::Handle<HFRecHitCollection> hfcoll;
    if ( ev.getByToken(tok_hf_, hfcoll) ) {
      
      for (HFRecHitCollection::const_iterator j = hfcoll->begin(); j != hfcoll->end(); j++) {
    HcalDetId cell(j->id());
    auto cellGeometry = (geometry->getSubdetectorGeometry(cell))->getGeometry (cell) ;
    double eta   = cellGeometry->getPosition().eta () ;
    double phi   = cellGeometry->getPosition().phi () ;
    double zc    = cellGeometry->getPosition().z ();
    int sub      = cell.subdet();
    int depth    = cell.depth();
    int inteta   = cell.ieta();
    int intphi   = cell.iphi();
    double en    = j->energy();
    double enM0  = 0.;
    double enM3  = 0.;
    double chi2  = 0.;
    double t     = j->time();
    int stwd     = j->flags();
    int auxstwd  = j->aux();
    
    int severityLevel = hcalSevLvl( (CaloRecHit*) &*j );
    if( cell.subdet()==HcalForward ){
      hcalHFSevLvlVec.push_back(severityLevel);
    } 
    
    if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
      
      csub.push_back(sub);
      cen.push_back(en);
      cenM0.push_back(enM0);
      cenM3.push_back(enM3);
      cchi2.push_back(chi2);
      ceta.push_back(eta);
      cphi.push_back(phi);
      ctime.push_back(t);
      cieta.push_back(inteta);
      ciphi.push_back(intphi);
      cdepth.push_back(depth);
      cz.push_back(zc);
      cstwd.push_back(stwd);
      cauxstwd.push_back(auxstwd);
    }
      }
    }
  }

  //HO
  if( subdet_ == 3 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
  
    edm::Handle<HORecHitCollection> hocoll;
    if ( ev.getByToken(tok_ho_, hocoll) ) {
      
      for (HORecHitCollection::const_iterator j = hocoll->begin(); j != hocoll->end(); j++) {
    HcalDetId cell(j->id());
    auto cellGeometry = (geometry->getSubdetectorGeometry(cell))->getGeometry (cell) ;
    double eta   = cellGeometry->getPosition().eta () ;
    double phi   = cellGeometry->getPosition().phi () ;
    double zc    = cellGeometry->getPosition().z ();
    int sub      = cell.subdet();
    int depth    = cell.depth();
    int inteta   = cell.ieta();
    int intphi   = cell.iphi();
    double t     = j->time();
    double en    = j->energy();
    double enM0  = 0.;
    double enM3  = 0.;
    double chi2 = 0.;
    int stwd     = j->flags();
    int auxstwd  = j->aux();
    
    int severityLevel = hcalSevLvl( (CaloRecHit*) &*j );
    if( cell.subdet()==HcalOuter ){
      hcalHOSevLvlVec.push_back(severityLevel);
    } 
    
    if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
      csub.push_back(sub);
      cen.push_back(en);
      cenM0.push_back(enM0);
      cenM3.push_back(enM3);
      cchi2.push_back(chi2);
      ceta.push_back(eta);
      cphi.push_back(phi);
      ctime.push_back(t);
      cieta.push_back(inteta);
      ciphi.push_back(intphi);
      cdepth.push_back(depth);
      cz.push_back(zc);
      cstwd.push_back(stwd);
      cauxstwd.push_back(auxstwd);
    }
      }
    }
  }
}

int HcalRecHitsAnalyzer::hcalSevLvl ( const CaloRecHit *  hit )

private

Definition at line 1128 of file HcalRecHitsAnalyzer.cc.

References DEFINE_FWK_MODULE, CaloRecHit::detid(), CaloRecHit::flags(), HcalTopology::getMergePositionFlag(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), HcalEndcap, HcalTopology::idFront(), interestingDetIdCollectionProducer_cfi::severityLevel, theHcalChStatus, theHcalSevLvlComputer, and theHcalTopology.

Referenced by fillRecHitsTmp().

                                                        {

   HcalDetId id = hit->detid();
   if (theHcalTopology->getMergePositionFlag() && id.subdet() == HcalEndcap) {
     id = theHcalTopology->idFront(id);
   }

   const uint32_t recHitFlag = hit->flags();
   const uint32_t dbStatusFlag = theHcalChStatus->getValues(id)->getValue();

   int severityLevel = theHcalSevLvlComputer->getSeverityLevel(id, recHitFlag, dbStatusFlag);

   return severityLevel;

} 

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

private

Definition at line 1097 of file HcalRecHitsAnalyzer.cc.

References PI, and tmp.

Referenced by analyze().

                                                                                  {
  // 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;

}

Member Data Documentation

std::vector<uint32_t> HcalRecHitsAnalyzer::cauxstwd

private

Definition at line 295 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cchi2

private

Definition at line 289 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::cdepth

private

Definition at line 285 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cen

private

Definition at line 286 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cenM0

private

Definition at line 287 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cenM3

private

Definition at line 288 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::ceta

private

Definition at line 290 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::cieta

private

Definition at line 283 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::ciphi

private

Definition at line 284 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cphi

private

Definition at line 291 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<uint32_t> HcalRecHitsAnalyzer::cstwd

private

Definition at line 294 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::csub

private

Definition at line 282 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::ctime

private

Definition at line 292 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cz

private

Definition at line 293 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

std::string HcalRecHitsAnalyzer::ecalselector_

private

Definition at line 77 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emap

private

Definition at line 139 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::emap_HO

private

Definition at line 140 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HB

private

Definition at line 142 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HBM0

private

Definition at line 143 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HBM3

private

Definition at line 144 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HE

private

Definition at line 145 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HEM0

private

Definition at line 146 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HEM3

private

Definition at line 147 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HEP17

private

Definition at line 148 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HEP17M0

private

Definition at line 149 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HEP17M3

private

Definition at line 150 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::emean_vs_ieta_HF

private

Definition at line 151 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::emean_vs_ieta_HO

private

Definition at line 152 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsAnalyzer::etype_

private

Definition at line 107 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

std::string HcalRecHitsAnalyzer::eventype_

private

Definition at line 78 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

bool HcalRecHitsAnalyzer::famos_

private

Definition at line 82 of file HcalRecHitsAnalyzer.h.

edm::ESHandle<CaloGeometry> HcalRecHitsAnalyzer::geometry

private

Definition at line 262 of file HcalRecHitsAnalyzer.h.

Referenced by dqmBeginRun(), and Vispa.Gui.ConnectableWidget.ConnectableWidget::leaveEvent().

std::vector<int> HcalRecHitsAnalyzer::hcalHBSevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHESevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHFSevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHOSevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::string HcalRecHitsAnalyzer::hcalselector_

private

Definition at line 76 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

const HcalDDDRecConstants* HcalRecHitsAnalyzer::hcons

private

Definition at line 84 of file HcalRecHitsAnalyzer.h.

Referenced by dqmBeginRun().

bool HcalRecHitsAnalyzer::hep17_

private

Definition at line 80 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and HcalRecHitsAnalyzer().

int HcalRecHitsAnalyzer::ieta_bins_

private

Definition at line 93 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::ieta_max_

private

Definition at line 94 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::ieta_min_

private

Definition at line 94 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::imc

private

Definition at line 109 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and HcalRecHitsAnalyzer().

int HcalRecHitsAnalyzer::iphi_bins_

private

Definition at line 90 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::iphi_max_

private

Definition at line 91 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::iphi_min_

private

Definition at line 91 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::iz

private

Definition at line 108 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

MonitorElement* HcalRecHitsAnalyzer::map_ecal

private

Definition at line 137 of file HcalRecHitsAnalyzer.h.

int HcalRecHitsAnalyzer::maxDepthAll_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::maxDepthHB_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::maxDepthHE_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::maxDepthHF_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::maxDepthHO_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

std::string HcalRecHitsAnalyzer::mc_

private

Definition at line 81 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meDeltaEta

private

Definition at line 175 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meDeltaPhi

private

Definition at line 174 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHB

private

Definition at line 249 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHE

private

Definition at line 250 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHF

private

Definition at line 252 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHO

private

Definition at line 251 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHB

private

Definition at line 246 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHE

private

Definition at line 247 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile

private

Definition at line 170 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_E

private

Definition at line 171 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_EH

private

Definition at line 172 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHB

private

Definition at line 255 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHE

private

Definition at line 256 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHB

private

Definition at line 212 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHE

private

Definition at line 219 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHB

private

Definition at line 259 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHE

private

Definition at line 260 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHB

private

Definition at line 184 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM0

private

Definition at line 185 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM3

private

Definition at line 186 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHE

private

Definition at line 192 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM0

private

Definition at line 193 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM3

private

Definition at line 194 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::meRecHitsEnergyHEP17

private

Definition at line 195 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::meRecHitsEnergyHEP17M0

private

Definition at line 196 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::meRecHitsEnergyHEP17M3

private

Definition at line 197 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHF

private

Definition at line 205 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHO

private

Definition at line 203 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HB

private

Definition at line 187 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HE

private

Definition at line 198 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HB

private

Definition at line 188 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HE

private

Definition at line 199 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HB

private

Definition at line 189 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HE

private

Definition at line 200 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HB

private

Definition at line 190 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HE

private

Definition at line 201 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHB

private

Definition at line 238 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHE

private

Definition at line 239 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHF

private

Definition at line 241 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFL

private

Definition at line 242 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFS

private

Definition at line 243 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHO

private

Definition at line 240 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHB

private

Definition at line 232 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHE

private

Definition at line 233 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHF

private

Definition at line 235 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHO

private

Definition at line 234 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meTE_HB

private

Definition at line 208 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HE

private

Definition at line 216 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HF

private

Definition at line 227 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HB

private

Definition at line 209 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HO

private

Definition at line 222 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HO

private

Definition at line 221 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HB

private

Definition at line 207 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HE

private

Definition at line 215 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HF

private

Definition at line 226 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB

private

Definition at line 211 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_High

private

Definition at line 213 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_Low

private

Definition at line 210 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE

private

Definition at line 218 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE_Low

private

Definition at line 217 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF

private

Definition at line 229 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF_Low

private

Definition at line 228 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO

private

Definition at line 223 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO_High

private

Definition at line 224 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHB

private

Definition at line 178 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHE

private

Definition at line 179 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHF

private

Definition at line 181 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHO

private

Definition at line 180 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsAnalyzer::nChannels_[5]

private

Definition at line 88 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::nevtot

private

Definition at line 298 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::Nhb

private

Definition at line 127 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::Nhe

private

Definition at line 128 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::Nhf

private

Definition at line 130 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::Nho

private

Definition at line 129 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::nrechits_vs_iphi_HBM

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::nrechits_vs_iphi_HBP

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::nrechits_vs_iphi_HEM

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::nrechits_vs_iphi_HEP

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::nrechits_vs_iphi_HFM

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::nrechits_vs_iphi_HFP

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * HcalRecHitsAnalyzer::nrechits_vs_iphi_HOM

private

Definition at line 167 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::nrechits_vs_iphi_HOP

private

Definition at line 167 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::occupancy_map_HB

private

Definition at line 154 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::occupancy_map_HE

private

Definition at line 155 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::occupancy_map_HF

private

Definition at line 156 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::occupancy_map_HO

private

Definition at line 157 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::occupancy_vs_ieta_HB

private

Definition at line 159 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::occupancy_vs_ieta_HE

private

Definition at line 160 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

std::vector<MonitorElement*> HcalRecHitsAnalyzer::occupancy_vs_ieta_HF

private

Definition at line 161 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::occupancy_vs_ieta_HO

private

Definition at line 162 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

std::string HcalRecHitsAnalyzer::outputFile_

private

Definition at line 75 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HB

private

Definition at line 276 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HE

private

Definition at line 277 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HF

private

Definition at line 278 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HO

private

Definition at line 279 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HB

private

Definition at line 265 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HE

private

Definition at line 266 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF

private

Definition at line 267 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF67

private

Definition at line 268 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HO

private

Definition at line 269 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HB

private

Definition at line 272 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HE

private

Definition at line 273 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HB

private

Definition at line 121 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HE

private

Definition at line 122 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HF

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HO

private

Definition at line 124 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string HcalRecHitsAnalyzer::sign_

private

Definition at line 79 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

int HcalRecHitsAnalyzer::subdet_

private

Definition at line 104 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), bookHistograms(), and HcalRecHitsAnalyzer().

const HcalChannelQuality* HcalRecHitsAnalyzer::theHcalChStatus

private

Definition at line 114 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

const HcalSeverityLevelComputer* HcalRecHitsAnalyzer::theHcalSevLvlComputer

private

Definition at line 116 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

const HcalTopology* HcalRecHitsAnalyzer::theHcalTopology

private

Definition at line 112 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

edm::EDGetTokenT<EBRecHitCollection> HcalRecHitsAnalyzer::tok_EB_

private

Definition at line 100 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<EERecHitCollection> HcalRecHitsAnalyzer::tok_EE_

private

Definition at line 101 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalRecHitsAnalyzer::tok_hbhe_

private

Definition at line 97 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HFRecHitCollection> HcalRecHitsAnalyzer::tok_hf_

private

Definition at line 99 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HORecHitCollection> HcalRecHitsAnalyzer::tok_ho_

private

Definition at line 98 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

std::string HcalRecHitsAnalyzer::topFolderName_

private

Definition at line 73 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().