#include <TopGenEventAnalyzer.h>

Inheritance diagram for TopGenEventAnalyzer:

Public Member Functions
	TopGenEventAnalyzer (const edm::ParameterSet &)

	~TopGenEventAnalyzer ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob ()

virtual void	endJob ()

Private Attributes
edm::EDGetTokenT< TtGenEvent >	inputGenEventToken_

TH1F *	nLep_

TH1F *	prodChan_

TH1F *	topBarEta_

TH1F *	topBarPhi_

TH1F *	topBarPt_

TH1F *	topEta_

TH1F *	topPhi_

TH1F *	topPt_

TH1F *	ttbarEta_

TH1F *	ttbarPhi_

TH1F *	ttbarPt_

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 16 of file TopGenEventAnalyzer.h.

Constructor & Destructor Documentation

TopGenEventAnalyzer::TopGenEventAnalyzer ( const edm::ParameterSet & cfg )

explicit

Definition at line 3 of file TopGenEventAnalyzer.cc.

References TFileService::make(), nLep_, prodChan_, topBarEta_, topBarPhi_, topBarPt_, topEta_, topPhi_, topPt_, ttbarEta_, ttbarPhi_, and ttbarPt_.

                                                                   :
   inputGenEventToken_(consumes<TtGenEvent>(cfg.getParameter<edm::InputTag>("genEvent")))
 {
   edm::Service<TFileService> fs;
   nLep_      = fs->make<TH1F>("nLep",      "N(Lepton)",     5,   0.,   5.);
   topPt_     = fs->make<TH1F>("topPt",     "pt (top)",    100,   0., 500.);
   topEta_    = fs->make<TH1F>("topEta",    "eta(top)",     40,  -5.,   5.);
   topPhi_    = fs->make<TH1F>("topPhi",    "phi(top)",     60, -3.5,  3.5);
   topBarPt_  = fs->make<TH1F>("topBarPt",  "pt (topBar)", 100,   0., 500.);
   topBarEta_ = fs->make<TH1F>("topBarEta", "eta(topBar)",  40,  -5.,   5.);
   topBarPhi_ = fs->make<TH1F>("topBarPhi", "phi(topBar)",  60, -3.5,  3.5);
   ttbarPt_   = fs->make<TH1F>("ttbarPt",   "pt (ttbar)",  100,   0., 500.);
   ttbarEta_  = fs->make<TH1F>("ttbarEta",  "eta(ttbar)",   40,  -5.,   5.);
   ttbarPhi_  = fs->make<TH1F>("ttbarPhi",  "phi(ttbar)",   60, -3.5,  3.5);
   prodChan_  = fs->make<TH1F>("prodChan",  "production mode", 3, 0, 3);
   prodChan_->GetXaxis()->SetBinLabel(1, "gg"   );
   prodChan_->GetXaxis()->SetBinLabel(2, "qqbar");
   prodChan_->GetXaxis()->SetBinLabel(3, "other");
 }

TopGenEventAnalyzer::~TopGenEventAnalyzer ( )

Definition at line 23 of file TopGenEventAnalyzer.cc.

24 {

25 }

Member Function Documentation

void TopGenEventAnalyzer::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	setup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 28 of file TopGenEventAnalyzer.cc.

References MCTruth::genEvent, edm::Event::getByToken(), inputGenEventToken_, nLep_, prodChan_, topBarEta_, topBarPhi_, topBarPt_, topEta_, topPhi_, topPt_, ttbarEta_, ttbarPhi_, and ttbarPt_.

 {
   edm::Handle<TtGenEvent> genEvent;
   evt.getByToken(inputGenEventToken_, genEvent);
 
   if(!genEvent->isTtBar())
     return;
 
   if(genEvent->fromGluonFusion())
     prodChan_->Fill("gg", 1);
   else if(genEvent->fromQuarkAnnihilation())
     prodChan_->Fill("qqbar", 1);
   else
     prodChan_->Fill("other", 1);
 
   // fill BR's
   nLep_  ->Fill(genEvent->numberOfLeptons());
 
   //fill top kinematic
   topPt_    ->Fill(genEvent->top   ()->pt ());
   topEta_   ->Fill(genEvent->top   ()->eta());
   topPhi_   ->Fill(genEvent->top   ()->phi());
   topBarPt_ ->Fill(genEvent->topBar()->pt ());
   topBarEta_->Fill(genEvent->topBar()->eta());
   topBarPhi_->Fill(genEvent->topBar()->phi());
 
   //fill ttbar kinematics
   ttbarPt_ ->Fill(genEvent->topPair()->pt() );
   ttbarEta_->Fill(genEvent->topPair()->eta());
   ttbarPhi_->Fill(genEvent->topPair()->phi());
 }