CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

JPTJetTester Class Reference

#include <JPTJetTester.h>

Inheritance diagram for JPTJetTester:
edm::EDAnalyzer edm::EDConsumerBase

List of all members.

Public Member Functions

virtual void analyze (const edm::Event &, const edm::EventSetup &)
virtual void beginJob ()
virtual void endJob ()
 JPTJetTester (const edm::ParameterSet &)
 ~JPTJetTester ()

Private Member Functions

void fillMatchHists (const reco::GenJet &fGenJet, const reco::JPTJet &fJPTJet, std::vector< reco::Vertex > goodVertices)
double getSumPt (const reco::TrackRefVector &tracks)

Private Attributes

std::string inputCaloMETLabel_
std::string inputGenMETLabel_
edm::InputTag inputMETLabel_
std::string JetCorrectionService
MonitorElementmAllGenJetsEta
MonitorElementmAllGenJetsPt
MonitorElementmCaloMET
MonitorElementmCaloMET_3000
MonitorElementmCaloMETPhi
MonitorElementmCaloMETSig
MonitorElementmCaloMETSig_3000
MonitorElementmCaloMEx
MonitorElementmCaloMEx_3000
MonitorElementmCaloMEy
MonitorElementmCaloMEy_3000
MonitorElementmCaloSumET
MonitorElementmCaloSumET_3000
MonitorElementmCorrJetEta
MonitorElementmCorrJetPhi
MonitorElementmCorrJetPt
MonitorElementmCorrJetPt_3000
MonitorElementmCorrJetPt_80
MonitorElementmDelEta
MonitorElementmDelPhi
MonitorElementmDelPt
MonitorElementmDeltaE
MonitorElementmDeltaEta
MonitorElementmDeltaPhi
MonitorElementmE
MonitorElementmE_3000
MonitorElementmE_80
MonitorElementmEBEne
MonitorElementmEBTime
MonitorElementmEEEne
MonitorElementmEETime
MonitorElementmEFirst
MonitorElementmEFirst_3000
MonitorElementmEFirst_80
MonitorElementmEmEnergyProfile
MonitorElementmEmTiming
MonitorElementmEScale
MonitorElementmEScale_pt10
MonitorElementmEScaleFineBin
MonitorElementmEta
MonitorElementmEtaFineBin
MonitorElementmEtaFineBin1m
MonitorElementmEtaFineBin1p
MonitorElementmEtaFineBin2m
MonitorElementmEtaFineBin2p
MonitorElementmEtaFineBin3m
MonitorElementmEtaFineBin3p
MonitorElementmEtaFirst
std::string METType_
double mGenEnergyFractionThreshold
MonitorElementmGenEta
MonitorElementmGenEtaFirst
MonitorElementmGenJetMatchEnergyFraction
MonitorElementmGenPhi
MonitorElementmGenPhiFirst
MonitorElementmGenPt
MonitorElementmGenPt_3000
MonitorElementmGenPt_80
MonitorElementmHadEnergyProfile
MonitorElementmHadTiming
MonitorElementmHBEne
MonitorElementmHBTime
MonitorElementmHEEne
MonitorElementmHETime
MonitorElementmHFEne
MonitorElementmHFTime
MonitorElementmHOEne
MonitorElementmHOTime
edm::InputTag mInputCollection
edm::InputTag mInputGenCollection
MonitorElementmjetArea
MonitorElementmJetEnergyProfile
MonitorElementmlinEScale
MonitorElementmMass
MonitorElementmMass_3000
MonitorElementmMass_80
MonitorElementmMatchedGenJetsEta
MonitorElementmMatchedGenJetsPt
double mMatchGenPtThreshold
MonitorElementmMjj
MonitorElementmMjj_3000
MonitorElementmMultiplicityElecInVertexInCalo
MonitorElementmMultiplicityElecInVertexOutCalo
MonitorElementmMultiplicityElecOutVertexInCalo
MonitorElementmMultiplicityMuonInVertexInCalo
MonitorElementmMultiplicityMuonInVertexOutCalo
MonitorElementmMultiplicityMuonOutVertexInCalo
MonitorElementmMultiplicityPionInVertexInCalo
MonitorElementmMultiplicityPionInVertexOutCalo
MonitorElementmMultiplicityPionOutVertexInCalo
MonitorElementmNJets1
MonitorElementmNJets2
MonitorElementmNJetsEtaC
MonitorElementmNJetsEtaF
MonitorElementmNJetsEtaF_30
std::string mOutputFile
MonitorElementmP
MonitorElementmP_3000
MonitorElementmP_80
MonitorElementmPhi
MonitorElementmPhiFineBin
MonitorElementmPhiFirst
MonitorElementmPt
MonitorElementmPt_3000
MonitorElementmPt_80
MonitorElementmPtFirst
MonitorElementmPtFirst_3000
MonitorElementmPtFirst_80
MonitorElementmPtFractionElecInVertexInCalo
MonitorElementmPtFractionElecInVertexOutCalo
MonitorElementmPtFractionElecOutVertexInCalo
MonitorElementmPtFractionMuonInVertexInCalo
MonitorElementmPtFractionMuonInVertexOutCalo
MonitorElementmPtFractionMuonOutVertexInCalo
MonitorElementmPtFractionPionInVertexInCalo
MonitorElementmPtFractionPionInVertexOutCalo
MonitorElementmPtFractionPionOutVertexInCalo
MonitorElementmPthat_3000
MonitorElementmPthat_80
MonitorElementmpTRatio
MonitorElementmpTRatio_1500_3500_d
MonitorElementmpTRatio_200_600_d
MonitorElementmpTRatio_30_200_d
MonitorElementmpTRatio_600_1500_d
MonitorElementmpTRatioB_d
MonitorElementmpTRatioE_d
MonitorElementmpTRatioF_d
MonitorElementmpTResponse
MonitorElementmpTResponse_1500_3500_d
MonitorElementmpTResponse_200_600_d
MonitorElementmpTResponse_30_200_d
MonitorElementmpTResponse_30_d
MonitorElementmpTResponse_600_1500_d
MonitorElementmpTResponse_nvtx_0_5
MonitorElementmpTResponse_nvtx_10_15
MonitorElementmpTResponse_nvtx_15_20
MonitorElementmpTResponse_nvtx_20_30
MonitorElementmpTResponse_nvtx_30_inf
MonitorElementmpTResponse_nvtx_5_10
MonitorElementmpTResponseB_d
MonitorElementmpTResponseE_d
MonitorElementmpTResponseF_d
MonitorElementmpTScale1D_1500_3500
MonitorElementmpTScale1D_200_600
MonitorElementmpTScale1D_30_200
MonitorElementmpTScale1D_600_1500
MonitorElementmpTScale1DB_1500_3500
MonitorElementmpTScale1DB_200_600
MonitorElementmpTScale1DB_30_200
MonitorElementmpTScale1DB_600_1500
MonitorElementmpTScale1DE_1500_3500
MonitorElementmpTScale1DE_200_600
MonitorElementmpTScale1DE_30_200
MonitorElementmpTScale1DE_600_1500
MonitorElementmpTScale1DF_1500_3500
MonitorElementmpTScale1DF_200_600
MonitorElementmpTScale1DF_30_200
MonitorElementmpTScale1DF_600_1500
MonitorElementmpTScale_1500_3500_d
MonitorElementmpTScale_1500_3500_s
MonitorElementmpTScale_200_600_d
MonitorElementmpTScale_200_600_s
MonitorElementmpTScale_30_200_d
MonitorElementmpTScale_30_200_s
MonitorElementmpTScale_600_1500_d
MonitorElementmpTScale_600_1500_s
MonitorElementmpTScale_a
MonitorElementmpTScale_a_nvtx_0_5
MonitorElementmpTScale_a_nvtx_10_15
MonitorElementmpTScale_a_nvtx_15_20
MonitorElementmpTScale_a_nvtx_20_30
MonitorElementmpTScale_a_nvtx_30_inf
MonitorElementmpTScale_a_nvtx_5_10
MonitorElementmpTScale_b
MonitorElementmpTScale_b_nvtx_0_5
MonitorElementmpTScale_b_nvtx_10_15
MonitorElementmpTScale_b_nvtx_15_20
MonitorElementmpTScale_b_nvtx_20_30
MonitorElementmpTScale_b_nvtx_30_inf
MonitorElementmpTScale_b_nvtx_5_10
MonitorElementmpTScale_c
MonitorElementmpTScale_c_nvtx_0_5
MonitorElementmpTScale_c_nvtx_10_15
MonitorElementmpTScale_c_nvtx_15_20
MonitorElementmpTScale_c_nvtx_20_30
MonitorElementmpTScale_c_nvtx_30_inf
MonitorElementmpTScale_c_nvtx_5_10
MonitorElementmpTScale_nvtx_0_5
MonitorElementmpTScale_nvtx_10_15
MonitorElementmpTScale_nvtx_15_20
MonitorElementmpTScale_nvtx_20_30
MonitorElementmpTScale_nvtx_30_inf
MonitorElementmpTScale_nvtx_5_10
MonitorElementmpTScale_pT
MonitorElementmpTScaleB_d
MonitorElementmpTScaleB_s
MonitorElementmpTScaleE_d
MonitorElementmpTScaleE_s
MonitorElementmpTScaleF_d
MonitorElementmpTScaleF_s
MonitorElementmpTScalePhiB_d
MonitorElementmpTScalePhiE_d
MonitorElementmpTScalePhiF_d
double mReverseEnergyFractionThreshold
MonitorElementmReverseMatchEnergyFraction
MonitorElementmRMatch
double mRThreshold
std::string mTurnOnEverything
MonitorElementnumberofevents
MonitorElementnvtx_0_30
MonitorElementnvtx_0_60

Detailed Description

Definition at line 24 of file JPTJetTester.h.


Constructor & Destructor Documentation

JPTJetTester::JPTJetTester ( const edm::ParameterSet iConfig)

Definition at line 46 of file JPTJetTester.cc.

References DQMStore::book1D(), DQMStore::bookProfile(), MuonErrorMatrixValues_cff::etaRange, edm::InputTag::label(), mCorrJetEta, mCorrJetPhi, mCorrJetPt, mCorrJetPt_80, mDelEta, mDelPhi, mDelPt, mDeltaEta, mDeltaPhi, mE, mE_80, mEmEnergyProfile, mEmTiming, mEScale_pt10, mEScaleFineBin, mEta, mEtaFineBin, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_80, mHadEnergyProfile, mHadTiming, mInputCollection, mJetEnergyProfile, mMass, mMass_80, mMjj, mMjj_3000, 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 
    = mEtaFirst = mPhiFirst  = mPtFirst = mPtFirst_80 = mPtFirst_3000
    = mMjj = mMjj_3000 = mDelEta = mDelPhi = mDelPt 
    = mHadTiming = mEmTiming 
    = mNJetsEtaC = mNJetsEtaF = mNJets1 = mNJets2
      = mDeltaEta = mDeltaPhi 
    = mHadEnergyProfile = mEmEnergyProfile = mJetEnergyProfile 
    = mEScale_pt10 = mEScaleFineBin
    = mpTScaleB_d = mpTScaleE_d = mpTScaleF_d
      = mpTScalePhiB_d = mpTScalePhiE_d = mpTScalePhiF_d
    = 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
    = 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 
      = 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/JPTJetTask_" + 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); 
    //
    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, 5000); 
    //
    mP                = dbe->book1D("P", "P", 100, 0, 500); 
    mP_80             = dbe->book1D("P_80", "P_80", 100, 0, 5000); 
    //
    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); 
    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); 
    //

    //
    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);
    //
    //
    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");
    //
    //
    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",25,0,2.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 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};


   // int log10PtFineBins = 50;
    //
    if (mTurnOnEverything.compare("yes")==0) {
    }

    //
    if (mTurnOnEverything.compare("yes")==0) {
  }
    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",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mpTScalePhiE_d = dbe->bookProfile("pTScalePhiE_d", "pTScalePhi_d_1.5<|eta|<3.0",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    mpTScalePhiF_d = dbe->bookProfile("pTScalePhiF_d", "pTScalePhi_d_3.0<|eta|<6.0",
                                   log10PtBins, log10PtMin, log10PtMax, 0, 2, " ");
    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);

  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 (dbe)

  if (mOutputFile.empty ()) {
    LogInfo("OutputInfo") << " JPTJet histograms will NOT be saved";
  } 
  else {
    LogInfo("OutputInfo") << " JPTJethistograms will be saved to file:" << mOutputFile;
  }
}
JPTJetTester::~JPTJetTester ( )

Definition at line 325 of file JPTJetTester.cc.

{
}

Member Function Documentation

void JPTJetTester::analyze ( const edm::Event mEvent,
const edm::EventSetup mSetup 
) [virtual]

Implements edm::EDAnalyzer.

Definition at line 337 of file JPTJetTester.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(), mCorrJetEta, mCorrJetPhi, mCorrJetPt, mCorrJetPt_80, mE, mE_80, mEmTiming, mEta, mEtaFineBin, mEtaFirst, mGenEta, mGenEtaFirst, mGenPhi, mGenPhiFirst, mGenPt, mGenPt_80, mHadTiming, mInputCollection, mInputGenCollection, mMass, mMass_80, mMatchGenPtThreshold, mMjj, mMjj_3000, 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, reco::LeafCandidate::phi(), position, reco::LeafCandidate::pt(), rho, pileupReCalc_HLTpaths::scale, reco::Jet::scaleEnergy(), and 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());
  }
*/
  // ***********************************
  // *** 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) {
  }

    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;
      }
    }
  } 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;
      }
    }
  } 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++) {
      }
    }
  } 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++) {
      }
    }
  } 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++) {
      }
    }
  } catch (...) {
    edm::LogInfo("OutputInfo") << " No EE RecHits.";
  }

  //***********************************
  //*** Get the Jet collection
  //***********************************
  math::XYZTLorentzVector p4tmp[2];
  Handle<JPTJetCollection> jptJets;
  mEvent.getByLabel(mInputCollection, jptJets);
  if (!jptJets.isValid()) return;
  JPTJetCollection::const_iterator jet = jptJets->begin ();
  int jetIndex = 0;
  int nJet = 0;
  int nJetF = 0;
  int nJetC = 0;
  int nJetF_30 =0;
  for (; jet != jptJets->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 (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 (jet == jptJets->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 (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() );
  }

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

  for (JPTJetCollection::const_iterator jet = jptJets->begin(); jet !=jptJets ->end(); jet++) 
  {
 
      JPTJet  correctedJet = *jet;
      //double scale = corrector->correction(jet->p4()); 
      //double scale = corrector->correction(*jet);
      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());
  }

      }
  }
  // Count Jets above Pt cut
  for (int istep = 0; istep < 100; ++istep) {
    int     njet = 0;
    float ptStep = (istep * (200./100.));

    for ( JPTJetCollection::const_iterator jpt = jptJets->begin(); jpt != jptJets->end(); ++ jpt ) {
      if ( jpt->pt() > ptStep ) njet++;
    }
    mNJets1->Fill( ptStep, njet );
  }

  for (int istep = 0; istep < 100; ++istep) {
    int     njet = 0;
    float ptStep = (istep * (4000./100.));
    for ( JPTJetCollection::const_iterator jpt = jptJets->begin(); jpt != jptJets->end(); ++ jpt ) {
      if ( jpt->pt() > ptStep ) njet++;
    }
    mNJets2->Fill( ptStep, njet );
  }

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 JPTJets to GenJets
  JetMatchingTools jetMatching (mEvent);
  if (!(mInputGenCollection.label().empty())) {
    //    Handle<GenJetCollection> genJets;
    //    mEvent.getByLabel(mInputGenCollection, genJets);

    std::vector <std::vector <const reco::GenParticle*> > genJetConstituents (genJets->size());
    std::vector <std::vector <const reco::GenParticle*> > jptJetConstituents (jptJets->size());
    
    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 (jptJets->size() <= 0) continue; // no JPTJets - nothing to match
      if (mRThreshold > 0) {
        unsigned iJPTJetBest = 0;
        double deltaRBest = 999.;
        for (unsigned iJPTJet = 0; iJPTJet < jptJets->size(); ++iJPTJet) {
          double dR = deltaR (genJet.eta(), genJet.phi(), (*jptJets) [iJPTJet].eta(), (*jptJets) [iJPTJet].phi());
          if (deltaRBest < mRThreshold && dR < mRThreshold && genJet.pt() > 5.) {
            /*
            std::cout << "Yet another matched jet for GenJet pt=" << genJet.pt()
                      << " previous JPTJet pt/dr: " << (*jptJets) [iJPTJetBest].pt() << '/' << deltaRBest
                      << " new JPTJet pt/dr: " << (*jptJets) [iJPTJet].pt() << '/' << dR
                      << std::endl;
            */
          }
          if (dR < deltaRBest) {
            iJPTJetBest = iJPTJet;
            deltaRBest = dR;
          }
        }
        if (mTurnOnEverything.compare("yes")==0) {
          //mRMatch->Fill (logPtGen, genJet.eta(), deltaRBest);
        }
        if (deltaRBest < mRThreshold) { // Matched
          fillMatchHists (genJet, (*jptJets) [iJPTJetBest],goodVertices);
        }

        double CorrdeltaRBest = 999.;
        double CorrJetPtBest = 0;
        for (JPTJetCollection::const_iterator jet = jptJets->begin(); jet !=jptJets ->end(); jet++) {
          JPTJet  correctedJet = *jet;
          //double scale = corrector->correction(jet->p4());
          //double scale = corrector->correction(*jet);
          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());


        }
      }
      
    }
  }
}

}
void JPTJetTester::beginJob ( void  ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 329 of file JPTJetTester.cc.

                           {
}
void JPTJetTester::endJob ( void  ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 332 of file JPTJetTester.cc.

References mOutputFile, and cuy::save.

void JPTJetTester::fillMatchHists ( const reco::GenJet fGenJet,
const reco::JPTJet fJPTJet,
std::vector< reco::Vertex goodVertices 
) [private]

Referenced by analyze().

double JPTJetTester::getSumPt ( const reco::TrackRefVector tracks) [private]

Definition at line 910 of file JPTJetTester.cc.

References edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::end(), and reco::TrackBase::pt().

                                                             {

  double sumpt = 0.;
  
  for (reco::TrackRefVector::const_iterator itrack = tracks.begin(); itrack != tracks.end(); ++itrack){
    const reco::Track& track = **itrack;
    sumpt += track.pt(); 
  }
  
  return sumpt;

}

Member Data Documentation

std::string JPTJetTester::inputCaloMETLabel_ [private]

Definition at line 45 of file JPTJetTester.h.

std::string JPTJetTester::inputGenMETLabel_ [private]

Definition at line 44 of file JPTJetTester.h.

Definition at line 42 of file JPTJetTester.h.

std::string JPTJetTester::JetCorrectionService [private]

Definition at line 301 of file JPTJetTester.h.

Referenced by analyze().

Definition at line 240 of file JPTJetTester.h.

Definition at line 238 of file JPTJetTester.h.

Definition at line 320 of file JPTJetTester.h.

Definition at line 321 of file JPTJetTester.h.

Definition at line 322 of file JPTJetTester.h.

Definition at line 325 of file JPTJetTester.h.

Definition at line 326 of file JPTJetTester.h.

Definition at line 316 of file JPTJetTester.h.

Definition at line 317 of file JPTJetTester.h.

Definition at line 318 of file JPTJetTester.h.

Definition at line 319 of file JPTJetTester.h.

Definition at line 323 of file JPTJetTester.h.

Definition at line 324 of file JPTJetTester.h.

Definition at line 83 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 84 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 80 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 82 of file JPTJetTester.h.

Definition at line 81 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 291 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 292 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 293 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 249 of file JPTJetTester.h.

Definition at line 245 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 246 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 61 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 63 of file JPTJetTester.h.

Definition at line 62 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 337 of file JPTJetTester.h.

Definition at line 338 of file JPTJetTester.h.

Definition at line 339 of file JPTJetTester.h.

Definition at line 340 of file JPTJetTester.h.

Definition at line 146 of file JPTJetTester.h.

Definition at line 148 of file JPTJetTester.h.

Definition at line 147 of file JPTJetTester.h.

Definition at line 308 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 77 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 247 of file JPTJetTester.h.

Definition at line 251 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 252 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 51 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 52 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 56 of file JPTJetTester.h.

Definition at line 53 of file JPTJetTester.h.

Definition at line 57 of file JPTJetTester.h.

Definition at line 54 of file JPTJetTester.h.

Definition at line 58 of file JPTJetTester.h.

Definition at line 55 of file JPTJetTester.h.

Definition at line 144 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

std::string JPTJetTester::METType_ [private]

Definition at line 43 of file JPTJetTester.h.

Definition at line 297 of file JPTJetTester.h.

Definition at line 227 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 234 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 242 of file JPTJetTester.h.

Definition at line 228 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 235 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 229 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 231 of file JPTJetTester.h.

Definition at line 230 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 307 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 76 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 329 of file JPTJetTester.h.

Definition at line 330 of file JPTJetTester.h.

Definition at line 331 of file JPTJetTester.h.

Definition at line 332 of file JPTJetTester.h.

Definition at line 335 of file JPTJetTester.h.

Definition at line 336 of file JPTJetTester.h.

Definition at line 333 of file JPTJetTester.h.

Definition at line 334 of file JPTJetTester.h.

Definition at line 39 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 40 of file JPTJetTester.h.

Referenced by analyze().

Definition at line 102 of file JPTJetTester.h.

Definition at line 309 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 248 of file JPTJetTester.h.

Definition at line 70 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 72 of file JPTJetTester.h.

Definition at line 71 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 241 of file JPTJetTester.h.

Definition at line 239 of file JPTJetTester.h.

Definition at line 296 of file JPTJetTester.h.

Referenced by analyze().

Definition at line 160 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 161 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 185 of file JPTJetTester.h.

Definition at line 186 of file JPTJetTester.h.

Definition at line 187 of file JPTJetTester.h.

Definition at line 181 of file JPTJetTester.h.

Definition at line 182 of file JPTJetTester.h.

Definition at line 183 of file JPTJetTester.h.

Definition at line 177 of file JPTJetTester.h.

Definition at line 178 of file JPTJetTester.h.

Definition at line 179 of file JPTJetTester.h.

Definition at line 156 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 157 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 153 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 154 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 137 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

std::string JPTJetTester::mOutputFile [private]

Definition at line 41 of file JPTJetTester.h.

Referenced by endJob(), and JPTJetTester().

Definition at line 64 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 66 of file JPTJetTester.h.

Definition at line 65 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 59 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 60 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 145 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 67 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 69 of file JPTJetTester.h.

Definition at line 68 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 149 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 151 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 150 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 172 of file JPTJetTester.h.

Definition at line 173 of file JPTJetTester.h.

Definition at line 174 of file JPTJetTester.h.

Definition at line 168 of file JPTJetTester.h.

Definition at line 169 of file JPTJetTester.h.

Definition at line 170 of file JPTJetTester.h.

Definition at line 164 of file JPTJetTester.h.

Definition at line 165 of file JPTJetTester.h.

Definition at line 166 of file JPTJetTester.h.

Definition at line 224 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 223 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 85 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 92 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 90 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 89 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 91 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 86 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 87 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 88 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 93 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 100 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 98 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 97 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 101 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 99 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 107 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 109 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 110 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 111 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 112 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 108 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 94 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 95 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 96 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 289 of file JPTJetTester.h.

Definition at line 287 of file JPTJetTester.h.

Definition at line 286 of file JPTJetTester.h.

Definition at line 288 of file JPTJetTester.h.

Definition at line 283 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 277 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 274 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 280 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 284 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 278 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 275 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 281 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 285 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 279 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 276 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 282 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 272 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 267 of file JPTJetTester.h.

Definition at line 270 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 265 of file JPTJetTester.h.

Definition at line 269 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 264 of file JPTJetTester.h.

Definition at line 271 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 266 of file JPTJetTester.h.

Definition at line 138 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 113 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 119 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 122 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 125 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 128 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 116 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 139 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 114 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 120 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 123 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 126 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 129 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 117 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 140 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 115 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 121 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 124 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 127 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 130 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 118 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 131 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 133 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 134 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 135 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 136 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 132 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 141 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 257 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 254 of file JPTJetTester.h.

Definition at line 258 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 255 of file JPTJetTester.h.

Definition at line 259 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 256 of file JPTJetTester.h.

Definition at line 260 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 261 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 262 of file JPTJetTester.h.

Referenced by JPTJetTester().

Definition at line 298 of file JPTJetTester.h.

Definition at line 243 of file JPTJetTester.h.

Definition at line 244 of file JPTJetTester.h.

double JPTJetTester::mRThreshold [private]

Definition at line 299 of file JPTJetTester.h.

Referenced by analyze().

std::string JPTJetTester::mTurnOnEverything [private]

Definition at line 304 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 48 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 105 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().

Definition at line 106 of file JPTJetTester.h.

Referenced by analyze(), and JPTJetTester().