#include <HcalRecHitsAnalyzer.h>

Inheritance diagram for HcalRecHitsAnalyzer:

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

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

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

	HcalRecHitsAnalyzer (edm::ParameterSet const &conf)

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

uint32_t	streamId () const

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

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

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

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

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

int	hcalSevLvl (const CaloRecHit *hit)

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

Private Attributes
std::vector< uint32_t >	cauxstwd

std::vector< double >	cchi2

std::vector< int >	cdepth

std::vector< double >	cen

std::vector< double >	cenM0

std::vector< double >	cenM3

std::vector< double >	ceta

std::vector< int >	cieta

std::vector< int >	ciphi

std::vector< double >	cphi

std::vector< uint32_t >	cstwd

std::vector< int >	csub

std::vector< double >	ctime

std::vector< double >	cz

std::string	ecalselector_

std::vector< MonitorElement * >	emap

std::vector< MonitorElement * >	emean_vs_ieta_HB

std::vector< MonitorElement * >	emean_vs_ieta_HBM0

std::vector< MonitorElement * >	emean_vs_ieta_HBM3

std::vector< MonitorElement * >	emean_vs_ieta_HE

std::vector< MonitorElement * >	emean_vs_ieta_HEM0

std::vector< MonitorElement * >	emean_vs_ieta_HEM3

std::vector< MonitorElement * >	emean_vs_ieta_HEP17

std::vector< MonitorElement * >	emean_vs_ieta_HEP17M0

std::vector< MonitorElement * >	emean_vs_ieta_HEP17M3

std::vector< MonitorElement * >	emean_vs_ieta_HF

MonitorElement *	emean_vs_ieta_HO

int	etype_

std::string	eventype_

bool	famos_

edm::ESHandle< CaloGeometry >	geometry

std::vector< int >	hcalHBSevLvlVec

std::vector< int >	hcalHESevLvlVec

std::vector< int >	hcalHFSevLvlVec

std::vector< int >	hcalHOSevLvlVec

std::string	hcalselector_

const HcalDDDRecConstants *	hcons

bool	hep17_

int	ieta_bins_

float	ieta_max_

float	ieta_min_

int	imc

int	iphi_bins_

float	iphi_max_

float	iphi_min_

int	iz

MonitorElement *	map_ecal

int	maxDepthAll_

int	maxDepthHB_

int	maxDepthHE_

int	maxDepthHF_

int	maxDepthHO_

std::string	mc_

MonitorElement *	meDeltaEta

MonitorElement *	meDeltaPhi

MonitorElement *	meEcalHcalEnergyConeHB

MonitorElement *	meEcalHcalEnergyConeHE

MonitorElement *	meEcalHcalEnergyConeHF

MonitorElement *	meEcalHcalEnergyConeHO

MonitorElement *	meEcalHcalEnergyHB

MonitorElement *	meEcalHcalEnergyHE

MonitorElement *	meEnConeEtaProfile

MonitorElement *	meEnConeEtaProfile_E

MonitorElement *	meEnConeEtaProfile_EH

MonitorElement *	meEnergyHcalVsEcalHB

MonitorElement *	meEnergyHcalVsEcalHE

MonitorElement *	meLog10Chi2profileHB

MonitorElement *	meLog10Chi2profileHE

MonitorElement *	meNumEcalRecHitsConeHB

MonitorElement *	meNumEcalRecHitsConeHE

MonitorElement *	meRecHitsEnergyHB

MonitorElement *	meRecHitsEnergyHBM0

MonitorElement *	meRecHitsEnergyHBM3

MonitorElement *	meRecHitsEnergyHE

MonitorElement *	meRecHitsEnergyHEM0

MonitorElement *	meRecHitsEnergyHEM3

std::vector< MonitorElement * >	meRecHitsEnergyHEP17

std::vector< MonitorElement * >	meRecHitsEnergyHEP17M0

std::vector< MonitorElement * >	meRecHitsEnergyHEP17M3

MonitorElement *	meRecHitsEnergyHF

MonitorElement *	meRecHitsEnergyHO

MonitorElement *	meRecHitsEnergyM2vM0HB

MonitorElement *	meRecHitsEnergyM2vM0HE

MonitorElement *	meRecHitsEnergyM3vM0HB

MonitorElement *	meRecHitsEnergyM3vM0HE

MonitorElement *	meRecHitsEnergyM3vM2HB

MonitorElement *	meRecHitsEnergyM3vM2HE

MonitorElement *	meRecHitsM2Chi2HB

MonitorElement *	meRecHitsM2Chi2HE

MonitorElement *	meSumRecHitsEnergyConeHB

MonitorElement *	meSumRecHitsEnergyConeHE

MonitorElement *	meSumRecHitsEnergyConeHF

MonitorElement *	meSumRecHitsEnergyConeHFL

MonitorElement *	meSumRecHitsEnergyConeHFS

MonitorElement *	meSumRecHitsEnergyConeHO

MonitorElement *	meSumRecHitsEnergyHB

MonitorElement *	meSumRecHitsEnergyHE

MonitorElement *	meSumRecHitsEnergyHF

MonitorElement *	meSumRecHitsEnergyHO

MonitorElement *	meTE_HB

MonitorElement *	meTE_HE

MonitorElement *	meTE_HF

MonitorElement *	meTE_High_HB

MonitorElement *	meTE_High_HO

MonitorElement *	meTE_HO

MonitorElement *	meTE_Low_HB

MonitorElement *	meTE_Low_HE

MonitorElement *	meTE_Low_HF

MonitorElement *	meTEprofileHB

MonitorElement *	meTEprofileHB_High

MonitorElement *	meTEprofileHB_Low

MonitorElement *	meTEprofileHE

MonitorElement *	meTEprofileHE_Low

MonitorElement *	meTEprofileHF

MonitorElement *	meTEprofileHF_Low

MonitorElement *	meTEprofileHO

MonitorElement *	meTEprofileHO_High

MonitorElement *	meTimeHB

MonitorElement *	meTimeHE

MonitorElement *	meTimeHF

MonitorElement *	meTimeHO

int	nChannels_ [5]

int	nevtot

std::vector< MonitorElement * >	Nhb

std::vector< MonitorElement * >	Nhe

std::vector< MonitorElement * >	Nhf

std::vector< MonitorElement * >	Nho

std::vector< MonitorElement * >	nrechits_vs_iphi_HBM

std::vector< MonitorElement * >	nrechits_vs_iphi_HBP

std::vector< MonitorElement * >	nrechits_vs_iphi_HEM

std::vector< MonitorElement * >	nrechits_vs_iphi_HEP

std::vector< MonitorElement * >	nrechits_vs_iphi_HFM

std::vector< MonitorElement * >	nrechits_vs_iphi_HFP

MonitorElement *	nrechits_vs_iphi_HOM

MonitorElement *	nrechits_vs_iphi_HOP

std::vector< MonitorElement * >	occupancy_map_HB

std::vector< MonitorElement * >	occupancy_map_HE

std::vector< MonitorElement * >	occupancy_map_HF

MonitorElement *	occupancy_map_HO

std::vector< MonitorElement * >	occupancy_vs_ieta_HB

std::vector< MonitorElement * >	occupancy_vs_ieta_HE

std::vector< MonitorElement * >	occupancy_vs_ieta_HF

MonitorElement *	occupancy_vs_ieta_HO

std::string	outputFile_

MonitorElement *	RecHit_Aux_StatusWord_HB

MonitorElement *	RecHit_Aux_StatusWord_HE

MonitorElement *	RecHit_Aux_StatusWord_HF

MonitorElement *	RecHit_Aux_StatusWord_HO

MonitorElement *	RecHit_StatusWord_HB

MonitorElement *	RecHit_StatusWord_HE

MonitorElement *	RecHit_StatusWord_HF

MonitorElement *	RecHit_StatusWord_HF67

MonitorElement *	RecHit_StatusWord_HO

MonitorElement *	RecHit_StatusWordCorr_HB

MonitorElement *	RecHit_StatusWordCorr_HE

MonitorElement *	sevLvl_HB

MonitorElement *	sevLvl_HE

MonitorElement *	sevLvl_HF

MonitorElement *	sevLvl_HO

std::string	sign_

int	subdet_

const HcalChannelQuality *	theHcalChStatus

const HcalSeverityLevelComputer *	theHcalSevLvlComputer

const HcalTopology *	theHcalTopology

edm::EDGetTokenT< EBRecHitCollection >	tok_EB_

edm::EDGetTokenT< EERecHitCollection >	tok_EE_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< HFRecHitCollection >	tok_hf_

edm::EDGetTokenT< HORecHitCollection >	tok_ho_

std::string	topFolderName_

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

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

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

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

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions 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 58 of file HcalRecHitsAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 6 of file HcalRecHitsAnalyzer.cc.

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

   : topFolderName_           ( conf.getParameter<std::string>("TopFolderName") )
 {
 
   // DQM ROOT output
   outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
   
   if ( !outputFile_.empty() ) {
     edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
   } else {
     edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
   }
   
   nevtot = 0;
  
   hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");
   ecalselector_ = conf.getUntrackedParameter<std::string>("ecalselector", "yes");
   eventype_     = conf.getUntrackedParameter<std::string>("eventype", "single");
   sign_         = conf.getUntrackedParameter<std::string>("sign", "*");
   //useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);
 
   //HEP17 configuration
   hep17_        = conf.getUntrackedParameter<bool>("hep17");
 
   //Collections
   tok_hbhe_ = consumes<HBHERecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HBHERecHitCollectionLabel"));
   tok_hf_  = consumes<HFRecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HFRecHitCollectionLabel"));
   tok_ho_ = consumes<HORecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HORecHitCollectionLabel"));
   edm::InputTag EBRecHitCollectionLabel = conf.getParameter<edm::InputTag>("EBRecHitCollectionLabel");
   tok_EB_ = consumes<EBRecHitCollection>(EBRecHitCollectionLabel);
   edm::InputTag EERecHitCollectionLabel = conf.getParameter<edm::InputTag>("EERecHitCollectionLabel");
   tok_EE_ = consumes<EERecHitCollection>(EERecHitCollectionLabel);
 
   subdet_ = 5;
   if (hcalselector_ == "noise") subdet_ = 0;
   if (hcalselector_ == "HB"   ) subdet_ = 1;
   if (hcalselector_ == "HE"   ) subdet_ = 2;
   if (hcalselector_ == "HO"   ) subdet_ = 3;
   if (hcalselector_ == "HF"   ) subdet_ = 4;
   if (hcalselector_ == "all"  ) subdet_ = 5;
   if (hcalselector_ == "ZS"   ) subdet_ = 6;
 
   etype_ = 1;
   if (eventype_ == "multi") etype_ = 2;
 
   iz = 1;
   if(sign_ == "-") iz = -1;
   if(sign_ == "*") iz = 0;
 
   imc = 0;
 
   }

Member Function Documentation

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

overridevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 505 of file HcalRecHitsAnalyzer.cc.

                                                                             {
 
   using namespace edm;
 
   // cuts for each subdet_ector mimiking  "Scheme B"
   //  double cutHB = 0.9, cutHE = 1.4, cutHO = 1.1, cutHFL = 1.2, cutHFS = 1.8; 
 
   // energy in HCAL
   double eHcal        = 0.;
   // Total numbet of RecHits in HCAL, in the cone, above 1 GeV theshold
   int nrechits       = 0;
   int nrechitsThresh = 0;
 
   // energy in ECAL
   double eEcal       = 0.;
   double eEcalB      = 0.;
   double eEcalE      = 0.;
   double eEcalCone   = 0.;
 
   // HCAL energy around MC eta-phi at all depths;
   double partR = 0.3;
 
   // Single particle samples: actual eta-phi position of cluster around
   // hottest cell
   double etaHot  = 99999.; 
   double phiHot  = 99999.; 
 
   //   previously was:  c.get<IdealGeometryRecord>().get (geometry);
   c.get<CaloGeometryRecord>().get (geometry);
 
   // HCAL Topology **************************************************
   edm::ESHandle<HcalTopology> topo;
   c.get<HcalRecNumberingRecord>().get(topo);
   theHcalTopology = topo.product();
 
   // HCAL channel status map ****************************************
   edm::ESHandle<HcalChannelQuality> hcalChStatus;
   c.get<HcalChannelQualityRcd>().get( "withTopo", hcalChStatus );
   theHcalChStatus = hcalChStatus.product();
 
   // Assignment of severity levels **********************************
   edm::ESHandle<HcalSeverityLevelComputer> hcalSevLvlComputerHndl;
   c.get<HcalSeverityLevelComputerRcd>().get(hcalSevLvlComputerHndl);
   theHcalSevLvlComputer = hcalSevLvlComputerHndl.product(); 
 
   // Fill working vectors of HCAL RecHits quantities (all of these are drawn)
   fillRecHitsTmp(subdet_, ev); 
 
   // HB   
   if( subdet_ ==5 || subdet_ == 1 ){ 
      for(unsigned int iv=0; iv<hcalHBSevLvlVec.size(); iv++){
         sevLvl_HB->Fill(hcalHBSevLvlVec[iv]);
      }    
   }
   // HE   
   if( subdet_ ==5 || subdet_ == 2 ){
      for(unsigned int iv=0; iv<hcalHESevLvlVec.size(); iv++){
         sevLvl_HE->Fill(hcalHESevLvlVec[iv]);
      }
   }
   // HO 
   if( subdet_ ==5 || subdet_ == 3 ){
      for(unsigned int iv=0; iv<hcalHOSevLvlVec.size(); iv++){
         sevLvl_HO->Fill(hcalHOSevLvlVec[iv]);
      }
   }
   // HF 
   if( subdet_ ==5 || subdet_ == 4 ){
      for(unsigned int iv=0; iv<hcalHFSevLvlVec.size(); iv++){
         sevLvl_HF->Fill(hcalHFSevLvlVec[iv]);
      }
   } 
 
   //===========================================================================
   // IN ALL other CASES : ieta-iphi maps 
   //===========================================================================
 
   // ECAL 
   if(ecalselector_ == "yes" && (subdet_ == 1 || subdet_ == 2 || subdet_ == 5)) {
 
     Handle<EBRecHitCollection> rhitEB;
     if ( ev.getByToken(tok_EB_, rhitEB) ) {
       
       for ( const auto & recHit : *(rhitEB.product()) ) {
     
     double en  = recHit.energy();
     eEcal  += en;
     eEcalB += en;
       }      
     }
   
     
     Handle<EERecHitCollection> rhitEE;
     if ( ev.getByToken(tok_EE_, rhitEE) ) {
     
       for ( const auto & recHit : *(rhitEE.product()) ) {
     
     double en   = recHit.energy();
     eEcal  += en;
     eEcalE += en;
     
       }
     }
   }     // end of ECAL selection 
 
   // Counting, including ZS items
   // Filling HCAL maps  ----------------------------------------------------
   //   double maxE = -99999.;
   
   // element 0: any depth. element 1,2,..: depth 1,2
   std::vector<int> nhb_v(maxDepthHB_+1,0);
   std::vector<int> nhe_v(maxDepthHE_+1,0);
   std::vector<int> nho_v(maxDepthHO_+1,0);
   std::vector<int> nhf_v(maxDepthHF_+1,0);
 
   for (unsigned int i = 0; i < cen.size(); i++) {
     
     int sub       = csub[i];
     int depth     = cdepth[i];
     int ieta      = cieta[i]; 
     int iphi      = ciphi[i]; 
     double en     = cen[i]; 
     double enM0   = cenM0[i];
     double enM3   = cenM3[i];
     //     double eta    = ceta[i]; 
     //     double phi    = cphi[i]; 
     uint32_t stwd = cstwd[i];
     uint32_t auxstwd = cauxstwd[i];
     //    double z   = cz[i];
 
     //This will be true if hep17 == "yes" and the rechit is in the hep17 wedge
     bool isHEP17 = (sub == 2) && (iphi >= 63) && (iphi <= 66) && (ieta > 0) && (hep17_);
 
     //Make sure that an invalid depth won't cause an error. We should probably report the problem as well.
     if( depth < 1 ) continue;
     if( sub == 1 && depth > maxDepthHB_ ) continue;
     if( sub == 2 && depth > maxDepthHE_ ) continue;
     if( sub == 3 && depth > maxDepthHO_ ) continue;
     if( sub == 4 && depth > maxDepthHF_ ) continue;
 
     if( sub ==1 ){ nhb_v[depth]++; nhb_v[0]++;} // element 0: any depth, element 1,2,..: depth 1,2,...
     if( sub ==2 ){ nhe_v[depth]++; nhe_v[0]++;} //
     if( sub ==3 ){ nho_v[depth]++; nho_v[0]++;} //
     if( sub ==4 ){ nhf_v[depth]++; nhf_v[0]++;} //
 
     if( subdet_ == 6) {                                    // ZS specific
     }
 
     if( subdet_ != 6) {  
       int ieta2 = ieta;
       int depth2 = depth;
       if(sub == 4){
     if (ieta2 < 0) ieta2--;
         else ieta2++;
       }
       if(sub == 3) depth2 = maxDepthAll_ - maxDepthHO_ + depth; //This will use the last depths for HO  
       emap[depth2-1]->Fill(double(ieta2),double(iphi),en);
 
       // to distinguish HE and HF
       if( depth == 1 || depth == 2 ) {
         int ieta1 =  ieta;
     if(sub == 4) { 
       if (ieta1 < 0) ieta1--;
           else  ieta1++;   
     }
       }
 
       if ( sub == 1){
      emean_vs_ieta_HB[depth-1]->Fill(double(ieta), en);
      emean_vs_ieta_HBM0[depth-1]->Fill(double(ieta), enM0);
      emean_vs_ieta_HBM3[depth-1]->Fill(double(ieta), enM3);
      occupancy_map_HB[depth-1]->Fill(double(ieta),double(iphi));
      if (ieta>0) nrechits_vs_iphi_HBP[depth-1]->Fill(double(iphi));
      else        nrechits_vs_iphi_HBM[depth-1]->Fill(double(iphi));
       }
       if ( sub == 2){
      if(!isHEP17){
         emean_vs_ieta_HE[depth-1]->Fill(double(ieta), en);
         emean_vs_ieta_HEM0[depth-1]->Fill(double(ieta), enM0);
         emean_vs_ieta_HEM3[depth-1]->Fill(double(ieta), enM3);
          }else{
         emean_vs_ieta_HEP17[depth-1]->Fill(double(ieta), en);
         emean_vs_ieta_HEP17M0[depth-1]->Fill(double(ieta), enM0);
         emean_vs_ieta_HEP17M3[depth-1]->Fill(double(ieta), enM3);
          }
      occupancy_map_HE[depth-1]->Fill(double(ieta),double(iphi));
      if (ieta>0) nrechits_vs_iphi_HEP[depth-1]->Fill(double(iphi));
      else        nrechits_vs_iphi_HEM[depth-1]->Fill(double(iphi));
       }
       if ( sub == 3){
      emean_vs_ieta_HO->Fill(double(ieta), en);
      occupancy_map_HO->Fill(double(ieta),double(iphi));
      if (ieta>0) nrechits_vs_iphi_HOP->Fill(double(iphi));
      else        nrechits_vs_iphi_HOM->Fill(double(iphi));
       }
       if ( sub == 4){
      emean_vs_ieta_HF[depth-1]->Fill(double(ieta), en);
      occupancy_map_HF[depth-1]->Fill(double(ieta),double(iphi));
      if (ieta>0) nrechits_vs_iphi_HFP[depth-1]->Fill(double(iphi));
      else        nrechits_vs_iphi_HFM[depth-1]->Fill(double(iphi));
       }
     }
 
 
     
     //32-bit status word  
     uint32_t statadd;
     unsigned int isw67 = 0;
 
     //Statusword correlation
     unsigned int sw27 = 27;
     unsigned int sw13 = 13;
 
     uint32_t statadd27 = 0x1<<sw27;
     uint32_t statadd13 = 0x1<<sw13;
 
     float status27 = 0;
     float status13 = 0;
 
     if(stwd & statadd27) status27 = 1;
     if(stwd & statadd13) status13 = 1;
 
     if        (sub == 1){
       RecHit_StatusWordCorr_HB->Fill(status13, status27);
     } else if (sub == 2){
       RecHit_StatusWordCorr_HE->Fill(status13, status27);
     }
 
 
     for (unsigned int isw = 0; isw < 32; isw++){
       statadd = 0x1<<(isw);
       if (stwd & statadd){
     if      (sub == 1) RecHit_StatusWord_HB->Fill(isw);
     else if (sub == 2) RecHit_StatusWord_HE->Fill(isw);
     else if (sub == 3) RecHit_StatusWord_HO->Fill(isw);
     else if (sub == 4){
       RecHit_StatusWord_HF->Fill(isw);
       if (isw == 6) isw67 += 1;
       if (isw == 7) isw67 += 2;
     }
       }
     }
 
     for (unsigned int isw =0; isw < 32; isw++){
       statadd = 0x1<<(isw);
       if( auxstwd & statadd ){
         if      (sub == 1) RecHit_Aux_StatusWord_HB->Fill(isw);
         else if (sub == 2) RecHit_Aux_StatusWord_HE->Fill(isw);
         else if (sub == 3) RecHit_Aux_StatusWord_HO->Fill(isw);
         else if (sub == 4) RecHit_Aux_StatusWord_HF->Fill(isw);
       }
 
     }
 
   } 
 
     for(int depth = 0; depth <= maxDepthHB_; depth++) Nhb[depth]->Fill(double(nhb_v[depth]));
     for(int depth = 0; depth <= maxDepthHE_; depth++) Nhe[depth]->Fill(double(nhe_v[depth]));
     for(int depth = 0; depth <= maxDepthHO_; depth++) Nho[depth]->Fill(double(nho_v[depth]));
     for(int depth = 0; depth <= maxDepthHF_; depth++) Nhf[depth]->Fill(double(nhf_v[depth]));
 
   //===========================================================================
   // SUBSYSTEMS,  
   //===========================================================================
   
   if ((subdet_ != 6) && (subdet_ != 0)) {
 
     double clusEta = 999.;
     double clusPhi = 999.; 
     double clusEn  = 0.;
     
     double HcalCone    = 0.;
 
     int ietaMax   =  9999;
     //     double enMax1 = -9999.;
     //     double enMax2 = -9999.;
     //     double enMax3 = -9999.;
     //     double enMax4 = -9999.;
     //     double enMax  = -9999.;
     //     double etaMax =  9999.;
 
     //   CYCLE over cells ====================================================
 
     for (unsigned int i = 0; i < cen.size(); i++) {
       int sub    = csub[i];
       double eta = ceta[i]; 
       double phi = cphi[i]; 
       double ieta = cieta[i]; 
       double iphi = ciphi[i]; 
       double en  = cen[i]; 
       double enM0  = cenM0[i]; 
       double enM3  = cenM3[i]; 
       double chi2  = cchi2[i];
       double chi2_log10=9.99; // initial value - stay with this value if chi2<0.
       if (chi2>0.) chi2_log10=log10(chi2);
       double t   = ctime[i];
       double depth = cdepth[i];
 
       bool isHEP17 = (sub == 2) && (iphi >= 63) && (iphi <= 66) && (ieta > 0) && (hep17_);
 
 //       int   ieta = cieta[i];
 
       double rhot = dR(etaHot, phiHot, eta, phi); 
       if(rhot < partR && en > 1.) { 
     clusEta = (clusEta * clusEn + eta * en)/(clusEn + en);
         clusPhi = phi12(clusPhi, clusEn, phi, en); 
         clusEn += en;
       }
 
       nrechits++;       
       eHcal += en;
 
       if(en > 1. ) nrechitsThresh++;
       
       //The energy and overall timing histos are drawn while
       //the ones split by depth are not
       if(sub == 1 && (subdet_ == 1 || subdet_ == 5)) {  
     meTimeHB->Fill(t);
     meRecHitsEnergyHB->Fill(en);
     meRecHitsEnergyHBM0->Fill(enM0);
     meRecHitsEnergyHBM3->Fill(enM3);
 
         meRecHitsEnergyM2vM0HB->Fill(enM0,en);
         meRecHitsEnergyM3vM0HB->Fill(enM0,enM3);
         meRecHitsEnergyM3vM2HB->Fill(en,enM3);
 
         meRecHitsM2Chi2HB->Fill(chi2_log10);
         meLog10Chi2profileHB->Fill(en,chi2_log10);
     
     meTE_Low_HB->Fill( en, t);
     meTE_HB->Fill( en, t);
     meTE_High_HB->Fill( en, t);
     meTEprofileHB_Low->Fill(en, t);
     meTEprofileHB->Fill(en, t);
     meTEprofileHB_High->Fill(en, t);
       }     
       if(sub == 2 && (subdet_ == 2 || subdet_ == 5)) {  
     meTimeHE->Fill(t);
         if(!isHEP17){
        meRecHitsEnergyHE->Fill(en);
        meRecHitsEnergyHEM0->Fill(enM0);
        meRecHitsEnergyHEM3->Fill(enM3);
         }else{
        meRecHitsEnergyHEP17[0]->Fill(en);
        meRecHitsEnergyHEP17M0[0]->Fill(enM0);
        meRecHitsEnergyHEP17M3[0]->Fill(enM3);
        meRecHitsEnergyHEP17[depth]->Fill(en);
        meRecHitsEnergyHEP17M0[depth]->Fill(enM0);
        meRecHitsEnergyHEP17M3[depth]->Fill(enM3);
         }
 
         meRecHitsEnergyM2vM0HE->Fill(enM0,en);
         meRecHitsEnergyM3vM0HE->Fill(enM0,enM3);
         meRecHitsEnergyM3vM2HE->Fill(en,enM3);
 
         meRecHitsM2Chi2HE->Fill(chi2_log10);
         meLog10Chi2profileHE->Fill(en,chi2_log10);  
 
     meTE_Low_HE->Fill( en, t);
     meTE_HE->Fill( en, t);
     meTEprofileHE_Low->Fill(en, t);
     meTEprofileHE->Fill(en, t);
       }
       if(sub == 4 && (subdet_ == 4 || subdet_ == 5)) {  
     meTimeHF->Fill(t);
     meRecHitsEnergyHF->Fill(en);      
 
     meTE_Low_HF->Fill(en, t);
     meTE_HF->Fill(en, t);
     meTEprofileHF_Low->Fill(en, t);
     meTEprofileHF->Fill(en, t);
 
       }
       if(sub == 3 && (subdet_ == 3 || subdet_ == 5)) {  
     meTimeHO->Fill(t);
     meRecHitsEnergyHO->Fill(en);
 
     meTE_HO->Fill( en, t);
     meTE_High_HO->Fill( en, t);
     meTEprofileHO->Fill(en, t);
     meTEprofileHO_High->Fill(en, t);
       }
     }
 
     if(imc != 0) {
       //Cone by depth are not drawn, the others are used for pion scan
       meEnConeEtaProfile       ->Fill(double(ietaMax),  HcalCone);   // 
       meEnConeEtaProfile_E     ->Fill(double(ietaMax), eEcalCone);   
       meEnConeEtaProfile_EH    ->Fill(double(ietaMax),  HcalCone+eEcalCone); 
     }
 
     // Single particle samples ONLY !  ======================================
     // Fill up some histos for "integrated" subsustems. 
     // These are not drawn
   }
 
   nevtot++;
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 120 of file HcalRecHitsAnalyzer.cc.

 { 
 
     Char_t histo[200];
 
     ibooker.setCurrentFolder(topFolderName_);
 
     // General counters (drawn)
 
     //Produce both a total per subdetector, and number of rechits per subdetector depth
     // The bins are 1 unit wide, and the range is determined by the number of channels per subdetector
 
     for(int depth = 0; depth <= maxDepthHB_; depth++){
       if(depth == 0){ sprintf  (histo, "N_HB" );}
       else{           sprintf  (histo, "N_HB_depth%d",depth );}
       int NBins = (int) (nChannels_[1] * 1.1);
       Nhb.push_back( ibooker.book1D(histo, histo, NBins, 0., (float)NBins) );
     } 
     for(int depth = 0; depth <= maxDepthHE_; depth++){
       if(depth == 0){ sprintf  (histo, "N_HE" );}
       else{           sprintf  (histo, "N_HE_depth%d",depth );}
       int NBins = (int) (nChannels_[2] * 1.1);
       Nhe.push_back( ibooker.book1D(histo, histo, NBins,0., (float)NBins) );
     } 
     for(int depth = 0; depth <= maxDepthHO_; depth++){
       if(depth == 0){ sprintf  (histo, "N_HO" );}
       else{           sprintf  (histo, "N_HO_depth%d",depth );}
       int NBins = (int) (nChannels_[3] * 1.1);
       Nho.push_back( ibooker.book1D(histo, histo, NBins,0., (float)NBins) );
     } 
     for(int depth = 0; depth <= maxDepthHF_; depth++){
       if(depth == 0){ sprintf  (histo, "N_HF" );}
       else{           sprintf  (histo, "N_HF_depth%d",depth );}
       int NBins = (int) (nChannels_[4] * 1.1);
       Nhf.push_back( ibooker.book1D(histo, histo, NBins,0., (float)NBins) );
     } 
 
     // ZS
     if(subdet_ == 6) {
 
     }
 
     // ALL others, except ZS
     else {  
       for(int depth = 1; depth <= maxDepthAll_; depth++){
         sprintf  (histo, "emap_depth%d",depth );
         emap.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
       } 
 
       //The mean energy histos are drawn, but not the RMS or emean seq
       
       for (int depth = 1; depth <= maxDepthHB_; depth++) {
     sprintf  (histo, "emean_vs_ieta_HB%d",depth );
     emean_vs_ieta_HB.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
 
     sprintf  (histo, "emean_vs_ieta_M0_HB%d",depth );
     emean_vs_ieta_HBM0.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
 
     sprintf  (histo, "emean_vs_ieta_M3_HB%d",depth );
     emean_vs_ieta_HBM3.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
       }
       for (int depth = 1; depth <= maxDepthHE_; depth++) {
     sprintf  (histo, "emean_vs_ieta_HE%d",depth );
     emean_vs_ieta_HE.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
 
     sprintf  (histo, "emean_vs_ieta_M0_HE%d",depth );
     emean_vs_ieta_HEM0.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
 
     sprintf  (histo, "emean_vs_ieta_M3_HE%d",depth );
     emean_vs_ieta_HEM3.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
       }
 
       if(hep17_){
          for (int depth = 1; depth <= maxDepthHE_; depth++) {
             sprintf  (histo, "emean_vs_ieta_HEP17_depth%d",depth );
             emean_vs_ieta_HEP17.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
 
         sprintf  (histo, "emean_vs_ieta_M0_HEP17_depth%d",depth );
         emean_vs_ieta_HEP17M0.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
 
         sprintf  (histo, "emean_vs_ieta_M3_HEP17_depth%d",depth );
         emean_vs_ieta_HEP17M3.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
          }
       }
 
       for (int depth = 1; depth <= maxDepthHF_; depth++) {
     sprintf  (histo, "emean_vs_ieta_HF%d",depth );
     emean_vs_ieta_HF.push_back( ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " ") );
       }
       sprintf  (histo, "emean_vs_ieta_HO" );
       emean_vs_ieta_HO = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -10., 2000., " " );
 
       //The only occupancy histos drawn are occupancy vs. ieta
       //but the maps are needed because this is where the latter are filled from
 
       for (int depth = 1; depth <= maxDepthHB_; depth++) {
          sprintf  (histo, "occupancy_map_HB%d",depth );
          occupancy_map_HB.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
       }
 
       for (int depth = 1; depth <= maxDepthHE_; depth++) {
          sprintf  (histo, "occupancy_map_HE%d",depth );
          occupancy_map_HE.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
       }
 
       sprintf  (histo, "occupancy_map_HO" );
       occupancy_map_HO = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);      
 
       for (int depth = 1; depth <= maxDepthHF_; depth++) {
          sprintf  (histo, "occupancy_map_HF%d",depth );
          occupancy_map_HF.push_back( ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_) );
       }
       
       // nrechits vs iphi
       for (int depth = 1; depth <= maxDepthHB_; depth++) {
     
      sprintf  (histo, "occupancy_vs_ieta_HB%d",depth );
      occupancy_vs_ieta_HB.push_back( ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_) );
          sprintf  (histo, "nrechits_vs_iphi_HBP_d%d",depth );
          nrechits_vs_iphi_HBP.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
          sprintf  (histo, "nrechits_vs_iphi_HBM_d%d",depth );
          nrechits_vs_iphi_HBM.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
       }
 
       for (int depth = 1; depth <= maxDepthHE_; depth++) {
          sprintf  (histo, "occupancy_vs_ieta_HE%d",depth );
          occupancy_vs_ieta_HE.push_back( ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_) );
          sprintf  (histo, "nrechits_vs_iphi_HEP_d%d",depth );
          nrechits_vs_iphi_HEP.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
          sprintf  (histo, "nrechits_vs_iphi_HEM_d%d",depth );
          nrechits_vs_iphi_HEM.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
       }
 
       sprintf  (histo, "occupancy_vs_ieta_HO" );
       occupancy_vs_ieta_HO = ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_);
       sprintf  (histo, "nrechits_vs_iphi_HOP" );
       nrechits_vs_iphi_HOP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);
       sprintf  (histo, "nrechits_vs_iphi_HOM" );
       nrechits_vs_iphi_HOM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);
 
       for (int depth = 1; depth <= maxDepthHF_; depth++) {
          sprintf  (histo, "occupancy_vs_ieta_HF%d",depth );
          occupancy_vs_ieta_HF.push_back( ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_) );
          sprintf  (histo, "nrechits_vs_iphi_HFP_d%d",depth );
          nrechits_vs_iphi_HFP.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
          sprintf  (histo, "nrechits_vs_iphi_HFM_d%d",depth );
          nrechits_vs_iphi_HFM.push_back( ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_) );
       }
 
       //All status word histos except HF67 are drawn
       sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HB" ) ;
       RecHit_StatusWord_HB = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
       
       sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HE" ) ;
       RecHit_StatusWord_HE = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
 
       sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HF" ) ;
       RecHit_StatusWord_HF = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
 
       sprintf (histo, "HcalRecHitTask_RecHit_StatusWord_HO" ) ;
       RecHit_StatusWord_HO = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
 
       //Aux status word histos
       sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HB" ) ;
       RecHit_Aux_StatusWord_HB = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
       
       sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HE" ) ;
       RecHit_Aux_StatusWord_HE = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
 
       sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HF" ) ;
       RecHit_Aux_StatusWord_HF = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
 
       sprintf (histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HO" ) ;
       RecHit_Aux_StatusWord_HO = ibooker.book1D(histo, histo, 32 , -0.5, 31.5); 
 
     }  // end-of (subdet_ =! 6)
 
       //Status word correlations
       sprintf (histo, "HcalRecHitTask_RecHit_StatusWordCorr_HB");
       RecHit_StatusWordCorr_HB = ibooker.book2D(histo, histo, 2, -0.5, 1.5, 2, -0.5, 1.5);
 
       sprintf (histo, "HcalRecHitTask_RecHit_StatusWordCorr_HE");
       RecHit_StatusWordCorr_HE = ibooker.book2D(histo, histo, 2, -0.5, 1.5, 2, -0.5, 1.5);
 
 
     //======================= Now various cases one by one ===================
 
     //Histograms drawn for single pion scan
     if(subdet_ != 0 && imc != 0) { // just not for noise  
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths");
       meEnConeEtaProfile = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -100., 2000., " ");  
       
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
       meEnConeEtaProfile_E = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -100., 2000., " ");  
       
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
       meEnConeEtaProfile_EH = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, -100., 2000., " ");  
     }
 
     // ************** HB **********************************
     if (subdet_ == 1 || subdet_ == 5 ){
 
       //Only severity level, energy of rechits and overall HB timing histos are drawn  
 
       sprintf(histo, "HcalRecHitTask_severityLevel_HB");
       sevLvl_HB = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HB" ) ;
       meRecHitsEnergyHB = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HB" ) ;
       meRecHitsEnergyHBM0 = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HB" ) ;
       meRecHitsEnergyHBM3 = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M2vM0_HB" ) ;
       meRecHitsEnergyM2vM0HB = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM0_HB" ) ;
       meRecHitsEnergyM3vM0HB = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM2_HB" ) ;
       meRecHitsEnergyM3vM2HB = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
       
       sprintf (histo, "HcalRecHitTask_M2Log10Chi2_of_rechits_HB" ) ;
       meRecHitsM2Chi2HB = ibooker.book1D(histo, histo, 120 , -2. , 10.); 
       
       sprintf (histo, "HcalRecHitTask_timing_HB" ) ;
       meTimeHB = ibooker.book1D(histo, histo, 70, -48., 92.); 
 
       //High, medium and low histograms to reduce RAM usage
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HB" ) ;
       meTE_Low_HB = ibooker.book2D(histo, histo, 50, -5., 45.,  70, -48., 92.);
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_HB" ) ;
       meTE_HB = ibooker.book2D(histo, histo, 150, -5., 295.,  70, -48., 92.);
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_High_HB" ) ;
       meTE_High_HB = ibooker.book2D(histo, histo, 150, -5., 2995.,  70, -48., 92.);
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HB" ) ;
       meTEprofileHB_Low = ibooker.bookProfile(histo, histo, 50, -5., 45., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HB" ) ;
       meTEprofileHB = ibooker.bookProfile(histo, histo, 150, -5., 295., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_Log10Chi2_vs_energy_profile_HB" ) ;
       meLog10Chi2profileHB = ibooker.bookProfile(histo, histo, 150, -5., 295., -2., 9.9, " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HB" ) ;
       meTEprofileHB_High = ibooker.bookProfile(histo, histo, 150, -5., 2995., -48., 92., " "); 
 
     }
     
     // ********************** HE ************************************
     if ( subdet_ == 2 || subdet_ == 5 ){
 
 
       //Only severity level, energy of rechits and overall HB timing histos are drawn  
       sprintf(histo, "HcalRecHitTask_severityLevel_HE");
       sevLvl_HE = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HE" ) ;
       meRecHitsEnergyHE = ibooker.book1D(histo, histo, 2010, -10., 2000.);
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HE" ) ;
       meRecHitsEnergyHEM0 = ibooker.book1D(histo, histo, 2010, -10., 2000.);
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HE" ) ;
       meRecHitsEnergyHEM3 = ibooker.book1D(histo, histo, 2010, -10., 2000.);
       
       if(hep17_){
          sprintf (histo, "HcalRecHitTask_energy_of_rechits_HEP17" ) ;
          meRecHitsEnergyHEP17.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
       
          sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HEP17" ) ;
          meRecHitsEnergyHEP17M0.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
       
          sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HEP17" ) ;
          meRecHitsEnergyHEP17M3.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
          for (int depth = 1; depth <= maxDepthHE_; depth++) {
             sprintf (histo, "HcalRecHitTask_energy_of_rechits_HEP17_depth%d",depth ) ;
             meRecHitsEnergyHEP17.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
       
             sprintf (histo, "HcalRecHitTask_energy_of_rechits_M0_HEP17_depth%d",depth ) ;
             meRecHitsEnergyHEP17M0.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
       
             sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3_HEP17_depth%d",depth ) ;
             meRecHitsEnergyHEP17M3.push_back(ibooker.book1D(histo, histo, 2010 , -10. , 2000.)); 
          }
       }
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M2vM0_HE" ) ;
       meRecHitsEnergyM2vM0HE = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM0_HE" ) ;
       meRecHitsEnergyM3vM0HE = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
       
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_M3vM2_HE" ) ;
       meRecHitsEnergyM3vM2HE = ibooker.book2D(histo, histo, 42 , -10. , 200., 42, -10., 200.); 
       
       sprintf (histo, "HcalRecHitTask_M2Log10Chi2_of_rechits_HE" ) ;
       meRecHitsM2Chi2HE = ibooker.book1D(histo, histo, 120 , -2. , 10.); 
       
       sprintf (histo, "HcalRecHitTask_timing_HE" ) ;
       meTimeHE = ibooker.book1D(histo, histo, 70, -48., 92.); 
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HE" ) ;
       meTE_Low_HE = ibooker.book2D(histo, histo, 80, -5., 75.,  70, -48., 92.);
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_HE" ) ;
       meTE_HE = ibooker.book2D(histo, histo, 200, -5., 395.,  70, -48., 92.);
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HE" ) ;
       meTEprofileHE_Low = ibooker.bookProfile(histo, histo, 80, -5., 75., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HE" ) ;
       meTEprofileHE = ibooker.bookProfile(histo, histo, 200, -5., 395., -48., 92., " "); 
       
       sprintf (histo, "HcalRecHitTask_Log10Chi2_vs_energy_profile_HE" ) ;
       meLog10Chi2profileHE = ibooker.bookProfile(histo, histo, 200, -5., 395., -2., 9.9, " ");       
 
     }
 
     // ************** HO ****************************************
     if ( subdet_ == 3 || subdet_ == 5  ){
       
       //Only severity level, energy of rechits and overall HB timing histos are drawn  
 
       sprintf(histo, "HcalRecHitTask_severityLevel_HO");
       sevLvl_HO = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HO" ) ;
       meRecHitsEnergyHO = ibooker.book1D(histo, histo, 2010 , -10. , 2000.);
       
       sprintf (histo, "HcalRecHitTask_timing_HO" ) ;
       meTimeHO = ibooker.book1D(histo, histo, 70, -48., 92.); 
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_HO" ) ;
       meTE_HO= ibooker.book2D(histo, histo, 60, -5., 55., 70, -48., 92.);
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_High_HO" ) ;
       meTE_High_HO= ibooker.book2D(histo, histo, 100, -5., 995., 70, -48., 92.);
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HO" ) ;
       meTEprofileHO = ibooker.bookProfile(histo, histo, 60, -5., 55., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HO" ) ;
       meTEprofileHO_High = ibooker.bookProfile(histo, histo, 100, -5., 995., -48., 92., " "); 
       
     }   
   
     // ********************** HF ************************************
     if ( subdet_ == 4 || subdet_ == 5 ){
 
       //Only severity level, energy of rechits and overall HB timing histos are drawn  
       
       sprintf(histo, "HcalRecHitTask_severityLevel_HF");
       sevLvl_HF = ibooker.book1D(histo, histo, 25, -0.5, 24.5); 
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HF" ) ;
       meRecHitsEnergyHF = ibooker.book1D(histo, histo, 2010 , -10. , 2000.); 
 
       sprintf (histo, "HcalRecHitTask_timing_HF" ) ;
       meTimeHF = ibooker.book1D(histo, histo, 70, -48., 92.); 
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_Low_HF" ) ;
       meTE_Low_HF = ibooker.book2D(histo, histo, 100, -5., 195., 70, -48., 92.);
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_HF" ) ;
       meTE_HF = ibooker.book2D(histo, histo, 200, -5., 995., 70, -48., 92.);
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HF" ) ;
       meTEprofileHF_Low = ibooker.bookProfile(histo, histo, 100, -5., 195., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HF" ) ;
       meTEprofileHF = ibooker.bookProfile(histo, histo, 200, -5., 995., -48., 92., " "); 
 
     }
 
 }

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

private

Definition at line 1111 of file HcalRecHitsAnalyzer.cc.

References PI, and tmp.

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

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

override

Definition at line 59 of file HcalRecHitsAnalyzer.cc.

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

                                                                                  {
   
     edm::ESHandle<HcalDDDRecConstants> pHRNDC;
     es.get<HcalRecNumberingRecord>().get( pHRNDC );
     hcons = &(*pHRNDC);
     maxDepthHB_ = hcons->getMaxDepth(0);
     maxDepthHE_ = hcons->getMaxDepth(1);
     maxDepthHF_ = hcons->getMaxDepth(2);
     maxDepthHO_ = hcons->getMaxDepth(3);
 
     edm::ESHandle<CaloGeometry> geometry;
 
     es.get<CaloGeometryRecord > ().get(geometry);
 
     const CaloGeometry* geo = geometry.product();
     const HcalGeometry* gHB = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel));
     const HcalGeometry* gHE = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap));
     const HcalGeometry* gHO = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalOuter));
     const HcalGeometry* gHF = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,HcalForward));
 
     nChannels_[1] = gHB->getHxSize(1); 
     nChannels_[2] = gHE->getHxSize(2); 
     nChannels_[3] = gHO->getHxSize(3); 
     nChannels_[4] = gHF->getHxSize(4); 
 
     nChannels_[0] = nChannels_[1] + nChannels_[2] + nChannels_[3] + nChannels_[4];
 
     //std::cout << "Channels HB:" << nChannels_[1] << " HE:" << nChannels_[2] << " HO:" << nChannels_[3] << " HF:" << nChannels_[4] << std::endl;
 
     //We hardcode the HF depths because in the dual readout configuration, rechits are not defined for depths 3&4
     maxDepthHF_ = (maxDepthHF_ > 2 ? 2 : maxDepthHF_); //We retain the dynamic possibility that HF might have 0 or 1 depths
 
     maxDepthAll_ = ( maxDepthHB_ + maxDepthHO_ > maxDepthHE_ ? maxDepthHB_ + maxDepthHO_ : maxDepthHE_ );
     maxDepthAll_ = ( maxDepthAll_ > maxDepthHF_ ? maxDepthAll_ : maxDepthHF_ );
 
     //Get Phi segmentation from geometry, use the max phi number so that all iphi values are included.
 
     int NphiMax = hcons->getNPhi(0);
 
     NphiMax = (hcons->getNPhi(1) > NphiMax ? hcons->getNPhi(1) : NphiMax);
     NphiMax = (hcons->getNPhi(2) > NphiMax ? hcons->getNPhi(2) : NphiMax);
     NphiMax = (hcons->getNPhi(3) > NphiMax ? hcons->getNPhi(3) : NphiMax);
 
     //Center the iphi bins on the integers
     iphi_min_ = 0.5;
     iphi_max_ = NphiMax + 0.5;
     iphi_bins_ = (int) (iphi_max_ - iphi_min_);
 
     //Retain classic behavior, all plots have same ieta range.
 
     int iEtaMax = (hcons->getEtaRange(0).second > hcons->getEtaRange(1).second ? hcons->getEtaRange(0).second : hcons->getEtaRange(1).second);
     iEtaMax = (iEtaMax > hcons->getEtaRange(2).second ? iEtaMax : hcons->getEtaRange(2).second);
     iEtaMax = (iEtaMax > hcons->getEtaRange(3).second ? iEtaMax : hcons->getEtaRange(3).second);
 
     //Give an empty bin around the subdet ieta range to make it clear that all ieta rings have been included
     ieta_min_ = -iEtaMax - 1.5;
     ieta_max_ = iEtaMax + 1.5;
     ieta_bins_ = (int) (ieta_max_ - ieta_min_);
 
   }

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

