CaloSimHitAnalysis Class Reference

Inheritance diagram for CaloSimHitAnalysis:

Public Member Functions
	CaloSimHitAnalysis (const edm::ParameterSet &ps)
	~CaloSimHitAnalysis () override
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 final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const 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
bool	wantsStreamRuns () const
	~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 > &, int)
void	analyzePassiveHits (std::vector< PassiveHit > &hits)
void	beginJob () override
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 int	allSteps_
const CaloGeometry *	caloGeometry_
const std::vector< std::string >	detNames_
std::map< int, unsigned int >	etaDepth_
const std::string	g4Label_
TH1F *	h_edep1_ [nCalo_]
TH1F *	h_edep_ [nCalo_]
TH1F *	h_edepEM_ [nCalo_]
TH1F *	h_edepHad_ [nCalo_]
TH1F *	h_edepp_
TH1F *	h_edepT1_ [nCalo_]
TH1F *	h_edepT_ [nCalo_]
std::vector< TH1F * >	h_edepTk_
TH1F *	h_eta_ [nCalo_]
TH2F *	h_etaphi_
TH1F *	h_etot_ [nCalo_]
TH1F *	h_etotg_ [nCalo_]
TH1F *	h_hit_ [nCalo_]
TH1F *	h_hitp_
TH1F *	h_phi_ [nCalo_]
TH1F *	h_rr_ [nCalo_]
TH2F *	h_rz1_
TH2F *	h_rz_
std::vector< TH2F * >	h_rzH_
TH1F *	h_stepp_
TH1F *	h_time_ [nCalo_]
TH1F *	h_timep_
TH1F *	h_timeT_ [nCalo_]
std::vector< TH1F * >	h_timeTk_
TH1F *	h_trackp_
TH1F *	h_zz_ [nCalo_]
const HcalGeometry *	hcalGeom_
const std::vector< std::string >	hitLab_
const double	maxEnergy_
const double	maxTime_
const bool	passive_
const double	tCut_
const bool	testNumber_
const std::vector< double >	timeSliceUnit_
const double	tMax_
const edm::EDGetTokenT< edm::PassiveHitContainer >	tok_passive_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tokGeom_
const std::vector< edm::EDGetTokenT< edm::PCaloHitContainer > >	toks_calo_
const double	tScale_

Static Private Attributes
static constexpr int	nCalo_ = 6

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 45 of file CaloSimHitAnalysis.cc.

Constructor & Destructor Documentation

CaloSimHitAnalysis()

CaloSimHitAnalysis::CaloSimHitAnalysis

(

const edm::ParameterSet &

ps

)

Definition at line 87 of file CaloSimHitAnalysis.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")),
      timeSliceUnit_(ps.getParameter<std::vector<double> >("timeSliceUnit")),
      maxEnergy_(ps.getUntrackedParameter<double>("maxEnergy", 250.0)),
      maxTime_(ps.getUntrackedParameter<double>("maxTime", 1000.0)),
      tMax_(ps.getUntrackedParameter<double>("timeCut", 100.0)),
      tScale_(ps.getUntrackedParameter<double>("timeScale", 1.0)),
      tCut_(ps.getUntrackedParameter<double>("timeThreshold", 15.0)),
      testNumber_(ps.getUntrackedParameter<bool>("testNumbering", false)),
      passive_(ps.getUntrackedParameter<bool>("passiveHits", false)),
      allSteps_(ps.getUntrackedParameter<int>("allSteps", 100)),
      detNames_(ps.getUntrackedParameter<std::vector<std::string> >("detNames")),
      tokGeom_(esConsumes<CaloGeometry, CaloGeometryRecord>()),
      toks_calo_{edm::vector_transform(hitLab_,
                                       [this](const std::string& name) {
                                         return consumes<edm::PCaloHitContainer>(edm::InputTag{g4Label_, name});
                                       })},
      tok_passive_(consumes<edm::PassiveHitContainer>(edm::InputTag(g4Label_, "AllPassiveHits"))) {
  usesResource(TFileService::kSharedResource);

  edm::LogVerbatim("HitStudy") << "Module Label: " << g4Label_ << "   Hits|timeSliceUnit:";
  for (unsigned int i = 0; i < hitLab_.size(); i++)
    edm::LogVerbatim("HitStudy") << "[" << i << "] " << hitLab_[i] << " " << timeSliceUnit_[i];
  edm::LogVerbatim("HitStudy") << "Passive Hits " << passive_ << " from AllPassiveHits";
  edm::LogVerbatim("HitStudy") << "maxEnergy: " << maxEnergy_ << " maxTime: " << maxTime_ << " tMax: " << tMax_;
  for (unsigned int k = 0; k < detNames_.size(); ++k)
    edm::LogVerbatim("HitStudy") << "Detector[" << k << "] " << detNames_[k];

  edm::Service<TFileService> tfile;
  if (!tfile.isAvailable())
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
  char name[29], title[120];
  std::string dets[nCalo_] = {"EB", "EE", "HB", "HE", "HO", "HF"};
  for (int i = 0; i < nCalo_; i++) {
    sprintf(name, "Hit%d", i);
    sprintf(title, "Number of hits in %s", dets[i].c_str());
    h_hit_[i] = tfile->make<TH1F>(name, title, 100, 0., 20000.);
    h_hit_[i]->GetXaxis()->SetTitle(title);
    h_hit_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "Time%d", i);
    sprintf(title, "Time of the hit (ns) in %s", dets[i].c_str());
    h_time_[i] = tfile->make<TH1F>(name, title, 100, 0., 200.);
    h_time_[i]->GetXaxis()->SetTitle(title);
    h_time_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "TimeT%d", i);
    sprintf(title, "Time of each hit (ns) in %s", dets[i].c_str());
    h_timeT_[i] = tfile->make<TH1F>(name, title, 100, 0., 200.);
    h_timeT_[i]->GetXaxis()->SetTitle(title);
    h_timeT_[i]->GetYaxis()->SetTitle("Hits");
    double ymax = (i > 1) ? 0.01 : 0.1;
    double ymx0 = (i > 1) ? 0.0025 : 0.025;
    sprintf(name, "Edep%d", i);
    sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
    h_edep_[i] = tfile->make<TH1F>(name, title, 100, 0., ymax);
    h_edep_[i]->GetXaxis()->SetTitle(title);
    h_edep_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "EdepT%d", i);
    sprintf(title, "Energy deposit (GeV) of each hit in %s", dets[i].c_str());
    h_edepT_[i] = tfile->make<TH1F>(name, title, 100, 0., ymx0);
    h_edepT_[i]->GetXaxis()->SetTitle(title);
    h_edepT_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "EdepEM%d", i);
    sprintf(title, "Energy deposit (GeV) by EM particles in %s", dets[i].c_str());
    h_edepEM_[i] = tfile->make<TH1F>(name, title, 100, 0., ymx0);
    h_edepEM_[i]->GetXaxis()->SetTitle(title);
    h_edepEM_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "EdepHad%d", i);
    sprintf(title, "Energy deposit (GeV) by hadrons in %s", dets[i].c_str());
    h_edepHad_[i] = tfile->make<TH1F>(name, title, 100, 0., ymx0);
    h_edepHad_[i]->GetXaxis()->SetTitle(title);
    h_edepHad_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "Edep15%d", i);
    sprintf(title, "Energy deposit (GeV) for T > %4.0f ns in %s", tCut_, dets[i].c_str());
    h_edep1_[i] = tfile->make<TH1F>(name, title, 100, 0., ymax);
    h_edep1_[i]->GetXaxis()->SetTitle(title);
    h_edep1_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "EdepT15%d", i);
    sprintf(title, "Energy deposit (GeV) of each hit for T > %4.0f in %s", tCut_, dets[i].c_str());
    h_edepT1_[i] = tfile->make<TH1F>(name, title, 100, 0., ymx0);
    h_edepT1_[i]->GetXaxis()->SetTitle(title);
    h_edepT1_[i]->GetYaxis()->SetTitle("Hits");
    ymax = (i > 1) ? 1.0 : maxEnergy_;
    sprintf(name, "Etot%d", i);
    sprintf(title, "Total energy deposit (GeV) in %s", dets[i].c_str());
    h_etot_[i] = tfile->make<TH1F>(name, title, 50, 0., ymax);
    h_etot_[i]->GetXaxis()->SetTitle(title);
    h_etot_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "EtotG%d", i);
    sprintf(title, "Total energy deposit (GeV) in %s (t < 100 ns)", dets[i].c_str());
    h_etotg_[i] = tfile->make<TH1F>(name, title, 50, 0., ymax);
    h_etotg_[i]->GetXaxis()->SetTitle(title);
    h_etotg_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "rr%d", i);
    sprintf(title, "R of hit point (cm) in %s", dets[i].c_str());
    h_rr_[i] = tfile->make<TH1F>(name, title, 100, 0., 250.);
    h_rr_[i]->GetXaxis()->SetTitle(title);
    h_rr_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "zz%d", i);
    sprintf(title, "z of hit point (cm) in %s", dets[i].c_str());
    h_zz_[i] = tfile->make<TH1F>(name, title, 240, -600., 600.);
    h_zz_[i]->GetXaxis()->SetTitle(title);
    h_zz_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "eta%d", i);
    sprintf(title, "#eta of hit point in %s", dets[i].c_str());
    h_eta_[i] = tfile->make<TH1F>(name, title, 100, -5.0, 5.0);
    h_eta_[i]->GetXaxis()->SetTitle(title);
    h_eta_[i]->GetYaxis()->SetTitle("Hits");
    sprintf(name, "phi%d", i);
    sprintf(title, "#phi of hit point in %s", dets[i].c_str());
    h_phi_[i] = tfile->make<TH1F>(name, title, 100, -M_PI, M_PI);
    h_phi_[i]->GetXaxis()->SetTitle(title);
    h_phi_[i]->GetYaxis()->SetTitle("Hits");
  }
  sprintf(title, "R vs Z of hit point");
  h_rz_ = tfile->make<TH2F>("rz", title, 120, 0., 600., 100, 0., 250.);
  h_rz_->GetXaxis()->SetTitle("z (cm)");
  h_rz_->GetYaxis()->SetTitle("R (cm)");
  sprintf(title, "R vs Z of hit point for hits with T > %4.0f ns", tCut_);
  h_rz1_ = tfile->make<TH2F>("rz2", title, 120, 0., 600., 100, 0., 250.);
  h_rz1_->GetXaxis()->SetTitle("z (cm)");
  h_rz1_->GetYaxis()->SetTitle("R (cm)");
  sprintf(title, "#phi vs #eta of hit point");
  h_etaphi_ = tfile->make<TH2F>("etaphi", title, 100, 0., 5., 100, 0., M_PI);
  h_etaphi_->GetXaxis()->SetTitle("#eta");
  h_etaphi_->GetYaxis()->SetTitle("#phi");

  if (passive_) {
    h_hitp_ = tfile->make<TH1F>("hitp", "All Steps", 100, 0.0, 20000.0);
    h_hitp_->GetXaxis()->SetTitle("Hits");
    h_hitp_->GetYaxis()->SetTitle("Events");
    h_trackp_ = tfile->make<TH1F>("trackp", "All Steps", 100, 0.0, 200000.0);
    h_hitp_->GetXaxis()->SetTitle("Tracks");
    h_hitp_->GetYaxis()->SetTitle("Events");
    h_edepp_ = tfile->make<TH1F>("edepp", "All Steps", 100, 0.0, 50.0);
    h_edepp_->GetXaxis()->SetTitle("Energy Deposit (MeV)");
    h_edepp_->GetYaxis()->SetTitle("Hits");
    h_timep_ = tfile->make<TH1F>("timep", "All Steps", 100, 0.0, 100.0);
    h_timep_->GetXaxis()->SetTitle("Hits");
    h_timep_->GetYaxis()->SetTitle("Hit Time (ns)");
    h_stepp_ = tfile->make<TH1F>("stepp", "All Steps", 1000, 0.0, 100.0);
    h_stepp_->GetXaxis()->SetTitle("Hits");
    h_stepp_->GetYaxis()->SetTitle("Step length (cm)");
    for (unsigned int k = 0; k < detNames_.size(); ++k) {
      sprintf(name, "edept%d", k);
      sprintf(title, "Energy Deposit (MeV) in %s", detNames_[k].c_str());
      h_edepTk_.emplace_back(tfile->make<TH1F>(name, title, 100, 0.0, 1.0));
      h_edepTk_.back()->GetYaxis()->SetTitle("Hits");
      h_edepTk_.back()->GetXaxis()->SetTitle(title);
      sprintf(name, "timet%d", k);
      sprintf(title, "Hit Time (ns) in %s", detNames_[k].c_str());
      h_timeTk_.emplace_back(tfile->make<TH1F>(name, title, 100, 0.0, 100.0));
      h_timeTk_.back()->GetYaxis()->SetTitle("Hits");
      h_timeTk_.back()->GetXaxis()->SetTitle(title);
    }
    if ((allSteps_ / 100) % 10 > 0) {
      for (int eta = 1; eta <= 29; ++eta) {
        int dmax = (eta < 16) ? 4 : 7;
        for (int depth = 1; depth <= dmax; ++depth) {
          sprintf(name, "Eta%dDepth%d", eta, depth);
          sprintf(title, "R vs Z (#eta = %d, depth = %d)", eta, depth);
          etaDepth_[eta * 100 + depth] = h_rzH_.size();
          h_rzH_.emplace_back(tfile->make<TH2F>(name, title, 120, 0., 600., 100, 0., 250.));
          h_rzH_.back()->GetXaxis()->SetTitle("z (cm)");
          h_rzH_.back()->GetYaxis()->SetTitle("R (cm)");
        }
      }
    }
  }
}

~CaloSimHitAnalysis()

CaloSimHitAnalysis::~CaloSimHitAnalysis ( )

inlineoverride

Definition at line 48 of file CaloSimHitAnalysis.cc.

{}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 279 of file CaloSimHitAnalysis.cc.

References analyzeHits(), analyzePassiveHits(), caloGeometry_, MillePedeFileConverter_cfg::e, CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, hcalGeom_, mps_fire::i, edm::HandleBase::isValid(), passive_, tok_passive_, tokGeom_, and toks_calo_.

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

  caloGeometry_ = &set.getData(tokGeom_);
  hcalGeom_ = static_cast<const HcalGeometry*>(caloGeometry_->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));

  for (unsigned int i = 0; i < toks_calo_.size(); i++) {
    const edm::Handle<edm::PCaloHitContainer>& hitsCalo = e.getHandle(toks_calo_[i]);
    bool getHits = (hitsCalo.isValid());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis: Input flags Hits[" << i << "]: " << getHits;
#endif
    if (getHits) {
      std::vector<PCaloHit> caloHits;
      caloHits.insert(caloHits.end(), hitsCalo->begin(), hitsCalo->end());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis: Hit buffer [" << i << "] " << caloHits.size();
#endif
      analyzeHits(caloHits, i);
    }
  }

  if (passive_) {
    const edm::Handle<edm::PassiveHitContainer>& hitsPassive = e.getHandle(tok_passive_);
    bool getHits = (hitsPassive.isValid());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis: Passive: " << getHits;
#endif
    if (getHits) {
      std::vector<PassiveHit> passiveHits;
      passiveHits.insert(passiveHits.end(), hitsPassive->begin(), hitsPassive->end());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis: Passive Hit buffer  " << passiveHits.size();
#endif
      analyzePassiveHits(passiveHits);
    }
  }
}

analyzeHits()

void CaloSimHitAnalysis::analyzeHits ( std::vector< PCaloHit > &  hits, int  indx  )

protected

Definition at line 318 of file CaloSimHitAnalysis.cc.

References funct::abs(), caloGeometry_, TauDecayModes::dec, hcalRecHitTable_cff::depth, dqmdumpme::first, CaloGeometry::getPosition(), HcalGeometry::getPosition(), h_edep1_, h_edep_, h_edepEM_, h_edepHad_, h_edepT1_, h_edepT_, h_eta_, h_etaphi_, h_etot_, h_etotg_, h_hit_, h_phi_, h_rr_, h_rz1_, h_rz_, h_time_, h_timeT_, h_zz_, HcalBarrel, HcalEndcap, HcalForward, hcalGeom_, HcalOuter, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, hit::id, heavyIonCSV_trainingSettings::idx, hcalRecHitTable_cff::ieta, dqmdumpme::k, nCalo_, phi, point, DetId::rawId(), edm::second(), HcalDetId::subdet(), tCut_, testNumber_, hcalRecHitTable_cff::time, timeSliceUnit_, tMax_, tScale_, HcalTestNumbering::unpackHcalIndex(), and z.

Referenced by analyze().

                                                                        {
  int nHit = hits.size();
  int nHB = 0, nHE = 0, nHO = 0, nHF = 0, nEB = 0, nEE = 0, nBad = 0;
#ifdef EDM_ML_DEBUG
  int iHit = 0;
#endif

  std::map<CaloHitID, std::pair<double, double> > hitMap;
  double etot[nCalo_], etotG[nCalo_];
  for (unsigned int k = 0; k < nCalo_; ++k) {
    etot[k] = etotG[k] = 0;
  }
  for (const auto& hit : hits) {
    double edep = hit.energy();
    double time = tScale_ * hit.time();
    uint32_t id = hit.id();
    int itra = hit.geantTrackId();
    double edepEM = hit.energyEM();
    double edepHad = hit.energyHad();
    int idx(-1);
    if (indx != 2) {
      idx = indx;
      if (indx == 0)
        ++nEB;
      else
        ++nEE;
    } else {
      int subdet(0);
      if (testNumber_) {
        int ieta(0), phi(0), z(0), lay(0), depth(0);
        HcalTestNumbering::unpackHcalIndex(id, subdet, z, depth, ieta, phi, lay);
        id = HcalDetId(static_cast<HcalSubdetector>(subdet), z * ieta, phi, depth).rawId();
      } else {
        subdet = HcalDetId(id).subdet();
      }
      if (subdet == static_cast<int>(HcalBarrel)) {
        idx = indx;
        nHB++;
      } else if (subdet == static_cast<int>(HcalEndcap)) {
        idx = indx + 1;
        nHE++;
      } else if (subdet == static_cast<int>(HcalOuter)) {
        idx = indx + 2;
        nHO++;
      } else if (subdet == static_cast<int>(HcalForward)) {
        idx = indx + 3;
        nHF++;
      }
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "Hit[" << iHit << ":" << nHit << ":" << idx << "] E " << edep << ":" << edepEM
                                 << ":" << edepHad << " T " << time << " itra " << itra << " ID " << std::hex << id
                                 << std::dec;
    ++iHit;
#endif
    if (idx >= 0) {
      CaloHitID hid(id, time, itra, 0, timeSliceUnit_[indx]);
      auto itr = hitMap.find(hid);
      if (itr == hitMap.end())
        hitMap[hid] = std::make_pair(time, edep);
      else
        ((itr->second).second) += edep;
      h_edepT_[idx]->Fill(edep);
      h_timeT_[idx]->Fill(time);
      if (edepEM > 0)
        h_edepEM_[idx]->Fill(edepEM);
      if (edepHad > 0)
        h_edepHad_[idx]->Fill(edepHad);
      if (time > tCut_)
        h_edepT1_[idx]->Fill(edep);
    } else {
      ++nBad;
    }
  }

  //Now make plots
  for (auto itr = hitMap.begin(); itr != hitMap.end(); ++itr) {
    int idx = -1;
    GlobalPoint point;
    DetId id((itr->first).unitID());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "Index " << indx << " Geom " << caloGeometry_ << ":" << hcalGeom_ << "  "
                                 << std::hex << id.rawId() << std::dec;
#endif
    if (indx != 2) {
      idx = indx;
      point = caloGeometry_->getPosition(id);
    } else {
      int subdet = id.subdetId();
      if (subdet == static_cast<int>(HcalBarrel)) {
        idx = indx;
      } else if (subdet == static_cast<int>(HcalEndcap)) {
        idx = indx + 1;
      } else if (subdet == static_cast<int>(HcalOuter)) {
        idx = indx + 2;
      } else if (subdet == static_cast<int>(HcalForward)) {
        idx = indx + 3;
      }
      point = hcalGeom_->getPosition(id);
    }
    double edep = (itr->second).second;
    double time = (itr->second).first;
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "Index " << idx << ":" << nCalo_ << " Point " << point << " E " << edep << " T "
                                 << time;
#endif
    if (idx >= 0) {
      h_time_[idx]->Fill(time);
      h_edep_[idx]->Fill(edep);
      h_rr_[idx]->Fill(point.perp());
      h_zz_[idx]->Fill(point.z());
      h_eta_[idx]->Fill(point.eta());
      h_phi_[idx]->Fill(point.phi());
      h_rz_->Fill(std::abs(point.z()), point.perp());
      h_etaphi_->Fill(std::abs(point.eta()), std::abs(point.phi()));
      etot[idx] += edep;
      if (time < tMax_)
        etotG[idx] += edep;
      if (time > tCut_) {
        h_edep1_[idx]->Fill(edep);
        h_rz1_->Fill(std::abs(point.z()), point.perp());
      }
    }
  }

  if (indx < 2) {
    h_etot_[indx]->Fill(etot[indx]);
    h_etotg_[indx]->Fill(etotG[indx]);
    if (indx == 0)
      h_hit_[indx]->Fill(double(nEB));
    else
      h_hit_[indx]->Fill(double(nEE));
  } else {
    h_hit_[2]->Fill(double(nHB));
    h_hit_[3]->Fill(double(nHE));
    h_hit_[4]->Fill(double(nHO));
    h_hit_[5]->Fill(double(nHF));
    for (int idx = 2; idx < nCalo_; idx++) {
      h_etot_[idx]->Fill(etot[idx]);
      h_etotg_[idx]->Fill(etotG[idx]);
    }
  }

  edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis::analyzeHits: EB " << nEB << " EE " << nEE << " HB " << nHB
                               << " HE " << nHE << " HO " << nHO << " HF " << nHF << " Bad " << nBad << " All " << nHit
                               << " Reduced " << hitMap.size();
}

analyzePassiveHits()

void CaloSimHitAnalysis::analyzePassiveHits ( std::vector< PassiveHit > &  hits )

protected

Definition at line 466 of file CaloSimHitAnalysis.cc.

References allSteps_, HcalDetId::depth(), detNames_, etaDepth_, h_edepp_, h_edepTk_, h_hitp_, h_rzH_, h_stepp_, h_timep_, h_timeTk_, h_trackp_, DetId::Hcal, hfClusterShapes_cfi::hits, hit::id, HcalDetId::ietaAbs(), dqmdumpme::k, Skims_PA_cff::name, findQualityFiles::rr, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, DiMuonV_cfg::tracks, hit::x, hit::y, and hit::z.

Referenced by analyze().

                                                                       {
  const std::string active = "Active";
  const std::string sensor = "Sensor";
  std::map<std::pair<std::string, uint32_t>, int> hitx;
  std::map<int, int> tracks;
  unsigned int passive1(0), passive2(0);
  for (auto& hit : hits) {
    std::string name = hit.vname();
    std::pair<std::string, uint32_t> volume = std::make_pair(name, (hit.id() % 1000000));
    auto itr = hitx.find(volume);
    if (itr == hitx.end())
      hitx[volume] = 1;
    else
      ++(itr->second);
    auto ktr = tracks.find(hit.trackId());
    if (ktr == tracks.end())
      tracks[hit.trackId()] = 1;
    else
      ++(ktr->second);
    h_edepp_->Fill(hit.energy());
    h_timep_->Fill(hit.time());
    h_stepp_->Fill(hit.stepLength());
    if ((name.find(active) != std::string::npos) || (name.find(sensor) != std::string::npos)) {
      unsigned idet = detNames_.size();
      for (unsigned int k = 0; k < detNames_.size(); ++k) {
        if (name.find(detNames_[k]) != std::string::npos) {
          idet = k;
          break;
        }
      }
      if (idet < detNames_.size()) {
        h_edepTk_[idet]->Fill(hit.energy());
        h_timeTk_[idet]->Fill(hit.time());
      }
    }

    if ((allSteps_ / 100) % 10 > 0) {
      uint32_t id = hit.id();
      if (DetId(id).det() == DetId::Hcal) {
        HcalDetId hid = HcalDetId(id);
        int indx = (100 * hid.ietaAbs() + hid.depth());
        auto itr = etaDepth_.find(indx);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis::ID: " << hid << " Index " << indx << " Iterator "
                                     << (itr != etaDepth_.end());
#endif
        ++passive1;
        if (itr != etaDepth_.end()) {
          uint32_t ipos = itr->second;
          double rr = std::sqrt(hit.x() * hit.x() + hit.y() * hit.y());
          if (ipos < h_rzH_.size()) {
            h_rzH_[ipos]->Fill(hit.z(), rr);
            ++passive2;
          }
        }
      }
    }
  }
  h_hitp_->Fill(hitx.size());
  h_trackp_->Fill(tracks.size());
  edm::LogVerbatim("HitStudy") << "CaloSimHitAnalysis::analyzPassiveHits: Total " << hits.size() << " Cells "
                               << hitx.size() << " Tracks " << tracks.size() << " Passive " << passive1 << ":"
                               << passive2;
}

beginJob()

void CaloSimHitAnalysis::beginJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 52 of file CaloSimHitAnalysis.cc.

{}

beginRun()

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

inlineoverrideprotected

Definition at line 54 of file CaloSimHitAnalysis.cc.

{}

endRun()

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

inlineoverrideprotected

Definition at line 55 of file CaloSimHitAnalysis.cc.

{}

fillDescriptions()

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

static

Definition at line 259 of file CaloSimHitAnalysis.cc.

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

                                                                                    {
  edm::ParameterSetDescription desc;
  std::vector<std::string> labels = {"EcalHitsEB1", "EcalHitsEE1", "HcalHits1"};
  std::vector<double> times = {1, 1, 1};
  desc.addUntracked<std::string>("moduleLabel", "g4SimHits");
  desc.add<std::vector<std::string> >("hitCollection", labels);
  desc.add<std::vector<double> >("timeSliceUnit", times);
  desc.addUntracked<double>("maxEnergy", 250.0);
  desc.addUntracked<double>("maxTime", 1000.0);
  desc.addUntracked<double>("timeCut", 100.0);
  desc.addUntracked<double>("timeScale", 1.0);
  desc.addUntracked<double>("timeThreshold", 15.0);
  desc.addUntracked<bool>("testNumbering", false);
  desc.addUntracked<bool>("passiveHits", false);
  std::vector<std::string> names = {"PixelBarrel", "PixelForward", "TIB", "TID", "TOB", "TEC"};
  desc.addUntracked<std::vector<std::string> >("detNames", names);
  desc.addUntracked<int>("allStep", 100);
  descriptions.add("caloSimHitAnalysis", desc);
}

Member Data Documentation

allSteps_

const int CaloSimHitAnalysis::allSteps_

private

Definition at line 66 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

caloGeometry_

const CaloGeometry* CaloSimHitAnalysis::caloGeometry_

private

Definition at line 72 of file CaloSimHitAnalysis.cc.

Referenced by analyze(), and analyzeHits().

detNames_

const std::vector<std::string> CaloSimHitAnalysis::detNames_

private

Definition at line 67 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

etaDepth_

std::map<int, unsigned int> CaloSimHitAnalysis::etaDepth_

private

Definition at line 84 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

g4Label_

const std::string CaloSimHitAnalysis::g4Label_

private

Definition at line 61 of file CaloSimHitAnalysis.cc.

Referenced by CaloSimHitAnalysis().

h_edep1_

TH1F * CaloSimHitAnalysis::h_edep1_[nCalo_]

private

Definition at line 77 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_edep_

TH1F * CaloSimHitAnalysis::h_edep_[nCalo_]

private

Definition at line 76 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_edepEM_

TH1F* CaloSimHitAnalysis::h_edepEM_[nCalo_]

private

Definition at line 78 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_edepHad_

TH1F * CaloSimHitAnalysis::h_edepHad_[nCalo_]

private

Definition at line 78 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_edepp_

TH1F * CaloSimHitAnalysis::h_edepp_

private

Definition at line 81 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_edepT1_

TH1F * CaloSimHitAnalysis::h_edepT1_[nCalo_]

private

Definition at line 77 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_edepT_

TH1F * CaloSimHitAnalysis::h_edepT_[nCalo_]

private

Definition at line 76 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_edepTk_

std::vector<TH1F*> CaloSimHitAnalysis::h_edepTk_

private

Definition at line 82 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_eta_

TH1F* CaloSimHitAnalysis::h_eta_[nCalo_]

private

Definition at line 79 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_etaphi_

TH2F * CaloSimHitAnalysis::h_etaphi_

private

Definition at line 80 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_etot_

TH1F * CaloSimHitAnalysis::h_etot_[nCalo_]

private

Definition at line 79 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_etotg_

TH1F * CaloSimHitAnalysis::h_etotg_[nCalo_]

private

Definition at line 79 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_hit_

TH1F* CaloSimHitAnalysis::h_hit_[nCalo_]

private

Definition at line 76 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_hitp_

TH1F* CaloSimHitAnalysis::h_hitp_

private

Definition at line 81 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_phi_

TH1F * CaloSimHitAnalysis::h_phi_[nCalo_]

private

Definition at line 79 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_rr_

TH1F * CaloSimHitAnalysis::h_rr_[nCalo_]

private

Definition at line 78 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_rz1_

TH2F * CaloSimHitAnalysis::h_rz1_

private

Definition at line 80 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_rz_

TH2F* CaloSimHitAnalysis::h_rz_

private

Definition at line 80 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_rzH_

std::vector<TH2F*> CaloSimHitAnalysis::h_rzH_

private

Definition at line 83 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_stepp_

TH1F * CaloSimHitAnalysis::h_stepp_

private

Definition at line 81 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_time_

TH1F * CaloSimHitAnalysis::h_time_[nCalo_]

private

Definition at line 76 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_timep_

TH1F * CaloSimHitAnalysis::h_timep_

private

Definition at line 81 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_timeT_

TH1F* CaloSimHitAnalysis::h_timeT_[nCalo_]

private

Definition at line 77 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

h_timeTk_

std::vector<TH1F*> CaloSimHitAnalysis::h_timeTk_

private

Definition at line 82 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_trackp_

TH1F * CaloSimHitAnalysis::h_trackp_

private

Definition at line 81 of file CaloSimHitAnalysis.cc.

Referenced by analyzePassiveHits().

h_zz_

TH1F * CaloSimHitAnalysis::h_zz_[nCalo_]

private

Definition at line 78 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

hcalGeom_

const HcalGeometry* CaloSimHitAnalysis::hcalGeom_

private

Definition at line 73 of file CaloSimHitAnalysis.cc.

Referenced by analyze(), and analyzeHits().

hitLab_

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

private

Definition at line 62 of file CaloSimHitAnalysis.cc.

Referenced by CaloSimHitAnalysis().

maxEnergy_

const double CaloSimHitAnalysis::maxEnergy_

private

Definition at line 64 of file CaloSimHitAnalysis.cc.

maxTime_

const double CaloSimHitAnalysis::maxTime_

private

Definition at line 64 of file CaloSimHitAnalysis.cc.

nCalo_

constexpr int CaloSimHitAnalysis::nCalo_ = 6

staticprivate

Definition at line 75 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

passive_

const bool CaloSimHitAnalysis::passive_

private

Definition at line 65 of file CaloSimHitAnalysis.cc.

Referenced by analyze().

tCut_

const double CaloSimHitAnalysis::tCut_

private

Definition at line 64 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

testNumber_

const bool CaloSimHitAnalysis::testNumber_

private

Definition at line 65 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

timeSliceUnit_

const std::vector<double> CaloSimHitAnalysis::timeSliceUnit_

private

Definition at line 63 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

tMax_

const double CaloSimHitAnalysis::tMax_

private

Definition at line 64 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().

tok_passive_

const edm::EDGetTokenT<edm::PassiveHitContainer> CaloSimHitAnalysis::tok_passive_

private

Definition at line 70 of file CaloSimHitAnalysis.cc.

Referenced by analyze().

tokGeom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> CaloSimHitAnalysis::tokGeom_

private

Definition at line 68 of file CaloSimHitAnalysis.cc.

Referenced by analyze().

toks_calo_

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

private

Definition at line 69 of file CaloSimHitAnalysis.cc.

Referenced by analyze().

tScale_

const double CaloSimHitAnalysis::tScale_

private

Definition at line 64 of file CaloSimHitAnalysis.cc.

Referenced by analyzeHits().