JetTester_HeavyIons.h

Go to the documentation of this file.

#ifndef ValidationRecoJetsJetTester_HeavyIons_h
#define ValidationRecoJetsJetTester_HeavyIons_h

// Producer for validation histograms for Calo, JPT and PF jet objects
// F. Ratnikov, Sept. 7, 2006
// Modified by Chiyoung Jeong, Feb. 2, 2010
// Modified by J. Piedra, Sept. 11, 2013
// Rewritten by Viola Sordini, Matthias Artur Weber, Robert Schoefbeck Nov./Dez. 2013
// Modified by Raghav Kunnawalkam Elayavalli, Aug 18th 2014 to run in 72X 
//                                          , Oct 22nd 2014 to run in 73X
//                                          , Dec 10th 2014 74X and adding the PF candidates information to easily detect  
//                                                          the voronoi subtraction algorithm failure modes.  

#include <cmath>
#include <string>

#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
#include "DataFormats/JetReco/interface/JPTJet.h"
#include "DataFormats/JetReco/interface/JPTJetCollection.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/JetReco/interface/BasicJet.h"
#include "DataFormats/JetReco/interface/BasicJetCollection.h"
#include "DataFormats/CaloTowers/interface/CaloTower.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/Candidate/interface/CandidateWithRef.h"

// include the pf candidates 
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
// include the voronoi subtraction
#include "DataFormats/HeavyIonEvent/interface/VoronoiBackground.h"
#include "RecoHI/HiJetAlgos/interface/UEParameters.h"
// include the centrality variables
#include "DataFormats/HeavyIonEvent/interface/Centrality.h"

#include "DataFormats/Math/interface/deltaR.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
#include "DataFormats/METReco/interface/GenMET.h"
#include "DataFormats/METReco/interface/GenMETCollection.h"
#include "DataFormats/METReco/interface/MET.h"
#include "DataFormats/METReco/interface/METCollection.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "JetMETCorrections/Objects/interface/JetCorrector.h"
#include "RecoJets/JetProducers/interface/JetMatchingTools.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"

const Int_t MAXPARTICLE = 10000;
const Double_t BarrelEta = 2.0;
const Double_t EndcapEta = 3.0;
const Double_t ForwardEta = 5.0;

class MonitorElement;

class JetTester_HeavyIons : public DQMEDAnalyzer {
 public:

  explicit JetTester_HeavyIons (const edm::ParameterSet&);
  ~JetTester_HeavyIons() override;

  void analyze(const edm::Event&, const edm::EventSetup&) override; 
  //virtual void beginJob(); 
  //virtual void endJob();
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  //reco::Vertex::Point getVtx(const edm::Event& ev);

  //double getEt(const DetID )

 private:
  
  void fillMatchHists(const double GenEta,  const double GenPhi,  const double GenPt,
              const double RecoEta, const double RecoPhi, const double RecoPt);
  
  edm::InputTag   mInputCollection;
  edm::InputTag   mInputGenCollection;
  edm::InputTag   mInputPFCandCollection;
  // edm::InputTag   mInputCandCollection;
  // edm::InputTag   rhoTag;
  edm::InputTag centralityTag_;
  edm::EDGetTokenT<reco::Centrality> centralityToken;
  edm::Handle<reco::Centrality> centrality_;

  edm::InputTag centralityBinTag_;
  edm::EDGetTokenT<int> centralityBinToken;
  edm::Handle<int>centralityBin_;
  
  std::string     mOutputFile;
  std::string     JetType;
  std::string     UEAlgo;
  edm::InputTag   Background;
  double          mRecoJetPtThreshold;
  double          mMatchGenPtThreshold;
  double          mGenEnergyFractionThreshold;
  double          mReverseEnergyFractionThreshold;
  double          mRThreshold;
  std::string     JetCorrectionService;

  //Tokens
  edm::EDGetTokenT<std::vector<reco::Vertex> > pvToken_;
  edm::EDGetTokenT<CaloTowerCollection > caloTowersToken_;
  edm::EDGetTokenT<reco::CaloJetCollection> caloJetsToken_;
  edm::EDGetTokenT<reco::PFJetCollection> pfJetsToken_;
  edm::EDGetTokenT<reco::BasicJetCollection> basicJetsToken_;
  edm::EDGetTokenT<reco::JPTJetCollection> jptJetsToken_;
  edm::EDGetTokenT<reco::GenJetCollection> genJetsToken_;
  edm::EDGetTokenT<GenEventInfoProduct> evtToken_;
  edm::EDGetTokenT<reco::PFCandidateCollection> pfCandToken_; 
  edm::EDGetTokenT<reco::CandidateView> pfCandViewToken_;
  edm::EDGetTokenT<reco::CandidateView> caloCandViewToken_;
  //edm::EDGetTokenT<reco::VoronoiMap> backgrounds_;
  edm::EDGetTokenT<edm::ValueMap<reco::VoronoiBackground>> backgrounds_;
  edm::EDGetTokenT<std::vector<float>> backgrounds_value_;
  edm::EDGetTokenT<std::vector<reco::Vertex> > hiVertexToken_;
  
  //Include Particle flow variables 
  MonitorElement *mNPFpart;
  MonitorElement *mPFPt;
  MonitorElement *mPFEta;
  MonitorElement *mPFPhi;
  MonitorElement *mPFVsPt;
  MonitorElement *mPFVsPtInitial;
  MonitorElement *mPFArea;
  MonitorElement *mNCalopart;
  MonitorElement *mCaloPt;
  MonitorElement *mCaloEta;
  MonitorElement *mCaloPhi;
  MonitorElement *mCaloVsPt;
  MonitorElement *mCaloVsPtInitial;
  MonitorElement *mCaloArea;
  MonitorElement *mSumpt;

  MonitorElement *mSumPFVsPt;
  MonitorElement *mSumPFVsPtInitial;
  MonitorElement *mSumPFPt;

  MonitorElement *mSumCaloVsPt;
  MonitorElement *mSumCaloVsPtInitial;
  MonitorElement *mSumCaloPt;

  MonitorElement *mSumSquaredPFVsPt;
  MonitorElement *mSumSquaredPFVsPtInitial;
  MonitorElement *mSumSquaredPFPt;

  MonitorElement *mSumSquaredCaloVsPt;
  MonitorElement *mSumSquaredCaloVsPtInitial;
  MonitorElement *mSumSquaredCaloPt;

  // Event variables (including centrality)
  MonitorElement* mNvtx;
  MonitorElement* mHF;

  // new additions Jan 12th 2015
  MonitorElement *mSumPFVsPt_HF;
  MonitorElement *mSumPFVsPtInitial_HF;
  MonitorElement *mSumPFPt_HF;
  MonitorElement *mDeltapT;
  MonitorElement *mDeltapT_eta;

  MonitorElement *mSumCaloVsPt_HF;
  MonitorElement *mSumCaloVsPtInitial_HF;
  MonitorElement *mSumCaloPt_HF;

  
  MonitorElement *mSumPFVsPtInitial_n5p191_n2p650;
  MonitorElement *mSumPFVsPtInitial_n2p650_n2p043;
  MonitorElement *mSumPFVsPtInitial_n2p043_n1p740;
  MonitorElement *mSumPFVsPtInitial_n1p740_n1p479;
  MonitorElement *mSumPFVsPtInitial_n1p479_n1p131;
  MonitorElement *mSumPFVsPtInitial_n1p131_n0p783;
  MonitorElement *mSumPFVsPtInitial_n0p783_n0p522;
  MonitorElement *mSumPFVsPtInitial_n0p522_0p522;
  MonitorElement *mSumPFVsPtInitial_0p522_0p783;
  MonitorElement *mSumPFVsPtInitial_0p783_1p131;
  MonitorElement *mSumPFVsPtInitial_1p131_1p479;
  MonitorElement *mSumPFVsPtInitial_1p479_1p740;
  MonitorElement *mSumPFVsPtInitial_1p740_2p043;
  MonitorElement *mSumPFVsPtInitial_2p043_2p650;
  MonitorElement *mSumPFVsPtInitial_2p650_5p191;

  MonitorElement *mSumPFVsPt_n5p191_n2p650;
  MonitorElement *mSumPFVsPt_n2p650_n2p043;
  MonitorElement *mSumPFVsPt_n2p043_n1p740;
  MonitorElement *mSumPFVsPt_n1p740_n1p479;
  MonitorElement *mSumPFVsPt_n1p479_n1p131;
  MonitorElement *mSumPFVsPt_n1p131_n0p783;
  MonitorElement *mSumPFVsPt_n0p783_n0p522;
  MonitorElement *mSumPFVsPt_n0p522_0p522;
  MonitorElement *mSumPFVsPt_0p522_0p783;
  MonitorElement *mSumPFVsPt_0p783_1p131;
  MonitorElement *mSumPFVsPt_1p131_1p479;
  MonitorElement *mSumPFVsPt_1p479_1p740;
  MonitorElement *mSumPFVsPt_1p740_2p043;
  MonitorElement *mSumPFVsPt_2p043_2p650;
  MonitorElement *mSumPFVsPt_2p650_5p191;

  MonitorElement *mSumPFPt_n5p191_n2p650;
  MonitorElement *mSumPFPt_n2p650_n2p043;
  MonitorElement *mSumPFPt_n2p043_n1p740;
  MonitorElement *mSumPFPt_n1p740_n1p479;
  MonitorElement *mSumPFPt_n1p479_n1p131;
  MonitorElement *mSumPFPt_n1p131_n0p783;
  MonitorElement *mSumPFPt_n0p783_n0p522;
  MonitorElement *mSumPFPt_n0p522_0p522;
  MonitorElement *mSumPFPt_0p522_0p783;
  MonitorElement *mSumPFPt_0p783_1p131;
  MonitorElement *mSumPFPt_1p131_1p479;
  MonitorElement *mSumPFPt_1p479_1p740;
  MonitorElement *mSumPFPt_1p740_2p043;
  MonitorElement *mSumPFPt_2p043_2p650;
  MonitorElement *mSumPFPt_2p650_5p191;
  
  
  MonitorElement *mSumCaloVsPtInitial_n5p191_n2p650;
  MonitorElement *mSumCaloVsPtInitial_n2p650_n2p043;
  MonitorElement *mSumCaloVsPtInitial_n2p043_n1p740;
  MonitorElement *mSumCaloVsPtInitial_n1p740_n1p479;
  MonitorElement *mSumCaloVsPtInitial_n1p479_n1p131;
  MonitorElement *mSumCaloVsPtInitial_n1p131_n0p783;
  MonitorElement *mSumCaloVsPtInitial_n0p783_n0p522;
  MonitorElement *mSumCaloVsPtInitial_n0p522_0p522;
  MonitorElement *mSumCaloVsPtInitial_0p522_0p783;
  MonitorElement *mSumCaloVsPtInitial_0p783_1p131;
  MonitorElement *mSumCaloVsPtInitial_1p131_1p479;
  MonitorElement *mSumCaloVsPtInitial_1p479_1p740;
  MonitorElement *mSumCaloVsPtInitial_1p740_2p043;
  MonitorElement *mSumCaloVsPtInitial_2p043_2p650;
  MonitorElement *mSumCaloVsPtInitial_2p650_5p191;

  MonitorElement *mSumCaloVsPt_n5p191_n2p650;
  MonitorElement *mSumCaloVsPt_n2p650_n2p043;
  MonitorElement *mSumCaloVsPt_n2p043_n1p740;
  MonitorElement *mSumCaloVsPt_n1p740_n1p479;
  MonitorElement *mSumCaloVsPt_n1p479_n1p131;
  MonitorElement *mSumCaloVsPt_n1p131_n0p783;
  MonitorElement *mSumCaloVsPt_n0p783_n0p522;
  MonitorElement *mSumCaloVsPt_n0p522_0p522;
  MonitorElement *mSumCaloVsPt_0p522_0p783;
  MonitorElement *mSumCaloVsPt_0p783_1p131;
  MonitorElement *mSumCaloVsPt_1p131_1p479;
  MonitorElement *mSumCaloVsPt_1p479_1p740;
  MonitorElement *mSumCaloVsPt_1p740_2p043;
  MonitorElement *mSumCaloVsPt_2p043_2p650;
  MonitorElement *mSumCaloVsPt_2p650_5p191;

  MonitorElement *mSumCaloPt_n5p191_n2p650;
  MonitorElement *mSumCaloPt_n2p650_n2p043;
  MonitorElement *mSumCaloPt_n2p043_n1p740;
  MonitorElement *mSumCaloPt_n1p740_n1p479;
  MonitorElement *mSumCaloPt_n1p479_n1p131;
  MonitorElement *mSumCaloPt_n1p131_n0p783;
  MonitorElement *mSumCaloPt_n0p783_n0p522;
  MonitorElement *mSumCaloPt_n0p522_0p522;
  MonitorElement *mSumCaloPt_0p522_0p783;
  MonitorElement *mSumCaloPt_0p783_1p131;
  MonitorElement *mSumCaloPt_1p131_1p479;
  MonitorElement *mSumCaloPt_1p479_1p740;
  MonitorElement *mSumCaloPt_1p740_2p043;
  MonitorElement *mSumCaloPt_2p043_2p650;
  MonitorElement *mSumCaloPt_2p650_5p191;
  
  // Jet parameters
  MonitorElement* mEta;
  MonitorElement* mPhi;
  MonitorElement* mPt;
  MonitorElement* mP;
  MonitorElement* mEnergy;
  MonitorElement* mMass;
  MonitorElement* mConstituents;
  MonitorElement* mJetArea;
  MonitorElement* mjetpileup;
  MonitorElement* mNJets;
  MonitorElement* mNJets_40;

  // histograms added on Oct 27th to study the gen-reco
  MonitorElement* mGenEta;
  MonitorElement* mGenPhi;
  MonitorElement* mGenPt;
  MonitorElement* mPtHat;

  MonitorElement* mPtRecoOverGen_B_20_30_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_20_30_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_20_30_Cent_0_10;
  MonitorElement* mPtRecoOverGen_B_30_50_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_30_50_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_30_50_Cent_0_10;
  MonitorElement* mPtRecoOverGen_B_50_80_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_50_80_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_50_80_Cent_0_10;
  MonitorElement* mPtRecoOverGen_B_80_120_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_80_120_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_80_120_Cent_0_10;
  MonitorElement* mPtRecoOverGen_B_120_180_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_120_180_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_120_180_Cent_0_10;
  MonitorElement* mPtRecoOverGen_B_180_300_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_180_300_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_180_300_Cent_0_10;
  MonitorElement* mPtRecoOverGen_B_300_Inf_Cent_0_10;
  MonitorElement* mPtRecoOverGen_E_300_Inf_Cent_0_10;
  MonitorElement* mPtRecoOverGen_F_300_Inf_Cent_0_10;
  
  MonitorElement* mPtRecoOverGen_B_20_30_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_20_30_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_20_30_Cent_10_30;
  MonitorElement* mPtRecoOverGen_B_30_50_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_30_50_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_30_50_Cent_10_30;
  MonitorElement* mPtRecoOverGen_B_50_80_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_50_80_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_50_80_Cent_10_30;
  MonitorElement* mPtRecoOverGen_B_80_120_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_80_120_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_80_120_Cent_10_30;
  MonitorElement* mPtRecoOverGen_B_120_180_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_120_180_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_120_180_Cent_10_30;
  MonitorElement* mPtRecoOverGen_B_180_300_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_180_300_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_180_300_Cent_10_30;
  MonitorElement* mPtRecoOverGen_B_300_Inf_Cent_10_30;
  MonitorElement* mPtRecoOverGen_E_300_Inf_Cent_10_30;
  MonitorElement* mPtRecoOverGen_F_300_Inf_Cent_10_30;

  MonitorElement* mPtRecoOverGen_B_20_30_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_20_30_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_20_30_Cent_30_50;
  MonitorElement* mPtRecoOverGen_B_30_50_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_30_50_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_30_50_Cent_30_50;
  MonitorElement* mPtRecoOverGen_B_50_80_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_50_80_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_50_80_Cent_30_50;
  MonitorElement* mPtRecoOverGen_B_80_120_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_80_120_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_80_120_Cent_30_50;
  MonitorElement* mPtRecoOverGen_B_120_180_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_120_180_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_120_180_Cent_30_50;
  MonitorElement* mPtRecoOverGen_B_180_300_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_180_300_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_180_300_Cent_30_50;
  MonitorElement* mPtRecoOverGen_B_300_Inf_Cent_30_50;
  MonitorElement* mPtRecoOverGen_E_300_Inf_Cent_30_50;
  MonitorElement* mPtRecoOverGen_F_300_Inf_Cent_30_50;

  MonitorElement* mPtRecoOverGen_B_20_30_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_20_30_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_20_30_Cent_50_80;
  MonitorElement* mPtRecoOverGen_B_30_50_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_30_50_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_30_50_Cent_50_80;
  MonitorElement* mPtRecoOverGen_B_50_80_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_50_80_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_50_80_Cent_50_80;
  MonitorElement* mPtRecoOverGen_B_80_120_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_80_120_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_80_120_Cent_50_80;
  MonitorElement* mPtRecoOverGen_B_120_180_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_120_180_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_120_180_Cent_50_80;
  MonitorElement* mPtRecoOverGen_B_180_300_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_180_300_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_180_300_Cent_50_80;
  MonitorElement* mPtRecoOverGen_B_300_Inf_Cent_50_80;
  MonitorElement* mPtRecoOverGen_E_300_Inf_Cent_50_80;
  MonitorElement* mPtRecoOverGen_F_300_Inf_Cent_50_80;

