#include <HcalTB06Analysis.h>

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

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () 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

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::vector< ModuleDescription const * > &modules, 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

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

edm::EDGetTokenT< edm::PCaloHitContainer >	m_BeamToken

bool	m_ECAL

edm::EDGetTokenT< edm::PCaloHitContainer >	m_EcalToken

double	m_ener

double	m_eta

double	m_factEcal

double	m_factHcal

edm::EDGetTokenT< edm::PCaloHitContainer >	m_HcalToken

HcalTB06Histo *	m_histo

int	m_idxetaEcal

int	m_idxetaHcal

int	m_idxphiEcal

int	m_idxphiHcal

std::vector< int >	m_PDG

double	m_phi

double	m_timeLimit

double	m_widthEcal

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)

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

Detailed Description

Definition at line 26 of file HcalTB06Analysis.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 47 of file HcalTB06Analysis.cc.

References edm::ParameterSet::getParameter(), createfilelist::int, 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_PDG, m_phi, m_timeLimit, m_widthEcal, m_widthHcal, maxEta, trackingParticleSelector_cfi::maxPhi, trackingPOGFilters_cfi::minEta, and trackingParticleSelector_cfi::minPhi.

                                                            : count(0) {
 
   usesResource("TFileService");
 
   m_ECAL = p.getParameter<bool>("ECAL");
   if(m_ECAL) {
     m_EcalToken = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits","EcalHitsEB"));
   }
   m_HcalToken = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits","HcalHits"));
   m_BeamToken = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits","HcalTB06BeamHits"));
   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");
 
   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   = (int)(beamEta/0.087) + 1;
   m_idxphiHcal   = (int)(beamPhi/0.087) + 6;
   if(m_idxphiHcal > 72) { m_idxphiHcal -= 73; }
 
   edm::ParameterSet ptb = p.getParameter<edm::ParameterSet>("TestBeamAnalysis");
   m_timeLimit = ptb.getParameter<double>("TimeLimit");
   m_widthEcal = ptb.getParameter<double>("EcalWidth");
   m_widthHcal = ptb.getParameter<double>("HcalWidth");
   m_factEcal  = ptb.getParameter<double>("EcalFactor");
   m_factHcal  = ptb.getParameter<double>("HcalFactor");
   double eMIP = ptb.getParameter<double>("MIP");
 
   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=       " << eMIP << " GeV"
     << "\n        TimeLimit=  " << m_timeLimit << " ns" << "\n";
   m_histo = new HcalTB06Histo(ptb);
 } 

HcalTB06Analysis::~HcalTB06Analysis ( )

override

Definition at line 103 of file HcalTB06Analysis.cc.

References m_histo.

                                     {
   delete m_histo;
 }

HcalTB06Analysis::HcalTB06Analysis ( const HcalTB06Analysis & )

delete

Member Function Documentation

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

override

Definition at line 116 of file HcalTB06Analysis.cc.

References funct::abs(), count, gather_cfg::cout, ModuleToSequenceAssign::Ecal, HcalTB06Histo::fillEdep(), HcalTB06Histo::fillPrimary(), HcalTB06Histo::fillTree(), edm::Event::getByToken(), 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, Vispa.Plugins.EdmBrowser.EdmDataAccessor::ne(), edm::Handle< T >::product(), HcalDetId::subdet(), protons_cff::time, and HcalTestBeamNumbering::unpackIndex().

 {
   ++count;
 
   //Beam Information
   m_histo->fillPrimary(m_ener, m_eta, m_phi);
 
   edm::Handle<edm::PCaloHitContainer> Ecal;
   edm::Handle<edm::PCaloHitContainer> Hcal;
   edm::Handle<edm::PCaloHitContainer> Beam;
   std::vector<double> eCalo(6,0), eTrig(7,0);
 
   const std::vector<PCaloHit>* EcalHits = nullptr;
   if(m_ECAL) { 
     evt.getByToken(m_EcalToken, Ecal); 
     EcalHits = Ecal.product();
   }
   evt.getByToken(m_HcalToken, Hcal);
   const std::vector<PCaloHit>* HcalHits = Hcal.product();
   evt.getByToken(m_BeamToken, Beam);
   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
       std::cout << "EB " << i << " " << ecalid.ieta() << ":" << m_idxetaEcal 
         << "   " << ecalid.iphi() << ":" << m_idxphiEcal << "   " 
         << (*EcalHits)[i].time() << ":" << m_timeLimit << "   " 
         << (*EcalHits)[i].energy() << std::endl;
 #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
       std::cout << "HC " << i << " " << hcalid.subdet() << "  " 
         << hcalid.ieta() << ":" << m_idxetaHcal << "   " 
         << hcalid.iphi() << ":" << m_idxphiHcal << "   " 
         << (*HcalHits)[i].time() << ":" << m_timeLimit << "   " 
         << (*HcalHits)[i].energy() << std::endl;
 #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 107 of file HcalTB06Analysis.cc.

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

void HcalTB06Analysis::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 111 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 47 of file HcalTB06Analysis.h.

Referenced by analyze(), and endJob().

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

private

Definition at line 44 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

bool HcalTB06Analysis::m_ECAL

private

Definition at line 45 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

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

private

Definition at line 42 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_ener

private

Definition at line 55 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_eta

private

Definition at line 53 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_factEcal

private

Definition at line 59 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_factHcal

private

Definition at line 60 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

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

private

Definition at line 43 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

HcalTB06Histo* HcalTB06Analysis::m_histo

private

Definition at line 63 of file HcalTB06Analysis.h.

Referenced by analyze(), HcalTB06Analysis(), and ~HcalTB06Analysis().

int HcalTB06Analysis::m_idxetaEcal

private

Definition at line 48 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxetaHcal

private

Definition at line 50 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxphiEcal

private

Definition at line 49 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

int HcalTB06Analysis::m_idxphiHcal

private

Definition at line 51 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

std::vector<int> HcalTB06Analysis::m_PDG

private

Definition at line 61 of file HcalTB06Analysis.h.

Referenced by HcalTB06Analysis().

double HcalTB06Analysis::m_phi

private

Definition at line 54 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_timeLimit

private

Definition at line 56 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_widthEcal

private

Definition at line 57 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

double HcalTB06Analysis::m_widthHcal

private

Definition at line 58 of file HcalTB06Analysis.h.

Referenced by analyze(), and HcalTB06Analysis().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation