#include <HiggsValidation.h>

Inheritance diagram for HiggsValidation:

Classes
class	MonitoredDecays

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override

	HiggsValidation (const edm::ParameterSet &)

	~HiggsValidation () override

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

	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

Private Member Functions
std::string	convert (int)

int	findHiggsDecayChannel (const HepMC::GenParticle , std::vector< HepMC::GenParticle > &decayprod)

Private Attributes
edm::ESHandle< HepPDT::ParticleDataTable >	fPDGTable
	PDT table. More...

edm::InputTag	hepmcCollection_

edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_

MonitorElement *	Higgs_eta

MonitorElement *	Higgs_mass

MonitorElement *	Higgs_pt

MonitorElement *	HiggsDecayChannels

std::vector< MonitorElement * >	HiggsDecayProd_eta

std::vector< MonitorElement * >	HiggsDecayProd_pt

MonitoredDecays *	monitoredDecays

MonitorElement *	nEvt

int	particle_id

std::string	particle_name

WeightManager	wmanager_

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 Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 36 of file HiggsValidation.h.

Constructor & Destructor Documentation

◆ HiggsValidation()

HiggsValidation::HiggsValidation ( const edm::ParameterSet & iPSet )

explicit

Definition at line 21 of file HiggsValidation.cc.

     : wmanager_(iPSet, consumesCollector()),
       hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")),
       particle_id(iPSet.getParameter<int>("pdg_id")),
       particle_name(iPSet.getParameter<std::string>("particleName")) {
   monitoredDecays = new MonitoredDecays(iPSet);
   hepmcCollectionToken_ = consumes<HepMCProduct>(hepmcCollection_);
 }

References hepmcCollection_, hepmcCollectionToken_, and monitoredDecays.

◆ ~HiggsValidation()

HiggsValidation::~HiggsValidation ( )

override

Definition at line 30 of file HiggsValidation.cc.

30 {}

Member Function Documentation

◆ analyze()

void HiggsValidation::analyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 99 of file HiggsValidation.cc.

                                                                                  {
   double weight = wmanager_.weight(iEvent);
   nEvt->Fill(0.5, weight);
  
   //Gathering the HepMCProduct information
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(hepmcCollectionToken_, evt);
  
   //Get EVENT
   HepMC::GenEvent* myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
  
   // loop over all particles
   bool filled = false;
   for (HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin();
        iter != myGenEvent->particles_end() && !filled;
        ++iter) {
     if (particle_id == std::abs((*iter)->pdg_id())) {
       std::vector<HepMC::GenParticle*> decayprod;
       int channel = findHiggsDecayChannel(*iter, decayprod);
       HiggsDecayChannels->Fill(channel, weight);
       Higgs_pt->Fill((*iter)->momentum().perp(), weight);
       Higgs_eta->Fill((*iter)->momentum().eta(), weight);
       Higgs_mass->Fill((*iter)->momentum().m(), weight);
       for (unsigned int i = 0; i < decayprod.size(); i++) {
         int idx = monitoredDecays->isDecayParticle(decayprod.at(i)->pdg_id());
         if (0 <= idx && idx <= (int)HiggsDecayProd_pt.size()) {
           HiggsDecayProd_pt.at(idx)->Fill(decayprod.at(i)->momentum().perp(), weight);
           HiggsDecayProd_eta.at(idx)->Fill(decayprod.at(i)->momentum().eta(), weight);
         }
       }
       filled = true;
     }
   }
  
   delete myGenEvent;
  
 }  //analyze

References funct::abs(), dqm::impl::MonitorElement::Fill(), findHiggsDecayChannel(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, Higgs_eta, Higgs_mass, Higgs_pt, HiggsDecayChannels, HiggsDecayProd_eta, HiggsDecayProd_pt, mps_fire::i, training_settings::idx, iEvent, HiggsValidation::MonitoredDecays::isDecayParticle(), monitoredDecays, nEvt, particle_id, WeightManager::weight(), mps_merge::weight, and wmanager_.

◆ bookHistograms()

void HiggsValidation::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 34 of file HiggsValidation.cc.

                                                                                             {
   std::string dir = "Generator/";
   dir += particle_name;
   DQMHelper dqm(&i);
   i.setCurrentFolder(dir);
  
   // Number of analyzed events
   nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1., "bin", "Number of Events");
  
   //decay type
  
   std::string channel = particle_name + "_DecayChannels";
   HiggsDecayChannels = dqm.book1dHisto(channel,
                                        particle_name + " decay channels",
                                        monitoredDecays->size(),
                                        0,
                                        monitoredDecays->size(),
                                        "Decay Channels",
                                        "Number of Events");
  
   for (size_t j = 0; j < monitoredDecays->size(); ++j) {
     HiggsDecayChannels->setBinLabel(1 + j, monitoredDecays->channel(j));
   }
  
   //Kinematics
   Higgs_pt = dqm.book1dHisto((particle_name + "_pt"),
                              (particle_name + " p_{t}"),
                              50,
                              0,
                              250,
                              "P_{t}^{" + particle_name + "} (GeV)",
                              "Number of Events");
   Higgs_eta = dqm.book1dHisto((particle_name + "_eta"),
                               (particle_name + " #eta"),
                               50,
                               -5,
                               5,
                               "#eta^{" + particle_name + "}",
                               "Number of Events");
   Higgs_mass = dqm.book1dHisto(
       (particle_name + "_m"), (particle_name + " M"), 500, 0, 500, "M^{" + particle_name + "}", "Number of Events");
  
   int idx = 0;
   for (unsigned int j = 0; j < monitoredDecays->NDecayParticles(); j++) {
     HiggsDecayProd_pt.push_back(dqm.book1dHisto((monitoredDecays->ConvertIndex(idx) + "_pt"),
                                                 (monitoredDecays->ConvertIndex(idx) + " p_{t}"),
                                                 50,
                                                 0,
                                                 250,
                                                 "P_{t}^{" + monitoredDecays->ConvertIndex(idx) + "} (GeV)",
                                                 "Number of Events"));
     HiggsDecayProd_eta.push_back(dqm.book1dHisto((monitoredDecays->ConvertIndex(idx) + "_eta"),
                                                  (monitoredDecays->ConvertIndex(idx) + " #eta"),
                                                  50,
                                                  -5,
                                                  5,
                                                  "#eta^{" + monitoredDecays->ConvertIndex(idx) + "}",
                                                  "Number of Events"));
     idx++;
   }
  
   return;
 }

References HiggsValidation::MonitoredDecays::channel(), HiggsValidation::MonitoredDecays::ConvertIndex(), DeadROC_duringRun::dir, Higgs_eta, Higgs_mass, Higgs_pt, HiggsDecayChannels, HiggsDecayProd_eta, HiggsDecayProd_pt, mps_fire::i, training_settings::idx, dqmiolumiharvest::j, monitoredDecays, HiggsValidation::MonitoredDecays::NDecayParticles(), nEvt, particle_name, dqm::impl::MonitorElement::setBinLabel(), HiggsValidation::MonitoredDecays::size(), and AlCaHLTBitMon_QueryRunRegistry::string.

◆ convert()

std::string HiggsValidation::convert ( int )

private

◆ dqmBeginRun()

void HiggsValidation::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 32 of file HiggsValidation.cc.

32 { c.getData(fPDGTable); }

References HltBtagPostValidation_cff::c, and fPDGTable.

◆ findHiggsDecayChannel()

int HiggsValidation::findHiggsDecayChannel	(	const HepMC::GenParticle *	genParticle,
		std::vector< HepMC::GenParticle * > &	decayprod
	)

private

Definition at line 137 of file HiggsValidation.cc.

                                                                                     {
   if (genParticle->status() == 1)
     return monitoredDecays->stable();
   std::vector<int> children;
   if (genParticle->end_vertex()) {
     HepMC::GenVertex::particle_iterator des;
     for (des = genParticle->end_vertex()->particles_begin(HepMC::descendants);
          des != genParticle->end_vertex()->particles_end(HepMC::descendants);
          ++des) {
       if ((*des)->pdg_id() == genParticle->pdg_id())
         continue;
  
       HepMC::GenVertex::particle_iterator mother = (*des)->production_vertex()->particles_begin(HepMC::parents);
       if ((*mother)->pdg_id() == genParticle->pdg_id()) {
         children.push_back((*des)->pdg_id());
         decayprod.push_back((*des));
       }
     }
   }
  
   if (children.size() == 2 && children.at(0) != 0 && children.at(1) != 0) {
     return monitoredDecays->position(children.at(0), children.at(1));
   }
   return monitoredDecays->undetermined();
 }

References class-composition::children, monitoredDecays, parents, HiggsValidation::MonitoredDecays::position(), HiggsValidation::MonitoredDecays::stable(), and HiggsValidation::MonitoredDecays::undetermined().

Referenced by analyze().

Member Data Documentation

◆ fPDGTable

edm::ESHandle<HepPDT::ParticleDataTable> HiggsValidation::fPDGTable

private

PDT table.

Definition at line 172 of file HiggsValidation.h.

Referenced by dqmBeginRun().

◆ hepmcCollection_

edm::InputTag HiggsValidation::hepmcCollection_

private

Definition at line 164 of file HiggsValidation.h.

Referenced by HiggsValidation().

◆ hepmcCollectionToken_

edm::EDGetTokenT<edm::HepMCProduct> HiggsValidation::hepmcCollectionToken_

private

Definition at line 184 of file HiggsValidation.h.

Referenced by analyze(), and HiggsValidation().

◆ Higgs_eta

MonitorElement* HiggsValidation::Higgs_eta

private

Definition at line 178 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Higgs_mass

MonitorElement* HiggsValidation::Higgs_mass

private

Definition at line 179 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Higgs_pt

MonitorElement* HiggsValidation::Higgs_pt

private

Definition at line 177 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ HiggsDecayChannels

MonitorElement* HiggsValidation::HiggsDecayChannels

private

Definition at line 175 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ HiggsDecayProd_eta

std::vector<MonitorElement *> HiggsValidation::HiggsDecayProd_eta

private

Definition at line 182 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ HiggsDecayProd_pt

std::vector<MonitorElement *> HiggsValidation::HiggsDecayProd_pt

private

Definition at line 181 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ monitoredDecays

MonitoredDecays* HiggsValidation::monitoredDecays

private

Definition at line 169 of file HiggsValidation.h.

Referenced by analyze(), bookHistograms(), findHiggsDecayChannel(), and HiggsValidation().

◆ nEvt

MonitorElement* HiggsValidation::nEvt

private

Definition at line 174 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

◆ particle_id

int HiggsValidation::particle_id

private

Definition at line 166 of file HiggsValidation.h.

Referenced by analyze().

◆ particle_name

std::string HiggsValidation::particle_name

private

Definition at line 167 of file HiggsValidation.h.

Referenced by bookHistograms().

◆ wmanager_

WeightManager HiggsValidation::wmanager_

private

Definition at line 46 of file HiggsValidation.h.

Referenced by analyze().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HiggsValidation()

◆ ~HiggsValidation()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ convert()

◆ dqmBeginRun()

◆ findHiggsDecayChannel()

Member Data Documentation

◆ fPDGTable

◆ hepmcCollection_

◆ hepmcCollectionToken_

◆ Higgs_eta

◆ Higgs_mass

◆ Higgs_pt

◆ HiggsDecayChannels

◆ HiggsDecayProd_eta

◆ HiggsDecayProd_pt

◆ monitoredDecays

◆ nEvt

◆ particle_id

◆ particle_name

◆ wmanager_