Inheritance diagram for HcalTB06Analysis:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	beginJob () override

void	endJob () override

	HcalTB06Analysis (const edm::ParameterSet &p)

	HcalTB06Analysis (const HcalTB06Analysis &)=delete

const HcalTB06Analysis &	operator= (const HcalTB06Analysis &)=delete

	~HcalTB06Analysis () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
int	count

const edm::EDGetTokenT< edm::PCaloHitContainer >	m_BeamToken

const bool	m_ECAL

edm::EDGetTokenT< edm::PCaloHitContainer >	m_EcalToken

const double	m_eMIP

const double	m_ener

const double	m_eta

const double	m_factEcal

const double	m_factHcal

const edm::EDGetTokenT< edm::PCaloHitContainer >	m_HcalToken

std::unique_ptr< HcalTB06Histo >	m_histo

int	m_idxetaEcal

int	m_idxetaHcal

int	m_idxphiEcal

int	m_idxphiHcal

const std::vector< int >	m_PDG

const double	m_phi

const edm::ParameterSet	m_ptb

const double	m_timeLimit

const double	m_widthEcal

const double	m_widthHcal

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Detailed Description

Definition at line 58 of file HcalTB06Analysis.cc.

Constructor & Destructor Documentation

◆ HcalTB06Analysis() [1/2]

HcalTB06Analysis::HcalTB06Analysis ( const edm::ParameterSet & p )

explicit

Definition at line 98 of file HcalTB06Analysis.cc.

References ProducerED_cfi::InputTag, m_ECAL, m_EcalToken, m_eMIP, m_ener, m_eta, m_factEcal, m_factHcal, m_histo, m_idxetaEcal, m_idxetaHcal, m_idxphiEcal, m_idxphiHcal, m_PDG, m_phi, m_ptb, m_timeLimit, m_widthEcal, m_widthHcal, razorScouting_cff::maxEta, HLT_2024v14_cff::maxPhi, EgHLTOffEleSelection_cfi::minEta, HLT_2024v14_cff::minPhi, and AlCaHLTBitMon_ParallelJobs::p.

     : m_ECAL(p.getParameter<bool>("ECAL")),
       m_eta(p.getParameter<double>("MinEta")),
       m_phi(p.getParameter<double>("MinPhi")),
       m_ener(p.getParameter<double>("MinE")),
       m_PDG(p.getParameter<std::vector<int> >("PartID")),
       m_HcalToken(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalHits"))),
       m_BeamToken(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalTB06BeamHits"))),
       m_ptb(p.getParameter<edm::ParameterSet>("TestBeamAnalysis")),
       m_timeLimit(m_ptb.getParameter<double>("TimeLimit")),
       m_widthEcal(m_ptb.getParameter<double>("EcalWidth")),
       m_widthHcal(m_ptb.getParameter<double>("HcalWidth")),
       m_factEcal(m_ptb.getParameter<double>("EcalFactor")),
       m_factHcal(m_ptb.getParameter<double>("HcalFactor")),
       m_eMIP(m_ptb.getParameter<double>("MIP")),
       count(0) {
   usesResource("TFileService");
   if (m_ECAL)
     m_EcalToken = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEB"));
   double minEta = p.getParameter<double>("MinEta");
   double maxEta = p.getParameter<double>("MaxEta");
   double minPhi = p.getParameter<double>("MinPhi");
   double maxPhi = p.getParameter<double>("MaxPhi");
   double beamEta = (maxEta + minEta) * 0.5;
   double beamPhi = (maxPhi + minPhi) * 0.5;
   if (beamPhi < 0) {
     beamPhi += twopi;
   }
 
   m_idxetaEcal = 13;
   m_idxphiEcal = 13;
 
   m_idxetaHcal = static_cast<int>(beamEta / 0.087) + 1;
   m_idxphiHcal = static_cast<int>(beamPhi / 0.087) + 6;
   if (m_idxphiHcal > 72) {
     m_idxphiHcal -= 73;
   }
 
   edm::LogInfo("HcalTB06Analysis") << "Beam parameters: E(GeV)= " << m_ener << " pdgID= " << m_PDG[0]
                                    << "\n eta= " << m_eta << "  idx_etaEcal= " << m_idxetaEcal
                                    << "  idx_etaHcal= " << m_idxetaHcal << "  phi= " << m_phi
                                    << "  idx_phiEcal= " << m_idxphiEcal << "  idx_phiHcal= " << m_idxphiHcal
                                    << "\n        EcalFactor= " << m_factEcal << "  EcalWidth= " << m_widthEcal << " GeV"
                                    << "\n        HcalFactor= " << m_factHcal << "  HcalWidth= " << m_widthHcal << " GeV"
                                    << "  MIP=       " << m_eMIP << " GeV"
                                    << "\n        TimeLimit=  " << m_timeLimit << " ns"
                                    << "\n";
   m_histo = std::make_unique<HcalTB06Histo>(m_ptb);
 }

◆ ~HcalTB06Analysis()

HcalTB06Analysis::~HcalTB06Analysis ( )

overridedefault

◆ HcalTB06Analysis() [2/2]

HcalTB06Analysis::HcalTB06Analysis ( const HcalTB06Analysis & )

delete

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 154 of file HcalTB06Analysis.cc.

References funct::abs(), count, edm::Event::getHandle(), HcalOuter, mps_fire::i, EBDetId::ieta(), HcalDetId::ieta(), EBDetId::iphi(), HcalDetId::iphi(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), m_BeamToken, m_ECAL, m_EcalToken, m_ener, m_eta, m_factEcal, m_factHcal, m_HcalToken, m_histo, m_idxetaEcal, m_idxetaHcal, m_idxphiEcal, m_idxphiHcal, m_phi, m_timeLimit, m_widthEcal, m_widthHcal, compareTotals::ne, cms::cuda::nh, edm::Handle< T >::product(), HcalDetId::subdet(), hcalRecHitTable_cff::time, and HcalTestBeamNumbering::unpackIndex().

                                                                         {
   ++count;
 
   //Beam Information
   m_histo->fillPrimary(m_ener, m_eta, m_phi);
 
   std::vector<double> eCalo(6, 0), eTrig(7, 0);
 
   const std::vector<PCaloHit>* EcalHits = nullptr;
   if (m_ECAL) {
     const edm::Handle<edm::PCaloHitContainer>& Ecal = evt.getHandle(m_EcalToken);
     EcalHits = Ecal.product();
   }
   const edm::Handle<edm::PCaloHitContainer>& Hcal = evt.getHandle(m_HcalToken);
   const std::vector<PCaloHit>* HcalHits = Hcal.product();
   const edm::Handle<edm::PCaloHitContainer>& Beam = evt.getHandle(m_BeamToken);
   const std::vector<PCaloHit>* BeamHits = Beam.product();
 
   // Total Energy
   double eecals = 0.;
   double ehcals = 0.;
 
   unsigned int ne = 0;
   unsigned int nh = 0;
   if (m_ECAL) {
     ne = EcalHits->size();
     for (unsigned int i = 0; i < ne; ++i) {
       EBDetId ecalid((*EcalHits)[i].id());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalTBSim") << "EB " << i << " " << ecalid.ieta() << ":" << m_idxetaEcal << "   "
                                     << ecalid.iphi() << ":" << m_idxphiEcal << "   " << (*EcalHits)[i].time() << ":"
                                     << m_timeLimit << "   " << (*EcalHits)[i].energy();
 #endif
       // 7x7 crystal selection
       if (std::abs(m_idxetaEcal - ecalid.ieta()) <= 3 && std::abs(m_idxphiEcal - ecalid.iphi()) <= 3 &&
           (*EcalHits)[i].time() < m_timeLimit) {
         eCalo[0] += (*EcalHits)[i].energy();
       }
     }
     if (m_widthEcal > 0.0) {
       eCalo[1] = G4RandGauss::shoot(0.0, m_widthEcal);
     }
     eecals = m_factEcal * (eCalo[0] + eCalo[1]);
   }
   if (HcalHits) {
     nh = HcalHits->size();
     for (unsigned int i = 0; i < nh; ++i) {
       HcalDetId hcalid((*HcalHits)[i].id());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalTBSim") << "HC " << i << " " << hcalid.subdet() << "  " << hcalid.ieta() << ":"
                                     << m_idxetaHcal << "   " << hcalid.iphi() << ":" << m_idxphiHcal << "   "
                                     << (*HcalHits)[i].time() << ":" << m_timeLimit << "   " << (*HcalHits)[i].energy();
 #endif
       // 3x3 towers selection
       if (std::abs(m_idxetaHcal - hcalid.ieta()) <= 1 && std::abs(m_idxphiHcal - hcalid.iphi()) <= 1 &&
           (*HcalHits)[i].time() < m_timeLimit) {
         if (hcalid.subdet() != HcalOuter) {
           eCalo[2] += (*HcalHits)[i].energy();
         } else {
           eCalo[4] += (*HcalHits)[i].energy();
         }
       }
     }
     if (m_widthHcal > 0.0) {
       eCalo[3] = G4RandGauss::shoot(0.0, m_widthHcal);
       eCalo[5] = G4RandGauss::shoot(0.0, m_widthHcal);
     }
     ehcals = m_factHcal * eCalo[2] + eCalo[3];
   }
   double etots = eecals + ehcals;
 
   edm::LogInfo("HcalTBSim") << "HcalTB06Analysis:: Etot(MeV)= " << etots << "   E(Ecal)= " << eecals
                             << "   E(Hcal)= " << ehcals << "  Nhits(ECAL)= " << ne << "  Nhits(HCAL)= " << nh;
   m_histo->fillEdep(etots, eecals, ehcals);
 
   if (BeamHits) {
     for (unsigned int i = 0; i < BeamHits->size(); ++i) {
       unsigned int id = ((*BeamHits)[i].id());
       int det, lay, ix, iy;
       HcalTestBeamNumbering::unpackIndex(id, det, lay, ix, iy);
       if ((det == 1) && ((*BeamHits)[i].time() < m_timeLimit)) {
         if (lay > 0 && lay <= 4) {
           eTrig[lay - 1] += (*BeamHits)[i].energy();
         } else if (lay == 7 || lay == 8) {
           eTrig[lay - 2] += (*BeamHits)[i].energy();
         } else if (lay >= 11 && lay <= 14) {
           eTrig[4] += (*BeamHits)[i].energy();
         }
       }
     }
   }
 
   edm::LogInfo("HcalTBSim") << "HcalTB06Analysis:: Trigger Info: " << eTrig[0] << ":" << eTrig[1] << ":" << eTrig[2]
                             << ":" << eTrig[3] << ":" << eTrig[4] << ":" << eTrig[5] << ":" << eTrig[6];
 
   m_histo->fillTree(eCalo, eTrig);
 }

◆ beginJob()

void HcalTB06Analysis::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 148 of file HcalTB06Analysis.cc.

148 { edm::LogInfo("HcalTB06Analysis") << " =====> Begin of Run"; }

edm::LogInfo

Log< level::Info, false > LogInfo

Definition: MessageLogger.h:131

◆ endJob()

void HcalTB06Analysis::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 150 of file HcalTB06Analysis.cc.

References count.

                               {
   edm::LogInfo("HcalTB06Analysis") << " =====> End of Run; Total number of events: " << count;
 }

◆ operator=()

const HcalTB06Analysis& HcalTB06Analysis::operator= ( const HcalTB06Analysis & )

delete

Member Data Documentation

◆ count

int HcalTB06Analysis::count

private

Definition at line 89 of file HcalTB06Analysis.cc.

Referenced by analyze(), and endJob().

◆ m_BeamToken

const edm::EDGetTokenT<edm::PCaloHitContainer> HcalTB06Analysis::m_BeamToken

private

Definition at line 79 of file HcalTB06Analysis.cc.

Referenced by analyze().

◆ m_ECAL

const bool HcalTB06Analysis::m_ECAL

private

Definition at line 71 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_EcalToken

edm::EDGetTokenT<edm::PCaloHitContainer> HcalTB06Analysis::m_EcalToken

private

Definition at line 77 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_eMIP

const double HcalTB06Analysis::m_eMIP

private

Definition at line 87 of file HcalTB06Analysis.cc.

Referenced by HcalTB06Analysis().

◆ m_ener

const double HcalTB06Analysis::m_ener

private

Definition at line 74 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_eta

const double HcalTB06Analysis::m_eta

private

Definition at line 72 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_factEcal

const double HcalTB06Analysis::m_factEcal

private

Definition at line 85 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_factHcal

const double HcalTB06Analysis::m_factHcal

private

Definition at line 86 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_HcalToken

const edm::EDGetTokenT<edm::PCaloHitContainer> HcalTB06Analysis::m_HcalToken

private

Definition at line 78 of file HcalTB06Analysis.cc.

Referenced by analyze().

◆ m_histo

std::unique_ptr<HcalTB06Histo> HcalTB06Analysis::m_histo

private

Definition at line 95 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_idxetaEcal

int HcalTB06Analysis::m_idxetaEcal

private

Definition at line 90 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_idxetaHcal

int HcalTB06Analysis::m_idxetaHcal

private

Definition at line 92 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_idxphiEcal

int HcalTB06Analysis::m_idxphiEcal

private

Definition at line 91 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_idxphiHcal

int HcalTB06Analysis::m_idxphiHcal

private

Definition at line 93 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_PDG

const std::vector<int> HcalTB06Analysis::m_PDG

private

Definition at line 75 of file HcalTB06Analysis.cc.

Referenced by HcalTB06Analysis().

◆ m_phi

const double HcalTB06Analysis::m_phi

private

Definition at line 73 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_ptb

const edm::ParameterSet HcalTB06Analysis::m_ptb

private

Definition at line 81 of file HcalTB06Analysis.cc.

Referenced by HcalTB06Analysis().

◆ m_timeLimit

const double HcalTB06Analysis::m_timeLimit

private

Definition at line 82 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_widthEcal

const double HcalTB06Analysis::m_widthEcal

private

Definition at line 83 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

◆ m_widthHcal

const double HcalTB06Analysis::m_widthHcal

private

Definition at line 84 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalTB06Analysis() [1/2]

◆ ~HcalTB06Analysis()

◆ HcalTB06Analysis() [2/2]

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ operator=()

Member Data Documentation

◆ count

◆ m_BeamToken

◆ m_ECAL

◆ m_EcalToken

◆ m_eMIP

◆ m_ener

◆ m_eta

◆ m_factEcal

◆ m_factHcal

◆ m_HcalToken

◆ m_histo

◆ m_idxetaEcal

◆ m_idxetaHcal

◆ m_idxphiEcal

◆ m_idxphiHcal

◆ m_PDG

◆ m_phi

◆ m_ptb

◆ m_timeLimit

◆ m_widthEcal

◆ m_widthHcal