HcalRecHitsAnalyzer Class Reference

#include <HcalRecHitsAnalyzer.h>

Inheritance diagram for HcalRecHitsAnalyzer:

Public Member Functions
void	analyze (edm::Event const &ev, edm::EventSetup const &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &) override
	HcalRecHitsAnalyzer (edm::ParameterSet const &conf)
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryEventToken_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryRunToken_
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< int >	csevlev
std::vector< uint32_t >	cstwd
std::vector< int >	csub
std::vector< double >	ctime
std::vector< double >	cz
std::string	ecalselector_
std::vector< MonitorElement * >	emap
MonitorElement *	emap_HO
std::vector< MonitorElement * >	emean_vs_ieta_HB
std::vector< MonitorElement * >	emean_vs_ieta_HBM0
std::vector< MonitorElement * >	emean_vs_ieta_HBM3
std::vector< MonitorElement * >	emean_vs_ieta_HE
std::vector< MonitorElement * >	emean_vs_ieta_HEM0
std::vector< MonitorElement * >	emean_vs_ieta_HEM3
std::vector< MonitorElement * >	emean_vs_ieta_HEP17
std::vector< MonitorElement * >	emean_vs_ieta_HEP17M0
std::vector< MonitorElement * >	emean_vs_ieta_HEP17M3
std::vector< MonitorElement * >	emean_vs_ieta_HF
MonitorElement *	emean_vs_ieta_HO
int	etype_
std::string	eventype_
bool	famos_
CaloGeometry const *	geometry = nullptr
edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd >	hcalChannelQualityToken_
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	hcalDDDRecConstantsToken_
std::vector< int >	hcalHBSevLvlVec
std::vector< int >	hcalHESevLvlVec
std::vector< int >	hcalHFSevLvlVec
std::vector< int >	hcalHOSevLvlVec
std::string	hcalselector_
edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	hcalSeverityLevelComputerToken_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	hcalTopologyToken_
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 *	meRecHitsCleanedEnergyHB
MonitorElement *	meRecHitsCleanedEnergyHE
MonitorElement *	meRecHitsCleanedEnergyHF
MonitorElement *	meRecHitsCleanedEnergyHO
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 = nullptr
const HcalSeverityLevelComputer *	theHcalSevLvlComputer = nullptr
const HcalTopology *	theHcalTopology = nullptr
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 DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 56 of file HcalRecHitsAnalyzer.h.

Constructor & Destructor Documentation

HcalRecHitsAnalyzer()

HcalRecHitsAnalyzer::HcalRecHitsAnalyzer

(

edm::ParameterSet const &

conf

)

Definition at line 7 of file HcalRecHitsAnalyzer.cc.

    : topFolderName_(conf.getParameter<std::string>("TopFolderName")),
      hcalDDDRecConstantsToken_{esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()},
      caloGeometryRunToken_{esConsumes<CaloGeometry, CaloGeometryRecord, edm::Transition::BeginRun>()},
      caloGeometryEventToken_{esConsumes<CaloGeometry, CaloGeometryRecord>()},
      hcalTopologyToken_{esConsumes<HcalTopology, HcalRecNumberingRecord>()},
      hcalChannelQualityToken_{esConsumes<HcalChannelQuality, HcalChannelQualityRcd>(edm::ESInputTag("", "withTopo"))},
      hcalSeverityLevelComputerToken_{esConsumes<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd>()} {
  // 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

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 547 of file HcalRecHitsAnalyzer.cc.

References caloGeometryEventToken_, cauxstwd, cchi2, cdepth, cen, cenM0, cenM3, ceta, nano_mu_local_reco_cff::chi2, cieta, ciphi, cphi, csevlev, cstwd, csub, ctime, hcalRecHitTable_cff::depth, dR(), emap, emap_HO, emean_vs_ieta_HB, emean_vs_ieta_HBM0, emean_vs_ieta_HBM3, emean_vs_ieta_HE, emean_vs_ieta_HEM0, emean_vs_ieta_HEM3, emean_vs_ieta_HEP17, emean_vs_ieta_HEP17M0, emean_vs_ieta_HEP17M3, emean_vs_ieta_HF, emean_vs_ieta_HO, PVValHelper::eta, makeMEIFBenchmarkPlots::ev, dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), fillRecHitsTmp(), edm::EventSetup::getData(), hcalChannelQualityToken_, hcalHBSevLvlVec, hcalHESevLvlVec, hcalHFSevLvlVec, hcalHOSevLvlVec, hcalSeverityLevelComputerToken_, hcalTopologyToken_, hep17_, mps_fire::i, hcalRecHitTable_cff::ieta, hcaltpdigi_cfi::ieta1, hcaltpdigi_cfi::ieta2, imc, hcalRecHitTable_cff::iphi, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meLog10Chi2profileHB, meLog10Chi2profileHE, meRecHitsCleanedEnergyHB, meRecHitsCleanedEnergyHE, meRecHitsCleanedEnergyHF, meRecHitsCleanedEnergyHO, meRecHitsEnergyHB, meRecHitsEnergyHBM0, meRecHitsEnergyHBM3, meRecHitsEnergyHE, meRecHitsEnergyHEM0, meRecHitsEnergyHEM3, meRecHitsEnergyHEP17, meRecHitsEnergyHEP17M0, meRecHitsEnergyHEP17M3, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitsEnergyM2vM0HB, meRecHitsEnergyM2vM0HE, meRecHitsEnergyM3vM0HB, meRecHitsEnergyM3vM0HE, meRecHitsEnergyM3vM2HB, meRecHitsEnergyM3vM2HE, meRecHitsM2Chi2HB, meRecHitsM2Chi2HE, meTE_HB, meTE_HE, meTE_HF, meTE_High_HB, meTE_High_HO, meTE_HO, meTE_Low_HB, meTE_Low_HE, meTE_Low_HF, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, meTimeHB, meTimeHE, meTimeHF, meTimeHO, nevtot, Nhb, Nhe, Nhf, Nho, nrechits_vs_iphi_HBM, nrechits_vs_iphi_HBP, nrechits_vs_iphi_HEM, nrechits_vs_iphi_HEP, nrechits_vs_iphi_HFM, nrechits_vs_iphi_HFP, nrechits_vs_iphi_HOM, nrechits_vs_iphi_HOP, occupancy_map_HB, occupancy_map_HE, occupancy_map_HF, occupancy_map_HO, phi, phi12(), RecHit_Aux_StatusWord_HB, RecHit_Aux_StatusWord_HE, RecHit_Aux_StatusWord_HF, RecHit_Aux_StatusWord_HO, RecHit_StatusWord_HB, RecHit_StatusWord_HE, RecHit_StatusWord_HF, RecHit_StatusWord_HO, RecHit_StatusWordCorr_HB, RecHit_StatusWordCorr_HE, sevLvl_HB, sevLvl_HE, sevLvl_HF, sevLvl_HO, subdet_, submitPVValidationJobs::t, theHcalChStatus, theHcalSevLvlComputer, and theHcalTopology.

                                                                                 {
  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 ECAL
  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:  iSetup.get<IdealGeometryRecord>().get (geometry);
  geometry = &iSetup.getData(caloGeometryEventToken_);

  // HCAL Topology **************************************************
  theHcalTopology = &iSetup.getData(hcalTopologyToken_);

  // HCAL channel status map ****************************************
  theHcalChStatus = &iSetup.getData(hcalChannelQualityToken_);

  // Assignment of severity levels **********************************
  theHcalSevLvlComputer = &iSetup.getData(hcalSeverityLevelComputerToken_);

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

  // Counting, including ZS items
  // Filling HCAL maps  ----------------------------------------------------
  //   double maxE = -99999.;

  // element 0: any depth. element 1,2,..: depth 1,2
  std::vector<int> nhb_v(maxDepthHB_ + 1, 0);
  std::vector<int> nhe_v(maxDepthHE_ + 1, 0);
  std::vector<int> nho_v(maxDepthHO_ + 1, 0);
  std::vector<int> nhf_v(maxDepthHF_ + 1, 0);

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

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

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

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

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

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

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

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

    // 32-bit status word
    uint32_t statadd;

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

    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];
      int sevlev = csevlev[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;
      }

      // 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);
        if (sevlev <= 9)
          meRecHitsCleanedEnergyHB->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);
          if (sevlev <= 9)
            meRecHitsCleanedEnergyHE->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);
        if (sevlev <= 9)
          meRecHitsCleanedEnergyHF->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);
        if (sevlev <= 9)
          meRecHitsCleanedEnergyHO->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++;
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 133 of file HcalRecHitsAnalyzer.cc.

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

{
  Char_t histo[200];

  ibooker.setCurrentFolder(topFolderName_);

  // General counters (drawn)

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

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

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

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

    // The mean energy histos are drawn, but not the RMS or emean seq

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    sprintf(histo, "HcalRecHitTask_RecHit_StatusWord_HE");
    RecHit_StatusWord_HE = ibooker.book1DD(histo, histo, 32, -0.5, 31.5);

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

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

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

    sprintf(histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HE");
    RecHit_Aux_StatusWord_HE = ibooker.book1DD(histo, histo, 32, -0.5, 31.5);

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

    sprintf(histo, "HcalRecHitTask_RecHit_Aux_StatusWord_HO");
    RecHit_Aux_StatusWord_HO = ibooker.book1DD(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.book1DD(histo, histo, 25, -0.5, 24.5);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HB");
    meRecHitsEnergyHB = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_cleaned_energy_of_rechits_HB");
    meRecHitsCleanedEnergyHB = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_M0_HB");
    meRecHitsEnergyHBM0 = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_M3_HB");
    meRecHitsEnergyHBM3 = ibooker.book1DD(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.book1DD(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.book1DD(histo, histo, 25, -0.5, 24.5);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HE");
    meRecHitsEnergyHE = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_cleaned_energy_of_rechits_HE");
    meRecHitsCleanedEnergyHE = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_M0_HE");
    meRecHitsEnergyHEM0 = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_M3_HE");
    meRecHitsEnergyHEM3 = ibooker.book1DD(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.book1DD(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.book1DD(histo, histo, 25, -0.5, 24.5);

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HO");
    meRecHitsEnergyHO = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_cleaned_energy_of_rechits_HO");
    meRecHitsCleanedEnergyHO = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_timing_HO");
    meTimeHO = ibooker.book1DD(histo, histo, 80, -80., 80.);

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_HO");
    meTE_HO = ibooker.book2D(histo, histo, 60, -5., 55., 80, -80., 80.);

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_High_HO");
    meTE_High_HO = ibooker.book2D(histo, histo, 100, -5., 995., 80, -80., 80.);

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_profile_HO");
    meTEprofileHO = ibooker.bookProfile(histo, histo, 60, -5., 55., -80., 80., " ");

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_profile_High_HO");
    meTEprofileHO_High = ibooker.bookProfile(histo, histo, 100, -5., 995., -80., 80., " ");
  }

  // ********************** 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.book1DD(histo, histo, 25, -0.5, 24.5);

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

    sprintf(histo, "HcalRecHitTask_cleaned_energy_of_rechits_HF");
    meRecHitsCleanedEnergyHF = ibooker.book1DD(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_timing_HF");
    meTimeHF = ibooker.book1DD(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., " ");
  }
}

dPhiWsign()

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

private

Definition at line 1144 of file HcalRecHitsAnalyzer.cc.

References testProducerWithPsetDescEmpty_cfi::a2, PI, and createJobs::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;
}

dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 73 of file HcalRecHitsAnalyzer.cc.

References caloGeometryRunToken_, edm::EventSetup::getData(), HcalDDDRecConstants::getEtaRange(), HcalGeometry::getHxSize(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNPhi(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, hcalDDDRecConstantsToken_, HcalEndcap, HcalForward, HcalOuter, ieta_bins_, ieta_max_, ieta_min_, hcalTTPDigis_cfi::iEtaMax, createfilelist::int, iphi_bins_, iphi_max_, iphi_min_, SiStripPI::max, maxDepthAll_, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, and nChannels_.

                                                                                    {
  HcalDDDRecConstants const &hcons = iSetup.getData(hcalDDDRecConstantsToken_);
  maxDepthHB_ = hcons.getMaxDepth(0);
  maxDepthHE_ = hcons.getMaxDepth(1);
  maxDepthHF_ = std::max(hcons.getMaxDepth(2), 1);
  maxDepthHO_ = hcons.getMaxDepth(3);

  CaloGeometry const &geo = iSetup.getData(caloGeometryRunToken_);

  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] = std::max(int(gHE->getHxSize(2)), 1);
  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_);
}

dR()

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

private

Definition at line 1111 of file HcalRecHitsAnalyzer.cc.

References HLT_2024v14_cff::eta1, HLT_2024v14_cff::eta2, mathSSE::sqrt(), and createJobs::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;
}

fillRecHitsTmp()

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

privatevirtual

Definition at line 940 of file HcalRecHitsAnalyzer.cc.

References edm::SortedCollection< T, SORT >::begin(), cauxstwd, cchi2, cdepth, cen, cenM0, cenM3, ceta, nano_mu_local_reco_cff::chi2, cieta, ciphi, cphi, csevlev, cstwd, csub, ctime, cz, hcalRecHitTable_cff::depth, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, makeMEIFBenchmarkPlots::ev, HcalBarrel, HcalEndcap, HcalForward, hcalHBSevLvlVec, hcalHESevLvlVec, hcalHFSevLvlVec, hcalHOSevLvlVec, HcalOuter, hcalSevLvl(), iz, dqmiolumiharvest::j, phi, interestingDetIdCollectionProducer_cfi::severityLevel, subdet_, submitPVValidationJobs::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();
  csevlev.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);
          csevlev.push_back(severityLevel);
        }
      }
    }
  }

  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);
          csevlev.push_back(severityLevel);
        }
      }
    }
  }

  // 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);
          csevlev.push_back(severityLevel);
        }
      }
    }
  }
}

