HcalSimHitCheck Class Reference

Inheritance diagram for HcalSimHitCheck:

Public Member Functions
	HcalSimHitCheck (const edm::ParameterSet &ps)
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	analyzeHits (std::vector< PCaloHit > &)
void	beginRun (edm::Run const &, edm::EventSetup const &) override
void	endRun (edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Private Attributes
const bool	checkHit_
const std::string	g4Label
const std::string	hcalHits
const HcalDDDRecConstants *	hcons_
const bool	hep17_
int	ieta_bins_HB
int	ieta_bins_HE
int	ieta_bins_HF
int	ieta_bins_HO
float	ieta_max_HB
float	ieta_max_HE
float	ieta_max_HF
float	ieta_max_HO
float	ieta_min_HB
float	ieta_min_HE
float	ieta_min_HF
float	ieta_min_HO
int	iphi_bins
float	iphi_max
float	iphi_min
int	maxDepth_
int	maxDepthHB_
int	maxDepthHE_
int	maxDepthHF_
int	maxDepthHO_
TH1D *	meAllNHit_
TH1D *	meBadDetHit_
TH1D *	meBadIdHit_
TH1D *	meBadSubHit_
TH1D *	meDepthHit_
TH1D *	meDetectHit_
TH1D *	meEnergyHit_
TH1D *	meEtaHit_
TH2D *	meEtaPhiHit_
std::vector< TH2D * >	meEtaPhiHitDepth_
TH1D *	meHBDepHit_
TH1D *	meHBEneHit2_
TH1D *	meHBEneHit_
TH2D *	meHBEneMap_
TH1D *	meHBEneSum_
TProfile *	meHBEneSum_vs_ieta_
TH1D *	meHBEtaHit_
TH1D *	meHBL10Ene_
TProfile *	meHBL10EneP_
TH1D *	meHBNHit_
TH1D *	meHBPhiHit_
TH1D *	meHBTimHit_
TH1D *	meHEDepHit_
TH1D *	meHEEneHit2_
TH1D *	meHEEneHit_
TH2D *	meHEEneMap_
TH1D *	meHEEneSum_
TProfile *	meHEEneSum_vs_ieta_
TH1D *	meHEEtaHit_
TH1D *	meHEL10Ene_
TProfile *	meHEL10EneP_
TH1D *	meHENHit_
TH1D *	meHEP17EneHit2_
TH1D *	meHEP17EneHit_
TH1D *	meHEPhiHit_
TH1D *	meHETimHit_
TH1D *	meHFDepHit_
TH1D *	meHFDepHitw_
TH1D *	meHFEneHit2_
TH1D *	meHFEneHit_
TH2D *	meHFEneMap_
TH1D *	meHFEneSum_
TProfile *	meHFEneSum_vs_ieta_
TH1D *	meHFEtaHit_
TH1D *	meHFL10Ene_
TProfile *	meHFL10EneP_
TH1D *	meHFNHit_
TH1D *	meHFPhiHit_
TH1D *	meHFTimHit_
TH1D *	meHODepHit_
TH1D *	meHOEneHit2_
TH1D *	meHOEneHit_
TH2D *	meHOEneMap_
TH1D *	meHOEneSum_
TProfile *	meHOEneSum_vs_ieta_
TH1D *	meHOEtaHit_
TH1D *	meHOL10Ene_
TProfile *	meHOL10EneP_
TH1D *	meHONHit_
TH1D *	meHOPhiHit_
TH1D *	meHOTimHit_
TH1D *	mePhiHit_
TH1D *	mePhiHitb_
TH1D *	meSubdetHit_
TH1D *	meTimeHit_
TH1D *	meTimeWHit_
const std::string	outFile_
const bool	testNumber_
const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hits_
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_HRNDC_
const int	verbose_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 34 of file HcalSimHitCheck.cc.

Constructor & Destructor Documentation

HcalSimHitCheck()

HcalSimHitCheck::HcalSimHitCheck

(

const edm::ParameterSet &

ps

)

Definition at line 91 of file HcalSimHitCheck.cc.

References checkHit_, g4Label, hcalHits, and outFile_.

    : g4Label(ps.getParameter<std::string>("moduleLabel")),
      hcalHits(ps.getParameter<std::string>("HitCollection")),
      outFile_(ps.getParameter<std::string>("outputFile")),
      verbose_(ps.getParameter<int>("Verbose")),
      checkHit_(true),
      testNumber_(ps.getParameter<bool>("TestNumber")),
      hep17_(ps.getParameter<bool>("hep17")),
      tok_hits_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label, hcalHits))),
      tok_HRNDC_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()) {
  edm::LogVerbatim("HcalSim") << "Module Label: " << g4Label << "   Hits: " << hcalHits << " / " << checkHit_
                              << "   Output: " << outFile_;
}

Member Function Documentation

analyze()

void HcalSimHitCheck::analyze ( edm::Event const &  e, edm::EventSetup const &    )

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 303 of file HcalSimHitCheck.cc.

References analyzeHits(), checkHit_, MillePedeFileConverter_cfg::e, edm::HandleBase::isValid(), tok_hits_, and verbose_.

                                                                      {
  if (verbose_ > 0)
    edm::LogVerbatim("HcalSim") << "Run = " << e.id().run() << " Event = " << e.id().event();

  bool getHits = false;
  if (checkHit_) {
    const edm::Handle<edm::PCaloHitContainer> &hitsHcal = e.getHandle(tok_hits_);
    if (hitsHcal.isValid()) {
      getHits = true;
      if (verbose_ > 0)
        edm::LogVerbatim("HcalSim") << "HcalValidation: Input flags Hits " << getHits;
      std::vector<PCaloHit> caloHits;
      caloHits.insert(caloHits.end(), hitsHcal->begin(), hitsHcal->end());
      if (verbose_ > 0)
        edm::LogVerbatim("HcalSim") << "HcalValidation: Hit buffer " << caloHits.size();
      analyzeHits(caloHits);
    } else {
      if (verbose_ > 0)
        edm::LogVerbatim("HcalSim") << "HcalValidation: Input flags Hits " << getHits;
    }
  }
}

