PFJetAnalyzer Class Reference

#include <PFJetAnalyzer.h>

Inheritance diagram for PFJetAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &, const reco::PFJetCollection &pfJets, 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 ()
	PFJetAnalyzer (const edm::ParameterSet &)
	Constructor. More...
void	setJetHiPass (int pass)
void	setJetLoPass (int pass)
void	setLeadJetFlag (int flag)
void	setSource (std::string source)
virtual	~PFJetAnalyzer ()
	Destructor. More...
Public Member Functions inherited from PFJetAnalyzerBase
void	analyze (const edm::Event &, const edm::EventSetup &, reco::PFJet &jet)
	Get the analysis of the muon properties. More...
	PFJetAnalyzerBase ()
	Constructor. More...
virtual	~PFJetAnalyzerBase ()
	Destructor. More...

Private Attributes
double	_asymmetryThirdJetCut
double	_balanceThirdJetCut
int	_JetHiPass
int	_JetLoPass
int	_leadJetFlag
double	_LooseCEFMax
double	_LooseCHFMin
double	_LooseNEFMax
double	_LooseNHFMax
double	_ptThreshold
std::string	_source
double	_ThisCEFMax
double	_ThisCHFMin
double	_ThisNEFMax
double	_ThisNHFMax
double	_TightCEFMax
double	_TightCHFMin
double	_TightNEFMax
double	_TightNHFMax
int	eBin
double	eMax
double	eMin
int	etaBin
double	etaMax
double	etaMin
int	fillpfJIDPassFrac
MonitorElement *	jetME
int	makedijetselection
MonitorElement *	mChargedEmEnergy
MonitorElement *	mChargedEmEnergy_profile
MonitorElement *	mChargedHadronEnergy
MonitorElement *	mChargedHadronEnergy_profile
MonitorElement *	mChargedMuEnergy
MonitorElement *	mChargedMuEnergy_profile
MonitorElement *	mChargedMultiplicity
MonitorElement *	mChargedMultiplicity_profile
MonitorElement *	mCHEn_highPt_Barrel
MonitorElement *	mCHEn_highPt_EndCap
MonitorElement *	mCHEn_lowPt_Barrel
MonitorElement *	mCHEn_lowPt_EndCap
MonitorElement *	mCHEn_mediumPt_Barrel
MonitorElement *	mCHEn_mediumPt_EndCap
MonitorElement *	mCHFrac_highPt_Barrel
MonitorElement *	mCHFrac_highPt_EndCap
MonitorElement *	mCHFrac_lowPt_Barrel
MonitorElement *	mCHFrac_lowPt_EndCap
MonitorElement *	mCHFrac_mediumPt_Barrel
MonitorElement *	mCHFrac_mediumPt_EndCap
MonitorElement *	mCHFracVSeta_highPt
MonitorElement *	mCHFracVSeta_lowPt
MonitorElement *	mCHFracVSeta_mediumPt
MonitorElement *	mCHFracVSpT_Barrel
MonitorElement *	mCHFracVSpT_EndCap
MonitorElement *	mChMultiplicity_highPt_Barrel
MonitorElement *	mChMultiplicity_highPt_EndCap
MonitorElement *	mChMultiplicity_highPt_Forward
MonitorElement *	mChMultiplicity_lowPt_Barrel
MonitorElement *	mChMultiplicity_lowPt_EndCap
MonitorElement *	mChMultiplicity_lowPt_Forward
MonitorElement *	mChMultiplicity_mediumPt_Barrel
MonitorElement *	mChMultiplicity_mediumPt_EndCap
MonitorElement *	mChMultiplicity_mediumPt_Forward
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 *	mEEffChargedFraction
MonitorElement *	mEEffNeutralFraction
MonitorElement *	mEFrac
MonitorElement *	mEFrac_Barrel
MonitorElement *	mEFrac_EndCap
MonitorElement *	mEFrac_Forward
MonitorElement *	mEFrac_profile
MonitorElement *	mElEn_highPt_Barrel
MonitorElement *	mElEn_highPt_EndCap
MonitorElement *	mElEn_lowPt_Barrel
MonitorElement *	mElEn_lowPt_EndCap
MonitorElement *	mElEn_mediumPt_Barrel
MonitorElement *	mElEn_mediumPt_EndCap
MonitorElement *	mElFrac_highPt_Barrel
MonitorElement *	mElFrac_highPt_EndCap
MonitorElement *	mElFrac_lowPt_Barrel
MonitorElement *	mElFrac_lowPt_EndCap
MonitorElement *	mElFrac_mediumPt_Barrel
MonitorElement *	mElFrac_mediumPt_EndCap
MonitorElement *	mElFracVSeta_highPt
MonitorElement *	mElFracVSeta_lowPt
MonitorElement *	mElFracVSeta_mediumPt
MonitorElement *	mElFracVSpT_Barrel
MonitorElement *	mElFracVSpT_EndCap
MonitorElement *	mEResChargedFraction
MonitorElement *	mEResNeutralFraction
MonitorElement *	mEta
MonitorElement *	mEta_Hi
MonitorElement *	mEta_profile
MonitorElement *	mEtaFirst
std::string	metname
MonitorElement *	mHFEEn_highPt_Forward
MonitorElement *	mHFEEn_lowPt_Forward
MonitorElement *	mHFEEn_mediumPt_Forward
MonitorElement *	mHFEFrac_highPt_Forward
MonitorElement *	mHFEFrac_lowPt_Forward
MonitorElement *	mHFEFrac_mediumPt_Forward
MonitorElement *	mHFEFracVSpT_Forward
MonitorElement *	mHFHEn_highPt_Forward
MonitorElement *	mHFHEn_lowPt_Forward
MonitorElement *	mHFHEn_mediumPt_Forward
MonitorElement *	mHFHFrac_highPt_Forward
MonitorElement *	mHFHFrac_lowPt_Forward
MonitorElement *	mHFHFrac_mediumPt_Forward
MonitorElement *	mHFHFracVSpT_Forward
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 *	mMuEn_highPt_Barrel
MonitorElement *	mMuEn_highPt_EndCap
MonitorElement *	mMuEn_lowPt_Barrel
MonitorElement *	mMuEn_lowPt_EndCap
MonitorElement *	mMuEn_mediumPt_Barrel
MonitorElement *	mMuEn_mediumPt_EndCap
MonitorElement *	mMuFrac_highPt_Barrel
MonitorElement *	mMuFrac_highPt_EndCap
MonitorElement *	mMuFrac_lowPt_Barrel
MonitorElement *	mMuFrac_lowPt_EndCap
MonitorElement *	mMuFrac_mediumPt_Barrel
MonitorElement *	mMuFrac_mediumPt_EndCap
MonitorElement *	mMuFracVSeta_highPt
MonitorElement *	mMuFracVSeta_lowPt
MonitorElement *	mMuFracVSeta_mediumPt
MonitorElement *	mMuFracVSpT_Barrel
MonitorElement *	mMuFracVSpT_EndCap
MonitorElement *	mMuMultiplicity_highPt_Barrel
MonitorElement *	mMuMultiplicity_highPt_EndCap
MonitorElement *	mMuMultiplicity_highPt_Forward
MonitorElement *	mMuMultiplicity_lowPt_Barrel
MonitorElement *	mMuMultiplicity_lowPt_EndCap
MonitorElement *	mMuMultiplicity_lowPt_Forward
MonitorElement *	mMuMultiplicity_mediumPt_Barrel
MonitorElement *	mMuMultiplicity_mediumPt_EndCap
MonitorElement *	mMuMultiplicity_mediumPt_Forward
MonitorElement *	mMuonMultiplicity
MonitorElement *	mMuonMultiplicity_profile
MonitorElement *	mNeuMultiplicity_highPt_Barrel
MonitorElement *	mNeuMultiplicity_highPt_EndCap
MonitorElement *	mNeuMultiplicity_highPt_Forward
MonitorElement *	mNeuMultiplicity_lowPt_Barrel
MonitorElement *	mNeuMultiplicity_lowPt_EndCap
MonitorElement *	mNeuMultiplicity_lowPt_Forward
MonitorElement *	mNeuMultiplicity_mediumPt_Barrel
MonitorElement *	mNeuMultiplicity_mediumPt_EndCap
MonitorElement *	mNeuMultiplicity_mediumPt_Forward
MonitorElement *	mNeutralEmEnergy
MonitorElement *	mNeutralEmEnergy_profile
MonitorElement *	mNeutralFraction
MonitorElement *	mNeutralFraction2
MonitorElement *	mNeutralHadronEnergy
MonitorElement *	mNeutralHadronEnergy_profile
MonitorElement *	mNeutralMultiplicity
MonitorElement *	mNeutralMultiplicity_profile
MonitorElement *	mNHEn_highPt_Barrel
MonitorElement *	mNHEn_highPt_EndCap
MonitorElement *	mNHEn_lowPt_Barrel
MonitorElement *	mNHEn_lowPt_EndCap
MonitorElement *	mNHEn_mediumPt_Barrel
MonitorElement *	mNHEn_mediumPt_EndCap
MonitorElement *	mNHFrac_highPt_Barrel
MonitorElement *	mNHFrac_highPt_EndCap
MonitorElement *	mNHFrac_lowPt_Barrel
MonitorElement *	mNHFrac_lowPt_EndCap
MonitorElement *	mNHFrac_mediumPt_Barrel
MonitorElement *	mNHFrac_mediumPt_EndCap
MonitorElement *	mNHFracVSeta_highPt
MonitorElement *	mNHFracVSeta_lowPt
MonitorElement *	mNHFracVSeta_mediumPt
MonitorElement *	mNHFracVSpT_Barrel
MonitorElement *	mNHFracVSpT_EndCap
MonitorElement *	mNJets
MonitorElement *	mNJets_profile
MonitorElement *	mPhEn_highPt_Barrel
MonitorElement *	mPhEn_highPt_EndCap
MonitorElement *	mPhEn_lowPt_Barrel
MonitorElement *	mPhEn_lowPt_EndCap
MonitorElement *	mPhEn_mediumPt_Barrel
MonitorElement *	mPhEn_mediumPt_EndCap
MonitorElement *	mPhFrac_highPt_Barrel
MonitorElement *	mPhFrac_highPt_EndCap
MonitorElement *	mPhFrac_lowPt_Barrel
MonitorElement *	mPhFrac_lowPt_EndCap
MonitorElement *	mPhFrac_mediumPt_Barrel
MonitorElement *	mPhFrac_mediumPt_EndCap
MonitorElement *	mPhFracVSeta_highPt
MonitorElement *	mPhFracVSeta_lowPt
MonitorElement *	mPhFracVSeta_mediumPt
MonitorElement *	mPhFracVSpT_Barrel
MonitorElement *	mPhFracVSpT_EndCap
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 *	mTightJIDPassFractionVSeta
MonitorElement *	mTightJIDPassFractionVSpt
MonitorElement *	nEEff
edm::ParameterSet	parameters
int	pBin
int	phiBin
double	phiMax
double	phiMin
double	pMax
double	pMin
int	ptBin
double	ptMax
double	ptMin
edm::InputTag	thePFJetCollectionLabel

Detailed Description

DQM monitoring source for PFlow Jets

Date:

2012/03/23 18:24:43

Revision:

1.11

Author

F. Chlebana - Fermilab

Definition at line 32 of file PFJetAnalyzer.h.

Constructor & Destructor Documentation

PFJetAnalyzer::PFJetAnalyzer

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 20 of file PFJetAnalyzer.cc.

References Parameters::parameters.

                                                        {

  parameters = pSet;
  _leadJetFlag = 0;
  _JetLoPass   = 0;
  _JetHiPass   = 0;
  _ptThreshold = 5.;
  _asymmetryThirdJetCut = 30.;
  _balanceThirdJetCut   = 0.2;
  _LooseCHFMin = -999.;
  _LooseNHFMax = -999.;
  _LooseCEFMax = -999.;
  _LooseNEFMax = -999.;
  _TightCHFMin = -999.;
  _TightNHFMax = -999.;
  _TightCEFMax = -999.;
  _TightNEFMax = -999.;
  _ThisCHFMin = -999.;
  _ThisNHFMax = -999.;
  _ThisCEFMax = -999.;
  _ThisNEFMax = -999.;

}

PFJetAnalyzer::~PFJetAnalyzer ( )

virtual

Destructor.

Definition at line 45 of file PFJetAnalyzer.cc.

{ }

Member Function Documentation

void PFJetAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const reco::PFJetCollection &	pfJets,
		const int	numPV
	)

Get the analysis.

Definition at line 423 of file PFJetAnalyzer.cc.

References corr, diffTreeTool::diff, dPhi(), eta(), edm::EventID::event(), edm::EventBase::id(), metsig::jet, LogTrace, metname, phi, and rand().

{
  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 ThisCHFcleaned=false;
  bool LooseCHFcleaned=false;
  bool TightCHFcleaned=false;

  srand( iEvent.id().event() % 10000);


  if (makedijetselection == 1) {
    //Dijet selection - careful: the pT is uncorrected!
    //if(makedijetselection==1 && pfJets.size()>=2){
    if(pfJets.size()>=2){
      double  dphiDJ = -999. ;

      bool LoosecleanedFirstJet =false; 
      bool LoosecleanedSecondJet=false; 
      bool TightcleanedFirstJet =false; 
      bool TightcleanedSecondJet=false; 
      bool LooseCHFcleanedFirstJet =false;
      bool LooseCHFcleanedSecondJet=false;
      bool TightCHFcleanedFirstJet =false;
      bool TightCHFcleanedSecondJet=false;

      //both jets pass pt threshold
      if ((pfJets.at(0)).pt() > _ptThreshold && (pfJets.at(1)).pt() > _ptThreshold ) {
    if(fabs((pfJets.at(0)).eta())<3. && fabs((pfJets.at(1)).eta())<3. ){
      //calculate dphi
      dphiDJ = fabs((pfJets.at(0)).phi()-(pfJets.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){
        //first jet
        LooseCHFcleanedFirstJet=true;
        TightCHFcleanedFirstJet=true;
        if((pfJets.at(0).chargedHadronEnergy()/pfJets.at(0).energy())<=_LooseCHFMin && fabs(pfJets.at(0).eta())<2.4) LooseCHFcleanedFirstJet=false; //apply CHF>0 only if |eta|<2.4
        if((pfJets.at(0).chargedHadronEnergy()/pfJets.at(0).energy())<=_TightCHFMin && fabs(pfJets.at(0).eta())<2.4) TightCHFcleanedFirstJet=false; //apply CHF>0 only if |eta|<2.4
        if(LooseCHFcleanedFirstJet && (pfJets.at(0).neutralHadronEnergy()/pfJets.at(0).energy())<_LooseNHFMax && (pfJets.at(0).chargedEmEnergy()/pfJets.at(0).energy())<_LooseCEFMax && (pfJets.at(0).neutralEmEnergy()/pfJets.at(0).energy())<_LooseNEFMax) LoosecleanedFirstJet=true;
        if(TightCHFcleanedFirstJet && (pfJets.at(0).neutralHadronEnergy()/pfJets.at(0).energy())<_TightNHFMax && (pfJets.at(0).chargedEmEnergy()/pfJets.at(0).energy())<_TightCEFMax && (pfJets.at(0).neutralEmEnergy()/pfJets.at(0).energy())<_TightNEFMax) TightcleanedFirstJet=true;

        //second jet
        LooseCHFcleanedSecondJet=true;
        TightCHFcleanedSecondJet=true;
        if((pfJets.at(1).chargedHadronEnergy()/pfJets.at(1).energy())<=_LooseCHFMin && fabs(pfJets.at(1).eta())<2.4) LooseCHFcleanedSecondJet=false; //apply CHF>0 only if |eta|<2.4
        if((pfJets.at(1).chargedHadronEnergy()/pfJets.at(1).energy())<=_TightCHFMin && fabs(pfJets.at(1).eta())<2.4) TightCHFcleanedSecondJet=false; //apply CHF>0 only if |eta|<2.4
        if(LooseCHFcleanedSecondJet && (pfJets.at(1).neutralHadronEnergy()/pfJets.at(1).energy())<_LooseNHFMax && (pfJets.at(1).chargedEmEnergy()/pfJets.at(1).energy())<_LooseCEFMax && (pfJets.at(1).neutralEmEnergy()/pfJets.at(1).energy())<_LooseNEFMax) LoosecleanedSecondJet=true;
        if(TightCHFcleanedSecondJet && (pfJets.at(1).neutralHadronEnergy()/pfJets.at(1).energy())<_TightNHFMax && (pfJets.at(1).chargedEmEnergy()/pfJets.at(1).energy())<_TightCEFMax && (pfJets.at(1).neutralEmEnergy()/pfJets.at(1).energy())<_TightNEFMax) TightcleanedSecondJet=true;
      
        if(fillpfJIDPassFrac==1) {
          //fill the profile for jid efficiency 
          if(LoosecleanedFirstJet) {
        mLooseJIDPassFractionVSeta->Fill(pfJets.at(0).eta(),1.);
        mLooseJIDPassFractionVSpt->Fill(pfJets.at(0).pt(),1.);
          } else {
        mLooseJIDPassFractionVSeta->Fill(pfJets.at(0).eta(),0.);
        mLooseJIDPassFractionVSpt->Fill(pfJets.at(0).pt(),0.);
          }
          if(TightcleanedFirstJet) {
        mTightJIDPassFractionVSeta->Fill(pfJets.at(0).eta(),1.);
        mTightJIDPassFractionVSpt->Fill(pfJets.at(0).pt(),1.);
          } else {
        mTightJIDPassFractionVSeta->Fill(pfJets.at(0).eta(),0.);
        mTightJIDPassFractionVSpt->Fill(pfJets.at(0).pt(),0.);
          }

          if(LoosecleanedSecondJet) {
        mLooseJIDPassFractionVSeta->Fill(pfJets.at(1).eta(),1.);
        mLooseJIDPassFractionVSpt->Fill(pfJets.at(1).pt(),1.);
          } else {
        mLooseJIDPassFractionVSeta->Fill(pfJets.at(1).eta(),0.);
        mLooseJIDPassFractionVSpt->Fill(pfJets.at(1).pt(),0.);
          }
          if(TightcleanedSecondJet) {
        mTightJIDPassFractionVSeta->Fill(pfJets.at(1).eta(),1.);
        mTightJIDPassFractionVSpt->Fill(pfJets.at(1).pt(),1.);
          } else {
        mTightJIDPassFractionVSeta->Fill(pfJets.at(1).eta(),0.);
        mTightJIDPassFractionVSpt->Fill(pfJets.at(1).pt(),0.);
          }
        }
        
        if(LoosecleanedFirstJet && LoosecleanedSecondJet) {
          //Filling variables for first jet
          if (mPt)   mPt->Fill (pfJets.at(0).pt());
          if (mEta)  mEta->Fill (pfJets.at(0).eta());
          if (mPhi)  mPhi->Fill (pfJets.at(0).phi());
          if (mPhiVSEta) mPhiVSEta->Fill(pfJets.at(0).eta(),pfJets.at(0).phi());
          
          if (mConstituents) mConstituents->Fill (pfJets.at(0).nConstituents());
          if (mHFrac)        mHFrac->Fill (pfJets.at(0).chargedHadronEnergyFraction()+pfJets.at(0).neutralHadronEnergyFraction());
          if (mEFrac)        mEFrac->Fill (pfJets.at(0).chargedEmEnergyFraction() +pfJets.at(0).neutralEmEnergyFraction());
          
          //if (mE) mE->Fill (pfJets.at(0).energy());
          //if (mP) mP->Fill (pfJets.at(0).p());
          //if (mMass) mMass->Fill (pfJets.at(0).mass());
            
          if (mChargedHadronEnergy)  mChargedHadronEnergy->Fill (pfJets.at(0).chargedHadronEnergy());
          if (mNeutralHadronEnergy)  mNeutralHadronEnergy->Fill (pfJets.at(0).neutralHadronEnergy());
          if (mChargedEmEnergy) mChargedEmEnergy->Fill(pfJets.at(0).chargedEmEnergy());
          if (mChargedMuEnergy) mChargedMuEnergy->Fill (pfJets.at(0).chargedMuEnergy ());
          if (mNeutralEmEnergy) mNeutralEmEnergy->Fill(pfJets.at(0).neutralEmEnergy());
          if (mChargedMultiplicity ) mChargedMultiplicity->Fill(pfJets.at(0).chargedMultiplicity());
          if (mNeutralMultiplicity ) mNeutralMultiplicity->Fill(pfJets.at(0).neutralMultiplicity());
          if (mMuonMultiplicity )mMuonMultiplicity->Fill (pfJets.at(0). muonMultiplicity());
          //_______________________________________________________
          if (mNeutralFraction) mNeutralFraction->Fill (pfJets.at(0).neutralMultiplicity()/pfJets.at(0).nConstituents());


          //Filling variables for second jet
          if (mPt)   mPt->Fill (pfJets.at(1).pt());
          if (mEta)  mEta->Fill (pfJets.at(1).eta());
          if (mPhi)  mPhi->Fill (pfJets.at(1).phi());
          if (mPhiVSEta) mPhiVSEta->Fill(pfJets.at(1).eta(),pfJets.at(1).phi());
          
          if (mConstituents) mConstituents->Fill (pfJets.at(1).nConstituents());
          if (mHFrac)        mHFrac->Fill (pfJets.at(1).chargedHadronEnergyFraction()+pfJets.at(1).neutralHadronEnergyFraction());
          if (mEFrac)        mEFrac->Fill (pfJets.at(1).chargedEmEnergyFraction() +pfJets.at(1).neutralEmEnergyFraction());
          
          //if (mE) mE->Fill (pfJets.at(1).energy());
          //if (mP) mP->Fill (pfJets.at(1).p());
          //if (mMass) mMass->Fill (pfJets.at(1).mass());
            
          if (mChargedHadronEnergy)  mChargedHadronEnergy->Fill (pfJets.at(1).chargedHadronEnergy());
          if (mNeutralHadronEnergy)  mNeutralHadronEnergy->Fill (pfJets.at(1).neutralHadronEnergy());
          if (mChargedEmEnergy) mChargedEmEnergy->Fill(pfJets.at(1).chargedEmEnergy());
          if (mChargedMuEnergy) mChargedMuEnergy->Fill (pfJets.at(1).chargedMuEnergy ());
          if (mNeutralEmEnergy) mNeutralEmEnergy->Fill(pfJets.at(1).neutralEmEnergy());
          if (mChargedMultiplicity ) mChargedMultiplicity->Fill(pfJets.at(1).chargedMultiplicity());
          if (mNeutralMultiplicity ) mNeutralMultiplicity->Fill(pfJets.at(1).neutralMultiplicity());
          if (mMuonMultiplicity )mMuonMultiplicity->Fill (pfJets.at(1). muonMultiplicity());
          //_______________________________________________________
          if (mNeutralFraction) mNeutralFraction->Fill (pfJets.at(1).neutralMultiplicity()/pfJets.at(1).nConstituents());


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

        if (mPt_profile)           mPt_profile          ->Fill(numPV, pfJets.at(iJet).pt());
        if (mEta_profile)          mEta_profile         ->Fill(numPV, pfJets.at(iJet).eta());
        if (mPhi_profile)          mPhi_profile         ->Fill(numPV, pfJets.at(iJet).phi());
        if (mConstituents_profile) mConstituents_profile->Fill(numPV, pfJets.at(iJet).nConstituents());
        if (mHFrac_profile)        mHFrac_profile       ->Fill(numPV, pfJets.at(iJet).chargedHadronEnergyFraction() + pfJets.at(iJet).neutralHadronEnergyFraction());
        if (mEFrac_profile)        mEFrac_profile       ->Fill(numPV, pfJets.at(iJet).chargedEmEnergyFraction()     + pfJets.at(iJet).neutralEmEnergyFraction());
          
        if (mChargedHadronEnergy_profile) mChargedHadronEnergy_profile->Fill(numPV, pfJets.at(iJet).chargedHadronEnergy());
        if (mNeutralHadronEnergy_profile) mNeutralHadronEnergy_profile->Fill(numPV, pfJets.at(iJet).neutralHadronEnergy());
        if (mChargedEmEnergy_profile)     mChargedEmEnergy_profile    ->Fill(numPV, pfJets.at(iJet).chargedEmEnergy());
        if (mChargedMuEnergy_profile)     mChargedMuEnergy_profile    ->Fill(numPV, pfJets.at(iJet).chargedMuEnergy());
        if (mNeutralEmEnergy_profile)     mNeutralEmEnergy_profile    ->Fill(numPV, pfJets.at(iJet).neutralEmEnergy());
        if (mChargedMultiplicity_profile) mChargedMultiplicity_profile->Fill(numPV, pfJets.at(iJet).chargedMultiplicity());
        if (mNeutralMultiplicity_profile) mNeutralMultiplicity_profile->Fill(numPV, pfJets.at(iJet).neutralMultiplicity());
        if (mMuonMultiplicity_profile)    mMuonMultiplicity_profile   ->Fill(numPV, pfJets.at(iJet).muonMultiplicity());
          }


        }// loose cleaned jets 1 and 2
      }// fabs dphi < 2.1
    }// fabs eta < 3
      }// pt jets > threshold
      //now do the dijet balance and asymmetry calculations
      if (fabs(pfJets.at(0).eta() < 1.4)) {
    double pt_dijet = (pfJets.at(0).pt() + pfJets.at(1).pt())/2;

    double dPhi = fabs((pfJets.at(0)).phi()-(pfJets.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(pfJets.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 < pfJets.size(); ++third) 
          if (pfJets.at(third).pt() > _asymmetryThirdJetCut) 
        thirdJetCut = false;
        if (thirdJetCut) {
          double dijetAsymmetry = (pfJets.at(jet1).pt() - pfJets.at(jet2).pt()) / (pfJets.at(jet1).pt() + pfJets.at(jet2).pt());
          mDijetAsymmetry->Fill(dijetAsymmetry);
        }// end restriction on third jet pt in asymmetry calculation
          
      }
      else {
        jet1 = 0;
        jet2 = 1;
      }
      
      pt_barrel = pfJets.at(jet1).pt();
      pt_probe  = pfJets.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 < pfJets.size(); ++third) 
        if (pfJets.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::PFJetCollection::const_iterator jet = pfJets.begin(); jet!=pfJets.end(); ++jet){
      LogTrace(metname)<<"[JetAnalyzer] Analyze PFJet";
 
      Thiscleaned=false;
      Loosecleaned=false;
      Tightcleaned=false;

      if (jet == pfJets.begin()) {
    fstPhi = jet->phi();
    _leadJetFlag = 1;
      } else {
    _leadJetFlag = 0;
      }

      if (jet == (pfJets.begin()+1)) sndPhi = jet->phi();
      //  if (jet->pt() < _ptThreshold) return;
      if (jet->pt() > _ptThreshold) {
    numofjets++ ;
    jetME->Fill(2);
            
    //calculate the jetID
    ThisCHFcleaned=true;
    LooseCHFcleaned=true;
    TightCHFcleaned=true;
    if((jet->chargedHadronEnergy()/jet->energy())<=_ThisCHFMin && fabs(jet->eta())<2.4) ThisCHFcleaned=false; //apply CHF>0 only if |eta|<2.4
    if((jet->chargedHadronEnergy()/jet->energy())<=_LooseCHFMin && fabs(jet->eta())<2.4) LooseCHFcleaned=false; //apply CHF>0 only if |eta|<2.4
    if((jet->chargedHadronEnergy()/jet->energy())<=_TightCHFMin && fabs(jet->eta())<2.4) TightCHFcleaned=false; //apply CHF>0 only if |eta|<2.4
    if(ThisCHFcleaned && (jet->neutralHadronEnergy()/jet->energy())<_ThisNHFMax && (jet->chargedEmEnergy()/jet->energy())<_ThisCEFMax && (jet->neutralEmEnergy()/jet->energy())<_ThisNEFMax) Thiscleaned=true;
    if(LooseCHFcleaned && (jet->neutralHadronEnergy()/jet->energy())<_LooseNHFMax && (jet->chargedEmEnergy()/jet->energy())<_LooseCEFMax && (jet->neutralEmEnergy()/jet->energy())<_LooseNEFMax) Loosecleaned=true;
    if(TightCHFcleaned && (jet->neutralHadronEnergy()/jet->energy())<_TightNHFMax && (jet->chargedEmEnergy()/jet->energy())<_TightCEFMax && (jet->neutralEmEnergy()/jet->energy())<_TightNEFMax) Tightcleaned=true;
      
    if(fillpfJIDPassFrac==1) {
      //fill the profile for jid efficiency 
      if(Loosecleaned) {
        mLooseJIDPassFractionVSeta->Fill(jet->eta(),1.);
        mLooseJIDPassFractionVSpt->Fill(jet->pt(),1.);
      } else {
        mLooseJIDPassFractionVSeta->Fill(jet->eta(),0.);
        mLooseJIDPassFractionVSpt->Fill(jet->pt(),0.);
      }
      if(Tightcleaned) {
        mTightJIDPassFractionVSeta->Fill(jet->eta(),1.);
        mTightJIDPassFractionVSpt->Fill(jet->pt(),1.);
      } else {
        mTightJIDPassFractionVSeta->Fill(jet->eta(),0.);
        mTightJIDPassFractionVSpt->Fill(jet->pt(),0.);
      }
    }
      
    if(!Thiscleaned) continue;
      
    // 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 (no longer used)
    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->chargedHadronEnergyFraction()+jet->neutralHadronEnergyFraction());   
      }
      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->chargedHadronEnergyFraction()+jet->neutralHadronEnergyFraction());   
      }
      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->chargedHadronEnergyFraction()+jet->neutralHadronEnergyFraction()); 
      }
    
      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->chargedHadronEnergyFraction()+jet->neutralHadronEnergyFraction());
    if (mEFrac)        mEFrac->Fill (jet->chargedEmEnergyFraction() +jet->neutralEmEnergyFraction());
      
    if (jet->pt()<= 50) {
      if (mCHFracVSeta_lowPt) mCHFracVSeta_lowPt->Fill(jet->eta(),jet->chargedHadronEnergyFraction());
      if (mNHFracVSeta_lowPt) mNHFracVSeta_lowPt->Fill(jet->eta(),jet->neutralHadronEnergyFraction());
      if (mPhFracVSeta_lowPt) mPhFracVSeta_lowPt->Fill(jet->eta(),jet->neutralEmEnergyFraction());
      if (mElFracVSeta_lowPt) mElFracVSeta_lowPt->Fill(jet->eta(),jet->chargedEmEnergyFraction());
      if (mMuFracVSeta_lowPt) mMuFracVSeta_lowPt->Fill(jet->eta(),jet->chargedMuEnergyFraction());
    }
    if (jet->pt()>50. && jet->pt()<=140.) {
      if (mCHFracVSeta_mediumPt) mCHFracVSeta_mediumPt->Fill(jet->eta(),jet->chargedHadronEnergyFraction());
      if (mNHFracVSeta_mediumPt) mNHFracVSeta_mediumPt->Fill(jet->eta(),jet->neutralHadronEnergyFraction());
      if (mPhFracVSeta_mediumPt) mPhFracVSeta_mediumPt->Fill(jet->eta(),jet->neutralEmEnergyFraction());
      if (mElFracVSeta_mediumPt) mElFracVSeta_mediumPt->Fill(jet->eta(),jet->chargedEmEnergyFraction());
      if (mMuFracVSeta_mediumPt) mMuFracVSeta_mediumPt->Fill(jet->eta(),jet->chargedMuEnergyFraction());
    }
    if (jet->pt()>140.) {
      if (mCHFracVSeta_highPt) mCHFracVSeta_highPt->Fill(jet->eta(),jet->chargedHadronEnergyFraction());
      if (mNHFracVSeta_highPt) mNHFracVSeta_highPt->Fill(jet->eta(),jet->neutralHadronEnergyFraction());
      if (mPhFracVSeta_highPt) mPhFracVSeta_highPt->Fill(jet->eta(),jet->neutralEmEnergyFraction());
      if (mElFracVSeta_highPt) mElFracVSeta_highPt->Fill(jet->eta(),jet->chargedEmEnergyFraction());
      if (mMuFracVSeta_highPt) mMuFracVSeta_highPt->Fill(jet->eta(),jet->chargedMuEnergyFraction());
    }

    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->chargedHadronEnergyFraction() + jet->neutralHadronEnergyFraction() );
    if (mEFrac_Barrel)           mEFrac->Fill (jet->chargedEmEnergyFraction() + jet->neutralEmEnergyFraction());    
      //fractions
      if (jet->pt()<=50.) {
        if (mCHFrac_lowPt_Barrel) mCHFrac_lowPt_Barrel->Fill(jet->chargedHadronEnergyFraction());
        if (mNHFrac_lowPt_Barrel) mNHFrac_lowPt_Barrel->Fill(jet->neutralHadronEnergyFraction());
        if (mPhFrac_lowPt_Barrel) mPhFrac_lowPt_Barrel->Fill(jet->neutralEmEnergyFraction());
        if (mElFrac_lowPt_Barrel) mElFrac_lowPt_Barrel->Fill(jet->chargedEmEnergyFraction());
        if (mMuFrac_lowPt_Barrel) mMuFrac_lowPt_Barrel->Fill(jet->chargedMuEnergyFraction());
        //
        if (mCHEn_lowPt_Barrel) mCHEn_lowPt_Barrel->Fill(jet->chargedHadronEnergy());
        if (mNHEn_lowPt_Barrel) mNHEn_lowPt_Barrel->Fill(jet->neutralHadronEnergy());
        if (mPhEn_lowPt_Barrel) mPhEn_lowPt_Barrel->Fill(jet->neutralEmEnergy());
        if (mElEn_lowPt_Barrel) mElEn_lowPt_Barrel->Fill(jet->chargedEmEnergy());
        if (mMuEn_lowPt_Barrel) mMuEn_lowPt_Barrel->Fill(jet->chargedMuEnergy());
      }
      if (jet->pt()>50. && jet->pt()<=140.) {
        if (mCHFrac_mediumPt_Barrel) mCHFrac_mediumPt_Barrel->Fill(jet->chargedHadronEnergyFraction());
        if (mNHFrac_mediumPt_Barrel) mNHFrac_mediumPt_Barrel->Fill(jet->neutralHadronEnergyFraction());
        if (mPhFrac_mediumPt_Barrel) mPhFrac_mediumPt_Barrel->Fill(jet->neutralEmEnergyFraction());
        if (mElFrac_mediumPt_Barrel) mElFrac_mediumPt_Barrel->Fill(jet->chargedEmEnergyFraction());
        if (mMuFrac_mediumPt_Barrel) mMuFrac_mediumPt_Barrel->Fill(jet->chargedMuEnergyFraction());
        //
        if (mCHEn_mediumPt_Barrel) mCHEn_mediumPt_Barrel->Fill(jet->chargedHadronEnergy());
        if (mNHEn_mediumPt_Barrel) mNHEn_mediumPt_Barrel->Fill(jet->neutralHadronEnergy());
        if (mPhEn_mediumPt_Barrel) mPhEn_mediumPt_Barrel->Fill(jet->neutralEmEnergy());
        if (mElEn_mediumPt_Barrel) mElEn_mediumPt_Barrel->Fill(jet->chargedEmEnergy());
        if (mMuEn_mediumPt_Barrel) mMuEn_mediumPt_Barrel->Fill(jet->chargedMuEnergy());
      }
      if (jet->pt()>140.) {
        if (mCHFrac_highPt_Barrel) mCHFrac_highPt_Barrel->Fill(jet->chargedHadronEnergyFraction());
        if (mNHFrac_highPt_Barrel) mNHFrac_highPt_Barrel->Fill(jet->neutralHadronEnergyFraction());
        if (mPhFrac_highPt_Barrel) mPhFrac_highPt_Barrel->Fill(jet->neutralEmEnergyFraction());
        if (mElFrac_highPt_Barrel) mElFrac_highPt_Barrel->Fill(jet->chargedEmEnergyFraction());
        if (mMuFrac_highPt_Barrel) mMuFrac_highPt_Barrel->Fill(jet->chargedMuEnergyFraction());
        //
        if (mCHEn_highPt_Barrel) mCHEn_highPt_Barrel->Fill(jet->chargedHadronEnergy());
        if (mNHEn_highPt_Barrel) mNHEn_highPt_Barrel->Fill(jet->neutralHadronEnergy());
        if (mPhEn_highPt_Barrel) mPhEn_highPt_Barrel->Fill(jet->neutralEmEnergy());
        if (mElEn_highPt_Barrel) mElEn_highPt_Barrel->Fill(jet->chargedEmEnergy());
        if (mMuEn_highPt_Barrel) mMuEn_highPt_Barrel->Fill(jet->chargedMuEnergy());
      }
      if(mChMultiplicity_lowPt_Barrel)  mChMultiplicity_lowPt_Barrel->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_lowPt_Barrel)  mNeuMultiplicity_lowPt_Barrel->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_lowPt_Barrel)  mMuMultiplicity_lowPt_Barrel->Fill(jet->muonMultiplicity());
      if(mChMultiplicity_mediumPt_Barrel)  mChMultiplicity_mediumPt_Barrel->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_mediumPt_Barrel)  mNeuMultiplicity_mediumPt_Barrel->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_mediumPt_Barrel)  mMuMultiplicity_mediumPt_Barrel->Fill(jet->muonMultiplicity());
      if(mChMultiplicity_highPt_Barrel)  mChMultiplicity_highPt_Barrel->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_highPt_Barrel)  mNeuMultiplicity_highPt_Barrel->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_highPt_Barrel)  mMuMultiplicity_highPt_Barrel->Fill(jet->muonMultiplicity());
      //
      if (mCHFracVSpT_Barrel) mCHFracVSpT_Barrel->Fill(jet->pt(),jet->chargedHadronEnergyFraction());
      if (mNHFracVSpT_Barrel) mNHFracVSpT_Barrel->Fill(jet->pt(),jet->neutralHadronEnergyFraction());
      if (mPhFracVSpT_Barrel) mPhFracVSpT_Barrel->Fill(jet->pt(),jet->neutralEmEnergyFraction());
      if (mElFracVSpT_Barrel) mElFracVSpT_Barrel->Fill(jet->pt(),jet->chargedEmEnergyFraction());
      if (mMuFracVSpT_Barrel) mMuFracVSpT_Barrel->Fill(jet->pt(),jet->chargedMuEnergyFraction());
    }
    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->chargedHadronEnergyFraction() + jet->neutralHadronEnergyFraction());
    if (mEFrac_EndCap)           mEFrac->Fill (jet->chargedEmEnergyFraction() + jet->neutralEmEnergyFraction());
      //fractions
      if (jet->pt()<=50.) {
        if (mCHFrac_lowPt_EndCap) mCHFrac_lowPt_EndCap->Fill(jet->chargedHadronEnergyFraction());
        if (mNHFrac_lowPt_EndCap) mNHFrac_lowPt_EndCap->Fill(jet->neutralHadronEnergyFraction());
        if (mPhFrac_lowPt_EndCap) mPhFrac_lowPt_EndCap->Fill(jet->neutralEmEnergyFraction());
        if (mElFrac_lowPt_EndCap) mElFrac_lowPt_EndCap->Fill(jet->chargedEmEnergyFraction());
        if (mMuFrac_lowPt_EndCap) mMuFrac_lowPt_EndCap->Fill(jet->chargedMuEnergyFraction());
        //
        if (mCHEn_lowPt_EndCap) mCHEn_lowPt_EndCap->Fill(jet->chargedHadronEnergy());
        if (mNHEn_lowPt_EndCap) mNHEn_lowPt_EndCap->Fill(jet->neutralHadronEnergy());
        if (mPhEn_lowPt_EndCap) mPhEn_lowPt_EndCap->Fill(jet->neutralEmEnergy());
        if (mElEn_lowPt_EndCap) mElEn_lowPt_EndCap->Fill(jet->chargedEmEnergy());
        if (mMuEn_lowPt_EndCap) mMuEn_lowPt_EndCap->Fill(jet->chargedMuEnergy());
      }
      if (jet->pt()>50. && jet->pt()<=140.) {
        if (mCHFrac_mediumPt_EndCap) mCHFrac_mediumPt_EndCap->Fill(jet->chargedHadronEnergyFraction());
        if (mNHFrac_mediumPt_EndCap) mNHFrac_mediumPt_EndCap->Fill(jet->neutralHadronEnergyFraction());
        if (mPhFrac_mediumPt_EndCap) mPhFrac_mediumPt_EndCap->Fill(jet->neutralEmEnergyFraction());
        if (mElFrac_mediumPt_EndCap) mElFrac_mediumPt_EndCap->Fill(jet->chargedEmEnergyFraction());
        if (mMuFrac_mediumPt_EndCap) mMuFrac_mediumPt_EndCap->Fill(jet->chargedMuEnergyFraction());
        //
        if (mCHEn_mediumPt_EndCap) mCHEn_mediumPt_EndCap->Fill(jet->chargedHadronEnergy());
        if (mNHEn_mediumPt_EndCap) mNHEn_mediumPt_EndCap->Fill(jet->neutralHadronEnergy());
        if (mPhEn_mediumPt_EndCap) mPhEn_mediumPt_EndCap->Fill(jet->neutralEmEnergy());
        if (mElEn_mediumPt_EndCap) mElEn_mediumPt_EndCap->Fill(jet->chargedEmEnergy());
        if (mMuEn_mediumPt_EndCap) mMuEn_mediumPt_EndCap->Fill(jet->chargedMuEnergy());
      }
      if (jet->pt()>140.) {
        if (mCHFrac_highPt_EndCap) mCHFrac_highPt_EndCap->Fill(jet->chargedHadronEnergyFraction());
        if (mNHFrac_highPt_EndCap) mNHFrac_highPt_EndCap->Fill(jet->neutralHadronEnergyFraction());
        if (mPhFrac_highPt_EndCap) mPhFrac_highPt_EndCap->Fill(jet->neutralEmEnergyFraction());
        if (mElFrac_highPt_EndCap) mElFrac_highPt_EndCap->Fill(jet->chargedEmEnergyFraction());
        if (mMuFrac_highPt_EndCap) mMuFrac_highPt_EndCap->Fill(jet->chargedMuEnergyFraction());
        //
        if (mCHEn_highPt_EndCap) mCHEn_highPt_EndCap->Fill(jet->chargedHadronEnergy());
        if (mNHEn_highPt_EndCap) mNHEn_highPt_EndCap->Fill(jet->neutralHadronEnergy());
        if (mPhEn_highPt_EndCap) mPhEn_highPt_EndCap->Fill(jet->neutralEmEnergy());
        if (mElEn_highPt_EndCap) mElEn_highPt_EndCap->Fill(jet->chargedEmEnergy());
        if (mMuEn_highPt_EndCap) mMuEn_highPt_EndCap->Fill(jet->chargedMuEnergy());
      }
      if(mChMultiplicity_lowPt_EndCap)  mChMultiplicity_lowPt_EndCap->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_lowPt_EndCap)  mNeuMultiplicity_lowPt_EndCap->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_lowPt_EndCap)  mMuMultiplicity_lowPt_EndCap->Fill(jet->muonMultiplicity());
      if(mChMultiplicity_mediumPt_EndCap)  mChMultiplicity_mediumPt_EndCap->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_mediumPt_EndCap)  mNeuMultiplicity_mediumPt_EndCap->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_mediumPt_EndCap)  mMuMultiplicity_mediumPt_EndCap->Fill(jet->muonMultiplicity());
      if(mChMultiplicity_highPt_EndCap)  mChMultiplicity_highPt_EndCap->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_highPt_EndCap)  mNeuMultiplicity_highPt_EndCap->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_highPt_EndCap)  mMuMultiplicity_highPt_EndCap->Fill(jet->muonMultiplicity());
      //
      if (mCHFracVSpT_EndCap) mCHFracVSpT_EndCap->Fill(jet->pt(),jet->chargedHadronEnergyFraction());
      if (mNHFracVSpT_EndCap) mNHFracVSpT_EndCap->Fill(jet->pt(),jet->neutralHadronEnergyFraction());
      if (mPhFracVSpT_EndCap) mPhFracVSpT_EndCap->Fill(jet->pt(),jet->neutralEmEnergyFraction());
      if (mElFracVSpT_EndCap) mElFracVSpT_EndCap->Fill(jet->pt(),jet->chargedEmEnergyFraction());
      if (mMuFracVSpT_EndCap) mMuFracVSpT_EndCap->Fill(jet->pt(),jet->chargedMuEnergyFraction());
    }
    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->chargedHadronEnergyFraction() + jet->neutralHadronEnergyFraction());    
    if (mEFrac_Forward)           mEFrac->Fill (jet->chargedEmEnergyFraction() + jet->neutralEmEnergyFraction());
      //fractions
      if (jet->pt()<=50.) {
        if(mHFEFrac_lowPt_Forward) mHFEFrac_lowPt_Forward->Fill(jet->HFEMEnergyFraction());
        if(mHFHFrac_lowPt_Forward) mHFHFrac_lowPt_Forward->Fill(jet->HFHadronEnergyFraction());
        //
        if(mHFEEn_lowPt_Forward) mHFEEn_lowPt_Forward->Fill(jet->HFEMEnergy());
        if(mHFHEn_lowPt_Forward) mHFHEn_lowPt_Forward->Fill(jet->HFHadronEnergy());
      }
      if (jet->pt()>50. && jet->pt()<=140.) {
        if(mHFEFrac_mediumPt_Forward) mHFEFrac_mediumPt_Forward->Fill(jet->HFEMEnergyFraction());
        if(mHFHFrac_mediumPt_Forward) mHFHFrac_mediumPt_Forward->Fill(jet->HFHadronEnergyFraction());
        //
        if(mHFEEn_mediumPt_Forward) mHFEEn_mediumPt_Forward->Fill(jet->HFEMEnergy());
        if(mHFHEn_mediumPt_Forward) mHFHEn_mediumPt_Forward->Fill(jet->HFHadronEnergy());
      }
      if (jet->pt()>140.) {
        if(mHFEFrac_highPt_Forward) mHFEFrac_highPt_Forward->Fill(jet->HFEMEnergyFraction());
        if(mHFHFrac_highPt_Forward) mHFHFrac_highPt_Forward->Fill(jet->HFHadronEnergyFraction());
        //
        if(mHFEEn_highPt_Forward) mHFEEn_highPt_Forward->Fill(jet->HFEMEnergy());
        if(mHFHEn_highPt_Forward) mHFHEn_highPt_Forward->Fill(jet->HFHadronEnergy());
      }
      if(mChMultiplicity_lowPt_Forward)  mChMultiplicity_lowPt_Forward->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_lowPt_Forward)  mNeuMultiplicity_lowPt_Forward->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_lowPt_Forward)  mMuMultiplicity_lowPt_Forward->Fill(jet->muonMultiplicity());
      if(mChMultiplicity_mediumPt_Forward)  mChMultiplicity_mediumPt_Forward->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_mediumPt_Forward)  mNeuMultiplicity_mediumPt_Forward->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_mediumPt_Forward)  mMuMultiplicity_mediumPt_Forward->Fill(jet->muonMultiplicity());
      if(mChMultiplicity_highPt_Forward)  mChMultiplicity_highPt_Forward->Fill(jet->chargedMultiplicity());
      if(mNeuMultiplicity_highPt_Forward)  mNeuMultiplicity_highPt_Forward->Fill(jet->neutralMultiplicity());
      if(mMuMultiplicity_highPt_Forward)  mMuMultiplicity_highPt_Forward->Fill(jet->muonMultiplicity());
      if(mHFHFracVSpT_Forward) mHFHFracVSpT_Forward->Fill(jet->pt(),jet->HFHadronEnergyFraction());
      if(mHFEFracVSpT_Forward) mHFEFracVSpT_Forward->Fill(jet->pt(),jet->HFEMEnergyFraction());
    }
    //if (mE) mE->Fill (jet->energy());
    //if (mP) mP->Fill (jet->p());
    //if (mMass) mMass->Fill (jet->mass());
            
    if (mChargedHadronEnergy)  mChargedHadronEnergy->Fill (jet->chargedHadronEnergy());
    if (mNeutralHadronEnergy)  mNeutralHadronEnergy->Fill (jet->neutralHadronEnergy());
    if (mChargedEmEnergy) mChargedEmEnergy->Fill(jet->chargedEmEnergy());
    if (mChargedMuEnergy) mChargedMuEnergy->Fill (jet->chargedMuEnergy ());
    if (mNeutralEmEnergy) mNeutralEmEnergy->Fill(jet->neutralEmEnergy());
    if (mChargedMultiplicity ) mChargedMultiplicity->Fill(jet->chargedMultiplicity());
    if (mNeutralMultiplicity ) mNeutralMultiplicity->Fill(jet->neutralMultiplicity());
    if (mMuonMultiplicity )mMuonMultiplicity->Fill (jet-> muonMultiplicity());
    //_______________________________________________________
    if (mNeutralFraction) mNeutralFraction->Fill (jet->neutralMultiplicity()/jet->nConstituents());


        // 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->chargedHadronEnergyFraction() + jet->neutralHadronEnergyFraction());
        if (mEFrac_profile)        mEFrac_profile       ->Fill(numPV, jet->chargedEmEnergyFraction()     + jet->neutralEmEnergyFraction());

        if (mChargedHadronEnergy_profile) mChargedHadronEnergy_profile->Fill(numPV, jet->chargedHadronEnergy());
        if (mNeutralHadronEnergy_profile) mNeutralHadronEnergy_profile->Fill(numPV, jet->neutralHadronEnergy());
        if (mChargedEmEnergy_profile)     mChargedEmEnergy_profile    ->Fill(numPV, jet->chargedEmEnergy());
        if (mChargedMuEnergy_profile)     mChargedMuEnergy_profile    ->Fill(numPV, jet->chargedMuEnergy ());
        if (mNeutralEmEnergy_profile)     mNeutralEmEnergy_profile    ->Fill(numPV, jet->neutralEmEnergy());
        if (mChargedMultiplicity_profile) mChargedMultiplicity_profile->Fill(numPV, jet->chargedMultiplicity());
        if (mNeutralMultiplicity_profile) mNeutralMultiplicity_profile->Fill(numPV, jet->neutralMultiplicity());
        if (mMuonMultiplicity_profile)    mMuonMultiplicity_profile   ->Fill(numPV, jet->muonMultiplicity());


    //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;
      }
    } // numofjets>1
      } // JetPt>_ptThreshold
    } // PF jet loop
    if (mNJets)   mNJets->Fill (numofjets);
    if (mDPhi)    mDPhi->Fill (dphi);


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


  } // non dijet selection
}

void PFJetAnalyzer::beginJob ( DQMStore *  dbe )

virtual

Inizialize parameters for histo binning.

Implements PFJetAnalyzerBase.

Definition at line 48 of file PFJetAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), jptDQMConfig_cff::eMax, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, LogTrace, metname, 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().

                                           {

  metname = "pFJetAnalyzer";

  LogTrace(metname)<<"[PFJetAnalyzer] Parameters initialization";
  //dbe->setCurrentFolder("JetMET/Jet/PFJets");//old version, now name set to source, which 
  //can be set for each instance of PFJetAnalyzer called inside JetMETAnalyzer. Useful, e.g., to 
  //name differently the dir for all jets and cleaned jets 
  dbe->setCurrentFolder("JetMET/Jet/"+_source);
  // dbe->setCurrentFolder("JetMET/Jet/PFJets");
            
  jetME = dbe->book1D("jetReco", "jetReco", 3, 1, 4);
  jetME->setBinLabel(2,"PFJets",1);

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

  _TightCHFMin = parameters.getParameter<double>("TightCHFMin");
  _TightNHFMax = parameters.getParameter<double>("TightNHFMax");
  _TightCEFMax = parameters.getParameter<double>("TightCEFMax");
  _TightNEFMax = parameters.getParameter<double>("TightNEFMax");
  _LooseCHFMin = parameters.getParameter<double>("LooseCHFMin");
  _LooseNHFMax = parameters.getParameter<double>("LooseNHFMax");
  _LooseCEFMax = parameters.getParameter<double>("LooseCEFMax");
  _LooseNEFMax = parameters.getParameter<double>("LooseNEFMax");

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

  _ThisCHFMin = parameters.getParameter<double>("ThisCHFMin");
  _ThisNHFMax = parameters.getParameter<double>("ThisNHFMax");
  _ThisCEFMax = parameters.getParameter<double>("ThisCEFMax");
  _ThisNEFMax = parameters.getParameter<double>("ThisNEFMax");

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


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


  //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", 50, 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)", 100, 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);
    
    mCHFracVSeta_lowPt= dbe->bookProfile("CHFracVSeta_lowPt","CHFracVSeta_lowPt",etaBin, etaMin, etaMax,0.,1.2);
    mNHFracVSeta_lowPt= dbe->bookProfile("NHFracVSeta_lowPt","NHFracVSeta_lowPt",etaBin, etaMin, etaMax,0.,1.2);
    mPhFracVSeta_lowPt= dbe->bookProfile("PhFracVSeta_lowPt","PhFracVSeta_lowPt",etaBin, etaMin, etaMax,0.,1.2);
    mElFracVSeta_lowPt= dbe->bookProfile("ElFracVSeta_lowPt","ElFracVSeta_lowPt",etaBin, etaMin, etaMax,0.,1.2);
    mMuFracVSeta_lowPt= dbe->bookProfile("MuFracVSeta_lowPt","MuFracVSeta_lowPt",etaBin, etaMin, etaMax,0.,1.2);
    mCHFracVSeta_mediumPt= dbe->bookProfile("CHFracVSeta_mediumPt","CHFracVSeta_mediumPt",etaBin, etaMin, etaMax,0.,1.2);
    mNHFracVSeta_mediumPt= dbe->bookProfile("NHFracVSeta_mediumPt","NHFracVSeta_mediumPt",etaBin, etaMin, etaMax,0.,1.2);
    mPhFracVSeta_mediumPt= dbe->bookProfile("PhFracVSeta_mediumPt","PhFracVSeta_mediumPt",etaBin, etaMin, etaMax,0.,1.2);
    mElFracVSeta_mediumPt= dbe->bookProfile("ElFracVSeta_mediumPt","ElFracVSeta_mediumPt",etaBin, etaMin, etaMax,0.,1.2);
    mMuFracVSeta_mediumPt= dbe->bookProfile("MuFracVSeta_mediumPt","MuFracVSeta_mediumPt",etaBin, etaMin, etaMax,0.,1.2);
    mCHFracVSeta_highPt= dbe->bookProfile("CHFracVSeta_highPt","CHFracVSeta_highPt",etaBin, etaMin, etaMax,0.,1.2);
    mNHFracVSeta_highPt= dbe->bookProfile("NHFracVSeta_highPt","NHFracVSeta_highPt",etaBin, etaMin, etaMax,0.,1.2);
    mPhFracVSeta_highPt= dbe->bookProfile("PhFracVSeta_highPt","PhFracVSeta_highPt",etaBin, etaMin, etaMax,0.,1.2);
    mElFracVSeta_highPt= dbe->bookProfile("ElFracVSeta_highPt","ElFracVSeta_highPt",etaBin, etaMin, etaMax,0.,1.2);
    mMuFracVSeta_highPt= dbe->bookProfile("MuFracVSeta_highPt","MuFracVSeta_highPt",etaBin, etaMin, etaMax,0.,1.2);
    
    //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", 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);
    // energy fractions
    mCHFrac_lowPt_Barrel     = dbe->book1D("CHFrac_lowPt_Barrel", "CHFrac_lowPt_Barrel", 120, -0.1, 1.1);
    mNHFrac_lowPt_Barrel     = dbe->book1D("NHFrac_lowPt_Barrel", "NHFrac_lowPt_Barrel", 120, -0.1, 1.1);
    mPhFrac_lowPt_Barrel     = dbe->book1D("PhFrac_lowPt_Barrel", "PhFrac_lowPt_Barrel", 120, -0.1, 1.1);
    mElFrac_lowPt_Barrel     = dbe->book1D("ElFrac_lowPt_Barrel", "ElFrac_lowPt_Barrel", 120, -0.1, 1.1);
    mMuFrac_lowPt_Barrel     = dbe->book1D("MuFrac_lowPt_Barrel", "MuFrac_lowPt_Barrel", 120, -0.1, 1.1);
    mCHFrac_mediumPt_Barrel  = dbe->book1D("CHFrac_mediumPt_Barrel", "CHFrac_mediumPt_Barrel", 120, -0.1, 1.1);
    mNHFrac_mediumPt_Barrel  = dbe->book1D("NHFrac_mediumPt_Barrel", "NHFrac_mediumPt_Barrel", 120, -0.1, 1.1);
    mPhFrac_mediumPt_Barrel  = dbe->book1D("PhFrac_mediumPt_Barrel", "PhFrac_mediumPt_Barrel", 120, -0.1, 1.1);
    mElFrac_mediumPt_Barrel  = dbe->book1D("ElFrac_mediumPt_Barrel", "ElFrac_mediumPt_Barrel", 120, -0.1, 1.1);
    mMuFrac_mediumPt_Barrel  = dbe->book1D("MuFrac_mediumPt_Barrel", "MuFrac_mediumPt_Barrel", 120, -0.1, 1.1);
    mCHFrac_highPt_Barrel    = dbe->book1D("CHFrac_highPt_Barrel", "CHFrac_highPt_Barrel", 120, -0.1, 1.1);
    mNHFrac_highPt_Barrel    = dbe->book1D("NHFrac_highPt_Barrel", "NHFrac_highPt_Barrel", 120, -0.1, 1.1);
    mPhFrac_highPt_Barrel    = dbe->book1D("PhFrac_highPt_Barrel", "PhFrac_highPt_Barrel", 120, -0.1, 1.1);
    mElFrac_highPt_Barrel    = dbe->book1D("ElFrac_highPt_Barrel", "ElFrac_highPt_Barrel", 120, -0.1, 1.1);
    mMuFrac_highPt_Barrel    = dbe->book1D("MuFrac_highPt_Barrel", "MuFrac_highPt_Barrel", 120, -0.1, 1.1);
    //energies
    mCHEn_lowPt_Barrel     = dbe->book1D("CHEn_lowPt_Barrel", "CHEn_lowPt_Barrel", ptBin, ptMin, ptMax);
    mNHEn_lowPt_Barrel     = dbe->book1D("NHEn_lowPt_Barrel", "NHEn_lowPt_Barrel", ptBin, ptMin, ptMax);
    mPhEn_lowPt_Barrel     = dbe->book1D("PhEn_lowPt_Barrel", "PhEn_lowPt_Barrel", ptBin, ptMin, ptMax);
    mElEn_lowPt_Barrel     = dbe->book1D("ElEn_lowPt_Barrel", "ElEn_lowPt_Barrel", ptBin, ptMin, ptMax);
    mMuEn_lowPt_Barrel     = dbe->book1D("MuEn_lowPt_Barrel", "MuEn_lowPt_Barrel", ptBin, ptMin, ptMax);
    mCHEn_mediumPt_Barrel  = dbe->book1D("CHEn_mediumPt_Barrel", "CHEn_mediumPt_Barrel", ptBin, ptMin, ptMax);
    mNHEn_mediumPt_Barrel  = dbe->book1D("NHEn_mediumPt_Barrel", "NHEn_mediumPt_Barrel", ptBin, ptMin, ptMax);
    mPhEn_mediumPt_Barrel  = dbe->book1D("PhEn_mediumPt_Barrel", "PhEn_mediumPt_Barrel", ptBin, ptMin, ptMax);
    mElEn_mediumPt_Barrel  = dbe->book1D("ElEn_mediumPt_Barrel", "ElEn_mediumPt_Barrel", ptBin, ptMin, ptMax);
    mMuEn_mediumPt_Barrel  = dbe->book1D("MuEn_mediumPt_Barrel", "MuEn_mediumPt_Barrel", ptBin, ptMin, ptMax);
    mCHEn_highPt_Barrel    = dbe->book1D("CHEn_highPt_Barrel", "CHEn_highPt_Barrel", ptBin, ptMin, ptMax);
    mNHEn_highPt_Barrel    = dbe->book1D("NHEn_highPt_Barrel", "NHEn_highPt_Barrel", ptBin, ptMin, ptMax);
    mPhEn_highPt_Barrel    = dbe->book1D("PhEn_highPt_Barrel", "PhEn_highPt_Barrel", ptBin, ptMin, ptMax);
    mElEn_highPt_Barrel    = dbe->book1D("ElEn_highPt_Barrel", "ElEn_highPt_Barrel", ptBin, ptMin, ptMax);
    mMuEn_highPt_Barrel    = dbe->book1D("MuEn_highPt_Barrel", "MuEn_highPt_Barrel", ptBin, ptMin, ptMax);
    //multiplicities
    mChMultiplicity_lowPt_Barrel    = dbe->book1D("ChMultiplicity_lowPt_Barrel", "ChMultiplicity_lowPt_Barrel", 30,0,30);
    mNeuMultiplicity_lowPt_Barrel   = dbe->book1D("NeuMultiplicity_lowPt_Barrel", "NeuMultiplicity_lowPt_Barrel", 30,0,30);
    mMuMultiplicity_lowPt_Barrel    = dbe->book1D("MuMultiplicity_lowPt_Barrel", "MuMultiplicity_lowPt_Barrel", 30,0,30);
    mChMultiplicity_mediumPt_Barrel    = dbe->book1D("ChMultiplicity_mediumPt_Barrel", "ChMultiplicity_mediumPt_Barrel", 30,0,30);
    mNeuMultiplicity_mediumPt_Barrel   = dbe->book1D("NeuMultiplicity_mediumPt_Barrel", "NeuMultiplicity_mediumPt_Barrel", 30,0,30);
    mMuMultiplicity_mediumPt_Barrel    = dbe->book1D("MuMultiplicity_mediumPt_Barrel", "MuMultiplicity_mediumPt_Barrel", 30,0,30);
    mChMultiplicity_highPt_Barrel    = dbe->book1D("ChMultiplicity_highPt_Barrel", "ChMultiplicity_highPt_Barrel", 30,0,30);
    mNeuMultiplicity_highPt_Barrel   = dbe->book1D("NeuMultiplicity_highPt_Barrel", "NeuMultiplicity_highPt_Barrel", 30,0,30);
    mMuMultiplicity_highPt_Barrel    = dbe->book1D("MuMultiplicity_highPt_Barrel", "MuMultiplicity_highPt_Barrel", 30,0,30);
    //
    mCHFracVSpT_Barrel= dbe->bookProfile("CHFracVSpT_Barrel","CHFracVSpT_Barrel",ptBin, ptMin, ptMax,0.,1.2);
    mNHFracVSpT_Barrel= dbe->bookProfile("NHFracVSpT_Barrel","NHFracVSpT_Barrel",ptBin, ptMin, ptMax,0.,1.2);
    mPhFracVSpT_Barrel= dbe->bookProfile("PhFracVSpT_Barrel","PhFracVSpT_Barrel",ptBin, ptMin, ptMax,0.,1.2);
    mElFracVSpT_Barrel= dbe->bookProfile("ElFracVSpT_Barrel","ElFracVSpT_Barrel",ptBin, ptMin, ptMax,0.,1.2);
    mMuFracVSpT_Barrel= dbe->bookProfile("MuFracVSpT_Barrel","MuFracVSpT_Barrel",ptBin, ptMin, ptMax,0.,1.2);
    mCHFracVSpT_EndCap= dbe->bookProfile("CHFracVSpT_EndCap","CHFracVSpT_EndCap",ptBin, ptMin, ptMax,0.,1.2);
    mNHFracVSpT_EndCap= dbe->bookProfile("NHFracVSpT_EndCap","NHFracVSpT_EndCap",ptBin, ptMin, ptMax,0.,1.2);
    mPhFracVSpT_EndCap= dbe->bookProfile("PhFracVSpT_EndCap","PhFracVSpT_EndCap",ptBin, ptMin, ptMax,0.,1.2);
    mElFracVSpT_EndCap= dbe->bookProfile("ElFracVSpT_EndCap","ElFracVSpT_EndCap",ptBin, ptMin, ptMax,0.,1.2);
    mMuFracVSpT_EndCap= dbe->bookProfile("MuFracVSpT_EndCap","MuFracVSpT_EndCap",ptBin, ptMin, ptMax,0.,1.2);
    mHFHFracVSpT_Forward= dbe->bookProfile("HFHFracVSpT_Forward","HFHFracVSpT_Forward",ptBin, ptMin, ptMax,0.,1.2);
    mHFEFracVSpT_Forward= dbe->bookProfile("HFEFracVSpT_Forward","HFEFracVSpT_Forward",ptBin, ptMin, ptMax,0.,1.2);

    mPhi_EndCap              = dbe->book1D("Phi_EndCap", "Phi_EndCap", phiBin, phiMin, phiMax);
    //mE_EndCap                = dbe->book1D("E_EndCap", "E_EndCap", eBin, eMin, eMax);
    mPt_EndCap               = dbe->book1D("Pt_EndCap", "Pt_EndCap", ptBin, ptMin, ptMax);
    //energy fractions
    mCHFrac_lowPt_EndCap     = dbe->book1D("CHFrac_lowPt_EndCap", "CHFrac_lowPt_EndCap", 120, -0.1, 1.1);
    mNHFrac_lowPt_EndCap     = dbe->book1D("NHFrac_lowPt_EndCap", "NHFrac_lowPt_EndCap", 120, -0.1, 1.1);
    mPhFrac_lowPt_EndCap     = dbe->book1D("PhFrac_lowPt_EndCap", "PhFrac_lowPt_EndCap", 120, -0.1, 1.1);
    mElFrac_lowPt_EndCap     = dbe->book1D("ElFrac_lowPt_EndCap", "ElFrac_lowPt_EndCap", 120, -0.1, 1.1);
    mMuFrac_lowPt_EndCap     = dbe->book1D("MuFrac_lowPt_EndCap", "MuFrac_lowPt_EndCap", 120, -0.1, 1.1);
    mCHFrac_mediumPt_EndCap  = dbe->book1D("CHFrac_mediumPt_EndCap", "CHFrac_mediumPt_EndCap", 120, -0.1, 1.1);
    mNHFrac_mediumPt_EndCap  = dbe->book1D("NHFrac_mediumPt_EndCap", "NHFrac_mediumPt_EndCap", 120, -0.1, 1.1);
    mPhFrac_mediumPt_EndCap  = dbe->book1D("PhFrac_mediumPt_EndCap", "PhFrac_mediumPt_EndCap", 120, -0.1, 1.1);
    mElFrac_mediumPt_EndCap  = dbe->book1D("ElFrac_mediumPt_EndCap", "ElFrac_mediumPt_EndCap", 120, -0.1, 1.1);
    mMuFrac_mediumPt_EndCap  = dbe->book1D("MuFrac_mediumPt_EndCap", "MuFrac_mediumPt_EndCap", 120, -0.1, 1.1);
    mCHFrac_highPt_EndCap    = dbe->book1D("CHFrac_highPt_EndCap", "CHFrac_highPt_EndCap", 120, -0.1, 1.1);
    mNHFrac_highPt_EndCap    = dbe->book1D("NHFrac_highPt_EndCap", "NHFrac_highPt_EndCap", 120, -0.1, 1.1);
    mPhFrac_highPt_EndCap    = dbe->book1D("PhFrac_highPt_EndCap", "PhFrac_highPt_EndCap", 120, -0.1, 1.1);
    mElFrac_highPt_EndCap    = dbe->book1D("ElFrac_highPt_EndCap", "ElFrac_highPt_EndCap", 120, -0.1, 1.1);
    mMuFrac_highPt_EndCap    = dbe->book1D("MuFrac_highPt_EndCap", "MuFrac_highPt_EndCap", 120, -0.1, 1.1);
    //energies
    mCHEn_lowPt_EndCap     = dbe->book1D("CHEn_lowPt_EndCap", "CHEn_lowPt_EndCap", ptBin, ptMin, ptMax);
    mNHEn_lowPt_EndCap     = dbe->book1D("NHEn_lowPt_EndCap", "NHEn_lowPt_EndCap", ptBin, ptMin, ptMax);
    mPhEn_lowPt_EndCap     = dbe->book1D("PhEn_lowPt_EndCap", "PhEn_lowPt_EndCap", ptBin, ptMin, ptMax);
    mElEn_lowPt_EndCap     = dbe->book1D("ElEn_lowPt_EndCap", "ElEn_lowPt_EndCap", ptBin, ptMin, ptMax);
    mMuEn_lowPt_EndCap     = dbe->book1D("MuEn_lowPt_EndCap", "MuEn_lowPt_EndCap", ptBin, ptMin, ptMax);
    mCHEn_mediumPt_EndCap  = dbe->book1D("CHEn_mediumPt_EndCap", "CHEn_mediumPt_EndCap", ptBin, ptMin, ptMax);
    mNHEn_mediumPt_EndCap  = dbe->book1D("NHEn_mediumPt_EndCap", "NHEn_mediumPt_EndCap", ptBin, ptMin, ptMax);
    mPhEn_mediumPt_EndCap  = dbe->book1D("PhEn_mediumPt_EndCap", "PhEn_mediumPt_EndCap", ptBin, ptMin, ptMax);
    mElEn_mediumPt_EndCap  = dbe->book1D("ElEn_mediumPt_EndCap", "ElEn_mediumPt_EndCap", ptBin, ptMin, ptMax);
    mMuEn_mediumPt_EndCap  = dbe->book1D("MuEn_mediumPt_EndCap", "MuEn_mediumPt_EndCap", ptBin, ptMin, ptMax);
    mCHEn_highPt_EndCap    = dbe->book1D("CHEn_highPt_EndCap", "CHEn_highPt_EndCap", ptBin, ptMin, ptMax);
    mNHEn_highPt_EndCap    = dbe->book1D("NHEn_highPt_EndCap", "NHEn_highPt_EndCap", ptBin, ptMin, ptMax);
    mPhEn_highPt_EndCap    = dbe->book1D("PhEn_highPt_EndCap", "PhEn_highPt_EndCap", ptBin, ptMin, ptMax);
    mElEn_highPt_EndCap    = dbe->book1D("ElEn_highPt_EndCap", "ElEn_highPt_EndCap", ptBin, ptMin, ptMax);
    mMuEn_highPt_EndCap    = dbe->book1D("MuEn_highPt_EndCap", "MuEn_highPt_EndCap", ptBin, ptMin, ptMax);
    //multiplicities
    mChMultiplicity_lowPt_EndCap    = dbe->book1D("ChMultiplicity_lowPt_EndCap", "ChMultiplicity_lowPt_EndCap", 30,0,30);
    mNeuMultiplicity_lowPt_EndCap   = dbe->book1D("NeuMultiplicity_lowPt_EndCap", "NeuMultiplicity_lowPt_EndCap", 30,0,30);
    mMuMultiplicity_lowPt_EndCap    = dbe->book1D("MuMultiplicity_lowPt_EndCap", "MuMultiplicity_lowPt_EndCap", 30,0,30);
    mChMultiplicity_mediumPt_EndCap    = dbe->book1D("ChMultiplicity_mediumPt_EndCap", "ChMultiplicity_mediumPt_EndCap", 30,0,30);
    mNeuMultiplicity_mediumPt_EndCap   = dbe->book1D("NeuMultiplicity_mediumPt_EndCap", "NeuMultiplicity_mediumPt_EndCap", 30,0,30);
    mMuMultiplicity_mediumPt_EndCap    = dbe->book1D("MuMultiplicity_mediumPt_EndCap", "MuMultiplicity_mediumPt_EndCap", 30,0,30);
    mChMultiplicity_highPt_EndCap    = dbe->book1D("ChMultiplicity_highPt_EndCap", "ChMultiplicity_highPt_EndCap", 30,0,30);
    mNeuMultiplicity_highPt_EndCap   = dbe->book1D("NeuMultiplicity_highPt_EndCap", "NeuMultiplicity_highPt_EndCap", 30,0,30);
    mMuMultiplicity_highPt_EndCap    = dbe->book1D("MuMultiplicity_highPt_EndCap", "MuMultiplicity_highPt_EndCap", 30,0,30);

    mPhi_Forward             = dbe->book1D("Phi_Forward", "Phi_Forward", phiBin, phiMin, phiMax);
    //mE_Forward               = dbe->book1D("E_Forward", "E_Forward", eBin, eMin, eMax);
    mPt_Forward              = dbe->book1D("Pt_Forward", "Pt_Forward", ptBin, ptMin, ptMax);
    //energy fraction
    mHFEFrac_lowPt_Forward    = dbe->book1D("HFEFrac_lowPt_Forward", "HFEFrac_lowPt_Forward", 120, -0.1, 1.1);
    mHFHFrac_lowPt_Forward    = dbe->book1D("HFHFrac_lowPt_Forward", "HFHFrac_lowPt_Forward", 120, -0.1, 1.1);
    mHFEFrac_mediumPt_Forward = dbe->book1D("HFEFrac_mediumPt_Forward", "HFEFrac_mediumPt_Forward", 120, -0.1, 1.1);
    mHFHFrac_mediumPt_Forward = dbe->book1D("HFHFrac_mediumPt_Forward", "HFHFrac_mediumPt_Forward", 120, -0.1, 1.1);
    mHFEFrac_highPt_Forward   = dbe->book1D("HFEFrac_highPt_Forward", "HFEFrac_highPt_Forward", 120, -0.1, 1.1);
    mHFHFrac_highPt_Forward   = dbe->book1D("HFHFrac_highPt_Forward", "HFHFrac_highPt_Forward", 120, -0.1, 1.1);
    //energies
    mHFEEn_lowPt_Forward    = dbe->book1D("HFEEn_lowPt_Forward", "HFEEn_lowPt_Forward", ptBin, ptMin, ptMax);
    mHFHEn_lowPt_Forward    = dbe->book1D("HFHEn_lowPt_Forward", "HFHEn_lowPt_Forward", ptBin, ptMin, ptMax);
    mHFEEn_mediumPt_Forward = dbe->book1D("HFEEn_mediumPt_Forward", "HFEEn_mediumPt_Forward", ptBin, ptMin, ptMax);
    mHFHEn_mediumPt_Forward = dbe->book1D("HFHEn_mediumPt_Forward", "HFHEn_mediumPt_Forward", ptBin, ptMin, ptMax);
    mHFEEn_highPt_Forward   = dbe->book1D("HFEEn_highPt_Forward", "HFEEn_highPt_Forward", ptBin, ptMin, ptMax);
    mHFHEn_highPt_Forward   = dbe->book1D("HFHEn_highPt_Forward", "HFHEn_highPt_Forward", ptBin, ptMin, ptMax);
    //multiplicities
    mChMultiplicity_lowPt_Forward     = dbe->book1D("ChMultiplicity_lowPt_Forward", "ChMultiplicity_lowPt_Forward", 30,0,30);
    mNeuMultiplicity_lowPt_Forward    = dbe->book1D("NeuMultiplicity_lowPt_Forward", "NeuMultiplicity_lowPt_Forward", 30,0,30);
    mMuMultiplicity_lowPt_Forward     = dbe->book1D("MuMultiplicity_lowPt_Forward", "MuMultiplicity_lowPt_Forward", 30,0,30);
    mChMultiplicity_mediumPt_Forward  = dbe->book1D("ChMultiplicity_mediumPt_Forward", "ChMultiplicity_mediumPt_Forward", 30,0,30);
    mNeuMultiplicity_mediumPt_Forward = dbe->book1D("NeuMultiplicity_mediumPt_Forward", "NeuMultiplicity_mediumPt_Forward", 30,0,30);
    mMuMultiplicity_mediumPt_Forward  = dbe->book1D("MuMultiplicity_mediumPt_Forward", "MuMultiplicity_mediumPt_Forward", 30,0,30);
    mChMultiplicity_highPt_Forward    = dbe->book1D("ChMultiplicity_highPt_Forward", "ChMultiplicity_highPt_Forward", 30,0,30);
    mNeuMultiplicity_highPt_Forward   = dbe->book1D("NeuMultiplicity_highPt_Forward", "NeuMultiplicity_highPt_Forward", 30,0,30);
    mMuMultiplicity_highPt_Forward    = dbe->book1D("MuMultiplicity_highPt_Forward", "MuMultiplicity_highPt_Forward", 30,0,30);



    // 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(fillpfJIDPassFrac==1) {
      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);
    }
  }


  mChargedHadronEnergy = dbe->book1D("mChargedHadronEnergy", "charged HAD energy",    100, 0, 100);
  mNeutralHadronEnergy = dbe->book1D("mNeutralHadronEnergy", "neutral HAD energy",    100, 0, 100);
  mChargedEmEnergy     = dbe->book1D("mChargedEmEnergy ",    "charged EM energy ",    100, 0, 100);
  mChargedMuEnergy     = dbe->book1D("mChargedMuEnergy",     "charged Mu energy",     100, 0, 100);
  mNeutralEmEnergy     = dbe->book1D("mNeutralEmEnergy",     "neutral EM energy",     100, 0, 100);
  mChargedMultiplicity = dbe->book1D("mChargedMultiplicity", "charged multiplicity ", 100, 0, 100);
  mNeutralMultiplicity = dbe->book1D("mNeutralMultiplicity", "neutral multiplicity",  100, 0, 100);
  mMuonMultiplicity    = dbe->book1D("mMuonMultiplicity",    "muon multiplicity",     100, 0, 100);

  
  // Book NPV profiles
  //----------------------------------------------------------------------------
  mChargedHadronEnergy_profile = dbe->bookProfile("mChargedHadronEnergy_profile", "charged HAD energy",   nbinsPV, PVlow, PVup, 100, 0, 100);
  mNeutralHadronEnergy_profile = dbe->bookProfile("mNeutralHadronEnergy_profile", "neutral HAD energy",   nbinsPV, PVlow, PVup, 100, 0, 100);
  mChargedEmEnergy_profile     = dbe->bookProfile("mChargedEmEnergy_profile",     "charged EM energy",    nbinsPV, PVlow, PVup, 100, 0, 100);
  mChargedMuEnergy_profile     = dbe->bookProfile("mChargedMuEnergy_profile",     "charged Mu energy",    nbinsPV, PVlow, PVup, 100, 0, 100);
  mNeutralEmEnergy_profile     = dbe->bookProfile("mNeutralEmEnergy_profile",     "neutral EM energy",    nbinsPV, PVlow, PVup, 100, 0, 100);
  mChargedMultiplicity_profile = dbe->bookProfile("mChargedMultiplicity_profile", "charged multiplicity", nbinsPV, PVlow, PVup, 100, 0, 100);
  mNeutralMultiplicity_profile = dbe->bookProfile("mNeutralMultiplicity_profile", "neutral multiplicity", nbinsPV, PVlow, PVup, 100, 0, 100);
  mMuonMultiplicity_profile    = dbe->bookProfile("mMuonMultiplicity_profile",    "muon multiplicity",    nbinsPV, PVlow, PVup, 100, 0, 100);

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


  // Set NPV profiles x-axis title
  //----------------------------------------------------------------------------
  mChargedHadronEnergy_profile->setAxisTitle("nvtx",1);
  mNeutralHadronEnergy_profile->setAxisTitle("nvtx",1);
  mChargedEmEnergy_profile    ->setAxisTitle("nvtx",1);
  mChargedMuEnergy_profile    ->setAxisTitle("nvtx",1);
  mNeutralEmEnergy_profile    ->setAxisTitle("nvtx",1);
  mChargedMultiplicity_profile->setAxisTitle("nvtx",1);
  mNeutralMultiplicity_profile->setAxisTitle("nvtx",1);
  mMuonMultiplicity_profile   ->setAxisTitle("nvtx",1);

  if (makedijetselection != 1) {
    mNJets_profile->setAxisTitle("nvtx",1);
  }


  //__________________________________________________
  mNeutralFraction     = dbe->book1D("NeutralFraction","Neutral Fraction",100,0,1);
  //}
  
  mDPhi                = dbe->book1D("DPhi", "dPhi btw the two leading jets", 100, 0., acos(-1.));
  
  if(makedijetselection==1) {
    mDijetAsymmetry                   = dbe->book1D("DijetAsymmetry", "DijetAsymmetry", 100, -1., 1.);
    mDijetBalance                     = dbe->book1D("DijetBalance",   "DijetBalance",   100, -2., 2.);
    if (fillpfJIDPassFrac==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 PFJetAnalyzer::endJob

(

void

)

Finish up a job.

Definition at line 419 of file PFJetAnalyzer.cc.

                           {
  
}

int PFJetAnalyzer::getLeadJetFlag ( )

inline

Definition at line 57 of file PFJetAnalyzer.h.

References _leadJetFlag.

                       {
    return  _leadJetFlag;
  }

void PFJetAnalyzer::setJetHiPass ( int  pass )

inline

Definition at line 64 of file PFJetAnalyzer.h.

References _JetHiPass.

                              {
    _JetHiPass = pass;
  }

void PFJetAnalyzer::setJetLoPass ( int  pass )

inline

Definition at line 60 of file PFJetAnalyzer.h.

References _JetLoPass.

                              {
    _JetLoPass = pass;
  }

void PFJetAnalyzer::setLeadJetFlag ( int  flag )

inline

Definition at line 54 of file PFJetAnalyzer.h.

References _leadJetFlag.

                                {
    _leadJetFlag = flag;
  }

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

inline

Definition at line 50 of file PFJetAnalyzer.h.

References _source, and LaserTracksInput_cfi::source.

                                   {
    _source = source;
  }

Member Data Documentation

double PFJetAnalyzer::_asymmetryThirdJetCut

private

Definition at line 76 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_balanceThirdJetCut

private

Definition at line 77 of file PFJetAnalyzer.h.

int PFJetAnalyzer::_JetHiPass

private

Definition at line 72 of file PFJetAnalyzer.h.

Referenced by setJetHiPass().

int PFJetAnalyzer::_JetLoPass

private

Definition at line 71 of file PFJetAnalyzer.h.

Referenced by setJetLoPass().

int PFJetAnalyzer::_leadJetFlag

private

Definition at line 73 of file PFJetAnalyzer.h.

Referenced by getLeadJetFlag(), and setLeadJetFlag().

double PFJetAnalyzer::_LooseCEFMax

private

Definition at line 110 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_LooseCHFMin

private

Definition at line 108 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_LooseNEFMax

private

Definition at line 111 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_LooseNHFMax

private

Definition at line 109 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_ptThreshold

private

Definition at line 74 of file PFJetAnalyzer.h.

std::string PFJetAnalyzer::_source

private

Definition at line 125 of file PFJetAnalyzer.h.

Referenced by setSource().

double PFJetAnalyzer::_ThisCEFMax

private

Definition at line 106 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_ThisCHFMin

private

Definition at line 104 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_ThisNEFMax

private

Definition at line 107 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_ThisNHFMax

private

Definition at line 105 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_TightCEFMax

private

Definition at line 114 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_TightCHFMin

private

Definition at line 112 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_TightNEFMax

private

Definition at line 115 of file PFJetAnalyzer.h.

double PFJetAnalyzer::_TightNHFMax

private

Definition at line 113 of file PFJetAnalyzer.h.

int PFJetAnalyzer::eBin

private

Definition at line 94 of file PFJetAnalyzer.h.

double PFJetAnalyzer::eMax

private

Definition at line 96 of file PFJetAnalyzer.h.

double PFJetAnalyzer::eMin

private

Definition at line 95 of file PFJetAnalyzer.h.

int PFJetAnalyzer::etaBin

private

Definition at line 82 of file PFJetAnalyzer.h.

double PFJetAnalyzer::etaMax

private

Definition at line 84 of file PFJetAnalyzer.h.

double PFJetAnalyzer::etaMin

private

Definition at line 83 of file PFJetAnalyzer.h.

int PFJetAnalyzer::fillpfJIDPassFrac

private

Definition at line 102 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::jetME

private

Definition at line 118 of file PFJetAnalyzer.h.

int PFJetAnalyzer::makedijetselection

private

Definition at line 79 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedEmEnergy

private

Definition at line 354 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedEmEnergy_profile

private

Definition at line 394 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedHadronEnergy

private

Definition at line 352 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedHadronEnergy_profile

private

Definition at line 392 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedMuEnergy

private

Definition at line 355 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedMuEnergy_profile

private

Definition at line 395 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedMultiplicity

private

Definition at line 357 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChargedMultiplicity_profile

private

Definition at line 397 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHEn_highPt_Barrel

private

Definition at line 180 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHEn_highPt_EndCap

private

Definition at line 250 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHEn_lowPt_Barrel

private

Definition at line 170 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHEn_lowPt_EndCap

private

Definition at line 240 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHEn_mediumPt_Barrel

private

Definition at line 175 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHEn_mediumPt_EndCap

private

Definition at line 245 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFrac_highPt_Barrel

private

Definition at line 165 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFrac_highPt_EndCap

private

Definition at line 234 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFrac_lowPt_Barrel

private

Definition at line 155 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFrac_lowPt_EndCap

private

Definition at line 224 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFrac_mediumPt_Barrel

private

Definition at line 160 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFrac_mediumPt_EndCap

private

Definition at line 229 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFracVSeta_highPt

private

Definition at line 218 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFracVSeta_lowPt

private

Definition at line 208 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFracVSeta_mediumPt

private

Definition at line 213 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFracVSpT_Barrel

private

Definition at line 195 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mCHFracVSpT_EndCap

private

Definition at line 200 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_highPt_Barrel

private

Definition at line 191 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_highPt_EndCap

private

Definition at line 262 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_highPt_Forward

private

Definition at line 291 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_lowPt_Barrel

private

Definition at line 185 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_lowPt_EndCap

private

Definition at line 256 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_lowPt_Forward

private

Definition at line 285 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_mediumPt_Barrel

private

Definition at line 188 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_mediumPt_EndCap

private

Definition at line 259 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mChMultiplicity_mediumPt_Forward

private

Definition at line 288 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents

private

Definition at line 147 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_Barrel

private

Definition at line 297 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_Barrel_Hi

private

Definition at line 313 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_EndCap

private

Definition at line 302 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_EndCap_Hi

private

Definition at line 317 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_Forward

private

Definition at line 307 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_Forward_Hi

private

Definition at line 321 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mConstituents_profile

private

Definition at line 388 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mDijetAsymmetry

private

Definition at line 379 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mDijetBalance

private

Definition at line 378 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mDPhi

private

Definition at line 334 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEEffChargedFraction

private

Definition at line 367 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEEffNeutralFraction

private

Definition at line 366 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEFrac

private

Definition at line 149 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEFrac_Barrel

private

Definition at line 299 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEFrac_EndCap

private

Definition at line 304 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEFrac_Forward

private

Definition at line 309 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEFrac_profile

private

Definition at line 390 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElEn_highPt_Barrel

private

Definition at line 183 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElEn_highPt_EndCap

private

Definition at line 253 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElEn_lowPt_Barrel

private

Definition at line 173 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElEn_lowPt_EndCap

private

Definition at line 243 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElEn_mediumPt_Barrel

private

Definition at line 178 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElEn_mediumPt_EndCap

private

Definition at line 248 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFrac_highPt_Barrel

private

Definition at line 168 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFrac_highPt_EndCap

private

Definition at line 237 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFrac_lowPt_Barrel

private

Definition at line 158 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFrac_lowPt_EndCap

private

Definition at line 227 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFrac_mediumPt_Barrel

private

Definition at line 163 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFrac_mediumPt_EndCap

private

Definition at line 232 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFracVSeta_highPt

private

Definition at line 221 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFracVSeta_lowPt

private

Definition at line 211 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFracVSeta_mediumPt

private

Definition at line 216 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFracVSpT_Barrel

private

Definition at line 198 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mElFracVSpT_EndCap

private

Definition at line 203 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEResChargedFraction

private

Definition at line 369 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEResNeutralFraction

private

Definition at line 368 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEta

private

Definition at line 145 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEta_Hi

private

Definition at line 347 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEta_profile

private

Definition at line 386 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mEtaFirst

private

Definition at line 337 of file PFJetAnalyzer.h.

std::string PFJetAnalyzer::metname

private

Definition at line 122 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEEn_highPt_Forward

private

Definition at line 283 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEEn_lowPt_Forward

private

Definition at line 279 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEEn_mediumPt_Forward

private

Definition at line 281 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEFrac_highPt_Forward

private

Definition at line 277 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEFrac_lowPt_Forward

private

Definition at line 273 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEFrac_mediumPt_Forward

private

Definition at line 275 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFEFracVSpT_Forward

private

Definition at line 206 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHEn_highPt_Forward

private

Definition at line 284 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHEn_lowPt_Forward

private

Definition at line 280 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHEn_mediumPt_Forward

private

Definition at line 282 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHFrac_highPt_Forward

private

Definition at line 278 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHFrac_lowPt_Forward

private

Definition at line 274 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHFrac_mediumPt_Forward

private

Definition at line 276 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFHFracVSpT_Forward

private

Definition at line 205 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac

private

Definition at line 148 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_Barrel

private

Definition at line 298 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_Barrel_Hi

private

Definition at line 314 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_EndCap

private

Definition at line 303 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_EndCap_Hi

private

Definition at line 318 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_Forward

private

Definition at line 308 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_Forward_Hi

private

Definition at line 322 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mHFrac_profile

private

Definition at line 389 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mLooseJIDPassFractionVSeta

private

Definition at line 372 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mLooseJIDPassFractionVSpt

private

Definition at line 373 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuEn_highPt_Barrel

private

Definition at line 184 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuEn_highPt_EndCap

private

Definition at line 254 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuEn_lowPt_Barrel

private

Definition at line 174 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuEn_lowPt_EndCap

private

Definition at line 244 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuEn_mediumPt_Barrel

private

Definition at line 179 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuEn_mediumPt_EndCap

private

Definition at line 249 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFrac_highPt_Barrel

private

Definition at line 169 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFrac_highPt_EndCap

private

Definition at line 238 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFrac_lowPt_Barrel

private

Definition at line 159 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFrac_lowPt_EndCap

private

Definition at line 228 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFrac_mediumPt_Barrel

private

Definition at line 164 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFrac_mediumPt_EndCap

private

Definition at line 233 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFracVSeta_highPt

private

Definition at line 222 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFracVSeta_lowPt

private

Definition at line 212 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFracVSeta_mediumPt

private

Definition at line 217 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFracVSpT_Barrel

private

Definition at line 199 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuFracVSpT_EndCap

private

Definition at line 204 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_highPt_Barrel

private

Definition at line 193 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_highPt_EndCap

private

Definition at line 264 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_highPt_Forward

private

Definition at line 293 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_lowPt_Barrel

private

Definition at line 187 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_lowPt_EndCap

private

Definition at line 258 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_lowPt_Forward

private

Definition at line 287 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_mediumPt_Barrel

private

Definition at line 190 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_mediumPt_EndCap

private

Definition at line 261 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuMultiplicity_mediumPt_Forward

private

Definition at line 290 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuonMultiplicity

private

Definition at line 359 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mMuonMultiplicity_profile

private

Definition at line 399 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_highPt_Barrel

private

Definition at line 192 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_highPt_EndCap

private

Definition at line 263 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_highPt_Forward

private

Definition at line 292 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_lowPt_Barrel

private

Definition at line 186 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_lowPt_EndCap

private

Definition at line 257 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_lowPt_Forward

private

Definition at line 286 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_mediumPt_Barrel

private

Definition at line 189 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_mediumPt_EndCap

private

Definition at line 260 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeuMultiplicity_mediumPt_Forward

private

Definition at line 289 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralEmEnergy

private

Definition at line 356 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralEmEnergy_profile

private

Definition at line 396 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralFraction

private

Definition at line 363 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralFraction2

private

Definition at line 364 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralHadronEnergy

private

Definition at line 353 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralHadronEnergy_profile

private

Definition at line 393 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralMultiplicity

private

Definition at line 358 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNeutralMultiplicity_profile

private

Definition at line 398 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHEn_highPt_Barrel

private

Definition at line 181 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHEn_highPt_EndCap

private

Definition at line 251 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHEn_lowPt_Barrel

private

Definition at line 171 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHEn_lowPt_EndCap

private

Definition at line 241 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHEn_mediumPt_Barrel

private

Definition at line 176 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHEn_mediumPt_EndCap

private

Definition at line 246 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFrac_highPt_Barrel

private

Definition at line 166 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFrac_highPt_EndCap

private

Definition at line 235 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFrac_lowPt_Barrel

private

Definition at line 156 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFrac_lowPt_EndCap

private

Definition at line 225 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFrac_mediumPt_Barrel

private

Definition at line 161 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFrac_mediumPt_EndCap

private

Definition at line 230 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFracVSeta_highPt

private

Definition at line 219 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFracVSeta_lowPt

private

Definition at line 209 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFracVSeta_mediumPt

private

Definition at line 214 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFracVSpT_Barrel

private

Definition at line 196 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNHFracVSpT_EndCap

private

Definition at line 201 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNJets

private

Definition at line 333 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mNJets_profile

private

Definition at line 384 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhEn_highPt_Barrel

private

Definition at line 182 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhEn_highPt_EndCap

private

Definition at line 252 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhEn_lowPt_Barrel

private

Definition at line 172 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhEn_lowPt_EndCap

private

Definition at line 242 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhEn_mediumPt_Barrel

private

Definition at line 177 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhEn_mediumPt_EndCap

private

Definition at line 247 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFrac_highPt_Barrel

private

Definition at line 167 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFrac_highPt_EndCap

private

Definition at line 236 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFrac_lowPt_Barrel

private

Definition at line 157 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFrac_lowPt_EndCap

private

Definition at line 226 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFrac_mediumPt_Barrel

private

Definition at line 162 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFrac_mediumPt_EndCap

private

Definition at line 231 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFracVSeta_highPt

private

Definition at line 220 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFracVSeta_lowPt

private

Definition at line 210 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFracVSeta_mediumPt

private

Definition at line 215 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFracVSpT_Barrel

private

Definition at line 197 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhFracVSpT_EndCap

private

Definition at line 202 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi

private

Definition at line 146 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_Barrel

private

Definition at line 153 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_Barrel_Hi

private

Definition at line 312 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_EndCap

private

Definition at line 268 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_EndCap_Hi

private

Definition at line 316 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_Forward

private

Definition at line 271 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_Forward_Hi

private

Definition at line 320 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_Hi

private

Definition at line 348 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_Lo

private

Definition at line 344 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhi_profile

private

Definition at line 387 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhiFirst

private

Definition at line 338 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPhiVSEta

private

Definition at line 150 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt

private

Definition at line 141 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_1

private

Definition at line 142 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_2

private

Definition at line 143 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_3

private

Definition at line 144 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_Barrel

private

Definition at line 152 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_Barrel_Hi

private

Definition at line 311 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_EndCap

private

Definition at line 267 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_EndCap_Hi

private

Definition at line 315 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_Forward

private

Definition at line 270 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_Forward_Hi

private

Definition at line 319 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_Hi

private

Definition at line 349 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_Lo

private

Definition at line 345 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPt_profile

private

Definition at line 385 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mPtFirst

private

Definition at line 340 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mTightJIDPassFractionVSeta

private

Definition at line 374 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::mTightJIDPassFractionVSpt

private

Definition at line 375 of file PFJetAnalyzer.h.

MonitorElement* PFJetAnalyzer::nEEff

private

Definition at line 370 of file PFJetAnalyzer.h.

edm::ParameterSet PFJetAnalyzer::parameters

private

Definition at line 120 of file PFJetAnalyzer.h.

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

int PFJetAnalyzer::pBin

private

Definition at line 98 of file PFJetAnalyzer.h.

int PFJetAnalyzer::phiBin

private

Definition at line 86 of file PFJetAnalyzer.h.

double PFJetAnalyzer::phiMax

private

Definition at line 88 of file PFJetAnalyzer.h.

double PFJetAnalyzer::phiMin

private

Definition at line 87 of file PFJetAnalyzer.h.

double PFJetAnalyzer::pMax

private

Definition at line 100 of file PFJetAnalyzer.h.

double PFJetAnalyzer::pMin

private

Definition at line 99 of file PFJetAnalyzer.h.

int PFJetAnalyzer::ptBin

private

Definition at line 90 of file PFJetAnalyzer.h.

double PFJetAnalyzer::ptMax

private

Definition at line 92 of file PFJetAnalyzer.h.

double PFJetAnalyzer::ptMin

private

Definition at line 91 of file PFJetAnalyzer.h.

edm::InputTag PFJetAnalyzer::thePFJetCollectionLabel

private

Definition at line 124 of file PFJetAnalyzer.h.