hcalSevLvl()

int HcalRecHitsAnalyzer::hcalSevLvl ( const CaloRecHit *  hit )

private

Definition at line 1161 of file HcalRecHitsAnalyzer.cc.

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

Referenced by fillRecHitsTmp().

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

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

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

  return severityLevel;
}

phi12()

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

private

Definition at line 1125 of file HcalRecHitsAnalyzer.cc.

References testProducerWithPsetDescEmpty_cfi::a2, PI, and createJobs::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

caloGeometryEventToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalRecHitsAnalyzer::caloGeometryEventToken_

private

Definition at line 102 of file HcalRecHitsAnalyzer.h.

Referenced by analyze().

caloGeometryRunToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalRecHitsAnalyzer::caloGeometryRunToken_

private

Definition at line 101 of file HcalRecHitsAnalyzer.h.

Referenced by dqmBeginRun().

cauxstwd

std::vector<uint32_t> HcalRecHitsAnalyzer::cauxstwd

private

Definition at line 300 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cchi2

std::vector<double> HcalRecHitsAnalyzer::cchi2

private

Definition at line 294 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cdepth

std::vector<int> HcalRecHitsAnalyzer::cdepth

private

Definition at line 290 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cen

std::vector<double> HcalRecHitsAnalyzer::cen

private

Definition at line 291 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cenM0

std::vector<double> HcalRecHitsAnalyzer::cenM0

private

Definition at line 292 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cenM3

std::vector<double> HcalRecHitsAnalyzer::cenM3

private

Definition at line 293 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

ceta

std::vector<double> HcalRecHitsAnalyzer::ceta

private

Definition at line 295 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cieta

std::vector<int> HcalRecHitsAnalyzer::cieta

private

Definition at line 288 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

ciphi

std::vector<int> HcalRecHitsAnalyzer::ciphi

private

Definition at line 289 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cphi

