#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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...

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

Constructor & Destructor Documentation

◆ DrellYanValidation()

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

explicit

Definition at line 19 of file DrellYanValidation.cc.

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

References hepmcCollection_, and hepmcCollectionToken_.

◆ ~DrellYanValidation()

DrellYanValidation::~DrellYanValidation ( )

override

Definition at line 27 of file DrellYanValidation.cc.

27 {}

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 81 of file DrellYanValidation.cc.

                                                                                     {
   // 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

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.

◆ bookHistograms()

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

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 31 of file DrellYanValidation.cc.

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

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

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 29 of file DrellYanValidation.cc.

29 { c.getData(fPDGTable); }

References c, and fPDGTable.

Member Data Documentation

◆ _flavor

int DrellYanValidation::_flavor

private

decay flavor

Definition at line 55 of file DrellYanValidation.h.

Referenced by analyze().

◆ _name

std::string DrellYanValidation::_name

private

decay flavor name

Definition at line 57 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 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_mass

MonitorElement* DrellYanValidation::dilep_mass

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_massPeak

MonitorElement * DrellYanValidation::dilep_massPeak

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_pt

MonitorElement * DrellYanValidation::dilep_pt

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_ptLog

MonitorElement * DrellYanValidation::dilep_ptLog

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ dilep_rap

MonitorElement * DrellYanValidation::dilep_rap

private

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

◆ gamma_energy

MonitorElement* DrellYanValidation::gamma_energy

private

Definition at line 52 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 59 of file DrellYanValidation.h.

Referenced by analyze(), and DrellYanValidation().

◆ leadeta

MonitorElement * DrellYanValidation::leadeta

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ leadpt

MonitorElement* DrellYanValidation::leadpt

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ nEvt

MonitorElement* DrellYanValidation::nEvt

private

Definition at line 48 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ seceta

MonitorElement * DrellYanValidation::seceta

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ secpt

MonitorElement * DrellYanValidation::secpt

private

Definition at line 51 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 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Zmass

MonitorElement* DrellYanValidation::Zmass

private

Definition at line 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ ZmassPeak

MonitorElement * DrellYanValidation::ZmassPeak

private

Definition at line 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Zpt

MonitorElement * DrellYanValidation::Zpt

private

Definition at line 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ ZptLog

MonitorElement * DrellYanValidation::ZptLog

private

Definition at line 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Zrap

MonitorElement * DrellYanValidation::Zrap

private

Definition at line 49 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

◆ gamma_energy

◆ hepmcCollection_

◆ hepmcCollectionToken_

◆ leadeta

◆ leadpt

◆ nEvt

◆ seceta

◆ secpt

◆ wmanager_

◆ Zdaughters

◆ Zmass

◆ ZmassPeak

◆ Zpt

◆ ZptLog

◆ Zrap