#include <WValidation.h>

Inheritance diagram for WValidation:

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

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

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

	WValidation (const edm::ParameterSet &)

virtual	~WValidation ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 35 of file WValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file WValidation.cc.

References 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_);
 }

WValidation::~WValidation ( )

virtual

Definition at line 28 of file WValidation.cc.

28 {}

Member Function Documentation

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

overridevirtual

Gathering the HepMCProduct information

Implements edm::stream::EDAnalyzerBase.

Definition at line 66 of file WValidation.cc.

References _flavor, funct::abs(), funct::cos(), cos_theta_gamma_lepton, MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), gamma_energy, HepMCValidationHelper::genMet(), edm::Event::getByToken(), hepmcCollectionToken_, getDQMSummary::iter, lepeta, lepmet_mT, lepmet_mTPeak, lepmet_pt, lepmet_ptLog, leppt, met, nEvt, edm::es::products(), python.multivaluedict::sort(), HepMCValidationHelper::sortByPt(), edmStreamStallGrapher::t, Wdaughters, WeightManager::weight(), histoStyle::weight, wmanager_, Wmass, WmassPeak, Wpt, WptLog, Wrap, x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

 { 
   
   // 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 = (abs(_flavor) == 11 || abs(_flavor) ==13 ) ? 1 : 3;
 
   bool vetotau = true; //(abs(_flavor) == 11 || abs(_flavor) == 12 || abs(_flavor) ==13 || abs(_flavor) ==14 || 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 && 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 (abs((*iter)->pdg_id()) == 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.size() != 0 && 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

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

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 30 of file WValidation.cc.

References _name, DQMStore::IBooker::book1D(), cos_theta_gamma_lepton, gamma_energy, lepeta, lepmet_mT, lepmet_mTPeak, lepmet_pt, lepmet_ptLog, leppt, met, nEvt, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, Wdaughters, Wmass, WmassPeak, Wpt, WptLog, and Wrap.

                                                                                          {
     std::string folderName = "Generator/W";
     folderName+=_name;
     i.setCurrentFolder(folderName.c_str());
     
     // Number of analyzed events
     nEvt = i.book1D("nEvt", "n analyzed Events", 1, 0., 1.);
     
     //Kinematics
     Wmass = i.book1D("Wmass","inv. Mass W", 70 ,0,140);
     WmassPeak = i.book1D("WmassPeak","inv. Mass W", 80 ,80 ,100);
     Wpt = i.book1D("Wpt","W pt",100,0,200);
     WptLog = i.book1D("WptLog","log(W pt)",100,0.,5.);
     Wrap = i.book1D("Wrap", "W y", 100, -5, 5);
     Wdaughters = i.book1D("Wdaughters", "W daughters", 60, -30, 30);
 
     lepmet_mT = i.book1D("lepmet_mT","lepton-met transverse mass", 70 ,0,140);
     lepmet_mTPeak = i.book1D("lepmet_mTPeak","lepton-met transverse mass", 80 ,80 ,100);
     lepmet_pt = i.book1D("lepmet_pt","lepton-met",100,0,200);
     lepmet_ptLog = i.book1D("lepmet_ptLog","log(lepton-met pt)",100,0.,5.);
 
     gamma_energy = i.book1D("gamma_energy", "photon energy in W rest frame", 200, 0., 100.);
     cos_theta_gamma_lepton = i.book1D("cos_theta_gamma_lepton",      "cos_theta_gamma_lepton in W rest frame",      200, -1, 1);
 
     leppt = i.book1D("leadpt","lepton pt", 200, 0., 200.);    
     met   = i.book1D("met","met", 200, 0., 200.);    
     lepeta = i.book1D("leadeta","leading lepton eta", 100, -5., 5.);
 
   return;
 }

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 62 of file WValidation.cc.

References fPDGTable, and edm::EventSetup::getData().

                                                                      {
   c.getData( fPDGTable );
 }

Member Data Documentation

int WValidation::_flavor

private

decay flavor

Definition at line 58 of file WValidation.h.

Referenced by analyze().

std::string WValidation::_name

private

decay flavor name

Definition at line 60 of file WValidation.h.

Referenced by bookHistograms().

MonitorElement * WValidation::cos_theta_gamma_lepton

private

Definition at line 55 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

PDT table.

Definition at line 49 of file WValidation.h.

Referenced by dqmBeginRun().

MonitorElement* WValidation::gamma_energy

private

Definition at line 55 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag WValidation::hepmcCollection_

private

Definition at line 46 of file WValidation.h.

Referenced by WValidation().

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

private

Definition at line 62 of file WValidation.h.

Referenced by analyze(), and WValidation().

MonitorElement * WValidation::lepeta

private

Definition at line 54 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* WValidation::lepmet_mT

private

Definition at line 53 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::lepmet_mTPeak

private

Definition at line 53 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::lepmet_pt

private

Definition at line 53 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::lepmet_ptLog

private

Definition at line 53 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::lepmet_rap

private

Definition at line 53 of file WValidation.h.

MonitorElement* WValidation::leppt

private

Definition at line 54 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::met

private

Definition at line 54 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* WValidation::nEvt

private

Definition at line 51 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::Wdaughters

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

WeightManager WValidation::wmanager_

private

Definition at line 45 of file WValidation.h.

Referenced by analyze().

MonitorElement* WValidation::Wmass

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::WmassPeak

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::Wpt

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::WptLog

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * WValidation::Wrap

private

Definition at line 52 of file WValidation.h.

Referenced by analyze(), and bookHistograms().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation