HcalRecHitsAnalyzer Class Reference

#include <HcalRecHitsAnalyzer.h>

Inheritance diagram for HcalRecHitsAnalyzer:

Public Member Functions
void	analyze (edm::Event const &ev, edm::EventSetup const &c) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &c) override
	HcalRecHitsAnalyzer (edm::ParameterSet const &conf)
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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

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

Detailed Description

Definition at line 57 of file HcalRecHitsAnalyzer.h.

Constructor & Destructor Documentation

HcalRecHitsAnalyzer::HcalRecHitsAnalyzer

(

edm::ParameterSet const &

conf

)

Definition at line 6 of file HcalRecHitsAnalyzer.cc.

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

    : topFolderName_(conf.getParameter<std::string>("TopFolderName")) {
  // DQM ROOT output
  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");

  if (!outputFile_.empty()) {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
  } else {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
  }

  nevtot = 0;

  hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");
  ecalselector_ = conf.getUntrackedParameter<std::string>("ecalselector", "yes");
  eventype_ = conf.getUntrackedParameter<std::string>("eventype", "single");
  sign_ = conf.getUntrackedParameter<std::string>("sign", "*");
  // useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);

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

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

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

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

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

  imc = 0;
}

Member Function Documentation

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

override

Definition at line 533 of file HcalRecHitsAnalyzer.cc.

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

                                                                            {
  using namespace edm;

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

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

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

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

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

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

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

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

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

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

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

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

  // ECAL
  if (ecalselector_ == "yes" && (subdet_ == 1 || subdet_ == 2 || subdet_ == 5)) {
    Handle<EBRecHitCollection> rhitEB;
    if (ev.getByToken(tok_EB_, rhitEB)) {
      for (const auto &recHit : *(rhitEB.product())) {
        double en = recHit.energy();
        eEcal += en;
        eEcalB += en;
      }
    }

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    float status27 = 0;
    float status13 = 0;

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

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

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

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

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

  //===========================================================================
  // SUBSYSTEMS,
  //===========================================================================

  if ((subdet_ != 6) && (subdet_ != 0)) {
    double clusEta = 999.;
    double clusPhi = 999.;
    double clusEn = 0.;

    double HcalCone = 0.;

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

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

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

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

      //       int   ieta = cieta[i];

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

      nrechits++;
      eHcal += en;

      if (en > 1.)
        nrechitsThresh++;

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

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

        meRecHitsM2Chi2HB->Fill(chi2_log10);
        meLog10Chi2profileHB->Fill(en, chi2_log10);

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

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

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

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

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

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

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

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

  nevtot++;
}

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

override

Definition at line 131 of file HcalRecHitsAnalyzer.cc.

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

{
  Char_t histo[200];

  ibooker.setCurrentFolder(topFolderName_);

  // General counters (drawn)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    // All status word histos except HF67 are drawn
    sprintf(histo, "HcalRecHitTask_RecHit_StatusWord_HB");
    RecHit_StatusWord_HB = ibooker.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_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_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_timing_HO");
    meTimeHO = ibooker.book1DD(histo, histo, 70, -48., 92.);

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

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

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

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

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

    sprintf(histo, "HcalRecHitTask_severityLevel_HF");
    sevLvl_HF = ibooker.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_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., " ");
  }
}

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

private

Definition at line 1166 of file HcalRecHitsAnalyzer.cc.

References PI, and tmp.

                                                              {
  // clockwise      phi2 w.r.t phi1 means "+" phi distance
  // anti-clockwise phi2 w.r.t phi1 means "-" phi distance

  double PI = 3.1415926535898;
  double a1 = phi1;
  double a2 = phi2;
  double tmp = a2 - a1;
  if (a1 * a2 < 0.) {
    if (a1 > 0.5 * PI)
      tmp += 2. * PI;
    if (a2 > 0.5 * PI)
      tmp -= 2. * PI;
  }
  return tmp;
}

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

override

Definition at line 66 of file HcalRecHitsAnalyzer.cc.

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

                                                                                {
  edm::ESHandle<HcalDDDRecConstants> pHRNDC;
  es.get<HcalRecNumberingRecord>().get(pHRNDC);
  hcons = &(*pHRNDC);
  maxDepthHB_ = hcons->getMaxDepth(0);
  maxDepthHE_ = hcons->getMaxDepth(1);
  maxDepthHF_ = hcons->getMaxDepth(2);
  maxDepthHO_ = hcons->getMaxDepth(3);

  edm::ESHandle<CaloGeometry> geometry;

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

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

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

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

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

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

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

  // Get Phi segmentation from geometry, use the max phi number so that all iphi
  // values are included.

  int NphiMax = hcons->getNPhi(0);

  NphiMax = (hcons->getNPhi(1) > NphiMax ? hcons->getNPhi(1) : NphiMax);
  NphiMax = (hcons->getNPhi(2) > NphiMax ? hcons->getNPhi(2) : NphiMax);
  NphiMax = (hcons->getNPhi(3) > NphiMax ? hcons->getNPhi(3) : NphiMax);

  // Center the iphi bins on the integers
  iphi_min_ = 0.5;
  iphi_max_ = NphiMax + 0.5;
  iphi_bins_ = (int)(iphi_max_ - iphi_min_);

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

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

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

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

private

Definition at line 1133 of file HcalRecHitsAnalyzer.cc.

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

Referenced by analyze().

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

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

privatevirtual

Definition at line 966 of file HcalRecHitsAnalyzer.cc.

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

Referenced by analyze().

                                                                        {
  using namespace edm;

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

  if (subdet_ == 1 || subdet_ == 2 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    // HBHE
    edm::Handle<HBHERecHitCollection> hbhecoll;
    if (ev.getByToken(tok_hbhe_, hbhecoll)) {
      for (HBHERecHitCollection::const_iterator j = hbhecoll->begin(); j != hbhecoll->end(); j++) {
        HcalDetId cell(j->id());
        const HcalGeometry *cellGeometry = dynamic_cast<const HcalGeometry *>(geometry->getSubdetectorGeometry(cell));
        double eta = cellGeometry->getPosition(cell).eta();
        double phi = cellGeometry->getPosition(cell).phi();
        double zc = cellGeometry->getPosition(cell).z();
        int sub = cell.subdet();
        int depth = cell.depth();
        int inteta = cell.ieta();
        int intphi = cell.iphi();
        double en = j->energy();
        double enM0 = j->eraw();
        double enM3 = j->eaux();
        double chi2 = j->chi2();
        double t = j->time();
        int stwd = j->flags();
        int auxstwd = j->aux();

        int severityLevel = hcalSevLvl((CaloRecHit *)&*j);
        if (cell.subdet() == HcalBarrel) {
          hcalHBSevLvlVec.push_back(severityLevel);
        } else if (cell.subdet() == HcalEndcap) {
          hcalHESevLvlVec.push_back(severityLevel);
        }

        if ((iz > 0 && eta > 0.) || (iz < 0 && eta < 0.) || iz == 0) {
          csub.push_back(sub);
          cen.push_back(en);
          cenM0.push_back(enM0);
          cenM3.push_back(enM3);
          cchi2.push_back(chi2);
          ceta.push_back(eta);
          cphi.push_back(phi);
          ctime.push_back(t);
          cieta.push_back(inteta);
          ciphi.push_back(intphi);
          cdepth.push_back(depth);
          cz.push_back(zc);
          cstwd.push_back(stwd);
          cauxstwd.push_back(auxstwd);
        }
      }
    }
  }

  if (subdet_ == 4 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    // HF
    edm::Handle<HFRecHitCollection> hfcoll;
    if (ev.getByToken(tok_hf_, hfcoll)) {
      for (HFRecHitCollection::const_iterator j = hfcoll->begin(); j != hfcoll->end(); j++) {
        HcalDetId cell(j->id());
        auto cellGeometry = (geometry->getSubdetectorGeometry(cell))->getGeometry(cell);
        double eta = cellGeometry->getPosition().eta();
        double phi = cellGeometry->getPosition().phi();
        double zc = cellGeometry->getPosition().z();
        int sub = cell.subdet();
        int depth = cell.depth();
        int inteta = cell.ieta();
        int intphi = cell.iphi();
        double en = j->energy();
        double enM0 = 0.;
        double enM3 = 0.;
        double chi2 = 0.;
        double t = j->time();
        int stwd = j->flags();
        int auxstwd = j->aux();

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

        if ((iz > 0 && eta > 0.) || (iz < 0 && eta < 0.) || iz == 0) {
          csub.push_back(sub);
          cen.push_back(en);
          cenM0.push_back(enM0);
          cenM3.push_back(enM3);
          cchi2.push_back(chi2);
          ceta.push_back(eta);
          cphi.push_back(phi);
          ctime.push_back(t);
          cieta.push_back(inteta);
          ciphi.push_back(intphi);
          cdepth.push_back(depth);
          cz.push_back(zc);
          cstwd.push_back(stwd);
          cauxstwd.push_back(auxstwd);
        }
      }
    }
  }

  // HO
  if (subdet_ == 3 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    edm::Handle<HORecHitCollection> hocoll;
    if (ev.getByToken(tok_ho_, hocoll)) {
      for (HORecHitCollection::const_iterator j = hocoll->begin(); j != hocoll->end(); j++) {
        HcalDetId cell(j->id());
        auto cellGeometry = (geometry->getSubdetectorGeometry(cell))->getGeometry(cell);
        double eta = cellGeometry->getPosition().eta();
        double phi = cellGeometry->getPosition().phi();
        double zc = cellGeometry->getPosition().z();
        int sub = cell.subdet();
        int depth = cell.depth();
        int inteta = cell.ieta();
        int intphi = cell.iphi();
        double t = j->time();
        double en = j->energy();
        double enM0 = 0.;
        double enM3 = 0.;
        double chi2 = 0.;
        int stwd = j->flags();
        int auxstwd = j->aux();

        int severityLevel = hcalSevLvl((CaloRecHit *)&*j);
        if (cell.subdet() == HcalOuter) {
          hcalHOSevLvlVec.push_back(severityLevel);
        }

        if ((iz > 0 && eta > 0.) || (iz < 0 && eta < 0.) || iz == 0) {
          csub.push_back(sub);
          cen.push_back(en);
          cenM0.push_back(enM0);
          cenM3.push_back(enM3);
          cchi2.push_back(chi2);
          ceta.push_back(eta);
          cphi.push_back(phi);
          ctime.push_back(t);
          cieta.push_back(inteta);
          ciphi.push_back(intphi);
          cdepth.push_back(depth);
          cz.push_back(zc);
          cstwd.push_back(stwd);
          cauxstwd.push_back(auxstwd);
        }
      }
    }
  }
}

int HcalRecHitsAnalyzer::hcalSevLvl ( const CaloRecHit *  hit )

private

Definition at line 1183 of file HcalRecHitsAnalyzer.cc.

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

Referenced by fillRecHitsTmp().

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

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

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

  return severityLevel;
}

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

private

Definition at line 1147 of file HcalRecHitsAnalyzer.cc.

References PI, and tmp.

Referenced by analyze().

                                                                                  {
  // weighted mean value of phi1 and phi2

  double tmp;
  double PI = 3.1415926535898;
  double a1 = phi1;
  double a2 = phi2;

  if (a1 > 0.5 * PI && a2 < 0.)
    a2 += 2 * PI;
  if (a2 > 0.5 * PI && a1 < 0.)
    a1 += 2 * PI;
  tmp = (a1 * en1 + a2 * en2) / (en1 + en2);
  if (tmp > PI)
    tmp -= 2. * PI;

  return tmp;
}

Member Data Documentation

std::vector<uint32_t> HcalRecHitsAnalyzer::cauxstwd

private

Definition at line 291 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cchi2

private

Definition at line 285 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::cdepth

private

Definition at line 281 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cen

private

Definition at line 282 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cenM0

private

Definition at line 283 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cenM3

private

Definition at line 284 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::ceta

private

Definition at line 286 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::cieta

private

Definition at line 279 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::ciphi

private

Definition at line 280 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cphi

private

Definition at line 287 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<uint32_t> HcalRecHitsAnalyzer::cstwd

private

Definition at line 290 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::csub

private

Definition at line 278 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::ctime

private

Definition at line 288 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<double> HcalRecHitsAnalyzer::cz

private

Definition at line 289 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp().

std::string HcalRecHitsAnalyzer::ecalselector_

private

Definition at line 77 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

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

private

Definition at line 138 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::emap_HO

private

Definition at line 139 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 141 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 142 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 143 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 144 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 145 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 146 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 147 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 148 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 149 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 150 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::emean_vs_ieta_HO

private

Definition at line 151 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsAnalyzer::etype_

private

Definition at line 106 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

std::string HcalRecHitsAnalyzer::eventype_

private

Definition at line 78 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

bool HcalRecHitsAnalyzer::famos_

private

Definition at line 82 of file HcalRecHitsAnalyzer.h.

edm::ESHandle<CaloGeometry> HcalRecHitsAnalyzer::geometry

private

Definition at line 258 of file HcalRecHitsAnalyzer.h.

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

std::vector<int> HcalRecHitsAnalyzer::hcalHBSevLvlVec

private

Definition at line 118 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHESevLvlVec

private

Definition at line 118 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHFSevLvlVec

private

Definition at line 118 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::vector<int> HcalRecHitsAnalyzer::hcalHOSevLvlVec

private

Definition at line 118 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and fillRecHitsTmp().

std::string HcalRecHitsAnalyzer::hcalselector_

private

Definition at line 76 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

const HcalDDDRecConstants* HcalRecHitsAnalyzer::hcons

private

Definition at line 84 of file HcalRecHitsAnalyzer.h.

Referenced by dqmBeginRun().

bool HcalRecHitsAnalyzer::hep17_

private

Definition at line 80 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::ieta_bins_

private

Definition at line 92 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::ieta_max_

private

Definition at line 93 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::ieta_min_

private

Definition at line 93 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::imc

private

Definition at line 108 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::iphi_bins_

private

Definition at line 89 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::iphi_max_

private

Definition at line 90 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float HcalRecHitsAnalyzer::iphi_min_

private

Definition at line 90 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::iz

private

Definition at line 107 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

MonitorElement* HcalRecHitsAnalyzer::map_ecal

private

Definition at line 136 of file HcalRecHitsAnalyzer.h.

int HcalRecHitsAnalyzer::maxDepthAll_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::maxDepthHB_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHE_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHF_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

int HcalRecHitsAnalyzer::maxDepthHO_

private

Definition at line 85 of file HcalRecHitsAnalyzer.h.

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

std::string HcalRecHitsAnalyzer::mc_

private

Definition at line 81 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meDeltaEta

private

Definition at line 174 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meDeltaPhi

private

Definition at line 173 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHB

private

Definition at line 245 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHE

private

Definition at line 246 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHF

private

Definition at line 248 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyConeHO

private

Definition at line 247 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHB

private

Definition at line 242 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEcalHcalEnergyHE

private

Definition at line 243 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile

private

Definition at line 169 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_E

private

Definition at line 170 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnConeEtaProfile_EH

private

Definition at line 171 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHB

private

Definition at line 251 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meEnergyHcalVsEcalHE

private

Definition at line 252 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHB

private

Definition at line 211 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meLog10Chi2profileHE

private

Definition at line 218 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHB

private

Definition at line 255 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meNumEcalRecHitsConeHE

private

Definition at line 256 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHB

private

Definition at line 183 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM0

private

Definition at line 184 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHBM3

private

Definition at line 185 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHE

private

Definition at line 191 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM0

private

Definition at line 192 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHEM3

private

Definition at line 193 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 194 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 195 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 196 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHF

private

Definition at line 204 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyHO

private

Definition at line 202 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HB

private

Definition at line 186 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM2vM0HE

private

Definition at line 197 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HB

private

Definition at line 187 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM0HE

private

Definition at line 198 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HB

private

Definition at line 188 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsEnergyM3vM2HE

private

Definition at line 199 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HB

private

Definition at line 189 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meRecHitsM2Chi2HE

private

Definition at line 200 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHB

private

Definition at line 235 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHE

private

Definition at line 236 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHF

private

Definition at line 238 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFL

private

Definition at line 239 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHFS

private

Definition at line 240 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyConeHO

private

Definition at line 237 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHB

private

Definition at line 230 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHE

private

Definition at line 231 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHF

private

Definition at line 233 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meSumRecHitsEnergyHO

private

Definition at line 232 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::meTE_HB

private

Definition at line 207 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HE

private

Definition at line 215 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HF

private

Definition at line 226 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HB

private

Definition at line 208 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_High_HO

private

Definition at line 221 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_HO

private

Definition at line 220 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HB

private

Definition at line 206 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HE

private

Definition at line 214 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTE_Low_HF

private

Definition at line 225 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB

private

Definition at line 210 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_High

private

Definition at line 212 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHB_Low

private

Definition at line 209 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE

private

Definition at line 217 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHE_Low

private

Definition at line 216 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF

private

Definition at line 228 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHF_Low

private

Definition at line 227 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO

private

Definition at line 222 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTEprofileHO_High

private

Definition at line 223 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHB

private

Definition at line 177 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHE

private

Definition at line 178 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHF

private

Definition at line 180 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::meTimeHO

private

Definition at line 179 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int HcalRecHitsAnalyzer::nChannels_[5]

private

Definition at line 87 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int HcalRecHitsAnalyzer::nevtot

private

Definition at line 294 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

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

private

Definition at line 126 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 127 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 129 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 128 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 163 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 163 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 164 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 165 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * HcalRecHitsAnalyzer::nrechits_vs_iphi_HOM

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::nrechits_vs_iphi_HOP

private

Definition at line 166 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 153 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 154 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 155 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::occupancy_map_HO

private

Definition at line 156 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 158 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

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

private

Definition at line 159 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

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

private

Definition at line 160 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::occupancy_vs_ieta_HO

private

Definition at line 161 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().

std::string HcalRecHitsAnalyzer::outputFile_

private

Definition at line 75 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HB

private

Definition at line 272 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HE

private

Definition at line 273 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HF

private

Definition at line 274 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_Aux_StatusWord_HO

private

Definition at line 275 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HB

private

Definition at line 261 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HE

private

Definition at line 262 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF

private

Definition at line 263 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HF67

private

Definition at line 264 of file HcalRecHitsAnalyzer.h.

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWord_HO

private

Definition at line 265 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HB

private

Definition at line 268 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::RecHit_StatusWordCorr_HE

private

Definition at line 269 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HB

private

Definition at line 120 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HE

private

Definition at line 121 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HF

private

Definition at line 122 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* HcalRecHitsAnalyzer::sevLvl_HO

private

Definition at line 123 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string HcalRecHitsAnalyzer::sign_

private

Definition at line 79 of file HcalRecHitsAnalyzer.h.

Referenced by HcalRecHitsAnalyzer().

int HcalRecHitsAnalyzer::subdet_

private

Definition at line 103 of file HcalRecHitsAnalyzer.h.

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

const HcalChannelQuality* HcalRecHitsAnalyzer::theHcalChStatus

private

Definition at line 113 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

const HcalSeverityLevelComputer* HcalRecHitsAnalyzer::theHcalSevLvlComputer

private

Definition at line 115 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

const HcalTopology* HcalRecHitsAnalyzer::theHcalTopology

private

Definition at line 111 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and hcalSevLvl().

edm::EDGetTokenT<EBRecHitCollection> HcalRecHitsAnalyzer::tok_EB_

private

Definition at line 99 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<EERecHitCollection> HcalRecHitsAnalyzer::tok_EE_

private

Definition at line 100 of file HcalRecHitsAnalyzer.h.

Referenced by analyze(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalRecHitsAnalyzer::tok_hbhe_

private

Definition at line 96 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HFRecHitCollection> HcalRecHitsAnalyzer::tok_hf_

private

Definition at line 98 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

edm::EDGetTokenT<HORecHitCollection> HcalRecHitsAnalyzer::tok_ho_

private

Definition at line 97 of file HcalRecHitsAnalyzer.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsAnalyzer().

std::string HcalRecHitsAnalyzer::topFolderName_

private

Definition at line 73 of file HcalRecHitsAnalyzer.h.

Referenced by bookHistograms().