PFJetTester Class Reference

#include <PFJetTester.h>

Inheritance diagram for PFJetTester:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()
	PFJetTester (const edm::ParameterSet &)
	~PFJetTester ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
void	fillMatchHists (const reco::GenJet &fGenJet, const reco::PFJet &fPFJet, std::vector< reco::Vertex > goodVertices)

Private Attributes
std::string	inputCaloMETLabel_
std::string	inputGenMETLabel_
edm::InputTag	inputMETLabel_
std::string	JetCorrectionService
MonitorElement *	mAllGenJetsEta
MonitorElement *	mAllGenJetsPt
MonitorElement *	mCaloMET
MonitorElement *	mCaloMET_3000
MonitorElement *	mCaloMETPhi
MonitorElement *	mCaloMETSig
MonitorElement *	mCaloMETSig_3000
MonitorElement *	mCaloMEx
MonitorElement *	mCaloMEx_3000
MonitorElement *	mCaloMEy
MonitorElement *	mCaloMEy_3000
MonitorElement *	mCaloSumET
MonitorElement *	mCaloSumET_3000
MonitorElement *	mChargedEmEnergy_3000
MonitorElement *	mChargedEmEnergy_80
MonitorElement *	mChargedEmEnergyFraction_B
MonitorElement *	mChargedEmEnergyFraction_E
MonitorElement *	mChargedEmEnergyFraction_F
MonitorElement *	mChargedHadronEnergy_3000
MonitorElement *	mChargedHadronEnergy_80
MonitorElement *	mChargedHadronEnergyFraction_B
MonitorElement *	mChargedHadronEnergyFraction_E
MonitorElement *	mChargedHadronEnergyFraction_F
MonitorElement *	mConstituents
MonitorElement *	mConstituents_3000
MonitorElement *	mConstituents_80
MonitorElement *	mCorrJetEta
MonitorElement *	mCorrJetPhi
MonitorElement *	mCorrJetPt
MonitorElement *	mCorrJetPt_3000
MonitorElement *	mCorrJetPt_80
MonitorElement *	mDelEta
MonitorElement *	mDelPhi
MonitorElement *	mDelPt
MonitorElement *	mDeltaE
MonitorElement *	mDeltaEta
MonitorElement *	mDeltaPhi
MonitorElement *	mE
MonitorElement *	mE_3000
MonitorElement *	mE_80
MonitorElement *	mEBEne
MonitorElement *	mEBTime
MonitorElement *	mEEEne
MonitorElement *	mEETime
MonitorElement *	mEFirst
MonitorElement *	mEFirst_3000
MonitorElement *	mEFirst_80
MonitorElement *	mElectronEnergy
MonitorElement *	mElectronEnergy_3000
MonitorElement *	mElectronEnergy_80
MonitorElement *	mEmEnergyInHF
MonitorElement *	mEmEnergyInHF_3000
MonitorElement *	mEmEnergyInHF_80
MonitorElement *	mEmEnergyProfile
MonitorElement *	mEmTiming
MonitorElement *	mEScale
MonitorElement *	mEScale_pt10
MonitorElement *	mEScaleFineBin
MonitorElement *	mEta
MonitorElement *	mEtaFineBin
MonitorElement *	mEtaFineBin1m
MonitorElement *	mEtaFineBin1p
MonitorElement *	mEtaFineBin2m
MonitorElement *	mEtaFineBin2p
MonitorElement *	mEtaFineBin3m
MonitorElement *	mEtaFineBin3p
MonitorElement *	mEtaFirst
std::string	METType_
double	mGenEnergyFractionThreshold
MonitorElement *	mGenEta
MonitorElement *	mGenEtaFirst
MonitorElement *	mGenJetMatchEnergyFraction
MonitorElement *	mGenPhi
MonitorElement *	mGenPhiFirst
MonitorElement *	mGenPt
MonitorElement *	mGenPt_3000
MonitorElement *	mGenPt_80
MonitorElement *	mHadEnergyInHF
MonitorElement *	mHadEnergyInHF_3000
MonitorElement *	mHadEnergyInHF_80
MonitorElement *	mHadEnergyProfile
MonitorElement *	mHadTiming
MonitorElement *	mHBEne
MonitorElement *	mHBTime
MonitorElement *	mHEEne
MonitorElement *	mHETime
MonitorElement *	mHFEne
MonitorElement *	mHFTime
MonitorElement *	mHOEne
MonitorElement *	mHOTime
edm::InputTag	mInputCollection
edm::InputTag	mInputGenCollection
MonitorElement *	mjetArea
MonitorElement *	mJetEnergyProfile
MonitorElement *	mlinEScale
MonitorElement *	mMass
MonitorElement *	mMass_3000
MonitorElement *	mMass_80
MonitorElement *	mMatchedGenJetsEta
MonitorElement *	mMatchedGenJetsPt
double	mMatchGenPtThreshold
MonitorElement *	mMjj
MonitorElement *	mMjj_3000
MonitorElement *	mMuonEnergy
MonitorElement *	mMuonEnergy_3000
MonitorElement *	mMuonEnergy_80
MonitorElement *	mNeutralEmEnergy_3000
MonitorElement *	mNeutralEmEnergy_80
MonitorElement *	mNeutralEmEnergyFraction_B
MonitorElement *	mNeutralEmEnergyFraction_E
MonitorElement *	mNeutralEmEnergyFraction_F
MonitorElement *	mNeutralHadronEnergy_3000
MonitorElement *	mNeutralHadronEnergy_80
MonitorElement *	mNeutralHadronEnergyFraction_B
MonitorElement *	mNeutralHadronEnergyFraction_E
MonitorElement *	mNeutralHadronEnergyFraction_F
MonitorElement *	mNJets1
MonitorElement *	mNJets2
MonitorElement *	mNJetsEtaC
MonitorElement *	mNJetsEtaF
MonitorElement *	mNJetsEtaF_30
std::string	mOutputFile
MonitorElement *	mP
MonitorElement *	mP_3000
MonitorElement *	mP_80
MonitorElement *	mPhi
MonitorElement *	mPhiFineBin
MonitorElement *	mPhiFirst
MonitorElement *	mPhotonEnergy
MonitorElement *	mPhotonEnergy_3000
MonitorElement *	mPhotonEnergy_80
MonitorElement *	mPt
MonitorElement *	mPt_3000
MonitorElement *	mPt_80
MonitorElement *	mPtFirst
MonitorElement *	mPtFirst_3000
MonitorElement *	mPtFirst_80
MonitorElement *	mPthat_3000
MonitorElement *	mPthat_80
MonitorElement *	mpTRatio
MonitorElement *	mpTRatio_1500_3500_d
MonitorElement *	mpTRatio_200_600_d
MonitorElement *	mpTRatio_30_200_d
MonitorElement *	mpTRatio_600_1500_d
MonitorElement *	mpTRatioB_d
MonitorElement *	mpTRatioE_d
MonitorElement *	mpTRatioF_d
MonitorElement *	mpTResponse
MonitorElement *	mpTResponse_1500_3500_d
MonitorElement *	mpTResponse_200_600_d
MonitorElement *	mpTResponse_30_200_d
MonitorElement *	mpTResponse_30_d
MonitorElement *	mpTResponse_600_1500_d
MonitorElement *	mpTResponse_nvtx_0_5
MonitorElement *	mpTResponse_nvtx_10_15
MonitorElement *	mpTResponse_nvtx_15_20
MonitorElement *	mpTResponse_nvtx_20_30
MonitorElement *	mpTResponse_nvtx_30_inf
MonitorElement *	mpTResponse_nvtx_5_10
MonitorElement *	mpTResponseB_d
MonitorElement *	mpTResponseE_d
MonitorElement *	mpTResponseF_d
MonitorElement *	mpTScale1D_1500_3500
MonitorElement *	mpTScale1D_200_600
MonitorElement *	mpTScale1D_30_200
MonitorElement *	mpTScale1D_600_1500
MonitorElement *	mpTScale1DB_1500_3500
MonitorElement *	mpTScale1DB_200_600
MonitorElement *	mpTScale1DB_30_200
MonitorElement *	mpTScale1DB_600_1500
MonitorElement *	mpTScale1DE_1500_3500
MonitorElement *	mpTScale1DE_200_600
MonitorElement *	mpTScale1DE_30_200
MonitorElement *	mpTScale1DE_600_1500
MonitorElement *	mpTScale1DF_1500_3500
MonitorElement *	mpTScale1DF_200_600
MonitorElement *	mpTScale1DF_30_200
MonitorElement *	mpTScale1DF_600_1500
MonitorElement *	mpTScale_1500_3500_d
MonitorElement *	mpTScale_1500_3500_s
MonitorElement *	mpTScale_200_600_d
MonitorElement *	mpTScale_200_600_s
MonitorElement *	mpTScale_30_200_d
MonitorElement *	mpTScale_30_200_s
MonitorElement *	mpTScale_600_1500_d
MonitorElement *	mpTScale_600_1500_s
MonitorElement *	mpTScale_a
MonitorElement *	mpTScale_a_nvtx_0_5
MonitorElement *	mpTScale_a_nvtx_10_15
MonitorElement *	mpTScale_a_nvtx_15_20
MonitorElement *	mpTScale_a_nvtx_20_30
MonitorElement *	mpTScale_a_nvtx_30_inf
MonitorElement *	mpTScale_a_nvtx_5_10
MonitorElement *	mpTScale_b
MonitorElement *	mpTScale_b_nvtx_0_5
MonitorElement *	mpTScale_b_nvtx_10_15
MonitorElement *	mpTScale_b_nvtx_15_20
MonitorElement *	mpTScale_b_nvtx_20_30
MonitorElement *	mpTScale_b_nvtx_30_inf
MonitorElement *	mpTScale_b_nvtx_5_10
MonitorElement *	mpTScale_c
MonitorElement *	mpTScale_c_nvtx_0_5
MonitorElement *	mpTScale_c_nvtx_10_15
MonitorElement *	mpTScale_c_nvtx_15_20
MonitorElement *	mpTScale_c_nvtx_20_30
MonitorElement *	mpTScale_c_nvtx_30_inf
MonitorElement *	mpTScale_c_nvtx_5_10
MonitorElement *	mpTScale_nvtx_0_5
MonitorElement *	mpTScale_nvtx_10_15
MonitorElement *	mpTScale_nvtx_15_20
MonitorElement *	mpTScale_nvtx_20_30
MonitorElement *	mpTScale_nvtx_30_inf
MonitorElement *	mpTScale_nvtx_5_10
MonitorElement *	mpTScale_pT
MonitorElement *	mpTScaleB_d
MonitorElement *	mpTScaleB_s
MonitorElement *	mpTScaleE_d
MonitorElement *	mpTScaleE_s
MonitorElement *	mpTScaleF_d
MonitorElement *	mpTScaleF_s
MonitorElement *	mpTScalePhiB_d
MonitorElement *	mpTScalePhiE_d
MonitorElement *	mpTScalePhiF_d
double	mReverseEnergyFractionThreshold
MonitorElement *	mReverseMatchEnergyFraction
MonitorElement *	mRMatch
double	mRThreshold
std::string	mTurnOnEverything
MonitorElement *	numberofevents
MonitorElement *	nvtx_0_30
MonitorElement *	nvtx_0_60

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 23 of file PFJetTester.h.

Constructor & Destructor Documentation