analyzeHits()

void HcalSimHitCheck::analyzeHits ( std::vector< PCaloHit > &  hits )

protected

Definition at line 326 of file HcalSimHitCheck.cc.

References TauDecayModes::dec, hcalRecHitTable_cff::depth, HcalDetId::depth(), DetId::det(), hcalRecHitTable_cff::energy, PVValHelper::eta, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons_, hep17_, hfClusterShapes_cfi::hits, mps_fire::i, HcalDetId::ieta(), ieta_bins_HB, ieta_bins_HE, ieta_bins_HF, ieta_bins_HO, ieta_min_HB, ieta_min_HE, ieta_min_HF, ieta_min_HO, createfilelist::int, HcalDetId::iphi(), iphi_bins, dqmiolumiharvest::j, meAllNHit_, meBadDetHit_, meBadIdHit_, meBadSubHit_, meDepthHit_, meDetectHit_, meEnergyHit_, meEtaHit_, meEtaPhiHit_, meEtaPhiHitDepth_, meHBDepHit_, meHBEneHit2_, meHBEneHit_, meHBEneMap_, meHBEneSum_, meHBEneSum_vs_ieta_, meHBEtaHit_, meHBL10Ene_, meHBL10EneP_, meHBNHit_, meHBPhiHit_, meHBTimHit_, meHEDepHit_, meHEEneHit2_, meHEEneHit_, meHEEneMap_, meHEEneSum_, meHEEneSum_vs_ieta_, meHEEtaHit_, meHEL10Ene_, meHEL10EneP_, meHENHit_, meHEP17EneHit2_, meHEP17EneHit_, meHEPhiHit_, meHETimHit_, meHFDepHit_, meHFDepHitw_, meHFEneHit2_, meHFEneHit_, meHFEneMap_, meHFEneSum_, meHFEneSum_vs_ieta_, meHFEtaHit_, meHFL10Ene_, meHFL10EneP_, meHFNHit_, meHFPhiHit_, meHFTimHit_, meHODepHit_, meHOEneHit2_, meHOEneHit_, meHOEneMap_, meHOEneSum_, meHOEneSum_vs_ieta_, meHOEtaHit_, meHOL10Ene_, meHOL10EneP_, meHONHit_, meHOPhiHit_, meHOTimHit_, mePhiHit_, mePhiHitb_, meSubdetHit_, meTimeHit_, meTimeWHit_, phi, HcalHitRelabeller::relabel(), HcalDetId::subdet(), testNumber_, hcalRecHitTable_cff::time, and verbose_.

