HcalRecHitsValidation Class Reference

#include <HcalRecHitsValidation.h>

Inheritance diagram for HcalRecHitsValidation:

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
	HcalRecHitsValidation (edm::ParameterSet const &conf)
	~HcalRecHitsValidation () override=default
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
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	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)
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 >	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
const edm::InputTag	EBRecHitCollectionLabel_
const std::string	ecalselector_
const edm::InputTag	EERecHitCollectionLabel_
const CaloGeometry *	geometry_
const std::string	hcalselector_
int	imc
int	iz
const std::string	mc_
MonitorElement *	meEnConeEtaProfile
MonitorElement *	meEnConeEtaProfile_E
MonitorElement *	meEnConeEtaProfile_EH
MonitorElement *	meEnergyHcalVsEcalHB
MonitorElement *	meEnergyHcalVsEcalHE
MonitorElement *	meLog10Chi2profileHB
MonitorElement *	meLog10Chi2profileHE
MonitorElement *	meRecHitsEnergyHB
MonitorElement *	meRecHitsEnergyHE
MonitorElement *	meRecHitsEnergyHF
MonitorElement *	meRecHitsEnergyHO
MonitorElement *	meRecHitSimHitHB
MonitorElement *	meRecHitSimHitHE
MonitorElement *	meRecHitSimHitHF
MonitorElement *	meRecHitSimHitHFL
MonitorElement *	meRecHitSimHitHFS
MonitorElement *	meRecHitSimHitHO
MonitorElement *	meRecHitSimHitProfileHB
MonitorElement *	meRecHitSimHitProfileHE
MonitorElement *	meRecHitSimHitProfileHF
MonitorElement *	meRecHitSimHitProfileHFL
MonitorElement *	meRecHitSimHitProfileHFS
MonitorElement *	meRecHitSimHitProfileHO
MonitorElement *	meRecHitsM2Chi2HB
MonitorElement *	meRecHitsM2Chi2HE
MonitorElement *	meTEprofileHB
MonitorElement *	meTEprofileHB_High
MonitorElement *	meTEprofileHB_Low
MonitorElement *	meTEprofileHE
MonitorElement *	meTEprofileHE_Low
MonitorElement *	meTEprofileHF
MonitorElement *	meTEprofileHF_Low
MonitorElement *	meTEprofileHO
MonitorElement *	meTEprofileHO_High
int	nevtot
const std::string	outputFile_
const std::string	sign_
int	subdet_
const bool	testNumber_
const edm::EDGetTokenT< EBRecHitCollection >	tok_EB_
const edm::EDGetTokenT< EERecHitCollection >	tok_EE_
const edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_Geom_
const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_
const edm::EDGetTokenT< HFRecHitCollection >	tok_hf_
const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hh_
const edm::EDGetTokenT< HORecHitCollection >	tok_ho_
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_HRNDC_
const 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 57 of file HcalRecHitsValidation.h.

Constructor & Destructor Documentation

HcalRecHitsValidation()

HcalRecHitsValidation::HcalRecHitsValidation

(

edm::ParameterSet const &

conf

)

Definition at line 8 of file HcalRecHitsValidation.cc.

References hcalselector_, imc, iz, mc_, nevtot, outputFile_, sign_, and subdet_.

    : topFolderName_(conf.getParameter<std::string>("TopFolderName")),
      outputFile_(conf.getUntrackedParameter<std::string>("outputFile", "myfile.root")),
      hcalselector_(conf.getUntrackedParameter<std::string>("hcalselector", "all")),
      ecalselector_(conf.getUntrackedParameter<std::string>("ecalselector", "yes")),
      sign_(conf.getUntrackedParameter<std::string>("sign", "*")),
      mc_(conf.getUntrackedParameter<std::string>("mc", "yes")),
      testNumber_(conf.getParameter<bool>("TestNumber")),
      EBRecHitCollectionLabel_(conf.getParameter<edm::InputTag>("EBRecHitCollectionLabel")),
      EERecHitCollectionLabel_(conf.getParameter<edm::InputTag>("EERecHitCollectionLabel")),
      tok_evt_(consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"))),
      tok_EB_(consumes<EBRecHitCollection>(EBRecHitCollectionLabel_)),
      tok_EE_(consumes<EERecHitCollection>(EERecHitCollectionLabel_)),
      tok_hh_(consumes<edm::PCaloHitContainer>(conf.getUntrackedParameter<edm::InputTag>("SimHitCollectionLabel"))),
      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"))),
      tok_HRNDC_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord>()),
      tok_Geom_(esConsumes<CaloGeometry, CaloGeometryRecord>()) {
  // DQM ROOT output
  if (!outputFile_.empty()) {
    edm::LogVerbatim("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
  } else {
    edm::LogVerbatim("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
  }

  nevtot = 0;

  // Collections

  // register for data access

  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;

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

  imc = 1;
  if (mc_ == "no")
    imc = 0;
}

~HcalRecHitsValidation()

HcalRecHitsValidation::~HcalRecHitsValidation ( )

overridedefault

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 186 of file HcalRecHitsValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), HltBtagPostValidation_cff::c, cchi2, cdepth, cen, ceta, hltPixelTracks_cff::chi2, cieta, cphi, csub, ctime, LEDCalibrationChannels::depth, dR(), ecalselector_, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), ALCARECOTkAlBeamHalo_cff::etaMax, makeMEIFBenchmarkPlots::ev, dqm::impl::MonitorElement::Fill(), fillRecHitsTmp(), geometry_, edm::HepMCProduct::GetEvent(), CaloSubdetectorGeometry::getGeometry(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, LEDCalibrationChannels::ieta, imc, edm::HandleBase::isValid(), PV_cfg::maxPt, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meLog10Chi2profileHB, meLog10Chi2profileHE, meRecHitsEnergyHB, meRecHitsEnergyHE, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitSimHitHB, meRecHitSimHitHE, meRecHitSimHitHF, meRecHitSimHitHFL, meRecHitSimHitHFS, meRecHitSimHitHO, meRecHitSimHitProfileHB, meRecHitSimHitProfileHE, meRecHitSimHitProfileHF, meRecHitSimHitProfileHFL, meRecHitSimHitProfileHFS, meRecHitSimHitProfileHO, meRecHitsM2Chi2HB, meRecHitsM2Chi2HE, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, nevtot, npart, AlCaHLTBitMon_ParallelJobs::p, phi, PV3DBase< T, PVType, FrameType >::phi(), edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, alignCSCRings::r, HcalHitRelabeller::relabel(), HcalDetId::subdet(), subdet_, submitPVValidationJobs::t, testNumber_, tok_EB_, tok_EE_, tok_evt_, tok_Geom_, tok_hh_, and tok_HRNDC_.

                                                                              {
  using namespace edm;

  const HcalDDDRecConstants *hcons = &c.getData(tok_HRNDC_);

  // 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 eHcalConeHB = 0.;
  double eHcalConeHE = 0.;
  double eHcalConeHO = 0.;
  double eHcalConeHF = 0.;
  double eHcalConeHFL = 0.;
  double eHcalConeHFS = 0.;

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

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

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

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

  if (imc != 0) {
    const edm::Handle<edm::HepMCProduct> &evtMC = ev.getHandle(tok_evt_);  // generator in late 310_preX
    if (!evtMC.isValid()) {
      edm::LogVerbatim("HcalRecHitsValidation") << "no HepMCProduct found";
    } else {
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HcalRecHitsValidation") << "*** source HepMCProduct found";
#endif
    }

    // MC particle with highest pt is taken as a direction reference
    double maxPt = -99999.;
    int npart = 0;
    const HepMC::GenEvent *myGenEvent = evtMC->GetEvent();
    for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
         ++p) {
      double phip = (*p)->momentum().phi();
      double etap = (*p)->momentum().eta();
      //    phi_MC = phip;
      //    eta_MC = etap;
      double pt = (*p)->momentum().perp();
      if (pt > maxPt) {
        npart++;
        maxPt = pt;
        phi_MC = phip;
        eta_MC = etap;
      }
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HcalRecHitsValidation") << "*** Max pT = " << maxPt;
#endif
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalRecHitsValidation") << "*** 2";
#endif

  geometry_ = &c.getData(tok_Geom_);

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

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalRecHitsValidation") << "*** 3";
#endif

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

  // ECAL
  if (ecalselector_ == "yes" && (subdet_ == 1 || subdet_ == 2 || subdet_ == 5)) {
    const edm::Handle<EBRecHitCollection> &rhitEB = ev.getHandle(tok_EB_);

    EcalRecHitCollection::const_iterator RecHit;
    EcalRecHitCollection::const_iterator RecHitEnd;

    if (rhitEB.isValid()) {
      RecHit = rhitEB.product()->begin();
      RecHitEnd = rhitEB.product()->end();

      for (; RecHit != RecHitEnd; ++RecHit) {
        EBDetId EBid = EBDetId(RecHit->id());

        auto cellGeometry = geometry_->getSubdetectorGeometry(EBid)->getGeometry(EBid);
        double eta = cellGeometry->getPosition().eta();
        double phi = cellGeometry->getPosition().phi();
        double en = RecHit->energy();

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

    const edm::Handle<EERecHitCollection> &rhitEE = ev.getHandle(tok_EE_);

    if (rhitEE.isValid()) {
      RecHit = rhitEE.product()->begin();
      RecHitEnd = rhitEE.product()->end();

      for (; RecHit != RecHitEnd; ++RecHit) {
        EEDetId EEid = EEDetId(RecHit->id());

        auto cellGeometry = geometry_->getSubdetectorGeometry(EEid)->getGeometry(EEid);
        double eta = cellGeometry->getPosition().eta();
        double phi = cellGeometry->getPosition().phi();
        double en = RecHit->energy();

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

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalRecHitsValidation") << "*** 4";
#endif

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

  if ((subdet_ != 6) && (subdet_ != 0)) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HcalRecHitsValidation") << "*** 6";
#endif

    double HcalCone = 0.;

    int ietaMax = 9999;
    double etaMax = 9999.;

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

    for (unsigned int i = 0; i < cen.size(); i++) {
      int sub = csub[i];
      int depth = cdepth[i];
      double eta = ceta[i];
      double phi = cphi[i];
      double en = cen[i];
      double t = ctime[i];
      int ieta = cieta[i];
      double chi2 = cchi2[i];

      double chi2_log10 = 9.99;  // initial value above histos limits , keep it if chi2 <= 0.
      if (chi2 > 0.)
        chi2_log10 = log10(chi2);

      nrechits++;
      if (en > 1.)
        nrechitsThresh++;

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

        HcalCone += en;
        // alternative: ietamax -> closest to MC eta  !!!
        float eta_diff = fabs(eta_MC - eta);
        if (eta_diff < etaMax) {
          etaMax = eta_diff;
          ietaMax = ieta;
        }
      }

      // The energy and overall timing histos are drawn while
      // the ones split by depth are not
      if (sub == 1 && (subdet_ == 1 || subdet_ == 5)) {
        meRecHitsM2Chi2HB->Fill(chi2_log10);
        meLog10Chi2profileHB->Fill(en, chi2_log10);

        meRecHitsEnergyHB->Fill(en);

        meTEprofileHB_Low->Fill(en, t);
        meTEprofileHB->Fill(en, t);
        meTEprofileHB_High->Fill(en, t);
      }
      if (sub == 2 && (subdet_ == 2 || subdet_ == 5)) {
        meRecHitsM2Chi2HE->Fill(chi2_log10);
        meLog10Chi2profileHE->Fill(en, chi2_log10);

        meRecHitsEnergyHE->Fill(en);

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

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

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

    if (imc != 0) {
      meEnConeEtaProfile->Fill(double(ietaMax), HcalCone);  //
      meEnConeEtaProfile_E->Fill(double(ietaMax), eEcalCone);
      meEnConeEtaProfile_EH->Fill(double(ietaMax), HcalCone + eEcalCone);
    }

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HcalRecHitsValidation") << "*** 7";
#endif
  }

  // SimHits vs. RecHits
  if (subdet_ > 0 && subdet_ < 6 && imc != 0) {  // not noise

    const edm::Handle<PCaloHitContainer> &hcalHits = ev.getHandle(tok_hh_);
    if (hcalHits.isValid()) {
      const PCaloHitContainer *SimHitResult = hcalHits.product();

      double enSimHitsHB = 0.;
      double enSimHitsHE = 0.;
      double enSimHitsHO = 0.;
      double enSimHitsHF = 0.;
      double enSimHitsHFL = 0.;
      double enSimHitsHFS = 0.;
      // sum of SimHits in the cone

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

        if (testNumber_)
          cell = HcalHitRelabeller::relabel(SimHits->id(), hcons);
        else
          cell = HcalDetId(SimHits->id());

        sub = cell.subdet();
        depth = cell.depth();

        if (sub != subdet_ && subdet_ != 5)
          continue;  // If we are not looking at all of the subdetectors and the
                     // simhit doesn't come from the specific subdetector of
                     // interest, then we won't do any thing with it

        const HcalGeometry *cellGeometry =
            dynamic_cast<const HcalGeometry *>(geometry_->getSubdetectorGeometry(DetId::Hcal, cell.subdet()));
        double etaS = cellGeometry->getPosition(cell).eta();
        double phiS = cellGeometry->getPosition(cell).phi();
        double en = SimHits->energy();

        double r = dR(eta_MC, phi_MC, etaS, phiS);

        if (r < partR) {  // just energy in the small cone

          if (sub == static_cast<int>(HcalBarrel))
            enSimHitsHB += en;
          if (sub == static_cast<int>(HcalEndcap))
            enSimHitsHE += en;
          if (sub == static_cast<int>(HcalOuter))
            enSimHitsHO += en;
          if (sub == static_cast<int>(HcalForward)) {
            enSimHitsHF += en;
            if (depth == 1)
              enSimHitsHFL += en;
            else
              enSimHitsHFS += en;
          }
        }
      }

      // Now some histos with SimHits

      if (subdet_ == 4 || subdet_ == 5) {
        meRecHitSimHitHF->Fill(enSimHitsHF, eHcalConeHF);
        meRecHitSimHitProfileHF->Fill(enSimHitsHF, eHcalConeHF);

        meRecHitSimHitHFL->Fill(enSimHitsHFL, eHcalConeHFL);
        meRecHitSimHitProfileHFL->Fill(enSimHitsHFL, eHcalConeHFL);
        meRecHitSimHitHFS->Fill(enSimHitsHFS, eHcalConeHFS);
        meRecHitSimHitProfileHFS->Fill(enSimHitsHFS, eHcalConeHFS);
      }
      if (subdet_ == 1 || subdet_ == 5) {
        meRecHitSimHitHB->Fill(enSimHitsHB, eHcalConeHB);
        meRecHitSimHitProfileHB->Fill(enSimHitsHB, eHcalConeHB);
      }
      if (subdet_ == 2 || subdet_ == 5) {
        meRecHitSimHitHE->Fill(enSimHitsHE, eHcalConeHE);
        meRecHitSimHitProfileHE->Fill(enSimHitsHE, eHcalConeHE);
      }
      if (subdet_ == 3 || subdet_ == 5) {
        meRecHitSimHitHO->Fill(enSimHitsHO, eHcalConeHO);
        meRecHitSimHitProfileHO->Fill(enSimHitsHO, eHcalConeHO);
      }
    }
  }

  nevtot++;
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 67 of file HcalRecHitsValidation.cc.

References timingPdfMaker::histo, cuy::ib, imc, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meLog10Chi2profileHB, meLog10Chi2profileHE, meRecHitsEnergyHB, meRecHitsEnergyHE, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitSimHitHB, meRecHitSimHitHE, meRecHitSimHitHF, meRecHitSimHitHFL, meRecHitSimHitHFS, meRecHitSimHitHO, meRecHitSimHitProfileHB, meRecHitSimHitProfileHE, meRecHitSimHitProfileHF, meRecHitSimHitProfileHFL, meRecHitSimHitProfileHFS, meRecHitSimHitProfileHO, meRecHitsM2Chi2HB, meRecHitsM2Chi2HE, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, subdet_, and topFolderName_.

                                                                                                          {
  Char_t histo[200];

  ib.setCurrentFolder(topFolderName_);

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

  // Histograms drawn for single pion scan
  if (subdet_ != 0 && imc != 0) {  // just not for noise
    sprintf(histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths");
    meEnConeEtaProfile = ib.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");

    sprintf(histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
    meEnConeEtaProfile_E = ib.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");

    sprintf(histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
    meEnConeEtaProfile_EH = ib.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");
  }

  // ************** HB **********************************
  if (subdet_ == 1 || subdet_ == 5) {
    sprintf(histo, "HcalRecHitTask_M2Log10Chi2_of_rechits_HB");  //   Chi2
    meRecHitsM2Chi2HB = ib.book1D(histo, histo, 120, -2., 10.);

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

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HB");
    meRecHitsEnergyHB = ib.book1D(histo, histo, 2010, -10., 2000.);

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

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HB");
    meTEprofileHB_Low = ib.bookProfile(histo, histo, 150, -5., 295., -48., 92., " ");

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

    if (imc != 0) {
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_HB");
      meRecHitSimHitHB = ib.book2D(histo, histo, 120, 0., 1.2, 300, 0., 150.);
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HB");
      meRecHitSimHitProfileHB = ib.bookProfile(histo, histo, 120, 0., 1.2, 0., 500., " ");
    }
  }

  // ********************** HE ************************************
  if (subdet_ == 2 || subdet_ == 5) {
    sprintf(histo, "HcalRecHitTask_M2Log10Chi2_of_rechits_HE");  //   Chi2
    meRecHitsM2Chi2HE = ib.book1D(histo, histo, 120, -2., 10.);

    sprintf(histo, "HcalRecHitTask_Log10Chi2_vs_energy_profile_HE");
    meLog10Chi2profileHE = ib.bookProfile(histo, histo, 1000, -5., 995., -2., 9.9, " ");

    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HE");
    meRecHitsEnergyHE = ib.book1D(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HE");
    meTEprofileHE_Low = ib.bookProfile(histo, histo, 80, -5., 75., -48., 92., " ");

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

    if (imc != 0) {
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_HE");
      meRecHitSimHitHE = ib.book2D(histo, histo, 120, 0., 0.6, 300, 0., 150.);
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HE");
      meRecHitSimHitProfileHE = ib.bookProfile(histo, histo, 120, 0., 0.6, 0., 500., " ");
    }
  }

  // ************** HO ****************************************
  if (subdet_ == 3 || subdet_ == 5) {
    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HO");
    meRecHitsEnergyHO = ib.book1D(histo, histo, 2010, -10., 2000.);

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

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

    if (imc != 0) {
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_HO");
      meRecHitSimHitHO = ib.book2D(histo, histo, 150, 0., 1.5, 350, 0., 350.);
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HO");
      meRecHitSimHitProfileHO = ib.bookProfile(histo, histo, 150, 0., 1.5, 0., 500., " ");
    }
  }

  // ********************** HF ************************************
  if (subdet_ == 4 || subdet_ == 5) {
    sprintf(histo, "HcalRecHitTask_energy_of_rechits_HF");
    meRecHitsEnergyHF = ib.book1D(histo, histo, 2010, -10., 2000.);

    sprintf(histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HF");
    meTEprofileHF_Low = ib.bookProfile(histo, histo, 100, -5., 195., -48., 92., " ");

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

    if (imc != 0) {
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_HF");
      meRecHitSimHitHF = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_HFL");
      meRecHitSimHitHFL = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_HFS");
      meRecHitSimHitHFS = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HF");
      meRecHitSimHitProfileHF = ib.bookProfile(histo, histo, 50, 0., 50., 0., 500., " ");
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HFL");
      meRecHitSimHitProfileHFL = ib.bookProfile(histo, histo, 50, 0., 50., 0., 500., " ");
      sprintf(histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HFS");
      meRecHitSimHitProfileHFS = ib.bookProfile(histo, histo, 50, 0., 50., 0., 500., " ");
    }
  }
}

dPhiWsign()

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

private

Definition at line 669 of file HcalRecHitsValidation.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;
}

dR()

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

private

Definition at line 636 of file HcalRecHitsValidation.cc.

References HLT_2023v12_cff::eta1, HLT_2023v12_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 HcalRecHitsValidation::fillRecHitsTmp ( int  subdet_, edm::Event const &  ev  )

privatevirtual

Definition at line 517 of file HcalRecHitsValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), cchi2, cdepth, cen, ceta, hltPixelTracks_cff::chi2, cieta, ciphi, cphi, csub, ctime, cz, LEDCalibrationChannels::depth, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, makeMEIFBenchmarkPlots::ev, geometry_, CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, edm::HandleBase::isValid(), iz, dqmiolumiharvest::j, phi, subdet_, submitPVValidationJobs::t, tok_hbhe_, tok_hf_, and tok_ho_.

Referenced by analyze().

                                                                          {
  using namespace edm;

  // initialize data vectors
  csub.clear();
  cen.clear();
  ceta.clear();
  cphi.clear();
  ctime.clear();
  cieta.clear();
  ciphi.clear();
  cdepth.clear();
  cz.clear();
  cchi2.clear();

  if (subdet_ == 1 || subdet_ == 2 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    // HBHE
    const edm::Handle<HBHERecHitCollection> &hbhecoll = ev.getHandle(tok_hbhe_);
    if (hbhecoll.isValid()) {
      for (HBHERecHitCollection::const_iterator j = hbhecoll->begin(); j != hbhecoll->end(); j++) {
        HcalDetId cell(j->id());
        const HcalGeometry *cellGeometry =
            dynamic_cast<const HcalGeometry *>(geometry_->getSubdetectorGeometry(DetId::Hcal, cell.subdet()));
        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() - 1;
        double en = j->energy();
        double t = j->time();
        double chi2 = j->chi2();

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

  if (subdet_ == 4 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    // HF
    const edm::Handle<HFRecHitCollection> &hfcoll = ev.getHandle(tok_hf_);
    if (hfcoll.isValid()) {
      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() - 1;
        double en = j->energy();
        double t = j->time();

        if ((iz > 0 && eta > 0.) || (iz < 0 && eta < 0.) || iz == 0) {
          csub.push_back(sub);
          cen.push_back(en);
          ceta.push_back(eta);
          cphi.push_back(phi);
          ctime.push_back(t);
          cieta.push_back(inteta);
          ciphi.push_back(intphi);
          cdepth.push_back(depth);
          cz.push_back(zc);
          cchi2.push_back(0.);
        }
      }
    }
  }

  // HO
  if (subdet_ == 3 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
    const edm::Handle<HORecHitCollection> &hocoll = ev.getHandle(tok_ho_);
    if (hocoll.isValid()) {
      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() - 1;
        double t = j->time();
        double en = j->energy();

        if ((iz > 0 && eta > 0.) || (iz < 0 && eta < 0.) || iz == 0) {
          csub.push_back(sub);
          cen.push_back(en);
          ceta.push_back(eta);
          cphi.push_back(phi);
          ctime.push_back(t);
          cieta.push_back(inteta);
          ciphi.push_back(intphi);
          cdepth.push_back(depth);
          cz.push_back(zc);
          cchi2.push_back(0.);
        }
      }
    }
  }
}

phi12()

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

private

Definition at line 650 of file HcalRecHitsValidation.cc.

References testProducerWithPsetDescEmpty_cfi::a2, PI, and createJobs::tmp.

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

cauxstwd

std::vector<uint32_t> HcalRecHitsValidation::cauxstwd

private

Definition at line 167 of file HcalRecHitsValidation.h.

cchi2

std::vector<double> HcalRecHitsValidation::cchi2

private

Definition at line 168 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

cdepth

std::vector<int> HcalRecHitsValidation::cdepth

private

Definition at line 160 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

cen

std::vector<double> HcalRecHitsValidation::cen

private

Definition at line 161 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

ceta

std::vector<double> HcalRecHitsValidation::ceta

private

Definition at line 162 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

cieta

std::vector<int> HcalRecHitsValidation::cieta

private

Definition at line 158 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

ciphi

std::vector<int> HcalRecHitsValidation::ciphi

private

Definition at line 159 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp().

cphi

std::vector<double> HcalRecHitsValidation::cphi

private

Definition at line 163 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

cstwd

std::vector<uint32_t> HcalRecHitsValidation::cstwd

private

Definition at line 166 of file HcalRecHitsValidation.h.

csub

std::vector<int> HcalRecHitsValidation::csub

private

Definition at line 157 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

ctime

std::vector<double> HcalRecHitsValidation::ctime

private

Definition at line 164 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

cz

std::vector<double> HcalRecHitsValidation::cz

private

Definition at line 165 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp().

EBRecHitCollectionLabel_

const edm::InputTag HcalRecHitsValidation::EBRecHitCollectionLabel_

private

Definition at line 80 of file HcalRecHitsValidation.h.

ecalselector_

const std::string HcalRecHitsValidation::ecalselector_

private

Definition at line 74 of file HcalRecHitsValidation.h.

Referenced by analyze().

EERecHitCollectionLabel_

const edm::InputTag HcalRecHitsValidation::EERecHitCollectionLabel_

private

Definition at line 81 of file HcalRecHitsValidation.h.

geometry_

const CaloGeometry* HcalRecHitsValidation::geometry_

private

Definition at line 154 of file HcalRecHitsValidation.h.

Referenced by analyze(), and fillRecHitsTmp().

hcalselector_

const std::string HcalRecHitsValidation::hcalselector_

private

Definition at line 73 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

imc

int HcalRecHitsValidation::imc

private

Definition at line 98 of file HcalRecHitsValidation.h.

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

iz

int HcalRecHitsValidation::iz

private

Definition at line 97 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp(), and HcalRecHitsValidation().

mc_

const std::string HcalRecHitsValidation::mc_

private

Definition at line 76 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

meEnConeEtaProfile

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile

private

Definition at line 105 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meEnConeEtaProfile_E

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_E

private

Definition at line 106 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meEnConeEtaProfile_EH

MonitorElement* HcalRecHitsValidation::meEnConeEtaProfile_EH

private

Definition at line 107 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meEnergyHcalVsEcalHB

MonitorElement* HcalRecHitsValidation::meEnergyHcalVsEcalHB

private

Definition at line 144 of file HcalRecHitsValidation.h.

meEnergyHcalVsEcalHE

MonitorElement* HcalRecHitsValidation::meEnergyHcalVsEcalHE

private

Definition at line 145 of file HcalRecHitsValidation.h.

meLog10Chi2profileHB

MonitorElement* HcalRecHitsValidation::meLog10Chi2profileHB

private

Definition at line 151 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meLog10Chi2profileHE

MonitorElement* HcalRecHitsValidation::meLog10Chi2profileHE

private

Definition at line 152 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHB

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHB

private

Definition at line 110 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHE

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHE

private

Definition at line 111 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHF

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHF

private

Definition at line 113 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitsEnergyHO

MonitorElement* HcalRecHitsValidation::meRecHitsEnergyHO

private

Definition at line 112 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitHB

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHB

private

Definition at line 129 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitHE

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHE

private

Definition at line 130 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitHF

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHF

private

Definition at line 132 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitHFL

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHFL

private

Definition at line 133 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitHFS

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHFS

private

Definition at line 134 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitHO

MonitorElement* HcalRecHitsValidation::meRecHitSimHitHO

private

Definition at line 131 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitProfileHB

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHB

private

Definition at line 136 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitProfileHE

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHE

private

Definition at line 137 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitProfileHF

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHF

private

Definition at line 139 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitProfileHFL

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHFL

private

Definition at line 140 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitProfileHFS

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHFS

private

Definition at line 141 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitSimHitProfileHO

MonitorElement* HcalRecHitsValidation::meRecHitSimHitProfileHO

private

Definition at line 138 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitsM2Chi2HB

MonitorElement* HcalRecHitsValidation::meRecHitsM2Chi2HB

private

Definition at line 148 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meRecHitsM2Chi2HE

MonitorElement* HcalRecHitsValidation::meRecHitsM2Chi2HE

private

Definition at line 149 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHB

MonitorElement* HcalRecHitsValidation::meTEprofileHB

private

Definition at line 116 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHB_High

MonitorElement* HcalRecHitsValidation::meTEprofileHB_High

private

Definition at line 117 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHB_Low

MonitorElement* HcalRecHitsValidation::meTEprofileHB_Low

private

Definition at line 115 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHE

MonitorElement* HcalRecHitsValidation::meTEprofileHE

private

Definition at line 120 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHE_Low

MonitorElement* HcalRecHitsValidation::meTEprofileHE_Low

private

Definition at line 119 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHF

MonitorElement* HcalRecHitsValidation::meTEprofileHF

private

Definition at line 126 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHF_Low

MonitorElement* HcalRecHitsValidation::meTEprofileHF_Low

private

Definition at line 125 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHO

MonitorElement* HcalRecHitsValidation::meTEprofileHO

private

Definition at line 122 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meTEprofileHO_High

MonitorElement* HcalRecHitsValidation::meTEprofileHO_High

private

Definition at line 123 of file HcalRecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

nevtot

int HcalRecHitsValidation::nevtot

private

Definition at line 171 of file HcalRecHitsValidation.h.

Referenced by analyze(), and HcalRecHitsValidation().

outputFile_

const std::string HcalRecHitsValidation::outputFile_

private

Definition at line 72 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

sign_

const std::string HcalRecHitsValidation::sign_

private

Definition at line 75 of file HcalRecHitsValidation.h.

Referenced by HcalRecHitsValidation().

subdet_

int HcalRecHitsValidation::subdet_

private

Definition at line 94 of file HcalRecHitsValidation.h.

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

testNumber_

const bool HcalRecHitsValidation::testNumber_

private

Definition at line 77 of file HcalRecHitsValidation.h.

Referenced by analyze().

tok_EB_

const edm::EDGetTokenT<EBRecHitCollection> HcalRecHitsValidation::tok_EB_

private

Definition at line 83 of file HcalRecHitsValidation.h.

Referenced by analyze().

tok_EE_

const edm::EDGetTokenT<EERecHitCollection> HcalRecHitsValidation::tok_EE_

private

Definition at line 84 of file HcalRecHitsValidation.h.

Referenced by analyze().

tok_evt_

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

private

Definition at line 82 of file HcalRecHitsValidation.h.

Referenced by analyze().

tok_Geom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalRecHitsValidation::tok_Geom_

private

Definition at line 91 of file HcalRecHitsValidation.h.

Referenced by analyze().

tok_hbhe_

const edm::EDGetTokenT<HBHERecHitCollection> HcalRecHitsValidation::tok_hbhe_

private

Definition at line 86 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp().

tok_hf_

const edm::EDGetTokenT<HFRecHitCollection> HcalRecHitsValidation::tok_hf_

private

Definition at line 87 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp().

tok_hh_

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

private

Definition at line 85 of file HcalRecHitsValidation.h.

Referenced by analyze().

tok_ho_

const edm::EDGetTokenT<HORecHitCollection> HcalRecHitsValidation::tok_ho_

private

Definition at line 88 of file HcalRecHitsValidation.h.

Referenced by fillRecHitsTmp().

tok_HRNDC_

const edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalRecHitsValidation::tok_HRNDC_

private

Definition at line 90 of file HcalRecHitsValidation.h.

Referenced by analyze().

topFolderName_

const std::string HcalRecHitsValidation::topFolderName_

private

Definition at line 71 of file HcalRecHitsValidation.h.

Referenced by bookHistograms().