|
|
|
#include <CMSDAS11DijetAnalyzer.h>
Public Member Functions | |
| void | analyze (const edm::Event &, const edm::EventSetup &) |
| virtual void | beginJob () |
| CMSDAS11DijetAnalyzer (const edm::ParameterSet &) | |
| virtual void | endJob (void) |
| virtual | ~CMSDAS11DijetAnalyzer () |
Static Public Member Functions | |
| static bool | compare_JetPt (const reco::CaloJet &jet1, const reco::CaloJet &jet2) |
Private Attributes | |
| TH1D * | hCorDijetMass |
| TH1D * | hCorDijetXsec |
| TH1D * | hDijetDeltaEta |
| TH1D * | hDijetDeltaPhi |
| TH1D * | hInnerDijetMass |
| TH1D * | hJet1EMF |
| TH1D * | hJet1Eta |
| TH1D * | hJet1Phi |
| TH1D * | hJet1Pt |
| TH1D * | hJet2EMF |
| TH1D * | hJet2Eta |
| TH1D * | hJet2Phi |
| TH1D * | hJet2Pt |
| TH1D * | hJetCorrPt |
| TH1D * | hJetEMF |
| TH1D * | hJetEta |
| TH1D * | hJetPhi |
| TH1D * | hJetRawPt |
| TH1D * | hOuterDijetMass |
| TH1D * | hRawDijetMass |
| TH1D * | hVertexZ |
| double | innerDeltaEta |
| double | JESbias |
| std::string | jetCorrections |
| edm::InputTag | jetSrc |
| double | outerDeltaEta |
| edm::InputTag | vertexSrc |
Definition at line 20 of file CMSDAS11DijetAnalyzer.h.
| CMSDAS11DijetAnalyzer::CMSDAS11DijetAnalyzer | ( | const edm::ParameterSet & | params | ) |
Definition at line 20 of file CMSDAS11DijetAnalyzer.cc.
References hCorDijetMass, hDijetDeltaEta, hDijetDeltaPhi, hInnerDijetMass, hJet1EMF, hJet1Eta, hJet1Phi, hJet1Pt, hJet2EMF, hJet2Eta, hJet2Phi, hJet2Pt, hJetCorrPt, hJetEMF, hJetEta, hJetPhi, hJetRawPt, hOuterDijetMass, hVertexZ, and NBINS.
: edm::EDAnalyzer(), jetSrc(params.getParameter<edm::InputTag>("jetSrc")), vertexSrc(params.getParameter<edm::InputTag>("vertexSrc")), jetCorrections(params.getParameter<std::string>("jetCorrections")), innerDeltaEta(params.getParameter<double>("innerDeltaEta")), outerDeltaEta(params.getParameter<double>("outerDeltaEta")), JESbias(params.getParameter<double>("JESbias")) { // setup file service edm::Service<TFileService> fs; const int NBINS=36; Double_t BOUNDARIES[NBINS] = { 220, 244, 270, 296, 325, 354, 386, 419, 453, 489, 526, 565, 606, 649, 693, 740, 788, 838, 890, 944, 1000, 1058, 1118, 1181, 1246, 1313, 1383, 1455, 1530, 1607, 1687, 1770, 1856, 1945, 2037, 2132 }; // setup histograms hVertexZ = fs->make<TH1D>("hVertexZ", "Z position of the Vertex",50,-20,20); hJetRawPt = fs->make<TH1D>("hJetRawPt","Raw Jet Pt",50,0,1000); hJetCorrPt = fs->make<TH1D>("hJetCorrPt","Corrected Jet Pt",50,0,1000); hJet1Pt = fs->make<TH1D>("hJet1Pt","Corrected Jet1 Pt",50,0,1000); hJet2Pt = fs->make<TH1D>("hJet2Pt","Corrected Jet2 Pt",50,0,1000); hJetEta = fs->make<TH1D>("hJetEta","Corrected Jet Eta", 10,-5,5); hJet1Eta = fs->make<TH1D>("hJet1Eta","Corrected Jet1 Eta",10,-5,5); hJet2Eta = fs->make<TH1D>("hJet2Eta","Corrected Jet2 Eta",10,-5,5); hJetPhi = fs->make<TH1D>("hJetPhi","Corrected Jet Phi", 10,-3.1415,3.1415); hJet1Phi = fs->make<TH1D>("hJet1Phi","Corrected Jet1 Phi",10,-3.1415,3.1415); hJet2Phi = fs->make<TH1D>("hJet2Phi","Corrected Jet2 Phi",10,-3.1415,3.1415); hJetEMF = fs->make<TH1D>("hJetEMF","EM Fraction of Jets",50,0,1); hJet1EMF = fs->make<TH1D>("hJet1EMF","EM Fraction of Jet1",50,0,1); hJet2EMF = fs->make<TH1D>("hJet2EMF","EM Fraction of Jet2",50,0,1); hCorDijetMass = fs->make<TH1D>("hCorDijetMass","Corrected Dijet Mass",NBINS-1,BOUNDARIES); hDijetDeltaPhi= fs->make<TH1D>("hDijetDeltaPhi","Dijet |#Delta #phi|",50,0,3.1415); hDijetDeltaEta= fs->make<TH1D>("hDijetDeltaEta","Dijet |#Delta #eta|",50,0,6); hInnerDijetMass = fs->make<TH1D>("hInnerDijetMass","Corrected Inner Dijet Mass",NBINS-1,BOUNDARIES); hOuterDijetMass = fs->make<TH1D>("hOuterDijetMass","Corrected Outer Dijet Mass",NBINS-1,BOUNDARIES); }
| virtual CMSDAS11DijetAnalyzer::~CMSDAS11DijetAnalyzer | ( | ) | [inline, virtual] |
Definition at line 24 of file CMSDAS11DijetAnalyzer.h.
{}
| void CMSDAS11DijetAnalyzer::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 69 of file CMSDAS11DijetAnalyzer.cc.
References compare_JetPt(), JetCorrector::correction(), gather_cfg::cout, reco::CaloJet::emEnergyFraction(), eta(), reco::LeafCandidate::eta(), edm::Event::getByLabel(), JetCorrector::getJetCorrector(), hCorDijetMass, hDijetDeltaEta, hDijetDeltaPhi, hInnerDijetMass, hJet1EMF, hJet1Eta, hJet1Phi, hJet1Pt, hJet2EMF, hJet2Eta, hJet2Phi, hJet2Pt, hJetCorrPt, hJetEMF, hJetEta, hJetPhi, hJetRawPt, hOuterDijetMass, hVertexZ, innerDeltaEta, edm::HandleBase::isValid(), JESbias, metsig::jet, jetCorrections, jetSrc, reco::Vertex::ndof(), outerDeltaEta, reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), phi, reco::Vertex::position(), reco::LeafCandidate::pt(), reco::Jet::scaleEnergy(), python::multivaluedict::sort(), vertexSrc, and reco::Vertex::z().
{
double mWeight;
edm::Handle<GenEventInfoProduct> hEventInfo;
iEvent.getByLabel("generator", hEventInfo);
if (hEventInfo.isValid())
mWeight = hEventInfo->weight();
else
mWeight = 1.;
// Get event ID information
//int nrun=iEvent.id().run();
//int nlumi=iEvent.luminosityBlock();
//int nevent=iEvent.id().event();
// Get Primary Vertex Information
// magic to get the vertices from EDM
edm::Handle< std::vector<reco::Vertex> > vertices_h;
iEvent.getByLabel(vertexSrc, vertices_h);
if (!vertices_h.isValid()) {
std::cout<<"Didja hear the one about the empty vertex collection?\n";
return;
}
// require in the event that there is at least one reconstructed vertex
if(vertices_h->size()<=0) return;
// pick the first (i.e. highest sum pt) verte
const reco::Vertex* theVertex=&(vertices_h->front());
// require that the vertex meets certain criteria
if(theVertex->ndof()<5) return;
if(fabs(theVertex->z())>24.0) return;
if(fabs(theVertex->position().rho())>2.0) return;
// Get Jet Information
// magic to get the jets from EDM
edm::Handle<reco::CaloJetCollection> jets_h;
iEvent.getByLabel(jetSrc, jets_h);
if (!jets_h.isValid()) {
std::cout<<"Didja hear the one about the empty jet collection?\n";
return;
}
// magic to get the jet energy corrections
const JetCorrector* corrector = JetCorrector::getJetCorrector(jetCorrections,iSetup);
// collection of selected jets
std::vector<reco::CaloJet> selectedJets;
// loop over the jet collection
for(reco::CaloJetCollection::const_iterator j_it = jets_h->begin(); j_it!=jets_h->end(); j_it++) {
reco::CaloJet jet = *j_it;
// calculate and apply the correction
double scale = corrector->correction(jet.p4());
// Introduce a purposeful bias to the correction, to show what happens
scale *= JESbias;
// fill the histograms
hJetRawPt->Fill(jet.pt(),mWeight);
hJetEta->Fill(jet.eta(),mWeight);
hJetPhi->Fill(jet.phi(),mWeight);
hJetEMF->Fill(jet.emEnergyFraction(),mWeight);
// correct the jet energy
jet.scaleEnergy(scale);
// now fill the correct jet pt after the correction
hJetCorrPt->Fill(jet.pt(),mWeight);
// put the selected jets into a collection
selectedJets.push_back(jet);
}
// require at least two jets to continue
if(selectedJets.size()<2) return;
//sort by corrected pt (not the same order as raw pt, sometimes)
sort(selectedJets.begin(), selectedJets.end(), compare_JetPt);
// select high pt, central, non-noise-like jets
if (selectedJets[0].pt()<50.0) return;
if (fabs(selectedJets[0].eta())>2.5) return;
if (selectedJets[0].emEnergyFraction()<0.01) return;
if (selectedJets[1].pt()<50.0) return;
if (fabs(selectedJets[1].eta())>2.5) return;
if (selectedJets[1].emEnergyFraction()<0.01) return;
// fill histograms for the jets in our dijets, only
hJet1Pt ->Fill(selectedJets[0].pt(),mWeight);
hJet1Eta->Fill(selectedJets[0].eta(),mWeight);
hJet1Phi->Fill(selectedJets[0].phi(),mWeight);
hJet1EMF->Fill(selectedJets[0].emEnergyFraction(),mWeight);
hJet2Pt ->Fill(selectedJets[1].pt(),mWeight);
hJet2Eta->Fill(selectedJets[1].eta(),mWeight);
hJet2Phi->Fill(selectedJets[1].phi(),mWeight);
hJet2EMF->Fill(selectedJets[1].emEnergyFraction(),mWeight);
//Get the mass of the two leading jets. Needs their 4-vectors...
double corMass = (selectedJets[0].p4()+selectedJets[1].p4()).M();
double deltaEta = fabs(selectedJets[0].eta()-selectedJets[1].eta());
if (corMass < 489) return;
if (deltaEta > 1.3) return;
hCorDijetMass->Fill(corMass,mWeight);
hDijetDeltaPhi->Fill(fabs(selectedJets[0].phi()-selectedJets[1].phi()) ,mWeight);
hDijetDeltaEta->Fill(deltaEta ,mWeight);
//Fill the inner and outer dijet mass spectra, to make the ratio from
if (deltaEta < innerDeltaEta) hInnerDijetMass->Fill(corMass,mWeight);
else if (deltaEta < outerDeltaEta) hOuterDijetMass->Fill(corMass,mWeight);
hVertexZ->Fill(theVertex->z(),mWeight);
return;
}
| virtual void CMSDAS11DijetAnalyzer::beginJob | ( | void | ) | [inline, virtual] |
| static bool CMSDAS11DijetAnalyzer::compare_JetPt | ( | const reco::CaloJet & | jet1, |
| const reco::CaloJet & | jet2 | ||
| ) | [inline, static] |
Definition at line 28 of file CMSDAS11DijetAnalyzer.h.
References reco::LeafCandidate::pt().
Referenced by analyze().
| void CMSDAS11DijetAnalyzer::endJob | ( | void | ) | [virtual] |
TH1D* CMSDAS11DijetAnalyzer::hCorDijetMass [private] |
Definition at line 50 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hCorDijetXsec [private] |
Definition at line 51 of file CMSDAS11DijetAnalyzer.h.
TH1D* CMSDAS11DijetAnalyzer::hDijetDeltaEta [private] |
Definition at line 61 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hDijetDeltaPhi [private] |
Definition at line 60 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hInnerDijetMass [private] |
Definition at line 63 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet1EMF [private] |
Definition at line 55 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet1Eta [private] |
Definition at line 53 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet1Phi [private] |
Definition at line 54 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet1Pt [private] |
Definition at line 52 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet2EMF [private] |
Definition at line 59 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet2Eta [private] |
Definition at line 57 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet2Phi [private] |
Definition at line 58 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJet2Pt [private] |
Definition at line 56 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJetCorrPt [private] |
Definition at line 43 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJetEMF [private] |
Definition at line 47 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJetEta [private] |
Definition at line 45 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJetPhi [private] |
Definition at line 46 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hJetRawPt [private] |
Definition at line 44 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hOuterDijetMass [private] |
Definition at line 64 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
TH1D* CMSDAS11DijetAnalyzer::hRawDijetMass [private] |
Definition at line 49 of file CMSDAS11DijetAnalyzer.h.
TH1D* CMSDAS11DijetAnalyzer::hVertexZ [private] |
Definition at line 42 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze(), and CMSDAS11DijetAnalyzer().
double CMSDAS11DijetAnalyzer::innerDeltaEta [private] |
Definition at line 37 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze().
double CMSDAS11DijetAnalyzer::JESbias [private] |
Definition at line 39 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze().
std::string CMSDAS11DijetAnalyzer::jetCorrections [private] |
Definition at line 36 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze().
edm::InputTag CMSDAS11DijetAnalyzer::jetSrc [private] |
Definition at line 34 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze().
double CMSDAS11DijetAnalyzer::outerDeltaEta [private] |
Definition at line 38 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze().
Definition at line 35 of file CMSDAS11DijetAnalyzer.h.
Referenced by analyze().
1.7.3