HcalTB06Analysis Class Reference

#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
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)
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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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, heepElectronID_HEEPV50_cff::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(), ntuplemaker::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().