|
|
|
#include <CaloJetTester.h>
Definition at line 21 of file CaloJetTester.h.
| CaloJetTester::CaloJetTester | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 48 of file CaloJetTester.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::book3D(), DQMStore::bookProfile(), DQMStore::bookProfile2D(), jptDQMConfig_cff::etaBins, ExpressReco_HICollisions_FallBack::etaMax, ExpressReco_HICollisions_FallBack::etaMin, edm::InputTag::label(), mAllGenJetsEta, mAllGenJetsPt, mCaloMET, mCaloMET_3000, mCaloMETPhi, mCaloMETSig, mCaloMETSig_3000, mCaloMEx, mCaloMEx_3000, mCaloMEy, mCaloMEy_3000, mCaloSumET, mCaloSumET_3000, mConstituents, mConstituents_80, mDelEta, mDelPhi, mDelPt, mDeltaE, mDeltaEta, mDeltaPhi, mE, mE_3000, mE_80, mEBEne, mEBTime, mEEEne, mEETime, mEFirst, mEFirst_3000, mEFirst_80, mEmEnergyInEB, mEmEnergyInEB_3000, mEmEnergyInEB_80, mEmEnergyInEE, mEmEnergyInEE_3000, mEmEnergyInEE_80, mEmEnergyInHF, mEmEnergyProfile, mEMJetEnergyProfile, mEmTiming, mEnergyFractionEm, mEnergyFractionHadronic, mEScale, mEScale_pt10, mEScaleFineBin, mEta, mEtaFineBin, mEtaFineBin1m, mEtaFineBin1p, mEtaFineBin2m, mEtaFineBin2p, mEtaFineBin3m, mEtaFineBin3p, mEtaFirst, mGenEta, mGenEtaFirst, mGenJetMatchEnergyFraction, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_3000, mGenPt_80, mHadEnergyInHB, mHadEnergyInHB_3000, mHadEnergyInHB_80, mHadEnergyInHE, mHadEnergyInHE_3000, mHadEnergyInHE_80, mHadEnergyInHF, mHadEnergyInHO, mHadEnergyInHO_3000, mHadEnergyInHO_80, mHadEnergyProfile, mHadJetEnergyProfile, mHadTiming, mHBEne, mHBTime, mHEEne, mHETime, mHFEne, mHFLong, mHFLong_3000, mHFLong_80, mHFShort, mHFShort_3000, mHFShort_80, mHFTime, mHFTotal, mHFTotal_3000, mHFTotal_80, mHOEne, mHOTime, mInputCollection, mJetEnergyProfile, mlinEScale, mMass, mMass_3000, mMass_80, mMatchedGenJetsEta, mMatchedGenJetsPt, mMaxEInEmTowers, mMaxEInHadTowers, mMjj, mMjj_3000, mN90, mNJets1, mNJets2, mNJetsEtaC, mNJetsEtaF, mOutputFile, mP, mP_3000, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_3000, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPthat_3000, mPthat_80, mpTScale1D_200_300, mpTScale1D_2700_3500, mpTScale1D_600_900, mpTScale1D_60_120, mpTScale1DB_200_300, mpTScale1DB_2700_3500, mpTScale1DB_600_900, mpTScale1DB_60_120, mpTScale1DE_200_300, mpTScale1DE_2700_3500, mpTScale1DE_600_900, mpTScale1DE_60_120, mpTScale1DF_200_300, mpTScale1DF_2700_3500, mpTScale1DF_600_900, mpTScale1DF_60_120, mpTScale_200_300_d, mpTScale_200_300_s, mpTScale_2700_3500_d, mpTScale_2700_3500_s, mpTScale_600_900_d, mpTScale_600_900_s, mpTScale_60_120_d, mpTScale_60_120_s, mpTScaleB_d, mpTScaleB_s, mpTScaleE_d, mpTScaleE_s, mpTScaleF_d, mpTScaleF_s, mReverseMatchEnergyFraction, mRMatch, mTurnOnEverything, numberofevents, and DQMStore::setCurrentFolder().
: mInputCollection (iConfig.getParameter<edm::InputTag>( "src" )), mInputGenCollection (iConfig.getParameter<edm::InputTag>( "srcGen" )), mOutputFile (iConfig.getUntrackedParameter<std::string>("outputFile", "")), mMatchGenPtThreshold (iConfig.getParameter<double>("genPtThreshold")), mGenEnergyFractionThreshold (iConfig.getParameter<double>("genEnergyFractionThreshold")), mReverseEnergyFractionThreshold (iConfig.getParameter<double>("reverseEnergyFractionThreshold")), mRThreshold (iConfig.getParameter<double>("RThreshold")), mTurnOnEverything (iConfig.getUntrackedParameter<std::string>("TurnOnEverything","")) { numberofevents = mEta = mEtaFineBin = mPhi = mPhiFineBin = mE = mE_80 = mE_3000 = mP = mP_80 = mP_3000 = mPt = mPt_80 = mPt_3000 = mMass = mMass_80 = mMass_3000 = mConstituents = mConstituents_80 = mEtaFirst = mPhiFirst = mEFirst = mEFirst_80 = mEFirst_3000 = mPtFirst = mPtFirst_80 = mPtFirst_3000 = mMjj = mMjj_3000 = mDelEta = mDelPhi = mDelPt = mMaxEInEmTowers = mMaxEInHadTowers = mHadEnergyInHO = mHadEnergyInHB = mHadEnergyInHF = mHadEnergyInHE = mHadEnergyInHO_80 = mHadEnergyInHB_80 = mHadEnergyInHE_80 = mHadEnergyInHO_3000 = mHadEnergyInHB_3000 = mHadEnergyInHE_3000 = mEmEnergyInEB = mEmEnergyInEE = mEmEnergyInHF = mEmEnergyInEB_80 = mEmEnergyInEE_80 = mEmEnergyInEB_3000 = mEmEnergyInEE_3000 = mEnergyFractionHadronic = mEnergyFractionEm = mHFLong = mHFTotal = mHFLong_80 = mHFLong_3000 = mHFShort = mHFShort_80 = mHFShort_3000 = mN90 = mCaloMEx = mCaloMEx_3000 = mCaloMEy = mCaloMEy_3000 = mCaloMETSig = mCaloMETSig_3000 = mCaloMET = mCaloMET_3000 = mCaloMETPhi = mCaloSumET = mCaloSumET_3000 = mHadTiming = mEmTiming = mNJetsEtaC = mNJetsEtaF = mNJets1 = mNJets2 = mAllGenJetsPt = mMatchedGenJetsPt = mAllGenJetsEta = mMatchedGenJetsEta = mGenJetMatchEnergyFraction = mReverseMatchEnergyFraction = mRMatch = mDeltaEta = mDeltaPhi = mEScale = mlinEScale = mDeltaE = mHadEnergyProfile = mEmEnergyProfile = mJetEnergyProfile = mHadJetEnergyProfile = mEMJetEnergyProfile = mEScale_pt10 = mEScaleFineBin = mpTScaleB_s = mpTScaleE_s = mpTScaleF_s = mpTScaleB_d = mpTScaleE_d = mpTScaleF_d = mpTScale_60_120_s = mpTScale_200_300_s = mpTScale_600_900_s = mpTScale_2700_3500_s = mpTScale_60_120_d = mpTScale_200_300_d = mpTScale_600_900_d = mpTScale_2700_3500_d = mpTScale1DB_60_120 = mpTScale1DE_60_120 = mpTScale1DF_60_120 = mpTScale1DB_200_300 = mpTScale1DE_200_300 = mpTScale1DF_200_300 = mpTScale1DB_600_900 = mpTScale1DE_600_900 = mpTScale1DF_600_900 = mpTScale1DB_2700_3500 = mpTScale1DE_2700_3500 = mpTScale1DF_2700_3500 = mpTScale1D_60_120 = mpTScale1D_200_300 = mpTScale1D_600_900 = mpTScale1D_2700_3500 = mHBEne = mHBTime = mHEEne = mHETime = mHFEne = mHFTime = mHOEne = mHOTime = mEBEne = mEBTime = mEEEne = mEETime = mPthat_80 = mPthat_3000 = 0; DQMStore* dbe = &*edm::Service<DQMStore>(); if (dbe) { dbe->setCurrentFolder("JetMET/RecoJetsV/CaloJetTask_" + mInputCollection.label()); // numberofevents = dbe->book1D("numberofevents","numberofevents", 3, 0 , 2); // mEta = dbe->book1D("Eta", "Eta", 100, -5, 5); mEtaFineBin = dbe->book1D("EtaFineBin_Pt10", "EtaFineBin_Pt10", 500, -5, 5); mEtaFineBin1p = dbe->book1D("EtaFineBin1p_Pt10", "EtaFineBin1p_Pt10", 100, 0, 1.3); mEtaFineBin2p = dbe->book1D("EtaFineBin2p_Pt10", "EtaFineBin2p_Pt10", 100, 1.3, 3); mEtaFineBin3p = dbe->book1D("EtaFineBin3p_Pt10", "EtaFineBin3p_Pt10", 100, 3, 5); mEtaFineBin1m = dbe->book1D("EtaFineBin1m_Pt10", "EtaFineBin1m_Pt10", 100, -1.3, 0); mEtaFineBin2m = dbe->book1D("EtaFineBin2m_Pt10", "EtaFineBin2m_Pt10", 100, -3, -1.3); mEtaFineBin3m = dbe->book1D("EtaFineBin3m_Pt10", "EtaFineBin3m_Pt10", 100, -5, -3); // mPhi = dbe->book1D("Phi", "Phi", 70, -3.5, 3.5); mPhiFineBin = dbe->book1D("PhiFineBin_Pt10", "PhiFineBin_Pt10", 350, -3.5, 3.5); // mE = dbe->book1D("E", "E", 100, 0, 500); mE_80 = dbe->book1D("E_80", "E_80", 100, 0, 4500); mE_3000 = dbe->book1D("E_3000", "E_3000", 100, 0, 6000); // mP = dbe->book1D("P", "P", 100, 0, 500); mP_80 = dbe->book1D("P_80", "P_80", 100, 0, 4500); mP_3000 = dbe->book1D("P_3000", "P_3000", 100, 0, 6000); // mPt = dbe->book1D("Pt", "Pt", 100, 0, 50); mPt_80 = dbe->book1D("Pt_80", "Pt_80", 100, 0, 140); mPt_3000 = dbe->book1D("Pt_3000", "Pt_3000", 100, 0, 4000); // mMass = dbe->book1D("Mass", "Mass", 100, 0, 25); mMass_80 = dbe->book1D("Mass_80", "Mass_80", 100, 0, 120); mMass_3000 = dbe->book1D("Mass_3000", "Mass_3000", 100, 0, 500); // mConstituents = dbe->book1D("Constituents", "# of Constituents", 100, 0, 100); mConstituents_80 = dbe->book1D("Constituents_80", "# of Constituents_80", 40, 0, 40); // 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); mEFirst_80 = dbe->book1D("EFirst_80", "EFirst_80", 100, 0, 180); mEFirst_3000 = dbe->book1D("EFirst_3000", "EFirst_3000", 100, 0, 4000); mPtFirst = dbe->book1D("PtFirst", "PtFirst", 100, 0, 50); mPtFirst_80 = dbe->book1D("PtFirst_80", "PtFirst_80", 100, 0, 140); mPtFirst_3000 = dbe->book1D("PtFirst_3000", "PtFirst_3000", 100, 0, 4000); // mMjj = dbe->book1D("Mjj", "Mjj", 100, 0, 2000); mMjj_3000 = dbe->book1D("Mjj_3000", "Mjj_3000", 100, 0, 10000); mDelEta = dbe->book1D("DelEta", "DelEta", 100, -.5, .5); mDelPhi = dbe->book1D("DelPhi", "DelPhi", 100, -.5, .5); mDelPt = dbe->book1D("DelPt", "DelPt", 100, -1, 1); // mMaxEInEmTowers = dbe->book1D("MaxEInEmTowers", "MaxEInEmTowers", 100, 0, 100); mMaxEInHadTowers = dbe->book1D("MaxEInHadTowers", "MaxEInHadTowers", 100, 0, 100); 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); // mHadEnergyInHO_80 = dbe->book1D("HadEnergyInHO_80", "HadEnergyInHO_80", 100, 0, 50); mHadEnergyInHB_80 = dbe->book1D("HadEnergyInHB_80", "HadEnergyInHB_80", 100, 0, 200); mHadEnergyInHE_80 = dbe->book1D("HadEnergyInHE_80", "HadEnergyInHE_80", 100, 0, 1000); mHadEnergyInHO_3000 = dbe->book1D("HadEnergyInHO_3000", "HadEnergyInHO_3000", 100, 0, 500); mHadEnergyInHB_3000 = dbe->book1D("HadEnergyInHB_3000", "HadEnergyInHB_3000", 100, 0, 3000); mHadEnergyInHE_3000 = dbe->book1D("HadEnergyInHE_3000", "HadEnergyInHE_3000", 100, 0, 2000); // mEmEnergyInEB = dbe->book1D("EmEnergyInEB", "EmEnergyInEB", 100, 0, 50); mEmEnergyInEE = dbe->book1D("EmEnergyInEE", "EmEnergyInEE", 100, 0, 50); mEmEnergyInHF = dbe->book1D("EmEnergyInHF", "EmEnergyInHF", 120, -20, 100); mEmEnergyInEB_80 = dbe->book1D("EmEnergyInEB_80", "EmEnergyInEB_80", 100, 0, 200); mEmEnergyInEE_80 = dbe->book1D("EmEnergyInEE_80", "EmEnergyInEE_80", 100, 0, 1000); mEmEnergyInEB_3000= dbe->book1D("EmEnergyInEB_3000", "EmEnergyInEB_3000", 100, 0, 3000); mEmEnergyInEE_3000= dbe->book1D("EmEnergyInEE_3000", "EmEnergyInEE_3000", 100, 0, 2000); mEnergyFractionHadronic = dbe->book1D("EnergyFractionHadronic", "EnergyFractionHadronic", 120, -0.1, 1.1); mEnergyFractionEm = dbe->book1D("EnergyFractionEm", "EnergyFractionEm", 120, -0.1, 1.1); // mHFTotal = dbe->book1D("HFTotal", "HFTotal", 100, 0, 500); mHFTotal_80 = dbe->book1D("HFTotal_80", "HFTotal_80", 100, 0, 3000); mHFTotal_3000 = dbe->book1D("HFTotal_3000", "HFTotal_3000", 100, 0, 6000); mHFLong = dbe->book1D("HFLong", "HFLong", 100, 0, 500); mHFLong_80 = dbe->book1D("HFLong_80", "HFLong_80", 100, 0, 200); mHFLong_3000 = dbe->book1D("HFLong_3000", "HFLong_3000", 100, 0, 1500); mHFShort = dbe->book1D("HFShort", "HFShort", 100, 0, 500); mHFShort_80 = dbe->book1D("HFShort_80", "HFShort_80", 100, 0, 200); mHFShort_3000 = dbe->book1D("HFShort_3000", "HFShort_3000", 100, 0, 1500); // mN90 = dbe->book1D("N90", "N90", 50, 0, 50); // mGenEta = dbe->book1D("GenEta", "GenEta", 100, -5, 5); mGenPhi = dbe->book1D("GenPhi", "GenPhi", 70, -3.5, 3.5); mGenPt = dbe->book1D("GenPt", "GenPt", 100, 0, 50); mGenPt_80 = dbe->book1D("GenPt_80", "GenPt_80", 100, 0, 140); mGenPt_3000 = dbe->book1D("GenPt_3000", "GenPt_3000", 100, 0, 1500); // mGenEtaFirst = dbe->book1D("GenEtaFirst", "GenEtaFirst", 100, -5, 5); mGenPhiFirst = dbe->book1D("GenPhiFirst", "GenPhiFirst", 70, -3.5, 3.5); // mCaloMEx = dbe->book1D("CaloMEx","CaloMEx",200,-150,150); mCaloMEx_3000 = dbe->book1D("CaloMEx_3000","CaloMEx_3000",100,-500,500); mCaloMEy = dbe->book1D("CaloMEy","CaloMEy",200,-150,150); mCaloMEy_3000 = dbe->book1D("CaloMEy_3000","CaloMEy_3000",100,-500,500); mCaloMETSig = dbe->book1D("CaloMETSig","CaloMETSig",100,0,15); mCaloMETSig_3000 = dbe->book1D("CaloMETSig_3000","CaloMETSig_3000",100,0,50); mCaloMET = dbe->book1D("CaloMET","CaloMET",100,0,150); mCaloMET_3000 = dbe->book1D("CaloMET_3000","CaloMET_3000",100,0,1000); mCaloMETPhi = dbe->book1D("CaloMETPhi","CaloMETPhi",70, -3.5, 3.5); mCaloSumET = dbe->book1D("CaloSumET","CaloSumET",100,0,500); mCaloSumET_3000 = dbe->book1D("CaloSumET_3000","CaloSumET_3000",100,3000,8000); // mHadTiming = dbe->book1D("HadTiming", "HadTiming", 75, -50, 100); mEmTiming = dbe->book1D("EMTiming", "EMTiming", 75, -50, 100); // mNJetsEtaC = dbe->book1D("NJetsEtaC_Pt10", "NJetsEtaC_Pt10", 15, 0, 15); mNJetsEtaF = dbe->book1D("NJetsEtaF_Pt10", "NJetsEtaF_Pt10", 15, 0, 15); // mNJets1 = dbe->bookProfile("NJets1", "NJets1", 100, 0, 200, 100, 0, 50, "s"); mNJets2 = dbe->bookProfile("NJets2", "NJets2", 100, 0, 4000, 100, 0, 50, "s"); // mHBEne = dbe->book1D( "HBEne", "HBEne", 1000, -20, 100 ); mHBTime = dbe->book1D( "HBTime", "HBTime", 200, -200, 200 ); mHEEne = dbe->book1D( "HEEne", "HEEne", 1000, -20, 100 ); mHETime = dbe->book1D( "HETime", "HETime", 200, -200, 200 ); mHOEne = dbe->book1D( "HOEne", "HOEne", 1000, -20, 100 ); mHOTime = dbe->book1D( "HOTime", "HOTime", 200, -200, 200 ); mHFEne = dbe->book1D( "HFEne", "HFEne", 1000, -20, 100 ); mHFTime = dbe->book1D( "HFTime", "HFTime", 200, -200, 200 ); mEBEne = dbe->book1D( "EBEne", "EBEne", 1000, -20, 100 ); mEBTime = dbe->book1D( "EBTime", "EBTime", 200, -200, 200 ); mEEEne = dbe->book1D( "EEEne", "EEEne", 1000, -20, 100 ); mEETime = dbe->book1D( "EETime", "EETime", 200, -200, 200 ); // mPthat_80 = dbe->book1D("Pthat_80", "Pthat_80", 100, 40.0, 160.0); mPthat_3000 = dbe->book1D("Pthat_3000", "Pthat_3000", 100, 2500.0, 4000.0); // double log10PtMin = 0.5; //=3.1622766 double log10PtMax = 3.75; //=5623.41325 int log10PtBins = 26; double etaMin = -5.; double etaMax = 5.; int etaBins = 50; double linPtMin = 5; double linPtMax = 155; int linPtBins = 15; int log10PtFineBins = 50; mAllGenJetsPt = dbe->book1D("GenJetLOGpT", "GenJet LOG(pT_gen)", log10PtBins, log10PtMin, log10PtMax); mMatchedGenJetsPt = dbe->book1D("MatchedGenJetLOGpT", "MatchedGenJet LOG(pT_gen)", log10PtBins, log10PtMin, log10PtMax); mAllGenJetsEta = dbe->book2D("GenJetEta", "GenJet Eta vs LOG(pT_gen)", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax); mMatchedGenJetsEta = dbe->book2D("MatchedGenJetEta", "MatchedGenJet Eta vs LOG(pT_gen)", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax); // if (mTurnOnEverything.compare("yes")==0) { mHadEnergyProfile = dbe->bookProfile2D("HadEnergyProfile", "HadEnergyProfile", 82, -41, 41, 73, 0, 73, 100, 0, 10000, "s"); mEmEnergyProfile = dbe->bookProfile2D("EmEnergyProfile", "EmEnergyProfile", 82, -41, 41, 73, 0, 73, 100, 0, 10000, "s"); } mJetEnergyProfile = dbe->bookProfile2D("JetEnergyProfile", "JetEnergyProfile", 50, -5, 5, 36, -3.1415987, 3.1415987, 100, 0, 10000, "s"); mHadJetEnergyProfile = dbe->bookProfile2D("HadJetEnergyProfile", "HadJetEnergyProfile", 50, -5, 5, 36, -3.1415987, 3.1415987, 100, 0, 10000, "s"); mEMJetEnergyProfile = dbe->bookProfile2D("EMJetEnergyProfile", "EMJetEnergyProfile", 50, -5, 5, 36, -3.1415987, 3.1415987, 100, 0, 10000, "s"); // if (mTurnOnEverything.compare("yes")==0) { mGenJetMatchEnergyFraction = dbe->book3D("GenJetMatchEnergyFraction", "GenJetMatchEnergyFraction vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 101, 0, 1.01); mReverseMatchEnergyFraction = dbe->book3D("ReverseMatchEnergyFraction", "ReverseMatchEnergyFraction vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 101, 0, 1.01); mRMatch = dbe->book3D("RMatch", "delta(R)(Gen-Calo) vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 60, 0, 3); mDeltaEta = dbe->book3D("DeltaEta", "DeltaEta vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, -1, 1); mDeltaPhi = dbe->book3D("DeltaPhi", "DeltaPhi vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, -1, 1); mEScale = dbe->book3D("EScale", "EnergyScale vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, 0, 2); mlinEScale = dbe->book3D("linEScale", "EnergyScale vs LOG(pT_gen) vs eta", linPtBins, linPtMin, linPtMax, etaBins, etaMin, etaMax, 100, 0, 2); mDeltaE = dbe->book3D("DeltaE", "DeltaE vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 2000, -200, 200); // mEScale_pt10 = dbe->book3D("EScale_pt10", "EnergyScale vs LOG(pT_gen) vs eta", log10PtBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, 0, 2); mEScaleFineBin = dbe->book3D("EScaleFineBins", "EnergyScale vs LOG(pT_gen) vs eta", log10PtFineBins, log10PtMin, log10PtMax, etaBins, etaMin, etaMax, 100, 0, 2); } mpTScaleB_s = dbe->bookProfile("pTScaleB_s", "pTScale_s_0<|eta|<1.3", log10PtBins, log10PtMin, log10PtMax, 0, 2, "s"); mpTScaleE_s = dbe->bookProfile("pTScaleE_s", "pTScale_s_1.3<|eta|<3.0", log10PtBins, log10PtMin, log10PtMax, 0, 2, "s"); mpTScaleF_s = dbe->bookProfile("pTScaleF_s", "pTScale_s_3.0<|eta|<5.0", log10PtBins, log10PtMin, log10PtMax, 0, 2, "s"); mpTScaleB_d = dbe->bookProfile("pTScaleB_d", "pTScale_d_0<|eta|<1.3", log10PtBins, log10PtMin, log10PtMax, 0, 2, " "); mpTScaleE_d = dbe->bookProfile("pTScaleE_d", "pTScale_d_1.3<|eta|<3.0", log10PtBins, log10PtMin, log10PtMax, 0, 2, " "); mpTScaleF_d = dbe->bookProfile("pTScaleF_d", "pTScale_d_3.0<|eta|<5.0", log10PtBins, log10PtMin, log10PtMax, 0, 2, " "); mpTScale_60_120_s = dbe->bookProfile("pTScale_60_120_s", "pTScale_s_60<pT<120", etaBins, etaMin, etaMax, 0., 2., "s"); mpTScale_200_300_s = dbe->bookProfile("pTScale_200_300_s", "pTScale_s_200<pT<300", etaBins, etaMin, etaMax, 0., 2., "s"); mpTScale_600_900_s = dbe->bookProfile("pTScale_600_900_s", "pTScale_s_600<pT<900", etaBins, etaMin, etaMax, 0., 2., "s"); mpTScale_2700_3500_s = dbe->bookProfile("pTScale_2700_3500_s", "pTScale_s_2700<pt<3500", etaBins, etaMin, etaMax, 0., 2., "s"); mpTScale_60_120_d = dbe->bookProfile("pTScale_60_120_d", "pTScale_d_60<pT<120", etaBins, etaMin, etaMax, 0., 2., " "); mpTScale_200_300_d = dbe->bookProfile("pTScale_200_300_d", "pTScale_d_200<pT<300", etaBins, etaMin, etaMax, 0., 2., " "); mpTScale_600_900_d = dbe->bookProfile("pTScale_600_900_d", "pTScale_d_600<pT<900", etaBins, etaMin, etaMax, 0., 2., " "); mpTScale_2700_3500_d = dbe->bookProfile("pTScale_2700_3500_d", "pTScale_d_2700<pt<3500", etaBins, etaMin, etaMax, 0., 2., " "); mpTScale1DB_60_120 = dbe->book1D("pTScale1DB_60_120", "pTScale_distribution_for_0<|eta|<1.3_60_120", 100, 0, 2); mpTScale1DE_60_120 = dbe->book1D("pTScale1DE_60_120", "pTScale_distribution_for_1.3<|eta|<3.0_60_120", 100, 0, 2); mpTScale1DF_60_120 = dbe->book1D("pTScale1DF_60_120", "pTScale_distribution_for_3.0<|eta|<5.0_60_120", 100, 0, 2); mpTScale1DB_200_300 = dbe->book1D("pTScale1DB_200_300", "pTScale_distribution_for_0<|eta|<1.3_200_300", 100, 0, 2); mpTScale1DE_200_300 = dbe->book1D("pTScale1DE_200_300", "pTScale_distribution_for_1.3<|eta|<3.0_200_300", 100, 0, 2); mpTScale1DF_200_300 = dbe->book1D("pTScale1DF_200_300", "pTScale_distribution_for_3.0<|eta|<5.0_200_300", 100, 0, 2); mpTScale1DB_600_900 = dbe->book1D("pTScale1DB_600_900", "pTScale_distribution_for_0<|eta|<1.3_600_900", 100, 0, 2); mpTScale1DE_600_900 = dbe->book1D("pTScale1DE_600_900", "pTScale_distribution_for_1.3<|eta|<3.0_600_900", 100, 0, 2); mpTScale1DF_600_900 = dbe->book1D("pTScale1DF_600_900", "pTScale_distribution_for_3.0<|eta|<5.0_600_900", 100, 0, 2); mpTScale1DB_2700_3500 = dbe->book1D("pTScale1DB_2700_3500", "pTScale_distribution_for_0<|eta|<1.3_2700_3500", 100, 0, 2); mpTScale1DE_2700_3500 = dbe->book1D("pTScale1DE_2700_3500", "pTScale_distribution_for_1.3<|eta|<3.0_2700_3500", 100, 0, 2); mpTScale1DF_2700_3500 = dbe->book1D("pTScale1DF_2700_3500", "pTScale_distribution_for_3.0<|eta|<5.0_2700_3500", 100, 0, 2); mpTScale1D_60_120 = dbe->book1D("pTScale1D_60_120", "pTScale_distribution_for_60<pT<120", 100, 0, 2); mpTScale1D_200_300 = dbe->book1D("pTScale1D_200_300", "pTScale_distribution_for_200<pT<300", 100, 0, 2); mpTScale1D_600_900 = dbe->book1D("pTScale1D_600_900", "pTScale_distribution_for_600<pT<900", 100, 0, 2); mpTScale1D_2700_3500 = dbe->book1D("pTScale1D_2700_3500", "pTScale_distribution_for_2700<pt<3500", 100, 0, 2); } if (mOutputFile.empty ()) { LogInfo("OutputInfo") << " CaloJet histograms will NOT be saved"; } else { LogInfo("OutputInfo") << " CaloJethistograms will be saved to file:" << mOutputFile; } }
| CaloJetTester::~CaloJetTester | ( | ) |
Definition at line 362 of file CaloJetTester.cc.
{
}
| void CaloJetTester::analyze | ( | const edm::Event & | mEvent, |
| const edm::EventSetup & | mSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 374 of file CaloJetTester.cc.
References ExpressReco_HICollisions_FallBack::caloTowers, deltaR(), reco::LeafCandidate::eta(), MonitorElement::Fill(), fillMatchHists(), edm::Event::getByLabel(), JetMatchingTools::getGenParticles(), edm::Event::getManyByType(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, i, edm::HandleBase::isValid(), j, metsig::jet, edm::InputTag::label(), mAllGenJetsEta, mAllGenJetsPt, mCaloMET, mCaloMET_3000, mCaloMETPhi, mCaloMETSig, mCaloMETSig_3000, mCaloMEx, mCaloMEx_3000, mCaloMEy, mCaloMEy_3000, mCaloSumET, mCaloSumET_3000, mConstituents, mConstituents_80, mE, mE_3000, mE_80, mEBEne, mEBTime, mEEEne, mEETime, mEFirst, mEFirst_3000, mEFirst_80, mEmEnergyInEB, mEmEnergyInEB_3000, mEmEnergyInEB_80, mEmEnergyInEE, mEmEnergyInEE_3000, mEmEnergyInEE_80, mEmEnergyInHF, mEmEnergyProfile, mEMJetEnergyProfile, mEmTiming, mEnergyFractionEm, mEnergyFractionHadronic, mEta, mEtaFineBin, mEtaFineBin1m, mEtaFineBin1p, mEtaFineBin2m, mEtaFineBin2p, mEtaFineBin3m, mEtaFineBin3p, mEtaFirst, reco::MET::mEtSig(), mGenEnergyFractionThreshold, mGenEta, mGenEtaFirst, mGenJetMatchEnergyFraction, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_3000, mGenPt_80, mHadEnergyInHB, mHadEnergyInHB_3000, mHadEnergyInHB_80, mHadEnergyInHE, mHadEnergyInHE_3000, mHadEnergyInHE_80, mHadEnergyInHF, mHadEnergyInHO, mHadEnergyInHO_3000, mHadEnergyInHO_80, mHadEnergyProfile, mHadJetEnergyProfile, mHadTiming, mHBEne, mHBTime, mHEEne, mHETime, mHFEne, mHFLong, mHFLong_3000, mHFLong_80, mHFShort, mHFShort_3000, mHFShort_80, mHFTime, mHFTotal, mHFTotal_3000, mHFTotal_80, mHOEne, mHOTime, mInputCollection, mInputGenCollection, mJetEnergyProfile, mMass, mMass_3000, mMass_80, mMatchGenPtThreshold, mMaxEInEmTowers, mMaxEInHadTowers, mMjj, mMjj_3000, mN90, mNJets1, mNJets2, mNJetsEtaC, mNJetsEtaF, mP, mP_3000, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_3000, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPthat_3000, mPthat_80, mReverseEnergyFractionThreshold, mReverseMatchEnergyFraction, mRMatch, mRThreshold, mTurnOnEverything, njet, numberofevents, JetMatchingTools::overlapEnergyFraction(), reco::LeafCandidate::phi(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::MET::sumEt().
{
double countsfornumberofevents = 1;
numberofevents->Fill(countsfornumberofevents);
// *********************************
// *** Get pThat
// *********************************
edm::Handle<HepMCProduct> evt;
mEvent.getByLabel("generator", evt);
if (evt.isValid()) {
HepMC::GenEvent * myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
double pthat = myGenEvent->event_scale();
mPthat_80->Fill(pthat);
mPthat_3000->Fill(pthat);
delete myGenEvent;
}
// ***********************************
// *** Get CaloMET
// ***********************************
const CaloMET *calomet;
edm::Handle<CaloMETCollection> calo;
mEvent.getByLabel("met", calo);
if (!calo.isValid()) {
edm::LogInfo("OutputInfo") << " failed to retrieve data required by MET Task";
edm::LogInfo("OutputInfo") << " MET Task cannot continue...!";
} else {
const CaloMETCollection *calometcol = calo.product();
calomet = &(calometcol->front());
double caloSumET = calomet->sumEt();
double caloMETSig = calomet->mEtSig();
double caloMET = calomet->pt();
double caloMEx = calomet->px();
double caloMEy = calomet->py();
double caloMETPhi = calomet->phi();
mCaloMEx->Fill(caloMEx);
mCaloMEx_3000->Fill(caloMEx);
mCaloMEy->Fill(caloMEy);
mCaloMEy_3000->Fill(caloMEy);
mCaloMET->Fill(caloMET);
mCaloMET_3000->Fill(caloMET);
mCaloMETPhi->Fill(caloMETPhi);
mCaloSumET->Fill(caloSumET);
mCaloSumET_3000->Fill(caloSumET);
mCaloMETSig->Fill(caloMETSig);
mCaloMETSig_3000->Fill(caloMETSig);
}
// ***********************************
// *** Get the CaloTower collection
// ***********************************
Handle<CaloTowerCollection> caloTowers;
mEvent.getByLabel( "towerMaker", caloTowers );
if (caloTowers.isValid()) {
for( CaloTowerCollection::const_iterator cal = caloTowers->begin(); cal != caloTowers->end(); ++ cal ){
//To compensate for the index
if (mTurnOnEverything.compare("yes")==0) {
if (cal->ieta() >> 0 ){mHadEnergyProfile->Fill (cal->ieta()-1, cal->iphi(), cal->hadEnergy());
mEmEnergyProfile->Fill (cal->ieta()-1, cal->iphi(), cal->emEnergy());}
mHadEnergyProfile->Fill (cal->ieta(), cal->iphi(), cal->hadEnergy());
mEmEnergyProfile->Fill (cal->ieta(), cal->iphi(), cal->emEnergy());
}
mHadTiming->Fill (cal->hcalTime());
mEmTiming->Fill (cal->ecalTime());
}
}
// ***********************************
// *** Get the RecHits collection
// ***********************************
try {
std::vector<edm::Handle<HBHERecHitCollection> > colls;
mEvent.getManyByType(colls);
std::vector<edm::Handle<HBHERecHitCollection> >::iterator i;
for (i=colls.begin(); i!=colls.end(); i++) {
for (HBHERecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
// std::cout << *j << std::endl;
if (j->id().subdet() == HcalBarrel) {
mHBEne->Fill(j->energy());
mHBTime->Fill(j->time());
}
if (j->id().subdet() == HcalEndcap) {
mHEEne->Fill(j->energy());
mHETime->Fill(j->time());
}
}
}
} catch (...) {
edm::LogInfo("OutputInfo") << " No HB/HE RecHits.";
}
try {
std::vector<edm::Handle<HFRecHitCollection> > colls;
mEvent.getManyByType(colls);
std::vector<edm::Handle<HFRecHitCollection> >::iterator i;
for (i=colls.begin(); i!=colls.end(); i++) {
for (HFRecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
// std::cout << *j << std::endl;
if (j->id().subdet() == HcalForward) {
mHFEne->Fill(j->energy());
mHFTime->Fill(j->time());
}
}
}
} catch (...) {
edm::LogInfo("OutputInfo") << " No HF RecHits.";
}
try {
std::vector<edm::Handle<HORecHitCollection> > colls;
mEvent.getManyByType(colls);
std::vector<edm::Handle<HORecHitCollection> >::iterator i;
for (i=colls.begin(); i!=colls.end(); i++) {
for (HORecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
if (j->id().subdet() == HcalOuter) {
mHOEne->Fill(j->energy());
mHOTime->Fill(j->time());
}
}
}
} catch (...) {
edm::LogInfo("OutputInfo") << " No HO RecHits.";
}
try {
std::vector<edm::Handle<EBRecHitCollection> > colls;
mEvent.getManyByType(colls);
std::vector<edm::Handle<EBRecHitCollection> >::iterator i;
for (i=colls.begin(); i!=colls.end(); i++) {
for (EBRecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
// if (j->id() == EcalBarrel) {
mEBEne->Fill(j->energy());
mEBTime->Fill(j->time());
// }
// std::cout << *j << std::endl;
// std::cout << j->id() << std::endl;
}
}
} catch (...) {
edm::LogInfo("OutputInfo") << " No EB RecHits.";
}
try {
std::vector<edm::Handle<EERecHitCollection> > colls;
mEvent.getManyByType(colls);
std::vector<edm::Handle<EERecHitCollection> >::iterator i;
for (i=colls.begin(); i!=colls.end(); i++) {
for (EERecHitCollection::const_iterator j=(*i)->begin(); j!=(*i)->end(); j++) {
// if (j->id().subdet() == EcalEndcap) {
mEEEne->Fill(j->energy());
mEETime->Fill(j->time());
// }
// std::cout << *j << std::endl;
}
}
} catch (...) {
edm::LogInfo("OutputInfo") << " No EE RecHits.";
}
//***********************************
//*** Get the Jet collection
//***********************************
math::XYZTLorentzVector p4tmp[2];
Handle<CaloJetCollection> caloJets;
mEvent.getByLabel(mInputCollection, caloJets);
if (!caloJets.isValid()) return;
CaloJetCollection::const_iterator jet = caloJets->begin ();
int jetIndex = 0;
int nJet = 0;
int nJetF = 0;
int nJetC = 0;
for (; jet != caloJets->end (); jet++, jetIndex++) {
if (mEta) mEta->Fill (jet->eta());
if (jet->pt() > 10.) {
if (fabs(jet->eta()) > 1.3)
nJetF++;
else
nJetC++;
}
if (jet->pt() > 10.) {
if (mEtaFineBin) mEtaFineBin->Fill (jet->eta());
if (mEtaFineBin1p) mEtaFineBin1p->Fill (jet->eta());
if (mEtaFineBin2p) mEtaFineBin2p->Fill (jet->eta());
if (mEtaFineBin3p) mEtaFineBin3p->Fill (jet->eta());
if (mEtaFineBin1m) mEtaFineBin1m->Fill (jet->eta());
if (mEtaFineBin2m) mEtaFineBin2m->Fill (jet->eta());
if (mEtaFineBin3m) mEtaFineBin3m->Fill (jet->eta());
if (mPhiFineBin) mPhiFineBin->Fill (jet->phi());
}
if (mPhi) mPhi->Fill (jet->phi());
if (mE) mE->Fill (jet->energy());
if (mE_80) mE_80->Fill (jet->energy());
if (mE_3000) mE_3000->Fill (jet->energy());
if (mP) mP->Fill (jet->p());
if (mP_80) mP_80->Fill (jet->p());
if (mP_3000) mP_3000->Fill (jet->p());
if (mPt) mPt->Fill (jet->pt());
if (mPt_80) mPt_80->Fill (jet->pt());
if (mPt_3000) mPt_3000->Fill (jet->pt());
if (mMass) mMass->Fill (jet->mass());
if (mMass_80) mMass_80->Fill (jet->mass());
if (mMass_3000) mMass_3000->Fill (jet->mass());
if (mConstituents) mConstituents->Fill (jet->nConstituents());
if (mConstituents_80) mConstituents_80->Fill (jet->nConstituents());
if (jet == caloJets->begin ()) { // first jet
if (mEtaFirst) mEtaFirst->Fill (jet->eta());
if (mPhiFirst) mPhiFirst->Fill (jet->phi());
if (mEFirst) mEFirst->Fill (jet->energy());
if (mEFirst_80) mEFirst_80->Fill (jet->energy());
if (mEFirst_3000) mEFirst_3000->Fill (jet->energy());
if (mPtFirst) mPtFirst->Fill (jet->pt());
if (mPtFirst_80) mPtFirst_80->Fill (jet->pt());
if (mPtFirst_3000) mPtFirst_3000->Fill (jet->pt());
}
if (jetIndex == 0) {
nJet++;
p4tmp[0] = jet->p4();
}
if (jetIndex == 1) {
nJet++;
p4tmp[1] = jet->p4();
}
if (mMaxEInEmTowers) mMaxEInEmTowers->Fill (jet->maxEInEmTowers());
if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->maxEInHadTowers());
if (mHadEnergyInHO) mHadEnergyInHO->Fill (jet->hadEnergyInHO());
if (mHadEnergyInHO_80) mHadEnergyInHO_80->Fill (jet->hadEnergyInHO());
if (mHadEnergyInHO_3000) mHadEnergyInHO_3000->Fill (jet->hadEnergyInHO());
if (mHadEnergyInHB) mHadEnergyInHB->Fill (jet->hadEnergyInHB());
if (mHadEnergyInHB_80) mHadEnergyInHB_80->Fill (jet->hadEnergyInHB());
if (mHadEnergyInHB_3000) mHadEnergyInHB_3000->Fill (jet->hadEnergyInHB());
if (mHadEnergyInHF) mHadEnergyInHF->Fill (jet->hadEnergyInHF());
if (mHadEnergyInHE) mHadEnergyInHE->Fill (jet->hadEnergyInHE());
if (mHadEnergyInHE_80) mHadEnergyInHE_80->Fill (jet->hadEnergyInHE());
if (mHadEnergyInHE_3000) mHadEnergyInHE_3000->Fill (jet->hadEnergyInHE());
if (mEmEnergyInEB) mEmEnergyInEB->Fill (jet->emEnergyInEB());
if (mEmEnergyInEB_80) mEmEnergyInEB_80->Fill (jet->emEnergyInEB());
if (mEmEnergyInEB_3000) mEmEnergyInEB_3000->Fill (jet->emEnergyInEB());
if (mEmEnergyInEE) mEmEnergyInEE->Fill (jet->emEnergyInEE());
if (mEmEnergyInEE_80) mEmEnergyInEE_80->Fill (jet->emEnergyInEE());
if (mEmEnergyInEE_3000) mEmEnergyInEE_3000->Fill (jet->emEnergyInEE());
if (mEmEnergyInHF) mEmEnergyInHF->Fill (jet->emEnergyInHF());
if (mEnergyFractionHadronic) mEnergyFractionHadronic->Fill (jet->energyFractionHadronic());
if (mEnergyFractionEm) mEnergyFractionEm->Fill (jet->emEnergyFraction());
if (mHFTotal) mHFTotal->Fill (jet->hadEnergyInHF()+jet->emEnergyInHF());
if (mHFTotal_80) mHFTotal_80->Fill (jet->hadEnergyInHF()+jet->emEnergyInHF());
if (mHFTotal_3000) mHFTotal_3000->Fill (jet->hadEnergyInHF()+jet->emEnergyInHF());
if (mHFLong) mHFLong->Fill (jet->hadEnergyInHF()*0.5+jet->emEnergyInHF());
if (mHFLong_80) mHFLong_80->Fill (jet->hadEnergyInHF()*0.5+jet->emEnergyInHF());
if (mHFLong_3000) mHFLong_3000->Fill (jet->hadEnergyInHF()*0.5+jet->emEnergyInHF());
if (mHFShort) mHFShort->Fill (jet->hadEnergyInHF()*0.5);
if (mHFShort_80) mHFShort_80->Fill (jet->hadEnergyInHF()*0.5);
if (mHFShort_3000) mHFShort_3000->Fill (jet->hadEnergyInHF()*0.5);
if (mN90) mN90->Fill (jet->n90());
mJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->energy());
mHadJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->hadEnergyInHO()+jet->hadEnergyInHB()+jet->hadEnergyInHF()+jet->hadEnergyInHE());
mEMJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->emEnergyInEB()+jet->emEnergyInEE()+jet->emEnergyInHF());
}
if (mNJetsEtaC) mNJetsEtaC->Fill( nJetC );
if (mNJetsEtaF) mNJetsEtaF->Fill( nJetF );
if (nJet == 2) {
if (mMjj) mMjj->Fill( (p4tmp[0]+p4tmp[1]).mass() );
if (mMjj_3000) mMjj_3000->Fill( (p4tmp[0]+p4tmp[1]).mass() );
}
// Count Jets above Pt cut
for (int istep = 0; istep < 100; ++istep) {
int njet = 0;
float ptStep = (istep * (200./100.));
for ( CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
if ( cal->pt() > ptStep ) njet++;
}
mNJets1->Fill( ptStep, njet );
}
for (int istep = 0; istep < 100; ++istep) {
int njet = 0;
float ptStep = (istep * (4000./100.));
for ( CaloJetCollection::const_iterator cal = caloJets->begin(); cal != caloJets->end(); ++ cal ) {
if ( cal->pt() > ptStep ) njet++;
}
mNJets2->Fill( ptStep, njet );
}
// Gen jet analysis
Handle<GenJetCollection> genJets;
mEvent.getByLabel(mInputGenCollection, genJets);
if (!genJets.isValid()) return;
GenJetCollection::const_iterator gjet = genJets->begin ();
int gjetIndex = 0;
for (; gjet != genJets->end (); gjet++, gjetIndex++) {
if (mGenEta) mGenEta->Fill (gjet->eta());
if (mGenPhi) mGenPhi->Fill (gjet->phi());
if (mGenPt) mGenPt->Fill (gjet->pt());
if (mGenPt_80) mGenPt_80->Fill (gjet->pt());
if (mGenPt_3000) mGenPt_3000->Fill (gjet->pt());
if (gjet == genJets->begin ()) { // first jet
if (mGenEtaFirst) mGenEtaFirst->Fill (gjet->eta());
if (mGenPhiFirst) mGenPhiFirst->Fill (gjet->phi());
}
}
// now match CaloJets to GenJets
JetMatchingTools jetMatching (mEvent);
if (!(mInputGenCollection.label().empty())) {
// Handle<GenJetCollection> genJets;
// mEvent.getByLabel(mInputGenCollection, genJets);
std::vector <std::vector <const reco::GenParticle*> > genJetConstituents (genJets->size());
std::vector <std::vector <const reco::GenParticle*> > caloJetConstituents (caloJets->size());
if (mRThreshold > 0) {
}
else {
for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
genJetConstituents [iGenJet] = jetMatching.getGenParticles ((*genJets) [iGenJet]);
}
for (unsigned iCaloJet = 0; iCaloJet < caloJets->size(); ++iCaloJet) {
caloJetConstituents [iCaloJet] = jetMatching.getGenParticles ((*caloJets) [iCaloJet], false);
}
}
for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) { //****************************************************************
//for (unsigned iGenJet = 0; iGenJet < 1; ++iGenJet) { // only FIRST Jet !!!!
const GenJet& genJet = (*genJets) [iGenJet];
double genJetPt = genJet.pt();
//std::cout << iGenJet <<". Genjet: pT = " << genJetPt << "GeV" << std::endl; // *****************************************************
if (fabs(genJet.eta()) > 5.) continue; // out of detector
if (genJetPt < mMatchGenPtThreshold) continue; // no low momentum
double logPtGen = log10 (genJetPt);
mAllGenJetsPt->Fill (logPtGen);
mAllGenJetsEta->Fill (logPtGen, genJet.eta());
if (caloJets->size() <= 0) continue; // no CaloJets - nothing to match
if (mRThreshold > 0) {
unsigned iCaloJetBest = 0;
double deltaRBest = 999.;
for (unsigned iCaloJet = 0; iCaloJet < caloJets->size(); ++iCaloJet) {
double dR = deltaR (genJet.eta(), genJet.phi(), (*caloJets) [iCaloJet].eta(), (*caloJets) [iCaloJet].phi());
if (deltaRBest < mRThreshold && dR < mRThreshold && genJet.pt() > 5.) {
/*
std::cout << "Yet another matched jet for GenJet pt=" << genJet.pt()
<< " previous CaloJet pt/dr: " << (*caloJets) [iCaloJetBest].pt() << '/' << deltaRBest
<< " new CaloJet pt/dr: " << (*caloJets) [iCaloJet].pt() << '/' << dR
<< std::endl;
*/
}
if (dR < deltaRBest) {
iCaloJetBest = iCaloJet;
deltaRBest = dR;
}
}
if (mTurnOnEverything.compare("yes")==0) {
mRMatch->Fill (logPtGen, genJet.eta(), deltaRBest);
}
if (deltaRBest < mRThreshold) { // Matched
fillMatchHists (genJet, (*caloJets) [iCaloJetBest]);
}
}
else {
unsigned iCaloJetBest = 0;
double energyFractionBest = 0.;
for (unsigned iCaloJet = 0; iCaloJet < caloJets->size(); ++iCaloJet) {
double energyFraction = jetMatching.overlapEnergyFraction (genJetConstituents [iGenJet],
caloJetConstituents [iCaloJet]);
if (energyFraction > energyFractionBest) {
iCaloJetBest = iCaloJet;
energyFractionBest = energyFraction;
}
}
if (mTurnOnEverything.compare("yes")==0) {
mGenJetMatchEnergyFraction->Fill (logPtGen, genJet.eta(), energyFractionBest);
}
if (energyFractionBest > mGenEnergyFractionThreshold) { // good enough
double reverseEnergyFraction = jetMatching.overlapEnergyFraction (caloJetConstituents [iCaloJetBest],
genJetConstituents [iGenJet]);
if (mTurnOnEverything.compare("yes")==0) {
mReverseMatchEnergyFraction->Fill (logPtGen, genJet.eta(), reverseEnergyFraction);
}
if (reverseEnergyFraction > mReverseEnergyFractionThreshold) { // Matched
fillMatchHists (genJet, (*caloJets) [iCaloJetBest]);
}
}
}
}
}
}
| void CaloJetTester::beginJob | ( | void | ) | [virtual] |
| void CaloJetTester::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 369 of file CaloJetTester.cc.
References mOutputFile.
{
if (!mOutputFile.empty() && &*edm::Service<DQMStore>()) edm::Service<DQMStore>()->save (mOutputFile);
}
| void CaloJetTester::fillMatchHists | ( | const reco::GenJet & | fGenJet, |
| const reco::CaloJet & | fCaloJet | ||
| ) | [private] |
Definition at line 785 of file CaloJetTester.cc.
References reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), MonitorElement::Fill(), mDelEta, mDelPhi, mDelPt, mDeltaE, mDeltaEta, mDeltaPhi, mEScale, mEScale_pt10, mEScaleFineBin, mlinEScale, mMatchedGenJetsEta, mMatchedGenJetsPt, mpTScale1D_200_300, mpTScale1D_2700_3500, mpTScale1D_600_900, mpTScale1D_60_120, mpTScale1DB_200_300, mpTScale1DB_2700_3500, mpTScale1DB_600_900, mpTScale1DB_60_120, mpTScale1DE_200_300, mpTScale1DE_2700_3500, mpTScale1DE_600_900, mpTScale1DE_60_120, mpTScale1DF_200_300, mpTScale1DF_2700_3500, mpTScale1DF_600_900, mpTScale1DF_60_120, mpTScale_200_300_d, mpTScale_200_300_s, mpTScale_2700_3500_d, mpTScale_2700_3500_s, mpTScale_600_900_d, mpTScale_600_900_s, mpTScale_60_120_d, mpTScale_60_120_s, mpTScaleB_d, mpTScaleB_s, mpTScaleE_d, mpTScaleE_s, mpTScaleF_d, mpTScaleF_s, mTurnOnEverything, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().
Referenced by analyze().
{
double logPtGen = log10 (fGenJet.pt());
double PtGen = fGenJet.pt();
double PtCalo = fCaloJet.pt();
mMatchedGenJetsPt->Fill (logPtGen);
mMatchedGenJetsEta->Fill (logPtGen, fGenJet.eta());
double PtThreshold = 10.;
if (mTurnOnEverything.compare("yes")==0) {
mDeltaEta->Fill (logPtGen, fGenJet.eta(), fCaloJet.eta()-fGenJet.eta());
mDeltaPhi->Fill (logPtGen, fGenJet.eta(), fCaloJet.phi()-fGenJet.phi());
mEScale->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
mlinEScale->Fill (fGenJet.pt(), fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
mDeltaE->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()-fGenJet.energy());
mEScaleFineBin->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
if (fGenJet.pt()>PtThreshold) {
mEScale_pt10->Fill (logPtGen, fGenJet.eta(), fCaloJet.energy()/fGenJet.energy());
}
}
if (fCaloJet.pt() > PtThreshold) {
mDelEta->Fill (fGenJet.eta()-fCaloJet.eta());
mDelPhi->Fill (fGenJet.phi()-fCaloJet.phi());
mDelPt->Fill ((fGenJet.pt()-fCaloJet.pt())/fGenJet.pt());
}
if (fabs(fGenJet.eta())<1.3) {
mpTScaleB_s->Fill (log10(PtGen), PtCalo/PtGen);
mpTScaleB_d->Fill (log10(PtGen), PtCalo/PtGen);
if (PtGen>60.0 && PtGen<120.0) {
mpTScale1DB_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>200.0 && PtGen<300.0) {
mpTScale1DB_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>600.0 && PtGen<900.0) {
mpTScale1DB_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>2700.0 && PtGen<3500.0) {
mpTScale1DB_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
}
}
if (fabs(fGenJet.eta())>1.3 && fabs(fGenJet.eta())<3.0) {
mpTScaleE_s->Fill (log10(PtGen), PtCalo/PtGen);
mpTScaleE_d->Fill (log10(PtGen), PtCalo/PtGen);
if (PtGen>60.0 && PtGen<120.0) {
mpTScale1DE_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>200.0 && PtGen<300.0) {
mpTScale1DE_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>600.0 && PtGen<900.0) {
mpTScale1DE_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>2700.0 && PtGen<3500.0) {
mpTScale1DE_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
}
}
if (fabs(fGenJet.eta())>3.0 && fabs(fGenJet.eta())<5.0) {
mpTScaleF_s->Fill (log10(PtGen), PtCalo/PtGen);
mpTScaleF_d->Fill (log10(PtGen), PtCalo/PtGen);
if (PtGen>60.0 && PtGen<120.0) {
mpTScale1DF_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>200.0 && PtGen<300.0) {
mpTScale1DF_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>600.0 && PtGen<900.0) {
mpTScale1DF_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (PtGen>2700.0 && PtGen<3500.0) {
mpTScale1DF_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
}
}
if (fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
mpTScale_60_120_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale_60_120_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale1D_60_120->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
mpTScale_200_300_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale_200_300_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale1D_200_300->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
mpTScale_600_900_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale_600_900_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale1D_600_900->Fill (fCaloJet.pt()/fGenJet.pt());
}
if (fGenJet.pt()>2700.0 && fGenJet.pt()<3500.0) {
mpTScale_2700_3500_s->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale_2700_3500_d->Fill (fGenJet.eta(),fCaloJet.pt()/fGenJet.pt());
mpTScale1D_2700_3500->Fill (fCaloJet.pt()/fGenJet.pt());
}
}
std::string CaloJetTester::inputCaloMETLabel_ [private] |
Definition at line 41 of file CaloJetTester.h.
std::string CaloJetTester::inputGenMETLabel_ [private] |
Definition at line 40 of file CaloJetTester.h.
edm::InputTag CaloJetTester::inputMETLabel_ [private] |
Definition at line 38 of file CaloJetTester.h.
MonitorElement* CaloJetTester::mAllGenJetsEta [private] |
Definition at line 146 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mAllGenJetsPt [private] |
Definition at line 144 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMET [private] |
Definition at line 219 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMET_3000 [private] |
Definition at line 220 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMETPhi [private] |
Definition at line 221 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMETSig [private] |
Definition at line 224 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMETSig_3000 [private] |
Definition at line 225 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMEx [private] |
Definition at line 215 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMEx_3000 [private] |
Definition at line 216 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMEy [private] |
Definition at line 217 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloMEy_3000 [private] |
Definition at line 218 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloSumET [private] |
Definition at line 222 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mCaloSumET_3000 [private] |
Definition at line 223 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mConstituents [private] |
Definition at line 69 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mConstituents_3000 [private] |
Definition at line 71 of file CaloJetTester.h.
MonitorElement* CaloJetTester::mConstituents_80 [private] |
Definition at line 70 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mDelEta [private] |
Definition at line 194 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mDelPhi [private] |
Definition at line 195 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mDelPt [private] |
Definition at line 196 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mDeltaE [private] |
Definition at line 155 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mDeltaEta [private] |
Definition at line 151 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mDeltaPhi [private] |
Definition at line 152 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mE [private] |
Definition at line 57 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mE_3000 [private] |
Definition at line 59 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mE_80 [private] |
Definition at line 58 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEBEne [private] |
Definition at line 236 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEBTime [private] |
Definition at line 237 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEEEne [private] |
Definition at line 238 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEETime [private] |
Definition at line 239 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEFirst [private] |
Definition at line 78 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEFirst_3000 [private] |
Definition at line 80 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEFirst_80 [private] |
Definition at line 79 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInEB [private] |
Definition at line 108 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInEB_3000 [private] |
Definition at line 113 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInEB_80 [private] |
Definition at line 111 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInEE [private] |
Definition at line 109 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInEE_3000 [private] |
Definition at line 114 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInEE_80 [private] |
Definition at line 112 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyInHF [private] |
Definition at line 110 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmEnergyProfile [private] |
Definition at line 209 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEMJetEnergyProfile [private] |
Definition at line 212 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEmTiming [private] |
Definition at line 73 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEnergyFractionEm [private] |
Definition at line 116 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
Definition at line 115 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEScale [private] |
Definition at line 153 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mEScale_pt10 [private] |
Definition at line 157 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mEScaleFineBin [private] |
Definition at line 158 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mEta [private] |
Definition at line 47 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin [private] |
Definition at line 48 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin1m [private] |
Definition at line 52 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin1p [private] |
Definition at line 49 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin2m [private] |
Definition at line 53 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin2p [private] |
Definition at line 50 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin3m [private] |
Definition at line 54 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFineBin3p [private] |
Definition at line 51 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mEtaFirst [private] |
Definition at line 76 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
std::string CaloJetTester::METType_ [private] |
Definition at line 39 of file CaloJetTester.h.
double CaloJetTester::mGenEnergyFractionThreshold [private] |
Definition at line 200 of file CaloJetTester.h.
Referenced by analyze().
MonitorElement* CaloJetTester::mGenEta [private] |
Definition at line 133 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mGenEtaFirst [private] |
Definition at line 140 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
Definition at line 148 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mGenPhi [private] |
Definition at line 134 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mGenPhiFirst [private] |
Definition at line 141 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mGenPt [private] |
Definition at line 135 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mGenPt_3000 [private] |
Definition at line 137 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mGenPt_80 [private] |
Definition at line 136 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHB [private] |
Definition at line 99 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHB_3000 [private] |
Definition at line 106 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHB_80 [private] |
Definition at line 103 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHE [private] |
Definition at line 101 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHE_3000 [private] |
Definition at line 107 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHE_80 [private] |
Definition at line 104 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHF [private] |
Definition at line 100 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHO [private] |
Definition at line 98 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHO_3000 [private] |
Definition at line 105 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyInHO_80 [private] |
Definition at line 102 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadEnergyProfile [private] |
Definition at line 208 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
Definition at line 211 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHadTiming [private] |
Definition at line 72 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHBEne [private] |
Definition at line 228 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHBTime [private] |
Definition at line 229 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHEEne [private] |
Definition at line 230 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHETime [private] |
Definition at line 231 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFEne [private] |
Definition at line 234 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFLong [private] |
Definition at line 120 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFLong_3000 [private] |
Definition at line 122 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFLong_80 [private] |
Definition at line 121 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFShort [private] |
Definition at line 123 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFShort_3000 [private] |
Definition at line 125 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFShort_80 [private] |
Definition at line 124 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFTime [private] |
Definition at line 235 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFTotal [private] |
Definition at line 117 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFTotal_3000 [private] |
Definition at line 119 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHFTotal_80 [private] |
Definition at line 118 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHOEne [private] |
Definition at line 232 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mHOTime [private] |
Definition at line 233 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
edm::InputTag CaloJetTester::mInputCollection [private] |
Definition at line 35 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
Definition at line 36 of file CaloJetTester.h.
Referenced by analyze().
MonitorElement* CaloJetTester::mJetEnergyProfile [private] |
Definition at line 210 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mlinEScale [private] |
Definition at line 154 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mMass [private] |
Definition at line 66 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mMass_3000 [private] |
Definition at line 68 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mMass_80 [private] |
Definition at line 67 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mMatchedGenJetsEta [private] |
Definition at line 147 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mMatchedGenJetsPt [private] |
Definition at line 145 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
double CaloJetTester::mMatchGenPtThreshold [private] |
Definition at line 199 of file CaloJetTester.h.
Referenced by analyze().
MonitorElement* CaloJetTester::mMaxEInEmTowers [private] |
Definition at line 96 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mMaxEInHadTowers [private] |
Definition at line 97 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mMjj [private] |
Definition at line 92 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mMjj_3000 [private] |
Definition at line 93 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mN90 [private] |
Definition at line 126 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mNJets1 [private] |
Definition at line 88 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mNJets2 [private] |
Definition at line 89 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mNJetsEtaC [private] |
Definition at line 85 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mNJetsEtaF [private] |
Definition at line 86 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
std::string CaloJetTester::mOutputFile [private] |
Definition at line 37 of file CaloJetTester.h.
Referenced by CaloJetTester(), and endJob().
MonitorElement* CaloJetTester::mP [private] |
Definition at line 60 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mP_3000 [private] |
Definition at line 62 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mP_80 [private] |
Definition at line 61 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPhi [private] |
Definition at line 55 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPhiFineBin [private] |
Definition at line 56 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPhiFirst [private] |
Definition at line 77 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPt [private] |
Definition at line 63 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPt_3000 [private] |
Definition at line 65 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPt_80 [private] |
Definition at line 64 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPtFirst [private] |
Definition at line 81 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPtFirst_3000 [private] |
Definition at line 83 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPtFirst_80 [private] |
Definition at line 82 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPthat_3000 [private] |
Definition at line 130 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mPthat_80 [private] |
Definition at line 129 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mpTScale1D_200_300 [private] |
Definition at line 190 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
Definition at line 192 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1D_600_900 [private] |
Definition at line 191 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1D_60_120 [private] |
Definition at line 189 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DB_200_300 [private] |
Definition at line 180 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
Definition at line 186 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DB_600_900 [private] |
Definition at line 183 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DB_60_120 [private] |
Definition at line 177 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DE_200_300 [private] |
Definition at line 181 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
Definition at line 187 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DE_600_900 [private] |
Definition at line 184 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DE_60_120 [private] |
Definition at line 178 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DF_200_300 [private] |
Definition at line 182 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
Definition at line 188 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DF_600_900 [private] |
Definition at line 185 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale1DF_60_120 [private] |
Definition at line 179 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale_200_300_d [private] |
Definition at line 173 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale_200_300_s [private] |
Definition at line 168 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
Definition at line 175 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
Definition at line 170 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale_600_900_d [private] |
Definition at line 174 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale_600_900_s [private] |
Definition at line 169 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale_60_120_d [private] |
Definition at line 172 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScale_60_120_s [private] |
Definition at line 167 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScaleB_d [private] |
Definition at line 163 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScaleB_s [private] |
Definition at line 160 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScaleE_d [private] |
Definition at line 164 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScaleE_s [private] |
Definition at line 161 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScaleF_d [private] |
Definition at line 165 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::mpTScaleF_s [private] |
Definition at line 162 of file CaloJetTester.h.
Referenced by CaloJetTester(), and fillMatchHists().
double CaloJetTester::mReverseEnergyFractionThreshold [private] |
Definition at line 201 of file CaloJetTester.h.
Referenced by analyze().
Definition at line 149 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
MonitorElement* CaloJetTester::mRMatch [private] |
Definition at line 150 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
double CaloJetTester::mRThreshold [private] |
Definition at line 202 of file CaloJetTester.h.
Referenced by analyze().
std::string CaloJetTester::mTurnOnEverything [private] |
Definition at line 205 of file CaloJetTester.h.
Referenced by analyze(), CaloJetTester(), and fillMatchHists().
MonitorElement* CaloJetTester::numberofevents [private] |
Definition at line 44 of file CaloJetTester.h.
Referenced by analyze(), and CaloJetTester().
1.7.3