#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

Detailed Description

Definition at line 33 of file DrellYanValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file DrellYanValidation.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_);
 }

DrellYanValidation::~DrellYanValidation ( )

override

Definition at line 28 of file DrellYanValidation.cc.

28 {}

Member Function Documentation

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

override

Gathering the HepMCProduct information

Definition at line 69 of file DrellYanValidation.cc.

References _flavor, funct::abs(), funct::cos(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), fPDGTable, gamma_energy, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, mps_fire::i, leadeta, leadpt, RazorAnalyzer::lep1, RazorAnalyzer::lep2, nEvt, source_particleGun_cfi::ParticleID, HiggsValidation_cfi::pdg_id, perp(), edm::es::products(), seceta, secpt, jetUpdater_cfi::sort, HepMCValidationHelper::sortByPt(), protons_cff::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

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

override

Setting the DQM top directories

Definition at line 34 of file DrellYanValidation.cc.

References _name, DQMHelper::book1dHisto(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, gamma_energy, leadeta, leadpt, nEvt, seceta, secpt, DQMStore::IBooker::setCurrentFolder(), 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;
 }

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

override

Definition at line 30 of file DrellYanValidation.cc.

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

                                                                             {
   c.getData( fPDGTable );
 }

Member Data Documentation

int DrellYanValidation::_flavor

private

decay flavor

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze().

std::string DrellYanValidation::_name

private

decay flavor name

Definition at line 59 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(), and plotting.PlotGroup::remove().

MonitorElement * DrellYanValidation::cos_theta_gamma_lepton

private

Definition at line 54 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::dilep_mass

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_massPeak

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_pt

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_ptLog

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_rap

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

PDT table.

Definition at line 48 of file DrellYanValidation.h.

Referenced by analyze(), and dqmBeginRun().

MonitorElement* DrellYanValidation::gamma_energy

private

Definition at line 54 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag DrellYanValidation::hepmcCollection_

private

Definition at line 45 of file DrellYanValidation.h.

Referenced by DrellYanValidation().

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

private

Definition at line 61 of file DrellYanValidation.h.

Referenced by analyze(), and DrellYanValidation().

MonitorElement * DrellYanValidation::leadeta

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::leadpt

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::nEvt

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::seceta

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::secpt

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

WeightManager DrellYanValidation::wmanager_

private

Definition at line 44 of file DrellYanValidation.h.

Referenced by analyze().

MonitorElement * DrellYanValidation::Zdaughters

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::Zmass

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::ZmassPeak

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::Zpt

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::ZptLog

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::Zrap

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation