CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

JPTJetTesterUnCorr Class Reference

#include <JPTJetTesterUnCorr.h>

Inheritance diagram for JPTJetTesterUnCorr:
edm::EDAnalyzer

List of all members.

Public Member Functions

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

Private Member Functions

void fillMatchHists (const reco::GenJet &fGenJet, const reco::JPTJet &fJPTJet)
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
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
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
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

Detailed Description

Definition at line 22 of file JPTJetTesterUnCorr.h.


Constructor & Destructor Documentation

JPTJetTesterUnCorr::JPTJetTesterUnCorr ( const edm::ParameterSet iConfig)

Definition at line 49 of file JPTJetTesterUnCorr.cc.

References DQMStore::book1D(), DQMStore::bookProfile(), MuonErrorMatrixValues_cff::etaRange, edm::InputTag::label(), 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, mOutputFile, mP, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPthat_3000, mPthat_80, 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, mpTScaleB_d, mpTScaleE_d, mpTScaleF_d, mpTScalePhiB_d, mpTScalePhiE_d, mpTScalePhiF_d, mTurnOnEverything, numberofevents, and DQMStore::setCurrentFolder().

  : mInputCollection (iConfig.getParameter<edm::InputTag>( "src" )),
    mInputGenCollection (iConfig.getParameter<edm::InputTag>( "srcGen" )),
    mOutputFile (iConfig.getUntrackedParameter<std::string>("outputFile", "")),
    mMatchGenPtThreshold (iConfig.getParameter<double>("genPtThreshold")),
    mGenEnergyFractionThreshold (iConfig.getParameter<double>("genEnergyFractionThreshold")),
    mReverseEnergyFractionThreshold (iConfig.getParameter<double>("reverseEnergyFractionThreshold")),
    mRThreshold (iConfig.getParameter<double>("RThreshold")),
    mTurnOnEverything (iConfig.getUntrackedParameter<std::string>("TurnOnEverything",""))
{
    numberofevents
    = mEta = mEtaFineBin = mPhi = mPhiFineBin = mE = mE_80
    = 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 = mHadEnergyInHF = mHadEnergyInHE 
    = mHadEnergyInHO_80 = mHadEnergyInHB_80 = mHadEnergyInHE_80 
    = mHadEnergyInHO_3000 = mHadEnergyInHB_3000 = mHadEnergyInHE_3000 
    = mEmEnergyInEB = mEmEnergyInEE = mEmEnergyInHF 
    = mEmEnergyInEB_80 = mEmEnergyInEE_80
    = mEmEnergyInEB_3000 = mEmEnergyInEE_3000
    = mEnergyFractionHadronic = mEnergyFractionEm 
    = mHFLong = mHFTotal = mHFLong_80 = mHFLong_3000 = mHFShort = mHFShort_80 = mHFShort_3000
    = mN90
    */
    //
    //= mCaloMEx = mCaloMEx_3000 = mCaloMEy = mCaloMEy_3000 = mCaloMETSig = mCaloMETSig_3000
    //= mCaloMET = mCaloMET_3000 =  mCaloMETPhi = mCaloSumET  = mCaloSumET_3000   
    = mHadTiming = mEmTiming 
    = mNJetsEtaC = mNJetsEtaF = mNJets1 = mNJets2
      //= mAllGenJetsPt = mMatchedGenJetsPt = mAllGenJetsEta = mMatchedGenJetsEta 
      //= mGenJetMatchEnergyFraction = mReverseMatchEnergyFraction = mRMatch
      = mDeltaEta = mDeltaPhi //= mEScale = mlinEScale = mDeltaE
    = mHadEnergyProfile = mEmEnergyProfile = mJetEnergyProfile 
    /*
    = mHadJetEnergyProfile = mEMJetEnergyProfile
    */
    = mEScale_pt10 = mEScaleFineBin
      //= mpTScaleB_s = mpTScaleE_s = mpTScaleF_s 
    = mpTScaleB_d = mpTScaleE_d = mpTScaleF_d
      = 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 //=mjetArea
    = 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); 
    /*
    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, 50); 
    mHadEnergyInHE    = dbe->book1D("HadEnergyInHE", "HadEnergyInHE", 100, 0, 100); 
    //
    mHadEnergyInHO_80    = dbe->book1D("HadEnergyInHO_80", "HadEnergyInHO_80", 100, 0, 50); 
    mHadEnergyInHB_80    = dbe->book1D("HadEnergyInHB_80", "HadEnergyInHB_80", 100, 0, 200); 
    mHadEnergyInHE_80    = dbe->book1D("HadEnergyInHE_80", "HadEnergyInHE_80", 100, 0, 1000); 
    mHadEnergyInHO_3000  = dbe->book1D("HadEnergyInHO_3000", "HadEnergyInHO_3000", 100, 0, 500); 
    mHadEnergyInHB_3000  = dbe->book1D("HadEnergyInHB_3000", "HadEnergyInHB_3000", 100, 0, 3000); 
    mHadEnergyInHE_3000  = dbe->book1D("HadEnergyInHE_3000", "HadEnergyInHE_3000", 100, 0, 2000); 
    //
    mEmEnergyInEB     = dbe->book1D("EmEnergyInEB", "EmEnergyInEB", 100, 0, 50); 
    mEmEnergyInEE     = dbe->book1D("EmEnergyInEE", "EmEnergyInEE", 100, 0, 50); 
    mEmEnergyInHF     = dbe->book1D("EmEnergyInHF", "EmEnergyInHF", 120, -20, 100); 
    mEmEnergyInEB_80  = dbe->book1D("EmEnergyInEB_80", "EmEnergyInEB_80", 100, 0, 200); 
    mEmEnergyInEE_80  = dbe->book1D("EmEnergyInEE_80", "EmEnergyInEE_80", 100, 0, 1000); 
    mEmEnergyInEB_3000= dbe->book1D("EmEnergyInEB_3000", "EmEnergyInEB_3000", 100, 0, 3000); 
    mEmEnergyInEE_3000= dbe->book1D("EmEnergyInEE_3000", "EmEnergyInEE_3000", 100, 0, 2000); 
    mEnergyFractionHadronic = dbe->book1D("EnergyFractionHadronic", "EnergyFractionHadronic", 120, -0.1, 1.1); 
    mEnergyFractionEm = dbe->book1D("EnergyFractionEm", "EnergyFractionEm", 120, -0.1, 1.1); 
    //
    mHFTotal          = dbe->book1D("HFTotal", "HFTotal", 100, 0, 500);
    mHFTotal_80       = dbe->book1D("HFTotal_80", "HFTotal_80", 100, 0, 3000);
    mHFTotal_3000     = dbe->book1D("HFTotal_3000", "HFTotal_3000", 100, 0, 6000);
    mHFLong           = dbe->book1D("HFLong", "HFLong", 100, 0, 500);
    mHFLong_80        = dbe->book1D("HFLong_80", "HFLong_80", 100, 0, 200);
    mHFLong_3000      = dbe->book1D("HFLong_3000", "HFLong_3000", 100, 0, 1500);
    mHFShort          = dbe->book1D("HFShort", "HFShort", 100, 0, 500);
    mHFShort_80       = dbe->book1D("HFShort_80", "HFShort_80", 100, 0, 200);
    mHFShort_3000     = dbe->book1D("HFShort_3000", "HFShort_3000", 100, 0, 1500);
    //
    mN90              = dbe->book1D("N90", "N90", 50, 0, 50); 
    */
    //
    mGenEta           = dbe->book1D("GenEta", "GenEta", 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);
    //
    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); 
    //mjetArea = dbe->book1D("jetArea","jetArea",25,0,2.5);
    //
    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, " ");
    /*
    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);
    */
  } // if (dbe)

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

Definition at line 391 of file JPTJetTesterUnCorr.cc.

{
}

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 403 of file JPTJetTesterUnCorr.cc.

References HiRecoJets_cff::caloTowers, deltaR(), reco::LeafCandidate::eta(), MonitorElement::Fill(), fillMatchHists(), edm::Event::getByLabel(), edm::Event::getManyByType(), i, edm::EventBase::isRealData(), edm::HandleBase::isValid(), j, metsig::jet, edm::InputTag::label(), scaleCards::mass, 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, mP, mP_80, mPhi, mPhiFineBin, mPhiFirst, mPt, mPt_80, mPtFirst, mPtFirst_3000, mPtFirst_80, mPthat_3000, mPthat_80, mRThreshold, mTurnOnEverything, njet, numberofevents, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().

{
  double countsfornumberofevents = 1;
  numberofevents->Fill(countsfornumberofevents);
  // *********************************
  // *** 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<JPTJetCollection> jptJets;
  mEvent.getByLabel(mInputCollection, jptJets);
  if (!jptJets.isValid()) return;
  JPTJetCollection::const_iterator jet = jptJets->begin ();
  int jetIndex = 0;
  int nJet = 0;
  int nJetF = 0;
  int nJetC = 0;
  for (; jet != jptJets->end (); jet++, jetIndex++) {

    if (jet->pt() > 10.) {
      if (fabs(jet->eta()) > 1.5) 
        nJetF++;
      else 
        nJetC++;          
    }
    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 == 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 (mMaxEInEmTowers) mMaxEInEmTowers->Fill ( (jet->getCaloJetRef()->maxEInEmTowers());
    if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->getCaloJetRef()->maxEInHadTowers());
    if (mHadEnergyInHO) mHadEnergyInHO->Fill (jet->getCaloJetRef()->hadEnergyInHO());
    if (mHadEnergyInHO_80)   mHadEnergyInHO_80->Fill (jet->getCaloJetRef()->hadEnergyInHO());
    if (mHadEnergyInHO_3000) mHadEnergyInHO_3000->Fill (jet->getCaloJetRef()->hadEnergyInHO());
    if (mHadEnergyInHB) mHadEnergyInHB->Fill (jet->getCaloJetRef()->hadEnergyInHB());
    if (mHadEnergyInHB_80)   mHadEnergyInHB_80->Fill (jet->getCaloJetRef()->hadEnergyInHB());
    if (mHadEnergyInHB_3000) mHadEnergyInHB_3000->Fill (jet->getCaloJetRef()->hadEnergyInHB());
    if (mHadEnergyInHF) mHadEnergyInHF->Fill (jet->getCaloJetRef()->hadEnergyInHF());
    if (mHadEnergyInHE) mHadEnergyInHE->Fill (jet->getCaloJetRef()->hadEnergyInHE());
    if (mHadEnergyInHE_80)   mHadEnergyInHE_80->Fill (jet->getCaloJetRef()->hadEnergyInHE());
    if (mHadEnergyInHE_3000) mHadEnergyInHE_3000->Fill (jet->getCaloJetRef()->hadEnergyInHE());
    if (mEmEnergyInEB) mEmEnergyInEB->Fill (jet->getCaloJetRef()->emEnergyInEB());
    if (mEmEnergyInEB_80)   mEmEnergyInEB_80->Fill (jet->getCaloJetRef()->emEnergyInEB());
    if (mEmEnergyInEB_3000) mEmEnergyInEB_3000->Fill (jet->getCaloJetRef()->emEnergyInEB());
    if (mEmEnergyInEE) mEmEnergyInEE->Fill (jet->getCaloJetRef()->emEnergyInEE());
    if (mEmEnergyInEE_80)   mEmEnergyInEE_80->Fill (jet->getCaloJetRef()->emEnergyInEE());
    if (mEmEnergyInEE_3000) mEmEnergyInEE_3000->Fill (jet->getCaloJetRef()->emEnergyInEE());
    if (mEmEnergyInHF) mEmEnergyInHF->Fill (jet->getCaloJetRef()->emEnergyInHF());
    if (mEnergyFractionHadronic) mEnergyFractionHadronic->Fill (jet->getCaloJetRef()->energyFractionHadronic());
    if (mEnergyFractionEm) mEnergyFractionEm->Fill (jet->getCaloJetRef()->emEnergyFraction());

    if (mHFTotal)      mHFTotal->Fill (jet->getCaloJetRef()->hadEnergyInHF()+jet->getCaloJetRef()->emEnergyInHF());
    if (mHFTotal_80)   mHFTotal_80->Fill (jet->getCaloJetRef()->hadEnergyInHF()+jet->getCaloJetRef()->emEnergyInHF());
    if (mHFTotal_3000) mHFTotal_3000->Fill (jet->getCaloJetRef()->hadEnergyInHF()+jet->getCaloJetRef()->emEnergyInHF());
    if (mHFLong)       mHFLong->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5+jet->getCaloJetRef()->emEnergyInHF());
    if (mHFLong_80)    mHFLong_80->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5+jet->getCaloJetRef()->emEnergyInHF());
    if (mHFLong_3000)  mHFLong_3000->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5+jet->getCaloJetRef()->emEnergyInHF());
    if (mHFShort)      mHFShort->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5);
    if (mHFShort_80)   mHFShort_80->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5);
    if (mHFShort_3000) mHFShort_3000->Fill (jet->getCaloJetRef()->hadEnergyInHF()*0.5);
    */

    /*
    if (mN90) mN90->Fill (jet->getCaloJetRef()->n90());
    */
    //mJetEnergyProfile->Fill (jet->eta(), jet->phi(), jet->energy());
    /*
    mHadJetEnergyProfile->Fill (jet->eta(), jet->phi(), 
                                jet->getCaloJetRef()->hadEnergyInHO()+
                                jet->getCaloJetRef()->hadEnergyInHB()+
                                jet->getCaloJetRef()->hadEnergyInHF()+
                                jet->getCaloJetRef()->hadEnergyInHE());
    mEMJetEnergyProfile->Fill (jet->eta(), jet->phi(), 
                               jet->getCaloJetRef()->emEnergyInEB()+
                               jet->getCaloJetRef()->emEnergyInEE()+
                               jet->getCaloJetRef()->emEnergyInHF());
    */
  }

  if (mNJetsEtaC) mNJetsEtaC->Fill( nJetC );
  if (mNJetsEtaF) mNJetsEtaF->Fill( nJetF );

  if (nJet == 2) {
    if (mMjj) mMjj->Fill( (p4tmp[0]+p4tmp[1]).mass() );
    if (mMjj_3000) mMjj_3000->Fill( (p4tmp[0]+p4tmp[1]).mass() );
  }

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

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

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

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());
    //     if (mRThreshold > 0) { 
    //     }
    //     else {
    //       for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
    //  genJetConstituents [iGenJet] = jetMatching.getGenParticles ((*genJets) [iGenJet]);
    //       }
    
    //       for (unsigned iJPTJet = 0; iJPTJet < jptJets->size(); ++iJPTJet) {
    //  jptJetConstituents [iJPTJet] = jetMatching.getGenParticles ((*(*jptJets) [iJPTJet].getCaloJetRef()), false);
    //       }
    //     }
    
    for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {               //****************************************************************
    //for (unsigned iGenJet = 0; iGenJet < 1; ++iGenJet) {                           // only FIRST Jet !!!!
      const GenJet& genJet = (*genJets) [iGenJet];
      double genJetPt = genJet.pt();

      //std::cout << iGenJet <<". Genjet: pT = " << genJetPt << "GeV" << std::endl;  //  *****************************************************

      if (fabs(genJet.eta()) > 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]);
        }
      }
      //       else {
      //        unsigned iJPTJetBest = 0;
      //        double energyFractionBest = 0.;
      //        for (unsigned iJPTJet = 0; iJPTJet < jptJets->size(); ++iJPTJet) {
      //          double energyFraction = jetMatching.overlapEnergyFraction (genJetConstituents [iGenJet], 
      //                                                                     jptJetConstituents [iJPTJet]);
      //          if (energyFraction > energyFractionBest) {
      //            iJPTJetBest = iJPTJet;
      //            energyFractionBest = energyFraction;
      //          }
      //        }
      //        if (mTurnOnEverything.compare("yes")==0) {
      //          mGenJetMatchEnergyFraction->Fill (logPtGen, genJet.eta(), energyFractionBest);
      //        }
      //        if (energyFractionBest > mGenEnergyFractionThreshold) { // good enough
      //          double reverseEnergyFraction = jetMatching.overlapEnergyFraction (jptJetConstituents [iJPTJetBest], 
      //                                                                            genJetConstituents [iGenJet]);
      //          if (mTurnOnEverything.compare("yes")==0) {
      //            mReverseMatchEnergyFraction->Fill (logPtGen, genJet.eta(), reverseEnergyFraction);
      //          }
      //          if (reverseEnergyFraction > mReverseEnergyFractionThreshold) { // Matched
      //            fillMatchHists (genJet, (*jptJets) [iJPTJetBest]);
      //          }
      //        }
      //       } // else - mRThreshold
    }
  }
}

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

Reimplemented from edm::EDAnalyzer.

Definition at line 395 of file JPTJetTesterUnCorr.cc.

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

Reimplemented from edm::EDAnalyzer.

Definition at line 398 of file JPTJetTesterUnCorr.cc.

References mOutputFile.

void JPTJetTesterUnCorr::fillMatchHists ( const reco::GenJet fGenJet,
const reco::JPTJet fJPTJet 
) [private]

Definition at line 825 of file JPTJetTesterUnCorr.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, 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 PtJpt = fJPTJet.pt();
  //mMatchedGenJetsPt->Fill (logPtGen);
  //mMatchedGenJetsEta->Fill (logPtGen, fGenJet.eta());

  double PtThreshold = 10.;

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

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

    }

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

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

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

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

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

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

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

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

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

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

  if (fGenJet.pt()>1500.0 && fGenJet.pt()<3500.0) {
    //mpTScale_1500_3500_s->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
    mpTScale_1500_3500_d->Fill (fGenJet.eta(),fJPTJet.pt()/fGenJet.pt());
    //mpTScale1D_1500_3500->Fill (fJPTJet.pt()/fGenJet.pt());
  }



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

Definition at line 946 of file JPTJetTesterUnCorr.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

Definition at line 43 of file JPTJetTesterUnCorr.h.

Definition at line 42 of file JPTJetTesterUnCorr.h.

Definition at line 40 of file JPTJetTesterUnCorr.h.

Definition at line 251 of file JPTJetTesterUnCorr.h.

Definition at line 190 of file JPTJetTesterUnCorr.h.

Definition at line 188 of file JPTJetTesterUnCorr.h.

Definition at line 270 of file JPTJetTesterUnCorr.h.

Definition at line 271 of file JPTJetTesterUnCorr.h.

Definition at line 272 of file JPTJetTesterUnCorr.h.

Definition at line 275 of file JPTJetTesterUnCorr.h.

Definition at line 276 of file JPTJetTesterUnCorr.h.

Definition at line 266 of file JPTJetTesterUnCorr.h.

Definition at line 267 of file JPTJetTesterUnCorr.h.

Definition at line 268 of file JPTJetTesterUnCorr.h.

Definition at line 269 of file JPTJetTesterUnCorr.h.

Definition at line 273 of file JPTJetTesterUnCorr.h.

Definition at line 274 of file JPTJetTesterUnCorr.h.

Definition at line 81 of file JPTJetTesterUnCorr.h.

Definition at line 82 of file JPTJetTesterUnCorr.h.

Definition at line 78 of file JPTJetTesterUnCorr.h.

