---

PFJetBenchmark.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_Benchmark_PFJetBenchmark_h
#define RecoParticleFlow_Benchmark_PFJetBenchmark_h

#include "RecoParticleFlow/Benchmark/interface/PFBenchmarkAlgo.h"

#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/BasicJetCollection.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

//#include "FWCore/ServiceRegistry/interface/Service.h" 
#include "PhysicsTools/UtilAlgos/interface/TFileService.h"
#include "TH1F.h"
#include "TH2F.h"
#include <string>
#include <TFile.h>


class PFJetBenchmark;

class TH1F;
class TH2F;

//class DQMStore; // CMSSW_2_X_X not needed here?

class PFJetBenchmark {
        
 public:
        
  PFJetBenchmark();
  virtual ~PFJetBenchmark();
        
  void setup(
             std::string Filename,
             bool debug, 
             bool plotAgainstReco=0, 
             bool onlyTwoJets=1,
             double deltaRMax=0.1,
             std::string benchmarkLabel_ = "ParticleFlow", 
             double recPt = -1, 
             double maxEta = -1,
             DQMStore * dbe_store = NULL
             );
  void process(const reco::PFJetCollection& , const reco::GenJetCollection& );
  void gettrue (const reco::GenJet* truth, double& true_ChargedHadEnergy, 
                double& true_NeutralHadEnergy, double& true_NeutralEmEnergy);
  void printPFJet (const reco::PFJet*);
  void printGenJet (const reco::GenJet*);
  double resPtMax() const {return resPtMax_;};
  double resChargedHadEnergyMax() const {return resChargedHadEnergyMax_;};
  double resNeutralHadEnergyMax() const {return resNeutralHadEnergyMax_;};
  double resNeutralEmEnergyMax() const {return resNeutralEmEnergyMax_;};
//  void save();        
  void write();
        
 private:
                
  TFile *file_;
        
  // histograms
  // Jets inclusive  distributions  (Pt > 20 GeV)
  TH1F *hNjets;
  TH1F *hjetsPt;
  TH1F *hjetsEta;
        
  // delta Pt or E quantities for Barrel
  TH1F *hBRPt;
  TH1F *hBRPt20_40;
  TH1F *hBRPt40_60;
  TH1F *hBRPt60_80;
  TH1F *hBRPt80_100;
  TH1F *hBRPt100_150;
  TH1F *hBRPt150_200;
  TH1F *hBRPt200_250;
  TH1F *hBRPt250_300;
  TH1F *hBRPt300_400;
  TH1F *hBRPt400_500;
  TH1F *hBRPt500_750;
  TH1F *hBRCHE;
  TH1F *hBRNHE;
  TH1F *hBRNEE;
  TH1F *hBRneut;
  TH2F *hBRPtvsPt ;
  TH2F *hBRCHEvsPt ;
  TH2F *hBRNHEvsPt;
  TH2F *hBRNEEvsPt;
  TH2F *hBRneutvsPt;
        
  // delta Pt or E quantities for Endcap
  TH1F *hERPt ;
  TH1F *hERPt20_40;
  TH1F *hERPt40_60;
  TH1F *hERPt60_80;
  TH1F *hERPt80_100;
  TH1F *hERPt100_150;
  TH1F *hERPt150_200;
  TH1F *hERPt200_250;
  TH1F *hERPt250_300;
  TH1F *hERPt300_400;
  TH1F *hERPt400_500;
  TH1F *hERPt500_750;
  TH1F *hERCHE;
  TH1F *hERNHE;
  TH1F *hERNEE;
  TH1F *hERneut;
  TH2F *hERPtvsPt ;
  TH2F *hERCHEvsPt;
  TH2F *hERNHEvsPt ;
  TH2F *hERNEEvsPt;
  TH2F *hERneutvsPt;
        
  std::string outputFile_;      
 protected:
                
  PFBenchmarkAlgo *algo_;
  bool debug_;
  bool plotAgainstReco_;
  bool onlyTwoJets_;
  double deltaRMax_;
  double resPtMax_;
  double resChargedHadEnergyMax_;
  double resNeutralHadEnergyMax_;
  double resNeutralEmEnergyMax_; 
  double recPt_cut;
  double maxEta_cut;
  DQMStore *dbe_;
};

#endif // RecoParticleFlow_Benchmark_PFJetBenchmark_h

---