Inheritance diagram for CaloSimHitStudy:

Public Member Functions
	CaloSimHitStudy (const edm::ParameterSet &ps)

	~CaloSimHitStudy () 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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	analyzeHits (const edm::Handle< edm::PSimHitContainer > &, int)

void	analyzeHits (std::vector< PCaloHit > &, std::vector< PCaloHit > &, std::vector< PCaloHit > &)

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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
TH1F *	edep_ [9]

TH1F *	edepC_ [9]

TH1F *	edepEM_ [9]

TH1F *	edepHad_ [9]

TH1F *	edepT_ [9]

TH1F *	edepTk_ [16]

const double	eMIP_

TH1F *	eneInc_

TH1F *	etaInc_

TH1F *	etot_ [9]

TH1F *	etotg_ [9]

const std::string	g4Label_

TH1F *	hit_ [9]

TH1F *	hitHigh

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

TH1F *	hitLow

TH1F *	hitMu

TH1F *	hitTk_ [16]

const double	maxEnergy_

const std::vector< std::string >	muonLab_ = {"MuonRPCHits", "MuonCSCHits", "MuonDTHits", "MuonGEMHits"}

TH1F *	phiInc_

TH1F *	ptInc_

const bool	storeRL_

const bool	testNumber_

TH1F *	time_ [9]

TH1F *	timeAll_ [9]

const std::vector< std::string >	tkHighLab_

const std::vector< std::string >	tkLowLab_

const double	tmax_

TH1F *	tofTk_ [16]

edm::EDGetTokenT < edm::HepMCProduct >	tok_evt_

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

std::vector< edm::EDGetTokenT < edm::PSimHitContainer > >	toks_tkHigh_

std::vector< edm::EDGetTokenT < edm::PSimHitContainer > >	toks_tkLow_

std::vector< edm::EDGetTokenT < edm::PSimHitContainer > >	toks_track_

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 39 of file CaloSimHitStudy.cc.

Constructor & Destructor Documentation

CaloSimHitStudy::CaloSimHitStudy ( const edm::ParameterSet & ps )

Definition at line 86 of file CaloSimHitStudy.cc.

References edep_, edepC_, edepEM_, edepHad_, edepT_, edepTk_, eMIP_, eneInc_, etaInc_, etot_, etotg_, g4Label_, edm::ParameterSet::getUntrackedParameter(), hit_, hitHigh, hitLab_, hitLow, hitMu, hitTk_, mps_fire::i, HLT_FULL_cff::InputTag, edm::Service< T >::isAvailable(), TFileService::kSharedResource, TFileService::make(), maxEnergy_, muonLab_, mergeVDriftHistosByStation::name, phiInc_, ptInc_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, time_, timeAll_, runGCPTkAlMap::title, tkHighLab_, tkLowLab_, tmax_, tofTk_, tok_evt_, toks_calo_, toks_tkHigh_, toks_tkLow_, toks_track_, and SiStrip_OfflineMonitoring_cff::ymax.

     : g4Label_(ps.getUntrackedParameter<std::string>("ModuleLabel")),
       hitLab_(ps.getUntrackedParameter<std::vector<std::string>>("CaloCollection")),
       maxEnergy_(ps.getUntrackedParameter<double>("MaxEnergy", 200.0)),
       tmax_(ps.getUntrackedParameter<double>("TimeCut", 100.0)),
       eMIP_(ps.getUntrackedParameter<double>("MIPCut", 0.70)),
       storeRL_(ps.getUntrackedParameter<bool>("StoreRL", false)),
       testNumber_(ps.getUntrackedParameter<bool>("TestNumbering", true)) {
   usesResource(TFileService::kSharedResource);
 
   tok_evt_ =
       consumes<edm::HepMCProduct>(edm::InputTag(ps.getUntrackedParameter<std::string>("SourceLabel", "VtxSmeared")));
   for (unsigned i = 0; i != hitLab_.size(); i++)
     toks_calo_.emplace_back(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hitLab_[i])));
   for (unsigned i = 0; i != muonLab_.size(); i++)
     toks_track_.emplace_back(consumes<edm::PSimHitContainer>(edm::InputTag(g4Label_, muonLab_[i])));
   for (unsigned i = 0; i != tkHighLab_.size(); i++)
     toks_tkHigh_.emplace_back(consumes<edm::PSimHitContainer>(edm::InputTag(g4Label_, tkHighLab_[i])));
   for (unsigned i = 0; i != tkLowLab_.size(); i++)
     toks_tkLow_.emplace_back(consumes<edm::PSimHitContainer>(edm::InputTag(g4Label_, tkLowLab_[i])));
 
   edm::LogVerbatim("HitStudy") << "Module Label: " << g4Label_ << "   Hits: " << hitLab_[0] << ", " << hitLab_[1]
                                << ", " << hitLab_[2] << ", " << hitLab_[3] << "   MaxEnergy: " << maxEnergy_
                                << "  Tmax: " << tmax_ << "  MIP Cut: " << eMIP_;
 
   edm::Service<TFileService> tfile;
 
   if (!tfile.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
   char name[20], title[120];
   sprintf(title, "Incident PT (GeV)");
   ptInc_ = tfile->make<TH1F>("PtInc", title, 1000, 0., maxEnergy_);
   ptInc_->GetXaxis()->SetTitle(title);
   ptInc_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Incident Energy (GeV)");
   eneInc_ = tfile->make<TH1F>("EneInc", title, 1000, 0., maxEnergy_);
   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");
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Completed defining histos for incident particle";
 #endif
   std::string dets[9] = {"EB", "EB(APD)", "EB(ATJ)", "EE", "ES", "HB", "HE", "HO", "HF"};
   double nhcMax[9] = {40000., 2000., 2000., 40000., 10000., 10000., 10000., 2000., 10000.};
   for (int i = 0; i < 9; 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., nhcMax[i]);
     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, 1000, 0., 1000.);
     time_[i]->GetXaxis()->SetTitle(title);
     time_[i]->GetYaxis()->SetTitle("Hits");
     sprintf(name, "TimeAll%d", i);
     sprintf(title, "Hit time (ns) in %s with no check on Track ID", dets[i].c_str());
     timeAll_[i] = tfile->make<TH1F>(name, title, 1000, 0., 1000.);
     timeAll_[i]->GetXaxis()->SetTitle(title);
     timeAll_[i]->GetYaxis()->SetTitle("Hits");
     double ymax = maxEnergy_;
     if (i == 1 || i == 2 || i == 4)
       ymax = 1.0;
     else if (i > 4 && i < 8)
       ymax = 10.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, "EdepEM%d", i);
     sprintf(title, "Energy deposit (GeV) by EM particles in %s", dets[i].c_str());
     edepEM_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepEM_[i]->GetXaxis()->SetTitle(title);
     edepEM_[i]->GetYaxis()->SetTitle("Hits");
     sprintf(name, "EdepHad%d", i);
     sprintf(title, "Energy deposit (GeV) by hadrons in %s", dets[i].c_str());
     edepHad_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepHad_[i]->GetXaxis()->SetTitle(title);
     edepHad_[i]->GetYaxis()->SetTitle("Hits");
     sprintf(name, "Etot%d", i);
     sprintf(title, "Total energy deposit (GeV) in %s", dets[i].c_str());
     etot_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     etot_[i]->GetXaxis()->SetTitle(title);
     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());
     etotg_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     etotg_[i]->GetXaxis()->SetTitle(title);
     etotg_[i]->GetYaxis()->SetTitle("Events");
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Completed defining histos for first level of Calorimeter";
 #endif
   std::string detx[9] = {"EB/EE (MIP)",
                          "HB/HE (MIP)",
                          "HB/HE/HO (MIP)",
                          "EB/EE (no MIP)",
                          "HB/HE (no MIP)",
                          "HB/HE/HO (no MIP)",
                          "EB/EE (All)",
                          "HB/HE (All)",
                          "HB/HE/HO (All)"};
   for (int i = 0; i < 9; i++) {
     double ymax = 1.0;
     if (i == 0 || i == 3 || i == 6)
       ymax = maxEnergy_;
     sprintf(name, "EdepCal%d", i);
     sprintf(title, "Energy deposit in %s", detx[i].c_str());
     edepC_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepC_[i]->GetXaxis()->SetTitle(title);
     edepC_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EdepCalT%d", i);
     sprintf(title, "Energy deposit (t < %f ns) in %s", tmax_, detx[i].c_str());
     edepT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepT_[i]->GetXaxis()->SetTitle(title);
     edepT_[i]->GetYaxis()->SetTitle("Events");
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Completed defining histos for second level of Calorimeter";
 #endif
   hitLow = tfile->make<TH1F>("HitLow", "Number of hits in Track (Low)", 1000, 0, 20000.);
   hitLow->GetXaxis()->SetTitle("Number of hits in Track (Low)");
   hitLow->GetYaxis()->SetTitle("Events");
   hitHigh = tfile->make<TH1F>("HitHigh", "Number of hits in Track (High)", 1000, 0, 5000.);
   hitHigh->GetXaxis()->SetTitle("Number of hits in Track (High)");
   hitHigh->GetYaxis()->SetTitle("Events");
   hitMu = tfile->make<TH1F>("HitMu", "Number of hits in Track (Muon)", 1000, 0, 2000.);
   hitMu->GetXaxis()->SetTitle("Number of hits in Muon");
   hitMu->GetYaxis()->SetTitle("Events");
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Completed defining histos for general tracking hits";
 #endif
   std::string dett[16] = {"Pixel Barrel (High)",
                           "Pixel Endcap (High)",
                           "TEC (High)",
                           "TIB (High)",
                           "TID (High)",
                           "TOB (High)",
                           "Pixel Barrel (Low)",
                           "Pixel Endcap (Low)",
                           "TEC (Low)",
                           "TIB (Low)",
                           "TID (Low)",
                           "TOB (Low)",
                           "RPC",
                           "CSC",
                           "DT",
                           "GEM"};
   double nhtMax[16] = {
       500., 500., 1000., 1000., 500., 1000., 5000., 2000., 10000., 5000., 2000., 5000., 500., 1000., 1000., 500.};
   for (int i = 0; i < 16; i++) {
     sprintf(name, "HitTk%d", i);
     sprintf(title, "Number of hits in %s", dett[i].c_str());
     hitTk_[i] = tfile->make<TH1F>(name, title, 1000, 0., nhtMax[i]);
     hitTk_[i]->GetXaxis()->SetTitle(title);
     hitTk_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "TimeTk%d", i);
     sprintf(title, "Time of the hit (ns) in %s", dett[i].c_str());
     tofTk_[i] = tfile->make<TH1F>(name, title, 1000, 0., 200.);
     tofTk_[i]->GetXaxis()->SetTitle(title);
     tofTk_[i]->GetYaxis()->SetTitle("Hits");
     sprintf(name, "EdepTk%d", i);
     sprintf(title, "Energy deposit (GeV) in %s", dett[i].c_str());
     double ymax = (i < 12) ? 0.25 : 0.005;
     edepTk_[i] = tfile->make<TH1F>(name, title, 250, 0., ymax);
     edepTk_[i]->GetXaxis()->SetTitle(title);
     edepTk_[i]->GetYaxis()->SetTitle("Hits");
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Completed defining histos for SimHit objects";
 #endif
 }

CaloSimHitStudy::~CaloSimHitStudy ( )

inlineoverride

Definition at line 42 of file CaloSimHitStudy.cc.

42 {}

Member Function Documentation

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 282 of file CaloSimHitStudy.cc.

References analyzeHits(), eneInc_, etaInc_, edm::EventID::event(), edm::Event::getHandle(), hitHigh, hitLow, hitMu, mps_fire::i, edm::EventBase::id(), edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, phiInc_, ptInc_, edm::EventID::run(), tok_evt_, toks_calo_, toks_tkHigh_, toks_tkLow_, and toks_track_.

                                                                      {
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy:Run = " << e.id().run() << " Event = " << e.id().event();
 
   const edm::Handle<edm::HepMCProduct> EvtHandle = e.getHandle(tok_evt_);
   const HepMC::GenEvent* myGenEvent = EvtHandle->GetEvent();
 
   double eInc = 0, etaInc = 0, 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();
   }
   double ptInc = eInc / std::cosh(etaInc);
   ptInc_->Fill(ptInc);
   eneInc_->Fill(eInc);
   etaInc_->Fill(etaInc);
   phiInc_->Fill(phiInc);
 
   std::vector<PCaloHit> ebHits, eeHits, hcHits;
   for (unsigned int i = 0; i < toks_calo_.size(); i++) {
     bool getHits = false;
     const edm::Handle<edm::PCaloHitContainer>& hitsCalo = e.getHandle(toks_calo_[i]);
     if (hitsCalo.isValid())
       getHits = true;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Input flags Hits " << getHits;
 #endif
     if (getHits) {
       std::vector<PCaloHit> caloHits;
       caloHits.insert(caloHits.end(), hitsCalo->begin(), hitsCalo->end());
       if (i == 0)
         ebHits.insert(ebHits.end(), hitsCalo->begin(), hitsCalo->end());
       else if (i == 1)
         eeHits.insert(eeHits.end(), hitsCalo->begin(), hitsCalo->end());
       else if (i == 3)
         hcHits.insert(hcHits.end(), hitsCalo->begin(), hitsCalo->end());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HitStudy") << "CaloSimHitStudy: Hit buffer " << caloHits.size();
 #endif
       analyzeHits(caloHits, i);
     }
   }
   analyzeHits(ebHits, eeHits, hcHits);
 
   std::vector<PSimHit> muonHits;
   for (unsigned int i = 0; i < toks_track_.size(); i++) {
     const edm::Handle<edm::PSimHitContainer>& hitsTrack = e.getHandle(toks_track_[i]);
     if (hitsTrack.isValid()) {
       muonHits.insert(muonHits.end(), hitsTrack->begin(), hitsTrack->end());
       analyzeHits(hitsTrack, i + 12);
     }
   }
   unsigned int nhmu = muonHits.size();
   hitMu->Fill(double(nhmu));
   std::vector<PSimHit> tkHighHits;
   for (unsigned int i = 0; i < toks_tkHigh_.size(); i++) {
     const edm::Handle<edm::PSimHitContainer>& hitsTrack = e.getHandle(toks_tkHigh_[i]);
     if (hitsTrack.isValid()) {
       tkHighHits.insert(tkHighHits.end(), hitsTrack->begin(), hitsTrack->end());
       analyzeHits(hitsTrack, i);
     }
   }
   unsigned int nhtkh = tkHighHits.size();
   hitHigh->Fill(double(nhtkh));
   std::vector<PSimHit> tkLowHits;
   for (unsigned int i = 0; i < toks_tkLow_.size(); i++) {
     const edm::Handle<edm::PSimHitContainer>& hitsTrack = e.getHandle(toks_tkLow_[i]);
     if (hitsTrack.isValid()) {
       tkLowHits.insert(tkLowHits.end(), hitsTrack->begin(), hitsTrack->end());
       analyzeHits(hitsTrack, i + 6);
     }
   }
   unsigned int nhtkl = tkLowHits.size();
   hitLow->Fill(double(nhtkl));
 }

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

protected

Definition at line 359 of file CaloSimHitStudy.cc.

References HLT_FULL_cff::depth, edep_, edepEM_, edepHad_, etot_, etotg_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hit_, mps_fire::i, PCaloHit::kEcalDepthIdMask, phi, storeRL_, testNumber_, time_, timeAll_, HcalTestNumbering::unpackHcalIndex(), and z.

Referenced by analyze().

                                                                      {
   int nHit = hits.size();
   int nHB(0), nHE(0), nHO(0), nHF(0), nEB(0), nEBAPD(0), nEBATJ(0), nEE(0), nES(0);
 #ifdef EDM_ML_DEBUG
   int nBad(0);
 #endif
   std::map<unsigned int, double> hitMap;
   std::vector<double> etot(9, 0), etotG(9, 0);
   for (int i = 0; i < nHit; i++) {
     double edep = hits[i].energy();
     double time = hits[i].time();
     unsigned int id_ = hits[i].id();
     double edepEM = hits[i].energyEM();
     double edepHad = hits[i].energyHad();
     if (indx == 0) {
       int dep = (hits[i].depth()) & PCaloHit::kEcalDepthIdMask;
       if (dep == 1)
         id_ |= 0x20000;
       else if (dep == 2)
         id_ |= 0x40000;
     }
     std::map<unsigned int, double>::const_iterator it = hitMap.find(id_);
     if (it == hitMap.end()) {
       hitMap.insert(std::pair<unsigned int, double>(id_, time));
     }
     int idx = -1;
     if (indx != 3) {
       if (indx == 0) {
         if (storeRL_)
           idx = 0;
         else
           idx = ((hits[i].depth()) & PCaloHit::kEcalDepthIdMask);
       } else
         idx = indx + 2;
       time_[idx]->Fill(time);
       edep_[idx]->Fill(edep);
       edepEM_[idx]->Fill(edepEM);
       edepHad_[idx]->Fill(edepHad);
       if (idx == 0)
         nEB++;
       else if (idx == 1)
         nEBAPD++;
       else if (idx == 2)
         nEBATJ++;
       else if (idx == 3)
         nEE++;
       else if (idx == 4)
         nES++;
 #ifdef EDM_ML_DEBUG
       else
         nBad++;
 #endif
       if (indx >= 0 && indx < 3) {
         etot[idx] += edep;
         if (time < 100)
           etotG[idx] += edep;
       }
     } 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);
       } else {
         HcalDetId id = HcalDetId(id_);
         subdet = id.subdet();
       }
       if (subdet == static_cast<int>(HcalBarrel)) {
         idx = indx + 2;
         nHB++;
       } else if (subdet == static_cast<int>(HcalEndcap)) {
         idx = indx + 3;
         nHE++;
       } else if (subdet == static_cast<int>(HcalOuter)) {
         idx = indx + 4;
         nHO++;
       } else if (subdet == static_cast<int>(HcalForward)) {
         idx = indx + 5;
         nHF++;
       }
       if (idx > 0) {
         time_[idx]->Fill(time);
         edep_[idx]->Fill(edep);
         edepEM_[idx]->Fill(edepEM);
         edepHad_[idx]->Fill(edepHad);
         etot[idx] += edep;
         if (time < 100)
           etotG[idx] += edep;
       } else {
 #ifdef EDM_ML_DEBUG
         nBad++;
 #endif
       }
     }
   }
   if (indx < 3) {
     etot_[indx + 2]->Fill(etot[indx + 2]);
     etotg_[indx + 2]->Fill(etotG[indx + 2]);
     if (indx == 0) {
       etot_[indx]->Fill(etot[indx]);
       etotg_[indx]->Fill(etotG[indx]);
       etot_[indx + 1]->Fill(etot[indx + 1]);
       etotg_[indx + 1]->Fill(etotG[indx + 1]);
       hit_[indx]->Fill(double(nEB));
       hit_[indx + 1]->Fill(double(nEBAPD));
       hit_[indx + 2]->Fill(double(nEBATJ));
     } else {
       hit_[indx + 2]->Fill(double(nHit));
     }
   } else if (indx == 3) {
     hit_[5]->Fill(double(nHB));
     hit_[6]->Fill(double(nHE));
     hit_[7]->Fill(double(nHO));
     hit_[8]->Fill(double(nHF));
     for (int idx = 5; idx < 9; idx++) {
       etot_[idx]->Fill(etot[idx]);
       etotg_[idx]->Fill(etotG[idx]);
     }
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy::analyzeHits: EB " << nEB << ", " << nEBAPD << ", " << nEBATJ
                                << " EE " << nEE << " ES " << nES << " HB " << nHB << " HE " << nHE << " HO " << nHO
                                << " HF " << nHF << " Bad " << nBad << " All " << nHit << " Reduced " << hitMap.size();
 #endif
   std::map<unsigned int, double>::const_iterator it = hitMap.begin();
   for (; it != hitMap.end(); it++) {
     double time = it->second;
     unsigned int id_ = (it->first);
     int idx = -1;
     if (indx < 3) {
       if (indx == 0) {
         if ((id_ & 0x20000) != 0)
           idx = indx + 1;
         else if ((id_ & 0x40000) != 0)
           idx = indx + 1;
         else
           idx = indx;
       } else {
         idx = indx + 2;
       }
       if (idx >= 0 && idx < 5)
         timeAll_[idx]->Fill(time);
     } else if (indx == 3) {
       int idx(-1), subdet(0);
       if (testNumber_) {
         int ieta(0), phi(0), z(0), lay(0), depth(0);
         HcalTestNumbering::unpackHcalIndex(id_, subdet, z, depth, ieta, phi, lay);
       } else {
         HcalDetId id = HcalDetId(id_);
         subdet = id.subdet();
       }
       if (subdet == static_cast<int>(HcalBarrel)) {
         idx = indx + 2;
       } else if (subdet == static_cast<int>(HcalEndcap)) {
         idx = indx + 3;
       } else if (subdet == static_cast<int>(HcalOuter)) {
         idx = indx + 4;
       } else if (subdet == static_cast<int>(HcalForward)) {
         idx = indx + 5;
       }
       if (idx > 0) {
         timeAll_[idx]->Fill(time);
       }
     }
   }
 }

void CaloSimHitStudy::analyzeHits	(	const edm::Handle< edm::PSimHitContainer > &	hits,
		int	indx
	)

protected

Definition at line 526 of file CaloSimHitStudy.cc.

References edepTk_, hitTk_, label, muonLab_, AlCaHLTBitMon_QueryRunRegistry::string, tkHighLab_, tkLowLab_, and tofTk_.

                                                                                       {
   int nHit = 0;
   edm::PSimHitContainer::const_iterator ihit;
   std::string label(" ");
   if (indx >= 0 && indx < 6)
     label = tkHighLab_[indx];
   else if (indx >= 6 && indx < 12)
     label = tkLowLab_[indx - 6];
   else if (indx >= 12 && indx < 16)
     label = muonLab_[indx - 12];
   for (ihit = hits->begin(); ihit != hits->end(); ihit++) {
     edepTk_[indx]->Fill(ihit->energyLoss());
     tofTk_[indx]->Fill(ihit->timeOfFlight());
     nHit++;
   }
   hitTk_[indx]->Fill(float(nHit));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy::analyzeHits: for " << label << " Index " << indx << " # of Hits "
                                << nHit;
 #endif
 }

void CaloSimHitStudy::analyzeHits	(	std::vector< PCaloHit > &	ebHits,
		std::vector< PCaloHit > &	eeHits,
		std::vector< PCaloHit > &	hcHits
	)

protected

Definition at line 548 of file CaloSimHitStudy.cc.

References HLT_FULL_cff::depth, edepC_, edepT_, eMIP_, HcalBarrel, HcalEndcap, HcalOuter, mps_fire::i, PCaloHit::kEcalDepthIdMask, phi, testNumber_, tmax_, HcalTestNumbering::unpackHcalIndex(), and z.

                                                                {
   double edepEB = 0, edepEBT = 0;
   for (unsigned int i = 0; i < ebHits.size(); i++) {
     double edep = ebHits[i].energy();
     double time = ebHits[i].time();
     if (((ebHits[i].depth()) & PCaloHit::kEcalDepthIdMask) == 0) {
       edepEB += edep;
       if (time < tmax_)
         edepEBT += edep;
     }
   }
   double edepEE = 0, edepEET = 0;
   for (unsigned int i = 0; i < eeHits.size(); i++) {
     double edep = eeHits[i].energy();
     double time = eeHits[i].time();
     edepEE += edep;
     if (time < tmax_)
       edepEET += edep;
   }
   double edepH = 0, edepHT = 0, edepHO = 0, edepHOT = 0;
   for (unsigned int i = 0; i < hcHits.size(); i++) {
     double edep = hcHits[i].energy();
     double time = hcHits[i].time();
     int subdet(0);
     if (testNumber_) {
       int ieta(0), phi(0), z(0), lay(0), depth(0);
       HcalTestNumbering::unpackHcalIndex(hcHits[i].id(), subdet, z, depth, ieta, phi, lay);
     } else {
       HcalDetId id = HcalDetId(hcHits[i].id());
       subdet = id.subdet();
     }
     if (subdet == static_cast<int>(HcalBarrel) || subdet == static_cast<int>(HcalEndcap)) {
       edepH += edep;
       if (time < tmax_)
         edepHT += edep;
     } else if (subdet == static_cast<int>(HcalOuter)) {
       edepHO += edep;
       if (time < tmax_)
         edepHOT += edep;
     }
   }
   double edepE = edepEB + edepEE;
   double edepET = edepEBT + edepEET;
   double edepHC = edepH + edepHO;
   double edepHCT = edepHT + edepHOT;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "CaloSimHitStudy::energy in EB " << edepEB << " (" << edepEBT << ") from "
                                << ebHits.size() << " hits; "
                                << " energy in EE " << edepEE << " (" << edepEET << ") from " << eeHits.size()
                                << " hits; energy in HC " << edepH << ", " << edepHO << " (" << edepHT << ", " << edepHOT
                                << ") from " << hcHits.size() << " hits";
 #endif
   edepC_[6]->Fill(edepE);
   edepT_[6]->Fill(edepET);
   edepC_[7]->Fill(edepH);
   edepT_[7]->Fill(edepHT);
   edepC_[8]->Fill(edepHC);
   edepT_[8]->Fill(edepHCT);
   if (edepE < eMIP_) {
     edepC_[0]->Fill(edepE);
     edepC_[1]->Fill(edepH);
     edepC_[2]->Fill(edepHC);
   } else {
     edepC_[3]->Fill(edepE);
     edepC_[4]->Fill(edepH);
     edepC_[5]->Fill(edepHC);
   }
   if (edepET < eMIP_) {
     edepT_[0]->Fill(edepET);
     edepT_[1]->Fill(edepHT);
     edepT_[2]->Fill(edepHCT);
   } else {
     edepT_[3]->Fill(edepET);
     edepT_[4]->Fill(edepHT);
     edepT_[5]->Fill(edepHCT);
   }
 }

void CaloSimHitStudy::beginJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 46 of file CaloSimHitStudy.cc.

46 {}

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

inlineoverrideprotected

Definition at line 48 of file CaloSimHitStudy.cc.

48 {}

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

inlineoverrideprotected

Definition at line 49 of file CaloSimHitStudy.cc.

49 {}

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

static

Definition at line 268 of file CaloSimHitStudy.cc.

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

                                                                                  {
   edm::ParameterSetDescription desc;
   std::vector<std::string> calonames = {"EcalHitsEB", "EcalHitsEE", "EcalHitsES", "HcalHits"};
   desc.addUntracked<std::string>("SourceLabel", "generatorSmeared");
   desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
   desc.addUntracked<std::vector<std::string>>("CaloCollection", calonames);
   desc.addUntracked<double>("MaxEnergy", 200.0);
   desc.addUntracked<double>("TimeCut", 100.0);
   desc.addUntracked<double>("MIPCut", 0.70);
   desc.addUntracked<bool>("StoreRL", false);
   desc.addUntracked<bool>("TestNumbering", true);
   descriptions.add("CaloSimHitStudy", desc);
 }

Member Data Documentation

TH1F * CaloSimHitStudy::edep_[9]

private

Definition at line 80 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::edepC_[9]

private

Definition at line 83 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::edepEM_[9]

private

Definition at line 80 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::edepHad_[9]

private

Definition at line 80 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::edepT_[9]

private

Definition at line 83 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::edepTk_[16]

private

Definition at line 83 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

const double CaloSimHitStudy::eMIP_

private

Definition at line 58 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::eneInc_

private

Definition at line 82 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::etaInc_

private

Definition at line 82 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

TH1F* CaloSimHitStudy::etot_[9]

private

Definition at line 81 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::etotg_[9]

private

Definition at line 81 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

const std::string CaloSimHitStudy::g4Label_

private

Definition at line 56 of file CaloSimHitStudy.cc.

Referenced by CaloSimHitStudy().

TH1F* CaloSimHitStudy::hit_[9]

private

Definition at line 80 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::hitHigh

private

Definition at line 81 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

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

private

Definition at line 57 of file CaloSimHitStudy.cc.

Referenced by CaloSimHitStudy().

TH1F* CaloSimHitStudy::hitLow

private

Definition at line 82 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::hitMu

private

Definition at line 81 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

TH1F* CaloSimHitStudy::hitTk_[16]

private

Definition at line 83 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

const double CaloSimHitStudy::maxEnergy_

private

Definition at line 58 of file CaloSimHitStudy.cc.

Referenced by CaloSimHitStudy().

const std::vector<std::string> CaloSimHitStudy::muonLab_ = {"MuonRPCHits", "MuonCSCHits", "MuonDTHits", "MuonGEMHits"}

private

Definition at line 66 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::phiInc_

private

Definition at line 82 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::ptInc_

private

Definition at line 82 of file CaloSimHitStudy.cc.

Referenced by analyze(), and CaloSimHitStudy().

const bool CaloSimHitStudy::storeRL_

private

Definition at line 59 of file CaloSimHitStudy.cc.

Referenced by analyzeHits().

const bool CaloSimHitStudy::testNumber_

private

Definition at line 59 of file CaloSimHitStudy.cc.

Referenced by analyzeHits().

TH1F * CaloSimHitStudy::time_[9]

private

Definition at line 80 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

TH1F * CaloSimHitStudy::timeAll_[9]

private

Definition at line 81 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

const std::vector<std::string> CaloSimHitStudy::tkHighLab_

private

Initial value:

= {"TrackerHitsPixelBarrelHighTof",
                                               "TrackerHitsPixelEndcapHighTof",
                                               "TrackerHitsTECHighTof",
                                               "TrackerHitsTIBHighTof",
                                               "TrackerHitsTIDHighTof",
                                               "TrackerHitsTOBHighTof"}

Definition at line 67 of file CaloSimHitStudy.cc.

Referenced by analyzeHits(), and CaloSimHitStudy().

const std::vector<std::string> CaloSimHitStudy::tkLowLab_

private

Initial value:

= {"TrackerHitsPixelBarrelLowTof",
                                              "TrackerHitsPixelEndcapLowTof",
                                              "TrackerHitsTECLowTof",
                                              "TrackerHitsTIBLowTof",
                                              "TrackerHitsTIDLowTof",
                                              "TrackerHitsTOBLowTof"}