#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 DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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
	)

overridevirtual

Implements DQMEDAnalyzer.

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");
 
 }