Back to home page

Project CMSSW displayed by LXR

 
 

    


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