#include <WValidation.h>

Inheritance diagram for WValidation:

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

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

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

	WValidation (const edm::ParameterSet &)

	~WValidation () 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

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

MonitorElement *	gamma_energy

edm::InputTag	hepmcCollection_

edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_

MonitorElement *	lepeta

MonitorElement *	lepmet_mT

MonitorElement *	lepmet_mTPeak

MonitorElement *	lepmet_pt

MonitorElement *	lepmet_ptLog

MonitorElement *	lepmet_rap

MonitorElement *	leppt

MonitorElement *	met

MonitorElement *	nEvt

MonitorElement *	Wdaughters

WeightManager	wmanager_

MonitorElement *	Wmass

MonitorElement *	WmassPeak

MonitorElement *	Wpt

MonitorElement *	WptLog

MonitorElement *	Wrap

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 33 of file WValidation.h.

Constructor & Destructor Documentation

◆ WValidation()

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

explicit

Definition at line 18 of file WValidation.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_.

◆ ~WValidation()

WValidation::~WValidation ( )

override

Definition at line 26 of file WValidation.cc.

26 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Gathering the HepMCProduct information

Reimplemented from DQMEDAnalyzer.

Definition at line 89 of file WValidation.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);
  
   nEvt->Fill(0.5, weight);
  
   std::vector<const HepMC::GenParticle*> allleptons;
   std::vector<const HepMC::GenParticle*> allneutrinos;
  
   //requires status 1 for leptons and neutrinos (except tau)
   int requiredstatus = (std::abs(_flavor) == 11 || std::abs(_flavor) == 13) ? 1 : 3;
  
   bool vetotau = true;
   // alternatively (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) {
       //@todo: improve this selection
       if ((*iter)->pdg_id() == _flavor)
         allleptons.push_back(*iter);
       else if (std::abs((*iter)->pdg_id()) == std::abs(_flavor) + 1)
         allneutrinos.push_back(*iter);
     }
   }
  
   //nothing to do if we don't have 2 particles
   if (allleptons.empty() || allneutrinos.empty())
     return;
  
   //sort them in pt
   std::sort(allleptons.begin(), allleptons.end(), HepMCValidationHelper::sortByPt);
   std::sort(allneutrinos.begin(), allneutrinos.end(), HepMCValidationHelper::sortByPt);
  
   //get the first lepton and the first neutrino, and check that one is particle one is antiparticle (product of pdgids < 0)
   std::vector<const HepMC::GenParticle*> products;
   if (allleptons.front()->pdg_id() * allneutrinos.front()->pdg_id() > 0)
     return;
  
   //require at least 20 GeV on the lepton
   if (allleptons.front()->momentum().perp() < 20. || allneutrinos.front()->momentum().perp() < 20.)
     return;
  
   //find possible qed fsr photons
   std::vector<const HepMC::GenParticle*> selectedLepton;
   selectedLepton.push_back(allleptons.front());
   std::vector<const HepMC::GenParticle*> fsrphotons;
   HepMCValidationHelper::findFSRPhotons(selectedLepton, myGenEvent, 0.1, fsrphotons);
  
   Wdaughters->Fill(allleptons.front()->pdg_id(), weight);
   Wdaughters->Fill(allneutrinos.front()->pdg_id(), weight);
  
   //assemble FourMomenta
   TLorentzVector lep1(allleptons[0]->momentum().x(),
                       allleptons[0]->momentum().y(),
                       allleptons[0]->momentum().z(),
                       allleptons[0]->momentum().t());
   TLorentzVector lep2(allneutrinos[0]->momentum().x(),
                       allneutrinos[0]->momentum().y(),
                       allneutrinos[0]->momentum().z(),
                       allneutrinos[0]->momentum().t());
   TLorentzVector dilepton_mom = lep1 + lep2;
   TLorentzVector dilepton_andphoton_mom = dilepton_mom;
   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;
     Wdaughters->Fill(fsrphotons[ipho]->pdg_id(), weight);
     gammasMomenta.push_back(phomom);
   }
   //Fill "true" W histograms
   Wmass->Fill(dilepton_andphoton_mom.M(), weight);
   WmassPeak->Fill(dilepton_andphoton_mom.M(), weight);
   Wpt->Fill(dilepton_andphoton_mom.Pt(), weight);
   WptLog->Fill(log10(dilepton_andphoton_mom.Pt()), weight);
   Wrap->Fill(dilepton_andphoton_mom.Rapidity(), weight);
  
   TLorentzVector met_mom = HepMCValidationHelper::genMet(myGenEvent, -3., 3.);
   TLorentzVector lep1T(lep1.Px(), lep1.Py(), 0., lep1.Et());
   TLorentzVector lepmet_mom = lep1T + met_mom;
   //Fill lepmet histograms
   lepmet_mT->Fill(lepmet_mom.M(), weight);
   lepmet_mTPeak->Fill(lepmet_mom.M(), weight);
   lepmet_pt->Fill(lepmet_mom.Pt(), weight);
   lepmet_ptLog->Fill(log10(lepmet_mom.Pt()), weight);
  
   //Fill lepton histograms
   leppt->Fill(lep1.Pt(), weight);
   lepeta->Fill(lep1.Eta(), weight);
   met->Fill(met_mom.Pt(), weight);
  
   //boost everything in the W 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());
     cos_theta_gamma_lepton->Fill(cos(dphi), weight);
   }
  
 }  //analyze

References _flavor, funct::abs(), funct::cos(), cos_theta_gamma_lepton, dqm::impl::MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), gamma_energy, HepMCValidationHelper::genMet(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, iEvent, RazorAnalyzer::lep1, RazorAnalyzer::lep2, lepeta, lepmet_mT, lepmet_mTPeak, lepmet_pt, lepmet_ptLog, leppt, met, nEvt, HiggsValidation_cfi::pdg_id, edm::es::products(), jetUpdater_cfi::sort, HepMCValidationHelper::sortByPt(), submitPVValidationJobs::t, Wdaughters, WeightManager::weight(), mps_merge::weight, wmanager_, Wmass, WmassPeak, Wpt, WptLog, Wrap, x, y, and z.

◆ bookHistograms()

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

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 28 of file WValidation.cc.

                                                                                         {
   std::string folderName = "Generator/W";
   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
   Wmass = dqm.book1dHisto("Wmass", "inv. Mass W", 70, 0, 140, "M_{T}^{W} (GeV)", "Number of Events");
   WmassPeak = dqm.book1dHisto("WmassPeak", "inv. Mass W", 80, 80, 100, "M_{T}^{W} (GeV)", "Number of Events");
   Wpt = dqm.book1dHisto("Wpt", "W pt", 100, 0, 200, "P_{T}^{W} (GeV)", "Number of Events");
   WptLog = dqm.book1dHisto("WptLog", "log(W pt)", 100, 0., 5., "Log_{10}(P_{T}^{W}) (GeV)", "Number of Events");
   Wrap = dqm.book1dHisto("Wrap", "W y", 100, -5, 5, "Y^{W}", "Number of Events");
   Wdaughters = dqm.book1dHisto("Wdaughters", "W daughters", 60, -30, 30, "W daughters (PDG ID)", "Number of Events");
  
   lepmet_mT = dqm.book1dHisto("lepmet_mT",
                               "lepton-met transverse mass",
                               70,
                               0,
                               140,
                               "M_{T}^{Lepton_{T}+E_{T}^{Miss}} (GeV)",
                               "Number of Events");
   lepmet_mTPeak = dqm.book1dHisto("lepmet_mTPeak",
                                   "lepton-met transverse mass",
                                   80,
                                   80,
                                   100,
                                   "M_{T}^{Lepton_{T}+E_{T}^{Miss}} (GeV)",
                                   "Number of Events");
   lepmet_pt = dqm.book1dHisto(
       "lepmet_pt", "lepton-met", 100, 0, 200, "P_{T}^{Lepton_{T}+E_{T}^{Miss}} (GeV)", "Number of Events");
   lepmet_ptLog = dqm.book1dHisto("lepmet_ptLog",
                                  "log(lepton-met pt)",
                                  100,
                                  0.,
                                  5.,
                                  "log_{10}(P_{T}^{Lepton_{T}+E_{T}^{Miss}}) (Log_{10}(GeV))",
                                  "Number of Events");
  
   gamma_energy = dqm.book1dHisto(
       "gamma_energy", "photon energy in W rest frame", 200, 0., 100., "E_{#gamma}^{W rest-frame}", "Number of Events");
   cos_theta_gamma_lepton = dqm.book1dHisto("cos_theta_gamma_lepton",
                                            "cos_theta_gamma_lepton in W rest frame",
                                            200,
                                            -1,
                                            1,
                                            "cos(#theta_{#gamma-lepton}^{W rest-frame})",
                                            "Number of Events");
  
   leppt = dqm.book1dHisto("leadpt", "lepton pt", 200, 0., 200., "P_{t}^{Lead-Lepton} (GeV)", "Number of Events");
   met = dqm.book1dHisto("met", "met", 200, 0., 200., "E_{T}^{Miss} (GeV)", "Number of Events");
   lepeta = dqm.book1dHisto("leadeta", "leading lepton eta", 100, -5., 5., "#eta^{Lead-Lepton}", "Number of Events");
  
   return;
 }

References _name, cos_theta_gamma_lepton, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, gamma_energy, mps_fire::i, lepeta, lepmet_mT, lepmet_mTPeak, lepmet_pt, lepmet_ptLog, leppt, met, nEvt, AlCaHLTBitMon_QueryRunRegistry::string, Wdaughters, Wmass, WmassPeak, Wpt, WptLog, and Wrap.

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 87 of file WValidation.cc.

87 { c.getData(fPDGTable); }

References c, and fPDGTable.

Member Data Documentation

◆ _flavor

int WValidation::_flavor

private

decay flavor

Definition at line 55 of file WValidation.h.

Referenced by analyze().

◆ _name

std::string WValidation::_name

private

decay flavor name

Definition at line 57 of file WValidation.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 * WValidation::cos_theta_gamma_lepton

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ fPDGTable

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

private

PDT table.

Definition at line 46 of file WValidation.h.

Referenced by dqmBeginRun().

◆ gamma_energy

MonitorElement* WValidation::gamma_energy

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ hepmcCollection_

edm::InputTag WValidation::hepmcCollection_

private

Definition at line 43 of file WValidation.h.

Referenced by WValidation().

◆ hepmcCollectionToken_

edm::EDGetTokenT<edm::HepMCProduct> WValidation::hepmcCollectionToken_

private

Definition at line 59 of file WValidation.h.

Referenced by analyze(), and WValidation().

◆ lepeta

MonitorElement * WValidation::lepeta

private

Definition at line 51 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ lepmet_mT

MonitorElement* WValidation::lepmet_mT

private

Definition at line 50 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ lepmet_mTPeak

MonitorElement * WValidation::lepmet_mTPeak

private

Definition at line 50 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ lepmet_pt

MonitorElement * WValidation::lepmet_pt

private

Definition at line 50 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ lepmet_ptLog

MonitorElement * WValidation::lepmet_ptLog

private

Definition at line 50 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ lepmet_rap

MonitorElement * WValidation::lepmet_rap

private

Definition at line 50 of file WValidation.h.

◆ leppt

MonitorElement* WValidation::leppt

private

Definition at line 51 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ met

MonitorElement * WValidation::met

private

Definition at line 51 of file WValidation.h.

Referenced by analyze(), objects.METAnalyzer.METAnalyzer::applyDeltaMet(), and bookHistograms().

◆ nEvt

MonitorElement* WValidation::nEvt

private

Definition at line 48 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Wdaughters

MonitorElement * WValidation::Wdaughters

private

Definition at line 49 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ wmanager_

WeightManager WValidation::wmanager_

private

Definition at line 42 of file WValidation.h.

Referenced by analyze().

◆ Wmass

MonitorElement* WValidation::Wmass

private

Definition at line 49 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ WmassPeak

MonitorElement * WValidation::WmassPeak

private

Definition at line 49 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Wpt

MonitorElement * WValidation::Wpt

private

Definition at line 49 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ WptLog

MonitorElement * WValidation::WptLog

private

Definition at line 49 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

◆ Wrap

MonitorElement * WValidation::Wrap

private

Definition at line 49 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ WValidation()

◆ ~WValidation()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

Member Data Documentation

◆ _flavor

◆ _name

◆ cos_theta_gamma_lepton

◆ fPDGTable

◆ gamma_energy

◆ hepmcCollection_

◆ hepmcCollectionToken_

◆ lepeta

◆ lepmet_mT

◆ lepmet_mTPeak

◆ lepmet_pt

◆ lepmet_ptLog

◆ lepmet_rap

◆ leppt

◆ met

◆ nEvt

◆ Wdaughters

◆ wmanager_

◆ Wmass

◆ WmassPeak

◆ Wpt

◆ WptLog

◆ Wrap