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 _asymmetryThirdJetCut
double _balanceThirdJetCut
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_Barrel_Hi
MonitorElementmConstituents_EndCap
MonitorElementmConstituents_EndCap_Hi
MonitorElementmConstituents_Forward
MonitorElementmConstituents_Forward_Hi
MonitorElementmDijetAsymmetry
MonitorElementmDijetBalance
MonitorElementmDPhi
MonitorElementmEFrac
MonitorElementmEFrac_Barrel
MonitorElementmEFrac_EndCap
MonitorElementmEFrac_Forward
MonitorElementmEmEnergyInEB
MonitorElementmEmEnergyInEE
MonitorElementmEmEnergyInHF
MonitorElementmEta
MonitorElementmEta_Hi
MonitorElementmEtaFirst
MonitorElementmfHPD
MonitorElementmfRBX
MonitorElementmHadEnergyInHB
MonitorElementmHadEnergyInHE
MonitorElementmHadEnergyInHF
MonitorElementmHadEnergyInHO
MonitorElementmHFrac
MonitorElementmHFrac_Barrel
MonitorElementmHFrac_Barrel_Hi
MonitorElementmHFrac_EndCap
MonitorElementmHFrac_EndCap_Hi
MonitorElementmHFrac_Forward
MonitorElementmHFrac_Forward_Hi
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/07/20 13:59:17
Revision:
1.15
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 46 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 245 of file JetAnalyzer.cc.

References corr, diffTreeTool::diff, dPhi(), eta(), edm::EventID::event(), edm::EventBase::id(), 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;

  srand( iEvent.id().event() % 10000);
  
  if(makedijetselection==1){
    //Dijet selection - careful: the pT is uncorrected!
    //if(makedijetselection==1 && caloJets.size()>=2){
    if(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 for first jet
            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 for second jet
            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());
              //if (mE)    mE->Fill ((caloJets.at(0)).energy());
              //if (mP)    mP->Fill ((caloJets.at(0)).p());
              //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 
              //  if (msigmaEta)  msigmaEta->Fill(sqrt((caloJets.at(0)).etaetaMoment()));
              //  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());
              //if (mE)    mE->Fill ((caloJets.at(1)).energy());
              //if (mP)    mP->Fill ((caloJets.at(1)).p());
              //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 
              //  if (msigmaEta)  msigmaEta->Fill(sqrt((caloJets.at(1)).etaetaMoment()));
              //  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
          }// FABS DPHI < 2.1
        }// fabs eta < 3
      }// pt jets > threshold
      //now do the dijet balance and asymmetry calculations
      if (fabs(caloJets.at(0).eta() < 1.4)) {
        double pt_dijet = (caloJets.at(0).pt() + caloJets.at(1).pt())/2;
        
        double dPhi = fabs((caloJets.at(0)).phi()-(caloJets.at(1)).phi());
        if (dPhi > 3.14) dPhi=fabs(dPhi -6.28 );
        
        if (dPhi > 2.7) {
          double pt_probe;
          double pt_barrel;
          int jet1, jet2;

          int randJet = rand() % 2;

          if (fabs(caloJets.at(1).eta() < 1.4)) {
            if (randJet) {
              jet1 = 0;
              jet2 = 1;
            }
            else {
              jet1 = 1;
              jet2 = 0;
            }
          
            /***Di-Jet Asymmetry****
             * leading jets eta < 1.4
             * leading jets dphi > 2.7
             * pt_third jet < threshold
             * A = (pt_1 - pt_2)/(pt_1 + pt_2)
             * jets 1 and two are randomly ordered
             */
            bool thirdJetCut = true;
            for (unsigned int third = 2; third < caloJets.size(); ++third) 
              if (caloJets.at(third).pt() > _asymmetryThirdJetCut) 
                thirdJetCut = false;
            if (thirdJetCut) {
              double dijetAsymmetry = (caloJets.at(jet1).pt() - caloJets.at(jet2).pt()) / (caloJets.at(jet1).pt() + caloJets.at(jet2).pt());
              mDijetAsymmetry->Fill(dijetAsymmetry);
            }// end restriction on third jet pt in asymmetry calculation
              
          }
          else {
            jet1 = 0;
            jet2 = 1;
          }
          
          pt_barrel = caloJets.at(jet1).pt();
          pt_probe  = caloJets.at(jet2).pt();
          
          //dijet balance cuts
          /***Di-Jet Balance****
           * pt_dijet = (pt_probe+pt_barrel)/2
           * leading jets dphi > 2.7
           * reject evnets where pt_third/pt_dijet > 0.2
           * pv selection
           * B = (pt_probe - pt_barrel)/pt_dijet
           * select probe randomly from 2 jets if both leading jets are in the barrel
           */
          bool thirdJetCut = true;
          for (unsigned int third = 2; third < caloJets.size(); ++third) 
            if (caloJets.at(third).pt()/pt_dijet > _balanceThirdJetCut) 
              thirdJetCut = false;
          if (thirdJetCut) {
            double dijetBalance = (pt_probe - pt_barrel) / pt_dijet;
            mDijetBalance->Fill(dijetBalance);
          }// end restriction on third jet pt ratio in balance calculation
        }// dPhi > 2.7
      }// leading jet eta cut for asymmetry and balance calculations
    }//jet size >= 2
  }// do dijet selection
  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) {
        //  if (msigmaEta)  msigmaEta->Fill(sqrt(jet->etaetaMoment()));
        //  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());
            //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());
            }
            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());
            }
            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());
            } */
            //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 && jet->pt()>100.)           mPt_Barrel_Hi->Fill(jet->pt());
              if (mEta_Hi && jet->pt()>100.)          mEta_Hi->Fill(jet->eta());
              if (mPhi_Barrel_Hi)          mPhi_Barrel_Hi->Fill(jet->phi());
              //if (mConstituents_Barrel_Hi) mConstituents_Barrel_Hi->Fill(jet->nConstituents());       
              //if (mHFrac_Barrel_Hi)        mHFrac_Barrel_Hi->Fill(jet->energyFractionHadronic());     
            }
            if ( (fabs(jet->eta()) > 1.3) && (fabs(jet->eta()) <= 3) ) {
              if (mPt_EndCap_Hi && jet->pt()>100.)           mPt_EndCap_Hi->Fill(jet->pt());
              if (mEta_Hi && jet->pt()>100.)          mEta_Hi->Fill(jet->eta());
              if (mPhi_EndCap_Hi)          mPhi_EndCap_Hi->Fill(jet->phi());
              //if (mConstituents_EndCap_Hi) mConstituents_EndCap_Hi->Fill(jet->nConstituents());       
              //if (mHFrac_EndCap_Hi)        mHFrac_EndCap_Hi->Fill(jet->energyFractionHadronic());     
            }
            if (fabs(jet->eta()) > 3.0) {
              if (mPt_Forward_Hi && jet->pt()>100.)           mPt_Forward_Hi->Fill(jet->pt());
              if (mEta_Hi && jet->pt()>100.)          mEta_Hi->Fill(jet->eta());
              if (mPhi_Forward_Hi)          mPhi_Forward_Hi->Fill(jet->phi());
              //if (mConstituents_Forward_Hi) mConstituents_Forward_Hi->Fill(jet->nConstituents());     
              //if (mHFrac_Forward_Hi)        mHFrac_Forward_Hi->Fill(jet->energyFractionHadronic());   
            }
            
            if (mEta_Hi && jet->pt()>100.) 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());
            //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());
            //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());
            //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());      
          }
          
          //if (mE)    mE->Fill (jet->energy());
          //if (mP)    mP->Fill (jet->p());
          // 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 50 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");
  _asymmetryThirdJetCut = parameters.getParameter<double>("asymmetryThirdJetCut");
  _balanceThirdJetCut   = parameters.getParameter<double>("balanceThirdJetCut");
  _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);
  //
  //mE                       = dbe->book1D("E", "E", eBin, eMin, eMax);
  //mP                       = dbe->book1D("P", "P", pBin, pMin, pMax);
  //  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", 60, 0, 300);   
    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)", 60, 0, 300);   
    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);

    //mPt_Barrel_Lo            = dbe->book1D("Pt_Barrel_Lo", "Pt Barrel (Pass Low Pt Jet Trigger)", 20, 0, 100);   
    //mPhi_Barrel_Lo           = dbe->book1D("Phi_Barrel_Lo", "Phi Barrel (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    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);
    
    //mPt_EndCap_Lo            = dbe->book1D("Pt_EndCap_Lo", "Pt EndCap (Pass Low Pt Jet Trigger)", 20, 0, 100);   
    //mPhi_EndCap_Lo           = dbe->book1D("Phi_EndCap_Lo", "Phi EndCap (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    mConstituents_EndCap     = dbe->book1D("Constituents_EndCap", "Constituents EndCap", 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);
    
    //mPt_Forward_Lo           = dbe->book1D("Pt_Forward_Lo", "Pt Forward (Pass Low Pt Jet Trigger)", 20, 0, 100);  
    //mPhi_Forward_Lo          = dbe->book1D("Phi_Forward_Lo", "Phi Forward (Pass Low Pt Jet Trigger)", phiBin, phiMin, phiMax);
    mConstituents_Forward    = dbe->book1D("Constituents_Forward", "Constituents Forward", 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);
    
    mPt_Barrel_Hi            = dbe->book1D("Pt_Barrel_Hi", "Pt Barrel (Pass Hi Pt Jet Trigger)", 60, 0, 300);   
    mPhi_Barrel_Hi           = dbe->book1D("Phi_Barrel_Hi", "Phi Barrel (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //mConstituents_Barrel_Hi  = dbe->book1D("Constituents_Barrel_Hi", "Constituents Barrel (Pass Hi Pt Jet Trigger)", 50, 0, 100);
    //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)", 60, 0, 300);  
    mPhi_EndCap_Hi           = dbe->book1D("Phi_EndCap_Hi", "Phi EndCap (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //mConstituents_EndCap_Hi  = dbe->book1D("Constituents_EndCap_Hi", "Constituents EndCap (Pass Hi Pt Jet Trigger)", 50, 0, 100);
    //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)", 60, 0, 300);  
    mPhi_Forward_Hi          = dbe->book1D("Phi_Forward_Hi", "Phi Forward (Pass Hi Pt Jet Trigger)", phiBin, phiMin, phiMax);
    //mConstituents_Forward_Hi = dbe->book1D("Constituents_Forward_Hi", "Constituents Forward (Pass Hi Pt Jet Trigger)", 50, 0, 100);
    //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);
    //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);
    //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);
    //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);
    //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.);

  //  msigmaEta                   = dbe->book1D("sigmaEta", "sigmaEta", 50, 0., 1.);
  //  msigmaPhi                   = dbe->book1D("sigmaPhi", "sigmaPhi", 50, 0., 0.5);
  
  if(makedijetselection==1) {
    mDijetAsymmetry                   = dbe->book1D("DijetAsymmetry", "DijetAsymmetry", 100, -1., 1.);
    mDijetBalance                     = dbe->book1D("DijetBalance",   "DijetBalance",   100, -2., 2.);
    if (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 239 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

Definition at line 91 of file JetAnalyzer.h.

Definition at line 92 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMax [private]

Definition at line 97 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMaxLoose [private]

Definition at line 102 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMaxTight [private]

Definition at line 105 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 99 of file JetAnalyzer.h.

Definition at line 104 of file JetAnalyzer.h.

Definition at line 107 of file JetAnalyzer.h.

double JetAnalyzer::_ptThreshold [private]

Definition at line 89 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMin [private]

Definition at line 98 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMinLoose [private]

Definition at line 103 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMinTight [private]

Definition at line 106 of file JetAnalyzer.h.

Definition at line 108 of file JetAnalyzer.h.

Definition at line 109 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 124 of file JetAnalyzer.h.

double JetAnalyzer::eMax [private]

Definition at line 126 of file JetAnalyzer.h.

double JetAnalyzer::eMin [private]

Definition at line 125 of file JetAnalyzer.h.

int JetAnalyzer::etaBin [private]

Definition at line 112 of file JetAnalyzer.h.

double JetAnalyzer::etaMax [private]

Definition at line 114 of file JetAnalyzer.h.

double JetAnalyzer::etaMin [private]

Definition at line 113 of file JetAnalyzer.h.

Definition at line 101 of file JetAnalyzer.h.

Definition at line 136 of file JetAnalyzer.h.

Definition at line 133 of file JetAnalyzer.h.

std::string JetAnalyzer::jetname [private]

Definition at line 80 of file JetAnalyzer.h.

Definition at line 94 of file JetAnalyzer.h.

Definition at line 158 of file JetAnalyzer.h.

Definition at line 174 of file JetAnalyzer.h.

Definition at line 190 of file JetAnalyzer.h.

Definition at line 179 of file JetAnalyzer.h.

Definition at line 194 of file JetAnalyzer.h.

Definition at line 184 of file JetAnalyzer.h.

Definition at line 198 of file JetAnalyzer.h.

Definition at line 254 of file JetAnalyzer.h.

Definition at line 253 of file JetAnalyzer.h.

Definition at line 211 of file JetAnalyzer.h.

Definition at line 160 of file JetAnalyzer.h.

Definition at line 176 of file JetAnalyzer.h.

Definition at line 181 of file JetAnalyzer.h.

Definition at line 186 of file JetAnalyzer.h.

Definition at line 227 of file JetAnalyzer.h.

Definition at line 228 of file JetAnalyzer.h.

Definition at line 229 of file JetAnalyzer.h.

Definition at line 156 of file JetAnalyzer.h.

Definition at line 248 of file JetAnalyzer.h.

Definition at line 214 of file JetAnalyzer.h.

Definition at line 233 of file JetAnalyzer.h.

Definition at line 234 of file JetAnalyzer.h.

Definition at line 224 of file JetAnalyzer.h.

Definition at line 226 of file JetAnalyzer.h.

Definition at line 225 of file JetAnalyzer.h.

Definition at line 223 of file JetAnalyzer.h.

Definition at line 159 of file JetAnalyzer.h.

Definition at line 175 of file JetAnalyzer.h.

Definition at line 191 of file JetAnalyzer.h.

Definition at line 180 of file JetAnalyzer.h.

Definition at line 195 of file JetAnalyzer.h.

Definition at line 185 of file JetAnalyzer.h.

Definition at line 199 of file JetAnalyzer.h.

Definition at line 238 of file JetAnalyzer.h.

Definition at line 239 of file JetAnalyzer.h.

Definition at line 221 of file JetAnalyzer.h.

Definition at line 222 of file JetAnalyzer.h.

Definition at line 232 of file JetAnalyzer.h.

Definition at line 210 of file JetAnalyzer.h.

Definition at line 157 of file JetAnalyzer.h.

Definition at line 164 of file JetAnalyzer.h.

Definition at line 189 of file JetAnalyzer.h.

Definition at line 167 of file JetAnalyzer.h.

Definition at line 193 of file JetAnalyzer.h.

Definition at line 170 of file JetAnalyzer.h.

Definition at line 197 of file JetAnalyzer.h.

Definition at line 249 of file JetAnalyzer.h.

Definition at line 245 of file JetAnalyzer.h.

Definition at line 215 of file JetAnalyzer.h.

Definition at line 161 of file JetAnalyzer.h.

Definition at line 152 of file JetAnalyzer.h.

Definition at line 153 of file JetAnalyzer.h.

Definition at line 154 of file JetAnalyzer.h.

Definition at line 155 of file JetAnalyzer.h.

Definition at line 163 of file JetAnalyzer.h.

Definition at line 188 of file JetAnalyzer.h.

Definition at line 166 of file JetAnalyzer.h.

Definition at line 192 of file JetAnalyzer.h.

Definition at line 169 of file JetAnalyzer.h.

Definition at line 196 of file JetAnalyzer.h.

Definition at line 250 of file JetAnalyzer.h.

Definition at line 246 of file JetAnalyzer.h.

Definition at line 217 of file JetAnalyzer.h.

Definition at line 235 of file JetAnalyzer.h.

Definition at line 240 of file JetAnalyzer.h.

Definition at line 241 of file JetAnalyzer.h.

Definition at line 78 of file JetAnalyzer.h.

int JetAnalyzer::pBin [private]

Definition at line 128 of file JetAnalyzer.h.

int JetAnalyzer::phiBin [private]

Definition at line 116 of file JetAnalyzer.h.

double JetAnalyzer::phiMax [private]

Definition at line 118 of file JetAnalyzer.h.

double JetAnalyzer::phiMin [private]

Definition at line 117 of file JetAnalyzer.h.

double JetAnalyzer::pMax [private]

Definition at line 130 of file JetAnalyzer.h.

double JetAnalyzer::pMin [private]

Definition at line 129 of file JetAnalyzer.h.

int JetAnalyzer::ptBin [private]

Definition at line 120 of file JetAnalyzer.h.

double JetAnalyzer::ptMax [private]

Definition at line 122 of file JetAnalyzer.h.

double JetAnalyzer::ptMin [private]

Definition at line 121 of file JetAnalyzer.h.

Definition at line 83 of file JetAnalyzer.h.