#include <RazorVarAnalyzer.h>

Inheritance diagram for RazorVarAnalyzer:

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

virtual void	bookMEs ()

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

	RazorVarAnalyzer (const edm::ParameterSet &)

virtual	~RazorVarAnalyzer ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from ScoutingAnalyzerBase
void	beginJob ()

MonitorElement *	bookH1 (const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")

MonitorElement *	bookH1andDivide (const std::string &name, MonitorElement num, MonitorElement denom, const std::string &titleX, const std::string &titleY, const std::string &title="")

MonitorElement *	bookH1BinArray (const std::string &name, const std::string &title, int nchX, float xbinsize, const std::string &titleX="", const std::string &titleY="Events", Option_t option="E1 P")

MonitorElement *	bookH1withSumw2 (const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")

MonitorElement *	bookH1withSumw2BinArray (const std::string &name, const std::string &title, int nchX, float xbinsize, const std::string &titleX="", const std::string &titleY="Events", Option_t option="E1 P")

MonitorElement *	bookH2 (const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")

MonitorElement *	bookH2andDivide (const std::string &name, MonitorElement num, MonitorElement denom, const std::string &titleX, const std::string &titleY, const std::string &title="")

MonitorElement *	bookH2withSumw2 (const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")

MonitorElement *	bookP1 (const std::string &name, const std::string &title, int nchX, double lowX, double highX, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="E1 P")

virtual void	endJob ()

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

std::string	newName (const std::string &name)

MonitorElement *	profileX (MonitorElement *me2d, const std::string &title="", const std::string &titleX="", const std::string &titleY="", Double_t minimum=-1111, Double_t maximum=-1111)

MonitorElement *	profileY (MonitorElement *me2d, const std::string &title="", const std::string &titleX="", const std::string &titleY="", Double_t minimum=-1111, Double_t maximum=-1111)

	ScoutingAnalyzerBase (const edm::ParameterSet &conf)

virtual	~ScoutingAnalyzerBase ()

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

CurrentProcessingContext const *	currentContext () const

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

Constructor & Destructor Documentation

RazorVarAnalyzer::RazorVarAnalyzer ( const edm::ParameterSet & conf )

explicit

Definition at line 16 of file RazorVarAnalyzer.cc.

                                                                 :
   ScoutingAnalyzerBase(conf),
   m_eleCollectionTag(conf.getUntrackedParameter<edm::InputTag>("eleCollectionName",edm::InputTag("hltPixelMatchElectronsActivity"))),
   m_jetCollectionTag(conf.getUntrackedParameter<edm::InputTag>("jetCollectionName",edm::InputTag("hltCaloJetIDPassed"))),
   m_muCollectionTag(conf.getUntrackedParameter<edm::InputTag>("muCollectionName",edm::InputTag("hltL3MuonCandidates"))),
   m_razorVarCollectionTag(conf.getUntrackedParameter<edm::InputTag>("razorVarCollectionName")){
 }

RazorVarAnalyzer::~RazorVarAnalyzer ( )

virtual

Definition at line 26 of file RazorVarAnalyzer.cc.

26 {}

Member Function Documentation

void RazorVarAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from ScoutingAnalyzerBase.

Definition at line 30 of file RazorVarAnalyzer.cc.

References MonitorElement::Fill(), edm::Event::getByLabel(), edm::HandleBase::isValid(), m_eleCollectionTag, m_jetCollectionTag, m_muCollectionTag, m_razorVarCollectionTag, m_rsqMREle, m_rsqMREleEle, m_rsqMREleMJ, m_rsqMREleMu, m_rsqMRFullyInc, m_rsqMRHad, m_rsqMRHadMJ, m_rsqMRInc10J, m_rsqMRInc12J, m_rsqMRInc14J, m_rsqMRInc4J, m_rsqMRInc6J, m_rsqMRInc8J, m_rsqMRMu, m_rsqMRMuMJ, m_rsqMRMuMu, and dttmaxenums::R.

                                                                                 {
   
   //count the number of jets with a minimal selection
   edm::Handle<reco::CaloJetCollection> calojets_handle;
   iEvent.getByLabel(m_jetCollectionTag,calojets_handle);
   
   unsigned int njets = 0;
   for(reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it){
     if(it->pt() >= 30. && fabs(it->eta()) <= 3.0){
       njets++;
     }
   }
 
   //count the number of muons
   edm::Handle<std::vector<reco::RecoChargedCandidate> > muon_handle;
   iEvent.getByLabel(m_muCollectionTag,muon_handle);
 
   unsigned int nmu_loose = 0;
   unsigned int nmu_tight = 0;
   if(muon_handle.isValid()){
     for(std::vector<reco::RecoChargedCandidate>::const_iterator it = muon_handle->begin(); it != muon_handle->end(); ++it){
       if(it->pt() >= 15 && fabs(it->eta()) <= 2.1) nmu_tight++;
       if(it->pt() >= 10 && fabs(it->eta()) <= 2.4) nmu_loose++;
     }
   }
   
   //count the number of electrons
   edm::Handle<reco::ElectronCollection> ele_handle;
   iEvent.getByLabel(m_eleCollectionTag,ele_handle);
   
   unsigned int nele_loose = 0;
   unsigned int nele_tight = 0;
   if(ele_handle.isValid()){
     for(reco::ElectronCollection::const_iterator it = ele_handle->begin(); it != ele_handle->end(); ++it){
       if(it->pt() >= 20 && fabs(it->eta()) <= 2.5) nele_tight++;
       if(it->pt() >= 10 && fabs(it->eta()) <= 2.5) nele_loose++;
     }
   }
   
   //now get the box number: {'MuEle':0,'MuMu':1,'EleEle':2,'Mu':3,'Ele':4,'Had':5}
   unsigned int box_num = 5;
   if(nmu_tight > 0 && nele_tight > 0){
     box_num = 0;
   }else if(nmu_tight > 0 && nmu_loose > 1){
     box_num = 1;
   }else if(nele_tight > 0 && nele_loose > 1){
     box_num = 2;
   }else if(nmu_tight > 0){
     box_num = 3;
   }else if(nele_tight > 0){
     box_num = 4;
   }
 
   edm::Handle<std::vector<double> > razorvar_handle;
   iEvent.getByLabel(m_razorVarCollectionTag,razorvar_handle);
 
   if(razorvar_handle->size() > 1){
     const double MR = razorvar_handle->at(0);
     const double R = razorvar_handle->at(1);
     m_rsqMRFullyInc->Fill(MR,R*R);
     if(njets >= 4) m_rsqMRInc4J->Fill(MR,R*R);
     if(njets >= 6) m_rsqMRInc6J->Fill(MR,R*R);
     if(njets >= 8) m_rsqMRInc8J->Fill(MR,R*R);
     if(njets >= 10) m_rsqMRInc10J->Fill(MR,R*R);
     if(njets >= 12) m_rsqMRInc12J->Fill(MR,R*R);
     if(njets >= 14) m_rsqMRInc14J->Fill(MR,R*R);
     
     //now fill the boxes
     if(box_num == 0) m_rsqMREleMu->Fill(MR,R*R);
     if(box_num == 1) m_rsqMRMuMu->Fill(MR,R*R);
     if(box_num == 2) m_rsqMREleEle->Fill(MR,R*R);
     if(box_num == 3) m_rsqMRMu->Fill(MR,R*R);
     if(box_num == 4) m_rsqMREle->Fill(MR,R*R);
     if(box_num == 5) m_rsqMRHad->Fill(MR,R*R);
 
     //finally the multijet boxes - think ttbar
     //muon boxes: muons are not in jets
     if( box_num == 3 && njets >= 4) m_rsqMRMuMJ->Fill(MR,R*R);
     //ele boxes: electrons are in jets
     else if( box_num == 4 && njets >= 5) m_rsqMREleMJ->Fill(MR,R*R);
     //fill the Had box
     else if( box_num == 5 && njets >= 6) m_rsqMRHadMJ->Fill(MR,R*R);
   }
   
 }

void RazorVarAnalyzer::bookMEs ( )

virtual

Reimplemented from ScoutingAnalyzerBase.

Definition at line 121 of file RazorVarAnalyzer.cc.

References ScoutingAnalyzerBase::bookH2withSumw2(), m_rsqMREle, m_rsqMREleEle, m_rsqMREleMJ, m_rsqMREleMu, m_rsqMRFullyInc, m_rsqMRHad, m_rsqMRHadMJ, m_rsqMRInc10J, m_rsqMRInc12J, m_rsqMRInc14J, m_rsqMRInc4J, m_rsqMRInc6J, m_rsqMRInc8J, m_rsqMRMu, m_rsqMRMuMJ, and m_rsqMRMuMu.

                               {
   
   //the full inclusive histograms
   m_rsqMRFullyInc = bookH2withSumw2("rsqMRFullyInc",
                     "M_{R} vs R^{2} (All Events)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
   m_rsqMRInc4J = bookH2withSumw2("rsqMRInc4J",
                     "M_{R} vs R^{2} (>= 4j)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
   m_rsqMRInc6J = bookH2withSumw2("rsqMRInc6J",
                     "M_{R} vs R^{2} (>= 6j)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
   m_rsqMRInc8J = bookH2withSumw2("rsqMRInc8J",
                     "M_{R} vs R^{2} (>= 8j)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
   m_rsqMRInc10J = bookH2withSumw2("rsqMRInc10J",
                     "M_{R} vs R^{2} (>= 10j)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
   m_rsqMRInc12J = bookH2withSumw2("rsqMRInc12J",
                     "M_{R} vs R^{2} (>= 12j)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
   m_rsqMRInc14J = bookH2withSumw2("rsqMRInc14J",
                     "M_{R} vs R^{2} (>= 14j)",
                     400,0.,4000.,
                     50,0.,1.,
                     "M_{R} [GeV]",
                     "R^{2}");
 
   //the by box histograms
   m_rsqMREleMu = bookH2withSumw2("rsqMREleMu",
                  "M_{R} vs R^{2} (EleMu box)",
                  400,0.,4000.,
                  50,0.,1.,
                  "M_{R} [GeV]",
                  "R^{2}");
   m_rsqMRMuMu = bookH2withSumw2("rsqMRMuMu",
                 "M_{R} vs R^{2} (MuMu box)",
                 400,0.,4000.,
                 50,0.,1.,
                 "M_{R} [GeV]",
                 "R^{2}");
   m_rsqMREleEle = bookH2withSumw2("rsqMREleEle",
                   "M_{R} vs R^{2} (EleEle box)",
                   400,0.,4000.,
                   50,0.,1.,
                   "M_{R} [GeV]",
                   "R^{2}");
   m_rsqMRMu = bookH2withSumw2("rsqMRMu",
                   "M_{R} vs R^{2} (Mu box)",
                   400,0.,4000.,
                   50,0.,1.,
                   "M_{R} [GeV]",
                   "R^{2}");
   m_rsqMREle = bookH2withSumw2("rsqMREle",
                    "M_{R} vs R^{2} (Ele box)",
                    400,0.,4000.,
                    50,0.,1.,
                    "M_{R} [GeV]",
                    "R^{2}");
   m_rsqMRHad = bookH2withSumw2("rsqMRHad",
                   "M_{R} vs R^{2} (Had box)",
                   400,0.,4000.,
                   50,0.,1.,
                   "M_{R} [GeV]",
                   "R^{2}");
 
   //the by box histograms
   m_rsqMRMuMJ = bookH2withSumw2("rsqMRMuMJ",
                   "M_{R} vs R^{2} (Mu box >= 4j)",
                   400,0.,4000.,
                   50,0.,1.,
                   "M_{R} [GeV]",
                   "R^{2}");
   m_rsqMREleMJ = bookH2withSumw2("rsqMREleMJ",
                    "M_{R} vs R^{2} (Ele box >= 5j)",
                    400,0.,4000.,
                    50,0.,1.,
                    "M_{R} [GeV]",
                    "R^{2}");
   m_rsqMRHadMJ = bookH2withSumw2("rsqMRHadMJ",
                  "M_{R} vs R^{2} (Had box >= 6j)",
                  400,0.,4000.,
                  50,0.,1.,
                  "M_{R} [GeV]",
                  "R^{2}");
 
 
 }

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

virtual

Reimplemented from ScoutingAnalyzerBase.

Definition at line 116 of file RazorVarAnalyzer.cc.

116 {

117 }

Member Data Documentation

edm::InputTag RazorVarAnalyzer::m_eleCollectionTag

private

Definition at line 23 of file RazorVarAnalyzer.h.

Referenced by analyze().

edm::InputTag RazorVarAnalyzer::m_jetCollectionTag

private

Definition at line 24 of file RazorVarAnalyzer.h.

Referenced by analyze().

edm::InputTag RazorVarAnalyzer::m_muCollectionTag

private

Definition at line 25 of file RazorVarAnalyzer.h.

Referenced by analyze().

edm::InputTag RazorVarAnalyzer::m_razorVarCollectionTag

private

Definition at line 26 of file RazorVarAnalyzer.h.

Referenced by analyze().

MonitorElement* RazorVarAnalyzer::m_rsqMREle

private

Definition at line 42 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMREleEle

private

Definition at line 40 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMREleMJ

private

Definition at line 47 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMREleMu

private

Definition at line 38 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRFullyInc

private

Definition at line 29 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRHad

private

Definition at line 43 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRHadMJ

private

Definition at line 48 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRInc10J

private

Definition at line 33 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRInc12J

private

Definition at line 34 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRInc14J

private

Definition at line 35 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRInc4J

private

Definition at line 30 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRInc6J

private

Definition at line 31 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRInc8J

private

Definition at line 32 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRMu

private

Definition at line 41 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRMuMJ

private

Definition at line 46 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

MonitorElement* RazorVarAnalyzer::m_rsqMRMuMu

private

Definition at line 39 of file RazorVarAnalyzer.h.

Referenced by analyze(), and bookMEs().

Private Attributes
edm::InputTag	m_eleCollectionTag

edm::InputTag	m_jetCollectionTag

edm::InputTag	m_muCollectionTag

edm::InputTag	m_razorVarCollectionTag

MonitorElement *	m_rsqMREle

MonitorElement *	m_rsqMREleEle

MonitorElement *	m_rsqMREleMJ

MonitorElement *	m_rsqMREleMu

MonitorElement *	m_rsqMRFullyInc

MonitorElement *	m_rsqMRHad

MonitorElement *	m_rsqMRHadMJ

MonitorElement *	m_rsqMRInc10J

MonitorElement *	m_rsqMRInc12J

MonitorElement *	m_rsqMRInc14J

MonitorElement *	m_rsqMRInc4J

MonitorElement *	m_rsqMRInc6J

MonitorElement *	m_rsqMRInc8J

MonitorElement *	m_rsqMRMu

MonitorElement *	m_rsqMRMuMJ

MonitorElement *	m_rsqMRMuMu

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation