JetAnalyzer_HeavyIons.h

Go to the documentation of this file.

#ifndef JetAnalyzer_HeavyIons_H
#define JetAnalyzer_HeavyIons_H


//
// Jet Tester class for heavy ion jets. for DQM jet analysis monitoring 
// For CMSSW_7_4_X, especially reading background subtracted jets 
// author: Raghav Kunnawalkam Elayavalli,
//         Jan 12th 2015 
//         Rutgers University, email: raghav.k.e at CERN dot CH 
//
// this class will be very similar to the class available in the validation suite under RecoJets/JetTester_HeavyIons 
//


#include <memory>
#include <fstream>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
#include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
#include "FWCore/Utilities/interface/EDGetToken.h"

#include "DataFormats/JetReco/interface/Jet.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/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 basic jet for the PuPF jets. 
#include "DataFormats/JetReco/interface/BasicJet.h"
#include "DataFormats/JetReco/interface/BasicJetCollection.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 "JetMETCorrections/Objects/interface/JetCorrector.h"
#include "RecoJets/JetProducers/interface/JetIDHelper.h"
#include "DQMOffline/JetMET/interface/JetMETDQMDCSFilter.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "DataFormats/Scalers/interface/DcsStatus.h" 
#include "PhysicsTools/SelectorUtils/interface/JetIDSelectionFunctor.h"
#include "PhysicsTools/SelectorUtils/interface/PFJetIDSelectionFunctor.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include <map>
#include <string>


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 JetAnalyzer_HeavyIons : public DQMEDAnalyzer {

 public:

  explicit JetAnalyzer_HeavyIons (const edm::ParameterSet&);
  virtual ~JetAnalyzer_HeavyIons();
  
  void analyze(const edm::Event&, const edm::EventSetup&) override; 
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  
 private:
  static constexpr int fourierOrder_ = 5;
  static constexpr int etaBins_ = 7;
  static constexpr int ptBins_ = 7;

  //default values - these are changed by the etaMap values for the CS plots
  const double edge_pseudorapidity[etaBins_ +1] = {-5, -3, -2.1, -1.3, 1.3, 2.1, 3, 5};
  const int ptBin[ptBins_+1] = {0, 20, 40, 60, 100, 150, 300, 99999};

  static constexpr int nedge_pseudorapidity = etaBins_ + 1; 
 
  edm::InputTag   mInputCollection;
  edm::InputTag   mInputPFCandCollection;
  edm::InputTag   mInputCsCandCollection; 
 
  std::string     mOutputFile;
  std::string     JetType;
  std::string     UEAlgo;
  edm::InputTag   Background;
  double          mRecoJetPtThreshold;
  double          mReverseEnergyFractionThreshold;
  double          mRThreshold;
  std::string     JetCorrectionService;

  edm::InputTag centralityTag_;
  edm::EDGetTokenT<reco::Centrality> centralityToken;
  edm::Handle<reco::Centrality> centrality_;

  edm::InputTag centralityBinTag_;
  edm::EDGetTokenT<int> centralityBinToken;
  edm::Handle<int>centralityBin_;

  //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::PFCandidateCollection> pfCandToken_; 
  edm::EDGetTokenT<reco::PFCandidateCollection> csCandToken_;
  edm::EDGetTokenT<reco::CandidateView> pfCandViewToken_;
  edm::EDGetTokenT<reco::CandidateView> caloCandViewToken_;

  edm::EDGetTokenT<std::vector<reco::Vertex> > hiVertexToken_;

  edm::EDGetTokenT<std::vector<double>>                  etaToken_;
  edm::EDGetTokenT<std::vector<double>>                  rhoToken_;
  edm::EDGetTokenT<std::vector<double>>                  rhomToken_;

  MonitorElement *mNPFpart;
  MonitorElement *mPFPt;
  MonitorElement *mPFEta;
  MonitorElement *mPFPhi;

  MonitorElement *mPFArea;
  MonitorElement *mPFDeltaR; //MZ
  MonitorElement *mPFDeltaR_Scaled_R; //MZ
  
  MonitorElement *mNCalopart;
  MonitorElement *mCaloPt;
  MonitorElement *mCaloEta;
  MonitorElement *mCaloPhi;
  MonitorElement *mCaloArea;
  MonitorElement *mSumpt;
  MonitorElement *mSumPFPt;
  MonitorElement *mSumPFPt_eta;

  MonitorElement *mSumCaloPt;
  MonitorElement *mSumCaloPt_eta;

  MonitorElement *mSumSquaredPFPt;
  MonitorElement *mSumSquaredPFPt_eta;

  MonitorElement *mSumSquaredCaloPt;
  MonitorElement *mSumSquaredCaloPt_eta;  

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

  // new additions Jan 12th 2015
  MonitorElement *mSumPFPt_HF;
  MonitorElement *mSumCaloPt_HF;
  
  MonitorElement *mSumPFPtEtaDep[etaBins_];
  MonitorElement *mSumCaloPtEtaDep[etaBins_];
  
  // 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;

  MonitorElement* mRhoDist_vsEta;
  MonitorElement* mRhoMDist_vsEta;
  MonitorElement* mRhoDist_vsPt;
  MonitorElement* mRhoMDist_vsPt;
  MonitorElement* mRhoDist_vsCent;
  MonitorElement* mRhoMDist_vsCent;
  MonitorElement* rhoEtaRange;
  MonitorElement* mCSCandpT_vsPt[etaBins_];
  MonitorElement* mCSCand_corrPFcand[etaBins_];
  MonitorElement* mSubtractedEFrac[ptBins_][etaBins_];
  MonitorElement* mSubtractedE[ptBins_][etaBins_];

  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;

};


#endif