PFJetTester::PFJetTester

(

const edm::ParameterSet &

iConfig

)

Definition at line 49 of file PFJetTester.cc.

References DQMStore::book1D(), DQMStore::bookProfile(), MuonErrorMatrixValues_cff::etaRange, edm::InputTag::label(), mChargedEmEnergy_3000, mChargedEmEnergy_80, mChargedEmEnergyFraction_B, mChargedEmEnergyFraction_E, mChargedEmEnergyFraction_F, mChargedHadronEnergy_3000, mChargedHadronEnergy_80, mChargedHadronEnergyFraction_B, mChargedHadronEnergyFraction_E, mChargedHadronEnergyFraction_F, mConstituents, mConstituents_80, mCorrJetEta, mCorrJetPhi, mCorrJetPt, mCorrJetPt_80, mDelEta, mDelPhi, mDelPt, mDeltaEta, mDeltaPhi, mE, mE_80, mEmEnergyInHF, mEmEnergyInHF_3000, mEmEnergyInHF_80, mEmTiming, mEScale_pt10, mEScaleFineBin, mEta, mEtaFineBin, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_80, mHadEnergyInHF, mHadEnergyInHF_3000, mHadEnergyInHF_80, mHadTiming, mInputCollection, mjetArea, mMass, mMass_80, mMjj, mMjj_3000, mNeutralEmEnergy_3000, mNeutralEmEnergy_80, mNeutralEmEnergyFraction_B, mNeutralEmEnergyFraction_E, mNeutralEmEnergyFraction_F, mNeutralHadronEnergy_3000, mNeutralHadronEnergy_80, mNeutralHadronEnergyFraction_B, mNeutralHadronEnergyFraction_E, mNeutralHadronEnergyFraction_F, mNJets1, mNJets2, mNJetsEtaC, mNJetsEtaF, mNJetsEtaF_30, mOutputFile, mP, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPthat_3000, mPthat_80, mpTRatio, mpTRatio_1500_3500_d, mpTRatio_200_600_d, mpTRatio_30_200_d, mpTRatio_600_1500_d, mpTRatioB_d, mpTRatioE_d, mpTRatioF_d, mpTResponse, mpTResponse_1500_3500_d, mpTResponse_200_600_d, mpTResponse_30_200_d, mpTResponse_30_d, mpTResponse_600_1500_d, mpTResponse_nvtx_0_5, mpTResponse_nvtx_10_15, mpTResponse_nvtx_15_20, mpTResponse_nvtx_20_30, mpTResponse_nvtx_30_inf, mpTResponse_nvtx_5_10, mpTResponseB_d, mpTResponseE_d, mpTResponseF_d, mpTScale1DB_1500_3500, mpTScale1DB_200_600, mpTScale1DB_30_200, mpTScale1DB_600_1500, mpTScale1DE_1500_3500, mpTScale1DE_200_600, mpTScale1DE_30_200, mpTScale1DE_600_1500, mpTScale1DF_1500_3500, mpTScale1DF_200_600, mpTScale1DF_30_200, mpTScale1DF_600_1500, mpTScale_1500_3500_d, mpTScale_200_600_d, mpTScale_30_200_d, mpTScale_600_1500_d, mpTScale_a, mpTScale_a_nvtx_0_5, mpTScale_a_nvtx_10_15, mpTScale_a_nvtx_15_20, mpTScale_a_nvtx_20_30, mpTScale_a_nvtx_30_inf, mpTScale_a_nvtx_5_10, mpTScale_b, mpTScale_b_nvtx_0_5, mpTScale_b_nvtx_10_15, mpTScale_b_nvtx_15_20, mpTScale_b_nvtx_20_30, mpTScale_b_nvtx_30_inf, mpTScale_b_nvtx_5_10, mpTScale_c, mpTScale_c_nvtx_0_5, mpTScale_c_nvtx_10_15, mpTScale_c_nvtx_15_20, mpTScale_c_nvtx_20_30, mpTScale_c_nvtx_30_inf, mpTScale_c_nvtx_5_10, mpTScale_nvtx_0_5, mpTScale_nvtx_10_15, mpTScale_nvtx_15_20, mpTScale_nvtx_20_30, mpTScale_nvtx_30_inf, mpTScale_nvtx_5_10, mpTScale_pT, mpTScaleB_d, mpTScaleE_d, mpTScaleF_d, mpTScalePhiB_d, mpTScalePhiE_d, mpTScalePhiF_d, mTurnOnEverything, numberofevents, nvtx_0_30, nvtx_0_60, 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")),
    JetCorrectionService  (iConfig.getParameter<std::string>  ("JetCorrectionService"  )),
    mTurnOnEverything (iConfig.getUntrackedParameter<std::string>("TurnOnEverything",""))
{
    numberofevents
    = mEta = mEtaFineBin = mPhi = mPhiFineBin = mE = mE_80
    = mP = mP_80 = mPt = mPt_80
    = mMass = mMass_80 = mConstituents = mConstituents_80
    = mEtaFirst = mPhiFirst = mPtFirst = mPtFirst_80 = mPtFirst_3000
    = mMjj = mMjj_3000 = mDelEta = mDelPhi = mDelPt 
      //    = mMaxEInEmTowers = mMaxEInHadTowers 
      //    = mHadEnergyInHO = mHadEnergyInHB = mHadEnergyInHE 
      = mHadEnergyInHF = mHadEnergyInHF_80 = mHadEnergyInHF_3000
      //    = mHadEnergyInHO_80 = mHadEnergyInHB_80 = mHadEnergyInHE_80 
      //    = mHadEnergyInHO_3000 = mHadEnergyInHB_3000 = mHadEnergyInHE_3000 
      //    = mEmEnergyInEB = mEmEnergyInEE
      = mEmEnergyInHF = mEmEnergyInHF_80 = mEmEnergyInHF_3000
      //    = mEmEnergyInEB_80 = mEmEnergyInEE_80
      //    = mEmEnergyInEB_3000 = mEmEnergyInEE_3000
      //    = mEnergyFractionHadronic = mEnergyFractionEm 
      //    = mHFLong = mHFTotal = mHFLong_80 = mHFLong_3000 = mHFShort = mHFShort_80 = mHFShort_3000
      //    = mN90
    = mChargedEmEnergy_80 = mChargedHadronEnergy_80 = mNeutralEmEnergy_80 = mNeutralHadronEnergy_80
    = mChargedEmEnergy_3000 = mChargedHadronEnergy_3000 = mNeutralEmEnergy_3000 = mNeutralHadronEnergy_3000
    = mChargedEmEnergyFraction_B = mChargedEmEnergyFraction_E = mChargedEmEnergyFraction_F 
    = mChargedHadronEnergyFraction_B = mChargedHadronEnergyFraction_E = mChargedHadronEnergyFraction_F 
    = mNeutralEmEnergyFraction_B = mNeutralEmEnergyFraction_E = mNeutralEmEnergyFraction_F
    = mNeutralHadronEnergyFraction_B = mNeutralHadronEnergyFraction_E = mNeutralHadronEnergyFraction_F
      //= 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
     = mpTScalePhiB_d = mpTScalePhiE_d = mpTScalePhiF_d
      //= mpTScale_30_200_s = mpTScale_200_600_s = mpTScale_600_1500_s = mpTScale_1500_3500_s
    = mpTScale_30_200_d = mpTScale_200_600_d = mpTScale_600_1500_d = mpTScale_1500_3500_d
      
    = mpTScale1DB_30_200    = mpTScale1DE_30_200    = mpTScale1DF_30_200 
    = mpTScale1DB_200_600   = mpTScale1DE_200_600   = mpTScale1DF_200_600 
    = mpTScale1DB_600_1500   = mpTScale1DE_600_1500   = mpTScale1DF_600_1500 
    = mpTScale1DB_1500_3500 = mpTScale1DE_1500_3500 = mpTScale1DF_1500_3500
    /*
    = mpTScale1D_30_200 = mpTScale1D_200_600 = mpTScale1D_600_1500 = mpTScale1D_1500_3500
    = mHBEne = mHBTime = mHEEne = mHETime = mHFEne = mHFTime = mHOEne = mHOTime
    = mEBEne = mEBTime = mEEEne = mEETime 
      */
    = mPthat_80 = mPthat_3000

      //Corr Jet
    = mCorrJetPt =mCorrJetPt_80 =mCorrJetEta =mCorrJetPhi =mpTRatio =mpTResponse
      = mpTRatioB_d = mpTRatioE_d = mpTRatioF_d
      = mpTRatio_30_200_d = mpTRatio_200_600_d = mpTRatio_600_1500_d = mpTRatio_1500_3500_d
      = mpTResponseB_d = mpTResponseE_d = mpTResponseF_d
      = mpTResponse_30_200_d = mpTResponse_200_600_d = mpTResponse_600_1500_d = mpTResponse_1500_3500_d
      = mpTResponse_30_d =mjetArea
      = nvtx_0_30 = nvtx_0_60
      = mpTResponse_nvtx_0_5 = mpTResponse_nvtx_5_10 =mpTResponse_nvtx_10_15 
      = mpTResponse_nvtx_15_20 = mpTResponse_nvtx_20_30 = mpTResponse_nvtx_30_inf
      = mpTScale_a_nvtx_0_5 = mpTScale_b_nvtx_0_5 = mpTScale_c_nvtx_0_5 
      = mpTScale_a_nvtx_5_10 = mpTScale_b_nvtx_5_10 = mpTScale_c_nvtx_5_10
      = mpTScale_a_nvtx_10_15 = mpTScale_b_nvtx_10_15 = mpTScale_c_nvtx_10_15
      = mpTScale_a_nvtx_15_20 = mpTScale_b_nvtx_15_20 = mpTScale_c_nvtx_15_20
      = mpTScale_a_nvtx_20_30 = mpTScale_b_nvtx_20_30 = mpTScale_c_nvtx_20_30
      = mpTScale_a_nvtx_30_inf = mpTScale_b_nvtx_30_inf = mpTScale_c_nvtx_30_inf
      = mpTScale_nvtx_0_5 = mpTScale_nvtx_5_10 = mpTScale_nvtx_10_15 
      = mpTScale_nvtx_15_20 = mpTScale_nvtx_20_30 = mpTScale_nvtx_30_inf
      = mNJetsEtaF_30
      = mpTScale_a = mpTScale_b = mpTScale_c = mpTScale_pT
      = 0;
  
  DQMStore* dbe = &*edm::Service<DQMStore>();
  if (dbe) {
    dbe->setCurrentFolder("JetMET/RecoJetsV/PFJetTask_" + mInputCollection.label());
    //
    numberofevents    = dbe->book1D("numberofevents","numberofevents", 3, 0 , 2);
    //
    mEta              = dbe->book1D("Eta", "Eta", 120, -6, 6); 
    mEtaFineBin       = dbe->book1D("EtaFineBin_Pt10", "EtaFineBin_Pt10", 600, -6, 6);
    /* 
    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, 4000);  
    //
    mP                = dbe->book1D("P", "P", 100, 0, 500); 
    mP_80             = dbe->book1D("P_80", "P_80", 100, 0, 4000); 
    //
    mPt               = dbe->book1D("Pt", "Pt", 100, 0, 150); 
    mPt_80            = dbe->book1D("Pt_80", "Pt_80", 100, 0, 4000);  
    //
    mMass             = dbe->book1D("Mass", "Mass", 100, 0, 200); 
    mMass_80          = dbe->book1D("Mass_80", "Mass_80", 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", 120, -6, 6); 
    mPhiFirst         = dbe->book1D("PhiFirst", "PhiFirst", 70, -3.5, 3.5);      
    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, 2500); 
    mHadEnergyInHF_80    = dbe->book1D("HadEnergyInHF_80", "HadEnergyInHF_80", 100, 0, 3000); 
    mHadEnergyInHF_3000    = dbe->book1D("HadEnergyInHF_3000", "HadEnergyInHF_3000", 100, 0, 1800); 
    //    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", 100, -20, 450); 
    mEmEnergyInHF_80     = dbe->book1D("EmEnergyInHF_80", "EmEnergyInHF_80", 100, -20, 440); 
    mEmEnergyInHF_3000     = dbe->book1D("EmEnergyInHF_3000", "EmEnergyInHF_3000", 100, -20, 190); 
    //    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); 
    //

    mChargedEmEnergy_80 = dbe->book1D("ChargedEmEnergy_80","ChargedEmEnergy_80",100,0.,500.);
    mChargedHadronEnergy_80 = dbe->book1D("ChargedHadronEnergy_80","ChargedHadronEnergy_80",100,0.,500.);
    mNeutralEmEnergy_80 = dbe->book1D("NeutralEmEnergy_80","NeutralEmEnergy_80",100,0.,500.);
    mNeutralHadronEnergy_80 = dbe->book1D("NeutralHadronEnergy_80","NeutralHadronEnergy_80",100,0.,500.);

    mChargedEmEnergy_3000 = dbe->book1D("ChargedEmEnergy_3000","ChargedEmEnergy_3000",100,0.,2000.);
    mChargedHadronEnergy_3000 = dbe->book1D("ChargedHadronEnergy_3000","ChargedHadronEnergy_3000",100,0.,2000.);
    mNeutralEmEnergy_3000 = dbe->book1D("NeutralEmEnergy_3000","NeutralEmEnergy_3000",100,0.,2000.);
    mNeutralHadronEnergy_3000 = dbe->book1D("NeutralHadronEnergy_3000","NeutralHadronEnergy_3000",100,0.,2000.);

    mChargedEmEnergyFraction_B = dbe->book1D("ChargedEmEnergyFraction_B","ChargedEmEnergyFraction_B",120,-0.1,1.1);
    mChargedEmEnergyFraction_E = dbe->book1D("ChargedEmEnergyFraction_E","ChargedEmEnergyFraction_E",120,-0.1,1.1);
    mChargedEmEnergyFraction_F = dbe->book1D("ChargedEmEnergyFraction_F","ChargedEmEnergyFraction_F",120,-0.1,1.1);
    mChargedHadronEnergyFraction_B = dbe->book1D("ChargedHadronEnergyFraction_B","ChargedHadronEnergyFraction_B",120,-0.1,1.1);
    mChargedHadronEnergyFraction_E = dbe->book1D("ChargedHadronEnergyFraction_E","ChargedHadronEnergyFraction_E",120,-0.1,1.1);
    mChargedHadronEnergyFraction_F = dbe->book1D("ChargedHadronEnergyFraction_F","ChargedHadronEnergyFraction_F",120,-0.1,1.1);
    mNeutralEmEnergyFraction_B = dbe->book1D("NeutralEmEnergyFraction_B","NeutralEmEnergyFraction_B",120,-0.1,1.1);
    mNeutralEmEnergyFraction_E = dbe->book1D("NeutralEmEnergyFraction_E","NeutralEmEnergyFraction_E",120,-0.1,1.1);
    mNeutralEmEnergyFraction_F = dbe->book1D("NeutralEmEnergyFraction_F","NeutralEmEnergyFraction_F",120,-0.1,1.1);
    mNeutralHadronEnergyFraction_B = dbe->book1D("NeutralHadronEnergyFraction_B","NeutralHadronEnergyFraction_B",120,-0.1,1.1);
    mNeutralHadronEnergyFraction_E = dbe->book1D("NeutralHadronEnergyFraction_E","NeutralHadronEnergyFraction_E",120,-0.1,1.1);
    mNeutralHadronEnergyFraction_F = dbe->book1D("NeutralHadronEnergyFraction_F","NeutralHadronEnergyFraction_F",120,-0.1,1.1);

    mGenEta           = dbe->book1D("GenEta", "GenEta", 120, -6, 6);
    mGenPhi           = dbe->book1D("GenPhi", "GenPhi", 70, -3.5, 3.5);
    mGenPt            = dbe->book1D("GenPt", "GenPt", 100, 0, 150);
    mGenPt_80         = dbe->book1D("GenPt_80", "GenPt_80", 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);
    mNJetsEtaF_30        = dbe->book1D("NJetsEtaF_Pt30", "NJetsEtaF_Pt30", 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, 0.0, 1000.0); 
    mPthat_3000          = dbe->book1D("Pthat_3000", "Pthat_3000", 100, 1000.0, 4000.0); 

    //Corr
    mCorrJetPt  = dbe->book1D("CorrPt", "CorrPt", 100, 0, 150);
    mCorrJetPt_80 = dbe->book1D("CorrPt_80", "CorrPt_80", 100, 0, 4000);
    mCorrJetEta = dbe->book1D("CorrEta", "CorrEta", 120, -6, 6);
    mCorrJetPhi = dbe->book1D("CorrPhi", "CorrPhi", 70, -3.5, 3.5);

    mjetArea = dbe->book1D("jetArea","jetArea",26,-0.5,12.5);

    //nvtx
    nvtx_0_30 = dbe->book1D("nvtx_0_30","nvtx_0_30",31,-0.5,30.5);
    nvtx_0_60 = dbe->book1D("nvtx_0_60","nvtx_0_60",61,-0.5,60.5);

    //pT scale with nvtx
    mpTScale_a_nvtx_0_5 = dbe->book1D("mpTScale_a_nvtx_0_5", "pTScale_a_nvtx_0_5_0<|eta|<1.3_60_120",100, 0, 2);
    mpTScale_b_nvtx_0_5 = dbe->book1D("mpTScale_b_nvtx_0_5", "pTScale_b_nvtx_0_5_0<|eta|<1.3_200_300",100, 0, 2);
    mpTScale_c_nvtx_0_5 = dbe->book1D("mpTScale_c_nvtx_0_5", "pTScale_c_nvtx_0_5_0<|eta|<1.3_600_900",100, 0, 2);
    mpTScale_a_nvtx_5_10 = dbe->book1D("mpTScale_a_nvtx_5_10", "pTScale_a_nvtx_5_10_0<|eta|<1.3_60_120",100, 0, 2);
    mpTScale_b_nvtx_5_10 = dbe->book1D("mpTScale_b_nvtx_5_10", "pTScale_b_nvtx_5_10_0<|eta|<1.3_200_300",100, 0, 2);
    mpTScale_c_nvtx_5_10 = dbe->book1D("mpTScale_c_nvtx_5_10", "pTScale_c_nvtx_5_10_0<|eta|<1.3_600_900",100, 0, 2);
    mpTScale_a_nvtx_10_15 = dbe->book1D("mpTScale_a_nvtx_10_15", "pTScale_a_nvtx_10_15_0<|eta|<1.3_60_120",100, 0, 2);
    mpTScale_b_nvtx_10_15 = dbe->book1D("mpTScale_b_nvtx_10_15", "pTScale_b_nvtx_10_15_0<|eta|<1.3_200_300",100, 0, 2);
    mpTScale_c_nvtx_10_15 = dbe->book1D("mpTScale_c_nvtx_10_15", "pTScale_c_nvtx_10_15_0<|eta|<1.3_600_900",100, 0, 2);
    mpTScale_a_nvtx_15_20 = dbe->book1D("mpTScale_a_nvtx_15_20", "pTScale_a_nvtx_15_20_0<|eta|<1.3_60_120",100, 0, 2);
    mpTScale_b_nvtx_15_20 = dbe->book1D("mpTScale_b_nvtx_15_20", "pTScale_b_nvtx_15_20_0<|eta|<1.3_200_300",100, 0, 2);
    mpTScale_c_nvtx_15_20 = dbe->book1D("mpTScale_c_nvtx_15_20", "pTScale_c_nvtx_15_20_0<|eta|<1.3_600_900",100, 0, 2);
    mpTScale_a_nvtx_20_30 = dbe->book1D("mpTScale_a_nvtx_20_30", "pTScale_a_nvtx_20_30_0<|eta|<1.3_60_120",100, 0, 2);
    mpTScale_b_nvtx_20_30 = dbe->book1D("mpTScale_b_nvtx_20_30", "pTScale_b_nvtx_20_30_0<|eta|<1.3_200_300",100, 0, 2);
    mpTScale_c_nvtx_20_30 = dbe->book1D("mpTScale_c_nvtx_20_30", "pTScale_c_nvtx_20_30_0<|eta|<1.3_600_900",100, 0, 2);
    mpTScale_a_nvtx_30_inf = dbe->book1D("mpTScale_a_nvtx_30_inf", "pTScale_a_nvtx_30_inf_0<|eta|<1.3_60_120",100, 0, 2);
    mpTScale_b_nvtx_30_inf = dbe->book1D("mpTScale_b_nvtx_30_inf", "pTScale_b_nvtx_30_inf_0<|eta|<1.3_200_300",100, 0, 2);
    mpTScale_c_nvtx_30_inf = dbe->book1D("mpTScale_c_nvtx_30_inf", "pTScale_c_nvtx_30_inf_0<|eta|<1.3_600_900",100, 0, 2);
    mpTScale_a = dbe->book1D("mpTScale_a", "pTScale_a_60_120",100, 0, 2);
    mpTScale_b = dbe->book1D("mpTScale_b", "pTScale_b_200_300",100, 0, 2);
    mpTScale_c = dbe->book1D("mpTScale_c", "pTScale_c_600_900",100, 0, 2);
    //
    double log10PtMin = 0.5; //=3.1622766
    double log10PtMax = 3.75; //=5623.41325
    int log10PtBins = 26; 
    //double etaMin = -6.;
    //double etaMax = 6.;
    //int etaBins = 50;
    double etaRange[91] = {-6.0,-5.8,-5.6,-5.4,-5.2,-5.0,-4.8,-4.6,-4.4,-4.2,-4.0,-3.8,-3.6,-3.4,-3.2,-3.0,-2.9,-2.8,-2.7,-2.6,-2.5,-2.4,-2.3,-2.2,-2.1,-2.0,-1.9,-1.8,-1.7,-1.6,-1.5,-1.4,-1.3,-1.2,-1.1,-1.0,-0.9,-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.2,3.4,3.6,3.8,4.0,4.2,4.4,4.6,4.8,5.0,5.2,5.4,5.6,5.8,6.0};

    //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, 30, 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.5",
                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mpTScaleE_d = dbe->bookProfile("pTScaleE_d", "pTScale_d_1.5<|eta|<3.0",
                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mpTScaleF_d = dbe->bookProfile("pTScaleF_d", "pTScale_d_3.0<|eta|<6.0",
                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mpTScalePhiB_d = dbe->bookProfile("pTScalePhiB_d", "pTScalePhi_d_0<|eta|<1.5",
                                   70, -3.5, 3.5, 0, 2, " ");
    mpTScalePhiE_d = dbe->bookProfile("pTScalePhiE_d", "pTScalePhi_d_1.5<|eta|<3.0",
                                   70, -3.5, 3.5, 0, 2, " ");
    mpTScalePhiF_d = dbe->bookProfile("pTScalePhiF_d", "pTScalePhi_d_3.0<|eta|<6.0",
                      70, -3.5, 3.5, 0, 2, " ");
    /*
    mpTScale_30_200_s    = dbe->bookProfile("pTScale_30_200_s", "pTScale_s_30<pT<200",
                      etaBins, etaMin, etaMax, 0., 2., "s");
    mpTScale_200_600_s   = dbe->bookProfile("pTScale_200_600_s", "pTScale_s_200<pT<600",
                      etaBins, etaMin, etaMax, 0., 2., "s");
    mpTScale_600_1500_s   = dbe->bookProfile("pTScale_600_1500_s", "pTScale_s_600<pT<1500",
                      etaBins, etaMin, etaMax, 0., 2., "s");
    mpTScale_1500_3500_s = dbe->bookProfile("pTScale_1500_3500_s", "pTScale_s_1500<pt<3500",
                                          etaBins, etaMin, etaMax, 0., 2., "s");
    */
    mpTScale_30_200_d    = dbe->bookProfile("pTScale_30_200_d", "pTScale_d_30<pT<200",
                      90,etaRange, 0., 2., " ");
    mpTScale_200_600_d   = dbe->bookProfile("pTScale_200_600_d", "pTScale_d_200<pT<600",
                      90,etaRange, 0., 2., " ");
    mpTScale_600_1500_d   = dbe->bookProfile("pTScale_600_1500_d", "pTScale_d_600<pT<1500",
                      90,etaRange, 0., 2., " ");
    mpTScale_1500_3500_d = dbe->bookProfile("pTScale_1500_3500_d", "pTScale_d_1500<pt<3500",
                                          90,etaRange, 0., 2., " ");
    
    mpTScale1DB_30_200 = dbe->book1D("pTScale1DB_30_200", "pTScale_distribution_for_0<|eta|<1.5_30_200",
                   100, 0, 2);
    mpTScale1DE_30_200 = dbe->book1D("pTScale1DE_30_200", "pTScale_distribution_for_1.5<|eta|<3.0_30_200",
                   50, 0, 2);
    mpTScale1DF_30_200 = dbe->book1D("pTScale1DF_30_200", "pTScale_distribution_for_3.0<|eta|<6.0_30_200",
                   50, 0, 2);

    mpTScale1DB_200_600 = dbe->book1D("pTScale1DB_200_600", "pTScale_distribution_for_0<|eta|<1.5_200_600",
                   100, 0, 2);
    mpTScale1DE_200_600 = dbe->book1D("pTScale1DE_200_600", "pTScale_distribution_for_1.5<|eta|<3.0_200_600",
                   50, 0, 2);
    mpTScale1DF_200_600 = dbe->book1D("pTScale1DF_200_600", "pTScale_distribution_for_3.0<|eta|<6.0_200_600",
                   50, 0, 2);

    mpTScale1DB_600_1500 = dbe->book1D("pTScale1DB_600_1500", "pTScale_distribution_for_0<|eta|<1.5_600_1500",
                   100, 0, 2);
    mpTScale1DE_600_1500 = dbe->book1D("pTScale1DE_600_1500", "pTScale_distribution_for_1.5<|eta|<3.0_600_1500",
                   50, 0, 2);
    mpTScale1DF_600_1500 = dbe->book1D("pTScale1DF_600_1500", "pTScale_distribution_for_3.0<|eta|<6.0_600_1500",
                   50, 0, 2);

    mpTScale1DB_1500_3500 = dbe->book1D("pTScale1DB_1500_3500", "pTScale_distribution_for_0<|eta|<1.5_1500_3500",
                   100, 0, 2);
    mpTScale1DE_1500_3500 = dbe->book1D("pTScale1DE_1500_3500", "pTScale_distribution_for_1.5<|eta|<3.0_1500_3500",
                   50, 0, 2);
    mpTScale1DF_1500_3500 = dbe->book1D("pTScale1DF_1500_3500", "pTScale_distribution_for_3.0<|eta|<6.0_1500_3500",
                   50, 0, 2);
    /*
    mpTScale1D_30_200    = dbe->book1D("pTScale1D_30_200", "pTScale_distribution_for_30<pT<200",
                        100, 0, 2);
    mpTScale1D_200_600    = dbe->book1D("pTScale1D_200_600", "pTScale_distribution_for_200<pT<600",
                        100, 0, 2);
    mpTScale1D_600_1500    = dbe->book1D("pTScale1D_600_1500", "pTScale_distribution_for_600<pT<1500",
                        100, 0, 2);
    mpTScale1D_1500_3500 = dbe->book1D("pTScale1D_1500_3500", "pTScale_distribution_for_1500<pt<3500",
                        100, 0, 2);
    */

  mpTScale_nvtx_0_5  = dbe->bookProfile("pTScale_nvtx_0_5", "pTScale_nvtx_0_5_0<|eta|<1.3",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
   mpTScale_nvtx_5_10  = dbe->bookProfile("pTScale_nvtx_5_10", "pTScale_nvtx_5_10_0<|eta|<1.3",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
   mpTScale_nvtx_10_15  = dbe->bookProfile("pTScale_nvtx_10_15", "pTScale_nvtx_10_15_0<|eta|<1.3",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
   mpTScale_nvtx_15_20  = dbe->bookProfile("pTScale_nvtx_15_20", "pTScale_nvtx_15_20_0<|eta|<1.3",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
   mpTScale_nvtx_20_30  = dbe->bookProfile("pTScale_nvtx_20_30", "pTScale_nvtx_20_30_0<|eta|<1.3",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
   mpTScale_nvtx_30_inf  = dbe->bookProfile("pTScale_nvtx_30_inf", "pTScale_nvtx_30_inf_0<|eta|<1.3",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
   mpTScale_pT = dbe->bookProfile("pTScale_pT", "pTScale_vs_pT",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mpTRatio = dbe->bookProfile("pTRatio", "pTRatio",
                                log10PtBins, log10PtMin, log10PtMax, 100, 0.,5., " ");
    mpTRatioB_d = dbe->bookProfile("pTRatioB_d", "pTRatio_d_0<|eta|<1.5",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 5, " ");
    mpTRatioE_d = dbe->bookProfile("pTRatioE_d", "pTRatio_d_1.5<|eta|<3.0",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 5, " ");
    mpTRatioF_d = dbe->bookProfile("pTRatioF_d", "pTRatio_d_3.0<|eta|<6.0",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 5, " ");
    mpTRatio_30_200_d    = dbe->bookProfile("pTRatio_30_200_d", "pTRatio_d_30<pT<200",
                                          90,etaRange, 0., 5., " ");
    mpTRatio_200_600_d   = dbe->bookProfile("pTRatio_200_600_d", "pTRatio_d_200<pT<600",
                                          90,etaRange, 0., 5., " ");
    mpTRatio_600_1500_d   = dbe->bookProfile("pTRatio_600_1500_d", "pTRatio_d_600<pT<1500",
                                          90,etaRange, 0., 5., " ");    mpTRatio_1500_3500_d = dbe->bookProfile("pTRatio_1500_3500_d", "pTRatio_d_1500<pt<3500",
                                          90,etaRange, 0., 5., " ");
    mpTResponse = dbe->bookProfile("pTResponse", "pTResponse",
                log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
    mpTResponseB_d = dbe->bookProfile("pTResponseB_d", "pTResponse_d_0<|eta|<1.5",
                                   log10PtBins, log10PtMin, log10PtMax, 0.8, 1.2, " ");
    mpTResponseE_d = dbe->bookProfile("pTResponseE_d", "pTResponse_d_1.5<|eta|<3.0",
                                   log10PtBins, log10PtMin, log10PtMax, 0.8, 1.2, " ");
    mpTResponseF_d = dbe->bookProfile("pTResponseF_d", "pTResponse_d_3.0<|eta|<6.0",
                                   log10PtBins, log10PtMin, log10PtMax, 0.8, 1.2, " ");
    mpTResponse_30_200_d    = dbe->bookProfile("pTResponse_30_200_d", "pTResponse_d_30<pT<200",
                                          90,etaRange, 0.8, 1.2, " ");
    mpTResponse_200_600_d   = dbe->bookProfile("pTResponse_200_600_d", "pTResponse_d_200<pT<600",
                                          90,etaRange, 0.8, 1.2, " ");
    mpTResponse_600_1500_d   = dbe->bookProfile("pTResponse_600_1500_d", "pTResponse_d_600<pT<1500",
                                          90,etaRange, 0.8, 1.2, " ");
    mpTResponse_1500_3500_d = dbe->bookProfile("pTResponse_1500_3500_d", "pTResponse_d_1500<pt<3500",
                           90,etaRange, 0.8, 1.2, " ");
    mpTResponse_30_d = dbe->bookProfile("pTResponse_30_d", "pTResponse_d_pt>30",
                           90,etaRange, 0.8, 1.2, " ");
    mpTResponse_nvtx_0_5 = dbe->bookProfile("pTResponse_nvtx_0_5", "pTResponse_nvtx_0_5", log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
   mpTResponse_nvtx_5_10 = dbe->bookProfile("pTResponse_nvtx_5_10", "pTResponse_nvtx_5_10", log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
   mpTResponse_nvtx_10_15 = dbe->bookProfile("pTResponse_nvtx_10_15", "pTResponse_nvtx_10_15", log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
   mpTResponse_nvtx_15_20 = dbe->bookProfile("pTResponse_nvtx_15_20", "pTResponse_nvtx_15_20", log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
   mpTResponse_nvtx_20_30 = dbe->bookProfile("pTResponse_nvtx_20_30", "pTResponse_nvtx_20_30", log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
   mpTResponse_nvtx_30_inf = dbe->bookProfile("pTResponse_nvtx_30_inf", "pTResponse_nvtx_30_inf", log10PtBins, log10PtMin, log10PtMax, 100, 0.8,1.2, " ");
  }

  if (mOutputFile.empty ()) {
    LogInfo("OutputInfo") << " PFJet histograms will NOT be saved";
  } 
  else {
    LogInfo("OutputInfo") << " PFJethistograms will be saved to file:" << mOutputFile;
  }
}

PFJetTester::~PFJetTester

(

)

Definition at line 534 of file PFJetTester.cc.

{
}

Member Function Documentation

void PFJetTester::analyze ( const edm::Event &  mEvent, const edm::EventSetup &  mSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 546 of file PFJetTester.cc.

References HiRecoJets_cff::caloTowers, JetCorrector::correction(), deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, alignCSCRings::e, reco::LeafCandidate::eta(), MonitorElement::Fill(), fillMatchHists(), edm::Event::getByLabel(), JetCorrector::getJetCorrector(), edm::Event::getManyByType(), metFilters_cff::goodVertices, i, edm::EventBase::isRealData(), edm::HandleBase::isValid(), j, metsig::jet, JetCorrectionService, edm::InputTag::label(), scaleCards::mass, mChargedEmEnergy_3000, mChargedEmEnergy_80, mChargedEmEnergyFraction_B, mChargedEmEnergyFraction_E, mChargedEmEnergyFraction_F, mChargedHadronEnergy_3000, mChargedHadronEnergy_80, mChargedHadronEnergyFraction_B, mChargedHadronEnergyFraction_E, mChargedHadronEnergyFraction_F, mConstituents, mConstituents_80, mCorrJetEta, mCorrJetPhi, mCorrJetPt, mCorrJetPt_80, mE, mE_80, mEmEnergyInHF, mEmEnergyInHF_3000, mEmEnergyInHF_80, mEmTiming, mEta, mEtaFineBin, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_80, mHadEnergyInHF, mHadEnergyInHF_3000, mHadEnergyInHF_80, mHadTiming, mInputCollection, mInputGenCollection, mjetArea, mMass, mMass_80, mMatchGenPtThreshold, mMjj, mMjj_3000, mNeutralEmEnergy_3000, mNeutralEmEnergy_80, mNeutralEmEnergyFraction_B, mNeutralEmEnergyFraction_E, mNeutralEmEnergyFraction_F, mNeutralHadronEnergy_3000, mNeutralHadronEnergy_80, mNeutralHadronEnergyFraction_B, mNeutralHadronEnergyFraction_E, mNeutralHadronEnergyFraction_F, mNJets1, mNJets2, mNJetsEtaC, mNJetsEtaF, mNJetsEtaF_30, mP, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPthat_3000, mPthat_80, mpTRatio, mpTRatio_1500_3500_d, mpTRatio_200_600_d, mpTRatio_30_200_d, mpTRatio_600_1500_d, mpTRatioB_d, mpTRatioE_d, mpTRatioF_d, mpTResponse, mpTResponse_1500_3500_d, mpTResponse_200_600_d, mpTResponse_30_200_d, mpTResponse_30_d, mpTResponse_600_1500_d, mpTResponse_nvtx_0_5, mpTResponse_nvtx_10_15, mpTResponse_nvtx_15_20, mpTResponse_nvtx_20_30, mpTResponse_nvtx_30_inf, mpTResponse_nvtx_5_10, mpTResponseB_d, mpTResponseE_d, mpTResponseF_d, mRThreshold, mTurnOnEverything, njet, numberofevents, nvtx_0_30, nvtx_0_60, pfJets_cff::pfJets, reco::LeafCandidate::phi(), position, reco::LeafCandidate::pt(), rho, pileupReCalc_HLTpaths::scale, reco::Jet::scaleEnergy(), and detailsBasic3DVector::z.

{
  double countsfornumberofevents = 1;
  numberofevents->Fill(countsfornumberofevents);
 
  //get primary vertices
  edm::Handle<vector<reco::Vertex> >pvHandle;
  try {
    mEvent.getByLabel( "offlinePrimaryVertices", pvHandle );
  } catch ( cms::Exception & e ) {
    //cout <<prefix<<"error: " << e.what() << endl;
  }
  vector<reco::Vertex> goodVertices;
  for (unsigned i = 0; i < pvHandle->size(); i++) {
    if ( (*pvHandle)[i].ndof() > 4 &&
       ( fabs((*pvHandle)[i].z()) <= 24. ) &&
       ( fabs((*pvHandle)[i].position().rho()) <= 2.0 ) )
       goodVertices.push_back((*pvHandle)[i]);
  }
  
  nvtx_0_30->Fill(goodVertices.size());
  nvtx_0_60->Fill(goodVertices.size());
  
  // *********************************
  // *** Get pThat
  // *********************************
if (!mEvent.isRealData()){
  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<PFJetCollection> pfJets;
  mEvent.getByLabel(mInputCollection, pfJets);
  if (!pfJets.isValid()) return;
  PFJetCollection::const_iterator jet = pfJets->begin ();
  int jetIndex = 0;
  int nJet = 0;
  int nJetF = 0;
  int nJetC = 0;
  int nJetF_30 =0;
  for (; jet != pfJets->end (); jet++, jetIndex++) {

    if (jet->pt() > 10.) {
      if (fabs(jet->eta()) > 1.5) 
    nJetF++;
      else 
    nJetC++;      
    }
    if (jet->pt() > 30.) nJetF_30++;
    if (jet->pt() > 10.) {
      if (mEta) mEta->Fill (jet->eta());
      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 (mjetArea) mjetArea->Fill(jet->jetArea());
    if (mPhi) mPhi->Fill (jet->phi());
    if (mE) mE->Fill (jet->energy());
    if (mE_80) mE_80->Fill (jet->energy());
    if (mP) mP->Fill (jet->p());
    if (mP_80) mP_80->Fill (jet->p());
    if (mPt) mPt->Fill (jet->pt());
    if (mPt_80) mPt_80->Fill (jet->pt());
    if (mMass) mMass->Fill (jet->mass());
    if (mMass_80) mMass_80->Fill (jet->mass());
    if (mConstituents) mConstituents->Fill (jet->nConstituents());
    if (mConstituents_80) mConstituents_80->Fill (jet->nConstituents());
    if (jet == pfJets->begin ()) { // first jet
      if (mEtaFirst) mEtaFirst->Fill (jet->eta());
      if (mPhiFirst) mPhiFirst->Fill (jet->phi());
      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->HFHadronEnergy());
    if (mHadEnergyInHF_80) mHadEnergyInHF_80->Fill (jet->HFHadronEnergy());
    if (mHadEnergyInHF_3000) mHadEnergyInHF_3000->Fill (jet->HFHadronEnergy());
    //    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->HFEMEnergy());
    if (mEmEnergyInHF_80) mEmEnergyInHF_80->Fill (jet->HFEMEnergy());
    if (mEmEnergyInHF_3000) mEmEnergyInHF_3000->Fill (jet->HFEMEnergy());
    //    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 (mChargedEmEnergy_80) mChargedEmEnergy_80->Fill (jet->chargedEmEnergy());
    if (mChargedHadronEnergy_80) mChargedHadronEnergy_80->Fill (jet->chargedHadronEnergy());
    if (mNeutralEmEnergy_80) mNeutralEmEnergy_80->Fill (jet->neutralEmEnergy());
    if (mNeutralHadronEnergy_80) mNeutralHadronEnergy_80->Fill (jet->neutralHadronEnergy());

    if (mChargedEmEnergy_3000) mChargedEmEnergy_3000->Fill (jet->chargedEmEnergy());
    if (mChargedHadronEnergy_3000) mChargedHadronEnergy_3000->Fill (jet->chargedHadronEnergy());
    if (mNeutralEmEnergy_3000) mNeutralEmEnergy_3000->Fill (jet->neutralEmEnergy());
    if (mNeutralHadronEnergy_3000) mNeutralHadronEnergy_3000->Fill (jet->neutralHadronEnergy());

    if (fabs(jet->eta())<1.5) mChargedEmEnergyFraction_B->Fill (jet->chargedEmEnergyFraction());
    if (fabs(jet->eta())>1.5 && fabs(jet->eta())<3.0) mChargedEmEnergyFraction_E->Fill (jet->chargedEmEnergyFraction());
    if (fabs(jet->eta())>3.0 && fabs(jet->eta())<6.0) mChargedEmEnergyFraction_F->Fill (jet->chargedEmEnergyFraction());
    if (fabs(jet->eta())<1.5) mChargedHadronEnergyFraction_B->Fill (jet->chargedHadronEnergyFraction());
    if (fabs(jet->eta())>1.5 && fabs(jet->eta())<3.0) mChargedHadronEnergyFraction_E->Fill (jet->chargedHadronEnergyFraction());
    if (fabs(jet->eta())>3.0 && fabs(jet->eta())<6.0) mChargedHadronEnergyFraction_F->Fill (jet->chargedHadronEnergyFraction());
    if (fabs(jet->eta())<1.5) mNeutralEmEnergyFraction_B->Fill (jet->neutralEmEnergyFraction());
    if (fabs(jet->eta())>1.5 && fabs(jet->eta())<3.0) mNeutralEmEnergyFraction_E->Fill (jet->neutralEmEnergyFraction());
    if (fabs(jet->eta())>3.0 && fabs(jet->eta())<6.0) mNeutralEmEnergyFraction_F->Fill (jet->neutralEmEnergyFraction());
    if (fabs(jet->eta())<1.5) mNeutralHadronEnergyFraction_B->Fill (jet->neutralHadronEnergyFraction());
    if (fabs(jet->eta())>1.5 && fabs(jet->eta())<3.0) mNeutralHadronEnergyFraction_E->Fill (jet->neutralHadronEnergyFraction());
    if (fabs(jet->eta())>3.0 && fabs(jet->eta())<6.0) mNeutralHadronEnergyFraction_F->Fill (jet->neutralHadronEnergyFraction());

    //    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 (mNJetsEtaF_30) mNJetsEtaF_30->Fill( nJetF_30 );  

  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 ( PFJetCollection::const_iterator cal = pfJets->begin(); cal != pfJets->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 ( PFJetCollection::const_iterator cal = pfJets->begin(); cal != pfJets->end(); ++ cal ) {
      if ( cal->pt() > ptStep ) njet++;
    }
    mNJets2->Fill( ptStep, njet );
  }

 // Correction jets
  const JetCorrector* corrector = JetCorrector::getJetCorrector (JetCorrectionService,mSetup);

  for (PFJetCollection::const_iterator jet = pfJets->begin(); jet !=pfJets ->end(); jet++) 
  {
 
      PFJet  correctedJet = *jet;
      //double scale = corrector->correction(jet->p4());
      double scale = corrector->correction(*jet,mEvent,mSetup); 
      correctedJet.scaleEnergy(scale); 
      if(correctedJet.pt()>30){
      mCorrJetPt->Fill(correctedJet.pt());
      mCorrJetPt_80->Fill(correctedJet.pt());
      if(correctedJet.pt()>10) mCorrJetEta->Fill(correctedJet.eta());
      mCorrJetPhi->Fill(correctedJet.phi());
      mpTRatio->Fill(log10(jet->pt()),correctedJet.pt()/jet->pt());

      if (fabs(jet->eta())<1.5) {
           mpTRatioB_d->Fill(log10(jet->pt()), correctedJet.pt()/jet->pt());
      }

     if (fabs(jet->eta())>1.5 && fabs(jet->eta())<3.0) {
        mpTRatioE_d->Fill (log10(jet->pt()), correctedJet.pt()/jet->pt());
     }
     if (fabs(jet->eta())>3.0 && fabs(jet->eta())<6.0) {
        mpTRatioF_d->Fill (log10(jet->pt()), correctedJet.pt()/jet->pt());
    }
     if (jet->pt()>30.0 && jet->pt()<200.0) {
    mpTRatio_30_200_d->Fill (jet->eta(),correctedJet.pt()/jet->pt());
  }
   if (jet->pt()>200.0 && jet->pt()<600.0) {
    mpTRatio_200_600_d->Fill (jet->eta(),correctedJet.pt()/jet->pt());
  }
if (jet->pt()>600.0 && jet->pt()<1500.0) {
    mpTRatio_600_1500_d->Fill (jet->eta(),correctedJet.pt()/jet->pt());
  }
if (jet->pt()>1500.0 && jet->pt()<3500.0) {
    mpTRatio_1500_3500_d->Fill (jet->eta(),correctedJet.pt()/jet->pt());
  }
      }



  }

if (!mEvent.isRealData()){
  // 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 (gjet == genJets->begin ()) { // first jet
      if (mGenEtaFirst) mGenEtaFirst->Fill (gjet->eta());
      if (mGenPhiFirst) mGenPhiFirst->Fill (gjet->phi());
    }
  }


  // now match PFJets 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*> > pfJetConstituents (pfJets->size());
    /*
    if (mRThreshold > 0) { 
    }
    else {
      for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
    genJetConstituents [iGenJet] = jetMatching.getGenParticles ((*genJets) [iGenJet]);
      }
      
      for (unsigned iPFJet = 0; iPFJet < pfJets->size(); ++iPFJet) {
    pfJetConstituents [iPFJet] = jetMatching.getGenParticles ((*pfJets) [iPFJet], 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()) > 6.) continue; // out of detector 
      if (genJetPt < mMatchGenPtThreshold) continue; // no low momentum 
      //double logPtGen = log10 (genJetPt);
      //mAllGenJetsPt->Fill (logPtGen);
      //mAllGenJetsEta->Fill (logPtGen, genJet.eta());
      if (pfJets->size() <= 0) continue; // no PFJets - nothing to match
      if (mRThreshold > 0) {
    unsigned iPFJetBest = 0;
    double deltaRBest = 999.;
    for (unsigned iPFJet = 0; iPFJet < pfJets->size(); ++iPFJet) {
      double dR = deltaR (genJet.eta(), genJet.phi(), (*pfJets) [iPFJet].eta(), (*pfJets) [iPFJet].phi());
      if (deltaRBest < mRThreshold && dR < mRThreshold && genJet.pt() > 5.) {
        /*
        std::cout << "Yet another matched jet for GenJet pt=" << genJet.pt()
              << " previous PFJet pt/dr: " << (*pfJets) [iPFJetBest].pt() << '/' << deltaRBest
              << " new PFJet pt/dr: " << (*pfJets) [iPFJet].pt() << '/' << dR
              << std::endl;
        */
      }
      if (dR < deltaRBest) {
        iPFJetBest = iPFJet;
        deltaRBest = dR;
      }
    }
    if (mTurnOnEverything.compare("yes")==0) {
      //mRMatch->Fill (logPtGen, genJet.eta(), deltaRBest);
    }
    if (deltaRBest < mRThreshold) { // Matched
      fillMatchHists (genJet, (*pfJets) [iPFJetBest],goodVertices);
    }

    double CorrdeltaRBest = 999.;
    double CorrJetPtBest = 0;
    for (PFJetCollection::const_iterator jet = pfJets->begin(); jet !=pfJets ->end(); jet++) {
      PFJet  correctedJet = *jet;
      //double scale = corrector->correction(jet->p4());
          double scale = corrector->correction(*jet,mEvent,mSetup); 
      correctedJet.scaleEnergy(scale);
      double CorrJetPt = correctedJet.pt();
      if(CorrJetPt>30){
      double CorrdR = deltaR (genJet.eta(), genJet.phi(), correctedJet.eta(), correctedJet.phi());
      if (CorrdR < CorrdeltaRBest) {
        CorrdeltaRBest = CorrdR;
        CorrJetPtBest = CorrJetPt;
      }
    }
    }
    if (deltaRBest < mRThreshold) { // Matched
      mpTResponse->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());
      
      if (fabs(genJet.eta())<1.5) {
        mpTResponseB_d->Fill(log10(genJet.pt()), CorrJetPtBest/genJet.pt());
      } 
      
      if (fabs(genJet.eta())>1.5 && fabs(genJet.eta())<3.0) {
        mpTResponseE_d->Fill (log10(genJet.pt()), CorrJetPtBest/genJet.pt());   
      }
      if (fabs(genJet.eta())>3.0 && fabs(genJet.eta())<6.0) {
        mpTResponseF_d->Fill (log10(genJet.pt()), CorrJetPtBest/genJet.pt());
      }
      if (genJet.pt()>30.0 && genJet.pt()<200.0) {
        mpTResponse_30_200_d->Fill (genJet.eta(),CorrJetPtBest/genJet.pt());
      }
      if (genJet.pt()>200.0 && genJet.pt()<600.0) {
        mpTResponse_200_600_d->Fill (genJet.eta(),CorrJetPtBest/genJet.pt());
      }
      if (genJet.pt()>600.0 && genJet.pt()<1500.0) {
        mpTResponse_600_1500_d->Fill (genJet.eta(),CorrJetPtBest/genJet.pt());
      }
      if (genJet.pt()>1500.0 && genJet.pt()<3500.0) {
        mpTResponse_1500_3500_d->Fill (genJet.eta(),CorrJetPtBest/genJet.pt());
      }
      if (genJet.pt()>30.0) {
        mpTResponse_30_d->Fill (genJet.eta(),CorrJetPtBest/genJet.pt());
      }
      if(goodVertices.size()<=5) mpTResponse_nvtx_0_5->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());
      if(goodVertices.size()>5 && goodVertices.size()<=10) mpTResponse_nvtx_5_10->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());
      if(goodVertices.size()>10 && goodVertices.size()<=15) mpTResponse_nvtx_10_15->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());   
      if(goodVertices.size()>15 && goodVertices.size()<=20) mpTResponse_nvtx_15_20->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());
      if(goodVertices.size()>20 && goodVertices.size()<=30) mpTResponse_nvtx_20_30->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());
      if(goodVertices.size()>30) mpTResponse_nvtx_30_inf->Fill(log10(genJet.pt()),CorrJetPtBest/genJet.pt());
    }

      }
      /*
      else {
    unsigned iPFJetBest = 0;
    double energyFractionBest = 0.;
    for (unsigned iPFJet = 0; iPFJet < pfJets->size(); ++iPFJet) {
      double energyFraction = jetMatching.overlapEnergyFraction (genJetConstituents [iGenJet], 
                                     pfJetConstituents [iPFJet]);
      if (energyFraction > energyFractionBest) {
        iPFJetBest = iPFJet;
        energyFractionBest = energyFraction;
      }
    }
    if (mTurnOnEverything.compare("yes")==0) {
      mGenJetMatchEnergyFraction->Fill (logPtGen, genJet.eta(), energyFractionBest);
    }
    if (energyFractionBest > mGenEnergyFractionThreshold) { // good enough
      double reverseEnergyFraction = jetMatching.overlapEnergyFraction (pfJetConstituents [iPFJetBest], 
                                        genJetConstituents [iGenJet]);
      if (mTurnOnEverything.compare("yes")==0) {
        mReverseMatchEnergyFraction->Fill (logPtGen, genJet.eta(), reverseEnergyFraction);
      }
      if (reverseEnergyFraction > mReverseEnergyFractionThreshold) { // Matched
        fillMatchHists (genJet, (*pfJets) [iPFJetBest]);
      }
    }
      }
      */
    }
  }
}

}

void PFJetTester::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 538 of file PFJetTester.cc.

                          {
}

void PFJetTester::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 541 of file PFJetTester.cc.

References mOutputFile.

                         {
 if (!mOutputFile.empty() && &*edm::Service<DQMStore>()) edm::Service<DQMStore>()->save (mOutputFile);
}

void PFJetTester::fillMatchHists ( const reco::GenJet &  fGenJet, const reco::PFJet &  fPFJet, std::vector< reco::Vertex >  goodVertices  )

private

Definition at line 1109 of file PFJetTester.cc.

References reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), MonitorElement::Fill(), mDelEta, mDelPhi, mDelPt, mDeltaEta, mDeltaPhi, mEScale_pt10, mEScaleFineBin, mpTScale1DB_1500_3500, mpTScale1DB_200_600, mpTScale1DB_30_200, mpTScale1DB_600_1500, mpTScale1DE_1500_3500, mpTScale1DE_200_600, mpTScale1DE_30_200, mpTScale1DE_600_1500, mpTScale1DF_1500_3500, mpTScale1DF_200_600, mpTScale1DF_30_200, mpTScale1DF_600_1500, mpTScale_1500_3500_d, mpTScale_200_600_d, mpTScale_30_200_d, mpTScale_600_1500_d, mpTScale_a, mpTScale_a_nvtx_0_5, mpTScale_a_nvtx_10_15, mpTScale_a_nvtx_15_20, mpTScale_a_nvtx_20_30, mpTScale_a_nvtx_30_inf, mpTScale_a_nvtx_5_10, mpTScale_b, mpTScale_b_nvtx_0_5, mpTScale_b_nvtx_10_15, mpTScale_b_nvtx_15_20, mpTScale_b_nvtx_20_30, mpTScale_b_nvtx_30_inf, mpTScale_b_nvtx_5_10, mpTScale_c, mpTScale_c_nvtx_0_5, mpTScale_c_nvtx_10_15, mpTScale_c_nvtx_15_20, mpTScale_c_nvtx_20_30, mpTScale_c_nvtx_30_inf, mpTScale_c_nvtx_5_10, mpTScale_nvtx_0_5, mpTScale_nvtx_10_15, mpTScale_nvtx_15_20, mpTScale_nvtx_20_30, mpTScale_nvtx_30_inf, mpTScale_nvtx_5_10, mpTScale_pT, mpTScaleB_d, mpTScaleE_d, mpTScaleF_d, mpTScalePhiB_d, mpTScalePhiE_d, mpTScalePhiF_d, mTurnOnEverything, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().

Referenced by analyze().

                                                                                                                       {
  double logPtGen = log10 (fGenJet.pt());
  double PtGen = fGenJet.pt();
  double PtPF = fPFJet.pt();
  //mMatchedGenJetsPt->Fill (logPtGen);
  //mMatchedGenJetsEta->Fill (logPtGen, fGenJet.eta());

  double PtThreshold = 10.;

  if (mTurnOnEverything.compare("yes")==0) {
    mDeltaEta->Fill (logPtGen, fGenJet.eta(), fPFJet.eta()-fGenJet.eta());
    mDeltaPhi->Fill (logPtGen, fGenJet.eta(), fPFJet.phi()-fGenJet.phi());
    //mEScale->Fill (logPtGen, fGenJet.eta(), fPFJet.energy()/fGenJet.energy());
    //mlinEScale->Fill (fGenJet.pt(), fGenJet.eta(), fPFJet.energy()/fGenJet.energy());
    //mDeltaE->Fill (logPtGen, fGenJet.eta(), fPFJet.energy()-fGenJet.energy());

    mEScaleFineBin->Fill (logPtGen, fGenJet.eta(), fPFJet.energy()/fGenJet.energy());
  
    if (fGenJet.pt()>PtThreshold) {
      mEScale_pt10->Fill (logPtGen, fGenJet.eta(), fPFJet.energy()/fGenJet.energy());

    }

  }
  if (fPFJet.pt() > PtThreshold) {
    mDelEta->Fill (fGenJet.eta()-fPFJet.eta());
    mDelPhi->Fill (fGenJet.phi()-fPFJet.phi());
    mDelPt->Fill  ((fGenJet.pt()-fPFJet.pt())/fGenJet.pt());
  }

  if (fabs(fGenJet.eta())<1.5) {

    //mpTScaleB_s->Fill (log10(PtGen), PtPF/PtGen);
    mpTScaleB_d->Fill (log10(PtGen), PtPF/PtGen);
    mpTScalePhiB_d->Fill (fGenJet.phi(), PtPF/PtGen);
    
    if (PtGen>30.0 && PtGen<200.0) {
      mpTScale1DB_30_200->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>200.0 && PtGen<600.0) {
      mpTScale1DB_200_600->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>600.0 && PtGen<1500.0) {
      mpTScale1DB_600_1500->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>1500.0 && PtGen<3500.0) {
      mpTScale1DB_1500_3500->Fill (fPFJet.pt()/fGenJet.pt());
    }
    
  }

  if (fabs(fGenJet.eta())>1.5 && fabs(fGenJet.eta())<3.0) {

    //mpTScaleE_s->Fill (log10(PtGen), PtPF/PtGen);
    mpTScaleE_d->Fill (log10(PtGen), PtPF/PtGen);
    mpTScalePhiE_d->Fill (fGenJet.phi(), PtPF/PtGen);
    
    if (PtGen>30.0 && PtGen<200.0) {
      mpTScale1DE_30_200->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>200.0 && PtGen<600.0) {
      mpTScale1DE_200_600->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>600.0 && PtGen<1500.0) {
      mpTScale1DE_600_1500->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>1500.0 && PtGen<3500.0) {
      mpTScale1DE_1500_3500->Fill (fPFJet.pt()/fGenJet.pt());
    }
    
  }

  if (fabs(fGenJet.eta())>3.0 && fabs(fGenJet.eta())<6.0) {

    //mpTScaleF_s->Fill (log10(PtGen), PtPF/PtGen);
    mpTScaleF_d->Fill (log10(PtGen), PtPF/PtGen);
    mpTScalePhiF_d->Fill (fGenJet.phi(), PtPF/PtGen);
    
    if (PtGen>30.0 && PtGen<200.0) {
      mpTScale1DF_30_200->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>200.0 && PtGen<600.0) {
      mpTScale1DF_200_600->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>600.0 && PtGen<1500.0) {
      mpTScale1DF_600_1500->Fill (fPFJet.pt()/fGenJet.pt());
    }
    if (PtGen>1500.0 && PtGen<3500.0) {
      mpTScale1DF_1500_3500->Fill (fPFJet.pt()/fGenJet.pt());
    }
    
  }

  if (fGenJet.pt()>30.0 && fGenJet.pt()<200.0) {
    //mpTScale_30_200_s->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    mpTScale_30_200_d->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    //mpTScale1D_30_200->Fill (fPFJet.pt()/fGenJet.pt());
  }

  if (fGenJet.pt()>200.0 && fGenJet.pt()<600.0) {
    //mpTScale_200_600_s->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    mpTScale_200_600_d->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    //mpTScale1D_200_600->Fill (fPFJet.pt()/fGenJet.pt());
  }

  if (fGenJet.pt()>600.0 && fGenJet.pt()<1500.0) {
    //mpTScale_600_1500_s->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    mpTScale_600_1500_d->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    //mpTScale1D_600_1500->Fill (fPFJet.pt()/fGenJet.pt());
  }

  if (fGenJet.pt()>1500.0 && fGenJet.pt()<3500.0) {
    //mpTScale_1500_3500_s->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    mpTScale_1500_3500_d->Fill (fGenJet.eta(),fPFJet.pt()/fGenJet.pt());
    //mpTScale1D_1500_3500->Fill (fPFJet.pt()/fGenJet.pt());
  }
  if (fabs(fGenJet.eta())<1.3) {
    if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
     if(goodVertices.size()<=5) mpTScale_a_nvtx_0_5->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
     if(goodVertices.size()<=5) mpTScale_b_nvtx_0_5->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
     if(goodVertices.size()<=5) mpTScale_c_nvtx_0_5->Fill( PtPF/PtGen);
    }

    if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
     if(goodVertices.size()>5 && goodVertices.size()<=10) mpTScale_a_nvtx_5_10->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
     if(goodVertices.size()>5 && goodVertices.size()<=10) mpTScale_b_nvtx_5_10->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
     if(goodVertices.size()>5 && goodVertices.size()<=10) mpTScale_c_nvtx_5_10->Fill( PtPF/PtGen);
    }

    if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
     if(goodVertices.size()>10 && goodVertices.size()<=15) mpTScale_a_nvtx_10_15->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
     if(goodVertices.size()>10 && goodVertices.size()<=15) mpTScale_b_nvtx_10_15->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
     if(goodVertices.size()>10 && goodVertices.size()<=15) mpTScale_c_nvtx_10_15->Fill( PtPF/PtGen);
    }

    if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
     if(goodVertices.size()>15 && goodVertices.size()<=20) mpTScale_a_nvtx_15_20->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
     if(goodVertices.size()>15 && goodVertices.size()<=20) mpTScale_b_nvtx_15_20->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
     if(goodVertices.size()>15 && goodVertices.size()<=20) mpTScale_c_nvtx_15_20->Fill( PtPF/PtGen);
    }

if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
     if(goodVertices.size()>20 && goodVertices.size()<=30) mpTScale_a_nvtx_20_30->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
     if(goodVertices.size()>20 && goodVertices.size()<=30) mpTScale_b_nvtx_20_30->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
     if(goodVertices.size()>20 && goodVertices.size()<=30) mpTScale_c_nvtx_20_30->Fill( PtPF/PtGen);
    }

if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) {
     if(goodVertices.size()>30) mpTScale_a_nvtx_30_inf->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) {
     if(goodVertices.size()>30) mpTScale_b_nvtx_30_inf->Fill( PtPF/PtGen);
    }
    if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) {
     if(goodVertices.size()>30) mpTScale_c_nvtx_30_inf->Fill( PtPF/PtGen);
    }

   if(goodVertices.size()<=5) mpTScale_nvtx_0_5->Fill(log10(PtGen),PtPF/PtGen);
    if(goodVertices.size()>5 && goodVertices.size()<=10) mpTScale_nvtx_5_10->Fill(log10(PtGen),PtPF/PtGen);
    if(goodVertices.size()>10 && goodVertices.size()<=15) mpTScale_nvtx_10_15->Fill(log10(PtGen),PtPF/PtGen);
    if(goodVertices.size()>15 && goodVertices.size()<=20) mpTScale_nvtx_15_20->Fill(log10(PtGen), PtPF/PtGen);
    if(goodVertices.size()>20 && goodVertices.size()<=30) mpTScale_nvtx_20_30->Fill(log10(PtGen), PtPF/PtGen);
    if(goodVertices.size()>30) mpTScale_nvtx_30_inf->Fill(log10(PtGen), PtPF/PtGen);
}
  if (fabs(fGenJet.eta())<1.3) {
  if(fGenJet.pt()>60.0 && fGenJet.pt()<120.0) mpTScale_a->Fill(PtPF/PtGen);
  if(fGenJet.pt()>200.0 && fGenJet.pt()<300.0) mpTScale_b->Fill(PtPF/PtGen);
  if(fGenJet.pt()>600.0 && fGenJet.pt()<900.0) mpTScale_c->Fill(PtPF/PtGen);
  }
  mpTScale_pT->Fill (log10(PtGen), PtPF/PtGen);
}

Member Data Documentation

std::string PFJetTester::inputCaloMETLabel_

private

Definition at line 43 of file PFJetTester.h.

std::string PFJetTester::inputGenMETLabel_

private

Definition at line 42 of file PFJetTester.h.

edm::InputTag PFJetTester::inputMETLabel_

private

Definition at line 40 of file PFJetTester.h.

std::string PFJetTester::JetCorrectionService

private

Definition at line 310 of file PFJetTester.h.

Referenced by analyze().

MonitorElement* PFJetTester::mAllGenJetsEta

private

Definition at line 249 of file PFJetTester.h.

MonitorElement* PFJetTester::mAllGenJetsPt

private

Definition at line 247 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMET

private

Definition at line 327 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMET_3000

private

Definition at line 328 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMETPhi

private

Definition at line 329 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMETSig

private

Definition at line 332 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMETSig_3000

private

Definition at line 333 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMEx

private

Definition at line 323 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMEx_3000

private

Definition at line 324 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMEy

private

Definition at line 325 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloMEy_3000

private

Definition at line 326 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloSumET

private

Definition at line 330 of file PFJetTester.h.

MonitorElement* PFJetTester::mCaloSumET_3000

private

Definition at line 331 of file PFJetTester.h.

MonitorElement* PFJetTester::mChargedEmEnergy_3000

private

Definition at line 167 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedEmEnergy_80

private

Definition at line 162 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedEmEnergyFraction_B

private

Definition at line 172 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedEmEnergyFraction_E

private

Definition at line 173 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedEmEnergyFraction_F

private

Definition at line 174 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedHadronEnergy_3000

private

Definition at line 168 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedHadronEnergy_80

private

Definition at line 163 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedHadronEnergyFraction_B

private

Definition at line 175 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedHadronEnergyFraction_E

private

Definition at line 176 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mChargedHadronEnergyFraction_F

private

Definition at line 177 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mConstituents

private

Definition at line 71 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mConstituents_3000

private

Definition at line 73 of file PFJetTester.h.

MonitorElement* PFJetTester::mConstituents_80

private

Definition at line 72 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mCorrJetEta

private

Definition at line 81 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mCorrJetPhi

private

Definition at line 82 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mCorrJetPt

private

Definition at line 78 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mCorrJetPt_3000

private

Definition at line 80 of file PFJetTester.h.

MonitorElement* PFJetTester::mCorrJetPt_80

private

Definition at line 79 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mDelEta

private

Definition at line 300 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mDelPhi

private

Definition at line 301 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mDelPt

private

Definition at line 302 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mDeltaE

private

Definition at line 258 of file PFJetTester.h.

MonitorElement* PFJetTester::mDeltaEta

private

Definition at line 254 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mDeltaPhi

private

Definition at line 255 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mE

private

Definition at line 59 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mE_3000

private

Definition at line 61 of file PFJetTester.h.

MonitorElement* PFJetTester::mE_80

private

Definition at line 60 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEBEne

private

Definition at line 344 of file PFJetTester.h.

MonitorElement* PFJetTester::mEBTime

private

Definition at line 345 of file PFJetTester.h.

MonitorElement* PFJetTester::mEEEne

private

Definition at line 346 of file PFJetTester.h.

MonitorElement* PFJetTester::mEETime

private

Definition at line 347 of file PFJetTester.h.

MonitorElement* PFJetTester::mEFirst

private

Definition at line 144 of file PFJetTester.h.

MonitorElement* PFJetTester::mEFirst_3000

private

Definition at line 146 of file PFJetTester.h.

MonitorElement* PFJetTester::mEFirst_80

private

Definition at line 145 of file PFJetTester.h.

MonitorElement* PFJetTester::mElectronEnergy

private

Definition at line 221 of file PFJetTester.h.

MonitorElement* PFJetTester::mElectronEnergy_3000

private

Definition at line 223 of file PFJetTester.h.

MonitorElement* PFJetTester::mElectronEnergy_80

