Inheritance diagram for EcalSimHitStudy:

Classes
struct	EcalHit

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

	~EcalSimHitStudy () 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	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
TH1F *	edep_ [ndets_]

TH1F *	edepAll_ [ndets_]

TH1F *	edepEM_ [ndets_]

TH1F *	edepHad_ [ndets_]

TH1F *	eneInc_

TH1F *	etaInc_

TH1F *	etot_ [ndets_]

TH1F *	etotg_ [ndets_]

const std::string	g4Label_

const CaloGeometry *	geometry_

TH1F *	hit_ [ndets_]

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

const double	maxEnergy_

TH1F *	phiInc_

TH2F *	poszn_ [ndets_]

TH2F *	poszp_ [ndets_]

TH1F *	ptInc_

TH1F *	r1by25_ [ndets_]

TH1F *	r1by9_ [ndets_]

TH1F *	r9by25_ [ndets_]

const int	selX_

TH1F *	sEtaEta_ [ndets_]

TH1F *	sEtaPhi_ [ndets_]

TH1F *	sPhiPhi_ [ndets_]

TH1F *	time_ [ndets_]

TH1F *	timeAll_ [ndets_]

const double	tmax_

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

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tokGeom_

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

const double	w0_

Static Private Attributes
static const int	ndets_ = 2

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 46 of file EcalSimHitStudy.cc.

Constructor & Destructor Documentation

◆ EcalSimHitStudy()

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

Definition at line 83 of file EcalSimHitStudy.cc.

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

     : tokGeom_(esConsumes<CaloGeometry, CaloGeometryRecord>()),
       g4Label_(ps.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits")),
       hitLab_(ps.getUntrackedParameter<std::vector<std::string>>("CaloCollection")),
       maxEnergy_(ps.getUntrackedParameter<double>("MaxEnergy", 200.0)),
       tmax_(ps.getUntrackedParameter<double>("TimeCut", 100.0)),
       w0_(ps.getUntrackedParameter<double>("W0", 4.7)),
       selX_(ps.getUntrackedParameter<int>("SelectX", -1)),
       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") << "Module Label: " << g4Label_ << "   Hits: " << hitLab_[0] << ", " << hitLab_[1]
                                << "   MaxEnergy: " << maxEnergy_ << "  Tmax: " << tmax_ << " Select " << selX_;
 }

◆ ~EcalSimHitStudy()

EcalSimHitStudy::~EcalSimHitStudy ( )

inlineoverride

Definition at line 49 of file EcalSimHitStudy.cc.

49 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 256 of file EcalSimHitStudy.cc.

References funct::abs(), analyzeHits(), MillePedeFileConverter_cfg::e, eneInc_, etaInc_, geometry_, edm::EventSetup::getData(), edm::HepMCProduct::GetEvent(), edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, phiInc_, ptInc_, tok_evt_, tokGeom_, and toks_calo_.

                                                                         {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "Run = " << e.id().run() << " Event = " << e.id().event();
 #endif
   // get handles to calogeometry
   geometry_ = &iS.getData(tokGeom_);
 
   double eInc = 0, etaInc = 0, phiInc = 0;
   int type(-1);
   edm::Handle<edm::HepMCProduct> EvtHandle;
   e.getByToken(tok_evt_, EvtHandle);
   if (EvtHandle.isValid()) {
     const HepMC::GenEvent* myGenEvent = EvtHandle->GetEvent();
 
     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);
 
     if (std::abs(etaInc) < 1.46)
       type = 0;
     else if (std::abs(etaInc) > 1.49 && std::abs(etaInc) < 3.0)
       type = 1;
   }
 
   int typeMin = (type < 0) ? 0 : type;
   int typeMax = (type < 0) ? 1 : type;
   for (int type = typeMin; type <= typeMax; ++type) {
     edm::Handle<edm::PCaloHitContainer> hitsCalo;
     e.getByToken(toks_calo_[type], hitsCalo);
     bool getHits = (hitsCalo.isValid());
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HitStudy") << "EcalSimHitStudy: Input flags Hits " << getHits << " with " << hitsCalo->size()
                                  << " hits";
 #endif
     if (getHits) {
       std::vector<PCaloHit> caloHits;
       caloHits.insert(caloHits.end(), hitsCalo->begin(), hitsCalo->end());
       if (!caloHits.empty())
         analyzeHits(caloHits, type);
     }
   }
 }

◆ analyzeHits()

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

protected

Definition at line 307 of file EcalSimHitStudy.cc.

References funct::abs(), TauDecayModes::dec, makePileupJSON::denom, HLT_2024v12_cff::dEta, HLT_2024v12_cff::dPhi, PVValHelper::dz, edep_, edepAll_, edepEM_, edepHad_, hcalRecHitTable_cff::energy, etot_, etotg_, geometry_, CaloGeometry::getGeometry(), hit_, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, hit::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), hcalRecHitTable_cff::iphi, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), dqmdumpme::k, dqm-mbProfile::log, M_PI, SiStripPI::max, position, poszn_, poszp_, r1by25_, r1by9_, r9by25_, singleTopDQM_cfi::select, selX_, sEtaEta_, sEtaPhi_, sPhiPhi_, hcalRecHitTable_cff::time, time_, timeAll_, tmax_, w(), w0_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and ecaldqm::zside().

Referenced by analyze().

                                                                      {
   unsigned int nEC(0);
   std::map<unsigned int, EcalHit> hitMap;
   double etot(0), etotG(0);
   for (auto hit : hits) {
     double edep = hit.energy();
     double time = hit.time();
     unsigned int id_ = hit.id();
     double edepEM = hit.energyEM();
     double edepHad = hit.energyHad();
     if (indx == 0) {
       if ((hit.depth() == 1) || (hit.depth() == 2))
         continue;
     }
     if (time <= tmax_) {
       auto it = hitMap.find(id_);
       if (it == hitMap.end()) {
         uint16_t dep = hit.depth();
         hitMap[id_] = EcalHit(dep, time, edep);
       } else {
         (it->second).energy += edep;
       }
       etotG += edep;
       ++nEC;
     }
     time_[indx]->Fill(time);
     edep_[indx]->Fill(edep);
     edepEM_[indx]->Fill(edepEM);
     edepHad_[indx]->Fill(edepHad);
     etot += edep;
   }
 
   double edepM(0);
   unsigned int idM(0);
   uint16_t depM(0);
   for (auto it : hitMap) {
     if (it.second.energy > edepM) {
       idM = it.first;
       edepM = it.second.energy;
       depM = it.second.id;
     }
   }
 
   bool select(true);
   if (selX_ >= 0) {
     if ((depM & 0X4) != 0)
       select = (selX_ > 0);
     else
       select = (selX_ == 0);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HitStudy") << "EcalSimHitStudy::analyzeHits: Index " << indx << " Emax " << edepM << " IDMax "
                                << std::hex << idM << ":" << depM << std::dec << " Select " << select << ":" << selX_
                                << " Hits " << hits.size() << ":" << nEC << ":" << hitMap.size() << " ETotal " << etot
                                << ":" << etotG;
 #endif
   if (select) {
     etot_[indx]->Fill(etot);
     etotg_[indx]->Fill(etotG);
     hit_[indx]->Fill(double(nEC));
     for (auto it : hitMap) {
       timeAll_[indx]->Fill((it.second).time);
       edepAll_[indx]->Fill((it.second).energy);
       DetId id(it.first);
       if (indx == 0) {
         if (EBDetId(id).zside() >= 0)
           poszp_[indx]->Fill(EBDetId(id).ietaAbs(), EBDetId(id).iphi());
         else
           poszn_[indx]->Fill(EBDetId(id).ietaAbs(), EBDetId(id).iphi());
       } else {
         GlobalPoint gpos = geometry_->getGeometry(id)->getPosition();
         if (EEDetId(id).zside() >= 0)
           poszp_[indx]->Fill(gpos.x(), gpos.y());
         else
           poszn_[indx]->Fill(gpos.x(), gpos.y());
       }
     }
 
     math::XYZVector meanPosition(0.0, 0.0, 0.0);
     std::vector<math::XYZVector> position;
     std::vector<double> energy;
     double e9(0), e25(0);
     for (auto it : hitMap) {
       DetId id(it.first);
       int deta(99), dphi(99), dz(0);
       if (indx == 0) {
         deta = std::abs(EBDetId(id).ietaAbs() - EBDetId(idM).ietaAbs());
         dphi = std::abs(EBDetId(id).iphi() - EBDetId(idM).iphi());
         if (dphi > 180)
           dphi = std::abs(dphi - 360);
         dz = std::abs(EBDetId(id).zside() - EBDetId(idM).zside());
       } else {
         deta = std::abs(EEDetId(id).ix() - EEDetId(idM).ix());
         dphi = std::abs(EEDetId(id).iy() - EEDetId(idM).iy());
         dz = std::abs(EEDetId(id).zside() - EEDetId(idM).zside());
       }
       if (deta <= 1 && dphi <= 1 && dz < 1)
         e9 += (it.second).energy;
       if (deta <= 2 && dphi <= 2 && dz < 1) {
         e25 += (it.second).energy;
         GlobalPoint gpos = geometry_->getGeometry(id)->getPosition();
         math::XYZVector pos(gpos.x(), gpos.y(), gpos.z());
         meanPosition += (it.second).energy * pos;
         position.push_back(pos);
         energy.push_back((it.second).energy);
       }
     }
     double r1by9 = (e9 > 0) ? (edepM / e9) : -1;
     double r1by25 = (e25 > 0) ? (edepM / e25) : -1;
     double r9by25 = (e25 > 0) ? (e9 / e25) : -1;
 
     meanPosition /= e25;
     double denom(0), numEtaEta(0), numEtaPhi(0), numPhiPhi(0);
     for (unsigned int k = 0; k < position.size(); ++k) {
       double dEta = position[k].eta() - meanPosition.eta();
       double dPhi = position[k].phi() - meanPosition.phi();
       if (dPhi > +M_PI) {
         dPhi = 2 * M_PI - dPhi;
       }
       if (dPhi < -M_PI) {
         dPhi = 2 * M_PI + dPhi;
       }
 
       double w = std::max(0.0, (w0_ + std::log(energy[k] / e25)));
       denom += w;
       numEtaEta += std::abs(w * dEta * dEta);
       numEtaPhi += std::abs(w * dEta * dPhi);
       numPhiPhi += std::abs(w * dPhi * dPhi);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HitStudy") << "[" << k << "] dEta " << dEta << " dPhi " << dPhi << " Wt " << energy[k] / e25
                                    << ":" << std::log(energy[k] / e25) << ":" << w;
 #endif
     }
     double sEtaEta = (denom > 0) ? (numEtaEta / denom) : -1.0;
     double sEtaPhi = (denom > 0) ? (numEtaPhi / denom) : -1.0;
     double sPhiPhi = (denom > 0) ? (numPhiPhi / denom) : -1.0;
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HitStudy") << "EcalSimHitStudy::Ratios " << r1by9 << " : " << r1by25 << " : " << r9by25
                                  << " Covariances " << sEtaEta << " : " << sEtaPhi << " : " << sPhiPhi;
 #endif
     r1by9_[indx]->Fill(r1by9);
     r1by25_[indx]->Fill(r1by25);
     r9by25_[indx]->Fill(r9by25);
     sEtaEta_[indx]->Fill(sEtaEta);
     sEtaPhi_[indx]->Fill(sEtaPhi);
     sPhiPhi_[indx]->Fill(sPhiPhi);
   }
 }

◆ beginJob()

void EcalSimHitStudy::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 114 of file EcalSimHitStudy.cc.

References edep_, edepAll_, edepEM_, edepHad_, eneInc_, etaInc_, etot_, etotg_, hit_, mps_fire::i, maxEnergy_, Skims_PA_cff::name, ndets_, phiInc_, poszn_, poszp_, ptInc_, r1by25_, r1by9_, r9by25_, sEtaEta_, sEtaPhi_, sPhiPhi_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, time_, timeAll_, 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";
   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");
   std::string dets[ndets_] = {"EB", "EE"};
   for (int i = 0; i < ndets_; i++) {
     sprintf(name, "Hit%d", i);
     sprintf(title, "Number of hits in %s", dets[i].c_str());
     hit_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0., 20000.);
     hit_[i]->GetXaxis()->SetTitle(title);
     hit_[i]->GetYaxis()->SetTitle("Events");
     hit_[i]->Sumw2();
     sprintf(name, "Time%d", i);
     sprintf(title, "Time of the hit (ns) in %s", dets[i].c_str());
     time_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0., 1000.);
     time_[i]->GetXaxis()->SetTitle(title);
     time_[i]->GetYaxis()->SetTitle("Hits");
     time_[i]->Sumw2();
     sprintf(name, "TimeAll%d", i);
     sprintf(title, "Hit time (ns) in %s (for first hit in the cell)", dets[i].c_str());
     timeAll_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0., 1000.);
     timeAll_[i]->GetXaxis()->SetTitle(title);
     timeAll_[i]->GetYaxis()->SetTitle("Hits");
     timeAll_[i]->Sumw2();
     sprintf(name, "Edep%d", i);
     sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
     edep_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
     edep_[i]->GetXaxis()->SetTitle(title);
     edep_[i]->GetYaxis()->SetTitle("Hits");
     edep_[i]->Sumw2();
     sprintf(name, "EdepAll%d", i);
     sprintf(title, "Total Energy deposit in the cell (GeV) in %s", dets[i].c_str());
     edepAll_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
     edepAll_[i]->GetXaxis()->SetTitle(title);
     edepAll_[i]->GetYaxis()->SetTitle("Hits");
     edepAll_[i]->Sumw2();
     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, dets[i].c_str(), 5000, 0., maxEnergy_);
     edepEM_[i]->GetXaxis()->SetTitle(title);
     edepEM_[i]->GetYaxis()->SetTitle("Hits");
     edepEM_[i]->Sumw2();
     sprintf(name, "EdepHad%d", i);
     sprintf(title, "Energy deposit (GeV) by hadrons in %s", dets[i].c_str());
     edepHad_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
     edepHad_[i]->GetXaxis()->SetTitle(title);
     edepHad_[i]->GetYaxis()->SetTitle("Hits");
     edepHad_[i]->Sumw2();
     sprintf(name, "Etot%d", i);
     sprintf(title, "Total energy deposit (GeV) in %s", dets[i].c_str());
     etot_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
     etot_[i]->GetXaxis()->SetTitle(title);
     etot_[i]->GetYaxis()->SetTitle("Events");
     etot_[i]->Sumw2();
     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, dets[i].c_str(), 5000, 0., maxEnergy_);
     etotg_[i]->GetXaxis()->SetTitle(title);
     etotg_[i]->GetYaxis()->SetTitle("Events");
     etotg_[i]->Sumw2();
     sprintf(name, "r1by9%d", i);
     sprintf(title, "E1/E9 in %s", dets[i].c_str());
     r1by9_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 100, 0.0, 1.0);
     r1by9_[i]->GetXaxis()->SetTitle(title);
     r1by9_[i]->GetYaxis()->SetTitle("Events");
     r1by9_[i]->Sumw2();
     sprintf(name, "r1by25%d", i);
     sprintf(title, "E1/E25 in %s", dets[i].c_str());
     r1by25_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 100, 0.0, 1.0);
     r1by25_[i]->GetXaxis()->SetTitle(title);
     r1by25_[i]->GetYaxis()->SetTitle("Events");
     r1by25_[i]->Sumw2();
     sprintf(name, "r9by25%d", i);
     sprintf(title, "E9/E25 in %s", dets[i].c_str());
     r9by25_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 100, 0.0, 1.0);
     r9by25_[i]->GetXaxis()->SetTitle(title);
     r9by25_[i]->GetYaxis()->SetTitle("Events");
     r9by25_[i]->Sumw2();
     double ymax = (i == 0) ? 0.0005 : 0.005;
     sprintf(name, "sEtaEta%d", i);
     sprintf(title, "Cov(#eta,#eta) in %s", dets[i].c_str());
     sEtaEta_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0.0, ymax);
     sEtaEta_[i]->GetXaxis()->SetTitle(title);
     sEtaEta_[i]->GetYaxis()->SetTitle("Events");
     sEtaEta_[i]->Sumw2();
     sprintf(name, "sEtaPhi%d", i);
     sprintf(title, "Cov(#eta,#phi) in %s", dets[i].c_str());
     sEtaPhi_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0.0, ymax);
     sEtaPhi_[i]->GetXaxis()->SetTitle(title);
     sEtaPhi_[i]->GetYaxis()->SetTitle("Events");
     sEtaPhi_[i]->Sumw2();
     ymax = (i == 0) ? 0.001 : 0.01;
     sprintf(name, "sPhiPhi%d", i);
     sprintf(title, "Cov(#phi,#phi) in %s", dets[i].c_str());
     sPhiPhi_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0.0, ymax);
     sPhiPhi_[i]->GetXaxis()->SetTitle(title);
     sPhiPhi_[i]->GetYaxis()->SetTitle("Events");
     sPhiPhi_[i]->Sumw2();
     if (i == 0) {
       sprintf(title, "%s+", dets[i].c_str());
       poszp_[i] = tfile->make<TH2F>("poszp0", title, 100, 0, 100, 360, 0, 360);
       poszp_[i]->GetXaxis()->SetTitle("i#eta");
       poszp_[i]->GetYaxis()->SetTitle("i#phi");
       sprintf(title, "%s-", dets[i].c_str());
       poszn_[i] = tfile->make<TH2F>("poszn0", title, 100, 0, 100, 360, 0, 360);
       poszn_[i]->GetXaxis()->SetTitle("i#eta");
       poszn_[i]->GetYaxis()->SetTitle("i#phi");
     } else {
       sprintf(title, "%s+", dets[i].c_str());
       poszp_[i] = tfile->make<TH2F>("poszp1", title, 100, -200, 200, 100, -200, 200);
       poszp_[i]->GetXaxis()->SetTitle("x (cm)");
       poszp_[i]->GetYaxis()->SetTitle("y (cm)");
       sprintf(title, "%s-", dets[i].c_str());
       poszn_[i] = tfile->make<TH2F>("poszn1", title, 100, -200, 200, 100, -200, 200);
       poszn_[i]->GetXaxis()->SetTitle("x (cm)");
       poszn_[i]->GetYaxis()->SetTitle("y (cm)");
     }
     poszp_[i]->GetYaxis()->SetTitleOffset(1.2);
     poszp_[i]->Sumw2();
     poszn_[i]->GetYaxis()->SetTitleOffset(1.2);
     poszn_[i]->Sumw2();
   }
 }

◆ beginRun()

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

inlineoverrideprotected

Definition at line 55 of file EcalSimHitStudy.cc.

55 {}

◆ endRun()

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

inlineoverrideprotected

Definition at line 56 of file EcalSimHitStudy.cc.

56 {}

◆ fillDescriptions()

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

static

Definition at line 102 of file EcalSimHitStudy.cc.

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

                                                                                  {
   edm::ParameterSetDescription desc;
   std::vector<std::string> calonames;
   desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
   desc.addUntracked<std::vector<std::string>>("CaloCollection", calonames);
   desc.addUntracked<std::string>("SourceLabel", "VtxSmeared");
   desc.addUntracked<double>("MaxEnergy", 200.0);
   desc.addUntracked<double>("TimeCut", 100.0);
   desc.addUntracked<int>("SelectX", -1);
   descriptions.add("EcalSimHitStudy", desc);
 }

Member Data Documentation

◆ edep_

TH1F * EcalSimHitStudy::edep_[ndets_]

private

Definition at line 76 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepAll_

TH1F * EcalSimHitStudy::edepAll_[ndets_]

private

Definition at line 77 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepEM_

TH1F* EcalSimHitStudy::edepEM_[ndets_]

private

Definition at line 76 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepHad_

TH1F * EcalSimHitStudy::edepHad_[ndets_]

private

Definition at line 76 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eneInc_

TH1F * EcalSimHitStudy::eneInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

◆ etaInc_

TH1F * EcalSimHitStudy::etaInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

◆ etot_

TH1F* EcalSimHitStudy::etot_[ndets_]

private

Definition at line 77 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ etotg_

TH1F * EcalSimHitStudy::etotg_[ndets_]

private

Definition at line 77 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ g4Label_

const std::string EcalSimHitStudy::g4Label_

private

Definition at line 67 of file EcalSimHitStudy.cc.

Referenced by EcalSimHitStudy().

◆ geometry_

const CaloGeometry* EcalSimHitStudy::geometry_

private

Definition at line 73 of file EcalSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

◆ hit_

TH1F* EcalSimHitStudy::hit_[ndets_]

private

Definition at line 75 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ hitLab_

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

private

Definition at line 68 of file EcalSimHitStudy.cc.

Referenced by EcalSimHitStudy().

◆ maxEnergy_

const double EcalSimHitStudy::maxEnergy_

private

Definition at line 69 of file EcalSimHitStudy.cc.

Referenced by beginJob().

◆ ndets_

const int EcalSimHitStudy::ndets_ = 2

staticprivate

Definition at line 65 of file EcalSimHitStudy.cc.

Referenced by beginJob().

◆ phiInc_

TH1F * EcalSimHitStudy::phiInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

◆ poszn_

TH2F * EcalSimHitStudy::poszn_[ndets_]

private

Definition at line 80 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ poszp_

TH2F* EcalSimHitStudy::poszp_[ndets_]

private

Definition at line 80 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ ptInc_

TH1F* EcalSimHitStudy::ptInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

◆ r1by25_

TH1F * EcalSimHitStudy::r1by25_[ndets_]

private

Definition at line 78 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ r1by9_

TH1F* EcalSimHitStudy::r1by9_[ndets_]

private

Definition at line 78 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ r9by25_

TH1F * EcalSimHitStudy::r9by25_[ndets_]

private

Definition at line 78 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ selX_

const int EcalSimHitStudy::selX_

private

Definition at line 70 of file EcalSimHitStudy.cc.

Referenced by analyzeHits().

◆ sEtaEta_

TH1F* EcalSimHitStudy::sEtaEta_[ndets_]

private

Definition at line 79 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ sEtaPhi_

TH1F * EcalSimHitStudy::sEtaPhi_[ndets_]

private

Definition at line 79 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ sPhiPhi_

TH1F * EcalSimHitStudy::sPhiPhi_[ndets_]

private

Definition at line 79 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ time_

TH1F * EcalSimHitStudy::time_[ndets_]

private

Definition at line 75 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ timeAll_

TH1F * EcalSimHitStudy::timeAll_[ndets_]

private

Definition at line 75 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ tmax_

const double EcalSimHitStudy::tmax_

private

Definition at line 69 of file EcalSimHitStudy.cc.

Referenced by analyzeHits().

◆ tok_evt_

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

private

Definition at line 71 of file EcalSimHitStudy.cc.

Referenced by analyze().

◆ tokGeom_

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

private

Definition at line 66 of file EcalSimHitStudy.cc.

Referenced by analyze().

◆ toks_calo_

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

private

Definition at line 72 of file EcalSimHitStudy.cc.

Referenced by analyze().

◆ w0_

const double EcalSimHitStudy::w0_

private

Definition at line 69 of file EcalSimHitStudy.cc.

Referenced by analyzeHits().

Classes

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalSimHitStudy()

◆ ~EcalSimHitStudy()

Member Function Documentation

◆ analyze()

◆ analyzeHits()

◆ beginJob()

◆ beginRun()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ edep_

◆ edepAll_

◆ edepEM_

◆ edepHad_

◆ eneInc_

◆ etaInc_

◆ etot_

◆ etotg_

◆ g4Label_

◆ geometry_

◆ hit_

◆ hitLab_

◆ maxEnergy_

◆ ndets_

◆ phiInc_

◆ poszn_

◆ poszp_

◆ ptInc_

◆ r1by25_

◆ r1by9_

◆ r9by25_

◆ selX_

◆ sEtaEta_

◆ sEtaPhi_

◆ sPhiPhi_

◆ time_

◆ timeAll_

◆ tmax_

◆ tok_evt_

◆ tokGeom_

◆ toks_calo_

◆ w0_