Referenced by analyze().

                                                           {
  int nHit = hits.size();
  int nHB = 0, nHE = 0, nHO = 0, nHF = 0, nBad1 = 0, nBad2 = 0, nBad = 0;
  std::vector<double> encontHB(140, 0.);
  std::vector<double> encontHE(140, 0.);
  std::vector<double> encontHF(140, 0.);
  std::vector<double> encontHO(140, 0.);
  double entotHB = 0, entotHE = 0, entotHF = 0, entotHO = 0;

  double HBEneMap[ieta_bins_HB][iphi_bins];
  double HEEneMap[ieta_bins_HE][iphi_bins];
  double HOEneMap[ieta_bins_HO][iphi_bins];
  double HFEneMap[ieta_bins_HF][iphi_bins];

  // Works in ieta_min_Hx is < 0
  int eta_offset_HB = -(int)ieta_min_HB;
  int eta_offset_HE = -(int)ieta_min_HE;
  int eta_offset_HO = -(int)ieta_min_HO;
  int eta_offset_HF = -(int)ieta_min_HF;

  for (int i = 0; i < ieta_bins_HB; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      HBEneMap[i][j] = 0.;
    }
  }

  for (int i = 0; i < ieta_bins_HE; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      HEEneMap[i][j] = 0.;
    }
  }

  for (int i = 0; i < ieta_bins_HO; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      HOEneMap[i][j] = 0.;
    }
  }

  for (int i = 0; i < ieta_bins_HF; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      HFEneMap[i][j] = 0.;
    }
  }

  for (int i = 0; i < nHit; i++) {
    double energy = hits[i].energy();
    double log10en = log10(energy);
    int log10i = int((log10en + 10.) * 10.);
    double time = hits[i].time();
    unsigned int id = hits[i].id();
    int det, subdet, depth, eta, phi;
    HcalDetId hid;
    if (testNumber_)
      hid = HcalHitRelabeller::relabel(id, hcons_);
    else
      hid = HcalDetId(id);
    det = hid.det();
    subdet = hid.subdet();
    depth = hid.depth();
    eta = hid.ieta();
    phi = hid.iphi();

    if (verbose_ > 1)
      edm::LogVerbatim("HcalSim") << "Hit[" << i << "] ID " << std::hex << id << std::dec << " Det " << det << " Sub "
                                  << subdet << " depth " << depth << " Eta " << eta << " Phi " << phi << " E " << energy
                                  << " time " << time;
    if (det == 4) {  // Check DetId.h
      if (subdet == static_cast<int>(HcalBarrel))
        nHB++;
      else if (subdet == static_cast<int>(HcalEndcap))
        nHE++;
      else if (subdet == static_cast<int>(HcalOuter))
        nHO++;
      else if (subdet == static_cast<int>(HcalForward))
        nHF++;
      else {
        nBad++;
        nBad2++;
      }
    } else {
      nBad++;
      nBad1++;
    }

    meDetectHit_->Fill(double(det));
    if (det == 4) {
      meSubdetHit_->Fill(double(subdet));
      meDepthHit_->Fill(double(depth));
      meEtaHit_->Fill(double(eta));
      meEtaPhiHit_->Fill(double(eta), double(phi));
      meEtaPhiHitDepth_[depth - 1]->Fill(double(eta), double(phi));

      // We will group the phi plots by HB,HO and HE,HF since these groups share
      // similar segmentation schemes
      if (subdet == static_cast<int>(HcalBarrel))
        mePhiHit_->Fill(double(phi));
      else if (subdet == static_cast<int>(HcalEndcap))
        mePhiHitb_->Fill(double(phi));
      else if (subdet == static_cast<int>(HcalOuter))
        mePhiHit_->Fill(double(phi));
      else if (subdet == static_cast<int>(HcalForward))
        mePhiHitb_->Fill(double(phi));

      // KC: HF energy is in photoelectrons rather than eV, so it will not be
      // included in total HCal energy
      if (subdet != static_cast<int>(HcalForward)) {
        meEnergyHit_->Fill(energy);

        // Since the HF energy is a different scale it does not make sense to
        // include it in the Energy Weighted Plot
        meTimeWHit_->Fill(double(time), energy);
      }
      meTimeHit_->Fill(time);

      if (subdet == static_cast<int>(HcalBarrel)) {
        meHBDepHit_->Fill(double(depth));
        meHBEtaHit_->Fill(double(eta));
        meHBPhiHit_->Fill(double(phi));
        meHBEneHit_->Fill(energy);
        meHBEneHit2_->Fill(energy);
        meHBTimHit_->Fill(time);
        meHBL10Ene_->Fill(log10en);
        if (log10i >= 0 && log10i < 140)
          encontHB[log10i] += energy;
        entotHB += energy;

        HBEneMap[eta + eta_offset_HB][phi - 1] += energy;

      } else if (subdet == static_cast<int>(HcalEndcap)) {
        meHEDepHit_->Fill(double(depth));
        meHEEtaHit_->Fill(double(eta));
        meHEPhiHit_->Fill(double(phi));

        bool isHEP17 = (phi >= 63) && (phi <= 66) && (eta > 0);
        if (hep17_) {
          if (!isHEP17) {
            meHEEneHit_->Fill(energy);
            meHEEneHit2_->Fill(energy);
          } else {
            meHEP17EneHit_->Fill(energy);
            meHEP17EneHit2_->Fill(energy);
          }
        } else {
          meHEEneHit_->Fill(energy);
          meHEEneHit2_->Fill(energy);
        }

        meHETimHit_->Fill(time);
        meHEL10Ene_->Fill(log10en);
        if (log10i >= 0 && log10i < 140)
          encontHE[log10i] += energy;
        entotHE += energy;

        HEEneMap[eta + eta_offset_HE][phi - 1] += energy;

      } else if (subdet == static_cast<int>(HcalOuter)) {
        meHODepHit_->Fill(double(depth));
        meHOEtaHit_->Fill(double(eta));
        meHOPhiHit_->Fill(double(phi));
        meHOEneHit_->Fill(energy);
        meHOEneHit2_->Fill(energy);
        meHOTimHit_->Fill(time);
        meHOL10Ene_->Fill(log10en);
        if (log10i >= 0 && log10i < 140)
          encontHO[log10i] += energy;
        entotHO += energy;

        HOEneMap[eta + eta_offset_HO][phi - 1] += energy;

      } else if (subdet == static_cast<int>(HcalForward)) {
        meHFDepHit_->Fill(double(depth));
        meHFDepHitw_->Fill(double(depth), energy);
        meHFEtaHit_->Fill(double(eta));
        meHFPhiHit_->Fill(double(phi));
        meHFEneHit_->Fill(energy);
        meHFEneHit2_->Fill(energy);
        meHFTimHit_->Fill(time);
        meHFL10Ene_->Fill(log10en);
        if (log10i >= 0 && log10i < 140)
          encontHF[log10i] += energy;
        entotHF += energy;

        HFEneMap[eta + eta_offset_HF][phi - 1] += energy;
      }
    }
  }
  if (entotHB != 0)
    for (int i = 0; i < 140; i++)
      meHBL10EneP_->Fill(-10. + (float(i) + 0.5) / 10., encontHB[i] / entotHB);
  if (entotHE != 0)
    for (int i = 0; i < 140; i++)
      meHEL10EneP_->Fill(-10. + (float(i) + 0.5) / 10., encontHE[i] / entotHE);
  if (entotHF != 0)
    for (int i = 0; i < 140; i++)
      meHFL10EneP_->Fill(-10. + (float(i) + 0.5) / 10., encontHF[i] / entotHF);
  if (entotHO != 0)
    for (int i = 0; i < 140; i++)
      meHOL10EneP_->Fill(-10. + (float(i) + 0.5) / 10., encontHO[i] / entotHO);

  meAllNHit_->Fill(double(nHit));
  meBadDetHit_->Fill(double(nBad1));
  meBadSubHit_->Fill(double(nBad2));
  meBadIdHit_->Fill(double(nBad));
  meHBNHit_->Fill(double(nHB));
  meHENHit_->Fill(double(nHE));
  meHONHit_->Fill(double(nHO));
  meHFNHit_->Fill(double(nHF));

  for (int i = 0; i < ieta_bins_HB; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      if (HBEneMap[i][j] != 0) {
        meHBEneSum_->Fill(HBEneMap[i][j]);
        meHBEneSum_vs_ieta_->Fill((i - eta_offset_HB), HBEneMap[i][j]);
        meHBEneMap_->Fill((i - eta_offset_HB), j + 1, HBEneMap[i][j]);
      }
    }
  }

  for (int i = 0; i < ieta_bins_HE; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      if (HEEneMap[i][j] != 0) {
        meHEEneSum_->Fill(HEEneMap[i][j]);
        meHEEneSum_vs_ieta_->Fill((i - eta_offset_HE), HEEneMap[i][j]);
        meHEEneMap_->Fill((i - eta_offset_HE), j + 1, HEEneMap[i][j]);
      }
    }
  }

  for (int i = 0; i < ieta_bins_HO; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      if (HOEneMap[i][j] != 0) {
        meHOEneSum_->Fill(HOEneMap[i][j]);
        meHOEneSum_vs_ieta_->Fill((i - eta_offset_HO), HOEneMap[i][j]);
        meHOEneMap_->Fill((i - eta_offset_HO), j + 1, HOEneMap[i][j]);
      }
    }
  }

  for (int i = 0; i < ieta_bins_HF; i++) {
    for (int j = 0; j < iphi_bins; j++) {
      if (HFEneMap[i][j] != 0) {
        meHFEneSum_->Fill(HFEneMap[i][j]);
        meHFEneSum_vs_ieta_->Fill((i - eta_offset_HF), HFEneMap[i][j]);
        meHFEneMap_->Fill((i - eta_offset_HF), j + 1, HFEneMap[i][j]);
      }
    }
  }

  if (verbose_ > 0)
    edm::LogVerbatim("HcalSim") << "HcalSimHitCheck::analyzeHits: HB " << nHB << " HE " << nHE << " HO " << nHO
                                << " HF " << nHF << " Bad " << nBad << " All " << nHit;
}

beginRun()

void HcalSimHitCheck::beginRun ( edm::Run const &  , edm::EventSetup const &  es  )

overrideprotected

Definition at line 116 of file HcalSimHitCheck.cc.

References checkHit_, hcalRecHitTable_cff::depth, compareTotals::fs, edm::EventSetup::getData(), HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNPhi(), hcons_, ieta_bins_HB, ieta_bins_HE, ieta_bins_HF, ieta_bins_HO, ieta_max_HB, ieta_max_HE, ieta_max_HF, ieta_max_HO, ieta_min_HB, ieta_min_HE, ieta_min_HF, ieta_min_HO, hcalTTPDigis_cfi::iEtaMax, createfilelist::int, iphi_bins, iphi_max, iphi_min, SiStripPI::max, maxDepth_, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meAllNHit_, meBadDetHit_, meBadIdHit_, meBadSubHit_, meDepthHit_, meDetectHit_, meEnergyHit_, meEtaHit_, meEtaPhiHit_, meEtaPhiHitDepth_, meHBDepHit_, meHBEneHit2_, meHBEneHit_, meHBEneMap_, meHBEneSum_, meHBEneSum_vs_ieta_, meHBEtaHit_, meHBL10Ene_, meHBL10EneP_, meHBNHit_, meHBPhiHit_, meHBTimHit_, meHEDepHit_, meHEEneHit2_, meHEEneHit_, meHEEneMap_, meHEEneSum_, meHEEneSum_vs_ieta_, meHEEtaHit_, meHEL10Ene_, meHEL10EneP_, meHENHit_, meHEP17EneHit2_, meHEP17EneHit_, meHEPhiHit_, meHETimHit_, meHFDepHit_, meHFDepHitw_, meHFEneHit2_, meHFEneHit_, meHFEneMap_, meHFEneSum_, meHFEneSum_vs_ieta_, meHFEtaHit_, meHFL10Ene_, meHFL10EneP_, meHFNHit_, meHFPhiHit_, meHFTimHit_, meHODepHit_, meHOEneHit2_, meHOEneHit_, meHOEneMap_, meHOEneSum_, meHOEneSum_vs_ieta_, meHOEtaHit_, meHOL10Ene_, meHOL10EneP_, meHONHit_, meHOPhiHit_, meHOTimHit_, mePhiHit_, mePhiHitb_, meSubdetHit_, meTimeHit_, meTimeWHit_, and tok_HRNDC_.

                                                                      {
  hcons_ = &es.getData(tok_HRNDC_);
  maxDepthHB_ = hcons_->getMaxDepth(0);
  maxDepthHE_ = hcons_->getMaxDepth(1);
  maxDepthHF_ = hcons_->getMaxDepth(2);
  maxDepthHO_ = hcons_->getMaxDepth(3);
  maxDepth_ = (maxDepthHB_ > maxDepthHE_ ? maxDepthHB_ : maxDepthHE_);
  maxDepth_ = (maxDepth_ > maxDepthHF_ ? maxDepth_ : maxDepthHF_);
  maxDepth_ = (maxDepth_ > maxDepthHO_ ? maxDepth_ : maxDepthHO_);

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

  int iEtaHBMax = hcons_->getEtaRange(0).second;
  int iEtaHEMax = std::max(hcons_->getEtaRange(1).second, 1);
  int iEtaHFMax = hcons_->getEtaRange(2).second;
  int iEtaHOMax = hcons_->getEtaRange(3).second;

  // Retain classic behavior, all plots have same ieta range.
  // Comment out code to allow each subdetector to have its on range

  int iEtaMax = (iEtaHBMax > iEtaHEMax ? iEtaHBMax : iEtaHEMax);
  iEtaMax = (iEtaMax > iEtaHFMax ? iEtaMax : iEtaHFMax);
  iEtaMax = (iEtaMax > iEtaHOMax ? iEtaMax : iEtaHOMax);

  iEtaHBMax = iEtaMax;
  iEtaHEMax = iEtaMax;
  iEtaHFMax = iEtaMax;
  iEtaHOMax = iEtaMax;

  // Give an empty bin around the subdet ieta range to make it clear that all
  // ieta rings have been included
  ieta_min_HB = -iEtaHBMax - 1.5;
  ieta_max_HB = iEtaHBMax + 1.5;
  ieta_bins_HB = (int)(ieta_max_HB - ieta_min_HB);

  ieta_min_HE = -iEtaHEMax - 1.5;
  ieta_max_HE = iEtaHEMax + 1.5;
  ieta_bins_HE = (int)(ieta_max_HE - ieta_min_HE);

  ieta_min_HF = -iEtaHFMax - 1.5;
  ieta_max_HF = iEtaHFMax + 1.5;
  ieta_bins_HF = (int)(ieta_max_HF - ieta_min_HF);

  ieta_min_HO = -iEtaHOMax - 1.5;
  ieta_max_HO = iEtaHOMax + 1.5;
  ieta_bins_HO = (int)(ieta_max_HO - ieta_min_HO);

  Char_t hname[100];
  Char_t htitle[100];

  edm::Service<TFileService> fs;

  // Histograms for Hits
  if (checkHit_) {
    meAllNHit_ = fs->make<TH1D>("Hit01", "Number of Hits in HCal", 20000, 0., 20000.);
    meBadDetHit_ = fs->make<TH1D>("Hit02", "Hits with wrong Det", 100, 0., 100.);
    meBadSubHit_ = fs->make<TH1D>("Hit03", "Hits with wrong Subdet", 100, 0., 100.);
    meBadIdHit_ = fs->make<TH1D>("Hit04", "Hits with wrong ID", 100, 0., 100.);
    meHBNHit_ = fs->make<TH1D>("Hit05", "Number of Hits in HB", 20000, 0., 20000.);
    meHENHit_ = fs->make<TH1D>("Hit06", "Number of Hits in HE", 10000, 0., 10000.);
    meHONHit_ = fs->make<TH1D>("Hit07", "Number of Hits in HO", 10000, 0., 10000.);
    meHFNHit_ = fs->make<TH1D>("Hit08", "Number of Hits in HF", 10000, 0., 10000.);
    meHBNHit_->Sumw2();
    meHENHit_->Sumw2();
    meHONHit_->Sumw2();
    meHFNHit_->Sumw2();
    meDetectHit_ = fs->make<TH1D>("Hit09", "Detector ID", 50, 0., 50.);
    meSubdetHit_ = fs->make<TH1D>("Hit10", "Subdetectors in HCal", 50, 0., 50.);
    meDepthHit_ = fs->make<TH1D>("Hit11", "Depths in HCal", 20, 0., 20.);
    meEtaHit_ = fs->make<TH1D>("Hit12", "Eta in HCal", ieta_bins_HF, ieta_min_HF, ieta_max_HF);
    meEtaPhiHit_ = fs->make<TH2D>(
        "Hit12b", "Eta-phi in HCal", ieta_bins_HF, ieta_min_HF, ieta_max_HF, iphi_bins, iphi_min, iphi_max);
    for (int depth = 1; depth <= maxDepth_; depth++) {
      sprintf(hname, "Hit12bd%d", depth);
      sprintf(htitle, "Eta-phi in HCal d%d", depth);
      meEtaPhiHitDepth_.emplace_back(
          fs->make<TH2D>(hname, htitle, ieta_bins_HF, ieta_min_HF, ieta_max_HF, iphi_bins, iphi_min, iphi_max));
    }
    // KC: There are different phi segmentation schemes, this plot uses wider
    // bins to represent the most sparse segmentation
    mePhiHit_ = fs->make<TH1D>("Hit13", "Phi in HCal (HB,HO)", iphi_bins, iphi_min, iphi_max);
    mePhiHitb_ = fs->make<TH1D>("Hit13b", "Phi in HCal (HE,HF)", iphi_bins, iphi_min, iphi_max);
    meEnergyHit_ = fs->make<TH1D>("Hit14", "Energy in HCal", 2000, 0., 20.);
    meTimeHit_ = fs->make<TH1D>("Hit15", "Time in HCal", 528, 0., 528.);
    meTimeWHit_ = fs->make<TH1D>("Hit16", "Time in HCal (E wtd)", 528, 0., 528.);
    meHBDepHit_ = fs->make<TH1D>("Hit17", "Depths in HB", 20, 0., 20.);
    meHEDepHit_ = fs->make<TH1D>("Hit18", "Depths in HE", 20, 0., 20.);
    meHODepHit_ = fs->make<TH1D>("Hit19", "Depths in HO", 20, 0., 20.);
    meHFDepHit_ = fs->make<TH1D>("Hit20", "Depths in HF", 20, 0., 20.);
    meHBDepHit_->Sumw2();
    meHEDepHit_->Sumw2();
    meHODepHit_->Sumw2();
    meHFDepHit_->Sumw2();
    meHFDepHitw_ = fs->make<TH1D>("Hit20b", "Depths in HF (p.e. weighted)", 20, 0., 20.);
    meHBEtaHit_ = fs->make<TH1D>("Hit21", "Eta in HB", ieta_bins_HB, ieta_min_HB, ieta_max_HB);
    meHEEtaHit_ = fs->make<TH1D>("Hit22", "Eta in HE", ieta_bins_HE, ieta_min_HE, ieta_max_HE);
    meHOEtaHit_ = fs->make<TH1D>("Hit23", "Eta in HO", ieta_bins_HO, ieta_min_HO, ieta_max_HO);
    meHFEtaHit_ = fs->make<TH1D>("Hit24", "Eta in HF", ieta_bins_HF, ieta_min_HF, ieta_max_HF);
    meHBEtaHit_->Sumw2();
    meHEEtaHit_->Sumw2();
    meHOEtaHit_->Sumw2();
    meHFEtaHit_->Sumw2();
    meHBPhiHit_ = fs->make<TH1D>("Hit25", "Phi in HB", iphi_bins, iphi_min, iphi_max);
    meHEPhiHit_ = fs->make<TH1D>("Hit26", "Phi in HE", iphi_bins, iphi_min, iphi_max);
    meHOPhiHit_ = fs->make<TH1D>("Hit27", "Phi in HO", iphi_bins, iphi_min, iphi_max);
    meHFPhiHit_ = fs->make<TH1D>("Hit28", "Phi in HF", iphi_bins, iphi_min, iphi_max);
    meHBPhiHit_->Sumw2();
    meHEPhiHit_->Sumw2();
    meHOPhiHit_->Sumw2();
    meHFPhiHit_->Sumw2();
    meHBEneHit_ = fs->make<TH1D>("Hit29", "Energy in HB", 2000, 0., 20.);
    meHEEneHit_ = fs->make<TH1D>("Hit30", "Energy in HE", 500, 0., 5.);
    meHEP17EneHit_ = fs->make<TH1D>("Hit30b", "Energy in HEP17", 500, 0., 5.);
    meHOEneHit_ = fs->make<TH1D>("Hit31", "Energy in HO", 500, 0., 5.);
    meHFEneHit_ = fs->make<TH1D>("Hit32", "Energy in HF", 1001, -0.5, 1000.5);
    meHBEneHit_->Sumw2();
    meHEEneHit_->Sumw2();
    meHOEneHit_->Sumw2();
    meHFEneHit_->Sumw2();

    // HxEneMap, HxEneSum, HxEneSum_vs_ieta plot the sum of the simhits energy
    // within a single ieta-iphi tower.

    meHBEneMap_ =
        fs->make<TH2D>("HBEneMap", "HBEneMap", ieta_bins_HB, ieta_min_HB, ieta_max_HB, iphi_bins, iphi_min, iphi_max);
    meHEEneMap_ =
        fs->make<TH2D>("HEEneMap", "HEEneMap", ieta_bins_HE, ieta_min_HE, ieta_max_HE, iphi_bins, iphi_min, iphi_max);
    meHOEneMap_ =
        fs->make<TH2D>("HOEneMap", "HOEneMap", ieta_bins_HO, ieta_min_HO, ieta_max_HO, iphi_bins, iphi_min, iphi_max);
    meHFEneMap_ =
        fs->make<TH2D>("HFEneMap", "HFEneMap", ieta_bins_HF, ieta_min_HF, ieta_max_HF, iphi_bins, iphi_min, iphi_max);

    meHBEneSum_ = fs->make<TH1D>("HBEneSum", "HBEneSum", 2000, 0., 20.);
    meHEEneSum_ = fs->make<TH1D>("HEEneSum", "HEEneSum", 500, 0., 5.);
    meHOEneSum_ = fs->make<TH1D>("HOEneSum", "HOEneSum", 500, 0., 5.);
    meHFEneSum_ = fs->make<TH1D>("HFEneSum", "HFEneSum", 1001, -0.5, 1000.5);
    meHBEneSum_->Sumw2();
    meHEEneSum_->Sumw2();
    meHOEneSum_->Sumw2();
    meHFEneSum_->Sumw2();

    meHBEneSum_vs_ieta_ = fs->make<TProfile>(
        "HBEneSum_vs_ieta", "HBEneSum_vs_ieta", ieta_bins_HB, ieta_min_HB, ieta_max_HB, -10.5, 2000.5);
    meHEEneSum_vs_ieta_ = fs->make<TProfile>(
        "HEEneSum_vs_ieta", "HEEneSum_vs_ieta", ieta_bins_HE, ieta_min_HE, ieta_max_HE, -10.5, 2000.5);
    meHOEneSum_vs_ieta_ = fs->make<TProfile>(
        "HOEneSum_vs_ieta", "HOEneSum_vs_ieta", ieta_bins_HO, ieta_min_HO, ieta_max_HO, -10.5, 2000.5);
    meHFEneSum_vs_ieta_ = fs->make<TProfile>(
        "HFEneSum_vs_ieta", "HFEneSum_vs_ieta", ieta_bins_HF, ieta_min_HF, ieta_max_HF, -10.5, 2000.5);

    meHBTimHit_ = fs->make<TH1D>("Hit33", "Time in HB", 528, 0., 528.);
    meHETimHit_ = fs->make<TH1D>("Hit34", "Time in HE", 528, 0., 528.);
    meHOTimHit_ = fs->make<TH1D>("Hit35", "Time in HO", 528, 0., 528.);
    meHFTimHit_ = fs->make<TH1D>("Hit36", "Time in HF", 528, 0., 528.);
    meHBTimHit_->Sumw2();
    meHETimHit_->Sumw2();
    meHOTimHit_->Sumw2();
    meHFTimHit_->Sumw2();
    // These are the zoomed in energy ranges
    meHBEneHit2_ = fs->make<TH1D>("Hit37", "Energy in HB 2", 100, 0., 0.0001);
    meHEEneHit2_ = fs->make<TH1D>("Hit38", "Energy in HE 2", 100, 0., 0.0001);
    meHEP17EneHit2_ = fs->make<TH1D>("Hit38b", "Energy in HEP17 2", 100, 0., 0.0001);
    meHOEneHit2_ = fs->make<TH1D>("Hit39", "Energy in HO 2", 100, 0., 0.0001);
    meHFEneHit2_ = fs->make<TH1D>("Hit40", "Energy in HF 2", 100, 0.5, 100.5);
    meHBL10Ene_ = fs->make<TH1D>("Hit41", "Log10Energy in HB", 140, -10., 4.);
    meHEL10Ene_ = fs->make<TH1D>("Hit42", "Log10Energy in HE", 140, -10., 4.);
    meHFL10Ene_ = fs->make<TH1D>("Hit43", "Log10Energy in HF", 50, -1., 4.);
    meHOL10Ene_ = fs->make<TH1D>("Hit44", "Log10Energy in HO", 140, -10., 4.);
    meHBL10EneP_ = fs->make<TProfile>("Hit45", "Log10Energy in HB vs Hit contribution", 140, -10., 4., 0., 1.);
    meHEL10EneP_ = fs->make<TProfile>("Hit46", "Log10Energy in HE vs Hit contribution", 140, -10., 4., 0., 1.);
    meHFL10EneP_ = fs->make<TProfile>("Hit47", "Log10Energy in HF vs Hit contribution", 140, -10., 4., 0., 1.);
    meHOL10EneP_ = fs->make<TProfile>("Hit48", "Log10Energy in HO vs Hit contribution", 140, -10., 4., 0., 1.);
  }
}

endRun()

void HcalSimHitCheck::endRun ( edm::Run const &  , edm::EventSetup const &    )

inlineoverrideprotected

Definition at line 43 of file HcalSimHitCheck.cc.

{}

fillDescriptions()

void HcalSimHitCheck::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 105 of file HcalSimHitCheck.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
  edm::ParameterSetDescription desc;
  desc.add<std::string>("moduleLabel", "g4SimHits");
  desc.add<std::string>("HitCollection", "HcalHits");
  desc.add<std::string>("outputFile", "hcHit.root");
  desc.add<int>("Verbose", 0);
  desc.add<bool>("TestNumber", true);
  desc.add<bool>("hep17", false);
  descriptions.add("hcalSimHitCheck", desc);
}

Member Data Documentation

checkHit_

const bool HcalSimHitCheck::checkHit_

private

Definition at line 66 of file HcalSimHitCheck.cc.

Referenced by analyze(), beginRun(), and HcalSimHitCheck().

g4Label

const std::string HcalSimHitCheck::g4Label

private

Definition at line 64 of file HcalSimHitCheck.cc.

Referenced by HcalSimHitCheck().

hcalHits

const std::string HcalSimHitCheck::hcalHits

private

Definition at line 64 of file HcalSimHitCheck.cc.

Referenced by HcalSimHitCheck().

hcons_

const HcalDDDRecConstants* HcalSimHitCheck::hcons_

private

Definition at line 48 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

hep17_

const bool HcalSimHitCheck::hep17_

private

Definition at line 66 of file HcalSimHitCheck.cc.

Referenced by analyzeHits().

ieta_bins_HB

int HcalSimHitCheck::ieta_bins_HB

private

Definition at line 55 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_bins_HE

int HcalSimHitCheck::ieta_bins_HE

private

Definition at line 57 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_bins_HF

int HcalSimHitCheck::ieta_bins_HF

private

Definition at line 61 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_bins_HO

int HcalSimHitCheck::ieta_bins_HO

private

Definition at line 59 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_max_HB

float HcalSimHitCheck::ieta_max_HB

private

Definition at line 56 of file HcalSimHitCheck.cc.

Referenced by beginRun().

ieta_max_HE

float HcalSimHitCheck::ieta_max_HE

private

Definition at line 58 of file HcalSimHitCheck.cc.

Referenced by beginRun().

ieta_max_HF

float HcalSimHitCheck::ieta_max_HF

private

Definition at line 62 of file HcalSimHitCheck.cc.

Referenced by beginRun().

ieta_max_HO

float HcalSimHitCheck::ieta_max_HO

private

Definition at line 60 of file HcalSimHitCheck.cc.

Referenced by beginRun().

ieta_min_HB

float HcalSimHitCheck::ieta_min_HB

private

Definition at line 56 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_min_HE

float HcalSimHitCheck::ieta_min_HE

private

Definition at line 58 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_min_HF

float HcalSimHitCheck::ieta_min_HF

private

Definition at line 62 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

ieta_min_HO

float HcalSimHitCheck::ieta_min_HO

private

Definition at line 60 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

iphi_bins

int HcalSimHitCheck::iphi_bins

private

Definition at line 53 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

iphi_max

float HcalSimHitCheck::iphi_max

private

Definition at line 54 of file HcalSimHitCheck.cc.

Referenced by beginRun().

iphi_min

float HcalSimHitCheck::iphi_min

private

Definition at line 54 of file HcalSimHitCheck.cc.

Referenced by beginRun().

maxDepth_

int HcalSimHitCheck::maxDepth_

private

Definition at line 51 of file HcalSimHitCheck.cc.

Referenced by beginRun().

maxDepthHB_

int HcalSimHitCheck::maxDepthHB_

private

Definition at line 49 of file HcalSimHitCheck.cc.

Referenced by beginRun().

maxDepthHE_

int HcalSimHitCheck::maxDepthHE_

private

Definition at line 49 of file HcalSimHitCheck.cc.

Referenced by beginRun().

maxDepthHF_

int HcalSimHitCheck::maxDepthHF_

private

Definition at line 50 of file HcalSimHitCheck.cc.

Referenced by beginRun().

maxDepthHO_

int HcalSimHitCheck::maxDepthHO_

private

Definition at line 50 of file HcalSimHitCheck.cc.

Referenced by beginRun().

meAllNHit_

TH1D* HcalSimHitCheck::meAllNHit_

private

Definition at line 71 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meBadDetHit_

TH1D * HcalSimHitCheck::meBadDetHit_

private

Definition at line 71 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meBadIdHit_

TH1D * HcalSimHitCheck::meBadIdHit_

private

Definition at line 71 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meBadSubHit_

TH1D * HcalSimHitCheck::meBadSubHit_

private

Definition at line 71 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meDepthHit_

TH1D * HcalSimHitCheck::meDepthHit_

private

Definition at line 73 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meDetectHit_

TH1D* HcalSimHitCheck::meDetectHit_

private

Definition at line 73 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meEnergyHit_

TH1D * HcalSimHitCheck::meEnergyHit_

private

Definition at line 76 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meEtaHit_

TH1D * HcalSimHitCheck::meEtaHit_

private

Definition at line 73 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meEtaPhiHit_

TH2D* HcalSimHitCheck::meEtaPhiHit_

private

Definition at line 74 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meEtaPhiHitDepth_

std::vector<TH2D *> HcalSimHitCheck::meEtaPhiHitDepth_

private

Definition at line 75 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBDepHit_

TH1D* HcalSimHitCheck::meHBDepHit_

private

Definition at line 77 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBEneHit2_

TH1D* HcalSimHitCheck::meHBEneHit2_

private

Definition at line 85 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBEneHit_

TH1D* HcalSimHitCheck::meHBEneHit_

private

Definition at line 80 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBEneMap_

TH2D* HcalSimHitCheck::meHBEneMap_

private

Definition at line 81 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBEneSum_

TH1D* HcalSimHitCheck::meHBEneSum_

private

Definition at line 82 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBEneSum_vs_ieta_

TProfile* HcalSimHitCheck::meHBEneSum_vs_ieta_

private

Definition at line 83 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBEtaHit_

TH1D* HcalSimHitCheck::meHBEtaHit_

private

Definition at line 78 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBL10Ene_

TH1D* HcalSimHitCheck::meHBL10Ene_

private

Definition at line 86 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBL10EneP_

TProfile* HcalSimHitCheck::meHBL10EneP_

private

Definition at line 87 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBNHit_

TH1D* HcalSimHitCheck::meHBNHit_

private

Definition at line 72 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBPhiHit_

TH1D* HcalSimHitCheck::meHBPhiHit_

private

