#include <DrellYanValidation.h>

Inheritance diagram for DrellYanValidation:

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

virtual void	beginJob ()

virtual void	beginRun (const edm::Run &, const edm::EventSetup &)

	DrellYanValidation (const edm::ParameterSet &)

virtual void	endJob ()

virtual void	endRun (const edm::Run &, const edm::EventSetup &)

virtual	~DrellYanValidation ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Attributes
int	_flavor
	decay flavor More...

std::string	_name
	decay flavor name More...

MonitorElement *	cos_theta_gamma_lepton

DQMStore *	dbe
	ME's "container". More...

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_

MonitorElement *	leadeta

MonitorElement *	leadpt

MonitorElement *	nEvt

MonitorElement *	seceta

MonitorElement *	secpt

MonitorElement *	Zdaughters

MonitorElement *	Zmass

MonitorElement *	ZmassPeak

MonitorElement *	Zpt

MonitorElement *	ZptLog

MonitorElement *	Zrap

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 34 of file DrellYanValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file DrellYanValidation.cc.

References dbe, and cmsCodeRules.cppFunctionSkipper::operator.

                                                                   :  
   hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")),
   _flavor(iPSet.getParameter<int>("decaysTo")),
   _name(iPSet.getParameter<std::string>("name")) 
 {    
   dbe = 0;
   dbe = edm::Service<DQMStore>().operator->();
 }

DrellYanValidation::~DrellYanValidation ( )

virtual

Definition at line 29 of file DrellYanValidation.cc.

29 {}

Member Function Documentation

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

virtual

Gathering the HepMCProduct information

Implements edm::EDAnalyzer.

Definition at line 76 of file DrellYanValidation.cc.

References _flavor, 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::getByLabel(), hepmcCollection_, i, leadeta, leadpt, nEvt, RecoTau_DiTaus_pt_20-420_cfg::ParticleID, perp(), edm::es::products(), seceta, secpt, python.multivaluedict::sort(), HepMCValidationHelper::sortByPt(), matplotRender::t, x, detailsBasic3DVector::y, detailsBasic3DVector::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.getByLabel(hepmcCollection_, evt);
 
   //Get EVENT
   const HepMC::GenEvent *myGenEvent = evt->GetEvent();
 
   nEvt->Fill(0.5);
 
   std::vector<const HepMC::GenParticle*> allproducts; 
 
   //requires status 1 for leptons and neutrinos (except tau)
   int requiredstatus = (abs(_flavor) == 11 || abs(_flavor) == 12 || abs(_flavor) ==13 || abs(_flavor) ==14 || abs(_flavor) ==16) ? 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) {
       if(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()); 
   Zdaughters->Fill(products[1]->pdg_id()); 
 
   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());
     gammasMomenta.push_back(phomom);
   }  
   //Fill Z histograms
   Zmass->Fill(dilepton_andphoton_mom.M());
   ZmassPeak->Fill(dilepton_andphoton_mom.M());
   Zpt->Fill(dilepton_andphoton_mom.Pt());
   ZptLog->Fill(log10(dilepton_andphoton_mom.Pt())); 
   Zrap->Fill(dilepton_andphoton_mom.Rapidity());
 
   //Fill dilepton histograms
   dilep_mass->Fill(dilepton_mom.M());
   dilep_massPeak->Fill(dilepton_mom.M());
   dilep_pt->Fill(dilepton_mom.Pt());
   dilep_ptLog->Fill(log10(dilepton_mom.Pt()));
   dilep_rap->Fill(dilepton_mom.Rapidity()); 
 
   //Fill lepton histograms 
   leadpt->Fill(lep1.Pt());
   secpt->Fill(lep2.Pt());
   leadeta->Fill(lep1.Eta());
   seceta->Fill(lep2.Eta());
 
   //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.size() != 0 && dilepton_andphoton_mom.M() > 50.) {
     gamma_energy->Fill(gammasMomenta.front().E());
     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));
   } 
 
 }//analyze

void DrellYanValidation::beginJob ( void )

virtual

Setting the DQM top directories

Reimplemented from edm::EDAnalyzer.

Definition at line 31 of file DrellYanValidation.cc.

References _name, DQMStore::book1D(), cos_theta_gamma_lepton, dbe, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, gamma_energy, leadeta, leadpt, nEvt, seceta, secpt, DQMStore::setCurrentFolder(), Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.

 {
   if(dbe){
     std::string folderName = "Generator/DrellYan";
     folderName+=_name;
     dbe->setCurrentFolder(folderName.c_str());
     
     // Number of analyzed events
     nEvt = dbe->book1D("nEvt", "n analyzed Events", 1, 0., 1.);
     
     //Kinematics
     Zmass = dbe->book1D("Zmass","inv. Mass Z", 70 ,0,140);
     ZmassPeak = dbe->book1D("ZmassPeak","inv. Mass Z", 80 ,80 ,100);
     Zpt = dbe->book1D("Zpt","Z pt",100,0,200);
     ZptLog = dbe->book1D("ZptLog","log(Z pt)",100,0.,5.);
     Zrap = dbe->book1D("Zrap", "Z y", 100, -5, 5);
     Zdaughters = dbe->book1D("Zdaughters", "Z daughters", 60, -30, 30);
 
     dilep_mass = dbe->book1D("dilep_mass","inv. Mass dilepton", 70 ,0,140);
     dilep_massPeak = dbe->book1D("dilep_massPeak","inv. Mass dilepton", 80 ,80 ,100);
     dilep_pt = dbe->book1D("dilep_pt","dilepton pt",100,0,200);
     dilep_ptLog = dbe->book1D("dilep_ptLog","log(dilepton pt)",100,0.,5.);
     dilep_rap = dbe->book1D("dilep_rap", "dilepton y", 100, -5, 5);
 
     gamma_energy = dbe->book1D("gamma_energy", "photon energy in Z rest frame", 200, 0., 100.);
     cos_theta_gamma_lepton = dbe->book1D("cos_theta_gamma_lepton",      "cos_theta_gamma_lepton in Z rest frame",      200, -1, 1);
 
     leadpt = dbe->book1D("leadpt","leading lepton pt", 200, 0., 200.);    
     secpt  = dbe->book1D("secpt","second lepton pt", 200, 0., 200.);    
     leadeta = dbe->book1D("leadeta","leading lepton eta", 100, -5., 5.);
     seceta  = dbe->book1D("seceta","second lepton eta", 100, -5., 5.);
 
   }
   return;
 }

void DrellYanValidation::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

virtual

Get PDT Table

Reimplemented from edm::EDAnalyzer.

Definition at line 69 of file DrellYanValidation.cc.

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

 {
   iSetup.getData( fPDGTable );
   return;
 }

void DrellYanValidation::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 68 of file DrellYanValidation.cc.

68 {return;}

void DrellYanValidation::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 75 of file DrellYanValidation.cc.

75 {return;}

Member Data Documentation

int DrellYanValidation::_flavor

private

decay flavor

Definition at line 62 of file DrellYanValidation.h.

Referenced by analyze().

std::string DrellYanValidation::_name

private

decay flavor name

Definition at line 64 of file DrellYanValidation.h.

Referenced by beginJob().

MonitorElement * DrellYanValidation::cos_theta_gamma_lepton

private

Definition at line 59 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

DQMStore* DrellYanValidation::dbe

private

ME's "container".

Definition at line 53 of file DrellYanValidation.h.

Referenced by beginJob(), and DrellYanValidation().

MonitorElement* DrellYanValidation::dilep_mass

private

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::dilep_massPeak

private

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::dilep_pt

private

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::dilep_ptLog

private

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::dilep_rap

private

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

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

private

PDT table.

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and beginRun().

MonitorElement* DrellYanValidation::gamma_energy

private

Definition at line 59 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

edm::InputTag DrellYanValidation::hepmcCollection_

private

Definition at line 47 of file DrellYanValidation.h.

Referenced by analyze().

MonitorElement * DrellYanValidation::leadeta

private

Definition at line 58 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement* DrellYanValidation::leadpt

private

Definition at line 58 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement* DrellYanValidation::nEvt

private

Definition at line 55 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::seceta

private

Definition at line 58 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::secpt

private

Definition at line 58 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::Zdaughters

private

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement* DrellYanValidation::Zmass

private

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::ZmassPeak

private

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::Zpt

private

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::ZptLog

private

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

MonitorElement * DrellYanValidation::Zrap

private

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze(), and beginJob().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation