|
|
|
#include <JPTJetTester.h>
Definition at line 21 of file JPTJetTester.h.
| JPTJetTester::JPTJetTester | ( | const edm::ParameterSet & | iConfig | ) |
Definition at line 49 of file JPTJetTester.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::book3D(), DQMStore::bookProfile(), DQMStore::bookProfile2D(), jptDQMConfig_cff::etaBins, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, edm::InputTag::label(), mAllGenJetsEta, mAllGenJetsPt, mCaloMET, mCaloMET_3000, mCaloMETPhi, mCaloMETSig, mCaloMETSig_3000, mCaloMEx, mCaloMEx_3000, mCaloMEy, mCaloMEy_3000, mCaloSumET, mCaloSumET_3000, mDelEta, mDelPhi, mDelPt, mDeltaE, mDeltaEta, mDeltaPhi, mE, mE_3000, mE_80, mEBEne, mEBTime, mEEEne, mEETime, mEFirst, mEFirst_3000, mEFirst_80, mEmEnergyProfile, mEmTiming, mEScale, mEScale_pt10, mEScaleFineBin, mEta, mEtaFineBin, mEtaFineBin1m, mEtaFineBin1p, mEtaFineBin2m, mEtaFineBin2p, mEtaFineBin3m, mEtaFineBin3p, mEtaFirst, mGenEta, mGenEtaFirst, mGenJetMatchEnergyFraction, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_3000, mGenPt_80, mHadEnergyProfile, mHadTiming, mHBEne, mHBTime, mHEEne, mHETime, mHFEne, mHFTime, mHOEne, mHOTime, mInputCollection, mJetEnergyProfile, mlinEScale, mMass, mMass_3000, mMass_80, mMatchedGenJetsEta, mMatchedGenJetsPt, mMjj, mMjj_3000, mMultiplicityElecInVertexInCalo, mMultiplicityElecInVertexOutCalo, mMultiplicityElecOutVertexInCalo, mMultiplicityMuonInVertexInCalo, mMultiplicityMuonInVertexOutCalo, mMultiplicityMuonOutVertexInCalo, mMultiplicityPionInVertexInCalo, mMultiplicityPionInVertexOutCalo, mMultiplicityPionOutVertexInCalo, mNJets1, mNJets2, mNJetsEtaC, mNJetsEtaF, mOutputFile, mP, mP_3000, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_3000, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPtFractionElecInVertexInCalo, mPtFractionElecInVertexOutCalo, mPtFractionElecOutVertexInCalo, mPtFractionMuonInVertexInCalo, mPtFractionMuonInVertexOutCalo, mPtFractionMuonOutVertexInCalo, mPtFractionPionInVertexInCalo, mPtFractionPionInVertexOutCalo, mPtFractionPionOutVertexInCalo, 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 */ = mPtFractionPionInVertexInCalo = mPtFractionPionInVertexOutCalo = mPtFractionPionOutVertexInCalo = mPtFractionMuonInVertexInCalo = mPtFractionMuonInVertexOutCalo = mPtFractionMuonOutVertexInCalo = mPtFractionElecInVertexInCalo = mPtFractionElecInVertexOutCalo = mPtFractionElecOutVertexInCalo = mMultiplicityPionInVertexInCalo = mMultiplicityPionInVertexOutCalo = mMultiplicityPionOutVertexInCalo = mMultiplicityMuonInVertexInCalo = mMultiplicityMuonInVertexOutCalo = mMultiplicityMuonOutVertexInCalo = mMultiplicityElecInVertexInCalo = mMultiplicityElecInVertexOutCalo = mMultiplicityElecOutVertexInCalo // = 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/JPTJetTask_" + 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); // mPtFractionPionInVertexInCalo = dbe->book1D("PtFractionPionInVertexInCalo", "PtFractionPionInVertexInCalo",51,0.,1.02); mPtFractionPionInVertexOutCalo = dbe->book1D("PtFractionPionInVertexOutCalo","PtFractionPionInVertexOutCalo",51,0.,1.02); mPtFractionPionOutVertexInCalo = dbe->book1D("PtFractionPionOutVertexInCalo","PtFractionPionOutVertexInCalo",51,0.,1.02); mPtFractionMuonInVertexInCalo = dbe->book1D("PtFractionMuonInVertexInCalo", "PtFractionMuonInVertexInCalo",51,0.,1.02); mPtFractionMuonInVertexOutCalo = dbe->book1D("PtFractionMuonInVertexOutCalo","PtFractionMuonInVertexOutCalo",51,0.,1.02); mPtFractionMuonOutVertexInCalo = dbe->book1D("PtFractionMuonOutVertexInCalo","PtFractionMuonOutVertexInCalo",51,0.,1.02); mPtFractionElecInVertexInCalo = dbe->book1D("PtFractionElecInVertexInCalo", "PtFractionElecInVertexInCalo",51,0.,1.02); mPtFractionElecInVertexOutCalo = dbe->book1D("PtFractionElecInVertexOutCalo","PtFractionElecInVertexOutCalo",51,0.,1.02); mPtFractionElecOutVertexInCalo = dbe->book1D("PtFractionElecOutVertexInCalo","PtFractionElecOutVertexInCalo",51,0.,1.02); // mMultiplicityPionInVertexInCalo = dbe->book1D("MultiplicityPionInVertexInCalo","MultiplicityPionInVertexInCalo",100,0.,100.); mMultiplicityPionInVertexOutCalo = dbe->book1D("MultiplicityPionInVertexOutCalo","MultiplicityPionInVertexOutCalo",100,0.,100.); mMultiplicityPionOutVertexInCalo = dbe->book1D("MultiplicityPionOutVertexInCalo","MultiplicityPionOutVertexInCalo",100,0.,100.); mMultiplicityMuonInVertexInCalo = dbe->book1D("MultiplicityMuonInVertexInCalo","MultiplicityMuonInVertexInCalo",100,0.,100.); mMultiplicityMuonInVertexOutCalo = dbe->book1D("MultiplicityMuonInVertexOutCalo","MultiplicityMuonInVertexOutCalo",100,0.,100.); mMultiplicityMuonOutVertexInCalo = dbe->book1D("MultiplicityMuonOutVertexInCalo","MultiplicityMuonOutVertexInCalo",100,0.,100.); mMultiplicityElecInVertexInCalo = dbe->book1D("MultiplicityElecInVertexInCalo","MultiplicityElecInVertexInCalo",100,0.,100.); mMultiplicityElecInVertexOutCalo = dbe->book1D("MultiplicityElecInVertexOutCalo","MultiplicityElecInVertexOutCalo",100,0.,100.); mMultiplicityElecOutVertexInCalo = dbe->book1D("MultiplicityElecOutVertexInCalo","MultiplicityElecOutVertexInCalo",100,0.,100.); /* 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 (dbe) if (mOutputFile.empty ()) { LogInfo("OutputInfo") << " JPTJet histograms will NOT be saved"; } else { LogInfo("OutputInfo") << " JPTJethistograms will be saved to file:" << mOutputFile; } }
| JPTJetTester::~JPTJetTester | ( | ) |
Definition at line 403 of file JPTJetTester.cc.
{
}
| void JPTJetTester::analyze | ( | const edm::Event & | mEvent, |
| const edm::EventSetup & | mSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 415 of file JPTJetTester.cc.
References HiRecoJets_cff::caloTowers, deltaR(), reco::LeafCandidate::eta(), MonitorElement::Fill(), fillMatchHists(), edm::Event::getByLabel(), edm::Event::getManyByType(), getSumPt(), 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, mE, mE_3000, mE_80, mEBEne, mEBTime, mEEEne, mEETime, mEFirst, mEFirst_3000, mEFirst_80, mEmEnergyProfile, mEmTiming, mEta, mEtaFineBin, mEtaFineBin1m, mEtaFineBin1p, mEtaFineBin2m, mEtaFineBin2p, mEtaFineBin3m, mEtaFineBin3p, mEtaFirst, reco::MET::mEtSig(), mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_3000, mGenPt_80, mHadEnergyProfile, mHadTiming, mHBEne, mHBTime, mHEEne, mHETime, mHFEne, mHFTime, mHOEne, mHOTime, mInputCollection, mInputGenCollection, mJetEnergyProfile, mMass, mMass_3000, mMass_80, mMatchGenPtThreshold, mMjj, mMjj_3000, mMultiplicityElecInVertexInCalo, mMultiplicityElecInVertexOutCalo, mMultiplicityElecOutVertexInCalo, mMultiplicityMuonInVertexInCalo, mMultiplicityMuonInVertexOutCalo, mMultiplicityMuonOutVertexInCalo, mMultiplicityPionInVertexInCalo, mMultiplicityPionInVertexOutCalo, mMultiplicityPionOutVertexInCalo, mNJets1, mNJets2, mNJetsEtaC, mNJetsEtaF, mP, mP_3000, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_3000, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPtFractionElecInVertexInCalo, mPtFractionElecInVertexOutCalo, mPtFractionElecOutVertexInCalo, mPtFractionMuonInVertexInCalo, mPtFractionMuonInVertexOutCalo, mPtFractionMuonOutVertexInCalo, mPtFractionPionInVertexInCalo, mPtFractionPionInVertexOutCalo, mPtFractionPionOutVertexInCalo, mPthat_3000, mPthat_80, mRMatch, mRThreshold, mTurnOnEverything, njet, numberofevents, 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<JPTJetCollection> jptJets;
mEvent.getByLabel(mInputCollection, jptJets);
if (!jptJets.isValid()) return;
JPTJetCollection::const_iterator jet = jptJets->begin ();
int jetIndex = 0;
int nJet = 0;
int nJetF = 0;
int nJetC = 0;
for (; jet != jptJets->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 == jptJets->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();
}
double sumPtPionsInVertexInCalo = getSumPt(jet->getPionsInVertexInCalo());
double sumPtPionsInVertexOutCalo = getSumPt(jet->getPionsInVertexOutCalo());
double sumPtPionsOutVertexInCalo = getSumPt(jet->getPionsOutVertexInCalo());
double sumPtMuonsInVertexInCalo = getSumPt(jet->getMuonsInVertexInCalo());
double sumPtMuonsInVertexOutCalo = getSumPt(jet->getMuonsInVertexOutCalo());
double sumPtMuonsOutVertexInCalo = getSumPt(jet->getMuonsOutVertexInCalo());
double sumPtElecsInVertexInCalo = getSumPt(jet->getElecsInVertexInCalo());
double sumPtElecsInVertexOutCalo = getSumPt(jet->getElecsInVertexOutCalo());
double sumPtElecsOutVertexInCalo = getSumPt(jet->getElecsOutVertexInCalo());
if (jet->pt()!=0.){
mPtFractionPionInVertexInCalo->Fill(sumPtPionsInVertexInCalo/jet->pt());
mPtFractionPionInVertexOutCalo->Fill(sumPtPionsInVertexOutCalo/jet->pt());
mPtFractionPionOutVertexInCalo->Fill(sumPtPionsOutVertexInCalo/jet->pt());
mPtFractionMuonInVertexInCalo->Fill(sumPtMuonsInVertexInCalo/jet->pt());
mPtFractionMuonInVertexOutCalo->Fill(sumPtMuonsInVertexOutCalo/jet->pt());
mPtFractionMuonOutVertexInCalo->Fill(sumPtMuonsOutVertexInCalo/jet->pt());
mPtFractionElecInVertexInCalo->Fill(sumPtElecsInVertexInCalo/jet->pt());
mPtFractionElecInVertexOutCalo->Fill(sumPtElecsInVertexOutCalo/jet->pt());
mPtFractionElecOutVertexInCalo->Fill(sumPtElecsOutVertexInCalo/jet->pt());
//
mMultiplicityPionInVertexInCalo->Fill(jet->getPionsInVertexInCalo().size());
mMultiplicityPionInVertexOutCalo->Fill(jet->getPionsInVertexOutCalo().size());
mMultiplicityPionOutVertexInCalo->Fill(jet->getPionsOutVertexInCalo().size());
mMultiplicityMuonInVertexInCalo->Fill(jet->getMuonsInVertexInCalo().size());
mMultiplicityMuonInVertexOutCalo->Fill(jet->getMuonsInVertexOutCalo().size());
mMultiplicityMuonOutVertexInCalo->Fill(jet->getMuonsOutVertexInCalo().size());
mMultiplicityElecInVertexInCalo->Fill(jet->getElecsInVertexInCalo().size());
mMultiplicityElecInVertexOutCalo->Fill(jet->getElecsInVertexOutCalo().size());
mMultiplicityElecOutVertexInCalo->Fill(jet->getElecsOutVertexInCalo().size());
}
/*
if (mMaxEInEmTowers) mMaxEInEmTowers->Fill ( (jet->getCaloJetRef()->maxEInEmTowers());
if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->getCaloJetRef()->maxEInHadTowers());
if (mHadEnergyInHO) mHadEnergyInHO->Fill (jet->getCaloJetRef()->hadEnergyInHO());
if (mHadEnergyInHO_80) mHadEnergyInHO_80->Fill (jet->getCaloJetRef()->hadEnergyInHO());
if (mHadEnergyInHO_3000) mHadEnergyInHO_3000->Fill (jet->getCaloJetRef()->hadEnergyInHO());
if (mHadEnergyInHB) mHadEnergyInHB->Fill (jet->getCaloJetRef()->hadEnergyInHB());
if (mHadEnergyInHB_80) mHadEnergyInHB_80->Fill (jet->getCaloJetRef()->hadEnergyInHB());
if (mHadEnergyInHB_3000) mHadEnergyInHB_3000->Fill (jet->getCaloJetRef()->hadEnergyInHB());
if (mHadEnergyInHF) mHadEnergyInHF->Fill (jet->getCaloJetRef()->hadEnergyInHF());
if (mHadEnergyInHE) mHadEnergyInHE->Fill (jet->getCaloJetRef()->hadEnergyInHE());
if (mHadEnergyInHE_80) mHadEnergyInHE_80->Fill (jet->getCaloJetRef()->hadEnergyInHE());
if (mHadEnergyInHE_3000) mHadEnergyInHE_3000->Fill (jet->getCaloJetRef()->hadEnergyInHE());
if (mEmEnergyInEB) mEmEnergyInEB->Fill (jet->getCaloJetRef()->emEnergyInEB());
if (mEmEnergyInEB_80) mEmEnergyInEB_80->Fill (jet->getCaloJetRef()->emEnergyInEB());
if (mEmEnergyInEB_3000) mEmEnergyInEB_3000->Fill (jet->getCaloJetRef()->emEnergyInEB());
if (mEmEnergyInEE) mEmEnergyInEE->Fill (jet->getCaloJetRef()->emEnergyInEE());
if (mEmEnergyInEE_80) mEmEnergyInEE_80->Fill (jet->getCaloJetRef()->emEnergyInEE());
if (mEmEnergyInEE_3000) mEmEnergyInEE_3000->Fill (jet->getCaloJetRef()->emEnergyInEE());
if (mEmEnergyInHF) mEmEnergyInHF->Fill (jet->getCaloJetRef()->emEnergyInHF());
if (mEnergyFractionHadronic) mEnergyFractionHadronic->Fill (jet->getCaloJetRef()->energyFractionHadronic());
if (mEnergyFractionEm) mEnergyFractionEm->Fill (jet->getCaloJetRef()->emEnergyFraction());
if (mHFTotal) mHFTotal->Fill (jet->getCaloJetRef()->hadEnergyInHF()+jet->getCaloJetRef()->emEnergyInHF());
if (mHFTotal_80) mHFTotal_80->Fill (jet->getCaloJetRef()->hadEnergyInHF()+jet->getCaloJetRef()->emEnergyInHF());
if (mHFTotal_3000) mHFTotal_3000->Fill (jet->getCaloJetRef()->hadEnergyInHF()+jet->getCaloJetRef()->emEnergyInHF());
if (mHFLong) mHFLong->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5+jet->getCaloJetRef()->emEnergyInHF());
if (mHFLong_80) mHFLong_80->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5+jet->getCaloJetRef()->emEnergyInHF());
if (mHFLong_3000) mHFLong_3000->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5+jet->getCaloJetRef()->emEnergyInHF());
if (mHFShort) mHFShort->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5);
if (mHFShort_80) mHFShort_80->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5);
if (mHFShort_3000) mHFShort_3000->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5);
*/
/*
if (mN90) mN90->Fill (jet->getCaloJetRef()->n90());
*/
mJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->energy());
/*
mHadJetEnergyProfile->Fill (jet->eta(), jet->phi(),
jet->getCaloJetRef()->hadEnergyInHO()+
jet->getCaloJetRef()->hadEnergyInHB()+
jet->getCaloJetRef()->hadEnergyInHF()+
jet->getCaloJetRef()->hadEnergyInHE());
mEMJetEnergyProfile->Fill (jet->eta(), jet->phi(),
jet->getCaloJetRef()->emEnergyInEB()+
jet->getCaloJetRef()->emEnergyInEE()+
jet->getCaloJetRef()->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 ( JPTJetCollection::const_iterator jpt = jptJets->begin(); jpt != jptJets->end(); ++ jpt ) {
if ( jpt->pt() > ptStep ) njet++;
}
mNJets1->Fill( ptStep, njet );
}
for (int istep = 0; istep < 100; ++istep) {
int njet = 0;
float ptStep = (istep * (4000./100.));
for ( JPTJetCollection::const_iterator jpt = jptJets->begin(); jpt != jptJets->end(); ++ jpt ) {
if ( jpt->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 JPTJets 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*> > jptJetConstituents (jptJets->size());
// if (mRThreshold > 0) {
// }
// else {
// for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
// genJetConstituents [iGenJet] = jetMatching.getGenParticles ((*genJets) [iGenJet]);
// }
// for (unsigned iJPTJet = 0; iJPTJet < jptJets->size(); ++iJPTJet) {
// jptJetConstituents [iJPTJet] = jetMatching.getGenParticles ((*(*jptJets) [iJPTJet].getCaloJetRef()), 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 (jptJets->size() <= 0) continue; // no JPTJets - nothing to match
if (mRThreshold > 0) {
unsigned iJPTJetBest = 0;
double deltaRBest = 999.;
for (unsigned iJPTJet = 0; iJPTJet < jptJets->size(); ++iJPTJet) {
double dR = deltaR (genJet.eta(), genJet.phi(), (*jptJets) [iJPTJet].eta(), (*jptJets) [iJPTJet].phi());
if (deltaRBest < mRThreshold && dR < mRThreshold && genJet.pt() > 5.) {
/*
std::cout << "Yet another matched jet for GenJet pt=" << genJet.pt()
<< " previous JPTJet pt/dr: " << (*jptJets) [iJPTJetBest].pt() << '/' << deltaRBest
<< " new JPTJet pt/dr: " << (*jptJets) [iJPTJet].pt() << '/' << dR
<< std::endl;
*/
}
if (dR < deltaRBest) {
iJPTJetBest = iJPTJet;
deltaRBest = dR;
}
}
if (mTurnOnEverything.compare("yes")==0) {
mRMatch->Fill (logPtGen, genJet.eta(), deltaRBest);
}
if (deltaRBest < mRThreshold) { // Matched
fillMatchHists (genJet, (*jptJets) [iJPTJetBest]);
}
}
// else {
// unsigned iJPTJetBest = 0;
// double energyFractionBest = 0.;
// for (unsigned iJPTJet = 0; iJPTJet < jptJets->size(); ++iJPTJet) {
// double energyFraction = jetMatching.overlapEnergyFraction (genJetConstituents [iGenJet],
// jptJetConstituents [iJPTJet]);
// if (energyFraction > energyFractionBest) {
// iJPTJetBest = iJPTJet;
// energyFractionBest = energyFraction;
// }
// }
// if (mTurnOnEverything.compare("yes")==0) {
// mGenJetMatchEnergyFraction->Fill (logPtGen, genJet.eta(), energyFractionBest);
// }
// if (energyFractionBest > mGenEnergyFractionThreshold) { // good enough
// double reverseEnergyFraction = jetMatching.overlapEnergyFraction (jptJetConstituents [iJPTJetBest],
// genJetConstituents [iGenJet]);
// if (mTurnOnEverything.compare("yes")==0) {
// mReverseMatchEnergyFraction->Fill (logPtGen, genJet.eta(), reverseEnergyFraction);
// }
// if (reverseEnergyFraction > mReverseEnergyFractionThreshold) { // Matched
// fillMatchHists (genJet, (*jptJets) [iJPTJetBest]);
// }
// }
// } // else - mRThreshold
}
}
}
| void JPTJetTester::beginJob | ( | void | ) | [virtual] |
| void JPTJetTester::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 410 of file JPTJetTester.cc.
References mOutputFile.
{
if (!mOutputFile.empty() && &*edm::Service<DQMStore>()) edm::Service<DQMStore>()->save (mOutputFile);
}
| void JPTJetTester::fillMatchHists | ( | const reco::GenJet & | fGenJet, |
| const reco::JPTJet & | fJPTJet | ||
| ) | [private] |
Definition at line 869 of file JPTJetTester.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 PtJpt = fJPTJet.pt();
mMatchedGenJetsPt->Fill (logPtGen);
mMatchedGenJetsEta->Fill (logPtGen, fGenJet.eta());
double PtThreshold = 10.;
if (mTurnOnEverything.compare("yes")==0) {
mDeltaEta->Fill (logPtGen, fGenJet.eta(), fJPTJet.eta()-fGenJet.eta());
mDeltaPhi->Fill (logPtGen, fGenJet.eta(), fJPTJet.phi()-fGenJet.phi());
mEScale->Fill (logPtGen, fGenJet.eta(), fJPTJet.energy()/fGenJet.energy());
mlinEScale->Fill (fGenJet.pt(), fGenJet.eta(), fJPTJet.energy()/fGenJet.energy());
mDeltaE->Fill (logPtGen, fGenJet.eta(), fJPTJet.energy()-fGenJet.energy());
mEScaleFineBin->Fill (logPtGen, fGenJet.eta(), fJPTJet.energy()/fGenJet.energy());
if (fGenJet.pt()>PtThreshold) {
mEScale_pt10->Fill (logPtGen, fGenJet.eta(), fJPTJet.energy()/fGenJet.energy());
}
}
if (fJPTJet.pt() > PtThreshold) {
mDelEta->Fill (fGenJet.eta()-fJPTJet.eta());
mDelPhi->Fill (fGenJet.phi()-fJPTJet.phi());
mDelPt->Fill ((fGenJet.pt()-fJPTJet.pt())/fGenJet.pt());
}
if (fabs(fGenJet.eta())<1.3) {
mpTScaleB_s->Fill (log10(PtGen), PtJpt/PtGen);
mpTScaleB_d->Fill (log10(PtGen), PtJpt/PtGen);
if (PtGen>60.0 && PtGen<120.0) {
mpTScale1DB_60_120->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>200.0 && PtGen<300.0) {
mpTScale1DB_200_300->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>600.0 && PtGen<900.0) {
mpTScale1DB_600_900->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>2700.0 && PtGen<3500.0) {
mpTScale1DB_2700_3500->Fill (fJPTJet.pt()/fGenJet.pt());
}
}
if (fabs(fGenJet.eta())>1.3 && fabs(fGenJet.eta())<3.0) {
mpTScaleE_s->Fill (log10(PtGen), PtJpt/PtGen);
mpTScaleE_d->Fill (log10(PtGen), PtJpt/PtGen);
if (PtGen>60.0 && PtGen<120.0) {
mpTScale1DE_60_120->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>200.0 && PtGen<300.0) {
mpTScale1DE_200_300->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>600.0 && PtGen<900.0) {
mpTScale1DE_600_900->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>2700.0 && PtGen<3500.0) {
mpTScale1DE_2700_3500->Fill (fJPTJet.pt()/fGenJet.pt());
}
}
if (fabs(fGenJet.eta())>3.0 && fabs(fGenJet.eta())<5.0) {
mpTScaleF_s->Fill (log10(PtGen), PtJpt/PtGen);
mpTScaleF_d->Fill (log10(PtGen), PtJpt/PtGen);
if (PtGen>60.0 && PtGen<120.0) {
mpTScale1DF_60_120->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>200.0 && PtGen<300.0) {
mpTScale1DF_200_300->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>600.0 && PtGen<900.0) {
mpTScale1DF_600_900->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (PtGen>2700.0 && PtGen<3500.0) {
mpTScale1DF_2700_3500->Fill (fJPTJet.pt()/fGenJet.pt());
}
}
if (fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
mpTScale_60_120_s->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale_60_120_d->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale1D_60_120->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
mpTScale_200_300_s->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale_200_300_d->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale1D_200_300->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
mpTScale_600_900_s->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale_600_900_d->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale1D_600_900->Fill (fJPTJet.pt()/fGenJet.pt());
}
if (fGenJet.pt()>2700.0 && fGenJet.pt()<3500.0) {
mpTScale_2700_3500_s->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale_2700_3500_d->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
mpTScale1D_2700_3500->Fill (fJPTJet.pt()/fGenJet.pt());
}
}
| double JPTJetTester::getSumPt | ( | const reco::TrackRefVector & | tracks | ) | [private] |
Definition at line 981 of file JPTJetTester.cc.
References edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::end(), and reco::TrackBase::pt().
Referenced by analyze().
{
double sumpt = 0.;
for (reco::TrackRefVector::const_iterator itrack = tracks.begin(); itrack != tracks.end(); ++itrack){
const reco::Track& track = **itrack;
sumpt += track.pt();
}
return sumpt;
}
std::string JPTJetTester::inputCaloMETLabel_ [private] |
Definition at line 42 of file JPTJetTester.h.
std::string JPTJetTester::inputGenMETLabel_ [private] |
Definition at line 41 of file JPTJetTester.h.
edm::InputTag JPTJetTester::inputMETLabel_ [private] |
Definition at line 39 of file JPTJetTester.h.
MonitorElement* JPTJetTester::mAllGenJetsEta [private] |
Definition at line 173 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mAllGenJetsPt [private] |
Definition at line 171 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMET [private] |
Definition at line 248 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMET_3000 [private] |
Definition at line 249 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMETPhi [private] |
Definition at line 250 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMETSig [private] |
Definition at line 253 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMETSig_3000 [private] |
Definition at line 254 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMEx [private] |
Definition at line 244 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMEx_3000 [private] |
Definition at line 245 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMEy [private] |
Definition at line 246 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloMEy_3000 [private] |
Definition at line 247 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloSumET [private] |
Definition at line 251 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mCaloSumET_3000 [private] |
Definition at line 252 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mConstituents [private] |
Definition at line 70 of file JPTJetTester.h.
MonitorElement* JPTJetTester::mConstituents_3000 [private] |
Definition at line 72 of file JPTJetTester.h.
MonitorElement* JPTJetTester::mConstituents_80 [private] |
Definition at line 71 of file JPTJetTester.h.
MonitorElement* JPTJetTester::mDelEta [private] |
Definition at line 221 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mDelPhi [private] |
Definition at line 222 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mDelPt [private] |
Definition at line 223 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mDeltaE [private] |
Definition at line 182 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mDeltaEta [private] |
Definition at line 178 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mDeltaPhi [private] |
Definition at line 179 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mE [private] |
Definition at line 58 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mE_3000 [private] |
Definition at line 60 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mE_80 [private] |
Definition at line 59 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEBEne [private] |
Definition at line 265 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEBTime [private] |
Definition at line 266 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEEEne [private] |
Definition at line 267 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEETime [private] |
Definition at line 268 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEFirst [private] |
Definition at line 79 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEFirst_3000 [private] |
Definition at line 81 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEFirst_80 [private] |
Definition at line 80 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEmEnergyProfile [private] |
Definition at line 236 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEmTiming [private] |
Definition at line 74 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEScale [private] |
Definition at line 180 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mEScale_pt10 [private] |
Definition at line 184 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mEScaleFineBin [private] |
Definition at line 185 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mEta [private] |
Definition at line 48 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin [private] |
Definition at line 49 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin1m [private] |
Definition at line 53 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin1p [private] |
Definition at line 50 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin2m [private] |
Definition at line 54 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin2p [private] |
Definition at line 51 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin3m [private] |
Definition at line 55 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFineBin3p [private] |
Definition at line 52 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mEtaFirst [private] |
Definition at line 77 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
std::string JPTJetTester::METType_ [private] |
Definition at line 40 of file JPTJetTester.h.
double JPTJetTester::mGenEnergyFractionThreshold [private] |
Definition at line 227 of file JPTJetTester.h.
MonitorElement* JPTJetTester::mGenEta [private] |
Definition at line 160 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mGenEtaFirst [private] |
Definition at line 167 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 175 of file JPTJetTester.h.
Referenced by JPTJetTester().
MonitorElement* JPTJetTester::mGenPhi [private] |
Definition at line 161 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mGenPhiFirst [private] |
Definition at line 168 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mGenPt [private] |
Definition at line 162 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mGenPt_3000 [private] |
Definition at line 164 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mGenPt_80 [private] |
Definition at line 163 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHadEnergyProfile [private] |
Definition at line 235 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHadTiming [private] |
Definition at line 73 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHBEne [private] |
Definition at line 257 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHBTime [private] |
Definition at line 258 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHEEne [private] |
Definition at line 259 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHETime [private] |
Definition at line 260 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHFEne [private] |
Definition at line 263 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHFTime [private] |
Definition at line 264 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHOEne [private] |
Definition at line 261 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mHOTime [private] |
Definition at line 262 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
edm::InputTag JPTJetTester::mInputCollection [private] |
Definition at line 36 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 37 of file JPTJetTester.h.
Referenced by analyze().
MonitorElement* JPTJetTester::mJetEnergyProfile [private] |
Definition at line 237 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mlinEScale [private] |
Definition at line 181 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mMass [private] |
Definition at line 67 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mMass_3000 [private] |
Definition at line 69 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mMass_80 [private] |
Definition at line 68 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mMatchedGenJetsEta [private] |
Definition at line 174 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mMatchedGenJetsPt [private] |
Definition at line 172 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
double JPTJetTester::mMatchGenPtThreshold [private] |
Definition at line 226 of file JPTJetTester.h.
Referenced by analyze().
MonitorElement* JPTJetTester::mMjj [private] |
Definition at line 93 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mMjj_3000 [private] |
Definition at line 94 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 118 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 119 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 120 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 114 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 115 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 116 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 110 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 111 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 112 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mNJets1 [private] |
Definition at line 89 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mNJets2 [private] |
Definition at line 90 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mNJetsEtaC [private] |
Definition at line 86 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mNJetsEtaF [private] |
Definition at line 87 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
std::string JPTJetTester::mOutputFile [private] |
Definition at line 38 of file JPTJetTester.h.
Referenced by endJob(), and JPTJetTester().
MonitorElement* JPTJetTester::mP [private] |
Definition at line 61 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mP_3000 [private] |
Definition at line 63 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mP_80 [private] |
Definition at line 62 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPhi [private] |
Definition at line 56 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPhiFineBin [private] |
Definition at line 57 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPhiFirst [private] |
Definition at line 78 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPt [private] |
Definition at line 64 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPt_3000 [private] |
Definition at line 66 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPt_80 [private] |
Definition at line 65 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPtFirst [private] |
Definition at line 82 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPtFirst_3000 [private] |
Definition at line 84 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPtFirst_80 [private] |
Definition at line 83 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 105 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 106 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 107 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 101 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 102 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 103 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 97 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 98 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
Definition at line 99 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPthat_3000 [private] |
Definition at line 157 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mPthat_80 [private] |
Definition at line 156 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1D_200_300 [private] |
Definition at line 217 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1D_2700_3500 [private] |
Definition at line 219 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1D_600_900 [private] |
Definition at line 218 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1D_60_120 [private] |
Definition at line 216 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DB_200_300 [private] |
Definition at line 207 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
Definition at line 213 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DB_600_900 [private] |
Definition at line 210 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DB_60_120 [private] |
Definition at line 204 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DE_200_300 [private] |
Definition at line 208 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
Definition at line 214 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DE_600_900 [private] |
Definition at line 211 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DE_60_120 [private] |
Definition at line 205 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DF_200_300 [private] |
Definition at line 209 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
Definition at line 215 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DF_600_900 [private] |
Definition at line 212 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale1DF_60_120 [private] |
Definition at line 206 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_200_300_d [private] |
Definition at line 200 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_200_300_s [private] |
Definition at line 195 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_2700_3500_d [private] |
Definition at line 202 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_2700_3500_s [private] |
Definition at line 197 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_600_900_d [private] |
Definition at line 201 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_600_900_s [private] |
Definition at line 196 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_60_120_d [private] |
Definition at line 199 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScale_60_120_s [private] |
Definition at line 194 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScaleB_d [private] |
Definition at line 190 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScaleB_s [private] |
Definition at line 187 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScaleE_d [private] |
Definition at line 191 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScaleE_s [private] |
Definition at line 188 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScaleF_d [private] |
Definition at line 192 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::mpTScaleF_s [private] |
Definition at line 189 of file JPTJetTester.h.
Referenced by fillMatchHists(), and JPTJetTester().
double JPTJetTester::mReverseEnergyFractionThreshold [private] |
Definition at line 228 of file JPTJetTester.h.
Definition at line 176 of file JPTJetTester.h.
Referenced by JPTJetTester().
MonitorElement* JPTJetTester::mRMatch [private] |
Definition at line 177 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
double JPTJetTester::mRThreshold [private] |
Definition at line 229 of file JPTJetTester.h.
Referenced by analyze().
std::string JPTJetTester::mTurnOnEverything [private] |
Definition at line 232 of file JPTJetTester.h.
Referenced by analyze(), fillMatchHists(), and JPTJetTester().
MonitorElement* JPTJetTester::numberofevents [private] |
Definition at line 45 of file JPTJetTester.h.
Referenced by analyze(), and JPTJetTester().
1.7.3