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

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

Detailed Description

Definition at line 55 of file HcalTB06Analysis.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 95 of file HcalTB06Analysis.cc.

References edm::ParameterSet::getParameter(), HLT_FULL_cff::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, maxEta, HLT_FULL_cff::maxPhi, and HLT_FULL_cff::minPhi.

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

overridedefault

HcalTB06Analysis::HcalTB06Analysis ( const HcalTB06Analysis & )

delete

Member Function Documentation

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 151 of file HcalTB06Analysis.cc.

References funct::abs(), count, edm::Event::getHandle(), HcalOuter, mps_fire::i, EBDetId::ieta(), HcalDetId::ieta(), EBDetId::iphi(), HcalDetId::iphi(), 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(), 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);
 }

void HcalTB06Analysis::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 145 of file HcalTB06Analysis.cc.

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

edm::LogInfo

Log< level::Info, false > LogInfo

Definition: MessageLogger.h:125

void HcalTB06Analysis::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 147 of file HcalTB06Analysis.cc.

References count.

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

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

delete

Member Data Documentation

int HcalTB06Analysis::count

private

Definition at line 86 of file HcalTB06Analysis.cc.

Referenced by analyze(), and endJob().

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

private

Definition at line 76 of file HcalTB06Analysis.cc.

Referenced by analyze().

const bool HcalTB06Analysis::m_ECAL

private

Definition at line 68 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

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

private

Definition at line 74 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const double HcalTB06Analysis::m_eMIP

private

Definition at line 84 of file HcalTB06Analysis.cc.

Referenced by HcalTB06Analysis().

const double HcalTB06Analysis::m_ener

private

Definition at line 71 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const double HcalTB06Analysis::m_eta

private

Definition at line 69 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const double HcalTB06Analysis::m_factEcal

private

Definition at line 82 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const double HcalTB06Analysis::m_factHcal

private

Definition at line 83 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

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

private

Definition at line 75 of file HcalTB06Analysis.cc.

Referenced by analyze().

std::unique_ptr<HcalTB06Histo> HcalTB06Analysis::m_histo

private

Definition at line 92 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxetaEcal

private

Definition at line 87 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxetaHcal

private

Definition at line 89 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxphiEcal

private

Definition at line 88 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxphiHcal

private

Definition at line 90 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

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

private

Definition at line 72 of file HcalTB06Analysis.cc.

Referenced by HcalTB06Analysis().

const double HcalTB06Analysis::m_phi

private

Definition at line 70 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const edm::ParameterSet HcalTB06Analysis::m_ptb

private

Definition at line 78 of file HcalTB06Analysis.cc.

Referenced by HcalTB06Analysis().

const double HcalTB06Analysis::m_timeLimit

private

Definition at line 79 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const double HcalTB06Analysis::m_widthEcal

private

Definition at line 80 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

const double HcalTB06Analysis::m_widthHcal

private

Definition at line 81 of file HcalTB06Analysis.cc.

Referenced by analyze(), and HcalTB06Analysis().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation