#include <RazorVarAnalyzer.h>

Inheritance diagram for RazorVarAnalyzer:

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

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

	RazorVarAnalyzer (const edm::ParameterSet &)

	~RazorVarAnalyzer () 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	m_eleCollectionTag

edm::EDGetTokenT< reco::ElectronCollection >	m_eleCollectionTagToken_

edm::InputTag	m_jetCollectionTag

edm::EDGetTokenT< reco::CaloJetCollection >	m_jetCollectionTagToken_

edm::InputTag	m_muCollectionTag

edm::EDGetTokenT< std::vector< reco::RecoChargedCandidate > >	m_muCollectionTagToken_

edm::InputTag	m_razorVarCollectionTag

edm::EDGetTokenT< std::vector< double > >	m_razorVarCollectionTagToken_

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

Additional Inherited Members
Protected Member Functions inherited from ScoutingAnalyzerBase
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

MonitorElement *	bookH1 (DQMStore::IBooker &, 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 (DQMStore::IBooker &, const std::string &name, MonitorElement num, MonitorElement denom, const std::string &titleX, const std::string &titleY, const std::string &title="")

MonitorElement *	bookH1BinArray (DQMStore::IBooker &, 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 (DQMStore::IBooker &, 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 (DQMStore::IBooker &, 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 (DQMStore::IBooker &, 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 (DQMStore::IBooker &, const std::string &name, MonitorElement num, MonitorElement denom, const std::string &titleX, const std::string &titleY, const std::string &title="")

MonitorElement *	bookH2withSumw2 (DQMStore::IBooker &, 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 (DQMStore::IBooker &, 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")

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

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

void	prepareBooking (DQMStore::IBooker &)

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

MonitorElement *	profileY (DQMStore::IBooker &, 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)

	~ScoutingAnalyzerBase () override

Detailed Description

Definition at line 10 of file RazorVarAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 15 of file RazorVarAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), m_eleCollectionTagToken_, m_jetCollectionTagToken_, m_muCollectionTagToken_, and m_razorVarCollectionTagToken_.

                                                                 :
   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")) {
     //set Token(-s)
     m_jetCollectionTagToken_ = consumes<reco::CaloJetCollection>(conf.getUntrackedParameter<edm::InputTag>("jetCollectionName",edm::InputTag("hltCaloJetIDPassed")));
     m_muCollectionTagToken_ = consumes<std::vector<reco::RecoChargedCandidate> >(conf.getUntrackedParameter<edm::InputTag>("muCollectionName",edm::InputTag("hltL3MuonCandidates")));
     m_eleCollectionTagToken_ = consumes<reco::ElectronCollection>(conf.getUntrackedParameter<edm::InputTag>("eleCollectionName",edm::InputTag("hltPixelMatchElectronsActivity")));
     m_razorVarCollectionTagToken_ = consumes<std::vector<double> >(conf.getUntrackedParameter<edm::InputTag>("razorVarCollectionName"));
 }

RazorVarAnalyzer::~RazorVarAnalyzer ( )

override

Definition at line 28 of file RazorVarAnalyzer.cc.

28 {}

Member Function Documentation

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

override

Definition at line 169 of file RazorVarAnalyzer.cc.

References MonitorElement::Fill(), edm::Event::getByToken(), edm::HandleBase::isValid(), m_eleCollectionTagToken_, m_jetCollectionTagToken_, m_muCollectionTagToken_, m_razorVarCollectionTagToken_, 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.getByToken(m_jetCollectionTagToken_, 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.getByToken(m_muCollectionTagToken_, 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.getByToken(m_eleCollectionTagToken_, 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.getByToken(m_razorVarCollectionTagToken_, 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::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

override

Definition at line 31 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, m_rsqMRMuMu, and ScoutingAnalyzerBase::prepareBooking().

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