CMS 3D CMS Logo

Classes | Functions | Variables
DijetRatio.h File Reference
#include <memory>
#include <string>
#include <iostream>
#include <map>
#include <algorithm>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "CLHEP/Vector/LorentzVector.h"
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"

Go to the source code of this file.

Classes

class  DijetRatio< Jet >
 

Functions

void hInit (TH1F *hJet[], const char *name)
 
template<class R >
void histoFill (TH1F *jetHisto[], edm::Handle< R > jetsRec, double eta1, double eta2)
 

Variables

const int histoSize = 5
 

Function Documentation

◆ hInit()

void hInit ( TH1F *  hJet[],
const char *  name 
)
inline

Definition at line 57 of file DijetRatio.h.

References mergeVDriftHistosByStation::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by DijetRatio< Jet >::beginJob().

57  {
58  int const binSize = 35;
59  float massBin[binSize + 1] = {100, 113, 132, 153, 176, 201, 229, 259, 292, 327, 366, 400,
60  453, 501, 553, 609, 669, 733, 802, 875, 954, 1038, 1127, 1222,
61  1323, 1431, 1546, 1667, 1796, 1934, 2079, 2233, 2396, 2569, 2752, 3000};
62 
63  // (jetEta1 > 0 && jetEta1 < 0.7), (jetEta2 > 0 && jetEta2 < 0.7 )
64  std::string tit = std::string(name) + "_Eta_innerEtaCut_outerEtaCut";
65  hJet[0] = new TH1F(tit.c_str(), "DiJet Mass", binSize, massBin);
66 
67  // (jetEta1 > 0.7 && jetEta1 < 1.3), (jetEta2 > 0.7 && jetEta2 < 1.3 )
68  tit = std::string(name) + "_Eta_0_innerEtaCut";
69  hJet[1] = new TH1F(tit.c_str(), "DiJet Mass", binSize, massBin);
70 
71  tit = std::string(name) + "_LeadJetEta";
72  hJet[2] = new TH1F(tit.c_str(), "1^{st} Leading Jet #eta", 120, -6., 6.);
73  tit = std::string(name) + "_SecondJetEta";
74  hJet[3] = new TH1F(tit.c_str(), "2^{nd} Leading Jet #eta", 120, -6., 6.);
75  tit = std::string(name) + "_numEvents";
76  hJet[4] = new TH1F(tit.c_str(), "No. of events", 10, 0., 10.);
77 
78  return;
79 }

◆ histoFill()

template<class R >
void histoFill ( TH1F *  jetHisto[],
edm::Handle< R >  jetsRec,
double  eta1,
double  eta2 
)

Definition at line 82 of file DijetRatio.h.

References StorageManager_cfg::e1, HLT_2024v14_cff::eta1, HLT_2024v14_cff::eta2, EgHLTOffHistBins_cfi::mass, and findQualityFiles::v.

Referenced by DijetRatio< Jet >::analyze().

82  {
83  //For no. of events
84  jetHisto[4]->Fill(1.);
85 
86  if ((*jetsRec).size() >= 2) {
87  double px1 = (*jetsRec)[0].px();
88  double py1 = (*jetsRec)[0].py();
89  double pz1 = (*jetsRec)[0].pz();
90  double e1 = (*jetsRec)[0].energy();
91  double jetEta1 = (*jetsRec)[0].eta();
92  jetHisto[2]->Fill(jetEta1);
93 
94  double px2 = (*jetsRec)[1].px();
95  double py2 = (*jetsRec)[1].py();
96  double pz2 = (*jetsRec)[1].pz();
97  double e2 = (*jetsRec)[1].energy();
98  double jetEta2 = (*jetsRec)[1].eta();
99  jetHisto[3]->Fill(jetEta2);
100 
101  CLHEP::HepLorentzVector v1(px1, py1, pz1, e1);
102  CLHEP::HepLorentzVector v2(px2, py2, pz2, e2);
103  CLHEP::HepLorentzVector v(0., 0., 0., 0.);
104  v = v1 + v2;
105  float mass = v.m();
106 
107  if (fabs(jetEta1) > 0.0 && fabs(jetEta1) < eta1)
108  if (fabs(jetEta2) > 0.0 && fabs(jetEta2) < eta1)
109  jetHisto[0]->Fill(mass);
110 
111  if (fabs(jetEta1) > eta1 && fabs(jetEta1) < eta2)
112  if (fabs(jetEta2) > eta1 && fabs(jetEta2) < eta2)
113  jetHisto[1]->Fill(mass);
114  }
115 } //histoFill

Variable Documentation

◆ histoSize

const int histoSize = 5

Definition at line 54 of file DijetRatio.h.

Referenced by DijetRatio< Jet >::endJob().