JetAnalyzer Class Reference

#include <JetAnalyzer.h>

Inheritance diagram for JetAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &, const reco::CaloJetCollection &caloJets, const int numPV)
	Get the analysis. More...
void	beginJob (DQMStore *dbe)
	Inizialize parameters for histo binning. More...
void	endJob ()
	Finish up a job. More...
int	getLeadJetFlag ()
	JetAnalyzer (const edm::ParameterSet &)
	Constructor. More...
void	setJetHiPass (int pass)
void	setJetLoPass (int pass)
void	setLeadJetFlag (int flag)
void	setSource (std::string source)
virtual	~JetAnalyzer ()
	Destructor. More...
Public Member Functions inherited from JetAnalyzerBase
void	analyze (const edm::Event &, const edm::EventSetup &, reco::CaloJet &caloJet)
	Get the analysis of the muon properties. More...
	JetAnalyzerBase ()
	Constructor. More...
virtual	~JetAnalyzerBase ()
	Destructor. More...

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::JetIDHelper *	jetID
MonitorElement *	jetME
std::string	jetname
int	makedijetselection
MonitorElement *	mConstituents
MonitorElement *	mConstituents_Barrel
MonitorElement *	mConstituents_Barrel_Hi
MonitorElement *	mConstituents_EndCap
MonitorElement *	mConstituents_EndCap_Hi
MonitorElement *	mConstituents_Forward
MonitorElement *	mConstituents_Forward_Hi
MonitorElement *	mConstituents_profile
MonitorElement *	mDijetAsymmetry
MonitorElement *	mDijetBalance
MonitorElement *	mDPhi
MonitorElement *	mEFrac
MonitorElement *	mEFrac_Barrel
MonitorElement *	mEFrac_EndCap
MonitorElement *	mEFrac_Forward
MonitorElement *	mEFrac_profile
MonitorElement *	mEmEnergyInEB
MonitorElement *	mEmEnergyInEE
MonitorElement *	mEmEnergyInHF
MonitorElement *	mEta
MonitorElement *	mEta_Hi
MonitorElement *	mEta_profile
MonitorElement *	mEtaFirst
MonitorElement *	mfHPD
MonitorElement *	mfRBX
MonitorElement *	mHadEnergyInHB
MonitorElement *	mHadEnergyInHE
MonitorElement *	mHadEnergyInHF
MonitorElement *	mHadEnergyInHO
MonitorElement *	mHFrac
MonitorElement *	mHFrac_Barrel
MonitorElement *	mHFrac_Barrel_Hi
MonitorElement *	mHFrac_EndCap
MonitorElement *	mHFrac_EndCap_Hi
MonitorElement *	mHFrac_Forward
MonitorElement *	mHFrac_Forward_Hi
MonitorElement *	mHFrac_profile
MonitorElement *	mLooseJIDPassFractionVSeta
MonitorElement *	mLooseJIDPassFractionVSpt
MonitorElement *	mMaxEInEmTowers
MonitorElement *	mMaxEInHadTowers
MonitorElement *	mN90Hits
MonitorElement *	mNJets
MonitorElement *	mNJets_profile
MonitorElement *	mPhi
MonitorElement *	mPhi_Barrel
MonitorElement *	mPhi_Barrel_Hi
MonitorElement *	mPhi_EndCap
MonitorElement *	mPhi_EndCap_Hi
MonitorElement *	mPhi_Forward
MonitorElement *	mPhi_Forward_Hi
MonitorElement *	mPhi_Hi
MonitorElement *	mPhi_Lo
MonitorElement *	mPhi_profile
MonitorElement *	mPhiFirst
MonitorElement *	mPhiVSEta
MonitorElement *	mPt
MonitorElement *	mPt_1
MonitorElement *	mPt_2
MonitorElement *	mPt_3
MonitorElement *	mPt_Barrel
MonitorElement *	mPt_Barrel_Hi
MonitorElement *	mPt_EndCap
MonitorElement *	mPt_EndCap_Hi
MonitorElement *	mPt_Forward
MonitorElement *	mPt_Forward_Hi
MonitorElement *	mPt_Hi
MonitorElement *	mPt_Lo
MonitorElement *	mPt_profile
MonitorElement *	mPtFirst
MonitorElement *	mresEMF
MonitorElement *	mTightJIDPassFractionVSeta
MonitorElement *	mTightJIDPassFractionVSpt
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

Author

F. Chlebana - Fermilab

Definition at line 39 of file JetAnalyzer.h.

Constructor & Destructor Documentation

JetAnalyzer::JetAnalyzer

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 20 of file JetAnalyzer.cc.

References Parameters::parameters.

                                                    {
  
  parameters   = pSet;
  _leadJetFlag = 0;
  _JetLoPass   = 0;
  _JetHiPass   = 0;
  _ptThreshold = 20.;
  _asymmetryThirdJetCut = 5.;
  _balanceThirdJetCut   = 0.2; 
  _n90HitsMin =0;
  _fHPDMax=1.;
  _resEMFMin=0.;
  _n90HitsMinLoose =0;
  _fHPDMaxLoose=1.;
  _resEMFMinLoose=0.;
  _n90HitsMinTight =0;
  _fHPDMaxTight=1.;
  _resEMFMinTight=0.;
  _sigmaEtaMinTight=-999.;
  _sigmaPhiMinTight=-999.;

} 

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,
		const int	numPV
	)

Get the analysis.

Definition at line 271 of file JetAnalyzer.cc.

References alignCSCRings::corr, diffTreeTool::diff, dPhi(), eta(), edm::EventID::event(), edm::EventBase::id(), metsig::jet, bTagSequences_cff::jetID, LogTrace, phi, RecoTauCleanerPlugins::pt, rand(), 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()));


          // Fill NPV profiles
          //----------------------------------------------------------------
          for (int ijet=0; ijet<2; ijet++) {

        if (mPt_profile)           mPt_profile          ->Fill(numPV, (caloJets.at(ijet)).pt());
        if (mEta_profile)          mEta_profile         ->Fill(numPV, (caloJets.at(ijet)).eta());
        if (mPhi_profile)          mPhi_profile         ->Fill(numPV, (caloJets.at(ijet)).phi());
        if (mConstituents_profile) mConstituents_profile->Fill(numPV, (caloJets.at(ijet)).nConstituents());
        if (mHFrac_profile)        mHFrac_profile       ->Fill(numPV, (caloJets.at(ijet)).energyFractionHadronic());
        if (mEFrac_profile)        mEFrac_profile       ->Fill(numPV, (caloJets.at(ijet)).emEnergyFraction());
          }
        }


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

      // Fill NPV profiles
      //--------------------------------------------------------------------
      if (mPt_profile)           mPt_profile          ->Fill(numPV, jet->pt());
      if (mEta_profile)          mEta_profile         ->Fill(numPV, jet->eta());
      if (mPhi_profile)          mPhi_profile         ->Fill(numPV, jet->phi());
      if (mConstituents_profile) mConstituents_profile->Fill(numPV, jet->nConstituents());
      if (mHFrac_profile)        mHFrac_profile       ->Fill(numPV, jet->energyFractionHadronic());
      if (mEFrac_profile)        mEFrac_profile       ->Fill(numPV, 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);


    if (mNJets_profile) mNJets_profile->Fill(numPV, numofjets);


  }//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, bTagSequences_cff::jetID, LogTrace, nbinsPV, Parameters::parameters, jptDQMConfig_cff::phiMax, jptDQMConfig_cff::phiMin, jptDQMConfig_cff::pMax, jptDQMConfig_cff::ptMax, PtMinSelector_cfg::ptMin, PVlow, PVup, MonitorElement::setAxisTitle(), 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);


  // Book NPV profiles
  //----------------------------------------------------------------------------
  mPt_profile           = dbe->bookProfile("Pt_profile",           "pt",                nbinsPV, PVlow, PVup,   ptBin,  ptMin,  ptMax);
  mEta_profile          = dbe->bookProfile("Eta_profile",          "eta",               nbinsPV, PVlow, PVup,  etaBin, etaMin, etaMax);
  mPhi_profile          = dbe->bookProfile("Phi_profile",          "phi",               nbinsPV, PVlow, PVup,  phiBin, phiMin, phiMax);
  mConstituents_profile = dbe->bookProfile("Constituents_profile", "# of constituents", nbinsPV, PVlow, PVup,      50,      0,    100);
  mHFrac_profile        = dbe->bookProfile("HFrac_profile",        "HFrac",             nbinsPV, PVlow, PVup,     120,   -0.1,    1.1);
  mEFrac_profile        = dbe->bookProfile("EFrac_profile",        "EFrac",             nbinsPV, PVlow, PVup,     120,   -0.1,    1.1);

  if (makedijetselection != 1)
    mNJets_profile = dbe->bookProfile("NJets_profile", "number of jets", nbinsPV, PVlow, PVup, 100, 0, 100);


  // Set NPV profiles x-axis title
  //----------------------------------------------------------------------------
  mPt_profile          ->setAxisTitle("nvtx",1);
  mEta_profile         ->setAxisTitle("nvtx",1);
  mPhi_profile         ->setAxisTitle("nvtx",1);
  mConstituents_profile->setAxisTitle("nvtx",1);
  mHFrac_profile       ->setAxisTitle("nvtx",1);
  mEFrac_profile       ->setAxisTitle("nvtx",1);

  if (makedijetselection != 1) {
    mNJets_profile->setAxisTitle("nvtx",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 265 of file JetAnalyzer.cc.

References bTagSequences_cff::jetID.

                         {
  delete jetID;
}

int JetAnalyzer::getLeadJetFlag ( )

inline

Definition at line 67 of file JetAnalyzer.h.

References _leadJetFlag.

                       {
    return  _leadJetFlag;
  }

void JetAnalyzer::setJetHiPass ( int  pass )

inline

Definition at line 74 of file JetAnalyzer.h.

References _JetHiPass.

                              {
    _JetHiPass = pass;
  }

void JetAnalyzer::setJetLoPass ( int  pass )

inline

Definition at line 70 of file JetAnalyzer.h.

References _JetLoPass.

                              {
    _JetLoPass = pass;
  }

void JetAnalyzer::setLeadJetFlag ( int  flag )

inline

Definition at line 64 of file JetAnalyzer.h.

References _leadJetFlag.

                                {
    _leadJetFlag = flag;
  }

void JetAnalyzer::setSource ( std::string  source )

inline

Definition at line 60 of file JetAnalyzer.h.

References _source, and source.

                                   {
    _source = source;
  }

Member Data Documentation

double JetAnalyzer::_asymmetryThirdJetCut

private

Definition at line 94 of file JetAnalyzer.h.

double JetAnalyzer::_balanceThirdJetCut

private

Definition at line 95 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMax

private

Definition at line 100 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMaxLoose

private

Definition at line 105 of file JetAnalyzer.h.

double JetAnalyzer::_fHPDMaxTight

private

Definition at line 108 of file JetAnalyzer.h.

int JetAnalyzer::_JetHiPass

private

Definition at line 89 of file JetAnalyzer.h.

Referenced by setJetHiPass().

int JetAnalyzer::_JetLoPass

private

Definition at line 88 of file JetAnalyzer.h.

Referenced by setJetLoPass().

int JetAnalyzer::_leadJetFlag

private

Definition at line 90 of file JetAnalyzer.h.

Referenced by getLeadJetFlag(), and setLeadJetFlag().

int JetAnalyzer::_n90HitsMin

private

Definition at line 102 of file JetAnalyzer.h.

int JetAnalyzer::_n90HitsMinLoose

private

Definition at line 107 of file JetAnalyzer.h.

int JetAnalyzer::_n90HitsMinTight

private

Definition at line 110 of file JetAnalyzer.h.

double JetAnalyzer::_ptThreshold

private

Definition at line 92 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMin

private

Definition at line 101 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMinLoose

private

Definition at line 106 of file JetAnalyzer.h.

double JetAnalyzer::_resEMFMinTight

private

Definition at line 109 of file JetAnalyzer.h.

double JetAnalyzer::_sigmaEtaMinTight

private

Definition at line 111 of file JetAnalyzer.h.

double JetAnalyzer::_sigmaPhiMinTight

private

Definition at line 112 of file JetAnalyzer.h.

std::string JetAnalyzer::_source

private

Definition at line 84 of file JetAnalyzer.h.

Referenced by setSource().

int JetAnalyzer::_theend

private

Definition at line 91 of file JetAnalyzer.h.

int JetAnalyzer::eBin

private

Definition at line 127 of file JetAnalyzer.h.

double JetAnalyzer::eMax

private

Definition at line 129 of file JetAnalyzer.h.

double JetAnalyzer::eMin

private

Definition at line 128 of file JetAnalyzer.h.

int JetAnalyzer::etaBin

private

Definition at line 115 of file JetAnalyzer.h.

double JetAnalyzer::etaMax

private

Definition at line 117 of file JetAnalyzer.h.

double JetAnalyzer::etaMin

private

Definition at line 116 of file JetAnalyzer.h.

int JetAnalyzer::fillJIDPassFrac

private

Definition at line 104 of file JetAnalyzer.h.

reco::helper::JetIDHelper* JetAnalyzer::jetID

private

Definition at line 139 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::jetME

private

Definition at line 136 of file JetAnalyzer.h.

std::string JetAnalyzer::jetname

private

Definition at line 83 of file JetAnalyzer.h.

int JetAnalyzer::makedijetselection

private

Definition at line 97 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents

private

Definition at line 161 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_Barrel

private

Definition at line 177 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_Barrel_Hi

private

Definition at line 193 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_EndCap

private

Definition at line 182 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_EndCap_Hi

private

Definition at line 197 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_Forward

private

Definition at line 187 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_Forward_Hi

private

Definition at line 201 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mConstituents_profile

private

Definition at line 266 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mDijetAsymmetry

private

Definition at line 257 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mDijetBalance

private

Definition at line 256 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mDPhi

private

Definition at line 214 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEFrac

private

Definition at line 163 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEFrac_Barrel

private

Definition at line 179 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEFrac_EndCap

private

Definition at line 184 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEFrac_Forward

private

Definition at line 189 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEFrac_profile

private

Definition at line 268 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEmEnergyInEB

private

Definition at line 230 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEmEnergyInEE

private

Definition at line 231 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEmEnergyInHF

private

Definition at line 232 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEta

private

Definition at line 159 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEta_Hi

private

Definition at line 251 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEta_profile

private

Definition at line 264 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mEtaFirst

private

Definition at line 217 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mfHPD

private

Definition at line 236 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mfRBX

private

Definition at line 237 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHadEnergyInHB

private

Definition at line 227 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHadEnergyInHE

private

Definition at line 229 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHadEnergyInHF

private

Definition at line 228 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHadEnergyInHO

private

Definition at line 226 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac

private

Definition at line 162 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_Barrel

private

Definition at line 178 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_Barrel_Hi

private

Definition at line 194 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_EndCap

private

Definition at line 183 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_EndCap_Hi

private

Definition at line 198 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_Forward

private

Definition at line 188 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_Forward_Hi

private

Definition at line 202 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mHFrac_profile

private

Definition at line 267 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mLooseJIDPassFractionVSeta

private

Definition at line 241 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mLooseJIDPassFractionVSpt

private

Definition at line 242 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mMaxEInEmTowers

private

Definition at line 224 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mMaxEInHadTowers

private

Definition at line 225 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mN90Hits

private

Definition at line 235 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mNJets

private

Definition at line 213 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mNJets_profile

private

Definition at line 262 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi

private

Definition at line 160 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_Barrel

private

Definition at line 167 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_Barrel_Hi

private

Definition at line 192 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_EndCap

private

Definition at line 170 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_EndCap_Hi

private

Definition at line 196 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_Forward

private

Definition at line 173 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_Forward_Hi

private

Definition at line 200 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_Hi

private

Definition at line 252 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_Lo

private

Definition at line 248 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhi_profile

private

Definition at line 265 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhiFirst

private

Definition at line 218 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPhiVSEta

private

Definition at line 164 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt

private

Definition at line 155 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_1

private

Definition at line 156 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_2

private

Definition at line 157 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_3

private

Definition at line 158 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_Barrel

private

Definition at line 166 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_Barrel_Hi

private

Definition at line 191 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_EndCap

private

Definition at line 169 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_EndCap_Hi

private

Definition at line 195 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_Forward

private

Definition at line 172 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_Forward_Hi

private

Definition at line 199 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_Hi

private

Definition at line 253 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_Lo

private

Definition at line 249 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPt_profile

private

Definition at line 263 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mPtFirst

private

Definition at line 220 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mresEMF

private

Definition at line 238 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mTightJIDPassFractionVSeta

private

Definition at line 243 of file JetAnalyzer.h.

MonitorElement* JetAnalyzer::mTightJIDPassFractionVSpt

private

Definition at line 244 of file JetAnalyzer.h.

edm::ParameterSet JetAnalyzer::parameters

private

Definition at line 81 of file JetAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

int JetAnalyzer::pBin

private

Definition at line 131 of file JetAnalyzer.h.

int JetAnalyzer::phiBin

private

Definition at line 119 of file JetAnalyzer.h.

double JetAnalyzer::phiMax

private

Definition at line 121 of file JetAnalyzer.h.

double JetAnalyzer::phiMin

private

Definition at line 120 of file JetAnalyzer.h.

double JetAnalyzer::pMax

private

Definition at line 133 of file JetAnalyzer.h.

double JetAnalyzer::pMin

private

Definition at line 132 of file JetAnalyzer.h.

int JetAnalyzer::ptBin

private

Definition at line 123 of file JetAnalyzer.h.

double JetAnalyzer::ptMax

private

Definition at line 125 of file JetAnalyzer.h.

double JetAnalyzer::ptMin

private

Definition at line 124 of file JetAnalyzer.h.

edm::InputTag JetAnalyzer::theCaloJetCollectionLabel

private

Definition at line 86 of file JetAnalyzer.h.