EcalSimHitsValidation Class Reference

#include <EcalSimHitsValidation.h>

Inheritance diagram for EcalSimHitsValidation:

Public Member Functions
	EcalSimHitsValidation (const edm::ParameterSet &ps)
	Constructor. More...
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
	Analyze. More...
void	bookHistograms (DQMStore::IBooker &ib, edm::Run const &, edm::EventSetup const &c) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType

Private Attributes
edm::EDGetTokenT< edm::PCaloHitContainer >	EBHitsCollectionToken
edm::EDGetTokenT< edm::PCaloHitContainer >	EEHitsCollectionToken
edm::EDGetTokenT< edm::PCaloHitContainer >	ESHitsCollectionToken
std::string	g4InfoLabel
edm::EDGetTokenT< edm::HepMCProduct >	HepMCToken
MonitorElement *	meEBEnergyFraction_
MonitorElement *	meEEEnergyFraction_
MonitorElement *	meESEnergyFraction_
MonitorElement *	meGunEnergy_
MonitorElement *	meGunEta_
MonitorElement *	meGunPhi_
bool	verbose_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 37 of file EcalSimHitsValidation.h.

Member Typedef Documentation

MapType

typedef std::map<uint32_t, float, std::less<uint32_t> > EcalSimHitsValidation::MapType

private

Definition at line 38 of file EcalSimHitsValidation.h.

Constructor & Destructor Documentation

EcalSimHitsValidation()

EcalSimHitsValidation::EcalSimHitsValidation

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 20 of file EcalSimHitsValidation.cc.

    : g4InfoLabel(ps.getParameter<std::string>("moduleLabelG4")),
      HepMCToken(consumes<edm::HepMCProduct>(ps.getParameter<std::string>("moduleLabelMC"))) {
  EBHitsCollectionToken =
      consumes<edm::PCaloHitContainer>(edm::InputTag(g4InfoLabel, ps.getParameter<std::string>("EBHitsCollection")));
  EEHitsCollectionToken =
      consumes<edm::PCaloHitContainer>(edm::InputTag(g4InfoLabel, ps.getParameter<std::string>("EEHitsCollection")));
  ESHitsCollectionToken =
      consumes<edm::PCaloHitContainer>(edm::InputTag(g4InfoLabel, ps.getParameter<std::string>("ESHitsCollection")));
 
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
}

References EBHitsCollectionToken, EEHitsCollectionToken, ESHitsCollectionToken, g4InfoLabel, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 57 of file EcalSimHitsValidation.cc.

                                                                             {
  edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
 
  std::vector<PCaloHit> theEBCaloHits;
  std::vector<PCaloHit> theEECaloHits;
  std::vector<PCaloHit> theESCaloHits;
 
  edm::Handle<edm::HepMCProduct> MCEvt;
  edm::Handle<edm::PCaloHitContainer> EcalHitsEB;
  edm::Handle<edm::PCaloHitContainer> EcalHitsEE;
  edm::Handle<edm::PCaloHitContainer> EcalHitsES;
 
  e.getByToken(HepMCToken, MCEvt);
  e.getByToken(EBHitsCollectionToken, EcalHitsEB);
  e.getByToken(EEHitsCollectionToken, EcalHitsEE);
  e.getByToken(ESHitsCollectionToken, EcalHitsES);
 
  for (HepMC::GenEvent::particle_const_iterator p = MCEvt->GetEvent()->particles_begin();
       p != MCEvt->GetEvent()->particles_end();
       ++p) {
    double htheta = (*p)->momentum().theta();
    double heta = -99999.;
    if (tan(htheta * 0.5) > 0) {
      heta = -log(tan(htheta * 0.5));
    }
    double hphi = (*p)->momentum().phi();
    hphi = (hphi >= 0) ? hphi : hphi + 2 * M_PI;
    hphi = hphi / M_PI * 180.;
 
    LogDebug("EventInfo") << "Particle gun type form MC = " << abs((*p)->pdg_id()) << "\n"
                          << "Energy = " << (*p)->momentum().e() << " Eta = " << heta << " Phi = " << hphi;
 
    meGunEnergy_->Fill((*p)->momentum().e());
    meGunEta_->Fill(heta);
    meGunPhi_->Fill(hphi);
  }
 
  double EBEnergy_ = 0.;
  if (EcalHitsEB.isValid()) {
    theEBCaloHits.insert(theEBCaloHits.end(), EcalHitsEB->begin(), EcalHitsEB->end());
    for (std::vector<PCaloHit>::iterator isim = theEBCaloHits.begin(); isim != theEBCaloHits.end(); ++isim) {
      EBEnergy_ += isim->energy();
    }
  }
 
  double EEEnergy_ = 0.;
  if (EcalHitsEE.isValid()) {
    theEECaloHits.insert(theEECaloHits.end(), EcalHitsEE->begin(), EcalHitsEE->end());
    for (std::vector<PCaloHit>::iterator isim = theEECaloHits.begin(); isim != theEECaloHits.end(); ++isim) {
      EEEnergy_ += isim->energy();
    }
  }
 
  double ESEnergy_ = 0.;
  if (EcalHitsES.isValid()) {
    theESCaloHits.insert(theESCaloHits.end(), EcalHitsES->begin(), EcalHitsES->end());
    for (std::vector<PCaloHit>::iterator isim = theESCaloHits.begin(); isim != theESCaloHits.end(); ++isim) {
      ESEnergy_ += isim->energy();
    }
  }
 
  double etot = EBEnergy_ + EEEnergy_ + ESEnergy_;
  double fracEB = 0.0;
  double fracEE = 0.0;
  double fracES = 0.0;
 
  if (etot > 0.0) {
    fracEB = EBEnergy_ / etot;
    fracEE = EEEnergy_ / etot;
    fracES = ESEnergy_ / etot;
  }
 
  meEBEnergyFraction_->Fill(fracEB);
  meEEEnergyFraction_->Fill(fracEE);
  meESEnergyFraction_->Fill(fracES);
}

References funct::abs(), MillePedeFileConverter_cfg::e, EBHitsCollectionToken, EEHitsCollectionToken, ESHitsCollectionToken, dqm::impl::MonitorElement::Fill(), edm::HepMCProduct::GetEvent(), HepMCToken, edm::HandleBase::isValid(), dqm-mbProfile::log, LogDebug, M_PI, meEBEnergyFraction_, meEEEnergyFraction_, meESEnergyFraction_, meGunEnergy_, meGunEta_, meGunPhi_, AlCaHLTBitMon_ParallelJobs::p, and funct::tan().

bookHistograms()

void EcalSimHitsValidation::bookHistograms ( DQMStore::IBooker &  ib, edm::Run const &  , edm::EventSetup const &  c  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 34 of file EcalSimHitsValidation.cc.

                                                                                                      {
  ib.setCurrentFolder("EcalHitsV/EcalSimHitsValidation");
  ib.setScope(MonitorElementData::Scope::RUN);
 
  std::string histo = "EcalSimHitsValidation Gun Momentum";
  meGunEnergy_ = ib.book1D(histo, histo, 100, 0., 1000.);
 
  histo = "EcalSimHitsValidation Gun Eta";
  meGunEta_ = ib.book1D(histo, histo, 700, -3.5, 3.5);
 
  histo = "EcalSimHitsValidation Gun Phi";
  meGunPhi_ = ib.book1D(histo, histo, 360, 0., 360.);
 
  histo = "EcalSimHitsValidation Barrel fraction of energy";
  meEBEnergyFraction_ = ib.book1D(histo, histo, 100, 0., 1.1);
 
  histo = "EcalSimHitsValidation Endcap fraction of energy";
  meEEEnergyFraction_ = ib.book1D(histo, histo, 100, 0., 1.1);
 
  histo = "EcalSimHitsValidation Preshower fraction of energy";
  meESEnergyFraction_ = ib.book1D(histo, histo, 60, 0., 0.001);
}

References timingPdfMaker::histo, cuy::ib, meEBEnergyFraction_, meEEEnergyFraction_, meESEnergyFraction_, meGunEnergy_, meGunEta_, meGunPhi_, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

EBHitsCollectionToken

edm::EDGetTokenT<edm::PCaloHitContainer> EcalSimHitsValidation::EBHitsCollectionToken

private

Definition at line 53 of file EcalSimHitsValidation.h.

Referenced by analyze(), and EcalSimHitsValidation().

EEHitsCollectionToken

edm::EDGetTokenT<edm::PCaloHitContainer> EcalSimHitsValidation::EEHitsCollectionToken

private

Definition at line 54 of file EcalSimHitsValidation.h.

Referenced by analyze(), and EcalSimHitsValidation().

ESHitsCollectionToken

edm::EDGetTokenT<edm::PCaloHitContainer> EcalSimHitsValidation::ESHitsCollectionToken

private

Definition at line 55 of file EcalSimHitsValidation.h.

Referenced by analyze(), and EcalSimHitsValidation().

g4InfoLabel

std::string EcalSimHitsValidation::g4InfoLabel

private

Definition at line 51 of file EcalSimHitsValidation.h.

Referenced by EcalSimHitsValidation().

HepMCToken

edm::EDGetTokenT<edm::HepMCProduct> EcalSimHitsValidation::HepMCToken

private

Definition at line 52 of file EcalSimHitsValidation.h.

Referenced by analyze().

meEBEnergyFraction_

MonitorElement* EcalSimHitsValidation::meEBEnergyFraction_

private

Definition at line 63 of file EcalSimHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meEEEnergyFraction_

MonitorElement* EcalSimHitsValidation::meEEEnergyFraction_

private

Definition at line 64 of file EcalSimHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meESEnergyFraction_

MonitorElement* EcalSimHitsValidation::meESEnergyFraction_

private

Definition at line 65 of file EcalSimHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meGunEnergy_

MonitorElement* EcalSimHitsValidation::meGunEnergy_

private

Definition at line 59 of file EcalSimHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meGunEta_

MonitorElement* EcalSimHitsValidation::meGunEta_

private

Definition at line 60 of file EcalSimHitsValidation.h.

Referenced by analyze(), and bookHistograms().

meGunPhi_

MonitorElement* EcalSimHitsValidation::meGunPhi_

private

Definition at line 61 of file EcalSimHitsValidation.h.

Referenced by analyze(), and bookHistograms().

verbose_

bool EcalSimHitsValidation::verbose_

private

Definition at line 57 of file EcalSimHitsValidation.h.

Referenced by EcalSimHitsValidation().