#include <DrellYanValidation.h>

Inheritance diagram for DrellYanValidation:

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

	DrellYanValidation (const edm::ParameterSet &)

	~DrellYanValidation () 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 Attributes
int	_flavor
	decay flavor More...

std::string	_name
	decay flavor name More...

MonitorElement *	cos_theta_gamma_lepton

MonitorElement *	dilep_mass

MonitorElement *	dilep_massPeak

MonitorElement *	dilep_pt

MonitorElement *	dilep_ptLog

MonitorElement *	dilep_rap

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

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

MonitorElement *	gamma_energy

edm::InputTag	hepmcCollection_

edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_

MonitorElement *	leadeta

MonitorElement *	leadpt

MonitorElement *	nEvt

MonitorElement *	seceta

MonitorElement *	secpt

WeightManager	wmanager_

MonitorElement *	Zdaughters

MonitorElement *	Zmass

MonitorElement *	ZmassPeak

MonitorElement *	Zpt

MonitorElement *	ZptLog

MonitorElement *	Zrap

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 32 of file DrellYanValidation.h.

Constructor & Destructor Documentation

◆ DrellYanValidation()

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

explicit

Definition at line 19 of file DrellYanValidation.cc.

References fPDGTableToken, hepmcCollection_, and hepmcCollectionToken_.

     : wmanager_(iPSet, consumesCollector()),
       hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")),
       _flavor(iPSet.getParameter<int>("decaysTo")),
       _name(iPSet.getParameter<std::string>("name")) {
   hepmcCollectionToken_ = consumes<HepMCProduct>(hepmcCollection_);
   fPDGTableToken = esConsumes<edm::Transition::BeginRun>();
 }

◆ ~DrellYanValidation()

DrellYanValidation::~DrellYanValidation ( )

override

Definition at line 28 of file DrellYanValidation.cc.

28 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Gathering the HepMCProduct information

Reimplemented from DQMEDAnalyzer.

Definition at line 84 of file DrellYanValidation.cc.

References _flavor, funct::abs(), cms::cuda::assert(), funct::cos(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, dqm::impl::MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), fPDGTable, gamma_energy, edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, mps_fire::i, iEvent, leadeta, leadpt, RazorAnalyzer::lep1, RazorAnalyzer::lep2, nEvt, LHEGenericFilter_cfi::ParticleID, HiggsValidation_cfi::pdg_id, perp(), edm::es::products(), seceta, secpt, jetUpdater_cfi::sort, HepMCValidationHelper::sortByPt(), submitPVValidationJobs::t, WeightManager::weight(), mps_merge::weight, wmanager_, x, y, z, Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.

                                                                                     {
   // we *DO NOT* rely on a Z entry in the particle listings!
 
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(hepmcCollectionToken_, evt);
 
   //Get EVENT
   const HepMC::GenEvent* myGenEvent = evt->GetEvent();
 
   double weight = wmanager_.weight(iEvent);
 
   //std::cout << "weight: " << weight << std::endl;
 
   nEvt->Fill(0.5, weight);
 
   std::vector<const HepMC::GenParticle*> allproducts;
 
   //requires status 1 for leptons and neutrinos (except tau)
   int requiredstatus = (std::abs(_flavor) == 11 || std::abs(_flavor) == 12 || std::abs(_flavor) == 13 ||
                         std::abs(_flavor) == 14 || std::abs(_flavor) == 16)
                            ? 1
                            : 3;
 
   bool vetotau =
       true;  //(std::abs(_flavor) == 11 || std::abs(_flavor) == 12 || std::abs(_flavor) ==13 || std::abs(_flavor) ==14 || std::abs(_flavor) ==16) ? true : false;
 
   for (HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin();
        iter != myGenEvent->particles_end();
        ++iter) {
     if (vetotau) {
       if ((*iter)->status() == 3 && std::abs((*iter)->pdg_id()) == 15)
         return;
     }
     if ((*iter)->status() == requiredstatus) {
       if (std::abs((*iter)->pdg_id()) == _flavor)
         allproducts.push_back(*iter);
     }
   }
 
   //nothing to do if we don't have 2 particles
   if (allproducts.size() < 2)
     return;
 
   //sort them in pt
   std::sort(allproducts.begin(), allproducts.end(), HepMCValidationHelper::sortByPt);
 
   //get the first element and the first following element with opposite charge
   std::vector<const HepMC::GenParticle*> products;
   products.push_back(allproducts.front());
   const HepPDT::ParticleData* PData1 = fPDGTable->particle(HepPDT::ParticleID(allproducts.front()->pdg_id()));
   double charge1 = PData1->charge();
   for (unsigned int i = 1; i < allproducts.size(); ++i) {
     const HepPDT::ParticleData* PData2 = fPDGTable->particle(HepPDT::ParticleID(allproducts[i]->pdg_id()));
     double charge2 = PData2->charge();
     if (charge1 * charge2 < 0)
       products.push_back(allproducts[i]);
   }
 
   //if we did not find any opposite charge pair there is nothing to do
   if (products.size() < 2)
     return;
 
   assert(products[0]->momentum().perp() >= products[1]->momentum().perp());
 
   //leading lepton with pt > 20.
   if (products[0]->momentum().perp() < 20.)
     return;
 
   //assemble FourMomenta
   TLorentzVector lep1(products[0]->momentum().x(),
                       products[0]->momentum().y(),
                       products[0]->momentum().z(),
                       products[0]->momentum().t());
   TLorentzVector lep2(products[1]->momentum().x(),
                       products[1]->momentum().y(),
                       products[1]->momentum().z(),
                       products[1]->momentum().t());
   TLorentzVector dilepton_mom = lep1 + lep2;
   TLorentzVector dilepton_andphoton_mom = dilepton_mom;
 
   //mass > 60.
   if (dilepton_mom.M() < 60.)
     return;
 
   //find possible qed fsr photons
   std::vector<const HepMC::GenParticle*> fsrphotons;
   HepMCValidationHelper::findFSRPhotons(products, myGenEvent, 0.1, fsrphotons);
 
   Zdaughters->Fill(products[0]->pdg_id(), weight);
   Zdaughters->Fill(products[1]->pdg_id(), weight);
 
   std::vector<TLorentzVector> gammasMomenta;
   for (unsigned int ipho = 0; ipho < fsrphotons.size(); ++ipho) {
     TLorentzVector phomom(fsrphotons[ipho]->momentum().x(),
                           fsrphotons[ipho]->momentum().y(),
                           fsrphotons[ipho]->momentum().z(),
                           fsrphotons[ipho]->momentum().t());
     dilepton_andphoton_mom += phomom;
     Zdaughters->Fill(fsrphotons[ipho]->pdg_id(), weight);
     gammasMomenta.push_back(phomom);
   }
   //Fill Z histograms
   Zmass->Fill(dilepton_andphoton_mom.M(), weight);
   ZmassPeak->Fill(dilepton_andphoton_mom.M(), weight);
   Zpt->Fill(dilepton_andphoton_mom.Pt(), weight);
   ZptLog->Fill(log10(dilepton_andphoton_mom.Pt()), weight);
   Zrap->Fill(dilepton_andphoton_mom.Rapidity(), weight);
 
   //Fill dilepton histograms
   dilep_mass->Fill(dilepton_mom.M(), weight);
   dilep_massPeak->Fill(dilepton_mom.M(), weight);
   dilep_pt->Fill(dilepton_mom.Pt(), weight);
   dilep_ptLog->Fill(log10(dilepton_mom.Pt()), weight);
   dilep_rap->Fill(dilepton_mom.Rapidity(), weight);
 
   //Fill lepton histograms
   leadpt->Fill(lep1.Pt(), weight);
   secpt->Fill(lep2.Pt(), weight);
   leadeta->Fill(lep1.Eta(), weight);
   seceta->Fill(lep2.Eta(), weight);
 
   //boost everything in the Z frame
   TVector3 boost = dilepton_andphoton_mom.BoostVector();
   boost *= -1.;
   lep1.Boost(boost);
   lep2.Boost(boost);
   for (unsigned int ipho = 0; ipho < gammasMomenta.size(); ++ipho) {
     gammasMomenta[ipho].Boost(boost);
   }
   std::sort(gammasMomenta.begin(), gammasMomenta.end(), HepMCValidationHelper::GreaterByE<TLorentzVector>);
 
   //fill gamma histograms
   if (!gammasMomenta.empty() && dilepton_andphoton_mom.M() > 50.) {
     gamma_energy->Fill(gammasMomenta.front().E(), weight);
     double dphi = lep1.DeltaR(gammasMomenta.front()) < lep2.DeltaR(gammasMomenta.front())
                       ? lep1.DeltaPhi(gammasMomenta.front())
                       : lep2.DeltaPhi(gammasMomenta.front());
     cos_theta_gamma_lepton->Fill(cos(dphi), weight);
   }
 
 }  //analyze

◆ bookHistograms()

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

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 34 of file DrellYanValidation.cc.

References _name, cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, ctppsRandomDQMSource_cfi::folderName, gamma_energy, mps_fire::i, leadeta, leadpt, nEvt, seceta, secpt, AlCaHLTBitMon_QueryRunRegistry::string, Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.

                                                                                                {
   std::string folderName = "Generator/DrellYan";
   folderName += _name;
   DQMHelper dqm(&i);
   i.setCurrentFolder(folderName);
 
   // Number of analyzed events
   nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1., "bin", "Number of Events");
 
   //Kinematics
   Zmass = dqm.book1dHisto("Zmass", "inv. Mass Z", 70, 0, 140, "M_{Z} (GeV)", "Number of Events");
   ZmassPeak = dqm.book1dHisto("ZmassPeak", "inv. Mass Z", 80, 80, 100, "M_{Z} (GeV)", "Number of Events");
   Zpt = dqm.book1dHisto("Zpt", "Z pt", 100, 0, 200, "P_{t}^{Z} (GeV)", "Number of Events");
   ZptLog =
       dqm.book1dHisto("ZptLog", "log(Z pt)", 100, 0., 5., "log_{10}(P_{t}^{Z}) (log_{10}(GeV))", "Number of Events");
   Zrap = dqm.book1dHisto("Zrap", "Z y", 100, -5, 5, "Y_{Z}", "Number of Events");
   Zdaughters = dqm.book1dHisto("Zdaughters", "Z daughters", 60, -30, 30, "Z daughters (PDG ID)", "Number of Events");
 
   dilep_mass = dqm.book1dHisto("dilep_mass", "inv. Mass dilepton", 70, 0, 140, "M_{ll} (GeV)", "Number of Events");
   dilep_massPeak =
       dqm.book1dHisto("dilep_massPeak", "inv. Mass dilepton", 80, 80, 100, "M_{ll} (GeV)", "Number of Events");
   dilep_pt = dqm.book1dHisto("dilep_pt", "dilepton pt", 100, 0, 200, "P_{t}^{ll} (GeV)", "Number of Events");
   dilep_ptLog = dqm.book1dHisto(
       "dilep_ptLog", "log(dilepton pt)", 100, 0., 5., "log_{10}(P_{t}^{ll}) (log_{10}(GeV))", "Number of Events");
   dilep_rap = dqm.book1dHisto("dilep_rap", "dilepton y", 100, -5, 5, "Y_{ll}", "Number of Events");
 
   gamma_energy = dqm.book1dHisto("gamma_energy",
                                  "photon energy in Z rest frame",
                                  200,
                                  0.,
                                  100.,
                                  "E_{#gamma}^{Z rest-frame} (GeV)",
                                  "Number of Events");
   cos_theta_gamma_lepton = dqm.book1dHisto("cos_theta_gamma_lepton",
                                            "cos_theta_gamma_lepton in Z rest frame",
                                            200,
                                            -1,
                                            1,
                                            "cos(#theta_{#gamma-lepton}^{Z rest-frame})",
                                            "Number of Events");
 
   leadpt = dqm.book1dHisto("leadpt", "leading lepton pt", 200, 0., 200., "P_{t}^{1st-lepton}", "Number of Events");
   secpt = dqm.book1dHisto("secpt", "second lepton pt", 200, 0., 200., "P_{t}^{2nd-lepton}", "Number of Events");
   leadeta = dqm.book1dHisto("leadeta", "leading lepton eta", 100, -5., 5., "#eta^{1st-lepton}", "Number of Events");
   seceta = dqm.book1dHisto("seceta", "second lepton eta", 100, -5., 5., "#eta^{2nd-lepton}", "Number of Events");
 
   return;
 }

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 30 of file DrellYanValidation.cc.

References HltBtagPostValidation_cff::c, fPDGTable, and fPDGTableToken.

                                                                             {
   fPDGTable = c.getHandle(fPDGTableToken);
 }

Member Data Documentation

◆ _flavor

int DrellYanValidation::_flavor

private

decay flavor

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze().

◆ _name

std::string DrellYanValidation::_name