private

Definition at line 222 of file PFJetTester.h.

MonitorElement* PFJetTester::mEmEnergyInHF

private

Definition at line 201 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEmEnergyInHF_3000

private

Definition at line 203 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEmEnergyInHF_80

private

Definition at line 202 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEmEnergyProfile

private

Definition at line 317 of file PFJetTester.h.

MonitorElement* PFJetTester::mEmTiming

private

Definition at line 75 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEScale

private

Definition at line 256 of file PFJetTester.h.

MonitorElement* PFJetTester::mEScale_pt10

private

Definition at line 260 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mEScaleFineBin

private

Definition at line 261 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mEta

private

Definition at line 49 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEtaFineBin

private

Definition at line 50 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mEtaFineBin1m

private

Definition at line 54 of file PFJetTester.h.

MonitorElement* PFJetTester::mEtaFineBin1p

private

Definition at line 51 of file PFJetTester.h.

MonitorElement* PFJetTester::mEtaFineBin2m

private

Definition at line 55 of file PFJetTester.h.

MonitorElement* PFJetTester::mEtaFineBin2p

private

Definition at line 52 of file PFJetTester.h.

MonitorElement* PFJetTester::mEtaFineBin3m

private

Definition at line 56 of file PFJetTester.h.

MonitorElement* PFJetTester::mEtaFineBin3p