std::vector<double> HcalRecHitsAnalyzer::cphi

private

Definition at line 296 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

csevlev

std::vector<int> HcalRecHitsAnalyzer::csevlev

private

Definition at line 301 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cstwd

std::vector<uint32_t> HcalRecHitsAnalyzer::cstwd

private

Definition at line 299 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

csub

std::vector<int> HcalRecHitsAnalyzer::csub

private

Definition at line 287 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

ctime

std::vector<double> HcalRecHitsAnalyzer::ctime

private

Definition at line 297 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

cz

std::vector<double> HcalRecHitsAnalyzer::cz

private

Definition at line 298 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

ecalselector_

std::string HcalRecHitsAnalyzer::ecalselector_

private

Definition at line 76 of file HcalRecHitsAnalyzer.h.

emap

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

private

Definition at line 143 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emap_HO

MonitorElement* HcalRecHitsAnalyzer::emap_HO

private

Definition at line 144 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HB

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

private

Definition at line 146 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HBM0

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

private

Definition at line 147 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HBM3

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

private

Definition at line 148 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HE

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

private

Definition at line 149 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HEM0

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

private

Definition at line 150 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HEM3

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

private

Definition at line 151 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HEP17

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

private

Definition at line 152 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HEP17M0

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

private

Definition at line 153 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HEP17M3

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

private

Definition at line 154 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HF

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

private

Definition at line 155 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

emean_vs_ieta_HO

MonitorElement* HcalRecHitsAnalyzer::emean_vs_ieta_HO

private

Definition at line 156 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

etype_

int HcalRecHitsAnalyzer::etype_

private

Definition at line 111 of file HcalRecHitsAnalyzer.h.

eventype_

std::string HcalRecHitsAnalyzer::eventype_

private

Definition at line 77 of file HcalRecHitsAnalyzer.h.

famos_

bool HcalRecHitsAnalyzer::famos_

private

Definition at line 81 of file HcalRecHitsAnalyzer.h.

geometry

CaloGeometry const* HcalRecHitsAnalyzer::geometry = nullptr

private

Definition at line 267 of file HcalRecHitsAnalyzer.h.

hcalChannelQualityToken_

edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> HcalRecHitsAnalyzer::hcalChannelQualityToken_

private

Definition at line 104 of file HcalRecHitsAnalyzer.h.

Referenced by analyze().

hcalDDDRecConstantsToken_

edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalRecHitsAnalyzer::hcalDDDRecConstantsToken_

private

Definition at line 100 of file HcalRecHitsAnalyzer.h.

Referenced by dqmBeginRun().

hcalHBSevLvlVec

std::vector<int> HcalRecHitsAnalyzer::hcalHBSevLvlVec

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

hcalHESevLvlVec

std::vector<int> HcalRecHitsAnalyzer::hcalHESevLvlVec

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

hcalHFSevLvlVec

std::vector<int> HcalRecHitsAnalyzer::hcalHFSevLvlVec

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

hcalHOSevLvlVec

std::vector<int> HcalRecHitsAnalyzer::hcalHOSevLvlVec

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

hcalselector_

std::string HcalRecHitsAnalyzer::hcalselector_

private

Definition at line 75 of file HcalRecHitsAnalyzer.h.

hcalSeverityLevelComputerToken_

edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> HcalRecHitsAnalyzer::hcalSeverityLevelComputerToken_

private

Definition at line 105 of file HcalRecHitsAnalyzer.h.

Referenced by analyze().

hcalTopologyToken_

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> HcalRecHitsAnalyzer::hcalTopologyToken_

private

Definition at line 103 of file HcalRecHitsAnalyzer.h.

Referenced by analyze().

hep17_

bool HcalRecHitsAnalyzer::hep17_

private

Definition at line 79 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

ieta_bins_

int HcalRecHitsAnalyzer::ieta_bins_

private

Definition at line 90 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

ieta_max_

float HcalRecHitsAnalyzer::ieta_max_

private

Definition at line 91 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

ieta_min_

float HcalRecHitsAnalyzer::ieta_min_

private

Definition at line 91 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

imc

int HcalRecHitsAnalyzer::imc

private

Definition at line 113 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

iphi_bins_

int HcalRecHitsAnalyzer::iphi_bins_

private

Definition at line 87 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

iphi_max_

float HcalRecHitsAnalyzer::iphi_max_

private

Definition at line 88 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

iphi_min_

float HcalRecHitsAnalyzer::iphi_min_

private

Definition at line 88 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

iz

int HcalRecHitsAnalyzer::iz

private

Definition at line 112 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

map_ecal

MonitorElement* HcalRecHitsAnalyzer::map_ecal

private

Definition at line 141 of file HcalRecHitsAnalyzer.h.

maxDepthAll_

int HcalRecHitsAnalyzer::maxDepthAll_

private

Definition at line 83 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

maxDepthHB_

int HcalRecHitsAnalyzer::maxDepthHB_

private

Definition at line 83 of file HcalRecHitsAnalyzer.h.

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

maxDepthHE_

int HcalRecHitsAnalyzer::maxDepthHE_

private

Definition at line 83 of file HcalRecHitsAnalyzer.h.

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

maxDepthHF_

int HcalRecHitsAnalyzer::maxDepthHF_

private

Definition at line 83 of file HcalRecHitsAnalyzer.h.

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

maxDepthHO_

int HcalRecHitsAnalyzer::maxDepthHO_

private

Definition at line 83 of file HcalRecHitsAnalyzer.h.

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

mc_

std::string HcalRecHitsAnalyzer::mc_

private

Definition at line 80 of file HcalRecHitsAnalyzer.h.

meDeltaEta

MonitorElement* HcalRecHitsAnalyzer::meDeltaEta

private

Definition at line 179 of file HcalRecHitsAnalyzer.h.

meDeltaPhi

MonitorElement* HcalRecHitsAnalyzer::meDeltaPhi

private

Definition at line 178 of file HcalRecHitsAnalyzer.h.

meEcalHcalEnergyConeHB

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHB

private

Definition at line 254 of file HcalRecHitsAnalyzer.h.

meEcalHcalEnergyConeHE

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHE

private

Definition at line 255 of file HcalRecHitsAnalyzer.h.

meEcalHcalEnergyConeHF

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHF

private

Definition at line 257 of file HcalRecHitsAnalyzer.h.

meEcalHcalEnergyConeHO

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHO

private

Definition at line 256 of file HcalRecHitsAnalyzer.h.

meEcalHcalEnergyHB

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHB

private

Definition at line 251 of file HcalRecHitsAnalyzer.h.

meEcalHcalEnergyHE

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHE

private

Definition at line 252 of file HcalRecHitsAnalyzer.h.

meEnConeEtaProfile

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile

private

Definition at line 174 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meEnConeEtaProfile_E

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_E

private

Definition at line 175 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meEnConeEtaProfile_EH

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_EH

private

Definition at line 176 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meEnergyHcalVsEcalHB

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHB

private

Definition at line 260 of file HcalRecHitsAnalyzer.h.

meEnergyHcalVsEcalHE

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHE

private

Definition at line 261 of file HcalRecHitsAnalyzer.h.

meLog10Chi2profileHB

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHB

private

Definition at line 220 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meLog10Chi2profileHE

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHE

private

Definition at line 227 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meNumEcalRecHitsConeHB

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHB

private

Definition at line 264 of file HcalRecHitsAnalyzer.h.

meNumEcalRecHitsConeHE

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHE

private

Definition at line 265 of file HcalRecHitsAnalyzer.h.

meRecHitsCleanedEnergyHB

MonitorElement* HcalRecHitsAnalyzer::meRecHitsCleanedEnergyHB

private

Definition at line 189 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsCleanedEnergyHE

MonitorElement* HcalRecHitsAnalyzer::meRecHitsCleanedEnergyHE

private

Definition at line 198 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsCleanedEnergyHF

MonitorElement* HcalRecHitsAnalyzer::meRecHitsCleanedEnergyHF

private

Definition at line 213 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsCleanedEnergyHO

MonitorElement* HcalRecHitsAnalyzer::meRecHitsCleanedEnergyHO

private

Definition at line 210 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHB

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHB

private

Definition at line 188 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHBM0

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM0

private

Definition at line 190 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHBM3

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM3

private

Definition at line 191 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHE

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHE

private

Definition at line 197 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHEM0

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM0

private

Definition at line 199 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHEM3

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM3

private

Definition at line 200 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHEP17

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

private

Definition at line 201 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHEP17M0

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

private

Definition at line 202 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHEP17M3

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

private

Definition at line 203 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHF

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHF

private

Definition at line 212 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHO

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHO

private

Definition at line 209 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyM2vM0HB

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HB

private

Definition at line 192 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyM2vM0HE

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HE

private

Definition at line 204 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyM3vM0HB

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HB

private

Definition at line 193 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyM3vM0HE

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HE

private

Definition at line 205 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyM3vM2HB

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HB

private

Definition at line 194 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyM3vM2HE

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HE

private

Definition at line 206 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsM2Chi2HB

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HB

private

Definition at line 195 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meRecHitsM2Chi2HE

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HE

private

Definition at line 207 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meSumRecHitsEnergyConeHB

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHB

private

Definition at line 244 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyConeHE

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHE

private

Definition at line 245 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyConeHF

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHF

private

Definition at line 247 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyConeHFL

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFL

private

Definition at line 248 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyConeHFS

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFS

private

Definition at line 249 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyConeHO

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHO

private

Definition at line 246 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyHB

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHB

private

Definition at line 239 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyHE

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHE

private

Definition at line 240 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyHF

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHF

private

Definition at line 242 of file HcalRecHitsAnalyzer.h.

meSumRecHitsEnergyHO

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHO

private

Definition at line 241 of file HcalRecHitsAnalyzer.h.

meTE_HB

MonitorElement* HcalRecHitsAnalyzer::meTE_HB

private

Definition at line 216 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_HE

MonitorElement* HcalRecHitsAnalyzer::meTE_HE

private

Definition at line 224 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_HF

MonitorElement* HcalRecHitsAnalyzer::meTE_HF

private

Definition at line 235 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_High_HB

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HB

private

Definition at line 217 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_High_HO

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HO

private

Definition at line 230 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_HO

MonitorElement* HcalRecHitsAnalyzer::meTE_HO

private

Definition at line 229 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_Low_HB

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HB

private

Definition at line 215 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_Low_HE

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HE

private

Definition at line 223 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTE_Low_HF

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HF

private

Definition at line 234 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHB

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB

private

Definition at line 219 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHB_High

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_High

private

Definition at line 221 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHB_Low

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_Low

private

Definition at line 218 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHE

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE

private

Definition at line 226 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHE_Low

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE_Low

private

Definition at line 225 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHF

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF

private

Definition at line 237 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHF_Low

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF_Low

private

Definition at line 236 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHO

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO

private

Definition at line 231 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHO_High

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO_High

private

Definition at line 232 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTimeHB

MonitorElement* HcalRecHitsAnalyzer::meTimeHB

private

Definition at line 182 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTimeHE

MonitorElement* HcalRecHitsAnalyzer::meTimeHE

private

Definition at line 183 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTimeHF

MonitorElement* HcalRecHitsAnalyzer::meTimeHF

private

Definition at line 185 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

meTimeHO

MonitorElement* HcalRecHitsAnalyzer::meTimeHO

private

Definition at line 184 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nChannels_

int HcalRecHitsAnalyzer::nChannels_[5]

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

nevtot

int HcalRecHitsAnalyzer::nevtot

private

Definition at line 304 of file HcalRecHitsAnalyzer.h.

Referenced by analyze().

Nhb

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

private

Definition at line 131 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

Nhe

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

private

Definition at line 132 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

Nhf

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

private

Definition at line 134 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

Nho

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

private

Definition at line 133 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HBM

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

private

Definition at line 168 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HBP

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

private

Definition at line 168 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HEM

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

private

Definition at line 169 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HEP

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

private

Definition at line 169 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HFM

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

private

Definition at line 170 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HFP

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

private

Definition at line 170 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HOM

MonitorElement * HcalRecHitsAnalyzer::nrechits_vs_iphi_HOM

private

Definition at line 171 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

nrechits_vs_iphi_HOP

MonitorElement* HcalRecHitsAnalyzer::nrechits_vs_iphi_HOP

private

Definition at line 171 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

occupancy_map_HB

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

private

Definition at line 158 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

occupancy_map_HE

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

private

Definition at line 159 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

occupancy_map_HF

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

private

Definition at line 160 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

occupancy_map_HO

MonitorElement* HcalRecHitsAnalyzer::occupancy_map_HO

private

Definition at line 161 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

occupancy_vs_ieta_HB

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

private

Definition at line 163 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

occupancy_vs_ieta_HE

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

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

occupancy_vs_ieta_HF

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

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

occupancy_vs_ieta_HO

MonitorElement* HcalRecHitsAnalyzer::occupancy_vs_ieta_HO

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

outputFile_

std::string HcalRecHitsAnalyzer::outputFile_

private

Definition at line 74 of file HcalRecHitsAnalyzer.h.

RecHit_Aux_StatusWord_HB

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HB

private

Definition at line 281 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_Aux_StatusWord_HE

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HE

private

Definition at line 282 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_Aux_StatusWord_HF

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HF

private

Definition at line 283 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_Aux_StatusWord_HO

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HO

private

Definition at line 284 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_StatusWord_HB

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HB

private

Definition at line 270 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_StatusWord_HE

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HE

private

Definition at line 271 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_StatusWord_HF

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF

private

Definition at line 272 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_StatusWord_HF67

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF67

private

Definition at line 273 of file HcalRecHitsAnalyzer.h.

RecHit_StatusWord_HO

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HO

private

Definition at line 274 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_StatusWordCorr_HB

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HB

private

Definition at line 277 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

RecHit_StatusWordCorr_HE

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HE

private

Definition at line 278 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

sevLvl_HB

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HB

private

Definition at line 125 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

sevLvl_HE

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HE

private

Definition at line 126 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

sevLvl_HF

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HF

private

Definition at line 127 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

sevLvl_HO

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HO

private

Definition at line 128 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

sign_

std::string HcalRecHitsAnalyzer::sign_

private

Definition at line 78 of file HcalRecHitsAnalyzer.h.

subdet_

int HcalRecHitsAnalyzer::subdet_

private

Definition at line 108 of file HcalRecHitsAnalyzer.h.

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

theHcalChStatus

const HcalChannelQuality* HcalRecHitsAnalyzer::theHcalChStatus = nullptr

private

Definition at line 118 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

theHcalSevLvlComputer

const HcalSeverityLevelComputer* HcalRecHitsAnalyzer::theHcalSevLvlComputer = nullptr

private

Definition at line 120 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

theHcalTopology

const HcalTopology* HcalRecHitsAnalyzer::theHcalTopology = nullptr

private

Definition at line 116 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

tok_EB_

edm::EDGetTokenT<EBRecHitCollection> HcalRecHitsAnalyzer::tok_EB_

private

Definition at line 97 of file HcalRecHitsAnalyzer.h.

tok_EE_

edm::EDGetTokenT<EERecHitCollection> HcalRecHitsAnalyzer::tok_EE_

private

Definition at line 98 of file HcalRecHitsAnalyzer.h.

tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalRecHitsAnalyzer::tok_hbhe_

private

Definition at line 94 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

tok_hf_

edm::EDGetTokenT<HFRecHitCollection> HcalRecHitsAnalyzer::tok_hf_

private

Definition at line 96 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

tok_ho_

edm::EDGetTokenT<HORecHitCollection> HcalRecHitsAnalyzer::tok_ho_

private

Definition at line 95 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

topFolderName_

std::string HcalRecHitsAnalyzer::topFolderName_

private

Definition at line 72 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().