private

decay flavor name

Definition at line 58 of file DrellYanValidation.h.

Referenced by plotting.Subtract::__str__(), plotting.Transform::__str__(), plotting.FakeDuplicate::__str__(), plotting.CutEfficiency::__str__(), plotting.AggregateBins::__str__(), plotting.AggregateHistos::__str__(), plotting.ROC::__str__(), trackingPlots.TimePerEventPlot::__str__(), trackingPlots.TimePerTrackPlot::__str__(), trackingPlots.TimePerEventPlot::_create(), plotting.PlotGroup::_drawSeparate(), plotting.PlotGroup::_save(), bookHistograms(), plotting.Subtract::create(), plotting.Transform::create(), plotting.FakeDuplicate::create(), plotting.CutEfficiency::create(), plotting.AggregateBins::create(), plotting.AggregateHistos::create(), trackingPlots.TimePerTrackPlot::create(), plotting.Plot::create(), plotting.PlotGroup::draw(), plotting.Plot::getName(), plotting.PlotGroup::getName(), plotting.PlotterFolder::getName(), plotting.PlotterItem::getName(), plotting.PlotterFolder::getSelectionNameIterator(), validation.SimpleSample::name(), trackingPlots.Iteration::name(), plotting.PlotterItem::readDirs(), plotting.PlotGroup::remove(), and plotting.Plot::setName().

◆ cos_theta_gamma_lepton

MonitorElement * DrellYanValidation::cos_theta_gamma_lepton

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_mass

MonitorElement* DrellYanValidation::dilep_mass

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_massPeak

MonitorElement * DrellYanValidation::dilep_massPeak

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_pt

MonitorElement * DrellYanValidation::dilep_pt

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_ptLog

MonitorElement * DrellYanValidation::dilep_ptLog

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_rap

MonitorElement * DrellYanValidation::dilep_rap

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ fPDGTable

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

private

PDT table.

Definition at line 46 of file DrellYanValidation.h.

Referenced by analyze(), and dqmBeginRun().

◆ fPDGTableToken

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

private

Definition at line 47 of file DrellYanValidation.h.

Referenced by dqmBeginRun(), and DrellYanValidation().

◆ gamma_energy

MonitorElement* DrellYanValidation::gamma_energy

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ hepmcCollection_

edm::InputTag DrellYanValidation::hepmcCollection_

private

Definition at line 43 of file DrellYanValidation.h.

Referenced by DrellYanValidation().

◆ hepmcCollectionToken_

edm::EDGetTokenT<edm::HepMCProduct> DrellYanValidation::hepmcCollectionToken_

private

Definition at line 60 of file DrellYanValidation.h.

Referenced by analyze(), and DrellYanValidation().

◆ leadeta

MonitorElement * DrellYanValidation::leadeta

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ leadpt

MonitorElement* DrellYanValidation::leadpt

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ nEvt

MonitorElement* DrellYanValidation::nEvt

private

Definition at line 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ seceta

MonitorElement * DrellYanValidation::seceta

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ secpt

MonitorElement * DrellYanValidation::secpt

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ wmanager_

WeightManager DrellYanValidation::wmanager_

private

Definition at line 42 of file DrellYanValidation.h.

Referenced by analyze().

◆ Zdaughters

MonitorElement * DrellYanValidation::Zdaughters

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Zmass

MonitorElement* DrellYanValidation::Zmass

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ ZmassPeak

MonitorElement * DrellYanValidation::ZmassPeak

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Zpt

MonitorElement * DrellYanValidation::Zpt

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ ZptLog

MonitorElement * DrellYanValidation::ZptLog

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Zrap

MonitorElement * DrellYanValidation::Zrap

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DrellYanValidation()

◆ ~DrellYanValidation()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

Member Data Documentation

◆ _flavor

◆ _name

◆ cos_theta_gamma_lepton

◆ dilep_mass

◆ dilep_massPeak

◆ dilep_pt

◆ dilep_ptLog

◆ dilep_rap

◆ fPDGTable

◆ fPDGTableToken

◆ gamma_energy

◆ hepmcCollection_

◆ hepmcCollectionToken_

◆ leadeta

◆ leadpt

◆ nEvt

◆ seceta

◆ secpt

◆ wmanager_

◆ Zdaughters

◆ Zmass

◆ ZmassPeak

◆ Zpt

◆ ZptLog

◆ Zrap