private

Definition at line 1085 of file HcalRecHitsAnalyzer.cc.

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

Referenced by analyze().

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

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

privatevirtual

Definition at line 906 of file HcalRecHitsAnalyzer.cc.

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

Referenced by analyze().

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

int HcalRecHitsAnalyzer::hcalSevLvl ( const CaloRecHit * hit )

private

Definition at line 1126 of file HcalRecHitsAnalyzer.cc.

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

Referenced by fillRecHitsTmp().

                                                         {
 
    HcalDetId id = hit->detid();
    if (theHcalTopology->withSpecialRBXHBHE() && id.subdet() == HcalEndcap) {
      id = theHcalTopology->idFront(id);
    }
 
    const uint32_t recHitFlag = hit->flags();
    const uint32_t dbStatusFlag = theHcalChStatus->getValues(id)->getValue();
 
    int severityLevel = theHcalSevLvlComputer->getSeverityLevel(id, recHitFlag, dbStatusFlag);
 
    return severityLevel;
 
 } 

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

private

Definition at line 1095 of file HcalRecHitsAnalyzer.cc.

References PI, and tmp.

Referenced by analyze().

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

Member Data Documentation

std::vector<uint32_t> HcalRecHitsAnalyzer::cauxstwd

private

Definition at line 294 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cchi2

private

Definition at line 288 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::cdepth

private

Definition at line 284 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cen

private

Definition at line 285 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cenM0

private

Definition at line 286 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cenM3

private

Definition at line 287 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::ceta

private

Definition at line 289 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::cieta

private

Definition at line 282 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::ciphi

private

Definition at line 283 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cphi

private

Definition at line 290 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<uint32_t> HcalRecHitsAnalyzer::cstwd

private

Definition at line 293 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::csub

private

Definition at line 281 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::ctime

private

Definition at line 291 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cz

private

Definition at line 292 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

std::string HcalRecHitsAnalyzer::ecalselector_

private

Definition at line 77 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

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

private

Definition at line 139 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 141 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 142 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 143 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 144 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 145 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 146 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 147 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 148 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 149 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 150 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::emean_vs_ieta_HO

private

Definition at line 151 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsAnalyzer::etype_

private

Definition at line 107 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

std::string HcalRecHitsAnalyzer::eventype_

private

Definition at line 78 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

bool HcalRecHitsAnalyzer::famos_

private

Definition at line 82 of file HcalRecHitsAnalyzer.h.

edm::ESHandle<CaloGeometry> HcalRecHitsAnalyzer::geometry

private

Definition at line 261 of file HcalRecHitsAnalyzer.h.

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

