CMS 3D CMS Logo

Public Member Functions | Private Attributes

JetAnalyzer Class Reference

#include <JetAnalyzer.h>

Inheritance diagram for JetAnalyzer:
JetAnalyzerBase

List of all members.

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &, const reco::CaloJetCollection &caloJets)
 Get the analysis.
void beginJob (DQMStore *dbe)
 Inizialize parameters for histo binning.
void endJob ()
 Finish up a job.
int getLeadJetFlag ()
 JetAnalyzer (const edm::ParameterSet &)
 Constructor.
void setJetHiPass (int pass)
void setJetLoPass (int pass)
void setLeadJetFlag (int flag)
void setSource (std::string source)
virtual ~JetAnalyzer ()
 Destructor.

Private Attributes

double _fHPDMax
double _fHPDMaxLoose
double _fHPDMaxTight
int _JetHiPass
int _JetLoPass
int _leadJetFlag
int _n90HitsMin
int _n90HitsMinLoose
int _n90HitsMinTight
double _ptThreshold
double _resEMFMin
double _resEMFMinLoose
double _resEMFMinTight
double _sigmaEtaMinTight
double _sigmaPhiMinTight
std::string _source
int _theend
int eBin
double eMax
double eMin
int etaBin
double etaMax
double etaMin
int fillJIDPassFrac
reco::helper::JetIDHelperjetID
MonitorElementjetME
std::string jetname
int makedijetselection
MonitorElementmConstituents
MonitorElementmConstituents_Barrel
MonitorElementmConstituents_EndCap
MonitorElementmConstituents_Forward
MonitorElementmDPhi
MonitorElementmEFrac
MonitorElementmEFrac_Barrel
MonitorElementmEFrac_EndCap
MonitorElementmEFrac_Forward
MonitorElementmEmEnergyInEB
MonitorElementmEmEnergyInEE
MonitorElementmEmEnergyInHF
MonitorElementmEta
MonitorElementmEta_Hi
MonitorElementmEta_Lo
MonitorElementmEtaFirst
MonitorElementmfHPD
MonitorElementmfRBX
MonitorElementmHadEnergyInHB
MonitorElementmHadEnergyInHE
MonitorElementmHadEnergyInHF
MonitorElementmHadEnergyInHO
MonitorElementmHFrac
MonitorElementmHFrac_Barrel
MonitorElementmHFrac_EndCap
MonitorElementmHFrac_Forward
MonitorElementmLooseJIDPassFractionVSeta
MonitorElementmLooseJIDPassFractionVSpt
MonitorElementmMaxEInEmTowers
MonitorElementmMaxEInHadTowers
MonitorElementmN90Hits
MonitorElementmNJets
MonitorElementmPhi
MonitorElementmPhi_Barrel
MonitorElementmPhi_Barrel_Hi
MonitorElementmPhi_EndCap
MonitorElementmPhi_EndCap_Hi
MonitorElementmPhi_Forward
MonitorElementmPhi_Forward_Hi
MonitorElementmPhi_Hi
MonitorElementmPhi_Lo
MonitorElementmPhiFirst
MonitorElementmPhiVSEta
MonitorElementmPt
MonitorElementmPt_1
MonitorElementmPt_2
MonitorElementmPt_3
MonitorElementmPt_Barrel
MonitorElementmPt_Barrel_Hi
MonitorElementmPt_EndCap
MonitorElementmPt_EndCap_Hi
MonitorElementmPt_Forward
MonitorElementmPt_Forward_Hi
MonitorElementmPt_Hi
MonitorElementmPt_Lo
MonitorElementmPtFirst
MonitorElementmresEMF
MonitorElementmTightJIDPassFractionVSeta
MonitorElementmTightJIDPassFractionVSpt
edm::ParameterSet parameters
int pBin
int phiBin
double phiMax
double phiMin
double pMax
double pMin
int ptBin
double ptMax
double ptMin
edm::InputTag theCaloJetCollectionLabel

Detailed Description

DQM monitoring source for Calo Jets

Date:
2011/08/12 15:30:37
Revision:
1.13.6.1
Author:
F. Chlebana - Fermilab

Definition at line 37 of file JetAnalyzer.h.


Constructor & Destructor Documentation

JetAnalyzer::JetAnalyzer ( const edm::ParameterSet pSet)
JetAnalyzer::~JetAnalyzer ( ) [virtual]

Destructor.

Definition at line 44 of file JetAnalyzer.cc.

{ }

Member Function Documentation

void JetAnalyzer::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup,
const reco::CaloJetCollection caloJets 
)

Get the analysis.

Definition at line 244 of file JetAnalyzer.cc.

References corr, diffTreeTool::diff, eta(), metsig::jet, LogTrace, phi, and mathSSE::sqrt().

                                                                 {
  int numofjets=0;
  double  fstPhi=0.;
  double  sndPhi=0.;
  double  diff = 0.;
  double  corr = 0.;
  double  dphi = -999. ;
  bool thiscleaned=false;
  bool Loosecleaned=false;
  bool Tightcleaned=false;
  bool thisemfclean=true;
  bool emfcleanLoose=true;
  bool emfcleanTight=true;
  if(makedijetselection==1){
    //Dijet selection - careful: the pT is uncorrected!
    if(makedijetselection==1 && caloJets.size()>=2){
      double  dphiDJ = -999. ;
      bool emfcleanLooseFirstJet=true;
      bool emfcleanLooseSecondJet=true;
      bool emfcleanTightFirstJet=true;
      bool emfcleanTightSecondJet=true;
      bool LoosecleanedFirstJet = false;
      bool LoosecleanedSecondJet = false;
      bool TightcleanedFirstJet = false;
      bool TightcleanedSecondJet = false;
      //both jets pass pt threshold
      if ((caloJets.at(0)).pt() > _ptThreshold && (caloJets.at(1)).pt() > _ptThreshold ) {
        if(fabs((caloJets.at(0)).eta())<3. && fabs((caloJets.at(1)).eta())<3. ){
          //calculate dphi
          dphiDJ = fabs((caloJets.at(0)).phi()-(caloJets.at(1)).phi());
          if (dphiDJ > 3.14) dphiDJ=fabs(dphiDJ -6.28 );
          //fill DPhi histo (before cutting)
          if (mDPhi) mDPhi->Fill (dphiDJ);
          //dphi cut
          if(fabs(dphiDJ)>2.1){
            //JetID 
            emfcleanLooseFirstJet=true;
            emfcleanTightFirstJet=true;
            emfcleanLooseSecondJet=true;
            emfcleanTightSecondJet=true;
            jetID->calculate(iEvent, (caloJets.at(0)));
            if(jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(0)).eta())<2.6) emfcleanLooseFirstJet=false;
            if(jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(0)).eta())<2.6) emfcleanTightFirstJet=false;
            if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseFirstJet) LoosecleanedFirstJet=true;
            if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt((caloJets.at(0)).etaetaMoment())>_sigmaEtaMinTight && sqrt((caloJets.at(0)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightFirstJet) TightcleanedFirstJet=true;
            //fill the JID variables histograms BEFORE you cut on them
            if (mN90Hits)         mN90Hits->Fill (jetID->n90Hits());
            if (mfHPD)            mfHPD->Fill (jetID->fHPD());
            if (mresEMF)         mresEMF->Fill (jetID->restrictedEMF());
            if (mfRBX)            mfRBX->Fill (jetID->fRBX());
            jetID->calculate(iEvent, (caloJets.at(1)));
            if(jetID->restrictedEMF()<_resEMFMinLoose && fabs((caloJets.at(1)).eta())<2.6) emfcleanLooseSecondJet=false;
            if(jetID->restrictedEMF()<_resEMFMinTight && fabs((caloJets.at(1)).eta())<2.6) emfcleanTightSecondJet=false;
            if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLooseSecondJet) LoosecleanedSecondJet=true;
            if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt((caloJets.at(1)).etaetaMoment())>_sigmaEtaMinTight && sqrt((caloJets.at(1)).phiphiMoment())>_sigmaPhiMinTight && emfcleanTightSecondJet) TightcleanedSecondJet=true;
            //fill the JID variables histograms BEFORE you cut on them
            if (mN90Hits)         mN90Hits->Fill (jetID->n90Hits());
            if (mfHPD)            mfHPD->Fill (jetID->fHPD());
            if (mresEMF)         mresEMF->Fill (jetID->restrictedEMF());
            if (mfRBX)            mfRBX->Fill (jetID->fRBX());
            if(LoosecleanedFirstJet && LoosecleanedSecondJet) { //only if both jets are (loose) cleaned
              //fill histos for first jet
              if (mPt)   mPt->Fill ((caloJets.at(0)).pt());
              if (mEta)  mEta->Fill ((caloJets.at(0)).eta());
              if (mPhi)  mPhi->Fill ((caloJets.at(0)).phi());
              if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(0)).eta(),(caloJets.at(0)).phi());
              if (mConstituents) mConstituents->Fill ((caloJets.at(0)).nConstituents());
              if (mHFrac)        mHFrac->Fill ((caloJets.at(0)).energyFractionHadronic());
              if (mEFrac)        mEFrac->Fill ((caloJets.at(0)).emEnergyFraction());
              //removed for optimizations//if (mE)    mE->Fill ((caloJets.at(0)).energy());
              //removed for optimizations//if (mP)    mP->Fill ((caloJets.at(0)).p());
              //removed for optimizations//if (mMass) mMass->Fill ((caloJets.at(0)).mass());
              if (mMaxEInEmTowers)  mMaxEInEmTowers->Fill ((caloJets.at(0)).maxEInEmTowers());
              if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(0)).maxEInHadTowers());
              if (mN90Hits)         mN90Hits->Fill (jetID->n90Hits());
              if (mfHPD)            mfHPD->Fill (jetID->fHPD());
              if (mresEMF)         mresEMF->Fill (jetID->restrictedEMF());
              if (mfRBX)            mfRBX->Fill (jetID->fRBX());
              //sigmaeta and sigmaphi only used in the tight selection.
              //fill the histos for them AFTER the loose selection 
              //removed for optimizations//if (msigmaEta)  msigmaEta->Fill(sqrt((caloJets.at(0)).etaetaMoment()));
              //removed for optimizations//if (msigmaPhi)  msigmaPhi->Fill(sqrt((caloJets.at(0)).phiphiMoment()));
              //fill histos for second jet
              if (mPt)   mPt->Fill ((caloJets.at(1)).pt());
              if (mEta)  mEta->Fill ((caloJets.at(1)).eta());
              if (mPhi)  mPhi->Fill ((caloJets.at(1)).phi());
              if (mPhiVSEta) mPhiVSEta->Fill((caloJets.at(1)).eta(),(caloJets.at(1)).phi());
              if (mConstituents) mConstituents->Fill ((caloJets.at(1)).nConstituents());
              if (mHFrac)        mHFrac->Fill ((caloJets.at(1)).energyFractionHadronic());
              if (mEFrac)        mEFrac->Fill ((caloJets.at(1)).emEnergyFraction());
              //removed for optimizations//if (mE)    mE->Fill ((caloJets.at(1)).energy());
              //removed for optimizations//if (mP)    mP->Fill ((caloJets.at(1)).p());
              //removed for optimizations//if (mMass) mMass->Fill ((caloJets.at(1)).mass());
              if (mMaxEInEmTowers)  mMaxEInEmTowers->Fill ((caloJets.at(1)).maxEInEmTowers());
              if (mMaxEInHadTowers) mMaxEInHadTowers->Fill ((caloJets.at(1)).maxEInHadTowers());
              //sigmaeta and sigmaphi only used in the tight selection.
              //fill the histos for them AFTER the loose selection 
              //removed for optimizations//if (msigmaEta)  msigmaEta->Fill(sqrt((caloJets.at(1)).etaetaMoment()));
              //removed for optimizations//if (msigmaPhi)  msigmaPhi->Fill(sqrt((caloJets.at(1)).phiphiMoment()));

            }
            //let's see how many of these jets passed the JetID cleaning
            if(fillJIDPassFrac==1) {
              if(LoosecleanedFirstJet) {
                mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),1.);
                mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.);
              } else  {
                mLooseJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),0.);
                mLooseJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.);
              }
              if(LoosecleanedSecondJet) {
                mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),1.);
                mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.);
              } else  {
                mLooseJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),0.);
                mLooseJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.);
              }
              //TIGHT JID
              if(TightcleanedFirstJet) {
                mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),1.);
                mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),1.);
              } else  {
                mTightJIDPassFractionVSeta->Fill((caloJets.at(0)).eta(),0.);
                mTightJIDPassFractionVSpt->Fill((caloJets.at(0)).pt(),0.);
              }
              if(TightcleanedSecondJet) {
                mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),1.);
                mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),1.);
              } else  {
                mTightJIDPassFractionVSeta->Fill((caloJets.at(1)).eta(),0.);
                mTightJIDPassFractionVSpt->Fill((caloJets.at(1)).pt(),0.);
              }

            }//if fillJIDPassFrac
          }
        }
      }
    }
  } else {
    for (reco::CaloJetCollection::const_iterator jet = caloJets.begin(); jet!=caloJets.end(); ++jet) {
      LogTrace(jetname)<<"[JetAnalyzer] Analyze Calo Jet";
      Loosecleaned=false;
      Tightcleaned=false;
      if (jet == caloJets.begin()) {
        fstPhi = jet->phi();
        _leadJetFlag = 1;
      } else {
        _leadJetFlag = 0;
      }
      if (jet == (caloJets.begin()+1)) sndPhi = jet->phi();
      //jetID
      jetID->calculate(iEvent, *jet);
      //minimal (uncorrected!) pT cut
      if (jet->pt() > _ptThreshold) {
        //removed for optimizations//if (msigmaEta)  msigmaEta->Fill(sqrt(jet->etaetaMoment()));
        //removed for optimizations//if (msigmaPhi)  msigmaPhi->Fill(sqrt(jet->phiphiMoment()));
        //cleaning to use for filling histograms
        thisemfclean=true;
        if(jetID->restrictedEMF()<_resEMFMin && fabs(jet->eta())<2.6) thisemfclean=false;
        if(jetID->n90Hits()>=_n90HitsMin && jetID->fHPD()<_fHPDMax && thisemfclean) thiscleaned=true;
        //loose and tight cleaning, used to fill the JetIDPAssFraction histos
        if(jetID->n90Hits()>=_n90HitsMinLoose && jetID->fHPD()<_fHPDMaxLoose && emfcleanLoose) Loosecleaned=true;
        if(jetID->n90Hits()>=_n90HitsMinTight && jetID->fHPD()<_fHPDMaxTight && sqrt(jet->etaetaMoment())>_sigmaEtaMinTight && sqrt(jet->phiphiMoment())>_sigmaPhiMinTight && emfcleanTight) Tightcleaned=true;

        if(fillJIDPassFrac==1) {
          if(Loosecleaned) {
            mLooseJIDPassFractionVSeta->Fill(jet->eta(),1.);
            mLooseJIDPassFractionVSpt->Fill(jet->pt(),1.);
          } else {
            mLooseJIDPassFractionVSeta->Fill(jet->eta(),0.);
            mLooseJIDPassFractionVSpt->Fill(jet->pt(),0.);
          }
          //TIGHT
          if(Tightcleaned) {
            mTightJIDPassFractionVSeta->Fill(jet->eta(),1.);
            mTightJIDPassFractionVSpt->Fill(jet->pt(),1.);
          } else {
            mTightJIDPassFractionVSeta->Fill(jet->eta(),0.);
            mTightJIDPassFractionVSpt->Fill(jet->pt(),0.);
          }
        }
        //eventually we could define the "cleaned" flag differently for e.g. HF
        if(thiscleaned) {
          numofjets++ ;
          jetME->Fill(1);      
          
          // Leading jet
          // Histograms are filled once per event
          if (_leadJetFlag == 1) { 
            if (mEtaFirst) mEtaFirst->Fill (jet->eta());
            if (mPhiFirst) mPhiFirst->Fill (jet->phi());
            //removed for optimizations//if (mEFirst)   mEFirst->Fill (jet->energy());
            if (mPtFirst)  mPtFirst->Fill (jet->pt());
          }
          // --- Passed the low pt jet trigger
          if (_JetLoPass == 1) {
            /*if (fabs(jet->eta()) <= 1.3) {
              if (mPt_Barrel_Lo)           mPt_Barrel_Lo->Fill(jet->pt());
              //if (mEta_Lo)                 mEta_Lo->Fill(jet->eta());
              if (mPhi_Barrel_Lo)          mPhi_Barrel_Lo->Fill(jet->phi());
              //removed for optimization//if (mConstituents_Barrel_Lo) mConstituents_Barrel_Lo->Fill(jet->nConstituents());     
              //removed for optimization//if (mHFrac_Barrel_Lo)        mHFrac_Barrel_Lo->Fill(jet->energyFractionHadronic());   
            }
            if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
              if (mPt_EndCap_Lo)           mPt_EndCap_Lo->Fill(jet->pt());
              //if (mEta_Lo)                 mEta_Lo->Fill(jet->eta());
              if (mPhi_EndCap_Lo)          mPhi_EndCap_Lo->Fill(jet->phi());
              //removed for optimization//if (mConstituents_EndCap_Lo) mConstituents_EndCap_Lo->Fill(jet->nConstituents());     
              //removed for optimization//if (mHFrac_EndCap_Lo)        mHFrac_EndCap_Lo->Fill(jet->energyFractionHadronic());   
            }
            if (fabs(jet->eta()) > 3.0) {
              if (mPt_Forward_Lo)           mPt_Forward_Lo->Fill(jet->pt());
              //if (mEta_Lo)                  mEta_Lo->Fill(jet->eta());
              if (mPhi_Forward_Lo)          mPhi_Forward_Lo->Fill(jet->phi());
              //removed for optimization//if (mConstituents_Forward_Lo) mConstituents_Forward_Lo->Fill(jet->nConstituents());   
              //removed for optimization//if (mHFrac_Forward_Lo)        mHFrac_Forward_Lo->Fill(jet->energyFractionHadronic()); 
              }*/
            if (mEta_Lo) mEta_Lo->Fill (jet->eta());
            if (mPhi_Lo) mPhi_Lo->Fill (jet->phi());
            if (mPt_Lo)  mPt_Lo->Fill (jet->pt());
          }
          
          // --- Passed the high pt jet trigger
          if (_JetHiPass == 1) {
            if (fabs(jet->eta()) <= 1.3) {
              if (mPt_Barrel_Hi)           mPt_Barrel_Hi->Fill(jet->pt());
              //if (mEta_Hi)                 mEta_Hi->Fill(jet->eta());
              if (mPhi_Barrel_Hi)          mPhi_Barrel_Hi->Fill(jet->phi());
              //removed for optimizations//if (mConstituents_Barrel_Hi) mConstituents_Barrel_Hi->Fill(jet->nConstituents());    
              //removed for optimizations//if (mHFrac_Barrel_Hi)        mHFrac_Barrel_Hi->Fill(jet->energyFractionHadronic());  
            }
            if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
              if (mPt_EndCap_Hi)           mPt_EndCap_Hi->Fill(jet->pt());
              //if (mEta_Hi)                 mEta_Hi->Fill(jet->eta());
              if (mPhi_EndCap_Hi)          mPhi_EndCap_Hi->Fill(jet->phi());
              //removed for optimizations//if (mConstituents_EndCap_Hi) mConstituents_EndCap_Hi->Fill(jet->nConstituents());    
              //removed for optimizations//if (mHFrac_EndCap_Hi)        mHFrac_EndCap_Hi->Fill(jet->energyFractionHadronic());  
            }
            if (fabs(jet->eta()) > 3.0) {
              if (mPt_Forward_Hi)           mPt_Forward_Hi->Fill(jet->pt());
              //if (mEta_Hi)                  mEta_Hi->Fill(jet->eta());
              if (mPhi_Forward_Hi)          mPhi_Forward_Hi->Fill(jet->phi());
              //removed for optimizations//if (mConstituents_Forward_Hi) mConstituents_Forward_Hi->Fill(jet->nConstituents());  
              //removed for optimizations//if (mHFrac_Forward_Hi)        mHFrac_Forward_Hi->Fill(jet->energyFractionHadronic());        
            }
            
            if (mEta_Hi) mEta_Hi->Fill (jet->eta());
            if (mPhi_Hi) mPhi_Hi->Fill (jet->phi());
            if (mPt_Hi)  mPt_Hi->Fill (jet->pt());
          }
          
          if (mPt)   mPt->Fill (jet->pt());
          if (mPt_1) mPt_1->Fill (jet->pt());
          if (mPt_2) mPt_2->Fill (jet->pt());
          if (mPt_3) mPt_3->Fill (jet->pt());
          if (mEta)  mEta->Fill (jet->eta());
          if (mPhi)  mPhi->Fill (jet->phi());
          
          if (mPhiVSEta) mPhiVSEta->Fill(jet->eta(),jet->phi());
          
          if (mConstituents) mConstituents->Fill (jet->nConstituents());
          if (mHFrac)        mHFrac->Fill (jet->energyFractionHadronic());
          if (mEFrac)        mEFrac->Fill (jet->emEnergyFraction());
          
          if (fabs(jet->eta()) <= 1.3) {
            if (mPt_Barrel)   mPt_Barrel->Fill (jet->pt());
            if (mPhi_Barrel)  mPhi_Barrel->Fill (jet->phi());
            //removed for optimizations//if (mE_Barrel)    mE_Barrel->Fill (jet->energy());
            if (mConstituents_Barrel)    mConstituents_Barrel->Fill(jet->nConstituents());      
            if (mHFrac_Barrel)           mHFrac_Barrel->Fill(jet->energyFractionHadronic());    
            if (mEFrac_Barrel)           mEFrac_Barrel->Fill(jet->emEnergyFraction());  
          }
          if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
            if (mPt_EndCap)   mPt_EndCap->Fill (jet->pt());
            if (mPhi_EndCap)  mPhi_EndCap->Fill (jet->phi());
            //removed for optimizations//if (mE_EndCap)    mE_EndCap->Fill (jet->energy());
            if (mConstituents_EndCap)    mConstituents_EndCap->Fill(jet->nConstituents());      
            if (mHFrac_EndCap)           mHFrac_EndCap->Fill(jet->energyFractionHadronic());    
            if (mEFrac_EndCap)           mEFrac_EndCap->Fill(jet->emEnergyFraction());  
          }
          if (fabs(jet->eta()) > 3.0) {
            if (mPt_Forward)   mPt_Forward->Fill (jet->pt());
            if (mPhi_Forward)  mPhi_Forward->Fill (jet->phi());
            //removed for optimizations//if (mE_Forward)    mE_Forward->Fill (jet->energy());
            if (mConstituents_Forward)    mConstituents_Forward->Fill(jet->nConstituents());    
            if (mHFrac_Forward)           mHFrac_Forward->Fill(jet->energyFractionHadronic());  
            if (mEFrac_Forward)           mEFrac_Forward->Fill(jet->emEnergyFraction());        
          }
          
          //removed for optimizations//if (mE)    mE->Fill (jet->energy());
          //removed for optimizations//if (mP)    mP->Fill (jet->p());
          //removed for optimizations//if (mMass) mMass->Fill (jet->mass());
          
          if (mMaxEInEmTowers)  mMaxEInEmTowers->Fill (jet->maxEInEmTowers());
          if (mMaxEInHadTowers) mMaxEInHadTowers->Fill (jet->maxEInHadTowers());
          
          if (mHadEnergyInHO)   mHadEnergyInHO->Fill (jet->hadEnergyInHO());
          if (mHadEnergyInHB)   mHadEnergyInHB->Fill (jet->hadEnergyInHB());
          if (mHadEnergyInHF)   mHadEnergyInHF->Fill (jet->hadEnergyInHF());
          if (mHadEnergyInHE)   mHadEnergyInHE->Fill (jet->hadEnergyInHE());
          if (mEmEnergyInEB)    mEmEnergyInEB->Fill (jet->emEnergyInEB());
          if (mEmEnergyInEE)    mEmEnergyInEE->Fill (jet->emEnergyInEE());
          if (mEmEnergyInHF)    mEmEnergyInHF->Fill (jet->emEnergyInHF());
          
          if (mN90Hits)         mN90Hits->Fill (jetID->n90Hits());
          if (mfHPD)            mfHPD->Fill (jetID->fHPD());
          if (mresEMF)         mresEMF->Fill (jetID->restrictedEMF());
          if (mfRBX)            mfRBX->Fill (jetID->fRBX());
          
          //calculate correctly the dphi
          if(numofjets>1) {
            diff = fabs(fstPhi - sndPhi);
            corr = 2*acos(-1.) - diff;
            if(diff < acos(-1.)) { 
              dphi = diff; 
            } else { 
              dphi = corr;
            }
          }
        }
      }//pt cut
    }
    if (mNJets) mNJets->Fill (numofjets);
    if (mDPhi && dphi>-998.) mDPhi->Fill (dphi);
  }//not dijet
}
void JetAnalyzer::beginJob ( DQMStore dbe) [virtual]

Inizialize parameters for histo binning.

Implements JetAnalyzerBase.

Definition at line 48 of file JetAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), jptDQMConfig_cff::eMax, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, LogTrace, Parameters::parameters, jptDQMConfig_cff::phiMax, jptDQMConfig_cff::phiMin, jptDQMConfig_cff::pMax, jptDQMConfig_cff::ptMax, PtMinSelector_cfg::ptMin, MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().

                                         {
  
  jetname = "jetAnalyzer";
  
  LogTrace(jetname)<<"[JetAnalyzer] Parameters initialization";
  dbe->setCurrentFolder("JetMET/Jet/"+_source);

  jetME = dbe->book1D("jetReco", "jetReco", 3, 1, 4);
  jetME->setBinLabel(1,"CaloJets",1);

  //
  jetID = new reco::helper::JetIDHelper(parameters.getParameter<ParameterSet>("JetIDParams"));
  //

  fillJIDPassFrac = parameters.getParameter<int>("fillJIDPassFrac");
  makedijetselection = parameters.getParameter<int>("makedijetselection");

  // monitoring of eta parameter
  etaBin = parameters.getParameter<int>("etaBin");
  etaMin = parameters.getParameter<double>("etaMin");
  etaMax = parameters.getParameter<double>("etaMax");

  // monitoring of phi paramater
  phiBin = parameters.getParameter<int>("phiBin");
  phiMin = parameters.getParameter<double>("phiMin");
  phiMax = parameters.getParameter<double>("phiMax");

  // monitoring of the transverse momentum
  ptBin = parameters.getParameter<int>("ptBin");
  ptMin = parameters.getParameter<double>("ptMin");
  ptMax = parameters.getParameter<double>("ptMax");

  // 
  eBin = parameters.getParameter<int>("eBin");
  eMin = parameters.getParameter<double>("eMin");
  eMax = parameters.getParameter<double>("eMax");

  // 
  pBin = parameters.getParameter<int>("pBin");
  pMin = parameters.getParameter<double>("pMin");
  pMax = parameters.getParameter<double>("pMax");

  //
  _ptThreshold = parameters.getParameter<double>("ptThreshold");
  _n90HitsMin = parameters.getParameter<int>("n90HitsMin");
  _fHPDMax = parameters.getParameter<double>("fHPDMax");
  _resEMFMin = parameters.getParameter<double>("resEMFMin");
  _sigmaEtaMinTight = parameters.getParameter<double>("sigmaEtaMinTight");
  _sigmaPhiMinTight = parameters.getParameter<double>("sigmaPhiMinTight");

  _n90HitsMinLoose = parameters.getParameter<int>("n90HitsMinLoose");
  _fHPDMaxLoose = parameters.getParameter<double>("fHPDMaxLoose");
  _resEMFMinLoose = parameters.getParameter<double>("resEMFMinLoose");
  _n90HitsMinTight = parameters.getParameter<int>("n90HitsMinTight");
  _fHPDMaxTight = parameters.getParameter<double>("fHPDMaxTight");
  _resEMFMinTight = parameters.getParameter<double>("resEMFMinTight");

  
  // Generic Jet Parameters
  mPt                      = dbe->book1D("Pt",  "Pt", ptBin, ptMin, ptMax);
  mEta                     = dbe->book1D("Eta", "Eta", etaBin, etaMin, etaMax);
  mPhi                     = dbe->book1D("Phi", "Phi", phiBin, phiMin, phiMax);
  mConstituents            = dbe->book1D("Constituents", "# of Constituents", 50, 0, 100);
  mHFrac                   = dbe->book1D("HFrac", "HFrac", 120, -0.1, 1.1);
  mEFrac                   = dbe->book1D("EFrac", "EFrac", 120, -0.1, 1.1);
  //
  //removed for optimizations//mE                       = dbe->book1D("E", "E", eBin, eMin, eMax);
  //removed for optimizations//mP                       = dbe->book1D("P", "P", pBin, pMin, pMax);
  //removed for optimizations//mMass                    = dbe->book1D("Mass", "Mass", 100, 0, 25);
  //
  mPhiVSEta                     = dbe->book2D("PhiVSEta", "PhiVSEta", 50, etaMin, etaMax, 24, phiMin, phiMax);
  if(makedijetselection!=1){
    mPt_1                    = dbe->book1D("Pt1", "Pt1", 20, 0, 100);
    mPt_2                    = dbe->book1D("Pt2", "Pt2", 100, 0, 500);
    mPt_3                    = dbe->book1D("Pt3", "Pt3", 100, 0, 5000);
    // Low and high pt trigger paths
    mPt_Lo                  = dbe->book1D("Pt_Lo", "Pt (Pass Low Pt Jet Trigger)", 20, 0, 100);
    mEta_Lo                 = dbe->book1D("Eta_Lo", "Eta (Pass Low Pt Jet Trigger)", etaBin, etaMin, etaMax);
    mPhi_Lo                 = dbe->book1D("Phi_Lo", "Phi (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    
    mPt_Hi                  = dbe->book1D("Pt_Hi", "Pt (Pass Hi Pt Jet Trigger)", 100, 0, 500);
    mEta_Hi                 = dbe->book1D("Eta_Hi", "Eta (Pass Hi Pt Jet Trigger)", etaBin, etaMin, etaMax);
    mPhi_Hi                 = dbe->book1D("Phi_Hi", "Phi (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    mNJets                   = dbe->book1D("NJets", "Number of Jets", 100, 0, 100);

    mConstituents_Barrel     = dbe->book1D("Constituents_Barrel", "Constituents Barrel above", 50, 0, 100);
    mHFrac_Barrel            = dbe->book1D("HFrac_Barrel", "HFrac Barrel", 100, 0, 1);
    mEFrac_Barrel            = dbe->book1D("EFrac_Barrel", "EFrac Barrel", 110, -0.05, 1.05);
    //removed for optimization//mPt_Barrel_Lo            = dbe->book1D("Pt_Barrel_Lo", "Pt Barrel (Pass Low Pt Jet Trigger)", 20, 0, 100);
    //removed for optimization//mPhi_Barrel_Lo           = dbe->book1D("Phi_Barrel_Lo", "Phi Barrel (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //removed for optimization//mConstituents_Barrel_Lo  = dbe->book1D("Constituents_Barrel_Lo", "Constituents Barrel (Pass Low Pt Jet Trigger)", 50, 0, 100);
    //removed for optimization//mHFrac_Barrel_Lo         = dbe->book1D("HFrac_Barrel_Lo", "HFrac Barrel (Pass Low Pt Jet Trigger)", 100, 0, 1);
    
    mConstituents_EndCap     = dbe->book1D("Constituents_EndCap", "Constituents EndCap above", 50, 0, 100);
    mHFrac_EndCap            = dbe->book1D("HFrac_EndCap", "HFrac EndCap", 100, 0, 1);
    mEFrac_EndCap            = dbe->book1D("EFrac_EndCap", "EFrac EndCap", 110, -0.05, 1.05);
    //removed for optimization//mPt_EndCap_Lo            = dbe->book1D("Pt_EndCap_Lo", "Pt EndCap (Pass Low Pt Jet Trigger)", 20, 0, 100);
    //removed for optimization//mPhi_EndCap_Lo           = dbe->book1D("Phi_EndCap_Lo", "Phi EndCap (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //removed for optimization//mConstituents_EndCap_Lo  = dbe->book1D("Constituents_EndCap_Lo", "Constituents EndCap (Pass Low Pt Jet Trigger)", 50, 0, 100);
    //removed for optimization//mHFrac_EndCap_Lo         = dbe->book1D("HFrac_Endcap_Lo", "HFrac EndCap (Pass Low Pt Jet Trigger)", 100, 0, 1);
    
    mConstituents_Forward     = dbe->book1D("Constituents_Forward", "Constituents Forward above", 50, 0, 100);
    mHFrac_Forward            = dbe->book1D("HFrac_Forward", "HFrac Forward", 100, 0, 1);
    mEFrac_Forward            = dbe->book1D("EFrac_Forward", "EFrac Forward", 110, -0.05, 1.05);
    //removed for optimization//mPt_Forward_Lo           = dbe->book1D("Pt_Forward_Lo", "Pt Forward (Pass Low Pt Jet Trigger)", 20, 0, 100);
    //removed for optimization//mPhi_Forward_Lo          = dbe->book1D("Phi_Forward_Lo", "Phi Forward (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //removed for optimization//mConstituents_Forward_Lo = dbe->book1D("Constituents_Forward_Lo", "Constituents Forward (Pass Low Pt Jet Trigger)", 50, 0, 100);
    //removed for optimization//mHFrac_Forward_Lo        = dbe->book1D("HFrac_Forward_Lo", "HFrac Forward (Pass Low Pt Jet Trigger)", 100, 0, 1);
    
    mPt_Barrel_Hi            = dbe->book1D("Pt_Barrel_Hi", "Pt Barrel (Pass Hi Pt Jet Trigger)", 100, 0, 500);
    mPhi_Barrel_Hi           = dbe->book1D("Phi_Barrel_Hi", "Phi Barrel (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //removed for optimization//mConstituents_Barrel_Hi  = dbe->book1D("Constituents_Barrel_Hi", "Constituents Barrel (Pass Hi Pt Jet Trigger)", 50, 0, 100);
    //removed for optimization//mHFrac_Barrel_Hi         = dbe->book1D("HFrac_Barrel_Hi", "HFrac Barrel (Pass Hi Pt Jet Trigger)", 100, 0, 1);
    
    mPt_EndCap_Hi            = dbe->book1D("Pt_EndCap_Hi", "Pt EndCap (Pass Hi Pt Jet Trigger)", 100, 0, 500);
    mPhi_EndCap_Hi           = dbe->book1D("Phi_EndCap_Hi", "Phi EndCap (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //removed for optimization//mConstituents_EndCap_Hi  = dbe->book1D("Constituents_EndCap_Hi", "Constituents EndCap (Pass Hi Pt Jet Trigger)", 50, 0, 100);
    //removed for optimization//mHFrac_EndCap_Hi         = dbe->book1D("HFrac_EndCap_Hi", "HFrac EndCap (Pass Hi Pt Jet Trigger)", 100, 0, 1);
    
    mPt_Forward_Hi           = dbe->book1D("Pt_Forward_Hi", "Pt Forward (Pass Hi Pt Jet Trigger)", 100, 0, 500);
    mPhi_Forward_Hi          = dbe->book1D("Phi_Forward_Hi", "Phi Forward (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //removed for optimization//mConstituents_Forward_Hi = dbe->book1D("Constituents_Forward_Hi", "Constituents Forward (Pass Hi Pt Jet Trigger)", 50, 0, 100);
    //removed for optimization//mHFrac_Forward_Hi        = dbe->book1D("HFrac_Forward_Hi", "HFrac Forward (Pass Hi Pt Jet Trigger)", 100, 0, 1);
    
    mPhi_Barrel              = dbe->book1D("Phi_Barrel", "Phi_Barrel", phiBin, phiMin, phiMax);
    //removed for optimization//mE_Barrel                = dbe->book1D("E_Barrel", "E_Barrel", eBin, eMin, eMax);
    mPt_Barrel               = dbe->book1D("Pt_Barrel", "Pt_Barrel", ptBin, ptMin, ptMax);
    
    mPhi_EndCap              = dbe->book1D("Phi_EndCap", "Phi_EndCap", phiBin, phiMin, phiMax);
    //removed for optimization//mE_EndCap                = dbe->book1D("E_EndCap", "E_EndCap", eBin, eMin, 2*eMax);
    mPt_EndCap               = dbe->book1D("Pt_EndCap", "Pt_EndCap", ptBin, ptMin, ptMax);
    
    mPhi_Forward             = dbe->book1D("Phi_Forward", "Phi_Forward", phiBin, phiMin, phiMax);
    //removed for optimization//mE_Forward               = dbe->book1D("E_Forward", "E_Forward", eBin, eMin, 4*eMax);
    mPt_Forward              = dbe->book1D("Pt_Forward", "Pt_Forward", ptBin, ptMin, ptMax);
    
    // Leading Jet Parameters
    mEtaFirst                = dbe->book1D("EtaFirst", "EtaFirst", 100, -5, 5);
    mPhiFirst                = dbe->book1D("PhiFirst", "PhiFirst", 70, -3.5, 3.5);
    //removed for optimizations//mEFirst                  = dbe->book1D("EFirst", "EFirst", 100, 0, 1000);
    mPtFirst                 = dbe->book1D("PtFirst", "PtFirst", 100, 0, 500);
    if(fillJIDPassFrac==1){//fillJIDPassFrac defines a collection of cleaned jets, for which we will want to fill the cleaning passing fraction
      mLooseJIDPassFractionVSeta      = dbe->bookProfile("LooseJIDPassFractionVSeta","LooseJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2);
      mLooseJIDPassFractionVSpt       = dbe->bookProfile("LooseJIDPassFractionVSpt","LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
      mTightJIDPassFractionVSeta      = dbe->bookProfile("TightJIDPassFractionVSeta","TightJIDPassFractionVSeta",etaBin, etaMin, etaMax,0.,1.2);
      mTightJIDPassFractionVSpt       = dbe->bookProfile("TightJIDPassFractionVSpt","TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);


    }
  }
  // CaloJet specific
  mMaxEInEmTowers         = dbe->book1D("MaxEInEmTowers", "MaxEInEmTowers", 100, 0, 100);
  mMaxEInHadTowers        = dbe->book1D("MaxEInHadTowers", "MaxEInHadTowers", 100, 0, 100);
  if(makedijetselection!=1) {
    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);
    mEmEnergyInEB           = dbe->book1D("EmEnergyInEB", "EmEnergyInEB", 100, 0, 50);
    mEmEnergyInEE           = dbe->book1D("EmEnergyInEE", "EmEnergyInEE", 100, 0, 50);
    mEmEnergyInHF           = dbe->book1D("EmEnergyInHF", "EmEnergyInHF", 120, -20, 100);
  }
  mDPhi                   = dbe->book1D("DPhi", "dPhi btw the two leading jets", 100, 0., acos(-1.));
  
  //JetID variables
  
  mresEMF                 = dbe->book1D("resEMF", "resEMF", 50, 0., 1.);
  mN90Hits                = dbe->book1D("N90Hits", "N90Hits", 50, 0., 50);
  mfHPD                   = dbe->book1D("fHPD", "fHPD", 50, 0., 1.);
  mfRBX                   = dbe->book1D("fRBX", "fRBX", 50, 0., 1.);
  
  //removed for optimizations//msigmaEta                   = dbe->book1D("sigmaEta", "sigmaEta", 50, 0., 0.5);
  //removed for optimizations//msigmaPhi                   = dbe->book1D("sigmaPhi", "sigmaPhi", 50, 0., 0.5);
  
  if(makedijetselection==1 && fillJIDPassFrac==1) {
    mLooseJIDPassFractionVSeta  = dbe->bookProfile("LooseJIDPassFractionVSeta","LooseJIDPassFractionVSeta",50, -3., 3.,0.,1.2);
    mLooseJIDPassFractionVSpt   = dbe->bookProfile("LooseJIDPassFractionVSpt","LooseJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);
    mTightJIDPassFractionVSeta  = dbe->bookProfile("TightJIDPassFractionVSeta","TightJIDPassFractionVSeta",50, -3., 3.,0.,1.2);
    mTightJIDPassFractionVSpt   = dbe->bookProfile("TightJIDPassFractionVSpt","TightJIDPassFractionVSpt",ptBin, ptMin, ptMax,0.,1.2);

  }
}
void JetAnalyzer::endJob ( void  )

Finish up a job.

Definition at line 238 of file JetAnalyzer.cc.

                         {
  delete jetID;
}
int JetAnalyzer::getLeadJetFlag ( ) [inline]

Definition at line 64 of file JetAnalyzer.h.

References _leadJetFlag.

                       {
    return  _leadJetFlag;
  }
void JetAnalyzer::setJetHiPass ( int  pass) [inline]

Definition at line 71 of file JetAnalyzer.h.

References _JetHiPass.

                              {
    _JetHiPass = pass;
  }
void JetAnalyzer::setJetLoPass ( int  pass) [inline]

Definition at line 67 of file JetAnalyzer.h.

References _JetLoPass.

                              {
    _JetLoPass = pass;
  }
void JetAnalyzer::setLeadJetFlag ( int  flag) [inline]

Definition at line 61 of file JetAnalyzer.h.

References _leadJetFlag.

void JetAnalyzer::setSource ( std::string  source) [inline]

Definition at line 57 of file JetAnalyzer.h.

References _source, and LaserTracksInput_cfi::source.

                                   {
    _source = source;
  }

Member Data Documentation

double JetAnalyzer::_fHPDMax [private]

Definition at line 94 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMaxLoose [private]

Definition at line 99 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMaxTight [private]

Definition at line 102 of file JetAnalyzer.h.

int JetAnalyzer::_JetHiPass [private]

Definition at line 86 of file JetAnalyzer.h.

Referenced by setJetHiPass().

int JetAnalyzer::_JetLoPass [private]

Definition at line 85 of file JetAnalyzer.h.

Referenced by setJetLoPass().

Definition at line 87 of file JetAnalyzer.h.

Referenced by getLeadJetFlag(), and setLeadJetFlag().

int JetAnalyzer::_n90HitsMin [private]

Definition at line 96 of file JetAnalyzer.h.

Definition at line 101 of file JetAnalyzer.h.

Definition at line 104 of file JetAnalyzer.h.

double JetAnalyzer::_ptThreshold [private]

Definition at line 89 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMin [private]

Definition at line 95 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMinLoose [private]

Definition at line 100 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMinTight [private]

Definition at line 103 of file JetAnalyzer.h.

Definition at line 105 of file JetAnalyzer.h.

Definition at line 106 of file JetAnalyzer.h.

std::string JetAnalyzer::_source [private]

Definition at line 81 of file JetAnalyzer.h.

Referenced by setSource().

int JetAnalyzer::_theend [private]

Definition at line 88 of file JetAnalyzer.h.

int JetAnalyzer::eBin [private]

Definition at line 121 of file JetAnalyzer.h.

double JetAnalyzer::eMax [private]

Definition at line 123 of file JetAnalyzer.h.

double JetAnalyzer::eMin [private]

Definition at line 122 of file JetAnalyzer.h.

int JetAnalyzer::etaBin [private]

Definition at line 109 of file JetAnalyzer.h.

double JetAnalyzer::etaMax [private]

Definition at line 111 of file JetAnalyzer.h.

double JetAnalyzer::etaMin [private]

Definition at line 110 of file JetAnalyzer.h.

Definition at line 98 of file JetAnalyzer.h.

Definition at line 133 of file JetAnalyzer.h.

Definition at line 130 of file JetAnalyzer.h.

std::string JetAnalyzer::jetname [private]

Definition at line 80 of file JetAnalyzer.h.

Definition at line 91 of file JetAnalyzer.h.

Definition at line 155 of file JetAnalyzer.h.

Definition at line 169 of file JetAnalyzer.h.

Definition at line 176 of file JetAnalyzer.h.

Definition at line 183 of file JetAnalyzer.h.

Definition at line 214 of file JetAnalyzer.h.

Definition at line 157 of file JetAnalyzer.h.

Definition at line 171 of file JetAnalyzer.h.

Definition at line 178 of file JetAnalyzer.h.

Definition at line 185 of file JetAnalyzer.h.

Definition at line 230 of file JetAnalyzer.h.

Definition at line 231 of file JetAnalyzer.h.

Definition at line 232 of file JetAnalyzer.h.

Definition at line 153 of file JetAnalyzer.h.

Definition at line 251 of file JetAnalyzer.h.

Definition at line 247 of file JetAnalyzer.h.

Definition at line 217 of file JetAnalyzer.h.

Definition at line 236 of file JetAnalyzer.h.

Definition at line 237 of file JetAnalyzer.h.

Definition at line 227 of file JetAnalyzer.h.

Definition at line 229 of file JetAnalyzer.h.

Definition at line 228 of file JetAnalyzer.h.

Definition at line 226 of file JetAnalyzer.h.

Definition at line 156 of file JetAnalyzer.h.

Definition at line 170 of file JetAnalyzer.h.

Definition at line 177 of file JetAnalyzer.h.

Definition at line 184 of file JetAnalyzer.h.

Definition at line 241 of file JetAnalyzer.h.

Definition at line 242 of file JetAnalyzer.h.

Definition at line 224 of file JetAnalyzer.h.

Definition at line 225 of file JetAnalyzer.h.

Definition at line 235 of file JetAnalyzer.h.

Definition at line 213 of file JetAnalyzer.h.

Definition at line 154 of file JetAnalyzer.h.

Definition at line 161 of file JetAnalyzer.h.

Definition at line 192 of file JetAnalyzer.h.

Definition at line 164 of file JetAnalyzer.h.

Definition at line 196 of file JetAnalyzer.h.

Definition at line 167 of file JetAnalyzer.h.

Definition at line 200 of file JetAnalyzer.h.

Definition at line 252 of file JetAnalyzer.h.

Definition at line 248 of file JetAnalyzer.h.

Definition at line 218 of file JetAnalyzer.h.

Definition at line 158 of file JetAnalyzer.h.

Definition at line 149 of file JetAnalyzer.h.

Definition at line 150 of file JetAnalyzer.h.

Definition at line 151 of file JetAnalyzer.h.

Definition at line 152 of file JetAnalyzer.h.

Definition at line 160 of file JetAnalyzer.h.

Definition at line 191 of file JetAnalyzer.h.

Definition at line 163 of file JetAnalyzer.h.

Definition at line 195 of file JetAnalyzer.h.

Definition at line 166 of file JetAnalyzer.h.

Definition at line 199 of file JetAnalyzer.h.

Definition at line 253 of file JetAnalyzer.h.

Definition at line 249 of file JetAnalyzer.h.

Definition at line 220 of file JetAnalyzer.h.

Definition at line 238 of file JetAnalyzer.h.

Definition at line 243 of file JetAnalyzer.h.

Definition at line 244 of file JetAnalyzer.h.

Definition at line 78 of file JetAnalyzer.h.

int JetAnalyzer::pBin [private]

Definition at line 125 of file JetAnalyzer.h.

int JetAnalyzer::phiBin [private]

Definition at line 113 of file JetAnalyzer.h.

double JetAnalyzer::phiMax [private]

Definition at line 115 of file JetAnalyzer.h.

double JetAnalyzer::phiMin [private]

Definition at line 114 of file JetAnalyzer.h.

double JetAnalyzer::pMax [private]

Definition at line 127 of file JetAnalyzer.h.

double JetAnalyzer::pMin [private]

Definition at line 126 of file JetAnalyzer.h.

int JetAnalyzer::ptBin [private]

Definition at line 117 of file JetAnalyzer.h.

double JetAnalyzer::ptMax [private]

Definition at line 119 of file JetAnalyzer.h.

double JetAnalyzer::ptMin [private]

Definition at line 118 of file JetAnalyzer.h.

Definition at line 83 of file JetAnalyzer.h.