  // generation profiles
  MonitorElement* mPtRecoOverGen_GenPt_B_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenPt_E_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenPt_F_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenPt_B_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenPt_E_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenPt_F_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenPt_B_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenPt_E_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenPt_F_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenPt_B_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenPt_E_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenPt_F_Cent_50_80;

  MonitorElement* mPtRecoOverGen_GenEta_20_30_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenEta_30_50_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenEta_50_80_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenEta_80_120_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenEta_120_180_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenEta_180_300_Cent_0_10;
  MonitorElement* mPtRecoOverGen_GenEta_300_Inf_Cent_0_10;

  MonitorElement* mPtRecoOverGen_GenEta_20_30_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenEta_30_50_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenEta_50_80_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenEta_80_120_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenEta_120_180_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenEta_180_300_Cent_10_30;
  MonitorElement* mPtRecoOverGen_GenEta_300_Inf_Cent_10_30;

  MonitorElement* mPtRecoOverGen_GenEta_20_30_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenEta_30_50_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenEta_50_80_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenEta_80_120_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenEta_120_180_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenEta_180_300_Cent_30_50;
  MonitorElement* mPtRecoOverGen_GenEta_300_Inf_Cent_30_50;
  
  MonitorElement* mPtRecoOverGen_GenEta_20_30_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenEta_30_50_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenEta_50_80_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenEta_80_120_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenEta_120_180_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenEta_180_300_Cent_50_80;
  MonitorElement* mPtRecoOverGen_GenEta_300_Inf_Cent_50_80;
  
  MonitorElement* mPFCandpT_vs_eta_Unknown; // pf id 0
  MonitorElement* mPFCandpT_vs_eta_ChargedHadron; // pf id - 1 
  MonitorElement* mPFCandpT_vs_eta_electron; // pf id - 2
  MonitorElement* mPFCandpT_vs_eta_muon; // pf id - 3
  MonitorElement* mPFCandpT_vs_eta_photon; // pf id - 4
  MonitorElement* mPFCandpT_vs_eta_NeutralHadron; // pf id - 5
  MonitorElement* mPFCandpT_vs_eta_HadE_inHF; // pf id - 6
  MonitorElement* mPFCandpT_vs_eta_EME_inHF; // pf id - 7
  
  MonitorElement* mPFCandpT_Barrel_Unknown; // pf id 0
  MonitorElement* mPFCandpT_Barrel_ChargedHadron; // pf id - 1 
  MonitorElement* mPFCandpT_Barrel_electron; // pf id - 2
  MonitorElement* mPFCandpT_Barrel_muon; // pf id - 3
  MonitorElement* mPFCandpT_Barrel_photon; // pf id - 4
  MonitorElement* mPFCandpT_Barrel_NeutralHadron; // pf id - 5
  MonitorElement* mPFCandpT_Barrel_HadE_inHF; // pf id - 6
  MonitorElement* mPFCandpT_Barrel_EME_inHF; // pf id - 7

  MonitorElement* mPFCandpT_Endcap_Unknown; // pf id 0
  MonitorElement* mPFCandpT_Endcap_ChargedHadron; // pf id - 1 
  MonitorElement* mPFCandpT_Endcap_electron; // pf id - 2
  MonitorElement* mPFCandpT_Endcap_muon; // pf id - 3
  MonitorElement* mPFCandpT_Endcap_photon; // pf id - 4
  MonitorElement* mPFCandpT_Endcap_NeutralHadron; // pf id - 5
  MonitorElement* mPFCandpT_Endcap_HadE_inHF; // pf id - 6
  MonitorElement* mPFCandpT_Endcap_EME_inHF; // pf id - 7

  MonitorElement* mPFCandpT_Forward_Unknown; // pf id 0
  MonitorElement* mPFCandpT_Forward_ChargedHadron; // pf id - 1 
  MonitorElement* mPFCandpT_Forward_electron; // pf id - 2
  MonitorElement* mPFCandpT_Forward_muon; // pf id - 3
  MonitorElement* mPFCandpT_Forward_photon; // pf id - 4
  MonitorElement* mPFCandpT_Forward_NeutralHadron; // pf id - 5
  MonitorElement* mPFCandpT_Forward_HadE_inHF; // pf id - 6
  MonitorElement* mPFCandpT_Forward_EME_inHF; // pf id - 7
  
 
  // Parameters

  bool            isCaloJet;
  bool            isJPTJet;
  bool            isPFJet;

  static const Int_t fourierOrder_ = 5;
  static const Int_t etaBins_ = 15;

  static const size_t nedge_pseudorapidity = etaBins_ + 1;


};

#endif