File indexing completed on 2024-04-06 12:11:26
0001 #ifndef JetComparison_H
0002 #define JetComparison_H
0003
0004 #include "FWCore/Framework/interface/Frameworkfwd.h"
0005 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
0006
0007 #include "FWCore/Framework/interface/Event.h"
0008 #include "FWCore/Framework/interface/EventSetup.h"
0009 #include "FWCore/Framework/interface/ESHandle.h"
0010 #include "FWCore/Framework/interface/MakerMacros.h"
0011 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0012 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0013
0014 #include "FWCore/ServiceRegistry/interface/Service.h"
0015 #include "CommonTools/UtilAlgos/interface/TFileService.h"
0016
0017 #include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
0018 #include "DataFormats/JetReco/interface/Jet.h"
0019 #include "DataFormats/JetReco/interface/CaloJetCollection.h"
0020 #include "DataFormats/JetReco/interface/CaloJet.h"
0021 #include "DataFormats/JetReco/interface/GenJet.h"
0022
0023 #include <vector>
0024 #include <utility>
0025 #include <ostream>
0026 #include <string>
0027 #include <algorithm>
0028 #include <cmath>
0029
0030 #include <TFile.h>
0031 #include <TH1.h>
0032 #include <TGraphErrors.h>
0033 #include <TROOT.h>
0034 #include <TFile.h>
0035 #include <TH1.h>
0036 #include <TH2.h>
0037
0038 class JetComparison : public edm::one::EDAnalyzer<edm::one::SharedResources> {
0039 public:
0040 JetComparison(edm::ParameterSet const& conf);
0041 ~JetComparison() override = default;
0042 void analyze(edm::Event const& e, edm::EventSetup const& c) override;
0043 virtual void beginRun(edm::Run const&, edm::EventSetup const&);
0044 void endJob() override;
0045 double deltaR2(double eta0, double phi0, double eta, double phi);
0046
0047 private:
0048 const edm::EDGetTokenT<reco::CaloJetCollection> jetToken_;
0049 const edm::EDGetTokenT<reco::GenJetCollection> genjetToken_;
0050 const float fMinEnergy;
0051
0052 TGraphErrors* gr;
0053
0054 int nEvent;
0055 TH1F* meEtJet;
0056 TH1F* meEtGen;
0057 TH1F* meEtJetMatched;
0058 TH1F* meEtGenMatched;
0059 TH1F* meEtaJet;
0060 TH1F* meEtaGen;
0061 TH2F* meRatio;
0062 TH2F* meEnergyHcalvsEcal;
0063 TH1F* meEnergyHO;
0064 TH1F* meEnergyHcal;
0065 TH1F* meEnergyEcal;
0066 TH1F* meNumFiredTowers;
0067 TH1F* meEnergyEcalTower;
0068 TH1F* meEnergyHcalTower;
0069 TH1F* meTotEnergy;
0070 TH1F* meNumberJet;
0071 TH1F* meDistR;
0072 TH2F* meDistR_vs_eta;
0073 TH2F* meHadronicFrac_vs_eta;
0074 TH2F* meNTowers90_vs_eta;
0075 TH2F* meNTowers60_vs_eta;
0076 TH2F* meNTowers_vs_eta;
0077 };
0078
0079 #endif