Definition at line 80 of file JPTJetTesterUnCorr.h.

Definition at line 79 of file JPTJetTesterUnCorr.h.

Definition at line 241 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 242 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 243 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 199 of file JPTJetTesterUnCorr.h.

Definition at line 195 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 196 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 59 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 61 of file JPTJetTesterUnCorr.h.

Definition at line 60 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 287 of file JPTJetTesterUnCorr.h.

Definition at line 288 of file JPTJetTesterUnCorr.h.

Definition at line 289 of file JPTJetTesterUnCorr.h.

Definition at line 290 of file JPTJetTesterUnCorr.h.

Definition at line 96 of file JPTJetTesterUnCorr.h.

Definition at line 98 of file JPTJetTesterUnCorr.h.

Definition at line 97 of file JPTJetTesterUnCorr.h.

Definition at line 258 of file JPTJetTesterUnCorr.h.

Referenced by JPTJetTesterUnCorr().

Definition at line 75 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 197 of file JPTJetTesterUnCorr.h.

Definition at line 201 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 202 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 49 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 50 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 54 of file JPTJetTesterUnCorr.h.

Definition at line 51 of file JPTJetTesterUnCorr.h.

Definition at line 55 of file JPTJetTesterUnCorr.h.

Definition at line 52 of file JPTJetTesterUnCorr.h.

Definition at line 56 of file JPTJetTesterUnCorr.h.

Definition at line 53 of file JPTJetTesterUnCorr.h.

Definition at line 94 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

std::string JPTJetTesterUnCorr::METType_ [private]

Definition at line 41 of file JPTJetTesterUnCorr.h.

Definition at line 247 of file JPTJetTesterUnCorr.h.

Definition at line 177 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 184 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 192 of file JPTJetTesterUnCorr.h.

Definition at line 178 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 185 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 179 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 181 of file JPTJetTesterUnCorr.h.

Definition at line 180 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 257 of file JPTJetTesterUnCorr.h.

Referenced by JPTJetTesterUnCorr().

Definition at line 74 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 279 of file JPTJetTesterUnCorr.h.

Definition at line 280 of file JPTJetTesterUnCorr.h.

Definition at line 281 of file JPTJetTesterUnCorr.h.

Definition at line 282 of file JPTJetTesterUnCorr.h.

Definition at line 285 of file JPTJetTesterUnCorr.h.

Definition at line 286 of file JPTJetTesterUnCorr.h.

Definition at line 283 of file JPTJetTesterUnCorr.h.

Definition at line 284 of file JPTJetTesterUnCorr.h.

Definition at line 37 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 38 of file JPTJetTesterUnCorr.h.

Referenced by analyze().

Definition at line 91 of file JPTJetTesterUnCorr.h.

Definition at line 259 of file JPTJetTesterUnCorr.h.

Referenced by JPTJetTesterUnCorr().

Definition at line 198 of file JPTJetTesterUnCorr.h.

Definition at line 68 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 70 of file JPTJetTesterUnCorr.h.

Definition at line 69 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 191 of file JPTJetTesterUnCorr.h.

Definition at line 189 of file JPTJetTesterUnCorr.h.

Definition at line 246 of file JPTJetTesterUnCorr.h.

Referenced by analyze().

Definition at line 110 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 111 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 135 of file JPTJetTesterUnCorr.h.

Definition at line 136 of file JPTJetTesterUnCorr.h.

Definition at line 137 of file JPTJetTesterUnCorr.h.

Definition at line 131 of file JPTJetTesterUnCorr.h.

Definition at line 132 of file JPTJetTesterUnCorr.h.

Definition at line 133 of file JPTJetTesterUnCorr.h.

Definition at line 127 of file JPTJetTesterUnCorr.h.

Definition at line 128 of file JPTJetTesterUnCorr.h.

Definition at line 129 of file JPTJetTesterUnCorr.h.

Definition at line 106 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 107 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 103 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 104 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

std::string JPTJetTesterUnCorr::mOutputFile [private]

Definition at line 39 of file JPTJetTesterUnCorr.h.

Referenced by endJob(), and JPTJetTesterUnCorr().

Definition at line 62 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 64 of file JPTJetTesterUnCorr.h.

Definition at line 63 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 57 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 58 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 95 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 65 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 67 of file JPTJetTesterUnCorr.h.

Definition at line 66 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 99 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 101 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 100 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 122 of file JPTJetTesterUnCorr.h.

Definition at line 123 of file JPTJetTesterUnCorr.h.

Definition at line 124 of file JPTJetTesterUnCorr.h.

Definition at line 118 of file JPTJetTesterUnCorr.h.

Definition at line 119 of file JPTJetTesterUnCorr.h.

Definition at line 120 of file JPTJetTesterUnCorr.h.

Definition at line 114 of file JPTJetTesterUnCorr.h.

Definition at line 115 of file JPTJetTesterUnCorr.h.

Definition at line 116 of file JPTJetTesterUnCorr.h.

Definition at line 174 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 173 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().

Definition at line 83 of file JPTJetTesterUnCorr.h.

Definition at line 90 of file JPTJetTesterUnCorr.h.

Definition at line 88 of file JPTJetTesterUnCorr.h.

Definition at line 87 of file JPTJetTesterUnCorr.h.

Definition at line 89 of file JPTJetTesterUnCorr.h.

Definition at line 84 of file JPTJetTesterUnCorr.h.

Definition at line 85 of file JPTJetTesterUnCorr.h.

Definition at line 86 of file JPTJetTesterUnCorr.h.

Definition at line 239 of file JPTJetTesterUnCorr.h.

Definition at line 237 of file JPTJetTesterUnCorr.h.

Definition at line 236 of file JPTJetTesterUnCorr.h.

Definition at line 238 of file JPTJetTesterUnCorr.h.

Definition at line 233 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 227 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 224 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 230 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 234 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 228 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 225 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 231 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 235 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 229 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 226 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 232 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 222 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 217 of file JPTJetTesterUnCorr.h.

Definition at line 220 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 215 of file JPTJetTesterUnCorr.h.

Definition at line 219 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 214 of file JPTJetTesterUnCorr.h.

Definition at line 221 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 216 of file JPTJetTesterUnCorr.h.

Definition at line 207 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 204 of file JPTJetTesterUnCorr.h.

Definition at line 208 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 205 of file JPTJetTesterUnCorr.h.

Definition at line 209 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 206 of file JPTJetTesterUnCorr.h.

Definition at line 210 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 211 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 212 of file JPTJetTesterUnCorr.h.

Referenced by fillMatchHists(), and JPTJetTesterUnCorr().

Definition at line 248 of file JPTJetTesterUnCorr.h.

Definition at line 193 of file JPTJetTesterUnCorr.h.

Definition at line 194 of file JPTJetTesterUnCorr.h.

Definition at line 249 of file JPTJetTesterUnCorr.h.

Referenced by analyze().

Definition at line 254 of file JPTJetTesterUnCorr.h.

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

Definition at line 46 of file JPTJetTesterUnCorr.h.

Referenced by analyze(), and JPTJetTesterUnCorr().