#include <TTbar_Kinematics.h>

Inheritance diagram for TTbar_Kinematics:

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

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

	TTbar_Kinematics (const edm::ParameterSet &)

	~TTbar_Kinematics () 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
edm::InputTag	genEventInfoProductTag_

edm::EDGetTokenT< GenEventInfoProduct >	genEventInfoProductTagToken_

edm::InputTag	genEvt_

MonitorElement *	hBottomE

MonitorElement *	hBottomEta

MonitorElement *	hBottomEtaPt

MonitorElement *	hBottomEtaPz

MonitorElement *	hBottomMass

MonitorElement *	hBottomMassDeltaY

MonitorElement *	hBottomMassEta

MonitorElement *	hBottomMassPz

MonitorElement *	hBottomMassY

MonitorElement *	hBottomPt

MonitorElement *	hBottomPtPz

MonitorElement *	hBottomPz

MonitorElement *	hBottomY

MonitorElement *	hBottomYPz

edm::InputTag	hepmcCollection_

edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_

MonitorElement *	hTopMass

MonitorElement *	hTopPt

MonitorElement *	hTopY

MonitorElement *	hTTbarMass

MonitorElement *	hTTbarPt

MonitorElement *	hTTbarY

MonitorElement *	hWminPz

MonitorElement *	hWplusPz

MonitorElement *	nEvt

TLorentzVector	tlv_Bottom

TLorentzVector	tlv_BottomBar

TLorentzVector	tlv_Top

TLorentzVector	tlv_TopBar

TLorentzVector	tlv_TTbar

TLorentzVector	tlv_Wmin

TLorentzVector	tlv_Wplus

double	weight

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 57 of file TTbar_Kinematics.h.

Constructor & Destructor Documentation

TTbar_Kinematics::TTbar_Kinematics ( const edm::ParameterSet & iConfig )

explicit

Definition at line 8 of file TTbar_Kinematics.cc.

References genEventInfoProductTag_, genEventInfoProductTagToken_, hepmcCollection_, and hepmcCollectionToken_.

                                                                  :
   hepmcCollection_(iConfig.getParameter<edm::InputTag>("hepmcCollection"))
   ,genEventInfoProductTag_(iConfig.getParameter<edm::InputTag>("genEventInfoProductTag"))
 {
 
   hepmcCollectionToken_=consumes<HepMCProduct>(hepmcCollection_);
   genEventInfoProductTagToken_=consumes<GenEventInfoProduct>(genEventInfoProductTag_);
 }

TTbar_Kinematics::~TTbar_Kinematics ( )

override

Definition at line 18 of file TTbar_Kinematics.cc.

18 {}

Member Function Documentation

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

override

Gathering the HepMCProduct information

Definition at line 26 of file TTbar_Kinematics.cc.

References PdtPdgMini::anti_b, PdtPdgMini::anti_t, PdtPdgMini::b, class-composition::children, MonitorElement::Fill(), genEventInfoProductTagToken_, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hBottomE, hBottomEta, hBottomEtaPt, hBottomEtaPz, hBottomMass, hBottomMassDeltaY, hBottomMassEta, hBottomMassPz, hBottomMassY, hBottomPt, hBottomPtPz, hBottomPz, hBottomY, hBottomYPz, hepmcCollectionToken_, hTopMass, hTopPt, hTopY, hTTbarMass, hTTbarPt, hTTbarY, hWminPz, hWplusPz, edm::HandleBase::isValid(), nEvt, PdtPdgMini::t, tlv_Bottom, tlv_BottomBar, tlv_Top, tlv_TopBar, tlv_TTbar, tlv_Wmin, tlv_Wplus, PdtPdgMini::W_minus, PdtPdgMini::W_plus, GenEventInfoProduct::weight(), and weight.

 {
 
   // --- the MC weights ---
   Handle<GenEventInfoProduct> evt_info;
   iEvent.getByToken(genEventInfoProductTagToken_, evt_info);
   if(!evt_info.isValid()) return;
   weight = evt_info->weight() ;
 
 
 
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(hepmcCollectionToken_, evt);
 
   //Get EVENT
   const HepMC::GenEvent *myGenEvent = evt->GetEvent();
 
   TLorentzVector tlv_Top, tlv_TopBar, tlv_Bottom, tlv_BottomBar ,tlv_Wplus ,tlv_Wmin , tlv_TTbar;
   bool top(false), antitop(false), antibottom(false), bottom(false), Wplus(false), Wmin(false);
   for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); iter != myGenEvent->particles_end(); iter++) {
     if((*iter)->pdg_id()==PdtPdgMini::t || (*iter)->pdg_id()==PdtPdgMini::anti_t){
       if( (*iter)->end_vertex()){
     HepMC::GenVertex::particle_iterator des;
     for(des = (*iter)->end_vertex()->particles_begin(HepMC::children);des!= (*iter)->end_vertex()->particles_end(HepMC::children);++des ){
       if((*des)->pdg_id()==PdtPdgMini::b){
         tlv_Bottom.SetPxPyPzE((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e());
         bottom=true;
       }
           if((*des)->pdg_id()==PdtPdgMini::anti_b){
         antibottom=true;
         tlv_BottomBar.SetPxPyPzE((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e());
          }
       if((*des)->pdg_id()==PdtPdgMini::W_plus){ 
         tlv_TopBar.SetPxPyPzE((*iter)->momentum().px(),(*iter)->momentum().py(),(*iter)->momentum().pz(),(*iter)->momentum().e()); antitop=true;
         tlv_Wplus.SetPxPyPzE((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e()); Wplus=true;
       }
       if((*des)->pdg_id()==PdtPdgMini::W_minus){ 
         tlv_Top.SetPxPyPzE((*iter)->momentum().px(),(*iter)->momentum().py(),(*iter)->momentum().pz(),(*iter)->momentum().e()); top=true;
         tlv_Wmin.SetPxPyPzE((*des)->momentum().px(),(*des)->momentum().py(),(*des)->momentum().pz(),(*des)->momentum().e()); Wmin=true;
       }
     }
       }
     }
   }
       
   tlv_TTbar = tlv_Top + tlv_TopBar ;
 
   //---topquarkquantities---
   nEvt->Fill(0.5,weight);
   if(top && antitop){
     hTopPt->Fill(tlv_Top.Pt(),weight);
     hTopPt->Fill(tlv_TopBar.Pt(),weight);
     
     hTopY->Fill(tlv_Top.Rapidity(),weight);
     hTopY->Fill(tlv_TopBar.Rapidity(),weight);
     
     hTopMass->Fill(tlv_Top.M(),weight);
     hTopMass->Fill(tlv_TopBar.M(),weight);
     
     //---ttbarpairquantities---
     hTTbarPt->Fill(tlv_TTbar.Pt(),weight);
     hTTbarY->Fill(tlv_TTbar.Rapidity(),weight);
     hTTbarMass->Fill(tlv_TTbar.M(),weight);
   }
   if(bottom && antibottom){
     hBottomPt->Fill(tlv_Bottom.Pt(),weight);
     hBottomPt->Fill(tlv_BottomBar.Pt(),weight);
     
     hBottomEta->Fill(tlv_Bottom.Eta(),weight);
     hBottomEta->Fill(tlv_BottomBar.Eta(),weight);
     
     //hBottomY->Fill(math::XYZTLorentzVector(bottom->momentum()).Rapidity(),weight);
     //hBottomY->Fill(math::XYZTLorentzVector(antibottom->momentum()).Rapidity(),weight);
     
     hBottomY->Fill(tlv_Bottom.Rapidity(),weight);
     hBottomY->Fill(tlv_BottomBar.Rapidity(),weight);
     
     hBottomPz->Fill(tlv_Bottom.Pz(),weight);
     hBottomPz->Fill(tlv_BottomBar.Pz(),weight);
     
     hBottomE->Fill(tlv_Bottom.E(),weight);
     hBottomE->Fill(tlv_BottomBar.E(),weight);
     
     hBottomMass->Fill(tlv_Bottom.M(),weight);
     hBottomMass->Fill(tlv_BottomBar.M(),weight);
     
     hBottomPtPz->Fill(tlv_Bottom.Pt(),tlv_Bottom.Pz(),weight);
     hBottomPtPz->Fill(tlv_BottomBar.Pt(),tlv_BottomBar.Pz(),weight);
     
     hBottomEtaPz->Fill(tlv_Bottom.Eta(),tlv_Bottom.Pz(),weight);
     hBottomEtaPz->Fill(tlv_BottomBar.Eta(),tlv_BottomBar.Pz(),weight);
     
     hBottomEtaPt->Fill(tlv_Bottom.Eta(),tlv_Bottom.Pt(),weight);
     hBottomEtaPt->Fill(tlv_BottomBar.Eta(),tlv_BottomBar.Pt(),weight);
     
     hBottomYPz->Fill(tlv_Bottom.Rapidity(),tlv_Bottom.Pz(),weight);
     hBottomYPz->Fill(tlv_BottomBar.Rapidity(),tlv_BottomBar.Pz(),weight);
     
     hBottomMassPz->Fill(tlv_Bottom.M(),tlv_Bottom.Pz(),weight);
     hBottomMassPz->Fill(tlv_BottomBar.M(),tlv_BottomBar.Pz(),weight);
     
     hBottomMassEta->Fill(tlv_Bottom.M(),tlv_Bottom.Eta(),weight);
     hBottomMassEta->Fill(tlv_BottomBar.M(),tlv_BottomBar.Eta(),weight);
     
     hBottomMassY->Fill(tlv_Bottom.M(),tlv_Bottom.Rapidity(),weight);
     hBottomMassY->Fill(tlv_BottomBar.M(),tlv_BottomBar.Rapidity(),weight);
     
     hBottomMassDeltaY->Fill(tlv_Bottom.M(),tlv_Bottom.Eta()-tlv_Bottom.Rapidity(),weight);
     hBottomMassDeltaY->Fill(tlv_BottomBar.M(),tlv_BottomBar.Eta()-tlv_BottomBar.Rapidity(),weight);
   }
   if(Wplus && Wmin){
     hWplusPz->Fill(tlv_Wplus.Pz(),weight);
     hWminPz->Fill(tlv_Wmin.Pz(),weight);
   }
 }

void TTbar_Kinematics::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	r,
		edm::EventSetup const &	e
	)

override

Definition at line 145 of file TTbar_Kinematics.cc.

References DQMHelper::book1dHisto(), DQMHelper::book2dHisto(), hBottomE, hBottomEta, hBottomEtaPt, hBottomEtaPz, hBottomMass, hBottomMassDeltaY, hBottomMassEta, hBottomMassPz, hBottomMassY, hBottomPt, hBottomPtPz, hBottomPz, hBottomY, hBottomYPz, hTopMass, hTopPt, hTopY, hTTbarMass, hTTbarPt, hTTbarY, hWminPz, hWplusPz, nEvt, and DQMStore::IBooker::setCurrentFolder().

                                                                                                 {
   DQMHelper dqm(&i); i.setCurrentFolder("Generator/TTbar");
 
   nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1.,"bins","Number of Events");
 
   hTopPt         = dqm.book1dHisto("TTbar_TopPt","t quark transverse momentum",1000,0.,1000.,"P_{t}^{t quark} (GeV)","Number of Events");
   hTopY          = dqm.book1dHisto("TTbar_TopY","t quark rapidity",200,-5.,5.,"Y_{t quark}","Number of Events");
   hTopMass       = dqm.book1dHisto("TTbar_TopMass","t quark mass",500,0.,500.,"M_{t quark} (GeV)","Number of Events");
  
   hTTbarPt       = dqm.book1dHisto("TTbar_TTbarPt","tt pair transverse momentum",1000,0.,1000.,"P_{t}^{tt pair} (GeV)","Number of Events");
   hTTbarY        = dqm.book1dHisto("TTbar_TTbarY","tt pair rapidity",200,-5.,5.,"Y_{tt pair}","Number of Events");
   hTTbarMass     = dqm.book1dHisto("TTbar_TTbarMass","tt pair mass",1000,0.,1000.,"M_{tt pair} (GeV)","Number of Events");
 
   hBottomPt      = dqm.book1dHisto("TTbar_BottomPt","b quark transverse momentum",1000,0.,1000.,"P_{t}^{b quark} (GeV)","Number of Events");
   hBottomEta     = dqm.book1dHisto("TTbar_BottomEta","b quark pseudo-rapidity",200,-5.,5.,"#eta_{b quark} (GeV)","Number of Events");
   hBottomY       = dqm.book1dHisto("TTbar_BottomY","b quark rapidity",200,-5.,5.,"M_{b quark} (GeV)","Number of Events");
   hBottomPz      = dqm.book1dHisto("TTbar_BottomPz","b quark longitudinal momentum",200,-100.,100.,"P_{z}^{b quark} (GeV)","Number of Events");
   hBottomE       = dqm.book1dHisto("TTbar_BottomE","b quark energy",1000,0.,1000.,"E_{b quark} (GeV)","Number of Events");
   hBottomMass    = dqm.book1dHisto("TTbar_BottomMass","b quark mass",50,0.,5.,"M_{b quark} (GeV)","Number of Events");
   
   hBottomPtPz    = dqm.book2dHisto("TTbar_BottomPtPz","b quark longitudinal vs transverse momentum",1000,0.,1000.,200,-100.,100.,"P_{z}^{b quark} (GeV)","P_{t}^{b quark} (GeV)");
   hBottomEtaPz   = dqm.book2dHisto("TTbar_BottomEtaPz","b quark longitudinal momentum vs pseudorapidity",200,-5.,5.,200,-100.,100.,"#eta_{b quark}","P_{z}^{b quark} (GeV)");
   hBottomEtaPt   = dqm.book2dHisto("TTbar_BottomEtaPt"," quark transveral   momentum vs pseudorapidity",200,-5.,5.,1000,0.,1000.,"#eta_{b quark}","P_{t}^{b quark} (GeV)");
   hBottomYPz     = dqm.book2dHisto("TTbar_BottomYPz","b quark longitudinal momentum vs rapidity",200,-5.,5.,200,-100.,100.,"Y_{b quark}","P_{z}^{b quark} (GeV)");
   hBottomMassPz  = dqm.book2dHisto("TTbar_BottomMassPz","b quark longitudinal momentum vs mass",50,0.,5.,200,-100.,100.,"M_{b quark} (GeV)","P_{z}^{b quark} (GeV)");
   hBottomMassEta = dqm.book2dHisto("TTbar_BottomMassEta","b quark pseudorapidity vs mass",50,0.,5.,200,-5.,5.,"M_{b quark} (GeV)","#eta_{b quark}");
   hBottomMassY   = dqm.book2dHisto("TTbar_BottomMassY","b quark rapidity vs mass",50,0.,5.,200,-5.,5.,"M_{b quark} (GeV)","Y_{b quark}");
   hBottomMassDeltaY = dqm.book2dHisto("TTbar_BottomMassDeltaY","b quark pseudorapidity - rapidity vs mass",50,0.,50.,2000,-5.,5.,"M_{b quark} (GeV)","Y_{b quark}");
 
   hWplusPz       = dqm.book1dHisto("TTbar_WplusPz","W+ boson longitudinal momentum",200,-100.,100.,"P_{z}^{W+} (GeV)","Number of Events");
   hWminPz        = dqm.book1dHisto("TTbar_WminPz","W- boson longitudinal momentum",200,-100.,100.,"P_{z}^{W-} (GeV)","Number of Events");
 
 }