HcalRecHitsValidation Class Reference

#include <HcalRecHitsValidation.h>

Inheritance diagram for HcalRecHitsValidation:

Public Member Functions
virtual void	analyze (edm::Event const &ev, edm::EventSetup const &c)
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
	HcalRecHitsValidation (edm::ParameterSet const &conf)
	~HcalRecHitsValidation ()
Public Member Functions inherited from thread_unsafe::DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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< int >	cdepth
std::vector< double >	cen
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
MonitorElement *	e_hb
MonitorElement *	e_he
MonitorElement *	e_hfl
MonitorElement *	e_hfs
MonitorElement *	e_ho
std::string	ecalselector_
MonitorElement *	emap_depth1
MonitorElement *	emap_depth2
MonitorElement *	emap_depth3
MonitorElement *	emap_depth4
double	emap_min [82][72][4][4]
MonitorElement *	emean_seqHB1
MonitorElement *	emean_seqHB2
MonitorElement *	emean_seqHE1
MonitorElement *	emean_seqHE2
MonitorElement *	emean_seqHE3
MonitorElement *	emean_seqHF1
MonitorElement *	emean_seqHF2
MonitorElement *	emean_seqHO
MonitorElement *	emean_vs_ieta_HB1
MonitorElement *	emean_vs_ieta_HB2
MonitorElement *	emean_vs_ieta_HE1
MonitorElement *	emean_vs_ieta_HE2
MonitorElement *	emean_vs_ieta_HE3
MonitorElement *	emean_vs_ieta_HF1
MonitorElement *	emean_vs_ieta_HF2
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_
int	imc
int	iz
MonitorElement *	map_depth1
MonitorElement *	map_depth2
MonitorElement *	map_depth3
MonitorElement *	map_depth4
MonitorElement *	map_ecal
MonitorElement *	map_econe_depth1
MonitorElement *	map_econe_depth2
MonitorElement *	map_econe_depth3
MonitorElement *	map_econe_depth4
std::string	mc_
MonitorElement *	meDeltaEta
MonitorElement *	meDeltaEtaS
MonitorElement *	meDeltaPhi
MonitorElement *	meDeltaPhiS
MonitorElement *	meEcalHcalEnergyConeHB
MonitorElement *	meEcalHcalEnergyConeHE
MonitorElement *	meEcalHcalEnergyConeHF
MonitorElement *	meEcalHcalEnergyConeHO
MonitorElement *	meEcalHcalEnergyHB
MonitorElement *	meEcalHcalEnergyHE
MonitorElement *	meEnConeEtaProfile
MonitorElement *	meEnConeEtaProfile_depth1
MonitorElement *	meEnConeEtaProfile_depth2
MonitorElement *	meEnConeEtaProfile_depth3
MonitorElement *	meEnConeEtaProfile_depth4
MonitorElement *	meEnConeEtaProfile_E
MonitorElement *	meEnConeEtaProfile_EH
MonitorElement *	meEnergyHcalVsEcalHB
MonitorElement *	meEnergyHcalVsEcalHE
MonitorElement *	meNumEcalRecHitsConeHB
MonitorElement *	meNumEcalRecHitsConeHE
MonitorElement *	meNumRecHitsConeHB
MonitorElement *	meNumRecHitsConeHE
MonitorElement *	meNumRecHitsConeHF
MonitorElement *	meNumRecHitsConeHO
MonitorElement *	meNumRecHitsThreshHB
MonitorElement *	meNumRecHitsThreshHE
MonitorElement *	meNumRecHitsThreshHO
MonitorElement *	meRecHitsEnergyHB
MonitorElement *	meRecHitsEnergyHE
MonitorElement *	meRecHitsEnergyHF
MonitorElement *	meRecHitsEnergyHO
MonitorElement *	meRecHitSimHitHB
MonitorElement *	meRecHitSimHitHE
MonitorElement *	meRecHitSimHitHF
MonitorElement *	meRecHitSimHitHFL
MonitorElement *	meRecHitSimHitHFS
MonitorElement *	meRecHitSimHitHO
MonitorElement *	meRecHitSimHitProfileHB
MonitorElement *	meRecHitSimHitProfileHE
MonitorElement *	meRecHitSimHitProfileHF
MonitorElement *	meRecHitSimHitProfileHFL
MonitorElement *	meRecHitSimHitProfileHFS
MonitorElement *	meRecHitSimHitProfileHO
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_HB1
MonitorElement *	meTE_HB2
MonitorElement *	meTE_HE
MonitorElement *	meTE_HE1
MonitorElement *	meTE_HE2
MonitorElement *	meTE_HF
MonitorElement *	meTE_HFL
MonitorElement *	meTE_HFS
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	nevtot
MonitorElement *	Nhb
MonitorElement *	Nhe
MonitorElement *	Nhf
MonitorElement *	Nho
MonitorElement *	occupancy_map_HB1
MonitorElement *	occupancy_map_HB2
MonitorElement *	occupancy_map_HE1
MonitorElement *	occupancy_map_HE2
MonitorElement *	occupancy_map_HE3
MonitorElement *	occupancy_map_HF1
MonitorElement *	occupancy_map_HF2
MonitorElement *	occupancy_map_HO
MonitorElement *	occupancy_seqHB1
MonitorElement *	occupancy_seqHB2
MonitorElement *	occupancy_seqHE1
MonitorElement *	occupancy_seqHE2
MonitorElement *	occupancy_seqHE3
MonitorElement *	occupancy_seqHF1
MonitorElement *	occupancy_seqHF2
MonitorElement *	occupancy_seqHO
MonitorElement *	occupancy_vs_ieta_HB1
MonitorElement *	occupancy_vs_ieta_HB2
MonitorElement *	occupancy_vs_ieta_HE1
MonitorElement *	occupancy_vs_ieta_HE2
MonitorElement *	occupancy_vs_ieta_HE3
MonitorElement *	occupancy_vs_ieta_HF1
MonitorElement *	occupancy_vs_ieta_HF2
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 *	RMS_seq_HB1
MonitorElement *	RMS_seq_HB2
MonitorElement *	RMS_seq_HE1
MonitorElement *	RMS_seq_HE2
MonitorElement *	RMS_seq_HE3
MonitorElement *	RMS_seq_HF1
MonitorElement *	RMS_seq_HF2
MonitorElement *	RMS_seq_HO
MonitorElement *	RMS_vs_ieta_HB1
MonitorElement *	RMS_vs_ieta_HB2
MonitorElement *	RMS_vs_ieta_HE1
MonitorElement *	RMS_vs_ieta_HE2
MonitorElement *	RMS_vs_ieta_HE3
MonitorElement *	RMS_vs_ieta_HF1
MonitorElement *	RMS_vs_ieta_HF2
MonitorElement *	RMS_vs_ieta_HO
MonitorElement *	sevLvl_HB
MonitorElement *	sevLvl_HE
MonitorElement *	sevLvl_HF
MonitorElement *	sevLvl_HO
std::string	sign_
int	subdet_
const HcalChannelQuality *	theHcalChStatus
const HcalSeverityLevelComputer *	theHcalSevLvlComputer
edm::EDGetTokenT < EBRecHitCollection >	tok_EB_
edm::EDGetTokenT < EERecHitCollection >	tok_EE_
edm::EDGetTokenT < edm::HepMCProduct >	tok_evt_
edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_
edm::EDGetTokenT < HFRecHitCollection >	tok_hf_
edm::EDGetTokenT < edm::PCaloHitContainer >	tok_hh_
edm::EDGetTokenT < HORecHitCollection >	tok_ho_
bool	useAllHistos_
MonitorElement *	ZS_HB1
MonitorElement *	ZS_HB2
MonitorElement *	ZS_HE1
MonitorElement *	ZS_HE2
MonitorElement *	ZS_HE3
MonitorElement *	ZS_HF1
MonitorElement *	ZS_HF2
MonitorElement *	ZS_HO
MonitorElement *	ZS_nHB1
MonitorElement *	ZS_nHB2
MonitorElement *	ZS_nHE1
MonitorElement *	ZS_nHE2
MonitorElement *	ZS_nHE3
MonitorElement *	ZS_nHF1
MonitorElement *	ZS_nHF2
MonitorElement *	ZS_nHO
MonitorElement *	ZS_seqHB1
MonitorElement *	ZS_seqHB2
MonitorElement *	ZS_seqHE1
MonitorElement *	ZS_seqHE2
MonitorElement *	ZS_seqHE3
MonitorElement *	ZS_seqHF1
MonitorElement *	ZS_seqHF2
MonitorElement *	ZS_seqHO

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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)

Detailed Description

Definition at line 56 of file HcalRecHitsValidation.h.

Constructor & Destructor Documentation

HcalRecHitsValidation::HcalRecHitsValidation

(

edm::ParameterSet const &

conf

)

Definition at line 5 of file HcalRecHitsValidation.cc.

References ecalselector_, etype_, eventype_, famos_, edm::ParameterSet::getUntrackedParameter(), hcalselector_, imc, iz, mc_, nevtot, outputFile_, sign_, AlCaHLTBitMon_QueryRunRegistry::string, subdet_, tok_EB_, tok_EE_, tok_evt_, tok_hbhe_, tok_hf_, tok_hh_, tok_ho_, and useAllHistos_.

                                                                        {
  // DQM ROOT output
  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
  
  if ( outputFile_.size() != 0 ) {
    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", "*");
  mc_           = conf.getUntrackedParameter<std::string>("mc", "yes");
  famos_        = conf.getUntrackedParameter<bool>("Famos", false);
  useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);

  //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"));

  // register for data access
  tok_evt_ = consumes<edm::HepMCProduct>(edm::InputTag("generator"));
  tok_EB_ = consumes<EBRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
  tok_EE_ = consumes<EERecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
  tok_hh_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits","HcalHits"));

  //  std::cout << "*** famos_ = " << famos_ << std::endl; 

  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 = 1;
  if(mc_ == "no") imc = 0;

}

HcalRecHitsValidation::~HcalRecHitsValidation

(

)

Definition at line 60 of file HcalRecHitsValidation.cc.

{ }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 701 of file HcalRecHitsValidation.cc.

References cauxstwd, cdepth, cen, ceta, cieta, ciphi, gather_cfg::cout, cphi, cstwd, csub, ctime, dPhiWsign(), dR(), e_hb, e_he, e_hfl, e_hfs, e_ho, ecalselector_, emap_depth1, emap_depth2, emap_depth3, emap_depth4, emap_min, emean_seqHB1, emean_seqHB2, emean_seqHE1, emean_seqHE2, emean_seqHE3, emean_seqHF1, emean_seqHF2, emean_seqHO, emean_vs_ieta_HB1, emean_vs_ieta_HB2, emean_vs_ieta_HE1, emean_vs_ieta_HE2, emean_vs_ieta_HE3, emean_vs_ieta_HF1, emean_vs_ieta_HF2, emean_vs_ieta_HO, PV3DBase< T, PVType, FrameType >::eta(), eta(), etype_, famos_, MonitorElement::Fill(), HcalObjRepresent::Fill(), fillRecHitsTmp(), geometry, edm::EventSetup::get(), edm::Event::getByToken(), CaloCellGeometry::getPosition(), hcalHBSevLvlVec, hcalHESevLvlVec, hcalHFSevLvlVec, hcalHOSevLvlVec, hcalselector_, i, imc, cmsHarvester::index, edm::HandleBase::isValid(), map_ecal, map_econe_depth1, map_econe_depth2, map_econe_depth3, map_econe_depth4, meDeltaEta, meDeltaEtaS, meDeltaPhi, meDeltaPhiS, meEcalHcalEnergyConeHB, meEcalHcalEnergyConeHE, meEcalHcalEnergyHB, meEcalHcalEnergyHE, meEnConeEtaProfile, meEnConeEtaProfile_depth1, meEnConeEtaProfile_depth2, meEnConeEtaProfile_depth3, meEnConeEtaProfile_depth4, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meEnergyHcalVsEcalHB, meEnergyHcalVsEcalHE, meNumEcalRecHitsConeHB, meNumEcalRecHitsConeHE, meNumRecHitsConeHB, meNumRecHitsConeHE, meNumRecHitsConeHF, meNumRecHitsConeHO, meNumRecHitsThreshHB, meNumRecHitsThreshHE, meNumRecHitsThreshHO, meRecHitsEnergyHB, meRecHitsEnergyHE, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitSimHitHB, meRecHitSimHitHE, meRecHitSimHitHF, meRecHitSimHitHFL, meRecHitSimHitHFS, meRecHitSimHitHO, meRecHitSimHitProfileHB, meRecHitSimHitProfileHE, meRecHitSimHitProfileHF, meRecHitSimHitProfileHFL, meRecHitSimHitProfileHFS, meRecHitSimHitProfileHO, meSumRecHitsEnergyConeHB, meSumRecHitsEnergyConeHE, meSumRecHitsEnergyConeHF, meSumRecHitsEnergyConeHFL, meSumRecHitsEnergyConeHFS, meSumRecHitsEnergyConeHO, meSumRecHitsEnergyHB, meSumRecHitsEnergyHE, meSumRecHitsEnergyHF, meSumRecHitsEnergyHO, meTE_HB, meTE_HB1, meTE_HB2, meTE_HE, meTE_HE1, meTE_HE2, meTE_HF, meTE_HFL, meTE_HFS, 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, npart, occupancy_map_HB1, occupancy_map_HB2, occupancy_map_HE1, occupancy_map_HE2, occupancy_map_HE3, occupancy_map_HF1, occupancy_map_HF2, occupancy_map_HO, AlCaHLTBitMon_ParallelJobs::p, phi, PV3DBase< T, PVType, FrameType >::phi(), phi12(), edm::Handle< T >::product(), RecoTauCleanerPlugins::pt, alignCSCRings::r, 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_HF67, RecHit_StatusWord_HO, HcalCondObjectContainerBase::setTopo(), sevLvl_HB, sevLvl_HE, sevLvl_HF, sevLvl_HO, HcalDetId::subdet(), subdet_, edmStreamStallGrapher::t, theHcalChStatus, theHcalSevLvlComputer, tok_EB_, tok_EE_, tok_evt_, tok_hh_, HcalCondObjectContainerBase::topo(), useAllHistos_, ZS_nHB1, ZS_nHB2, ZS_nHE1, ZS_nHE2, ZS_nHE3, ZS_nHF1, ZS_nHF2, and ZS_nHO.

                                                                              {

  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.;
  double eHcalCone    = 0.;  
  double eHcalConeHB  = 0.;  
  double eHcalConeHE  = 0.;  
  double eHcalConeHO  = 0.;  
  double eHcalConeHF  = 0.;  
  double eHcalConeHFL = 0.;  
  double eHcalConeHFS = 0.;  
  // Total numbet of RecHits in HCAL, in the cone, above 1 GeV theshold
  int nrechits       = 0;
  int nrechitsCone   = 0;
  int nrechitsThresh = 0;

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

  // MC info 
  double phi_MC = -999999.;  // phi of initial particle from HepMC
  double eta_MC = -999999.;  // eta of initial particle from HepMC

  // HCAL energy around MC eta-phi at all depths;
  double partR = 0.3;
  double ehcal_coneMC_1 = 0.;
  double ehcal_coneMC_2 = 0.;
  double ehcal_coneMC_3 = 0.;
  double ehcal_coneMC_4 = 0.;

  // Cone size for serach of the hottest HCAL cell around MC
  double searchR = 1.0; 
  double eps     = 0.001;

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

  // MC information

  //  std::cout << "*** 1" << std::endl; 


  if(imc != 0) { 

  edm::Handle<edm::HepMCProduct> evtMC;
  ev.getByToken(tok_evt_,evtMC);  // generator in late 310_preX
  if (!evtMC.isValid()) {
    std::cout << "no HepMCProduct found" << std::endl;    
  } else {
    //    std::cout << "*** source HepMCProduct found"<< std::endl;
  }  

  // 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();
    //    phi_MC = phip;
    //    eta_MC = etap;
    double pt  = (*p)->momentum().perp();
    if(pt > maxPt) { npart++; maxPt = pt; phi_MC = phip; eta_MC = etap; }
  }
  //  std::cout << "*** Max pT = " << maxPt <<  std::endl;  

  }

  //   std::cout << "*** 2" << std::endl; 
  //   previously was:  c.get<IdealGeometryRecord>().get (geometry);
  c.get<CaloGeometryRecord>().get (geometry);

  // HCAL channel status map ****************************************
  edm::ESHandle<HcalChannelQuality> hcalChStatus;
  c.get<HcalChannelQualityRcd>().get( hcalChStatus );

  edm::ESHandle<HcalTopology> topo;
  c.get<IdealGeometryRecord>().get(topo);

  theHcalChStatus = hcalChStatus.product();

  if( !theHcalChStatus->topo() ) theHcalChStatus->setTopo(topo.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]);
     }
  } 

  //  std::cout << "*** 3" << std::endl; 


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

  // ECAL 
  if(ecalselector_ == "yes" && (subdet_ == 1 || subdet_ == 2 || subdet_ == 5)) {
    Handle<EBRecHitCollection> rhitEB;


      ev.getByToken(tok_EB_, rhitEB);

    EcalRecHitCollection::const_iterator RecHit = rhitEB.product()->begin();  
    EcalRecHitCollection::const_iterator RecHitEnd = rhitEB.product()->end();  
    
    for (; RecHit != RecHitEnd ; ++RecHit) {
      EBDetId EBid = EBDetId(RecHit->id());
       
      const CaloCellGeometry* cellGeometry =
    geometry->getSubdetectorGeometry (EBid)->getGeometry (EBid) ;
      double eta = cellGeometry->getPosition ().eta () ;
      double phi = cellGeometry->getPosition ().phi () ;
      double en  = RecHit->energy();
      eEcal  += en;
      eEcalB += en;

      if (useAllHistos_) map_ecal->Fill(eta, phi, en);

      double r   = dR(eta_MC, phi_MC, eta, phi);
      if( r < partR)  {
    eEcalCone += en;
    numrechitsEcal++; 
      }
    }

    
    Handle<EERecHitCollection> rhitEE;
 
      ev.getByToken(tok_EE_, rhitEE);

    RecHit = rhitEE.product()->begin();  
    RecHitEnd = rhitEE.product()->end();  
    
    for (; RecHit != RecHitEnd ; ++RecHit) {
      EEDetId EEid = EEDetId(RecHit->id());
      
      const CaloCellGeometry* cellGeometry =
    geometry->getSubdetectorGeometry (EEid)->getGeometry (EEid) ;
      double eta = cellGeometry->getPosition ().eta () ;
      double phi = cellGeometry->getPosition ().phi () ;    
      double en   = RecHit->energy();
      eEcal  += en;
      eEcalE += en;

      if (useAllHistos_) map_ecal->Fill(eta, phi, en);

      double r   = dR(eta_MC, phi_MC, eta, phi);
      if( r < partR)  {
    eEcalCone += en;
    numrechitsEcal++; 
      }
    }
  }     // end of ECAL selection 


  //     std::cout << "*** 4" << std::endl; 


  // Counting, including ZS items
  // Filling HCAL maps  ----------------------------------------------------
  double maxE = -99999.;
  
  int nhb1 = 0;
  int nhb2 = 0;
  int nhe1 = 0;
  int nhe2 = 0;
  int nhe3 = 0;
  int nho  = 0;
  int nhf1 = 0;
  int nhf2 = 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 eta    = ceta[i]; 
    double phi    = cphi[i]; 
    uint32_t stwd = cstwd[i];
    uint32_t auxstwd = cauxstwd[i];
    //    double z   = cz[i];

    int index = ieta * 72 + iphi; //  for sequential histos
    
    /*   
     std::cout << "*** point 4-1" << " ieta, iphi, depth, sub = "
     << ieta << ", " << iphi << ", " << depth << ", " << sub  
     << std::endl;
    */
    
    
    if( sub == 1 && depth == 1)  nhb1++;
    if( sub == 1 && depth == 2)  nhb2++;
    if( sub == 2 && depth == 1)  nhe1++;
    if( sub == 2 && depth == 2)  nhe2++;
    if( sub == 2 && depth == 3)  nhe3++;
    if( sub == 3 && depth == 4)  nho++;
    if( sub == 4 && depth == 1)  nhf1++;
    if( sub == 4 && depth == 2)  nhf2++;
    
    if( subdet_ == 6) {                                    // ZS specific
      if( en < emap_min[ieta+41][iphi][depth-1][sub-1] )
    emap_min[ieta+41][iphi][depth-1][sub-1] = en;
    }
    
    double emin = 1.;
    if(fabs(eta) > 3.) emin = 5.; 
    
    double r  = dR(eta_MC, phi_MC, eta, phi);
    if( r < searchR ) { // search for hottest cell in a big cone
      if(maxE < en && en > emin) {
    maxE    = en;
    etaHot  = eta;
    phiHot  = phi;
      }
    }

    /*   
    if(ieta == 27 ) { 
      std::cout << "*** ieta=28, iphi = " << iphi << "  det = " 
        << sub << "  depth = " << depth << std::endl;
    }
    */

    if( subdet_ != 6) {  

      //      std::cout << "*** 4-1" << std::endl; 
      //The emean_vs_ieta histos are drawn as well as the e_maps


      // to distinguish HE and HF
      if( depth == 1 || depth == 2 ) {
        int ieta1 =  ieta;
    if(sub == 4) { 
      if (ieta1 < 0) ieta1--;
          else  ieta1++;   
    }
    if (depth == 1) emap_depth1->Fill(double(ieta1), double(iphi), en);
    if (depth == 2) emap_depth2->Fill(double(ieta1), double(iphi), en);
      }

      if( depth == 3) emap_depth3->Fill(double(ieta), double(iphi), en);
      if( depth == 4) emap_depth4->Fill(double(ieta), double(iphi), en);
      
      if (depth == 1 && sub == 1 ) {
    emean_vs_ieta_HB1->Fill(double(ieta), en);
    occupancy_map_HB1->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHB1->Fill(double(index), en);
    }
      }
      if (depth == 2  && sub == 1) {
    emean_vs_ieta_HB2->Fill(double(ieta), en);
    occupancy_map_HB2->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHB2->Fill(double(index), en);
    }
      }
      if (depth == 1 && sub == 2) {
    emean_vs_ieta_HE1->Fill(double(ieta), en);
    occupancy_map_HE1->Fill(double(ieta), double(iphi));   
    if(useAllHistos_){
      emean_seqHE1->Fill(double(index), en);
    }
      }
      if (depth == 2 && sub == 2) {
    emean_vs_ieta_HE2->Fill(double(ieta), en);
    occupancy_map_HE2->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHE2->Fill(double(index), en);
    }
      }
      if (depth == 3 && sub == 2) {
    emean_vs_ieta_HE3->Fill(double(ieta), en);
    occupancy_map_HE3->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHE3->Fill(double(index), en);
    }
      }
      if (depth == 4 ) {
    emean_vs_ieta_HO->Fill(double(ieta), en);
    occupancy_map_HO->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHO->Fill(double(index), en);
    }
      }
      if (depth == 1 && sub == 4) {
    emean_vs_ieta_HF1->Fill(double(ieta), en);
    occupancy_map_HF1->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHF1->Fill(double(index), en);
    }
      }
      if (depth == 2 && sub == 4) {
    emean_vs_ieta_HF2->Fill(double(ieta), en);
    occupancy_map_HF2->Fill(double(ieta), double(iphi));          
    if(useAllHistos_){
      emean_seqHF2->Fill(double(index), en);
    }
      }
    }
    

    if( r < partR ) {
      if (depth == 1) ehcal_coneMC_1 += en; 
      if (depth == 2) ehcal_coneMC_2 += en; 
      if (depth == 3) ehcal_coneMC_3 += en; 
      if (depth == 4) ehcal_coneMC_4 += en; 
    }
    
    //32-bit status word  
    uint32_t statadd;
    unsigned int isw67 = 0;
    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;
    }
      }
    }
    if (isw67 != 0 && useAllHistos_) RecHit_StatusWord_HF67->Fill(isw67); //This one is not drawn

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

    }

  } 
 
  //  std::cout << "*** 4-2" << std::endl; 
  
  if( subdet_ == 6 && useAllHistos_) {               // ZS plots; not drawn
    ZS_nHB1->Fill(double(nhb1));  
    ZS_nHB2->Fill(double(nhb2));  
    ZS_nHE1->Fill(double(nhe1));  
    ZS_nHE2->Fill(double(nhe2));  
    ZS_nHE3->Fill(double(nhe3));  
    ZS_nHO ->Fill(double(nho));  
    ZS_nHF1->Fill(double(nhf1));  
    ZS_nHF2->Fill(double(nhf2));  
  }
  else{ 
    Nhb->Fill(double(nhb1 + nhb2));
    Nhe->Fill(double(nhe1 + nhe2 + nhe3));
    Nho->Fill(double(nho));
    Nhf->Fill(double(nhf1 + nhf2));

    //These are not drawn
    if(imc != 0 && useAllHistos_) {
      map_econe_depth1->Fill(eta_MC, phi_MC, ehcal_coneMC_1);
      map_econe_depth2->Fill(eta_MC, phi_MC, ehcal_coneMC_2);
      map_econe_depth3->Fill(eta_MC, phi_MC, ehcal_coneMC_3);
      map_econe_depth4->Fill(eta_MC, phi_MC, ehcal_coneMC_4);
    }
  }

  //  std::cout << "*** 5" << std::endl; 
    

  //  NOISE ================================================================= 
  //Not drawn
  if (hcalselector_ == "noise" && useAllHistos_) {
    for (unsigned int i = 0; i < cen.size(); i++) {
      
      int sub   = csub[i];
      int depth = cdepth[i];
      double en = cen[i]; 
      
      if (sub == 1) e_hb->Fill(en);
      if (sub == 2) e_he->Fill(en);  
      if (sub == 3) e_ho->Fill(en);  
      if (sub == 4) {
    if(depth == 1)  
      e_hfl->Fill(en);  
    else 
      e_hfs->Fill(en);  
      }
    }
  }

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

    //       std::cout << "*** 6" << std::endl; 
    
    
    double clusEta = 999.;
    double clusPhi = 999.; 
    double clusEn  = 0.;
    
    double HcalCone_d1 = 0.;
    double HcalCone_d2 = 0.;
    double HcalCone_d3 = 0.;
    double HcalCone_d4 = 0.;
    double HcalCone    = 0.;

    int ietaMax1  =  9999;
    int ietaMax2  =  9999;
    int ietaMax3  =  9999;
    int ietaMax4  =  9999;
    int ietaMax   =  9999;
    double enMax1 = -9999.;
    double enMax2 = -9999.;
    double enMax3 = -9999.;
    double enMax4 = -9999.;
    //    double enMax  = -9999.;
    double etaMax =  9999.;

    /*
    std::cout << "*** point 5-1" << "  eta_MC, phi_MC    etaHot,  phiHot = "
          << eta_MC  << ", " << phi_MC << "   "
          << etaHot  << ", " << phiHot  
          << std::endl;
    */

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

    for (unsigned int i = 0; i < cen.size(); i++) {
      int sub    = csub[i];
      int depth  = cdepth[i];
      double eta = ceta[i]; 
      double phi = cphi[i]; 
      double en  = cen[i]; 
      double t   = ctime[i];
      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++;

      double r    = dR(eta_MC, phi_MC, eta, phi);
      if( r < partR ){
        if(sub == 1)   eHcalConeHB += en;
        if(sub == 2)   eHcalConeHE += en;
        if(sub == 3)   eHcalConeHO += en;
        if(sub == 4) {
      eHcalConeHF += en;
      if (depth == 1) eHcalConeHFL += en;
      else            eHcalConeHFS += en;
    }
    eHcalCone += en;
    nrechitsCone++;

    // search for most energetic cell at the given depth in the cone
        if(depth == 1) {
      HcalCone_d1 += en;
      if(enMax1 < en) {
        enMax1   = en;
        ietaMax1 = ieta;
      }
    }
        if(depth == 2) {
      HcalCone_d2 += en;
      if(enMax2 < en) {
        enMax2   = en;
        ietaMax2 = ieta;
      }
    }
        if(depth == 3) {
      HcalCone_d3 += en;
      if(enMax3 < en) {
        enMax3   = en;
        ietaMax3 = ieta;
      }
    }
        if(depth == 4) {
      HcalCone_d4 += en;
      if(enMax4 < en) {
        enMax4   = en;
        ietaMax4 = ieta;
      }
    }

    if(depth != 4) {
      HcalCone += en;
    }


    // regardless of the depths (but excluding HO), just hottest cell
    /*
    if(depth != 4) {
      if(enMax   < en) {
        enMax   = en;
        ietaMax = ieta;
      }
    }   
    */

        // alternative: ietamax -> closest to MC eta  !!!
    float eta_diff = fabs(eta_MC - eta);
    if(eta_diff < etaMax) {
      etaMax  = eta_diff; 
          ietaMax = ieta; 
    }
      }
      
      //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);
    
    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 (useAllHistos_){
      if      (depth == 1) meTE_HB1->Fill( en, t);
      else if (depth == 2) meTE_HB2->Fill( en, t);
    }
      }     
      if(sub == 2 && (subdet_ == 2 || subdet_ == 5)) {  
    meTimeHE->Fill(t);
    meRecHitsEnergyHE->Fill(en);

    meTE_Low_HE->Fill( en, t);
    meTE_HE->Fill( en, t);
    meTEprofileHE_Low->Fill(en, t);
    meTEprofileHE->Fill(en, t);

    if (useAllHistos_){
      if      (depth == 1) meTE_HE1->Fill( en, t);
      else if (depth == 2) meTE_HE2->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 (useAllHistos_){
      if   (depth == 1) meTE_HFL->Fill( en, t);
      else              meTE_HFS->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
      if (useAllHistos_){
    meEnConeEtaProfile_depth1->Fill(double(ietaMax1), HcalCone_d1);
    meEnConeEtaProfile_depth2->Fill(double(ietaMax2), HcalCone_d2);
    meEnConeEtaProfile_depth3->Fill(double(ietaMax3), HcalCone_d3);
    meEnConeEtaProfile_depth4->Fill(double(ietaMax4), HcalCone_d4);
      }
      meEnConeEtaProfile       ->Fill(double(ietaMax),  HcalCone);   // 
      meEnConeEtaProfile_E     ->Fill(double(ietaMax), eEcalCone);   
      meEnConeEtaProfile_EH    ->Fill(double(ietaMax),  HcalCone+eEcalCone); 
    }

    //     std::cout << "*** 7" << std::endl; 

    
    // Single particle samples ONLY !  ======================================
    // Fill up some histos for "integrated" subsustems. 
    // These are not drawn
    if(etype_ == 1 && useAllHistos_) {

      /*
      std::cout << "*** point 7-1" << "  eta_MC, phi_MC   clusEta, clusPhi = "
                << eta_MC  << ", " << phi_MC << "   "
        << clusEta << ", " << clusPhi 
        << std::endl;
      */    

      double phidev = dPhiWsign(clusPhi, phi_MC);
      meDeltaPhi->Fill(eta_MC, phidev);
      double etadev = clusEta - eta_MC;
      meDeltaEta->Fill(eta_MC, etadev);

      if(subdet_ == 1) {
    meSumRecHitsEnergyHB->Fill(eHcal);
    if(imc != 0) meSumRecHitsEnergyConeHB->Fill(eHcalConeHB);    
    if(imc != 0) meNumRecHitsConeHB->Fill(double(nrechitsCone));
    meNumRecHitsThreshHB->Fill(double(nrechitsThresh));
      }

      if(subdet_ == 2) {
        meSumRecHitsEnergyHE->Fill(eHcal);
    if(imc != 0) meSumRecHitsEnergyConeHE->Fill(eHcalConeHE);    
    if(imc != 0) meNumRecHitsConeHE->Fill(double(nrechitsCone));
    meNumRecHitsThreshHE->Fill(double(nrechitsThresh));
      }

      if(subdet_ == 3) {
    meSumRecHitsEnergyHO->Fill(eHcal);
    if(imc != 0) meSumRecHitsEnergyConeHO->Fill(eHcalConeHO);    
    if(imc != 0) meNumRecHitsConeHO->Fill(double(nrechitsCone));
    meNumRecHitsThreshHO->Fill(double(nrechitsThresh));
      }

      if(subdet_ == 4) {
        if(eHcalConeHF > eps ) {
      meSumRecHitsEnergyHF ->Fill(eHcal);
      if(imc != 0) { 
        meSumRecHitsEnergyConeHF ->Fill(eHcalConeHF);    
        meNumRecHitsConeHF->Fill(double(nrechitsCone));
        meSumRecHitsEnergyConeHFL ->Fill(eHcalConeHFL);    
        meSumRecHitsEnergyConeHFS ->Fill(eHcalConeHFS);    
      }
    }
      }

      //         std::cout << "*** 8" << std::endl; 


      // Also combine with ECAL if needed 
      if(subdet_ == 1  && ecalselector_ == "yes") {
    
    /*
      std::cout << "*** point 8-1" 
      << "  eEcalB " << eEcalB << "  eHcal " << eHcal
      << "  eEcalCone " <<  eEcalCone << "  eHcalCone " 
          << eHcalCone
          << "  numrechitsEcal " <<  numrechitsEcal
          << std::endl;
          
    */
    
        meEcalHcalEnergyHB->Fill(eEcalB+eHcal);
        meEcalHcalEnergyConeHB->Fill(eEcalCone+eHcalCone);
        meNumEcalRecHitsConeHB->Fill(double(numrechitsEcal));
    
      }
      
      if(subdet_ == 2  && ecalselector_ == "yes"){
    
    /*
      std::cout << "*** point 8-2a" 
      << "  eEcalE " << eEcalE << "  eHcal " << eHcal
      << "  eEcalCone " <<  eEcalCone << "  eHcalCone " 
      << eHcalCone
      << "  numrechitsEcal " <<  numrechitsEcal
      << std::endl;
    */
    
    meEcalHcalEnergyHE->Fill(eEcalE+eHcal);
    if(imc != 0) meEcalHcalEnergyConeHE->Fill(eEcalCone+eHcalCone);
    if(imc != 0) meNumEcalRecHitsConeHE->Fill(double(numrechitsEcal));
      } 

      // Banana plots finally
      if(imc != 0) {
    if(subdet_ == 1 && ecalselector_ == "yes")
      meEnergyHcalVsEcalHB -> Fill(eEcalCone,eHcalCone);
    if(subdet_ == 2 && ecalselector_ == "yes") 
      meEnergyHcalVsEcalHE -> Fill(eEcalCone,eHcalCone);
      }
    }
  }
  //  std::cout << "*** 9" << std::endl; 


  //===========================================================================
  // Getting SimHits
  //===========================================================================

  if(subdet_ > 0 && subdet_ < 6 && imc !=0 && !famos_  ) {  // not noise 

    double maxES = -9999.;
    double etaHotS = 1000.;
    double phiHotS = 1000.;
    
    edm::Handle<PCaloHitContainer> hcalHits;
    ev.getByToken(tok_hh_,hcalHits);
    const PCaloHitContainer * SimHitResult = hcalHits.product () ;
    
    double enSimHits    = 0.;
    double enSimHitsHB  = 0.;
    double enSimHitsHE  = 0.;
    double enSimHitsHO  = 0.;
    double enSimHitsHF  = 0.;
    double enSimHitsHFL = 0.;
    double enSimHitsHFS = 0.;
    // sum of SimHits in the cone 
    
    for (std::vector<PCaloHit>::const_iterator SimHits = SimHitResult->begin () ; SimHits != SimHitResult->end(); ++SimHits) {
      HcalDetId cell(SimHits->id());
      int sub =  cell.subdet();
      const CaloCellGeometry* cellGeometry =
    geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
      double etaS = cellGeometry->getPosition().eta () ;
      double phiS = cellGeometry->getPosition().phi () ;
      double en   = SimHits->energy();    

      double emin = 0.01;
      if(fabs(etaS) > 3.) emin = 1.;   

      double r  = dR(eta_MC, phi_MC, etaS, phiS);
      if( r < searchR ) { // search for hottest cell in a big cone
    if(maxES < en && en > emin ) {
      maxES    = en;
      etaHotS  = etaS;
      phiHotS  = phiS;
    }
      }
       
      if ( r < partR ){ // just energy in the small cone
    enSimHits += en;
    if(sub == 1) enSimHitsHB += en; 
    if(sub == 2) enSimHitsHE += en; 
    if(sub == 3) enSimHitsHO += en; 
    if(sub == 4) {
      enSimHitsHF += en;
      int depth = cell.depth();
      if(depth == 1) enSimHitsHFL += en;
      else           enSimHitsHFS += en;
    } 
      }
    }


    // Second look over SimHits: cluster finding

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

    for (std::vector<PCaloHit>::const_iterator SimHits = SimHitResult->begin () ; SimHits != SimHitResult->end(); ++SimHits) {
     HcalDetId cell(SimHits->id());

      const CaloCellGeometry* cellGeometry =
    geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
      double etaS = cellGeometry->getPosition().eta () ;
      double phiS = cellGeometry->getPosition().phi () ;
      double en   =  SimHits->energy();    

      double emin = 0.01;
      if(fabs(etaS) > 3.) emin = 1.; 

      double rhot = dR(etaHotS, phiHotS, etaS, phiS); 
      if(rhot < partR && en > emin) { 
    clusEta = (clusEta * clusEn + etaS * en)/(clusEn + en);
        clusPhi = phi12(clusPhi, clusEn, phiS, en); 
        clusEn += en;
      }
    }

    // SimHits cluster deviation from MC (eta, phi)
    // These are not drawn
    if (useAllHistos_){
      if(etype_ == 1) {
    double phidev = dPhiWsign(clusPhi, phi_MC);
    meDeltaPhiS->Fill(eta_MC, phidev);
    double etadev = clusEta - eta_MC;
    meDeltaEtaS->Fill(eta_MC, etadev);
      }
      // Now some histos with SimHits
    
      if(subdet_ == 4 || subdet_ == 5) {
    if(eHcalConeHF > eps) {
      meRecHitSimHitHF->Fill( enSimHitsHF, eHcalConeHF );
      meRecHitSimHitProfileHF->Fill( enSimHitsHF, eHcalConeHF);
      
      meRecHitSimHitHFL->Fill( enSimHitsHFL, eHcalConeHFL );
      meRecHitSimHitProfileHFL->Fill( enSimHitsHFL, eHcalConeHFL);
      meRecHitSimHitHFS->Fill( enSimHitsHFS, eHcalConeHFS );
      meRecHitSimHitProfileHFS->Fill( enSimHitsHFS, eHcalConeHFS);       
    }
      }
      if(subdet_ == 1  || subdet_ == 5) { 
    meRecHitSimHitHB->Fill( enSimHitsHB,eHcalConeHB );
    meRecHitSimHitProfileHB->Fill( enSimHitsHB,eHcalConeHB);
      }
      if(subdet_ == 2  || subdet_ == 5) { 
    meRecHitSimHitHE->Fill( enSimHitsHE,eHcalConeHE );
    meRecHitSimHitProfileHE->Fill( enSimHitsHE,eHcalConeHE);
      }
      if(subdet_ == 3  || subdet_ == 5) { 
    meRecHitSimHitHO->Fill( enSimHitsHO,eHcalConeHO );
    meRecHitSimHitProfileHO->Fill( enSimHitsHO,eHcalConeHO);
      }
    }
  }

  nevtot++;
}

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

virtual

Implements thread_unsafe::DQMEDAnalyzer.

Definition at line 62 of file HcalRecHitsValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), e_hb, e_he, e_hfl, e_hfs, e_ho, ecalselector_, emap_depth1, emap_depth2, emap_depth3, emap_depth4, emap_min, emean_seqHB1, emean_seqHB2, emean_seqHE1, emean_seqHE2, emean_seqHE3, emean_seqHF1, emean_seqHF2, emean_seqHO, emean_vs_ieta_HB1, emean_vs_ieta_HB2, emean_vs_ieta_HE1, emean_vs_ieta_HE2, emean_vs_ieta_HE3, emean_vs_ieta_HF1, emean_vs_ieta_HF2, emean_vs_ieta_HO, etype_, hcalselector_, timingPdfMaker::histo, imc, map_depth1, map_depth2, map_depth3, map_depth4, map_ecal, map_econe_depth1, map_econe_depth2, map_econe_depth3, map_econe_depth4, meDeltaEta, meDeltaEtaS, meDeltaPhi, meDeltaPhiS, meEcalHcalEnergyConeHB, meEcalHcalEnergyConeHE, meEcalHcalEnergyHB, meEcalHcalEnergyHE, meEnConeEtaProfile, meEnConeEtaProfile_depth1, meEnConeEtaProfile_depth2, meEnConeEtaProfile_depth3, meEnConeEtaProfile_depth4, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meEnergyHcalVsEcalHB, meEnergyHcalVsEcalHE, meNumEcalRecHitsConeHB, meNumEcalRecHitsConeHE, meNumRecHitsConeHB, meNumRecHitsConeHE, meNumRecHitsConeHF, meNumRecHitsConeHO, meNumRecHitsThreshHB, meNumRecHitsThreshHE, meNumRecHitsThreshHO, meRecHitsEnergyHB, meRecHitsEnergyHE, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitSimHitHB, meRecHitSimHitHE, meRecHitSimHitHF, meRecHitSimHitHFL, meRecHitSimHitHFS, meRecHitSimHitHO, meRecHitSimHitProfileHB, meRecHitSimHitProfileHE, meRecHitSimHitProfileHF, meRecHitSimHitProfileHFL, meRecHitSimHitProfileHFS, meRecHitSimHitProfileHO, meSumRecHitsEnergyConeHB, meSumRecHitsEnergyConeHE, meSumRecHitsEnergyConeHF, meSumRecHitsEnergyConeHFL, meSumRecHitsEnergyConeHFS, meSumRecHitsEnergyConeHO, meSumRecHitsEnergyHB, meSumRecHitsEnergyHE, meSumRecHitsEnergyHF, meSumRecHitsEnergyHO, meTE_HB, meTE_HB1, meTE_HB2, meTE_HE, meTE_HE1, meTE_HE2, meTE_HF, meTE_HFL, meTE_HFS, 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, Nhb, Nhe, Nhf, Nho, occupancy_map_HB1, occupancy_map_HB2, occupancy_map_HE1, occupancy_map_HE2, occupancy_map_HE3, occupancy_map_HF1, occupancy_map_HF2, occupancy_map_HO, occupancy_seqHB1, occupancy_seqHB2, occupancy_seqHE1, occupancy_seqHE2, occupancy_seqHE3, occupancy_seqHF1, occupancy_seqHF2, occupancy_seqHO, occupancy_vs_ieta_HB1, occupancy_vs_ieta_HB2, occupancy_vs_ieta_HE1, occupancy_vs_ieta_HE2, occupancy_vs_ieta_HE3, occupancy_vs_ieta_HF1, occupancy_vs_ieta_HF2, 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_HF67, RecHit_StatusWord_HO, RMS_seq_HB1, RMS_seq_HB2, RMS_seq_HE1, RMS_seq_HE2, RMS_seq_HE3, RMS_seq_HF1, RMS_seq_HF2, RMS_seq_HO, RMS_vs_ieta_HB1, RMS_vs_ieta_HB2, RMS_vs_ieta_HE1, RMS_vs_ieta_HE2, RMS_vs_ieta_HE3, RMS_vs_ieta_HF1, RMS_vs_ieta_HF2, RMS_vs_ieta_HO, DQMStore::IBooker::setCurrentFolder(), sevLvl_HB, sevLvl_HE, sevLvl_HF, sevLvl_HO, subdet_, useAllHistos_, ZS_HB1, ZS_HB2, ZS_HE1, ZS_HE2, ZS_HE3, ZS_HF1, ZS_HF2, ZS_HO, ZS_nHB1, ZS_nHB2, ZS_nHE1, ZS_nHE2, ZS_nHE3, ZS_nHF1, ZS_nHF2, ZS_nHO, ZS_seqHB1, ZS_seqHB2, ZS_seqHE1, ZS_seqHE2, ZS_seqHE3, ZS_seqHF1, ZS_seqHF2, and ZS_seqHO.

{

  Char_t histo[200];

    ib.setCurrentFolder("HcalRecHitsV/HcalRecHitTask");

    // General counters (drawn)
    sprintf  (histo, "N_HB" );
    Nhb = ib.book1D(histo, histo, 2600,0.,2600.);
    sprintf  (histo, "N_HE" );
    Nhe = ib.book1D(histo, histo, 2600,0.,2600.);
    sprintf  (histo, "N_HO" );
    Nho = ib.book1D(histo, histo, 2200,0.,2200.);
    sprintf  (histo, "N_HF" );
    Nhf = ib.book1D(histo, histo, 1800,0., 1800.);

    // ZS
    if(subdet_ == 6) {

      for (unsigned int i1 = 0;  i1 < 82; i1++) {
    for (unsigned int i2 = 0;  i2 < 72; i2++) {
      for (unsigned int i3 = 0;  i3 < 4;  i3++) {
        for (unsigned int i4 = 0;  i4 < 4;  i4++) {
          emap_min [i1][i2][i3][i4] = 99999.;     
        }
      }
    }
      }

      //None of the ZS histos are drawn
      if (useAllHistos_){
    sprintf  (histo, "ZSmin_map_depth1" );
    map_depth1 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
    sprintf  (histo, "ZSmin_map_depth2" );
    map_depth2 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
    sprintf  (histo, "ZSmin_map_depth3" );
    map_depth3 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
    sprintf  (histo, "ZSmin_map_depth4" );
    map_depth4 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
    
      
    sprintf  (histo, "ZS_Nreco_HB1" );
    ZS_nHB1 = ib.book1D(histo, histo, 2500, 0., 2500.);
    sprintf  (histo, "ZS_Nreco_HB2" );
    ZS_nHB2 = ib.book1D(histo, histo,  500, 0.,  500.);
    sprintf  (histo, "ZS_Nreco_HE1" );
    ZS_nHE1 = ib.book1D(histo, histo, 2000, 0., 2000.);
    sprintf  (histo, "ZS_Nreco_HE2" );
    ZS_nHE2 = ib.book1D(histo, histo, 2000, 0., 2000.);
    sprintf  (histo, "ZS_Nreco_HE3" );
    ZS_nHE3 = ib.book1D(histo, histo,  500, 0.,  500.);
    sprintf  (histo, "ZS_Nreco_HO" );
    ZS_nHO  = ib.book1D(histo, histo, 2500, 0., 2500.);
    sprintf  (histo, "ZS_Nreco_HF1" );
    ZS_nHF1 = ib.book1D(histo, histo, 1000, 0., 1000.);
    sprintf  (histo, "ZS_Nreco_HF2" );
    ZS_nHF2 = ib.book1D(histo, histo, 1000, 0., 1000.);
      
    sprintf  (histo, "ZSmin_simple1D_HB1" );
    ZS_HB1 = ib.book1D(histo, histo,120, -2., 10.);
    sprintf  (histo, "ZSmin_simple1D_HB2" );
    ZS_HB2 = ib.book1D(histo, histo,120, -2., 10.);
    sprintf  (histo, "ZSmin_simple1D_HE1" );
    ZS_HE1 = ib.book1D(histo, histo,120, -2., 10.);
    sprintf  (histo, "ZSmin_simple1D_HE2" );
    ZS_HE2 = ib.book1D(histo, histo,120, -2., 10.);
    sprintf  (histo, "ZSmin_simple1D_HE3" );
    ZS_HE3 = ib.book1D(histo, histo,120, -2., 10.);
    sprintf  (histo, "ZSmin_simple1D_HO" );
    ZS_HO = ib.book1D(histo, histo,120, -2., 10.);
    sprintf  (histo, "ZSmin_simple1D_HF1" );
    ZS_HF1 = ib.book1D(histo, histo,200, -10., 10.);
    sprintf  (histo, "ZSmin_simple1D_HF2" );
    ZS_HF2 = ib.book1D(histo, histo,200, -10., 10.);
    
    sprintf  (histo, "ZSmin_sequential1D_HB1" );
    ZS_seqHB1 = ib.book1D(histo, histo,2400, -1200., 1200.);
    sprintf  (histo, "ZSmin_sequential1D_HB2" );
    ZS_seqHB2 = ib.book1D(histo, histo,2400, -1200., 1200.);
    sprintf  (histo, "ZSmin_sequential1D_HE1" );
    ZS_seqHE1 = ib.book1D(histo, histo,4400, -2200., 2200.);
    sprintf  (histo, "ZSmin_sequential1D_HE2" );
    ZS_seqHE2 = ib.book1D(histo, histo,4400, -2200., 2200.);
    sprintf  (histo, "ZSmin_sequential1D_HE3" );
    ZS_seqHE3 = ib.book1D(histo, histo,4400, -2200., 2200.);
    sprintf  (histo, "ZSmin_sequential1D_HO" );
    ZS_seqHO  = ib.book1D(histo, histo,2400, -1200., 1200.);
    sprintf  (histo, "ZSmin_sequential1D_HF1" );
    ZS_seqHF1 = ib.book1D(histo, histo,6000, -3000., 3000.);
    sprintf  (histo, "ZSmin_sequential1D_HF2" );
    ZS_seqHF2 = ib.book1D(histo, histo,6000, -3000., 3000.);
      }
    }

    // ALL others, except ZS
    else {
  
      sprintf  (histo, "emap_depth1" );
      emap_depth1 = ib.book2D(histo, histo, 84, -42., 42., 72, 0., 72.);
      sprintf  (histo, "emap_depth2" );
      emap_depth2 = ib.book2D(histo, histo, 84, -42., 42., 72, 0., 72.);
      sprintf  (histo, "emap_depth3" );
      emap_depth3 = ib.book2D(histo, histo, 84, -42., 42., 72, 0., 72.);
      sprintf  (histo, "emap_depth4" );
      emap_depth4 = ib.book2D(histo, histo, 84, -42., 42., 72, 0., 72.);
      
      if (useAllHistos_){
    
    if (ecalselector_ == "yes") {
      sprintf  (histo, "map_ecal" );
      map_ecal = ib.book2D(histo, histo, 70, -3.045, 3.045, 72, -3.1415926536, 3.1415926536);
    }
      }
      
      //The mean energy histos are drawn, but not the RMS or emean seq
      sprintf  (histo, "emean_vs_ieta_HB1" );
      emean_vs_ieta_HB1 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s");
      sprintf  (histo, "emean_vs_ieta_HB2" );
      emean_vs_ieta_HB2 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s");
      sprintf  (histo, "emean_vs_ieta_HE1" );
      emean_vs_ieta_HE1 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10. ,2000., "s" );
      sprintf  (histo, "emean_vs_ieta_HE2" );
      emean_vs_ieta_HE2 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s");
      sprintf  (histo, "emean_vs_ieta_HE3" );
      emean_vs_ieta_HE3 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s" );
      sprintf  (histo, "emean_vs_ieta_HO" );
      emean_vs_ieta_HO = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s" );
      sprintf  (histo, "emean_vs_ieta_HF1" );
      emean_vs_ieta_HF1 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s" );
      sprintf  (histo, "emean_vs_ieta_HF2" );
      emean_vs_ieta_HF2 = ib.bookProfile(histo, histo, 82, -41., 41., 2010, -10., 2000., "s" );

      if (useAllHistos_){
    sprintf  (histo, "RMS_vs_ieta_HB1" );
    RMS_vs_ieta_HB1 = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HB2" );
    RMS_vs_ieta_HB2 = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HE1" );
    RMS_vs_ieta_HE1 = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HE2" );
    RMS_vs_ieta_HE2 = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HE3" );
    RMS_vs_ieta_HE3 = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HO" );
    RMS_vs_ieta_HO = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HF1" );
    RMS_vs_ieta_HF1 = ib.book1D(histo, histo, 82, -41., 41.);
    sprintf  (histo, "RMS_vs_ieta_HF2" );
    RMS_vs_ieta_HF2 = ib.book1D(histo, histo, 82, -41., 41.);
    
    // Sequential emean and RMS
    sprintf  (histo, "emean_seq_HB1" );
    emean_seqHB1 = ib.bookProfile(histo, histo, 2400, -1200., 1200.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HB2" );
    emean_seqHB2 = ib.bookProfile(histo, histo, 2400, -1200., 1200.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HE1" );
    emean_seqHE1 = ib.bookProfile(histo, histo, 4400, -2200., 2200.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HE2" );
    emean_seqHE2 = ib.bookProfile(histo, histo, 4400, -2200., 2200.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HE3" );
    emean_seqHE3 = ib.bookProfile(histo, histo, 4400, -2200., 2200.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HO" );
    emean_seqHO = ib.bookProfile(histo, histo,  2400, -1200., 1200.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HF1" );
    emean_seqHF1 = ib.bookProfile(histo, histo, 6000, -3000., 3000.,  2010, -10., 2000., "s" );
    sprintf  (histo, "emean_seq_HF2" );
    emean_seqHF2 = ib.bookProfile(histo, histo, 6000, -3000., 3000.,  2010, -10., 2000., "s" );
    
    sprintf  (histo, "RMS_seq_HB1" );
    RMS_seq_HB1 = ib.book1D(histo, histo, 2400, -1200., 1200.);
    sprintf  (histo, "RMS_seq_HB2" );
    RMS_seq_HB2 = ib.book1D(histo, histo, 2400, -1200., 1200.);
    sprintf  (histo, "RMS_seq_HE1" );
    RMS_seq_HE1 = ib.book1D(histo, histo, 4400, -2200., 2200.);
    sprintf  (histo, "RMS_seq_HE2" );
    RMS_seq_HE2 = ib.book1D(histo, histo, 4400, -2200., 2200.);
    sprintf  (histo, "RMS_seq_HE3" );
    RMS_seq_HE3 = ib.book1D(histo, histo, 4400, -2200., 2200.);
    sprintf  (histo, "RMS_seq_HO" );
    RMS_seq_HO = ib.book1D(histo, histo, 2400, -1200., 1200.);
    sprintf  (histo, "RMS_seq_HF1" );
    RMS_seq_HF1 = ib.book1D(histo, histo, 6000, -3000., 3000.);
    sprintf  (histo, "RMS_seq_HF2" );
    RMS_seq_HF2 = ib.book1D(histo, histo, 6000, -3000., 3000.);
      }
      // Occupancy
      //The only occupancy histos drawn are occupancy vs. ieta
      //but the maps are needed because this is where the latter are filled from
      sprintf  (histo, "occupancy_map_HB1" );
      occupancy_map_HB1 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
      sprintf  (histo, "occupancy_map_HB2" );
      occupancy_map_HB2 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
      sprintf  (histo, "occupancy_map_HE1" );
      occupancy_map_HE1 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
      sprintf  (histo, "occupancy_map_HE2" );
      occupancy_map_HE2 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);      
      sprintf  (histo, "occupancy_map_HE3" );
      occupancy_map_HE3 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
      sprintf  (histo, "occupancy_map_HO" );
      occupancy_map_HO = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);      
      sprintf  (histo, "occupancy_map_HF1" );
      occupancy_map_HF1 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
      sprintf  (histo, "occupancy_map_HF2" );
      occupancy_map_HF2 = ib.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
      
      //These are drawn
      sprintf  (histo, "occupancy_vs_ieta_HB1" );
      occupancy_vs_ieta_HB1 = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HB2" );
      occupancy_vs_ieta_HB2 = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HE1" );
      occupancy_vs_ieta_HE1 = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HE2" );
      occupancy_vs_ieta_HE2 = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HE3" );
      occupancy_vs_ieta_HE3 = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HO" );
      occupancy_vs_ieta_HO = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HF1" );
      occupancy_vs_ieta_HF1 = ib.book1D(histo, histo, 82, -41., 41.);
      sprintf  (histo, "occupancy_vs_ieta_HF2" );
      occupancy_vs_ieta_HF2 = ib.book1D(histo, histo, 82, -41., 41.);
      
      //These are not
      if (useAllHistos_){
    sprintf  (histo, "occ_sequential1D_HB1" );
    occupancy_seqHB1 = ib.book1D(histo, histo,2400, -1200., 1200.);
    sprintf  (histo, "occ_sequential1D_HB2" );
    occupancy_seqHB2 = ib.book1D(histo, histo,2400, -1200., 1200.);
    sprintf  (histo, "occ_sequential1D_HE1" );
    occupancy_seqHE1 = ib.book1D(histo, histo,4400, -2200., 2200.);
    sprintf  (histo, "occ_sequential1D_HE2" );
    occupancy_seqHE2 = ib.book1D(histo, histo,4400, -2200., 2200.);
    sprintf  (histo, "occ_sequential1D_HE3" );
    occupancy_seqHE3 = ib.book1D(histo, histo,4400, -2200., 2200.);
    sprintf  (histo, "occ_sequential1D_HO" );
    occupancy_seqHO  = ib.book1D(histo, histo,2400, -1200., 1200.);
    sprintf  (histo, "occ_sequential1D_HF1" );
    occupancy_seqHF1 = ib.book1D(histo, histo,6000, -3000., 3000.);
    sprintf  (histo, "occ_sequential1D_HF2" );
    occupancy_seqHF2 = ib.book1D(histo, histo,6000, -3000., 3000.);
      }

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

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

      if (useAllHistos_){
    sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HF67" ) ;
    RecHit_StatusWord_HF67 = ib.book1D(histo, histo, 3 , 0.5, 3.5); 
      }
      sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HO" ) ;
      RecHit_StatusWord_HO = ib.book1D(histo, histo, 32 , -0.5, 31.5); 

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

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

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

      //These are not drawn
      if(imc !=0 && useAllHistos_) { 
    sprintf  (histo, "map_econe_depth1" );
    map_econe_depth1 =
      ib.book2D(histo, histo, 520, -5.2, 5.2, 72, -3.1415926536, 3.1415926536);
    sprintf  (histo, "map_econe_depth2" );
    map_econe_depth2 =
      ib.book2D(histo, histo, 520, -5.2, 5.2, 72, -3.1415926536, 3.1415926536);
    sprintf  (histo, "map_econe_depth3" );
    map_econe_depth3 =
      ib.book2D(histo, histo, 520, -5.2, 5.2, 72, -3.1415926536, 3.1415926536);
    sprintf  (histo, "map_econe_depth4" );
    map_econe_depth4 =
      ib.book2D(histo, histo, 520, -5.2, 5.2, 72, -3.1415926536, 3.1415926536);
      }
    }  // end-of (subdet_ =! 6)

    //======================= 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 = ib.bookProfile(histo, histo, 82, -41., 41.,        2100, -100., 2000.);  
      
      sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
      meEnConeEtaProfile_E = ib.bookProfile(histo, histo, 82, -41., 41.,      2100, -100., 2000.);  
      
      sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
      meEnConeEtaProfile_EH = ib.bookProfile(histo, histo, 82, -41., 41.,     2100, -100., 2000.);  
    }
    //The other cone profile, delta ieta/phi and noise histos are not drawn
    if (useAllHistos_){
      if(subdet_ != 0 && imc != 0) { // just not for noise  
    
    //    meEnConeEtaProfiel_depth1->Fill(eta_RecHit, HcalCone_d1);
    
    sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_depth1");
    meEnConeEtaProfile_depth1 = ib.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);   
    
    sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_depth2");
    meEnConeEtaProfile_depth2 = ib.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);  
    
    sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_depth3");
    meEnConeEtaProfile_depth3 = ib.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);  
    
    sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_depth4");
    meEnConeEtaProfile_depth4 = ib.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);  
    
      }
      
      if(etype_ == 1 && subdet_ != 0) { // single part., not for noise
    
    sprintf  (histo, "Delta_phi_cluster-MC");
    meDeltaPhi =  ib.book2D(histo, histo, 520, -5.2, 5.2, 60, -0.6, 0.6);
    
    sprintf  (histo, "Delta_eta_cluster-MC");
    meDeltaEta =  ib.book2D(histo, histo, 520, -5.2, 5.2, 60, -0.6, 0.6);
      
    sprintf  (histo, "Delta_phi_simcluster-MC");
    meDeltaPhiS =  ib.book2D(histo, histo, 520, -5.2, 5.2, 60, -0.6, 0.6);
    
    sprintf  (histo, "Delta_eta_simcluster-MC");
    meDeltaEtaS =  ib.book2D(histo, histo, 520, -5.2, 5.2, 60, -0.6, 0.6);
      }
      // NOISE-specific
      
      if (hcalselector_ == "noise" ){
    
    sprintf  (histo, "e_hb" ) ;
    e_hb = ib.book1D(histo, histo,1000, -5., 5.);
    sprintf  (histo, "e_he" ) ;
    e_he = ib.book1D(histo, histo,1000, -5., 5.);
    sprintf  (histo, "e_ho" ) ;
    e_ho = ib.book1D(histo, histo,1000, -5., 5.);
    sprintf  (histo, "e_hfl" ) ;
    e_hfl = ib.book1D(histo, histo,2000, -10., 10.);
    sprintf  (histo, "e_hfs" ) ;
    e_hfs = ib.book1D(histo, histo,2000, -10., 10.);
      }
    }
    // ************** HB **********************************
    if (subdet_ == 1 || subdet_ == 5 ){

      //Only severity level, energy of rechits and overall HB timing histos are drawn  
      if (useAllHistos_){
    if(etype_ == 1 && subdet_ == 1 ) { 
      if(imc != 0) {
        sprintf (histo, "HcalRecHitTask_number_of_rechits_in_cone_HB" ) ;
        meNumRecHitsConeHB    = ib.book1D(histo, histo, 100, 0., 100.);
        
        sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_in_cone_HB" ) ;
        meSumRecHitsEnergyConeHB = ib.book1D(histo,histo, 60 ,-20., 280.);
      }
      
      sprintf (histo, "HcalRecHitTask_number_of_rechits_above_1GeV_HB");
      meNumRecHitsThreshHB = ib.book1D(histo, histo,  30, 0., 30.); 
      
      sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_HB" ) ;
      meSumRecHitsEnergyHB = ib.book1D(histo,histo, 60 , -20., 280.);
    
      if (ecalselector_ == "yes") {  
        if(imc != 0) {
          sprintf (histo, "HcalRecHitTask_number_of_ecalrechits_in_cone_HB");
          meNumEcalRecHitsConeHB = ib.book1D(histo, histo, 300, 0., 300.);      
          sprintf (histo, "HcalRecHitTask_energy_ecal_plus_hcal_in_cone_HB");
          meEcalHcalEnergyConeHB =  ib.book1D(histo,histo, 60 , -20., 280.);
        }
        
        sprintf (histo, "HcalRecHitTask_energy_hcal_vs_ecal_HB");
        meEnergyHcalVsEcalHB = ib.book2D(histo, histo, 300, 0., 150., 300, 0., 150.);   
        sprintf (histo, "HcalRecHitTask_energy_ecal_plus_hcal_HB" ) ;
        meEcalHcalEnergyHB = ib.book1D(histo,histo, 60 , -20., 280.);
      }
    }
      }
      
      sprintf(histo, "HcalRecHitTask_severityLevel_HB");
      sevLvl_HB = ib.book1D(histo, histo, 25, -0.5, 24.5); 

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HB" ) ;
      meRecHitsEnergyHB = ib.book1D(histo, histo, 2010 , -10. , 2000.); 
      
      sprintf (histo, "HcalRecHitTask_timing_HB" ) ;
      meTimeHB = ib.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 = ib.book2D(histo, histo, 50, -5., 45.,  70, -48., 92.);

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

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

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

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

      //Timing by depth and rechits vs simhits are not drawn
      if (useAllHistos_){
    sprintf (histo, "HcalRecHitTask_timing_vs_energy_HB_depth1" ) ;
    meTE_HB1 = ib.book2D(histo, histo, 3000, -5., 2995.,  70, -48., 92.);
    
    sprintf (histo, "HcalRecHitTask_timing_vs_energy_HB_depth2" ) ;
    meTE_HB2 = ib.book2D(histo, histo, 3000, -5., 2995.,  70, -48., 92.);
    
    if(imc != 0) {
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HB");
      meRecHitSimHitHB = ib.book2D(histo, histo, 120, 0., 1.2,  300, 0., 150.);
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HB");
      meRecHitSimHitProfileHB = ib.bookProfile(histo, histo, 120, 0., 1.2, 500, 0., 500.);  
    }      
      }
    }
    
    // ********************** HE ************************************
    if ( subdet_ == 2 || subdet_ == 5 ){

      //None of these are drawn
      if (useAllHistos_){
    if(etype_ == 1 && subdet_ == 2 ) { 
      
      if(imc != 0) {
        sprintf (histo, "HcalRecHitTask_number_of_rechits_in_cone_HE" ) ;
        meNumRecHitsConeHE    = ib.book1D(histo, histo, 100, 0., 100.);
        
        sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_in_cone_HE" ) ;
        meSumRecHitsEnergyConeHE = ib.book1D(histo,histo, 60 ,-20., 280.);
      } 
      
      sprintf (histo, "HcalRecHitTask_number_of_rechits_above_1GeV_HE");
      meNumRecHitsThreshHE = ib.book1D(histo, histo,  30, 0., 30.);  
      
      sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_HE" ) ;
      meSumRecHitsEnergyHE = ib.book1D(histo,histo, 60 , -20., 280.);
      
      if (ecalselector_ == "yes") {     
        sprintf (histo, "HcalRecHitTask_energy_ecal_plus_hcal_HE" ) ;
        meEcalHcalEnergyHE = ib.book1D(histo,histo, 80, -20., 380.);
        
        sprintf (histo, "HcalRecHitTask_energy_hcal_vs_ecal_HE");
        meEnergyHcalVsEcalHE = ib.book2D(histo, histo, 300, 0., 150., 300, 0., 150.);
        if(imc != 0) {
          sprintf (histo, "HcalRecHitTask_number_of_ecalrechits_in_cone_HE");
          meNumEcalRecHitsConeHE = ib.book1D(histo, histo, 300, 0., 300.);   
          sprintf (histo, "HcalRecHitTask_energy_ecal_plus_hcal_in_cone_HE");
          meEcalHcalEnergyConeHE =  ib.book1D(histo,histo, 60,-20., 280.);
        }
      }       
    }
      }
      
      //Only severity level, energy of rechits and overall HB timing histos are drawn  
      sprintf(histo, "HcalRecHitTask_severityLevel_HE");
      sevLvl_HE = ib.book1D(histo, histo, 25, -0.5, 24.5); 
      
      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HE" ) ;
      meRecHitsEnergyHE = ib.book1D(histo, histo, 510, -10., 500.); 
      
      sprintf (histo, "HcalRecHitTask_timing_HE" ) ;
      meTimeHE = ib.book1D(histo, histo, 70, -48., 92.); 
      
      sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HE" ) ;
      meTE_Low_HE = ib.book2D(histo, histo, 80, -5., 75.,  70, -48., 92.);

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

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HE" ) ;
      meTEprofileHE = ib.bookProfile(histo, histo, 200, -5., 395., 70, -48., 92.); 

      //Timing by depth and rechits vs simhits are not drawn
      if (useAllHistos_){
    sprintf (histo, "HcalRecHitTask_timing_vs_energy_HE_depth1" ) ;
    meTE_HE1 = ib.book2D(histo, histo, 1000, -5., 995., 70, -48., 92.);
    
    sprintf (histo, "HcalRecHitTask_timing_vs_energy_HE_depth2" ) ;
    meTE_HE2 = ib.book2D(histo, histo, 1000, -5., 995.,  70, -48., 92.);
    
    if(imc != 0) {
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HE");
      meRecHitSimHitHE = ib.book2D(histo, histo, 120, 0., 0.6,  300, 0., 150.);
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HE");
      meRecHitSimHitProfileHE = ib.bookProfile(histo, histo, 120, 0., 0.6, 500, 0., 500.);  
    }
      }
      
    }

    // ************** HO ****************************************
    if ( subdet_ == 3 || subdet_ == 5  ){
      
      //Only severity level, energy of rechits and overall HB timing histos are drawn  
      if (useAllHistos_){
    if(etype_ == 1 && subdet_ == 3) { 
      if (imc != 0) {
        sprintf (histo, "HcalRecHitTask_number_of_rechits_in_cone_HO" ) ;
        meNumRecHitsConeHO    = ib.book1D(histo, histo, 100, 0 , 100.);
        
        sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_in_cone_HO" ) ;
        meSumRecHitsEnergyConeHO = ib.book1D(histo,histo, 80 ,-20., 380.);
      }
      
      sprintf (histo, "HcalRecHitTask_number_of_rechits_above_1GeV_HO");
      meNumRecHitsThreshHO = ib.book1D(histo, histo,   100, 0., 100.);   
      
      sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_HO" ) ;
      meSumRecHitsEnergyHO = ib.book1D(histo,histo, 80 , -20., 380.);
    }
      }      
      
      sprintf(histo, "HcalRecHitTask_severityLevel_HO");
      sevLvl_HO = ib.book1D(histo, histo, 25, -0.5, 24.5); 

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

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

      sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HO" ) ;
      meTEprofileHO_High = ib.bookProfile(histo, histo, 100, -5., 995.,  70, -48., 92.); 
      
      //Rechits vs simhits are not drawn
      if (useAllHistos_){
    if(imc != 0) {
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HO");
      meRecHitSimHitHO = ib.book2D(histo, histo, 150, 0., 1.5,  350, 0., 350.);
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HO");
      meRecHitSimHitProfileHO = ib.bookProfile(histo, histo, 150, 0., 1.5, 500, 0., 500.);  
    }
      }
    }   
  
    // ********************** HF ************************************
    if ( subdet_ == 4 || subdet_ == 5 ){

      //Only severity level, energy of rechits and overall HB timing histos are drawn  
      if (useAllHistos_){
    if(etype_ == 1 &&  subdet_ == 4) { 
      
      if(imc != 0) {
        sprintf (histo, "HcalRecHitTask_number_of_rechits_in_cone_HF" ) ;
        meNumRecHitsConeHF    = ib.book1D(histo, histo, 30, 0 , 30.);
        
        sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_in_cone_HF" ) ;
        meSumRecHitsEnergyConeHF = ib.book1D(histo,histo,100, -20., 180.);
        
        sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_in_cone_HFL" );
        meSumRecHitsEnergyConeHFL = ib.book1D(histo,histo,100,-20., 180.);
        
        sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_in_cone_HFS");
        meSumRecHitsEnergyConeHFS = ib.book1D(histo,histo,100,-20., 180.);
      }
      sprintf (histo, "HcalRecHitTask_sum_of_rechits_energy_HF" ) ;
      meSumRecHitsEnergyHF = ib.book1D(histo,histo, 80 , -20., 380.);  
    }
      }
      
      sprintf(histo, "HcalRecHitTask_severityLevel_HF");
      sevLvl_HF = ib.book1D(histo, histo, 25, -0.5, 24.5); 

      sprintf (histo, "HcalRecHitTask_energy_of_rechits_HF" ) ;
      meRecHitsEnergyHF = ib.book1D(histo, histo, 1010 , -10. , 1000.); 

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

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

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

      //Timing by L/S and rechits vs simhits are not drawn
      if (useAllHistos_){
    sprintf (histo, "HcalRecHitTask_timing_vs_energy_HFL" ) ;
    meTE_HFL = ib.book2D(histo, histo, 1000, -5., 995., 70, -48., 92.);
    
    sprintf (histo, "HcalRecHitTask_timing_vs_energy_HFS" ) ;
    meTE_HFS = ib.book2D(histo, histo, 1000, -5., 995., 70, -48., 92.);
    
    if(imc != 0) {
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HF");
      meRecHitSimHitHF  = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);      
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HFL");
      meRecHitSimHitHFL = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);      
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HFS");
      meRecHitSimHitHFS = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HF");
      meRecHitSimHitProfileHF  = ib.bookProfile(histo, histo, 50, 0., 50., 500, 0., 500.);  
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HFL");
      meRecHitSimHitProfileHFL = ib.bookProfile(histo, histo, 50, 0., 50., 500, 0., 500.);  
      sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HFS");
      meRecHitSimHitProfileHFS = ib.bookProfile(histo, histo, 50, 0., 50., 500, 0., 500.);  
    }
      }
    }

}

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

private

Definition at line 1747 of file HcalRecHitsValidation.cc.

References PI, and tmp.

Referenced by analyze().

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

}

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

private

Definition at line 1721 of file HcalRecHitsValidation.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 HcalRecHitsValidation::fillRecHitsTmp ( int  subdet_, edm::Event const &  ev  )

privatevirtual

Definition at line 1561 of file HcalRecHitsValidation.cc.

References cauxstwd, cdepth, cen, ceta, cieta, ciphi, cphi, cstwd, csub, ctime, cz, eta(), geometry, edm::Event::getByToken(), HcalBarrel, HcalEndcap, HcalForward, hcalHBSevLvlVec, hcalHESevLvlVec, hcalHFSevLvlVec, hcalHOSevLvlVec, HcalOuter, hcalSevLvl(), iz, j, phi, edmStreamStallGrapher::t, tok_hbhe_, tok_hf_, and tok_ho_.

Referenced by analyze().

                                                                         {
  
  using namespace edm;
  
  
  // initialize data vectors
  csub.clear();
  cen.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;
    ev.getByToken(tok_hbhe_, hbhecoll);
    
    for (HBHERecHitCollection::const_iterator j=hbhecoll->begin(); j != hbhecoll->end(); j++) {
      HcalDetId cell(j->id());
      const CaloCellGeometry* 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();
      if(inteta > 0) inteta -= 1;
      int intphi  = cell.iphi()-1;
      double en   = j->energy();
      double t    = j->time();
      int stwd    = j->flags();
      int auxstwd = j->aux();
      
      int serivityLevel = hcalSevLvl( (CaloRecHit*) &*j );
      if( cell.subdet()==HcalBarrel ){
         hcalHBSevLvlVec.push_back(serivityLevel);
      }else if (cell.subdet()==HcalEndcap ){
         hcalHESevLvlVec.push_back(serivityLevel);
      } 
      
      if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
    
    csub.push_back(sub);
    cen.push_back(en);
    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;
    ev.getByToken(tok_hf_, hfcoll);

    for (HFRecHitCollection::const_iterator j = hfcoll->begin(); j != hfcoll->end(); j++) {
      HcalDetId cell(j->id());
      const CaloCellGeometry* 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();
      if(inteta > 0) inteta -= 1;
      int intphi   = cell.iphi()-1;
      double en    = j->energy();
      double t     = j->time();
      int stwd     = j->flags();
      int auxstwd  = j->aux();

      int serivityLevel = hcalSevLvl( (CaloRecHit*) &*j );
      if( cell.subdet()==HcalForward ){
         hcalHFSevLvlVec.push_back(serivityLevel);
      } 

      if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
    
    csub.push_back(sub);
    cen.push_back(en);
    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;
    ev.getByToken(tok_ho_, hocoll);
    
    for (HORecHitCollection::const_iterator j = hocoll->begin(); j != hocoll->end(); j++) {
      HcalDetId cell(j->id());
      const CaloCellGeometry* 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();
      if(inteta > 0) inteta -= 1;
      int intphi   = cell.iphi()-1;
      double t     = j->time();
      double en    = j->energy();
      int stwd     = j->flags();
      int auxstwd  = j->aux();

      int serivityLevel = hcalSevLvl( (CaloRecHit*) &*j );
      if( cell.subdet()==HcalOuter ){
         hcalHOSevLvlVec.push_back(serivityLevel);
      } 
      
      if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
    csub.push_back(sub);
    cen.push_back(en);
    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 HcalRecHitsValidation::hcalSevLvl ( const CaloRecHit *  hit )

private

Definition at line 1762 of file HcalRecHitsValidation.cc.

References CaloRecHit::detid(), CaloRecHit::flags(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), theHcalChStatus, and theHcalSevLvlComputer.

Referenced by fillRecHitsTmp().

                                                          {

   const DetId id = hit->detid();

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

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

   return severityLevel;

} 

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

private

Definition at line 1731 of file HcalRecHitsValidation.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> HcalRecHitsValidation::cauxstwd

private

Definition at line 376 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::cdepth

private

Definition at line 369 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsValidation::cen

private

Definition at line 370 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsValidation::ceta

private

Definition at line 371 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::cieta

private

Definition at line 367 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::ciphi

private

Definition at line 368 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsValidation::cphi

private

Definition at line 372 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<uint32_t> HcalRecHitsValidation::cstwd

private

Definition at line 375 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::csub

private

Definition at line 366 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsValidation::ctime

private

Definition at line 373 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsValidation::cz

private

Definition at line 374 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp().

MonitorElement* HcalRecHitsValidation::e_hb

private

Definition at line 247 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::e_he

private

Definition at line 248 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::e_hfl

private

Definition at line 250 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::e_hfs

private

Definition at line 251 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::e_ho

private

Definition at line 249 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

std::string HcalRecHitsValidation::ecalselector_

private

Definition at line 73 of file HcalRecHitsValidation.h.

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

MonitorElement* HcalRecHitsValidation::emap_depth1

private

Definition at line 157 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emap_depth2

private

Definition at line 158 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emap_depth3

private

Definition at line 159 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emap_depth4

private

Definition at line 160 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

double HcalRecHitsValidation::emap_min[82][72][4][4]

private

Definition at line 379 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHB1

private

Definition at line 180 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHB2

private

Definition at line 181 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHE1

private

Definition at line 182 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHE2

private

Definition at line 183 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHE3

private

Definition at line 184 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHF1

private

Definition at line 186 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHF2

private

Definition at line 187 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_seqHO

private

Definition at line 185 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HB1

private

Definition at line 162 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HB2

private

Definition at line 163 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HE1

private

Definition at line 164 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HE2

private

Definition at line 165 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HE3

private

Definition at line 166 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HF1

private

Definition at line 168 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HF2

private

Definition at line 169 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::emean_vs_ieta_HO

private

Definition at line 167 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsValidation::etype_

private

Definition at line 93 of file HcalRecHitsValidation.h.

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

std::string HcalRecHitsValidation::eventype_

private

Definition at line 74 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

bool HcalRecHitsValidation::famos_

private

Definition at line 77 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

edm::ESHandle<CaloGeometry> HcalRecHitsValidation::geometry

private

Definition at line 350 of file HcalRecHitsValidation.h.

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

std::vector<int> HcalRecHitsValidation::hcalHBSevLvlVec

private

Definition at line 103 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::hcalHESevLvlVec

private

Definition at line 103 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::hcalHFSevLvlVec

private

Definition at line 103 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsValidation::hcalHOSevLvlVec

private

Definition at line 103 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

std::string HcalRecHitsValidation::hcalselector_

private

Definition at line 72 of file HcalRecHitsValidation.h.

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

int HcalRecHitsValidation::imc

private

Definition at line 95 of file HcalRecHitsValidation.h.

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

int HcalRecHitsValidation::iz

private

Definition at line 94 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsValidation().

MonitorElement* HcalRecHitsValidation::map_depth1

private

Definition at line 118 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::map_depth2

private

Definition at line 119 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::map_depth3

private

Definition at line 120 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::map_depth4

private

Definition at line 121 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::map_ecal

private

Definition at line 155 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::map_econe_depth1

private

Definition at line 227 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::map_econe_depth2

private

Definition at line 228 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::map_econe_depth3

private

Definition at line 229 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::map_econe_depth4

private

Definition at line 230 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

std::string HcalRecHitsValidation::mc_

private

Definition at line 76 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

MonitorElement* HcalRecHitsValidation::meDeltaEta

private

Definition at line 242 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meDeltaEtaS

private

Definition at line 244 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meDeltaPhi

private

Definition at line 241 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meDeltaPhiS

private

Definition at line 243 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEcalHcalEnergyConeHB

private

Definition at line 322 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEcalHcalEnergyConeHE

private

Definition at line 323 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEcalHcalEnergyConeHF

private

Definition at line 325 of file HcalRecHitsValidation.h.

MonitorElement* HcalRecHitsValidation::meEcalHcalEnergyConeHO

private

Definition at line 324 of file HcalRecHitsValidation.h.

MonitorElement* HcalRecHitsValidation::meEcalHcalEnergyHB

private

Definition at line 319 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEcalHcalEnergyHE

private

Definition at line 320 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile

private

Definition at line 237 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_depth1

private

Definition at line 233 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_depth2

private

Definition at line 234 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_depth3

private

Definition at line 235 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_depth4

private

Definition at line 236 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_E

private

Definition at line 238 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_EH

private

Definition at line 239 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnergyHcalVsEcalHB

private

Definition at line 343 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meEnergyHcalVsEcalHE

private

Definition at line 344 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumEcalRecHitsConeHB

private

Definition at line 347 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumEcalRecHitsConeHE

private

Definition at line 348 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsConeHB

private

Definition at line 259 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsConeHE

private

Definition at line 260 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsConeHF

private

Definition at line 262 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsConeHO

private

Definition at line 261 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsThreshHB

private

Definition at line 254 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsThreshHE

private

Definition at line 255 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meNumRecHitsThreshHO

private

Definition at line 256 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHB

private

Definition at line 271 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHE

private

Definition at line 272 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHF

private

Definition at line 274 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHO

private

Definition at line 273 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHB

private

Definition at line 328 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHE

private

Definition at line 329 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHF

private

Definition at line 331 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHFL

private

Definition at line 332 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHFS

private

Definition at line 333 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHO

private

Definition at line 330 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHB

private

Definition at line 335 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHE

private

Definition at line 336 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHF

private

Definition at line 338 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHFL

private

Definition at line 339 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHFS

private

Definition at line 340 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHO

private

Definition at line 337 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyConeHB

private

Definition at line 311 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyConeHE

private

Definition at line 312 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyConeHF

private

Definition at line 314 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyConeHFL

private

Definition at line 315 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyConeHFS

private

Definition at line 316 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyConeHO

private

Definition at line 313 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyHB

private

Definition at line 305 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyHE

private

Definition at line 306 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyHF

private

Definition at line 308 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meSumRecHitsEnergyHO

private

Definition at line 307 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HB

private

Definition at line 277 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HB1

private

Definition at line 279 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HB2

private

Definition at line 280 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HE

private

Definition at line 286 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HE1

private

Definition at line 287 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HE2

private

Definition at line 288 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HF

private

Definition at line 298 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HFL

private

Definition at line 299 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HFS

private

Definition at line 300 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_High_HB

private

Definition at line 278 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_High_HO

private

Definition at line 293 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_HO

private

Definition at line 292 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_Low_HB

private

Definition at line 276 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_Low_HE

private

Definition at line 285 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTE_Low_HF

private

Definition at line 297 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHB

private

Definition at line 282 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHB_High

private

Definition at line 283 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHB_Low

private

Definition at line 281 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHE

private

Definition at line 290 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHE_Low

private

Definition at line 289 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHF

private

Definition at line 302 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHF_Low

private

Definition at line 301 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHO

private

Definition at line 294 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTEprofileHO_High

private

Definition at line 295 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTimeHB

private

Definition at line 265 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTimeHE

private

Definition at line 266 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTimeHF

private

Definition at line 268 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::meTimeHO

private

Definition at line 267 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsValidation::nevtot

private

Definition at line 382 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

MonitorElement* HcalRecHitsValidation::Nhb

private

Definition at line 111 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::Nhe

private

Definition at line 112 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::Nhf

private

Definition at line 114 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::Nho

private

Definition at line 113 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HB1

private

Definition at line 198 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HB2

private

Definition at line 199 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HE1

private

Definition at line 200 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HE2

private

Definition at line 201 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HE3

private

Definition at line 202 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HF1

private

Definition at line 204 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HF2

private

Definition at line 205 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_map_HO

private

Definition at line 203 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHB1

private

Definition at line 216 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHB2

private

Definition at line 217 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHE1

private

Definition at line 218 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHE2

private

Definition at line 219 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHE3

private

Definition at line 220 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHF1

private

Definition at line 222 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHF2

private

Definition at line 223 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_seqHO

private

Definition at line 221 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HB1

private

Definition at line 207 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HB2

private

Definition at line 208 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HE1

private

Definition at line 209 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HE2

private

Definition at line 210 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HE3

private

Definition at line 211 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HF1

private

Definition at line 213 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HF2

private

Definition at line 214 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::occupancy_vs_ieta_HO

private

Definition at line 212 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

std::string HcalRecHitsValidation::outputFile_

private

Definition at line 71 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

MonitorElement* HcalRecHitsValidation::RecHit_Aux_StatusWord_HB

private

Definition at line 360 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_Aux_StatusWord_HE

private

Definition at line 361 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_Aux_StatusWord_HF

private

Definition at line 362 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_Aux_StatusWord_HO

private

Definition at line 363 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_StatusWord_HB

private

Definition at line 353 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_StatusWord_HE

private

Definition at line 354 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_StatusWord_HF

private

Definition at line 355 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_StatusWord_HF67

private

Definition at line 356 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RecHit_StatusWord_HO

private

Definition at line 357 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HB1

private

Definition at line 189 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HB2

private

Definition at line 190 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HE1

private

Definition at line 191 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HE2

private

Definition at line 192 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HE3

private

Definition at line 193 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HF1

private

Definition at line 195 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HF2

private

Definition at line 196 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_seq_HO

private

Definition at line 194 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HB1

private

Definition at line 171 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HB2

private

Definition at line 172 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HE1

private

Definition at line 173 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HE2

private

Definition at line 174 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HE3

private

Definition at line 175 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HF1

private

Definition at line 177 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HF2

private

Definition at line 178 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::RMS_vs_ieta_HO

private

Definition at line 176 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::sevLvl_HB

private

Definition at line 105 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::sevLvl_HE

private

Definition at line 106 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::sevLvl_HF

private

Definition at line 107 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::sevLvl_HO

private

Definition at line 108 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

std::string HcalRecHitsValidation::sign_

private

Definition at line 75 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

int HcalRecHitsValidation::subdet_

private

Definition at line 90 of file HcalRecHitsValidation.h.

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

const HcalChannelQuality* HcalRecHitsValidation::theHcalChStatus

private

Definition at line 98 of file HcalRecHitsValidation.h.

Referenced by analyze(), and hcalSevLvl().

const HcalSeverityLevelComputer* HcalRecHitsValidation::theHcalSevLvlComputer

private

Definition at line 100 of file HcalRecHitsValidation.h.

Referenced by analyze(), and hcalSevLvl().

edm::EDGetTokenT<EBRecHitCollection> HcalRecHitsValidation::tok_EB_

private

Definition at line 82 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

edm::EDGetTokenT<EERecHitCollection> HcalRecHitsValidation::tok_EE_

private

Definition at line 83 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

edm::EDGetTokenT<edm::HepMCProduct> HcalRecHitsValidation::tok_evt_

private

Definition at line 81 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

edm::EDGetTokenT<HBHERecHitCollection> HcalRecHitsValidation::tok_hbhe_

private

Definition at line 85 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsValidation().

edm::EDGetTokenT<HFRecHitCollection> HcalRecHitsValidation::tok_hf_

private

Definition at line 86 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsValidation().

edm::EDGetTokenT<edm::PCaloHitContainer> HcalRecHitsValidation::tok_hh_

private

Definition at line 84 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

edm::EDGetTokenT<HORecHitCollection> HcalRecHitsValidation::tok_ho_

private

Definition at line 87 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsValidation().

bool HcalRecHitsValidation::useAllHistos_

private

Definition at line 78 of file HcalRecHitsValidation.h.

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

MonitorElement* HcalRecHitsValidation::ZS_HB1

private

Definition at line 123 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HB2

private

Definition at line 124 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HE1

private

Definition at line 125 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HE2

private

Definition at line 126 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HE3

private

Definition at line 127 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HF1

private

Definition at line 129 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HF2

private

Definition at line 130 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_HO

private

Definition at line 128 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHB1

private

Definition at line 132 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHB2

private

Definition at line 133 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHE1

private

Definition at line 134 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHE2

private

Definition at line 135 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHE3

private

Definition at line 136 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHF1

private

Definition at line 138 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHF2

private

Definition at line 139 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_nHO

private

Definition at line 137 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHB1

private

Definition at line 141 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHB2

private

Definition at line 142 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHE1

private

Definition at line 143 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHE2

private

Definition at line 144 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHE3

private

Definition at line 145 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHF1

private

Definition at line 147 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHF2

private

Definition at line 148 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsValidation::ZS_seqHO

private

Definition at line 146 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().