private

Definition at line 53 of file PFJetTester.h.

MonitorElement* PFJetTester::mEtaFirst

private

Definition at line 142 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

std::string PFJetTester::METType_

private

Definition at line 41 of file PFJetTester.h.

double PFJetTester::mGenEnergyFractionThreshold

private

Definition at line 306 of file PFJetTester.h.

MonitorElement* PFJetTester::mGenEta

private

Definition at line 236 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mGenEtaFirst

private

Definition at line 243 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mGenJetMatchEnergyFraction

private

Definition at line 251 of file PFJetTester.h.

MonitorElement* PFJetTester::mGenPhi

private

Definition at line 237 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mGenPhiFirst

private

Definition at line 244 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mGenPt

private

Definition at line 238 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mGenPt_3000

private

Definition at line 240 of file PFJetTester.h.

MonitorElement* PFJetTester::mGenPt_80

private

Definition at line 239 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mHadEnergyInHF

private

Definition at line 189 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mHadEnergyInHF_3000

private

Definition at line 191 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mHadEnergyInHF_80

private

Definition at line 190 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mHadEnergyProfile

private

Definition at line 316 of file PFJetTester.h.

MonitorElement* PFJetTester::mHadTiming

private

Definition at line 74 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mHBEne

private

Definition at line 336 of file PFJetTester.h.

MonitorElement* PFJetTester::mHBTime

private

Definition at line 337 of file PFJetTester.h.

MonitorElement* PFJetTester::mHEEne

private

Definition at line 338 of file PFJetTester.h.

MonitorElement* PFJetTester::mHETime

private

Definition at line 339 of file PFJetTester.h.

MonitorElement* PFJetTester::mHFEne

private

Definition at line 342 of file PFJetTester.h.

MonitorElement* PFJetTester::mHFTime

private

Definition at line 343 of file PFJetTester.h.

MonitorElement* PFJetTester::mHOEne

private

Definition at line 340 of file PFJetTester.h.

MonitorElement* PFJetTester::mHOTime

private

Definition at line 341 of file PFJetTester.h.

edm::InputTag PFJetTester::mInputCollection

private

Definition at line 37 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

edm::InputTag PFJetTester::mInputGenCollection

private

Definition at line 38 of file PFJetTester.h.

Referenced by analyze().

MonitorElement* PFJetTester::mjetArea

private

Definition at line 100 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mJetEnergyProfile

private

Definition at line 318 of file PFJetTester.h.

MonitorElement* PFJetTester::mlinEScale

private

Definition at line 257 of file PFJetTester.h.

MonitorElement* PFJetTester::mMass

private

Definition at line 68 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mMass_3000

private

Definition at line 70 of file PFJetTester.h.

MonitorElement* PFJetTester::mMass_80

private

Definition at line 69 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mMatchedGenJetsEta

private

Definition at line 250 of file PFJetTester.h.

MonitorElement* PFJetTester::mMatchedGenJetsPt

private

Definition at line 248 of file PFJetTester.h.

double PFJetTester::mMatchGenPtThreshold

private

Definition at line 305 of file PFJetTester.h.

Referenced by analyze().

MonitorElement* PFJetTester::mMjj

private

Definition at line 158 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mMjj_3000

private

Definition at line 159 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mMuonEnergy

private

Definition at line 224 of file PFJetTester.h.

MonitorElement* PFJetTester::mMuonEnergy_3000

private

Definition at line 226 of file PFJetTester.h.

MonitorElement* PFJetTester::mMuonEnergy_80

private

Definition at line 225 of file PFJetTester.h.

MonitorElement* PFJetTester::mNeutralEmEnergy_3000

private

Definition at line 169 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralEmEnergy_80

private

Definition at line 164 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralEmEnergyFraction_B

private

Definition at line 178 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralEmEnergyFraction_E

private

Definition at line 179 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralEmEnergyFraction_F

private

Definition at line 180 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralHadronEnergy_3000

private

Definition at line 170 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralHadronEnergy_80

private

Definition at line 165 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralHadronEnergyFraction_B

private

Definition at line 181 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralHadronEnergyFraction_E

private

Definition at line 182 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNeutralHadronEnergyFraction_F

private

Definition at line 183 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNJets1

private

Definition at line 154 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNJets2

private

Definition at line 155 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNJetsEtaC

private

Definition at line 151 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNJetsEtaF

private

Definition at line 152 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mNJetsEtaF_30

private

Definition at line 135 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

std::string PFJetTester::mOutputFile

private

Definition at line 39 of file PFJetTester.h.

Referenced by endJob(), and PFJetTester().

MonitorElement* PFJetTester::mP

private

Definition at line 62 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mP_3000

private

Definition at line 64 of file PFJetTester.h.

MonitorElement* PFJetTester::mP_80

private

Definition at line 63 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPhi

private

Definition at line 57 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPhiFineBin

private

Definition at line 58 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPhiFirst

private

Definition at line 143 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPhotonEnergy

private

Definition at line 227 of file PFJetTester.h.

MonitorElement* PFJetTester::mPhotonEnergy_3000

private

Definition at line 229 of file PFJetTester.h.

MonitorElement* PFJetTester::mPhotonEnergy_80

private

Definition at line 228 of file PFJetTester.h.

MonitorElement* PFJetTester::mPt

private

Definition at line 65 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPt_3000

private

Definition at line 67 of file PFJetTester.h.

MonitorElement* PFJetTester::mPt_80

private

Definition at line 66 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPtFirst

private

Definition at line 147 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPtFirst_3000

private

Definition at line 149 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPtFirst_80

private

Definition at line 148 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPthat_3000

private

Definition at line 233 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mPthat_80

private

Definition at line 232 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatio

private

Definition at line 83 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatio_1500_3500_d

private

Definition at line 90 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatio_200_600_d

private

Definition at line 88 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatio_30_200_d

private

Definition at line 87 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatio_600_1500_d

private

Definition at line 89 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatioB_d

private

Definition at line 84 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatioE_d

private

Definition at line 85 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTRatioF_d

private

Definition at line 86 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse

private

Definition at line 91 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_1500_3500_d

private

Definition at line 98 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_200_600_d

private

Definition at line 96 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_30_200_d

private

Definition at line 95 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_30_d

private

Definition at line 99 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_600_1500_d

private

Definition at line 97 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_nvtx_0_5

private

Definition at line 105 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_nvtx_10_15

private

Definition at line 107 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_nvtx_15_20

private

Definition at line 108 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_nvtx_20_30

private

Definition at line 109 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_nvtx_30_inf

private

Definition at line 110 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponse_nvtx_5_10

private

Definition at line 106 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponseB_d

private

Definition at line 92 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponseE_d

private

Definition at line 93 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTResponseF_d

private

Definition at line 94 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1D_1500_3500

private

Definition at line 298 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale1D_200_600

private

Definition at line 296 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale1D_30_200

private

Definition at line 295 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale1D_600_1500

private

Definition at line 297 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale1DB_1500_3500

private

Definition at line 292 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DB_200_600

private

Definition at line 286 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DB_30_200

private

Definition at line 283 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DB_600_1500

private

Definition at line 289 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DE_1500_3500

private

Definition at line 293 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DE_200_600

private

Definition at line 287 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DE_30_200

private

Definition at line 284 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DE_600_1500

private

Definition at line 290 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DF_1500_3500

private

Definition at line 294 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DF_200_600

private

Definition at line 288 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DF_30_200

private

Definition at line 285 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale1DF_600_1500

private

Definition at line 291 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_1500_3500_d

private

Definition at line 281 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_1500_3500_s

private

Definition at line 276 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale_200_600_d

private

Definition at line 279 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_200_600_s

private

Definition at line 274 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale_30_200_d

private

Definition at line 278 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_30_200_s

private

Definition at line 273 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale_600_1500_d

private

Definition at line 280 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_600_1500_s

private

Definition at line 275 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScale_a

private

Definition at line 136 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_a_nvtx_0_5

private

Definition at line 111 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_a_nvtx_10_15

private

Definition at line 117 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_a_nvtx_15_20

private

Definition at line 120 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_a_nvtx_20_30

private

Definition at line 123 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_a_nvtx_30_inf

private

Definition at line 126 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_a_nvtx_5_10

private

Definition at line 114 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b

private

Definition at line 137 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b_nvtx_0_5

private

Definition at line 112 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b_nvtx_10_15

private

Definition at line 118 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b_nvtx_15_20

private

Definition at line 121 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b_nvtx_20_30

private

Definition at line 124 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b_nvtx_30_inf

private

Definition at line 127 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_b_nvtx_5_10

private

Definition at line 115 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c

private

Definition at line 138 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c_nvtx_0_5

private

Definition at line 113 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c_nvtx_10_15

private

Definition at line 119 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c_nvtx_15_20

private

Definition at line 122 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c_nvtx_20_30

private

Definition at line 125 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c_nvtx_30_inf

private

Definition at line 128 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_c_nvtx_5_10

private

Definition at line 116 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_nvtx_0_5

private

Definition at line 129 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_nvtx_10_15

private

Definition at line 131 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_nvtx_15_20

private

Definition at line 132 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_nvtx_20_30

private

Definition at line 133 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_nvtx_30_inf

private

Definition at line 134 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_nvtx_5_10

private

Definition at line 130 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScale_pT

private

Definition at line 139 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScaleB_d

private

Definition at line 266 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScaleB_s

private

Definition at line 263 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScaleE_d

private

Definition at line 267 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScaleE_s

private

Definition at line 264 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScaleF_d

private

Definition at line 268 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScaleF_s

private

Definition at line 265 of file PFJetTester.h.

MonitorElement* PFJetTester::mpTScalePhiB_d

private

Definition at line 269 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScalePhiE_d

private

Definition at line 270 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::mpTScalePhiF_d

private

Definition at line 271 of file PFJetTester.h.

Referenced by fillMatchHists(), and PFJetTester().

double PFJetTester::mReverseEnergyFractionThreshold

private

Definition at line 307 of file PFJetTester.h.

MonitorElement* PFJetTester::mReverseMatchEnergyFraction

private

Definition at line 252 of file PFJetTester.h.

MonitorElement* PFJetTester::mRMatch

private

Definition at line 253 of file PFJetTester.h.

double PFJetTester::mRThreshold

private

Definition at line 308 of file PFJetTester.h.

Referenced by analyze().

std::string PFJetTester::mTurnOnEverything

private

Definition at line 313 of file PFJetTester.h.

Referenced by analyze(), fillMatchHists(), and PFJetTester().

MonitorElement* PFJetTester::numberofevents

private

Definition at line 46 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::nvtx_0_30

private

Definition at line 103 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().

MonitorElement* PFJetTester::nvtx_0_60

private

Definition at line 104 of file PFJetTester.h.

Referenced by analyze(), and PFJetTester().