HiggsValidation Class Reference

#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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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::ESGetToken< HepPDT::ParticleDataTable, edm::DefaultRecord >	fPDGTableToken
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 >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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.

References fPDGTableToken, hepmcCollection_, hepmcCollectionToken_, and monitoredDecays.

    : 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_);
  fPDGTableToken = esConsumes<edm::Transition::BeginRun>();
}

~HiggsValidation()

HiggsValidation::~HiggsValidation ( )

override

Definition at line 31 of file HiggsValidation.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 102 of file HiggsValidation.cc.

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, heavyIonCSV_trainingSettings::idx, iEvent, HiggsValidation::MonitoredDecays::isDecayParticle(), monitoredDecays, nEvt, particle_id, WeightManager::weight(), mps_merge::weight, and wmanager_.

                                                                                 {
  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

bookHistograms()

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

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 37 of file HiggsValidation.cc.

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

                                                                                            {
  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;
}

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 33 of file HiggsValidation.cc.

References c, fPDGTable, and fPDGTableToken.

                                                                         {
  fPDGTable = c.getHandle(fPDGTableToken);
}

findHiggsDecayChannel()

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

private

Definition at line 140 of file HiggsValidation.cc.

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

Referenced by analyze().

                                                                                    {
  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();
}

Member Data Documentation

fPDGTable

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

private

PDT table.

Definition at line 172 of file HiggsValidation.h.

Referenced by dqmBeginRun().

fPDGTableToken

edm::ESGetToken<HepPDT::ParticleDataTable, edm::DefaultRecord> HiggsValidation::fPDGTableToken

private

Definition at line 173 of file HiggsValidation.h.

Referenced by dqmBeginRun(), and HiggsValidation().

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 185 of file HiggsValidation.h.

Referenced by analyze(), and HiggsValidation().

Higgs_eta

MonitorElement* HiggsValidation::Higgs_eta

private

Definition at line 179 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

Higgs_mass

MonitorElement* HiggsValidation::Higgs_mass

private

Definition at line 180 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

Higgs_pt

MonitorElement* HiggsValidation::Higgs_pt

private

Definition at line 178 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

HiggsDecayChannels

MonitorElement* HiggsValidation::HiggsDecayChannels

private

Definition at line 176 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

HiggsDecayProd_eta

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

private

Definition at line 183 of file HiggsValidation.h.

Referenced by analyze(), and bookHistograms().

HiggsDecayProd_pt

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

private

Definition at line 182 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 175 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().