#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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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(), HepMCValidationHelper::sortByPt(), OrderedSet::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 HltBtagPostValidation_cff::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(), and plotting.PlotGroup::remove().

◆ 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