|
|
|
#include <RazorVarAnalyzer.h>
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 () |
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 |
Definition at line 7 of file RazorVarAnalyzer.h.
| 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.
{}
| 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] |
Definition at line 23 of file RazorVarAnalyzer.h.
Referenced by analyze().
Definition at line 24 of file RazorVarAnalyzer.h.
Referenced by analyze().
Definition at line 25 of file RazorVarAnalyzer.h.
Referenced by analyze().
Definition at line 26 of file RazorVarAnalyzer.h.
Referenced by analyze().
MonitorElement* RazorVarAnalyzer::m_rsqMREle [private] |
Definition at line 42 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMREleEle [private] |
Definition at line 40 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMREleMJ [private] |
Definition at line 47 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMREleMu [private] |
Definition at line 38 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRFullyInc [private] |
Definition at line 29 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRHad [private] |
Definition at line 43 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRHadMJ [private] |
Definition at line 48 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRInc10J [private] |
Definition at line 33 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRInc12J [private] |
Definition at line 34 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRInc14J [private] |
Definition at line 35 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRInc4J [private] |
Definition at line 30 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRInc6J [private] |
Definition at line 31 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRInc8J [private] |
Definition at line 32 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRMu [private] |
Definition at line 41 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRMuMJ [private] |
Definition at line 46 of file RazorVarAnalyzer.h.
MonitorElement* RazorVarAnalyzer::m_rsqMRMuMu [private] |
Definition at line 39 of file RazorVarAnalyzer.h.
1.7.3