std::vector<int> HcalRecHitsAnalyzer::hcalHBSevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHESevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHFSevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHOSevLvlVec

private

Definition at line 119 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::string HcalRecHitsAnalyzer::hcalselector_

private

Definition at line 76 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

const HcalDDDRecConstants* HcalRecHitsAnalyzer::hcons

private

Definition at line 84 of file HcalRecHitsAnalyzer.h.

Referenced by dqmBeginRun().

bool HcalRecHitsAnalyzer::hep17_

private

Definition at line 80 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::ieta_bins_

private

Definition at line 93 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::ieta_max_

private

Definition at line 94 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::ieta_min_

private

Definition at line 94 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::imc

private

Definition at line 109 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::iphi_bins_

private

Definition at line 90 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::iphi_max_

private

Definition at line 91 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::iphi_min_

private

Definition at line 91 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::iz

private

Definition at line 108 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

MonitorElement* HcalRecHitsAnalyzer::map_ecal

private

Definition at line 137 of file HcalRecHitsAnalyzer.h.

int HcalRecHitsAnalyzer::maxDepthAll_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHB_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHE_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHF_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHO_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

std::string HcalRecHitsAnalyzer::mc_

private

Definition at line 81 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meDeltaEta

private

Definition at line 174 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meDeltaPhi

private

Definition at line 173 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHB

private

Definition at line 248 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHE

private

Definition at line 249 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHF

private

Definition at line 251 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHO

private

Definition at line 250 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHB

private

Definition at line 245 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHE

private

Definition at line 246 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile

private

Definition at line 169 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_E

private

Definition at line 170 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_EH

private

Definition at line 171 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHB

private

Definition at line 254 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHE

private

Definition at line 255 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHB

private

Definition at line 211 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHE

private

Definition at line 218 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHB

private

Definition at line 258 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHE

private

Definition at line 259 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHB

private

Definition at line 183 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM0

private

Definition at line 184 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM3

private

Definition at line 185 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHE

private

Definition at line 191 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM0

private

Definition at line 192 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM3

private

Definition at line 193 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 194 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 195 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 196 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHF

private

Definition at line 204 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHO

private

Definition at line 202 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HB

private

Definition at line 186 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HE

private

Definition at line 197 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HB

private

Definition at line 187 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HE

private

Definition at line 198 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HB

private

Definition at line 188 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HE

private

Definition at line 199 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HB

private

Definition at line 189 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HE

private

Definition at line 200 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHB

private

Definition at line 237 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHE

private

Definition at line 238 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHF

private

Definition at line 240 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFL

private

Definition at line 241 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFS

private

Definition at line 242 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHO

private

Definition at line 239 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHB

private

Definition at line 231 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHE

private

Definition at line 232 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHF

private

Definition at line 234 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHO

private

Definition at line 233 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meTE_HB

private

Definition at line 207 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HE

private

Definition at line 215 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HF

private

Definition at line 226 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HB

private

Definition at line 208 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HO

private

Definition at line 221 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HO

private

Definition at line 220 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HB

private

Definition at line 206 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HE

private

Definition at line 214 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HF

private

Definition at line 225 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB

private

Definition at line 210 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_High

private

Definition at line 212 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_Low

private

Definition at line 209 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE

private

Definition at line 217 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE_Low

private

Definition at line 216 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF

private

Definition at line 228 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF_Low

private

Definition at line 227 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO

private

Definition at line 222 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO_High

private

Definition at line 223 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHB

private

Definition at line 177 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHE

private

Definition at line 178 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHF

private

Definition at line 180 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHO

private

Definition at line 179 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsAnalyzer::nChannels_[5]

private

Definition at line 88 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::nevtot

private

Definition at line 297 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

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

private

Definition at line 127 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 128 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 130 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 129 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 163 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 163 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * HcalRecHitsAnalyzer::nrechits_vs_iphi_HOM

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::nrechits_vs_iphi_HOP

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 153 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 154 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 155 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::occupancy_map_HO

private

Definition at line 156 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 158 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

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

private

Definition at line 159 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

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

private

Definition at line 160 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::occupancy_vs_ieta_HO

private

Definition at line 161 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

std::string HcalRecHitsAnalyzer::outputFile_

private

Definition at line 75 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HB

private

Definition at line 275 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HE

private

Definition at line 276 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HF

private

Definition at line 277 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HO

private

Definition at line 278 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HB

private

Definition at line 264 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HE

private

Definition at line 265 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF

private

Definition at line 266 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF67

private

Definition at line 267 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HO

private

Definition at line 268 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HB

private

Definition at line 271 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HE

private

Definition at line 272 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HB

private

Definition at line 121 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HE

private

Definition at line 122 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HF

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HO

private

Definition at line 124 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string HcalRecHitsAnalyzer::sign_

private

Definition at line 79 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

int HcalRecHitsAnalyzer::subdet_

private

Definition at line 104 of file HcalRecHitsAnalyzer.h.

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

const HcalChannelQuality* HcalRecHitsAnalyzer::theHcalChStatus

private

Definition at line 114 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

const HcalSeverityLevelComputer* HcalRecHitsAnalyzer::theHcalSevLvlComputer

private

Definition at line 116 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

const HcalTopology* HcalRecHitsAnalyzer::theHcalTopology

private

Definition at line 112 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

edm::EDGetTokenT<EBRecHitCollection> HcalRecHitsAnalyzer::tok_EB_

private

Definition at line 100 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<EERecHitCollection> HcalRecHitsAnalyzer::tok_EE_

private

Definition at line 101 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalRecHitsAnalyzer::tok_hbhe_

private

Definition at line 97 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HFRecHitCollection> HcalRecHitsAnalyzer::tok_hf_

private

Definition at line 99 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HORecHitCollection> HcalRecHitsAnalyzer::tok_ho_

private

Definition at line 98 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

std::string HcalRecHitsAnalyzer::topFolderName_

private

Definition at line 73 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation