HOSimHitStudy Class Reference

Inheritance diagram for HOSimHitStudy:

Public Member Functions
	HOSimHitStudy (const edm::ParameterSet &ps)
	~HOSimHitStudy () override=default
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 (const edm::Event &e, const edm::EventSetup &c) override
void	analyzeHits ()
void	beginJob () override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
void	endJob () 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
std::vector< PCaloHit >	ecalHits
TH1F *	edep_ [3]
TH1F *	edepT_ [3]
TH1F *	edepTW_ [3]
TH1F *	edepTWT_ [3]
TH1F *	edepZon_ [3]
TH1F *	edepZonT_ [3]
TH1F *	eEB_
TH1F *	eEBHB_
TH1F *	eEBHBHO_
TH1F *	eEBHBHOT_
TH1F *	eEBHBT_
TH1F *	eEBT_
TProfile *	eHO17_
TProfile *	eHO17T_
TProfile *	eHO18_
TProfile *	eHO18T_
TProfile *	eHO1_
TProfile *	eHO1T_
TProfile2D *	eHO2_
TProfile2D *	eHO2T_
TH1F *	eHOE17_ [15]
TH1F *	eHOE17T_ [15]
TH1F *	eHOE18_ [15]
TH1F *	eHOE18T_ [15]
TH1F *	eHOE_ [15]
TH1F *	eHOET_ [15]
TH1F *	eHOEta17_ [15]
TH1F *	eHOEta17T_ [15]
TH1F *	eHOEta18_ [15]
TH1F *	eHOEta18T_ [15]
TH1F *	eHOEta_ [15]
TH1F *	eHOEtaT_ [15]
double	eInc
TH1F *	eneInc_
double	etaInc
TH1F *	etaInc_
const std::string	g4Label_
std::vector< PCaloHit >	hcalHits
TH1F *	hit_ [3]
const std::vector< std::string >	hitLab_
TH1F *	hitTow_ [3]
const double	maxEnergy_
TProfile *	nHO1_
TProfile *	nHO1T_
TProfile2D *	nHO2_
TProfile2D *	nHO2T_
TH1F *	nHOE1_ [15]
TH1F *	nHOE1T_ [15]
TProfile *	nHOE2_ [15]
TProfile *	nHOE2T_ [15]
TH1F *	nHOEta1_ [15]
TH1F *	nHOEta1T_ [15]
TProfile *	nHOEta2_ [15]
TProfile *	nHOEta2T_ [15]
double	phiInc
TH1F *	phiInc_
const bool	print_
const double	scaleEB_
const double	scaleHB_
const double	scaleHO_
const bool	scheme_
const double	tcut_
TH1F *	time_ [3]
const edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_
const std::vector< edm::EDGetTokenT< edm::PCaloHitContainer > >	toks_calo_

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 41 of file HOSimHitStudy.cc.

Constructor & Destructor Documentation

HOSimHitStudy()

HOSimHitStudy::HOSimHitStudy

(

const edm::ParameterSet &

ps

)

Definition at line 81 of file HOSimHitStudy.cc.

References g4Label_, hitLab_, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, and edm::vector_transform().

    : g4Label_(ps.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits")),
      hitLab_(ps.getParameter<std::vector<std::string>>("HitCollection")),
      maxEnergy_(ps.getUntrackedParameter<double>("MaxEnergy", 200.0)),
      scaleEB_(ps.getUntrackedParameter<double>("ScaleEB", 1.0)),
      scaleHB_(ps.getUntrackedParameter<double>("ScaleHB", 100.0)),
      scaleHO_(ps.getUntrackedParameter<double>("ScaleHO", 2.0)),
      tcut_(ps.getUntrackedParameter<double>("TimeCut", 100.0)),
      scheme_(ps.getUntrackedParameter<bool>("TestNumbering", false)),
      print_(ps.getUntrackedParameter<bool>("PrintExcessEnergy", true)),
      tok_evt_(consumes<edm::HepMCProduct>(
          edm::InputTag(ps.getUntrackedParameter<std::string>("SourceLabel", "VtxSmeared")))),
      toks_calo_{edm::vector_transform(hitLab_, [this](const std::string &name) {
        return consumes<edm::PCaloHitContainer>(edm::InputTag{g4Label_, name});
      })} {
  usesResource(TFileService::kSharedResource);

  edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Module Label: " << g4Label_ << "   Hits: " << hitLab_[0] << ", "
                               << hitLab_[1] << "   MaxEnergy: " << maxEnergy_ << "   Scale factor for EB " << scaleEB_
                               << ", for HB " << scaleHB_ << " and for HO " << scaleHO_ << " time Cut " << tcut_;
}

~HOSimHitStudy()

HOSimHitStudy::~HOSimHitStudy ( )

overridedefault

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 364 of file HOSimHitStudy.cc.

References analyzeHits(), MillePedeFileConverter_cfg::e, ecalHits, eInc, etaInc, edm::HepMCProduct::GetEvent(), hcalHits, hitLab_, mps_fire::i, edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, phiInc, tok_evt_, and toks_calo_.

                                                                    {
  edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Run = " << e.id().run() << " Event = " << e.id().event();

  const edm::Handle<edm::HepMCProduct> &EvtHandle = e.getHandle(tok_evt_);
  const HepMC::GenEvent *myGenEvent = EvtHandle->GetEvent();

  eInc = etaInc = phiInc = 0;
  HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
  if (p != myGenEvent->particles_end()) {
    eInc = (*p)->momentum().e();
    etaInc = (*p)->momentum().eta();
    phiInc = (*p)->momentum().phi();
  }

  edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Energy = " << eInc << " Eta = " << etaInc
                               << " Phi = " << phiInc / CLHEP::deg;

  for (int i = 0; i < 2; i++) {
    bool getHits = false;
    if (i == 0)
      ecalHits.clear();
    else
      hcalHits.clear();
    const edm::Handle<edm::PCaloHitContainer> &hitsCalo = e.getHandle(toks_calo_[i]);
    if (hitsCalo.isValid())
      getHits = true;
    edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Input flag " << hitLab_[i] << " getHits flag " << getHits;

    if (getHits) {
      unsigned int isiz;
      if (i == 0) {
        ecalHits.insert(ecalHits.end(), hitsCalo->begin(), hitsCalo->end());
        isiz = ecalHits.size();
      } else {
        hcalHits.insert(hcalHits.end(), hitsCalo->begin(), hitsCalo->end());
        isiz = hcalHits.size();
      }
      edm::LogVerbatim("HitStudy") << "HOSimHitStudy:: Hit buffer for " << hitLab_[i] << " has " << isiz << " hits";
    }
  }
  analyzeHits();
}

analyzeHits()

void HOSimHitStudy::analyzeHits ( )

protected

Definition at line 407 of file HOSimHitStudy.cc.

References funct::abs(), TauDecayModes::dec, hcalRecHitTable_cff::depth, ecalHits, edep_, edepT_, edepTW_, edepTWT_, edepZon_, edepZonT_, eEB_, eEBHB_, eEBHBHO_, eEBHBHOT_, eEBHBT_, eEBT_, eHO17_, eHO17T_, eHO18_, eHO18T_, eHO1_, eHO1T_, eHO2_, eHO2T_, eHOE17_, eHOE17T_, eHOE18_, eHOE18T_, eHOE_, eHOET_, eHOEta17_, eHOEta18_, eHOEta_, eInc, eneInc_, PVValHelper::eta, etaInc, etaInc_, HcalBarrel, hcalHits, HcalOuter, hit_, hitTow_, mps_fire::i, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, dqmiolumiharvest::j, dqmdumpme::k, maxEnergy_, nHO1_, nHO1T_, nHO2_, nHO2T_, nHOE1_, nHOE1T_, nHOE2_, nHOE2T_, nHOEta1_, nHOEta1T_, nHOEta2_, nHOEta2T_, convertSQLiteXML::ok, phi, phiInc, phiInc_, print_, scaleEB_, scaleHB_, scaleHO_, scheme_, AlCaHLTBitMon_QueryRunRegistry::string, tcut_, hcalRecHitTable_cff::time, time_, HcalTestNumberingScheme::unpackHcalIndex(), and ecaldqm::zside().

Referenced by analyze().

                                {
  //initialize
  int nhit[3];
  double etot[3], etotT[3];
  std::vector<unsigned int> ebID, hbID, hoID;
  std::vector<double> ebEtow, hbEtow, hoEtow;
  std::vector<double> ebEtowT, hbEtowT, hoEtowT;
  for (int k = 0; k < 3; k++) {
    nhit[k] = 0;
    etot[k] = 0;
    etotT[k] = 0;
  }
  eneInc_->Fill(eInc);
  etaInc_->Fill(etaInc);
  phiInc_->Fill(phiInc);

  double eHO17 = 0, eHO18 = 0, eHO17T = 0, eHO18T = 0;
  double eHOE17[15], eHOE18[15], eHOE17T[15], eHOE18T[15];
  for (int k = 0; k < 15; k++) {
    eHOE17[k] = eHOE18[k] = eHOE17T[k] = eHOE18T[k] = 0;
  }
  // Loop over containers
  for (int k = 0; k < 2; k++) {
    int nHit;
    if (k == 0) {
      nHit = ecalHits.size();
    } else {
      nHit = hcalHits.size();
    }
    for (int i = 0; i < nHit; i++) {
      double edep, time;
      int indx = -1;
      unsigned int id_;
      if (k == 0) {
        indx = 0;
        edep = ecalHits[i].energy();
        time = ecalHits[i].time();
        id_ = ecalHits[i].id();
      } else {
        edep = hcalHits[i].energy();
        time = hcalHits[i].time();
        id_ = hcalHits[i].id();
        int subdet, zside, depth, eta, phi, lay;
        if (scheme_) {
          HcalTestNumberingScheme::unpackHcalIndex(id_, subdet, zside, depth, eta, phi, lay);
        } else {
          HcalDetId id = HcalDetId(id_);
          subdet = id.subdet();
          zside = id.zside();
          depth = id.depth();
          eta = id.ietaAbs();
          phi = id.iphi();
          lay = -1;
        }
        edm::LogVerbatim("HitStudy") << "HOSimHitStudy:: Hit " << k << " Subdet:" << subdet << " zside:" << zside
                                     << " depth:" << depth << " layer:" << lay << " eta:" << eta << " phi:" << phi;
        if (subdet == static_cast<int>(HcalBarrel))
          indx = 1;
        else if (subdet == static_cast<int>(HcalOuter)) {
          indx = 2;
          if (lay == 18) {
            eHO17 += edep;
            if (time < tcut_)
              eHO17T += edep;
            if (eta >= 0 && eta < 15) {
              eHOE17[eta] += edep;
              if (time < tcut_)
                eHOE17T[eta] += edep;
            }
          } else {
            eHO18 += edep;
            if (time < tcut_)
              eHO18T += edep;
            if (eta >= 0 && eta < 15) {
              eHOE18[eta] += edep;
              if (time < tcut_)
                eHOE18T[eta] += edep;
            }
          }
        }
      }
      if (indx >= 0) {
        double edepT = edep;
        time_[indx]->Fill(time);
        edep_[indx]->Fill(edep);
        etot[indx] += edep;
        if (time < tcut_) {
          etotT[indx] += edep;
          edepT_[indx]->Fill(edep);
          edepT = 0;
        }
        nhit[indx]++;
        if (indx == 0) {
          bool ok = false;
          for (unsigned int j = 0; j < ebID.size(); j++) {
            if (id_ == ebID[j]) {
              ebEtow[j] += edep;
              ebEtowT[j] += edepT;
              ok = true;
              break;
            }
          }
          if (!ok) {
            ebID.push_back(id_);
            ebEtow.push_back(edep);
            ebEtowT.push_back(edepT);
          }
        } else if (indx == 1) {
          bool ok = false;
          for (unsigned int j = 0; j < hbID.size(); j++) {
            if (id_ == hbID[j]) {
              hbEtow[j] += edep;
              hbEtowT[j] += edepT;
              ok = true;
              break;
            }
          }
          if (!ok) {
            hbID.push_back(id_);
            hbEtow.push_back(edep);
            hbEtowT.push_back(edepT);
          }
        } else {
          bool ok = false;
          for (unsigned int j = 0; j < hoID.size(); j++) {
            if (id_ == hoID[j]) {
              hoEtow[j] += edep;
              hoEtowT[j] += edepT;
              ok = true;
              break;
            }
          }
          if (!ok) {
            hoID.push_back(id_);
            hoEtow.push_back(edep);
            hoEtowT.push_back(edepT);
          }
        }
      }
    }
  }

  // Now for towers and total energy deposits
  for (int k = 0; k < 3; k++) {
    hit_[k]->Fill(double(nhit[k]));
    edepZon_[k]->Fill(etot[k]);
    edepZonT_[k]->Fill(etotT[k]);
  }
  hitTow_[0]->Fill(double(ebEtow.size()));
  for (unsigned int i = 0; i < ebEtow.size(); i++) {
    edepTW_[0]->Fill(ebEtow[i]);
    edepTWT_[0]->Fill(ebEtowT[i]);
  }
  hitTow_[1]->Fill(double(hbEtow.size()));
  for (unsigned int i = 0; i < hbEtow.size(); i++) {
    edepTW_[1]->Fill(hbEtow[i]);
    edepTWT_[1]->Fill(hbEtowT[i]);
  }
  hitTow_[2]->Fill(double(hoEtow.size()));
  for (unsigned int i = 0; i < hoEtow.size(); i++) {
    edepTW_[2]->Fill(hoEtow[i]);
    edepTWT_[2]->Fill(hoEtowT[i]);
  }
  double eEB = scaleEB_ * etot[0];
  double eEBHB = eEB + scaleHB_ * etot[1];
  double eEBHBHO = eEBHB + scaleHB_ * scaleHO_ * etot[2];
  eEB_->Fill(eEB);
  eEBHB_->Fill(eEBHB);
  eEBHBHO_->Fill(eEBHBHO);
  double eEBT = scaleEB_ * etotT[0];
  double eEBHBT = eEBT + scaleHB_ * etotT[1];
  double eEBHBHOT = eEBHBT + scaleHB_ * scaleHO_ * etotT[2];
  eEBT_->Fill(eEBT);
  eEBHBT_->Fill(eEBHBT);
  eEBHBHOT_->Fill(eEBHBHOT);
  eHO1_->Fill(etaInc, eHO17 + eHO18);
  eHO2_->Fill(etaInc, phiInc, eHO17 + eHO18);
  eHO17_->Fill(etaInc, eHO17);
  eHO18_->Fill(etaInc, eHO18);
  eHO1T_->Fill(etaInc, eHO17T + eHO18T);
  eHO2T_->Fill(etaInc, phiInc, eHO17T + eHO18T);
  eHO17T_->Fill(etaInc, eHO17T);
  eHO18T_->Fill(etaInc, eHO18T);
  int nHO = 0, nHOT = 0;
  if (eHO17 > 0)
    nHO++;
  if (eHO17T > 0)
    nHOT++;
  if (eHO18 > 0)
    nHO++;
  if (eHO18T > 0)
    nHOT++;
  nHO1_->Fill(etaInc, (double)(nHO));
  nHO2_->Fill(etaInc, phiInc, (double)(nHO));
  nHO1T_->Fill(etaInc, (double)(nHOT));
  nHO2T_->Fill(etaInc, phiInc, (double)(nHOT));
  int ieta = 15;
  for (int k = 0; k < 15; ++k) {
    if (std::abs(etaInc) < 0.087 * (k + 1)) {
      ieta = k;
      break;
    }
  }
  if (ieta >= 0 && ieta < 15) {
    eHOE_[ieta]->Fill(eHO17 + eHO18);
    eHOE17_[ieta]->Fill(eHO17);
    eHOE18_[ieta]->Fill(eHO18);
    eHOET_[ieta]->Fill(eHO17T + eHO18T);
    eHOE17T_[ieta]->Fill(eHO17T);
    eHOE18T_[ieta]->Fill(eHO18T);
    eHOEta_[ieta]->Fill(eHOE17[ieta] + eHOE18[ieta]);
    eHOEta17_[ieta]->Fill(eHOE17[ieta]);
    eHOEta18_[ieta]->Fill(eHOE18[ieta]);
    nHOE1_[ieta]->Fill((double)(nHO));
    nHOE2_[ieta]->Fill(phiInc, (double)(nHO));
    nHOE1T_[ieta]->Fill((double)(nHOT));
    nHOE2T_[ieta]->Fill(phiInc, (double)(nHOT));
    int nHOE = 0, nHOET = 0;
    if (eHOE17[ieta] > 0)
      nHOE++;
    if (eHOE17T[ieta] > 0)
      nHOET++;
    if (eHOE18[ieta] > 0)
      nHOE++;
    if (eHOE18T[ieta] > 0)
      nHOET++;
    nHOEta1_[ieta]->Fill((double)(nHOE));
    nHOEta2_[ieta]->Fill(phiInc, (double)(nHOE));
    nHOEta1T_[ieta]->Fill((double)(nHOET));
    nHOEta2T_[ieta]->Fill(phiInc, (double)(nHOET));
  }

  edm::LogVerbatim("HitStudy") << "HOSimHitStudy::analyzeHits: Hits in EB " << nhit[0] << " in " << ebEtow.size()
                               << " towers with total E " << etot[0] << "|" << etotT[0]
                               << "\n                            Hits in HB " << nhit[1] << " in " << hbEtow.size()
                               << " towers with total E " << etot[1] << "|" << etotT[1]
                               << "\n                            Hits in HO " << nhit[2] << " in " << hoEtow.size()
                               << " towers with total E " << etot[2] << "|" << etotT[2]
                               << "\n                            Energy in EB " << eEB << "|" << eEBT << " with HB "
                               << eEBHB << "|" << eEBHBT << " and with HO " << eEBHBHO << "|" << eEBHBHOT
                               << "\n                            E in HO layers " << eHO17 << "|" << eHO17T << " "
                               << eHO18 << "|" << eHO18T << " number of HO hits " << nHO << "|" << nHOT;

  if (eEBHBHO > 0.75 * maxEnergy_ && print_) {
    edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Event with excessive energy: EB = " << eEB << " EB+HB = " << eEBHB
                                 << " EB+HB+HO = " << eEBHBHO;
    const std::string Dets[3] = {"EB", "HB", "HO"};
    for (int k = 0; k < 2; k++) {
      int nHit;
      if (k == 0) {
        nHit = ecalHits.size();
      } else {
        nHit = hcalHits.size();
      }
      for (int i = 0; i < nHit; i++) {
        double edep, time;
        int indx = -1;
        unsigned int id_;
        int ieta, iphi, depth = 0;
        if (k == 0) {
          indx = 0;
          edep = ecalHits[i].energy();
          time = ecalHits[i].time();
          id_ = ecalHits[i].id();
          EBDetId id = EBDetId(id_);
          ieta = id.ieta();
          iphi = id.iphi();
        } else {
          indx = -1;
          edep = hcalHits[i].energy();
          time = hcalHits[i].time();
          id_ = hcalHits[i].id();
          HcalDetId id = HcalDetId(id_);
          int subdet = id.subdet();
          if (subdet == static_cast<int>(HcalBarrel))
            indx = 1;
          else if (subdet == static_cast<int>(HcalOuter))
            indx = 2;
          ieta = id.ieta();
          iphi = id.iphi();
          depth = id.depth();
        }
        if (indx >= 0) {
          edm::LogVerbatim("HitStudy") << "HOSimHitStudy::" << Dets[indx] << " " << i << std::hex << id_ << std::dec
                                       << " (" << ieta << "|" << iphi << "|" << depth << ") " << std::setw(8) << edep
                                       << " " << std::setw(8) << time;
        }
      }
    }
  }
}

beginJob()

void HOSimHitStudy::beginJob ( void  )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 119 of file HOSimHitStudy.cc.

References edep_, edepT_, edepTW_, edepTWT_, edepZon_, edepZonT_, eEB_, eEBHB_, eEBHBHO_, eEBHBHOT_, eEBHBT_, eEBT_, eHO17_, eHO17T_, eHO18_, eHO18T_, eHO1_, eHO1T_, eHO2_, eHO2T_, eHOE17_, eHOE17T_, eHOE18_, eHOE18T_, eHOE_, eHOET_, eHOEta17_, eHOEta17T_, eHOEta18_, eHOEta18T_, eHOEta_, eHOEtaT_, eneInc_, etaInc_, hit_, hitTow_, mps_fire::i, createfilelist::int, maxEnergy_, Skims_PA_cff::name, nHO1_, nHO1T_, nHO2_, nHO2T_, nHOE1_, nHOE1T_, nHOE2_, nHOE2T_, nHOEta1_, nHOEta1T_, nHOEta2_, nHOEta2T_, phiInc_, AlCaHLTBitMon_QueryRunRegistry::string, tcut_, compare::tfile, time_, runGCPTkAlMap::title, and L1TOccupancyClient_cfi::ymax.

                             {
  edm::Service<TFileService> tfile;

  if (!tfile.isAvailable())
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
  std::string dets[3] = {"EB", "HB", "HO"};
  char name[60], title[100];
  double ymax = maxEnergy_;
  sprintf(title, "Incident Energy (GeV)");
  eneInc_ = tfile->make<TH1F>("EneInc", title, 1000, 0., ymax);
  eneInc_->GetXaxis()->SetTitle(title);
  eneInc_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Incident #eta");
  etaInc_ = tfile->make<TH1F>("EtaInc", title, 200, -5., 5.);
  etaInc_->GetXaxis()->SetTitle(title);
  etaInc_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Incident #phi");
  phiInc_ = tfile->make<TH1F>("PhiInc", title, 200, -3.1415926, 3.1415926);
  phiInc_->GetXaxis()->SetTitle(title);
  phiInc_->GetYaxis()->SetTitle("Events");
  int itcut = (int)(tcut_);
  for (int i = 0; i < 3; i++) {
    sprintf(name, "Hit%d", i);
    sprintf(title, "Number of hits in %s", dets[i].c_str());
    hit_[i] = tfile->make<TH1F>(name, title, 1000, 0., 20000.);
    hit_[i]->GetXaxis()->SetTitle(title);
    hit_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "Time%d", i);
    sprintf(title, "Time of the hit (ns) in %s", dets[i].c_str());
    time_[i] = tfile->make<TH1F>(name, title, 1200, 0., 1200.);
    time_[i]->GetXaxis()->SetTitle(title);
    time_[i]->GetYaxis()->SetTitle("Hits");
    if (i > 0)
      ymax = 1.0;
    else
      ymax = 50.0;
    sprintf(name, "Edep%d", i);
    sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
    edep_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
    edep_[i]->GetXaxis()->SetTitle(title);
    edep_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "EdepT%d", i);
    sprintf(title, "Energy deposit (GeV) in %s for t < %d ns", dets[i].c_str(), itcut);
    edepT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
    edepT_[i]->GetXaxis()->SetTitle(title);
    edepT_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "HitTow%d", i);
    sprintf(title, "Number of towers with hits in %s", dets[i].c_str());
    hitTow_[i] = tfile->make<TH1F>(name, title, 1000, 0., 20000.);
    hitTow_[i]->GetXaxis()->SetTitle(title);
    hitTow_[i]->GetYaxis()->SetTitle("Events");
    if (i > 0)
      ymax = 0.05 * maxEnergy_;
    else
      ymax = maxEnergy_;
    sprintf(name, "EdepTW%d", i);
    sprintf(title, "Energy deposit (GeV) in %s Tower", dets[i].c_str());
    edepTW_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
    edepTW_[i]->GetXaxis()->SetTitle(title);
    edepTW_[i]->GetYaxis()->SetTitle("Towers");
    sprintf(name, "EdepTWT%d", i);
    sprintf(title, "Energy deposit (GeV) in %s Tower for t < %d ns", dets[i].c_str(), itcut);
    edepTWT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
    edepTWT_[i]->GetXaxis()->SetTitle(title);
    edepTWT_[i]->GetYaxis()->SetTitle("Towers");
    sprintf(name, "EdepZone%d", i);
    sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
    edepZon_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
    edepZon_[i]->GetXaxis()->SetTitle(title);
    edepZon_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EdepZoneT%d", i);
    sprintf(title, "Energy deposit (GeV) in %s for t < %d ns", dets[i].c_str(), itcut);
    edepZonT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
    edepZonT_[i]->GetXaxis()->SetTitle(title);
    edepZonT_[i]->GetYaxis()->SetTitle("Events");
  }
  sprintf(title, "Energy Measured in EB (GeV)");
  eEB_ = tfile->make<TH1F>("EEB", title, 5000, 0., maxEnergy_);
  eEB_->GetXaxis()->SetTitle(title);
  eEB_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Energy Measured in EB+HB (GeV)");
  eEBHB_ = tfile->make<TH1F>("EEBHB", title, 5000, 0., maxEnergy_);
  eEBHB_->GetXaxis()->SetTitle(title);
  eEBHB_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Energy Measured in EB+HB+HO (GeV)");
  eEBHBHO_ = tfile->make<TH1F>("EEBHBHO", title, 5000, 0., maxEnergy_);
  eEBHBHO_->GetXaxis()->SetTitle(title);
  eEBHBHO_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Energy Measured in EB (GeV) for t < %d ns", itcut);
  eEBT_ = tfile->make<TH1F>("EEBT", title, 5000, 0., maxEnergy_);
  eEBT_->GetXaxis()->SetTitle(title);
  eEBT_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Energy Measured in EB+HB (GeV) for t < %d ns", itcut);
  eEBHBT_ = tfile->make<TH1F>("EEBHBT", title, 5000, 0., maxEnergy_);
  eEBHBT_->GetXaxis()->SetTitle(title);
  eEBHBT_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Energy Measured in EB+HB+HO (GeV) for t < %d ns", itcut);
  eEBHBHOT_ = tfile->make<TH1F>("EEBHBHOT", title, 5000, 0., maxEnergy_);
  eEBHBHOT_->GetXaxis()->SetTitle(title);
  eEBHBHOT_->GetYaxis()->SetTitle("Events");
  sprintf(title, "SimHit energy in HO");
  eHO1_ = tfile->make<TProfile>("EHO1", title, 30, -1.305, 1.305);
  eHO1_->GetXaxis()->SetTitle(title);
  eHO1_->GetYaxis()->SetTitle("Events");
  eHO2_ = tfile->make<TProfile2D>("EHO2", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
  eHO2_->GetXaxis()->SetTitle(title);
  eHO2_->GetYaxis()->SetTitle("Events");
  sprintf(title, "SimHit energy in HO Layer 17");
  eHO17_ = tfile->make<TProfile>("EHO17", title, 30, -1.305, 1.305);
  eHO17_->GetXaxis()->SetTitle(title);
  eHO17_->GetYaxis()->SetTitle("Events");
  sprintf(title, "SimHit energy in HO Layer 18");
  eHO18_ = tfile->make<TProfile>("EHO18", title, 30, -1.305, 1.305);
  eHO18_->GetXaxis()->SetTitle(title);
  eHO18_->GetYaxis()->SetTitle("Events");
  sprintf(title, "SimHit energy in HO for t < %d ns", itcut);
  eHO1T_ = tfile->make<TProfile>("EHO1T", title, 30, -1.305, 1.305);
  eHO1T_->GetXaxis()->SetTitle(title);
  eHO1T_->GetYaxis()->SetTitle("Events");
  eHO2T_ = tfile->make<TProfile2D>("EHO2T", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
  eHO2T_->GetXaxis()->SetTitle(title);
  eHO2T_->GetYaxis()->SetTitle("Events");
  sprintf(title, "SimHit energy in HO Layer 17 for t < %d ns", itcut);
  eHO17T_ = tfile->make<TProfile>("EHO17T", title, 30, -1.305, 1.305);
  eHO17T_->GetXaxis()->SetTitle(title);
  eHO17T_->GetYaxis()->SetTitle("Events");
  sprintf(title, "SimHit energy in HO Layer 18 for t < %d ns", itcut);
  eHO18T_ = tfile->make<TProfile>("EHO18T", title, 30, -1.305, 1.305);
  eHO18T_->GetXaxis()->SetTitle(title);
  eHO18T_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Number of layers hit in HO");
  nHO1_ = tfile->make<TProfile>("NHO1", title, 30, -1.305, 1.305);
  nHO1_->GetXaxis()->SetTitle(title);
  nHO1_->GetYaxis()->SetTitle("Events");
  nHO2_ = tfile->make<TProfile2D>("NHO2", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
  nHO2_->GetXaxis()->SetTitle(title);
  nHO2_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Number of layers hit in HO for t < %d ns", itcut);
  nHO1T_ = tfile->make<TProfile>("NHO1T", title, 30, -1.305, 1.305);
  nHO1T_->GetXaxis()->SetTitle(title);
  nHO1T_->GetYaxis()->SetTitle("Events");
  nHO2T_ = tfile->make<TProfile2D>("NHO2T", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
  nHO2T_->GetXaxis()->SetTitle(title);
  nHO2T_->GetYaxis()->SetTitle("Events");
  for (int i = 0; i < 15; i++) {
    sprintf(name, "EHOE%d", i + 1);
    sprintf(title, "SimHit energy in HO (Beam in #eta=%d bin)", i + 1);
    eHOE_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOE_[i]->GetXaxis()->SetTitle(title);
    eHOE_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOE17%d", i + 1);
    sprintf(title, "SimHit energy in Layer 17 (Beam in #eta=%d bin)", i + 1);
    eHOE17_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOE17_[i]->GetXaxis()->SetTitle(title);
    eHOE17_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOE18%d", i + 1);
    sprintf(title, "SimHit energy in Layer 18 (Beam in #eta=%d bin)", i + 1);
    eHOE18_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOE18_[i]->GetXaxis()->SetTitle(title);
    eHOE18_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOE%dT", i + 1);
    sprintf(title, "SimHit energy in HO (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
    eHOET_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOET_[i]->GetXaxis()->SetTitle(title);
    eHOET_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOE17%dT", i + 1);
    sprintf(title, "SimHit energy in Layer 17 (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
    eHOE17T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOE17T_[i]->GetXaxis()->SetTitle(title);
    eHOE17T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOE18%dT", i + 1);
    sprintf(title, "SimHit energy in Layer 18 (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
    eHOE18T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOE18T_[i]->GetXaxis()->SetTitle(title);
    eHOE18T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOEta%d", i + 1);
    sprintf(title, "SimHit energy in HO #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
    eHOEta_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOEta_[i]->GetXaxis()->SetTitle(title);
    eHOEta_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOEta17%d", i + 1);
    sprintf(title, "SimHit energy in Layer 17 #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
    eHOEta17_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOEta17_[i]->GetXaxis()->SetTitle(title);
    eHOEta17_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOEta18%d", i + 1);
    sprintf(title, "SimHit energy in Layer 18 #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
    eHOEta18_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOEta18_[i]->GetXaxis()->SetTitle(title);
    eHOEta18_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOEta%dT", i + 1);
    sprintf(title, "SimHit energy in HO #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
    eHOEtaT_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOEtaT_[i]->GetXaxis()->SetTitle(title);
    eHOEtaT_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOEta17%dT", i + 1);
    sprintf(title, "SimHit energy in Layer 17 #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
    eHOEta17T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOEta17T_[i]->GetXaxis()->SetTitle(title);
    eHOEta17T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EHOEta18%dT", i + 1);
    sprintf(title, "SimHit energy in Layer 18 #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
    eHOEta18T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
    eHOEta18T_[i]->GetXaxis()->SetTitle(title);
    eHOEta18T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOE1%d", i + 1);
    sprintf(title, "Number of layers hit in HO  (Beam in #eta=%d bin)", i + 1);
    nHOE1_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
    nHOE1_[i]->GetXaxis()->SetTitle(title);
    nHOE1_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOE2%d", i + 1);
    nHOE2_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
    nHOE2_[i]->GetXaxis()->SetTitle(title);
    nHOE2_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOE1%dT", i + 1);
    sprintf(title, "Number of layers hit in HO (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
    nHOE1T_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
    nHOE1T_[i]->GetXaxis()->SetTitle(title);
    nHOE1T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOE2%dT", i + 1);
    nHOE2T_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
    nHOE2T_[i]->GetXaxis()->SetTitle(title);
    nHOE2T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOEta1%d", i + 1);
    sprintf(title, "Number of layers hit in HO #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
    nHOEta1_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
    nHOEta1_[i]->GetXaxis()->SetTitle(title);
    nHOEta1_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOEta2%d", i + 1);
    nHOEta2_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
    nHOEta2_[i]->GetXaxis()->SetTitle(title);
    nHOEta2_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOEta1%dT", i + 1);
    sprintf(title, "Number of layers hit in HO #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
    nHOEta1T_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
    nHOEta1T_[i]->GetXaxis()->SetTitle(title);
    nHOEta1T_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "NHOEta2%dT", i + 1);
    nHOEta2T_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
    nHOEta2T_[i]->GetXaxis()->SetTitle(title);
    nHOEta2T_[i]->GetYaxis()->SetTitle("Events");
  }
}

beginRun()

void HOSimHitStudy::beginRun ( edm::Run const &  , edm::EventSetup const &    )

inlineoverrideprotected

Definition at line 49 of file HOSimHitStudy.cc.

{}

endJob()

void HOSimHitStudy::endJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 50 of file HOSimHitStudy.cc.

{}

endRun()

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

inlineoverrideprotected

Definition at line 51 of file HOSimHitStudy.cc.

{}

fillDescriptions()

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

static

Definition at line 103 of file HOSimHitStudy.cc.

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

                                                                               {
  edm::ParameterSetDescription desc;
  std::vector<std::string> labels = {"EcalHitsEB", "HcalHits"};
  desc.addUntracked<std::string>("SourceLabel", "generatorSmeared");
  desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
  desc.addUntracked<std::vector<std::string>>("HitCollection", labels);
  desc.addUntracked<double>("MaxEnergy", 50.0);
  desc.addUntracked<double>("ScaleEB", 1.02);
  desc.addUntracked<double>("ScaleHB", 104.4);
  desc.addUntracked<double>("ScaleHO", 2.33);
  desc.addUntracked<double>("TimeCut", 100.0);
  desc.addUntracked<bool>("PrintExcessEnergy", true);
  desc.addUntracked<bool>("TestNumbering", false);
  descriptions.add("HOSimHitStudy", desc);
}

Member Data Documentation

ecalHits

std::vector<PCaloHit> HOSimHitStudy::ecalHits

private

Definition at line 64 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

edep_

TH1F* HOSimHitStudy::edep_[3]

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepT_

TH1F* HOSimHitStudy::edepT_[3]

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepTW_

TH1F * HOSimHitStudy::edepTW_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepTWT_

TH1F * HOSimHitStudy::edepTWT_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepZon_

TH1F* HOSimHitStudy::edepZon_[3]

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepZonT_

TH1F * HOSimHitStudy::edepZonT_[3]

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eEB_

TH1F * HOSimHitStudy::eEB_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eEBHB_

TH1F * HOSimHitStudy::eEBHB_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eEBHBHO_

TH1F * HOSimHitStudy::eEBHBHO_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eEBHBHOT_

TH1F * HOSimHitStudy::eEBHBHOT_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eEBHBT_

TH1F * HOSimHitStudy::eEBHBT_

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eEBT_

TH1F * HOSimHitStudy::eEBT_

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO17_

TProfile * HOSimHitStudy::eHO17_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO17T_

TProfile * HOSimHitStudy::eHO17T_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO18_

TProfile * HOSimHitStudy::eHO18_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO18T_

TProfile * HOSimHitStudy::eHO18T_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO1_

TProfile* HOSimHitStudy::eHO1_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO1T_

TProfile * HOSimHitStudy::eHO1T_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO2_

TProfile2D* HOSimHitStudy::eHO2_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHO2T_

TProfile2D * HOSimHitStudy::eHO2T_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOE17_

TH1F* HOSimHitStudy::eHOE17_[15]

private

Definition at line 69 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOE17T_

TH1F* HOSimHitStudy::eHOE17T_[15]

private

Definition at line 70 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOE18_

TH1F * HOSimHitStudy::eHOE18_[15]

private

Definition at line 69 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOE18T_

TH1F * HOSimHitStudy::eHOE18T_[15]

private

Definition at line 70 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOE_

TH1F * HOSimHitStudy::eHOE_[15]

private

Definition at line 69 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOET_

TH1F * HOSimHitStudy::eHOET_[15]

private

Definition at line 70 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOEta17_

TH1F* HOSimHitStudy::eHOEta17_[15]

private

Definition at line 71 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOEta17T_

TH1F* HOSimHitStudy::eHOEta17T_[15]

private

Definition at line 72 of file HOSimHitStudy.cc.

Referenced by beginJob().

eHOEta18_

TH1F * HOSimHitStudy::eHOEta18_[15]

private

Definition at line 71 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOEta18T_

TH1F * HOSimHitStudy::eHOEta18T_[15]

private

Definition at line 72 of file HOSimHitStudy.cc.

Referenced by beginJob().

eHOEta_

TH1F * HOSimHitStudy::eHOEta_[15]

private

Definition at line 71 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eHOEtaT_

TH1F * HOSimHitStudy::eHOEtaT_[15]

private

Definition at line 72 of file HOSimHitStudy.cc.

Referenced by beginJob().

eInc

double HOSimHitStudy::eInc

private

Definition at line 78 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

eneInc_

TH1F * HOSimHitStudy::eneInc_

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

etaInc

double HOSimHitStudy::etaInc

private

Definition at line 78 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

etaInc_

TH1F * HOSimHitStudy::etaInc_

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

g4Label_

const std::string HOSimHitStudy::g4Label_

private

Definition at line 57 of file HOSimHitStudy.cc.

Referenced by HOSimHitStudy().

hcalHits

std::vector<PCaloHit> HOSimHitStudy::hcalHits

private

Definition at line 64 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

hit_

TH1F* HOSimHitStudy::hit_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

hitLab_

const std::vector<std::string> HOSimHitStudy::hitLab_

private

Definition at line 58 of file HOSimHitStudy.cc.

Referenced by analyze(), and HOSimHitStudy().

hitTow_

TH1F * HOSimHitStudy::hitTow_[3]

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

maxEnergy_

const double HOSimHitStudy::maxEnergy_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHO1_

TProfile* HOSimHitStudy::nHO1_

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHO1T_

TProfile * HOSimHitStudy::nHO1T_

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHO2_

TProfile2D * HOSimHitStudy::nHO2_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHO2T_

TProfile2D * HOSimHitStudy::nHO2T_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOE1_

TH1F* HOSimHitStudy::nHOE1_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOE1T_

TH1F * HOSimHitStudy::nHOE1T_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOE2_

TProfile * HOSimHitStudy::nHOE2_[15]

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOE2T_

TProfile * HOSimHitStudy::nHOE2T_[15]

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOEta1_

TH1F * HOSimHitStudy::nHOEta1_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOEta1T_

TH1F * HOSimHitStudy::nHOEta1T_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOEta2_

TProfile* HOSimHitStudy::nHOEta2_[15]

private

Definition at line 76 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

nHOEta2T_

TProfile * HOSimHitStudy::nHOEta2T_[15]

private

Definition at line 76 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

phiInc

double HOSimHitStudy::phiInc

private

Definition at line 78 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

phiInc_

TH1F * HOSimHitStudy::phiInc_

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

print_

const bool HOSimHitStudy::print_

private

Definition at line 61 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

scaleEB_

const double HOSimHitStudy::scaleEB_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

scaleHB_

const double HOSimHitStudy::scaleHB_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

scaleHO_

const double HOSimHitStudy::scaleHO_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

scheme_

const bool HOSimHitStudy::scheme_

private

Definition at line 61 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

tcut_

const double HOSimHitStudy::tcut_

private

Definition at line 60 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

time_

TH1F * HOSimHitStudy::time_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

tok_evt_

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

private

Definition at line 62 of file HOSimHitStudy.cc.

Referenced by analyze().

toks_calo_

const std::vector<edm::EDGetTokenT<edm::PCaloHitContainer> > HOSimHitStudy::toks_calo_

private

Definition at line 63 of file HOSimHitStudy.cc.

Referenced by analyze().