#include <TTbar_GenJetAnalyzer.h>

Inheritance diagram for TTbar_GenJetAnalyzer:

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_GenJetAnalyzer (const edm::ParameterSet &)

	~TTbar_GenJetAnalyzer () 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_

std::map< std::string, MonitorElement * >	hists_

edm::InputTag	jets_

edm::EDGetTokenT< std::vector< reco::GenJet > >	jetsToken_

double	weight

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 55 of file TTbar_GenJetAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 4 of file TTbar_GenJetAnalyzer.cc.

References genEventInfoProductTag_, genEventInfoProductTagToken_, jets_, and jetsToken_.

                                                                         :
   jets_(iConfig.getParameter<edm::InputTag>("jets")),
   genEventInfoProductTag_(iConfig.getParameter<edm::InputTag>("genEventInfoProductTag"))
 {
 
   genEventInfoProductTagToken_=consumes<GenEventInfoProduct>(genEventInfoProductTag_);
   jetsToken_=consumes<std::vector<reco::GenJet> >(jets_);
 
 }

TTbar_GenJetAnalyzer::~TTbar_GenJetAnalyzer ( )

override

Definition at line 15 of file TTbar_GenJetAnalyzer.cc.

15 {}

Member Function Documentation

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

override

Definition at line 18 of file TTbar_GenJetAnalyzer.cc.

References genEventInfoProductTagToken_, edm::Event::getByToken(), hists_, edm::HandleBase::isValid(), fwrapper::jets, jetsToken_, GenEventInfoProduct::weight(), and weight.

 {
   using namespace edm;
  
   // --- the MC weights ---
   Handle<GenEventInfoProduct> evt_info;
   iEvent.getByToken(genEventInfoProductTagToken_, evt_info);
   if(!evt_info.isValid()) return;
   weight = evt_info->weight() ;
 
   // Gather information in the GenJet collection
   edm::Handle<std::vector<reco::GenJet> > jets;
   iEvent.getByToken(jetsToken_, jets );
 
   if(!jets.isValid()) return;
   // loop Jet collection and fill histograms
   int njets = 0;
   for(std::vector<reco::GenJet>::const_iterator jet_it=jets->begin(); jet_it!=jets->end(); ++jet_it){  
 
     ++njets; 
 
     hists_["jetPtAll" ]->Fill( jet_it->pt() , weight );
     hists_["jetEtaAll"]->Fill( jet_it->eta(), weight );
 
     if (njets == 1) { hists_["jetPt1" ]->Fill( jet_it->pt() , weight );
                       hists_["jetEta1"]->Fill( jet_it->eta(), weight );
                     }
     if (njets == 2) { hists_["jetPt2" ]->Fill( jet_it->pt() , weight );
                       hists_["jetEta2"]->Fill( jet_it->eta(), weight );
                     }
     if (njets == 3) { hists_["jetPt3" ]->Fill( jet_it->pt() , weight );
                       hists_["jetEta3"]->Fill( jet_it->eta(), weight );
                     }
     if (njets == 4) { hists_["jetPt4" ]->Fill( jet_it->pt() , weight );
                       hists_["jetEta4"]->Fill( jet_it->eta(), weight );
                     }
   }
 
 }

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

override

Definition at line 59 of file TTbar_GenJetAnalyzer.cc.

References DQMHelper::book1dHisto(), hists_, and DQMStore::IBooker::setCurrentFolder().

                                                                                                   {
   DQMHelper dqm(&i); i.setCurrentFolder("Generator/TTbar");
   hists_["jetPtAll" ] = dqm.book1dHisto("TTbar_jetPtAll" , "pt" , 1000,  0., 1000.,"P_{t}^{All-Jets} (GeV)","Number of Events"); 
   hists_["jetPt1"   ] = dqm.book1dHisto("TTbar_jetPt1"   , "pt" , 1000,  0., 1000.,"P_{t}^{1st-Jet} (GeV)","Number of Events"); 
   hists_["jetPt2"   ] = dqm.book1dHisto("TTbar_jetPt2"   , "pt" , 1000,  0., 1000.,"P_{t}^{2nd-Jet} (GeV)","Number of Events"); 
   hists_["jetPt3"   ] = dqm.book1dHisto("TTbar_jetPt3"   , "pt" , 1000,  0., 1000.,"P_{t}^{3rd-Jet} (GeV)","Number of Events"); 
   hists_["jetPt4"   ] = dqm.book1dHisto("TTbar_jetPt4"   , "pt" , 1000,  0., 1000.,"P_{t}^{4th-Jet} (GeV)","Number of Events"); 
                                                                                     
   hists_["jetEtaAll"] = dqm.book1dHisto("TTbar_jetEtaAll", "eta",  100, -5.,    5.,"#eta^{All-Jets}","Number of Events"); 
   hists_["jetEta1"  ] = dqm.book1dHisto("TTbar_jetEta1"  , "eta",  100, -5.,    5.,"#eta^{1st-Jet}","Number of Events"); 
   hists_["jetEta2"  ] = dqm.book1dHisto("TTbar_jetEta2"  , "eta",  100, -5.,    5.,"#eta^{2nd-Jet}","Number of Events"); 
   hists_["jetEta3"  ] = dqm.book1dHisto("TTbar_jetEta3"  , "eta",  100, -5.,    5.,"#eta^{3rd-Jet}","Number of Events"); 
   hists_["jetEta4"  ] = dqm.book1dHisto("TTbar_jetEta4"  , "eta",  100, -5.,    5.,"#eta^{4th-Jet}","Number of Events"); 
 }