#include <EcalPreshowerSimHitsValidation.h>

Inheritance diagram for EcalPreshowerSimHitsValidation:

Public Member Functions
	EcalPreshowerSimHitsValidation (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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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
std::string	EEHitsCollection

edm::EDGetTokenT< edm::PCaloHitContainer >	EEHitsToken

std::string	ESHitsCollection

edm::EDGetTokenT< edm::PCaloHitContainer >	ESHitsToken

std::string	g4InfoLabel

std::string	HepMCLabel

edm::EDGetTokenT< edm::HepMCProduct >	HepMCToken

MonitorElement *	me2eszmOver1eszm_

MonitorElement *	me2eszpOver1eszp_

MonitorElement *	meE1alphaE2zm_

MonitorElement *	meE1alphaE2zp_

MonitorElement *	meEEoverESzm_

MonitorElement *	meEEoverESzp_

MonitorElement *	meESEnergyHits1zm_

MonitorElement *	meESEnergyHits1zp_

MonitorElement *	meESEnergyHits2zm_

MonitorElement *	meESEnergyHits2zp_

MonitorElement *	meEShitLog10Energy_

MonitorElement *	meEShitLog10EnergyNorm_

MonitorElement *	menESHits1zm_

MonitorElement *	menESHits1zp_

MonitorElement *	menESHits2zm_

MonitorElement *	menESHits2zp_

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 34 of file EcalPreshowerSimHitsValidation.h.

Member Typedef Documentation

◆ MapType

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

private

Definition at line 35 of file EcalPreshowerSimHitsValidation.h.

Constructor & Destructor Documentation

◆ EcalPreshowerSimHitsValidation()

EcalPreshowerSimHitsValidation::EcalPreshowerSimHitsValidation ( const edm::ParameterSet & ps )

Constructor.

Definition at line 18 of file EcalPreshowerSimHitsValidation.cc.

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

References EEHitsCollection, EEHitsToken, ESHitsCollection, ESHitsToken, g4InfoLabel, edm::ParameterSet::getUntrackedParameter(), HepMCLabel, HepMCToken, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 87 of file EcalPreshowerSimHitsValidation.cc.

                                                                                       {
   edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
  
   edm::Handle<edm::HepMCProduct> MCEvt;
   e.getByToken(HepMCToken, MCEvt);
  
   edm::Handle<edm::PCaloHitContainer> EcalHitsEE;
   e.getByToken(EEHitsToken, EcalHitsEE);
  
   edm::Handle<edm::PCaloHitContainer> EcalHitsES;
   e.getByToken(ESHitsToken, EcalHitsES);
  
   std::vector<PCaloHit> theEECaloHits;
   if (EcalHitsEE.isValid()) {
     theEECaloHits.insert(theEECaloHits.end(), EcalHitsEE->begin(), EcalHitsEE->end());
   }
  
   std::vector<PCaloHit> theESCaloHits;
   if (EcalHitsES.isValid()) {
     theESCaloHits.insert(theESCaloHits.end(), EcalHitsES->begin(), EcalHitsES->end());
   }
  
   double ESEnergy_ = 0.;
   // std::map<unsigned int, std::vector<PCaloHit>,std::less<unsigned int> >
   // CaloHitMap;
  
   // endcap
   double EEetzp_ = 0.;
   double EEetzm_ = 0.;
   for (std::vector<PCaloHit>::iterator isim = theEECaloHits.begin(); isim != theEECaloHits.end(); ++isim) {
     EEDetId eeid(isim->id());
     if (eeid.zside() > 0)
       EEetzp_ += isim->energy();
     if (eeid.zside() < 0)
       EEetzm_ += isim->energy();
   }
  
   uint32_t nESHits1zp = 0;
   uint32_t nESHits1zm = 0;
   uint32_t nESHits2zp = 0;
   uint32_t nESHits2zm = 0;
   double ESet1zp_ = 0.;
   double ESet2zp_ = 0.;
   double ESet1zm_ = 0.;
   double ESet2zm_ = 0.;
   std::vector<double> econtr(140, 0.);
  
   for (std::vector<PCaloHit>::iterator isim = theESCaloHits.begin(); isim != theESCaloHits.end(); ++isim) {
     // CaloHitMap[(*isim).id()].push_back((*isim));
  
     ESDetId esid(isim->id());
  
     LogDebug("HitInfo") << " CaloHit " << isim->getName() << "\n"
                         << " DetID = " << isim->id() << " ESDetId: z side " << esid.zside() << "  plane "
                         << esid.plane() << esid.six() << ',' << esid.siy() << ':' << esid.strip() << "\n"
                         << " Time = " << isim->time() << "\n"
                         << " Track Id = " << isim->geantTrackId() << "\n"
                         << " Energy = " << isim->energy();
  
     ESEnergy_ += isim->energy();
     if (isim->energy() > 0) {
       meEShitLog10Energy_->Fill(log10(isim->energy()));
       int log10i = int((log10(isim->energy()) + 10.) * 10.);
       if (log10i >= 0 && log10i < 140)
         econtr[log10i] += isim->energy();
     }
  
     if (esid.plane() == 1) {
       if (esid.zside() > 0) {
         nESHits1zp++;
         ESet1zp_ += isim->energy();
         meESEnergyHits1zp_->Fill(isim->energy());
       } else if (esid.zside() < 0) {
         nESHits1zm++;
         ESet1zm_ += isim->energy();
         meESEnergyHits1zm_->Fill(isim->energy());
       }
     } else if (esid.plane() == 2) {
       if (esid.zside() > 0) {
         nESHits2zp++;
         ESet2zp_ += isim->energy();
         meESEnergyHits2zp_->Fill(isim->energy());
       } else if (esid.zside() < 0) {
         nESHits2zm++;
         ESet2zm_ += isim->energy();
         meESEnergyHits2zm_->Fill(isim->energy());
       }
     }
   }
  
   menESHits1zp_->Fill(nESHits1zp);
   menESHits1zm_->Fill(nESHits1zm);
  
   menESHits2zp_->Fill(nESHits2zp);
   menESHits2zm_->Fill(nESHits2zm);
  
   if (ESEnergy_ != 0) {
     for (int i = 0; i < 140; i++) {
       meEShitLog10EnergyNorm_->Fill(-10. + (float(i) + 0.5) / 10., econtr[i] / ESEnergy_);
     }
   }
  
   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));
     }
  
     if (heta > 1.653 && heta < 2.6) {
       meE1alphaE2zp_->Fill(ESet1zp_ + 0.7 * ESet2zp_);
       meEEoverESzp_->Fill((ESet1zp_ + 0.7 * ESet2zp_) / 0.00009, EEetzp_);
       if (ESet1zp_ != 0.)
         me2eszpOver1eszp_->Fill(ESet2zp_ / ESet1zp_);
     }
     if (heta < -1.653 && heta > -2.6) {
       meE1alphaE2zm_->Fill(ESet1zm_ + 0.7 * ESet2zm_);
       meEEoverESzm_->Fill((ESet1zm_ + 0.7 * ESet2zm_) / 0.00009, EEetzm_);
       if (ESet1zm_ != 0.)
         me2eszmOver1eszm_->Fill(ESet2zm_ / ESet1zm_);
     }
   }
 }

References MillePedeFileConverter_cfg::e, EEHitsToken, ESHitsToken, dqm::impl::MonitorElement::Fill(), edm::HepMCProduct::GetEvent(), HepMCToken, mps_fire::i, createfilelist::int, edm::HandleBase::isValid(), dqm-mbProfile::log, LogDebug, me2eszmOver1eszm_, me2eszpOver1eszp_, meE1alphaE2zm_, meE1alphaE2zp_, meEEoverESzm_, meEEoverESzp_, meESEnergyHits1zm_, meESEnergyHits1zp_, meESEnergyHits2zm_, meESEnergyHits2zp_, meEShitLog10Energy_, meEShitLog10EnergyNorm_, menESHits1zm_, menESHits1zp_, menESHits2zm_, menESHits2zp_, AlCaHLTBitMon_ParallelJobs::p, and funct::tan().

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 34 of file EcalPreshowerSimHitsValidation.cc.

                                                                                                                {
   ib.setCurrentFolder("EcalHitsV/EcalSimHitsValidation");
   ib.setScope(MonitorElementData::Scope::RUN);
  
   std::string histo = "ES hits layer 1 multiplicity z+";
   menESHits1zp_ = ib.book1D(histo, histo, 50, 0., 50.);
  
   histo = "ES hits layer 2 multiplicity z+";
   menESHits2zp_ = ib.book1D(histo, histo, 50, 0., 50.);
  
   histo = "ES hits layer 1 multiplicity z-";
   menESHits1zm_ = ib.book1D(histo, histo, 50, 0., 50.);
  
   histo = "ES hits layer 2 multiplicity z-";
   menESHits2zm_ = ib.book1D(histo, histo, 50, 0., 50.);
  
   histo = "ES hits energy layer 1 z+";
   meESEnergyHits1zp_ = ib.book1D(histo, histo, 100, 0., 0.05);
  
   histo = "ES hits energy layer 2 z+";
   meESEnergyHits2zp_ = ib.book1D(histo, histo, 100, 0., 0.05);
  
   histo = "ES hits energy layer 1 z-";
   meESEnergyHits1zm_ = ib.book1D(histo, histo, 100, 0., 0.05);
  
   histo = "ES hits energy layer 2 z-";
   meESEnergyHits2zm_ = ib.book1D(histo, histo, 100, 0., 0.05);
  
   histo = "ES hits log10energy spectrum";
   meEShitLog10Energy_ = ib.book1D(histo, histo, 140, -10., 4.);
  
   histo = "ES hits log10energy spectrum vs normalized energy";
   meEShitLog10EnergyNorm_ = ib.bookProfile(histo, histo, 140, -10., 4., 100, 0., 1.);
  
   histo = "ES E1+07E2 z+";
   meE1alphaE2zp_ = ib.book1D(histo, histo, 100, 0., 0.05);
  
   histo = "ES E1+07E2 z-";
   meE1alphaE2zm_ = ib.book1D(histo, histo, 100, 0., 0.05);
  
   histo = "EE vs ES z+";
   meEEoverESzp_ = ib.bookProfile(histo, histo, 250, 0., 500., 200, 0., 200.);
  
   histo = "EE vs ES z-";
   meEEoverESzm_ = ib.bookProfile(histo, histo, 250, 0., 500., 200, 0., 200.);
  
   histo = "ES ene2oEne1 z+";
   me2eszpOver1eszp_ = ib.book1D(histo, histo, 50, 0., 10.);
  
   histo = "ES ene2oEne1 z-";
   me2eszmOver1eszm_ = ib.book1D(histo, histo, 50, 0., 10.);
 }