#include <JetAnalyzer.h>
DQM monitoring source for Calo Jets
Definition at line 37 of file JetAnalyzer.h.
JetAnalyzer::JetAnalyzer | ( | const edm::ParameterSet & | pSet | ) |
Constructor.
Definition at line 22 of file JetAnalyzer.cc.
References ExpressReco_HICollisions_FallBack::parameters.
{ parameters = pSet; _leadJetFlag = 0; _JetLoPass = 0; _JetHiPass = 0; _ptThreshold = 5.; _n90HitsMin =0; _fHPDMax=1.; _resEMFMin=0.; _n90HitsMinLoose =0; _fHPDMaxLoose=1.; _resEMFMinLoose=0.; _n90HitsMinTight =0; _fHPDMaxTight=1.; _resEMFMinTight=0.; _sigmaEtaMinTight=-999.; _sigmaPhiMinTight=-999.; }
JetAnalyzer::~JetAnalyzer | ( | ) | [virtual] |
void JetAnalyzer::analyze | ( | const edm::Event & | iEvent, |
const edm::EventSetup & | iSetup, | ||
const reco::CaloJetCollection & | caloJets | ||
) |
Get the analysis.
Definition at line 234 of file JetAnalyzer.cc.
References corr, diffTreeTool::diff, relval_parameters_module::energy, eta(), metsig::jet, LogTrace, L1TEmulatorMonitor_cff::p, phi, ExpressReco_HICollisions_FallBack::pt, and mathSSE::sqrt().
{ int numofjets=0; double fstPhi=0.; double sndPhi=0.; double diff = 0.; double corr = 0.; double dphi = -999. ; bool thiscleaned=false; bool Loosecleaned=false; bool Tightcleaned=false; bool thisemfclean=true; bool emfcleanLoose=true; bool emfcleanTight=true; if(makedijetselection==1){ //Dijet selection - careful: the pT is uncorrected! if(makedijetselection==1 && caloJets.size()>=2){ double dphiDJ = -999. ; bool emfcleanLooseFirstJet=true; bool emfcleanLooseSecondJet=true; bool emfcleanTightFirstJet=true; bool emfcleanTightSecondJet=true; bool LoosecleanedFirstJet = false; bool LoosecleanedSecondJet = false; bool TightcleanedFirstJet = false; bool TightcleanedSecondJet = false; //both jets pass pt threshold if ((caloJets.at(0)).pt() > _ptThreshold && (caloJets.at(1)).pt() > _ptThreshold ) { if(fabs((caloJets.at(0)).eta())<3. && fabs((caloJets.at(1)).eta())<3. ){ //calculate dphi dphiDJ = fabs((caloJets.at(0)).phi()-(caloJets.at(1)).phi()); if (dphiDJ > 3.14) dphiDJ=fabs(dphiDJ -6.28 ); //fill DPhi histo (before cutting) if (mDPhi) mDPhi->Fill (dphiDJ); //dphi cut if(fabs(dphiDJ)>2.1){ //JetID emfcleanLooseFirstJet=true; emfcleanTightFirstJet=true; emfcleanLooseSecondJet=true; emfcleanTightSecondJet=true; jetID->calculate(iEvent, (caloJets.at(0))); if(jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(0)).eta())<2.6) emfcleanLooseFirstJet=false; if(jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(0)).eta())<2.6) emfcleanTightFirstJet=false; if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseFirstJet) LoosecleanedFirstJet=true; if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt((caloJets.at(0)).etaetaMoment())>_sigmaEtaMinTight && sqrt((caloJets.at(0)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightFirstJet) TightcleanedFirstJet=true; //fill the JID variables histograms BEFORE you cut on them if (mN90Hits) mN90Hits->Fill (jetID->n90Hits()); if (mfHPD) mfHPD->Fill (jetID->fHPD()); if (mresEMF) mresEMF->Fill (jetID->restrictedEMF()); if (mfRBX) mfRBX->Fill (jetID->fRBX()); jetID->calculate(iEvent, (caloJets.at(1))); if(jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(1)).eta())<2.6) emfcleanLooseSecondJet=false; if(jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(1)).eta())<2.6) emfcleanTightSecondJet=false; if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseSecondJet) LoosecleanedSecondJet=true; if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt((caloJets.at(1)).etaetaMoment())>_sigmaEtaMinTight && sqrt((caloJets.at(1)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightSecondJet) TightcleanedSecondJet=true; //fill the JID variables histograms BEFORE you cut on them if (mN90Hits) mN90Hits->Fill (jetID->n90Hits()); if (mfHPD) mfHPD->Fill (jetID->fHPD()); if (mresEMF) mresEMF->Fill (jetID->restrictedEMF()); if (mfRBX) mfRBX->Fill (jetID->fRBX()); if(LoosecleanedFirstJet && LoosecleanedSecondJet) { //only if both jets are (loose) cleaned //fill histos for first jet if (mPt) mPt->Fill ((caloJets.at(0)).pt()); if (mEta) mEta->Fill ((caloJets.at(0)).eta()); if (mPhi) mPhi->Fill ((caloJets.at(0)).phi()); if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(0)).eta(),(caloJets.at(0)).phi()); if (mConstituents) mConstituents->Fill ((caloJets.at(0)).nConstituents()); if (mHFrac) mHFrac->Fill ((caloJets.at(0)).energyFractionHadronic()); if (mEFrac) mEFrac->Fill ((caloJets.at(0)).emEnergyFraction()); if (mE) mE->Fill ((caloJets.at(0)).energy()); if (mP) mP->Fill ((caloJets.at(0)).p()); if (mMass) mMass->Fill ((caloJets.at(0)).mass()); if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ((caloJets.at(0)).maxEInEmTowers()); if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(0)).maxEInHadTowers()); if (mN90Hits) mN90Hits->Fill (jetID->n90Hits()); if (mfHPD) mfHPD->Fill (jetID->fHPD()); if (mresEMF) mresEMF->Fill (jetID->restrictedEMF()); if (mfRBX) mfRBX->Fill (jetID->fRBX()); //sigmaeta and sigmaphi only used in the tight selection. //fill the histos for them AFTER the loose selection if (msigmaEta) msigmaEta->Fill(sqrt((caloJets.at(0)).etaetaMoment())); if (msigmaPhi) msigmaPhi->Fill(sqrt((caloJets.at(0)).phiphiMoment())); //fill histos for second jet if (mPt) mPt->Fill ((caloJets.at(1)).pt()); if (mEta) mEta->Fill ((caloJets.at(1)).eta()); if (mPhi) mPhi->Fill ((caloJets.at(1)).phi()); if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(1)).eta(),(caloJets.at(1)).phi()); if (mConstituents) mConstituents->Fill ((caloJets.at(1)).nConstituents()); if (mHFrac) mHFrac->Fill ((caloJets.at(1)).energyFractionHadronic()); if (mEFrac) mEFrac->Fill ((caloJets.at(1)).emEnergyFraction()); if (mE) mE->Fill ((caloJets.at(1)).energy()); if (mP) mP->Fill ((caloJets.at(1)).p()); if (mMass) mMass->Fill ((caloJets.at(1)).mass()); if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ((caloJets.at(1)).maxEInEmTowers()); if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(1)).maxEInHadTowers()); //sigmaeta and sigmaphi only used in the tight selection. //fill the histos for them AFTER the loose selection if (msigmaEta) msigmaEta->Fill(sqrt((caloJets.at(1)).etaetaMoment())); if (msigmaPhi) msigmaPhi->Fill(sqrt((caloJets.at(1)).phiphiMoment())); } //let's see how many of these jets passed the JetID cleaning if(fillJIDPassFrac==1) { if(LoosecleanedFirstJet) { mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),1.); mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.); } else { mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),0.); mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.); } if(LoosecleanedSecondJet) { mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),1.); mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.); } else { mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),0.); mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.); } //TIGHT JID if(TightcleanedFirstJet) { mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),1.); mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.); } else { mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),0.); mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.); } if(TightcleanedSecondJet) { mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),1.); mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.); } else { mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),0.); mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.); } }//if fillJIDPassFrac } } } } } else { for (reco::CaloJetCollection::const_iterator jet = caloJets.begin(); jet!=caloJets.end(); ++jet) { LogTrace(jetname)<<"[JetAnalyzer] Analyze Calo Jet"; Loosecleaned=false; Tightcleaned=false; if (jet == caloJets.begin()) { fstPhi = jet->phi(); _leadJetFlag = 1; } else { _leadJetFlag = 0; } if (jet == (caloJets.begin()+1)) sndPhi = jet->phi(); //jetID jetID->calculate(iEvent, *jet); //minimal (uncorrected!) pT cut if (jet->pt() > _ptThreshold) { if (msigmaEta) msigmaEta->Fill(sqrt(jet->etaetaMoment())); if (msigmaPhi) msigmaPhi->Fill(sqrt(jet->phiphiMoment())); //cleaning to use for filling histograms thisemfclean=true; if(jetID->restrictedEMF()<_resEMFMin && fabs(jet->eta())<2.6) thisemfclean=false; if(jetID->n90Hits()>=_n90HitsMin && jetID->fHPD()<_fHPDMax && thisemfclean) thiscleaned=true; //loose and tight cleaning, used to fill the JetIDPAssFraction histos if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLoose) Loosecleaned=true; if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt(jet->etaetaMoment())>_sigmaEtaMinTight && sqrt(jet->phiphiMoment())>_sigmaPhiMinTight && emfcleanTight) Tightcleaned=true; if(fillJIDPassFrac==1) { if(Loosecleaned) { mLooseJIDPassFractionVSeta->Fill(jet->eta(),1.); mLooseJIDPassFractionVSpt->Fill(jet->pt(),1.); } else { mLooseJIDPassFractionVSeta->Fill(jet->eta(),0.); mLooseJIDPassFractionVSpt->Fill(jet->pt(),0.); } //TIGHT if(Tightcleaned) { mTightJIDPassFractionVSeta->Fill(jet->eta(),1.); mTightJIDPassFractionVSpt->Fill(jet->pt(),1.); } else { mTightJIDPassFractionVSeta->Fill(jet->eta(),0.); mTightJIDPassFractionVSpt->Fill(jet->pt(),0.); } } //eventually we could define the "cleaned" flag differently for e.g. HF if(thiscleaned) { numofjets++ ; jetME->Fill(1); // Leading jet // Histograms are filled once per event if (_leadJetFlag == 1) { if (mEtaFirst) mEtaFirst->Fill (jet->eta()); if (mPhiFirst) mPhiFirst->Fill (jet->phi()); if (mEFirst) mEFirst->Fill (jet->energy()); if (mPtFirst) mPtFirst->Fill (jet->pt()); } // --- Passed the low pt jet trigger if (_JetLoPass == 1) { if (fabs(jet->eta()) <= 1.3) { if (mPt_Barrel_Lo) mPt_Barrel_Lo->Fill(jet->pt()); if (mEta_Lo) mEta_Lo->Fill(jet->eta()); if (mPhi_Barrel_Lo) mPhi_Barrel_Lo->Fill(jet->phi()); if (mConstituents_Barrel_Lo) mConstituents_Barrel_Lo->Fill(jet->nConstituents()); if (mHFrac_Barrel_Lo) mHFrac_Barrel_Lo->Fill(jet->energyFractionHadronic()); } if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) { if (mPt_EndCap_Lo) mPt_EndCap_Lo->Fill(jet->pt()); if (mEta_Lo) mEta_Lo->Fill(jet->eta()); if (mPhi_EndCap_Lo) mPhi_EndCap_Lo->Fill(jet->phi()); if (mConstituents_EndCap_Lo) mConstituents_EndCap_Lo->Fill(jet->nConstituents()); if (mHFrac_EndCap_Lo) mHFrac_EndCap_Lo->Fill(jet->energyFractionHadronic()); } if (fabs(jet->eta()) > 3.0) { if (mPt_Forward_Lo) mPt_Forward_Lo->Fill(jet->pt()); if (mEta_Lo) mEta_Lo->Fill(jet->eta()); if (mPhi_Forward_Lo) mPhi_Forward_Lo->Fill(jet->phi()); if (mConstituents_Forward_Lo) mConstituents_Forward_Lo->Fill(jet->nConstituents()); if (mHFrac_Forward_Lo) mHFrac_Forward_Lo->Fill(jet->energyFractionHadronic()); } if (mEta_Lo) mEta_Lo->Fill (jet->eta()); if (mPhi_Lo) mPhi_Lo->Fill (jet->phi()); if (mPt_Lo) mPt_Lo->Fill (jet->pt()); } // --- Passed the high pt jet trigger if (_JetHiPass == 1) { if (fabs(jet->eta()) <= 1.3) { if (mPt_Barrel_Hi) mPt_Barrel_Hi->Fill(jet->pt()); if (mEta_Hi) mEta_Hi->Fill(jet->eta()); if (mPhi_Barrel_Hi) mPhi_Barrel_Hi->Fill(jet->phi()); if (mConstituents_Barrel_Hi) mConstituents_Barrel_Hi->Fill(jet->nConstituents()); if (mHFrac_Barrel_Hi) mHFrac_Barrel_Hi->Fill(jet->energyFractionHadronic()); } if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) { if (mPt_EndCap_Hi) mPt_EndCap_Hi->Fill(jet->pt()); if (mEta_Hi) mEta_Hi->Fill(jet->eta()); if (mPhi_EndCap_Hi) mPhi_EndCap_Hi->Fill(jet->phi()); if (mConstituents_EndCap_Hi) mConstituents_EndCap_Hi->Fill(jet->nConstituents()); if (mHFrac_EndCap_Hi) mHFrac_EndCap_Hi->Fill(jet->energyFractionHadronic()); } if (fabs(jet->eta()) > 3.0) { if (mPt_Forward_Hi) mPt_Forward_Hi->Fill(jet->pt()); if (mEta_Hi) mEta_Hi->Fill(jet->eta()); if (mPhi_Forward_Hi) mPhi_Forward_Hi->Fill(jet->phi()); if (mConstituents_Forward_Hi) mConstituents_Forward_Hi->Fill(jet->nConstituents()); if (mHFrac_Forward_Hi) mHFrac_Forward_Hi->Fill(jet->energyFractionHadronic()); } if (mEta_Hi) mEta_Hi->Fill (jet->eta()); if (mPhi_Hi) mPhi_Hi->Fill (jet->phi()); if (mPt_Hi) mPt_Hi->Fill (jet->pt()); } if (mPt) mPt->Fill (jet->pt()); if (mPt_1) mPt_1->Fill (jet->pt()); if (mPt_2) mPt_2->Fill (jet->pt()); if (mPt_3) mPt_3->Fill (jet->pt()); if (mEta) mEta->Fill (jet->eta()); if (mPhi) mPhi->Fill (jet->phi()); if (mPhiVSEta) mPhiVSEta->Fill(jet->eta(),jet->phi()); if (mConstituents) mConstituents->Fill (jet->nConstituents()); if (mHFrac) mHFrac->Fill (jet->energyFractionHadronic()); if (mEFrac) mEFrac->Fill (jet->emEnergyFraction()); if (fabs(jet->eta()) <= 1.3) { if (mPt_Barrel) mPt_Barrel->Fill (jet->pt()); if (mPhi_Barrel) mPhi_Barrel->Fill (jet->phi()); if (mE_Barrel) mE_Barrel->Fill (jet->energy()); } if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) { if (mPt_EndCap) mPt_EndCap->Fill (jet->pt()); if (mPhi_EndCap) mPhi_EndCap->Fill (jet->phi()); if (mE_EndCap) mE_EndCap->Fill (jet->energy()); } if (fabs(jet->eta()) > 3.0) { if (mPt_Forward) mPt_Forward->Fill (jet->pt()); if (mPhi_Forward) mPhi_Forward->Fill (jet->phi()); if (mE_Forward) mE_Forward->Fill (jet->energy()); } if (mE) mE->Fill (jet->energy()); if (mP) mP->Fill (jet->p()); if (mMass) mMass->Fill (jet->mass()); if (mMaxEInEmTowers) mMaxEInEmTowers->Fill (jet->maxEInEmTowers()); if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->maxEInHadTowers()); if (mHadEnergyInHO) mHadEnergyInHO->Fill (jet->hadEnergyInHO()); if (mHadEnergyInHB) mHadEnergyInHB->Fill (jet->hadEnergyInHB()); if (mHadEnergyInHF) mHadEnergyInHF->Fill (jet->hadEnergyInHF()); if (mHadEnergyInHE) mHadEnergyInHE->Fill (jet->hadEnergyInHE()); if (mEmEnergyInEB) mEmEnergyInEB->Fill (jet->emEnergyInEB()); if (mEmEnergyInEE) mEmEnergyInEE->Fill (jet->emEnergyInEE()); if (mEmEnergyInHF) mEmEnergyInHF->Fill (jet->emEnergyInHF()); if (mN90Hits) mN90Hits->Fill (jetID->n90Hits()); if (mfHPD) mfHPD->Fill (jetID->fHPD()); if (mresEMF) mresEMF->Fill (jetID->restrictedEMF()); if (mfRBX) mfRBX->Fill (jetID->fRBX()); //calculate correctly the dphi if(numofjets>1) { diff = fabs(fstPhi - sndPhi); corr = 2*acos(-1.) - diff; if(diff < acos(-1.)) { dphi = diff; } else { dphi = corr; } } } }//pt cut } if (mNJets) mNJets->Fill (numofjets); if (mDPhi && dphi>-998.) mDPhi->Fill (dphi); }//not dijet }
void JetAnalyzer::beginJob | ( | DQMStore * | dbe | ) | [virtual] |
Inizialize parameters for histo binning.
Implements JetAnalyzerBase.
Definition at line 48 of file JetAnalyzer.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), jptDQMConfig_cff::eMax, ExpressReco_HICollisions_FallBack::etaBin, ExpressReco_HICollisions_FallBack::etaMax, ExpressReco_HICollisions_FallBack::etaMin, LogTrace, ExpressReco_HICollisions_FallBack::parameters, ExpressReco_HICollisions_FallBack::phiBin, ExpressReco_HICollisions_FallBack::phiMax, ExpressReco_HICollisions_FallBack::phiMin, ExpressReco_HICollisions_FallBack::pMax, ExpressReco_HICollisions_FallBack::pMin, ExpressReco_HICollisions_FallBack::ptMax, ExpressReco_HICollisions_FallBack::ptMin, MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().
{ jetname = "jetAnalyzer"; LogTrace(jetname)<<"[JetAnalyzer] Parameters initialization"; dbe->setCurrentFolder("JetMET/Jet/"+_source); jetME = dbe->book1D("jetReco", "jetReco", 3, 1, 4); jetME->setBinLabel(1,"CaloJets",1); // jetID = new reco::helper::JetIDHelper(parameters.getParameter<ParameterSet>("JetIDParams")); // fillJIDPassFrac = parameters.getParameter<int>("fillJIDPassFrac"); makedijetselection = parameters.getParameter<int>("makedijetselection"); // monitoring of eta parameter etaBin = parameters.getParameter<int>("etaBin"); etaMin = parameters.getParameter<double>("etaMin"); etaMax = parameters.getParameter<double>("etaMax"); // monitoring of phi paramater phiBin = parameters.getParameter<int>("phiBin"); phiMin = parameters.getParameter<double>("phiMin"); phiMax = parameters.getParameter<double>("phiMax"); // monitoring of the transverse momentum ptBin = parameters.getParameter<int>("ptBin"); ptMin = parameters.getParameter<double>("ptMin"); ptMax = parameters.getParameter<double>("ptMax"); // eBin = parameters.getParameter<int>("eBin"); eMin = parameters.getParameter<double>("eMin"); eMax = parameters.getParameter<double>("eMax"); // pBin = parameters.getParameter<int>("pBin"); pMin = parameters.getParameter<double>("pMin"); pMax = parameters.getParameter<double>("pMax"); // _ptThreshold = parameters.getParameter<double>("ptThreshold"); _n90HitsMin = parameters.getParameter<int>("n90HitsMin"); _fHPDMax = parameters.getParameter<double>("fHPDMax"); _resEMFMin = parameters.getParameter<double>("resEMFMin"); _sigmaEtaMinTight = parameters.getParameter<double>("sigmaEtaMinTight"); _sigmaPhiMinTight = parameters.getParameter<double>("sigmaPhiMinTight"); _n90HitsMinLoose = parameters.getParameter<int>("n90HitsMinLoose"); _fHPDMaxLoose = parameters.getParameter<double>("fHPDMaxLoose"); _resEMFMinLoose = parameters.getParameter<double>("resEMFMinLoose"); _n90HitsMinTight = parameters.getParameter<int>("n90HitsMinTight"); _fHPDMaxTight = parameters.getParameter<double>("fHPDMaxTight"); _resEMFMinTight = parameters.getParameter<double>("resEMFMinTight"); // Generic Jet Parameters mPt = dbe->book1D("Pt", "Pt", ptBin, ptMin, ptMax); mEta = dbe->book1D("Eta", "Eta", etaBin, etaMin, etaMax); mPhi = dbe->book1D("Phi", "Phi", phiBin, phiMin, phiMax); mConstituents = dbe->book1D("Constituents", "# of Constituents", 100, 0, 100); mHFrac = dbe->book1D("HFrac", "HFrac", 120, -0.1, 1.1); mEFrac = dbe->book1D("EFrac", "EFrac", 120, -0.1, 1.1); // mE = dbe->book1D("E", "E", eBin, eMin, eMax); mP = dbe->book1D("P", "P", pBin, pMin, pMax); mMass = dbe->book1D("Mass", "Mass", 100, 0, 25); // mPhiVSEta = dbe->book2D("PhiVSEta", "PhiVSEta", 50, etaMin, etaMax, 24, phiMin, phiMax); if(makedijetselection!=1){ mPt_1 = dbe->book1D("Pt1", "Pt1", 100, 0, 100); mPt_2 = dbe->book1D("Pt2", "Pt2", 100, 0, 300); mPt_3 = dbe->book1D("Pt3", "Pt3", 100, 0, 5000); // Low and high pt trigger paths mPt_Lo = dbe->book1D("Pt_Lo", "Pt (Pass Low Pt Jet Trigger)", 100, 0, 100); mEta_Lo = dbe->book1D("Eta_Lo", "Eta (Pass Low Pt Jet Trigger)", etaBin, etaMin, etaMax); mPhi_Lo = dbe->book1D("Phi_Lo", "Phi (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax); mPt_Hi = dbe->book1D("Pt_Hi", "Pt (Pass Hi Pt Jet Trigger)", 100, 0, 300); mEta_Hi = dbe->book1D("Eta_Hi", "Eta (Pass Hi Pt Jet Trigger)", etaBin, etaMin, etaMax); mPhi_Hi = dbe->book1D("Phi_Hi", "Phi (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax); mNJets = dbe->book1D("NJets", "Number of Jets", 100, 0, 100); mPt_Barrel_Lo = dbe->book1D("Pt_Barrel_Lo", "Pt Barrel (Pass Low Pt Jet Trigger)", 100, 0, 100); mPhi_Barrel_Lo = dbe->book1D("Phi_Barrel_Lo", "Phi Barrel (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax); mConstituents_Barrel_Lo = dbe->book1D("Constituents_Barrel_Lo", "Constituents Barrel (Pass Low Pt Jet Trigger)", 100, 0, 100); mHFrac_Barrel_Lo = dbe->book1D("HFrac_Barrel_Lo", "HFrac Barrel (Pass Low Pt Jet Trigger)", 100, 0, 1); mPt_EndCap_Lo = dbe->book1D("Pt_EndCap_Lo", "Pt EndCap (Pass Low Pt Jet Trigger)", 100, 0, 100); mPhi_EndCap_Lo = dbe->book1D("Phi_EndCap_Lo", "Phi EndCap (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax); mConstituents_EndCap_Lo = dbe->book1D("Constituents_EndCap_Lo", "Constituents EndCap (Pass Low Pt Jet Trigger)", 100, 0, 100); mHFrac_EndCap_Lo = dbe->book1D("HFrac_Endcap_Lo", "HFrac EndCap (Pass Low Pt Jet Trigger)", 100, 0, 1); mPt_Forward_Lo = dbe->book1D("Pt_Forward_Lo", "Pt Forward (Pass Low Pt Jet Trigger)", 100, 0, 100); mPhi_Forward_Lo = dbe->book1D("Phi_Forward_Lo", "Phi Forward (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax); mConstituents_Forward_Lo = dbe->book1D("Constituents_Forward_Lo", "Constituents Forward (Pass Low Pt Jet Trigger)", 100, 0, 100); mHFrac_Forward_Lo = dbe->book1D("HFrac_Forward_Lo", "HFrac Forward (Pass Low Pt Jet Trigger)", 100, 0, 1); mPt_Barrel_Hi = dbe->book1D("Pt_Barrel_Hi", "Pt Barrel (Pass Hi Pt Jet Trigger)", 100, 0, 300); mPhi_Barrel_Hi = dbe->book1D("Phi_Barrel_Hi", "Phi Barrel (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax); mConstituents_Barrel_Hi = dbe->book1D("Constituents_Barrel_Hi", "Constituents Barrel (Pass Hi Pt Jet Trigger)", 100, 0, 100); mHFrac_Barrel_Hi = dbe->book1D("HFrac_Barrel_Hi", "HFrac Barrel (Pass Hi Pt Jet Trigger)", 100, 0, 1); mPt_EndCap_Hi = dbe->book1D("Pt_EndCap_Hi", "Pt EndCap (Pass Hi Pt Jet Trigger)", 100, 0, 300); mPhi_EndCap_Hi = dbe->book1D("Phi_EndCap_Hi", "Phi EndCap (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax); mConstituents_EndCap_Hi = dbe->book1D("Constituents_EndCap_Hi", "Constituents EndCap (Pass Hi Pt Jet Trigger)", 100, 0, 100); mHFrac_EndCap_Hi = dbe->book1D("HFrac_EndCap_Hi", "HFrac EndCap (Pass Hi Pt Jet Trigger)", 100, 0, 1); mPt_Forward_Hi = dbe->book1D("Pt_Forward_Hi", "Pt Forward (Pass Hi Pt Jet Trigger)", 100, 0, 300); mPhi_Forward_Hi = dbe->book1D("Phi_Forward_Hi", "Phi Forward (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax); mConstituents_Forward_Hi = dbe->book1D("Constituents_Forward_Hi", "Constituents Forward (Pass Hi Pt Jet Trigger)", 100, 0, 100); mHFrac_Forward_Hi = dbe->book1D("HFrac_Forward_Hi", "HFrac Forward (Pass Hi Pt Jet Trigger)", 100, 0, 1); mPhi_Barrel = dbe->book1D("Phi_Barrel", "Phi_Barrel", phiBin, phiMin, phiMax); mE_Barrel = dbe->book1D("E_Barrel", "E_Barrel", eBin, eMin, eMax); mPt_Barrel = dbe->book1D("Pt_Barrel", "Pt_Barrel", ptBin, ptMin, ptMax); mPhi_EndCap = dbe->book1D("Phi_EndCap", "Phi_EndCap", phiBin, phiMin, phiMax); mE_EndCap = dbe->book1D("E_EndCap", "E_EndCap", eBin, eMin, 2*eMax); mPt_EndCap = dbe->book1D("Pt_EndCap", "Pt_EndCap", ptBin, ptMin, ptMax); mPhi_Forward = dbe->book1D("Phi_Forward", "Phi_Forward", phiBin, phiMin, phiMax); mE_Forward = dbe->book1D("E_Forward", "E_Forward", eBin, eMin, 4*eMax); mPt_Forward = dbe->book1D("Pt_Forward", "Pt_Forward", ptBin, ptMin, ptMax); // Leading Jet Parameters mEtaFirst = dbe->book1D("EtaFirst", "EtaFirst", 100, -5, 5); mPhiFirst = dbe->book1D("PhiFirst", "PhiFirst", 70, -3.5, 3.5); mEFirst = dbe->book1D("EFirst", "EFirst", 100, 0, 1000); mPtFirst = dbe->book1D("PtFirst", "PtFirst", 100, 0, 500); if(fillJIDPassFrac==1){//fillJIDPassFrac defines a collection of cleaned jets, for which we will want to fill the cleaning passing fraction mLooseJIDPassFractionVSeta = dbe->bookProfile("LooseJIDPassFractionVSeta","LooseJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2); mLooseJIDPassFractionVSpt = dbe->bookProfile("LooseJIDPassFractionVSpt","LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2); mTightJIDPassFractionVSeta = dbe->bookProfile("TightJIDPassFractionVSeta","TightJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2); mTightJIDPassFractionVSpt = dbe->bookProfile("TightJIDPassFractionVSpt","TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2); } } // CaloJet specific mMaxEInEmTowers = dbe->book1D("MaxEInEmTowers", "MaxEInEmTowers", 100, 0, 100); mMaxEInHadTowers = dbe->book1D("MaxEInHadTowers", "MaxEInHadTowers", 100, 0, 100); if(makedijetselection!=1) { mHadEnergyInHO = dbe->book1D("HadEnergyInHO", "HadEnergyInHO", 100, 0, 10); mHadEnergyInHB = dbe->book1D("HadEnergyInHB", "HadEnergyInHB", 100, 0, 50); mHadEnergyInHF = dbe->book1D("HadEnergyInHF", "HadEnergyInHF", 100, 0, 50); mHadEnergyInHE = dbe->book1D("HadEnergyInHE", "HadEnergyInHE", 100, 0, 100); mEmEnergyInEB = dbe->book1D("EmEnergyInEB", "EmEnergyInEB", 100, 0, 50); mEmEnergyInEE = dbe->book1D("EmEnergyInEE", "EmEnergyInEE", 100, 0, 50); mEmEnergyInHF = dbe->book1D("EmEnergyInHF", "EmEnergyInHF", 120, -20, 100); } mDPhi = dbe->book1D("DPhi", "dPhi btw the two leading jets", 100, 0., acos(-1.)); //JetID variables mresEMF = dbe->book1D("resEMF", "resEMF", 50, 0., 1.); mN90Hits = dbe->book1D("N90Hits", "N90Hits", 50, 0., 50); mfHPD = dbe->book1D("fHPD", "fHPD", 50, 0., 1.); mfRBX = dbe->book1D("fRBX", "fRBX", 50, 0., 1.); msigmaEta = dbe->book1D("sigmaEta", "sigmaEta", 50, 0., 0.5); msigmaPhi = dbe->book1D("sigmaPhi", "sigmaPhi", 50, 0., 0.5); if(makedijetselection==1 && fillJIDPassFrac==1) { mLooseJIDPassFractionVSeta = dbe->bookProfile("LooseJIDPassFractionVSeta","LooseJIDPassFractionVSeta",50, -3., 3.,0.,1.2); mLooseJIDPassFractionVSpt = dbe->bookProfile("LooseJIDPassFractionVSpt","LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2); mTightJIDPassFractionVSeta = dbe->bookProfile("TightJIDPassFractionVSeta","TightJIDPassFractionVSeta",50, -3., 3.,0.,1.2); mTightJIDPassFractionVSpt = dbe->bookProfile("TightJIDPassFractionVSpt","TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2); } }
void JetAnalyzer::endJob | ( | void | ) |
int JetAnalyzer::getLeadJetFlag | ( | ) | [inline] |
void JetAnalyzer::setJetHiPass | ( | int | pass | ) | [inline] |
void JetAnalyzer::setJetLoPass | ( | int | pass | ) | [inline] |
void JetAnalyzer::setLeadJetFlag | ( | int | flag | ) | [inline] |
void JetAnalyzer::setSource | ( | std::string | source | ) | [inline] |
Definition at line 57 of file JetAnalyzer.h.
References _source, and LaserTracksInput_cfi::source.
double JetAnalyzer::_fHPDMax [private] |
Definition at line 94 of file JetAnalyzer.h.
double JetAnalyzer::_fHPDMaxLoose [private] |
Definition at line 99 of file JetAnalyzer.h.
double JetAnalyzer::_fHPDMaxTight [private] |
Definition at line 102 of file JetAnalyzer.h.
int JetAnalyzer::_JetHiPass [private] |
Definition at line 86 of file JetAnalyzer.h.
Referenced by setJetHiPass().
int JetAnalyzer::_JetLoPass [private] |
Definition at line 85 of file JetAnalyzer.h.
Referenced by setJetLoPass().
int JetAnalyzer::_leadJetFlag [private] |
Definition at line 87 of file JetAnalyzer.h.
Referenced by getLeadJetFlag(), and setLeadJetFlag().
int JetAnalyzer::_n90HitsMin [private] |
Definition at line 96 of file JetAnalyzer.h.
int JetAnalyzer::_n90HitsMinLoose [private] |
Definition at line 101 of file JetAnalyzer.h.
int JetAnalyzer::_n90HitsMinTight [private] |
Definition at line 104 of file JetAnalyzer.h.
double JetAnalyzer::_ptThreshold [private] |
Definition at line 89 of file JetAnalyzer.h.
double JetAnalyzer::_resEMFMin [private] |
Definition at line 95 of file JetAnalyzer.h.
double JetAnalyzer::_resEMFMinLoose [private] |
Definition at line 100 of file JetAnalyzer.h.
double JetAnalyzer::_resEMFMinTight [private] |
Definition at line 103 of file JetAnalyzer.h.
double JetAnalyzer::_sigmaEtaMinTight [private] |
Definition at line 105 of file JetAnalyzer.h.
double JetAnalyzer::_sigmaPhiMinTight [private] |
Definition at line 106 of file JetAnalyzer.h.
std::string JetAnalyzer::_source [private] |
Definition at line 81 of file JetAnalyzer.h.
Referenced by setSource().
int JetAnalyzer::_theend [private] |
Definition at line 88 of file JetAnalyzer.h.
int JetAnalyzer::eBin [private] |
Definition at line 121 of file JetAnalyzer.h.
double JetAnalyzer::eMax [private] |
Definition at line 123 of file JetAnalyzer.h.
double JetAnalyzer::eMin [private] |
Definition at line 122 of file JetAnalyzer.h.
int JetAnalyzer::etaBin [private] |
Definition at line 109 of file JetAnalyzer.h.
double JetAnalyzer::etaMax [private] |
Definition at line 111 of file JetAnalyzer.h.
double JetAnalyzer::etaMin [private] |
Definition at line 110 of file JetAnalyzer.h.
int JetAnalyzer::fillJIDPassFrac [private] |
Definition at line 98 of file JetAnalyzer.h.
reco::helper::JetIDHelper* JetAnalyzer::jetID [private] |
Definition at line 133 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::jetME [private] |
Definition at line 130 of file JetAnalyzer.h.
std::string JetAnalyzer::jetname [private] |
Definition at line 80 of file JetAnalyzer.h.
int JetAnalyzer::makedijetselection [private] |
Definition at line 91 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mConstituents [private] |
Definition at line 155 of file JetAnalyzer.h.
Definition at line 184 of file JetAnalyzer.h.
Definition at line 171 of file JetAnalyzer.h.
Definition at line 188 of file JetAnalyzer.h.
Definition at line 175 of file JetAnalyzer.h.
Definition at line 192 of file JetAnalyzer.h.
Definition at line 179 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mDPhi [private] |
Definition at line 205 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mE [private] |
Definition at line 201 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mE_Barrel [private] |
Definition at line 197 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mE_EndCap [private] |
Definition at line 198 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mE_Forward [private] |
Definition at line 199 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFirst [private] |
Definition at line 210 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEFrac [private] |
Definition at line 157 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEmEnergyInEB [private] |
Definition at line 221 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEmEnergyInEE [private] |
Definition at line 222 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEmEnergyInHF [private] |
Definition at line 223 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEta [private] |
Definition at line 153 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEta_Hi [private] |
Definition at line 242 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEta_Lo [private] |
Definition at line 238 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mEtaFirst [private] |
Definition at line 208 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mfHPD [private] |
Definition at line 227 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mfRBX [private] |
Definition at line 228 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHB [private] |
Definition at line 218 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHE [private] |
Definition at line 220 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHF [private] |
Definition at line 219 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHadEnergyInHO [private] |
Definition at line 217 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac [private] |
Definition at line 156 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Barrel_Hi [private] |
Definition at line 185 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Barrel_Lo [private] |
Definition at line 172 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_EndCap_Hi [private] |
Definition at line 189 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_EndCap_Lo [private] |
Definition at line 176 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Forward_Hi [private] |
Definition at line 193 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mHFrac_Forward_Lo [private] |
Definition at line 180 of file JetAnalyzer.h.
Definition at line 232 of file JetAnalyzer.h.
Definition at line 233 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mMass [private] |
Definition at line 203 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mMaxEInEmTowers [private] |
Definition at line 215 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mMaxEInHadTowers [private] |
Definition at line 216 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mN90Hits [private] |
Definition at line 226 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mNJets [private] |
Definition at line 204 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mP [private] |
Definition at line 202 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi [private] |
Definition at line 154 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Barrel [private] |
Definition at line 161 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Barrel_Hi [private] |
Definition at line 183 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Barrel_Lo [private] |
Definition at line 170 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_EndCap [private] |
Definition at line 164 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_EndCap_Hi [private] |
Definition at line 187 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_EndCap_Lo [private] |
Definition at line 174 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Forward [private] |
Definition at line 167 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Forward_Hi [private] |
Definition at line 191 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Forward_Lo [private] |
Definition at line 178 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Hi [private] |
Definition at line 243 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhi_Lo [private] |
Definition at line 239 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhiFirst [private] |
Definition at line 209 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPhiVSEta [private] |
Definition at line 158 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt [private] |
Definition at line 149 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_1 [private] |
Definition at line 150 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_2 [private] |
Definition at line 151 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_3 [private] |
Definition at line 152 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Barrel [private] |
Definition at line 160 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Barrel_Hi [private] |
Definition at line 182 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Barrel_Lo [private] |
Definition at line 169 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_EndCap [private] |
Definition at line 163 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_EndCap_Hi [private] |
Definition at line 186 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_EndCap_Lo [private] |
Definition at line 173 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Forward [private] |
Definition at line 166 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Forward_Hi [private] |
Definition at line 190 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Forward_Lo [private] |
Definition at line 177 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Hi [private] |
Definition at line 244 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPt_Lo [private] |
Definition at line 240 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mPtFirst [private] |
Definition at line 211 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::mresEMF [private] |
Definition at line 229 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::msigmaEta [private] |
Definition at line 230 of file JetAnalyzer.h.
MonitorElement* JetAnalyzer::msigmaPhi [private] |
Definition at line 231 of file JetAnalyzer.h.
Definition at line 234 of file JetAnalyzer.h.
Definition at line 235 of file JetAnalyzer.h.
edm::ParameterSet JetAnalyzer::parameters [private] |
Definition at line 78 of file JetAnalyzer.h.
int JetAnalyzer::pBin [private] |
Definition at line 125 of file JetAnalyzer.h.
int JetAnalyzer::phiBin [private] |
Definition at line 113 of file JetAnalyzer.h.
double JetAnalyzer::phiMax [private] |
Definition at line 115 of file JetAnalyzer.h.
double JetAnalyzer::phiMin [private] |
Definition at line 114 of file JetAnalyzer.h.
double JetAnalyzer::pMax [private] |
Definition at line 127 of file JetAnalyzer.h.
double JetAnalyzer::pMin [private] |
Definition at line 126 of file JetAnalyzer.h.
int JetAnalyzer::ptBin [private] |
Definition at line 117 of file JetAnalyzer.h.
double JetAnalyzer::ptMax [private] |
Definition at line 119 of file JetAnalyzer.h.
double JetAnalyzer::ptMin [private] |
Definition at line 118 of file JetAnalyzer.h.
Definition at line 83 of file JetAnalyzer.h.