Definition at line 79 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHBTimHit_

TH1D* HcalSimHitCheck::meHBTimHit_

private

Definition at line 84 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEDepHit_

TH1D * HcalSimHitCheck::meHEDepHit_

private

Definition at line 77 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEEneHit2_

TH1D * HcalSimHitCheck::meHEEneHit2_

private

Definition at line 85 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEEneHit_

TH1D * HcalSimHitCheck::meHEEneHit_

private

Definition at line 80 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEEneMap_

TH2D * HcalSimHitCheck::meHEEneMap_

private

Definition at line 81 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEEneSum_

TH1D * HcalSimHitCheck::meHEEneSum_

private

Definition at line 82 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEEneSum_vs_ieta_

TProfile * HcalSimHitCheck::meHEEneSum_vs_ieta_

private

Definition at line 83 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEEtaHit_

TH1D * HcalSimHitCheck::meHEEtaHit_

private

Definition at line 78 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEL10Ene_

TH1D * HcalSimHitCheck::meHEL10Ene_

private

Definition at line 86 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEL10EneP_

TProfile * HcalSimHitCheck::meHEL10EneP_

private

Definition at line 87 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHENHit_

TH1D * HcalSimHitCheck::meHENHit_

private

Definition at line 72 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEP17EneHit2_

TH1D * HcalSimHitCheck::meHEP17EneHit2_

private

Definition at line 88 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEP17EneHit_

TH1D* HcalSimHitCheck::meHEP17EneHit_

private

Definition at line 88 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHEPhiHit_

TH1D * HcalSimHitCheck::meHEPhiHit_

private

Definition at line 79 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHETimHit_

TH1D * HcalSimHitCheck::meHETimHit_

private

Definition at line 84 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFDepHit_

TH1D * HcalSimHitCheck::meHFDepHit_

private

Definition at line 77 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFDepHitw_

TH1D * HcalSimHitCheck::meHFDepHitw_

private

Definition at line 77 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFEneHit2_

TH1D * HcalSimHitCheck::meHFEneHit2_

private

Definition at line 85 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFEneHit_

TH1D * HcalSimHitCheck::meHFEneHit_

private

Definition at line 80 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFEneMap_

TH2D * HcalSimHitCheck::meHFEneMap_

private

Definition at line 81 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFEneSum_

TH1D * HcalSimHitCheck::meHFEneSum_

private

Definition at line 82 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFEneSum_vs_ieta_

TProfile * HcalSimHitCheck::meHFEneSum_vs_ieta_

private

Definition at line 83 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFEtaHit_

TH1D * HcalSimHitCheck::meHFEtaHit_

private

Definition at line 78 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFL10Ene_

TH1D * HcalSimHitCheck::meHFL10Ene_

private

Definition at line 86 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFL10EneP_

TProfile * HcalSimHitCheck::meHFL10EneP_

private

Definition at line 87 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFNHit_

TH1D * HcalSimHitCheck::meHFNHit_

private

Definition at line 72 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFPhiHit_

TH1D * HcalSimHitCheck::meHFPhiHit_

private

Definition at line 79 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHFTimHit_

TH1D * HcalSimHitCheck::meHFTimHit_

private

Definition at line 84 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHODepHit_

TH1D * HcalSimHitCheck::meHODepHit_

private

Definition at line 77 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOEneHit2_

TH1D * HcalSimHitCheck::meHOEneHit2_

private

Definition at line 85 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOEneHit_

TH1D * HcalSimHitCheck::meHOEneHit_

private

Definition at line 80 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOEneMap_

TH2D * HcalSimHitCheck::meHOEneMap_

private

Definition at line 81 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOEneSum_

TH1D * HcalSimHitCheck::meHOEneSum_

private

Definition at line 82 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOEneSum_vs_ieta_

TProfile * HcalSimHitCheck::meHOEneSum_vs_ieta_

private

Definition at line 83 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOEtaHit_

TH1D * HcalSimHitCheck::meHOEtaHit_

private

Definition at line 78 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOL10Ene_

TH1D * HcalSimHitCheck::meHOL10Ene_

private

Definition at line 86 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOL10EneP_

TProfile * HcalSimHitCheck::meHOL10EneP_

private

Definition at line 87 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHONHit_

TH1D * HcalSimHitCheck::meHONHit_

private

Definition at line 72 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOPhiHit_

TH1D * HcalSimHitCheck::meHOPhiHit_

private

Definition at line 79 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meHOTimHit_

TH1D * HcalSimHitCheck::meHOTimHit_

private

Definition at line 84 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

mePhiHit_

TH1D* HcalSimHitCheck::mePhiHit_

private

Definition at line 76 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

mePhiHitb_

TH1D * HcalSimHitCheck::mePhiHitb_

private

Definition at line 76 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meSubdetHit_

TH1D * HcalSimHitCheck::meSubdetHit_

private

Definition at line 73 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meTimeHit_

TH1D * HcalSimHitCheck::meTimeHit_

private

Definition at line 76 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

meTimeWHit_

TH1D * HcalSimHitCheck::meTimeWHit_

private

Definition at line 76 of file HcalSimHitCheck.cc.

Referenced by analyzeHits(), and beginRun().

outFile_

const std::string HcalSimHitCheck::outFile_

private

Definition at line 64 of file HcalSimHitCheck.cc.

Referenced by HcalSimHitCheck().

testNumber_

const bool HcalSimHitCheck::testNumber_

private

Definition at line 66 of file HcalSimHitCheck.cc.

Referenced by analyzeHits().

tok_hits_

const edm::EDGetTokenT<edm::PCaloHitContainer> HcalSimHitCheck::tok_hits_

private

Definition at line 68 of file HcalSimHitCheck.cc.

Referenced by analyze().

tok_HRNDC_

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

private

Definition at line 69 of file HcalSimHitCheck.cc.

Referenced by beginRun().

verbose_

const int HcalSimHitCheck::verbose_

private

Definition at line 65 of file HcalSimHitCheck.cc.

Referenced by analyze(), and analyzeHits().