CommonTools/AnalysisJet.h at master · pandolf/CommonTools · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
#ifndef AnalysisJet_h
#define AnalysisJet_h

#include "TLorentzVector.h"
#include <string>

using namespace std;

class AnalysisJet : public TLorentzVector {

 public:

// ---- Constructor & Destructor ----------------------------------------------------------------------------
  AnalysisJet( float x=0., float y=0., float z=0., float t=0.) : TLorentzVector( x, y, z, t ) {
    rmsCand=0.;
    ptD=0.;
    nCharged=0;
    nNeutral=0;
  }

  AnalysisJet( const TLorentzVector &v) : TLorentzVector( v ) {
    rmsCand=0.;
    ptD=0.;
    nCharged=0;
    nNeutral=0;
  }

  virtual ~AnalysisJet(){}
// ---- Member functions ------------------------------------------------------------------------
   Bool_t passedJetID( const std::string& strength ) const;

// ---- Data Members ----------------------------------------------------------------------------
  float rmsCand;
  float ptD;
  int nCharged; //DON'T USE
  int nNeutral; //DON'T USE
  float QGLikelihood;
  float QGLikelihood2012;
  float QGLikelihoodNoPU;


  //extra QGL variables
  float axis1;
  float axis2;
  float pull;
  float tana;

  float ptD_QC;
  float rmsCand_QC;
  float axis1_QC;
  float axis2_QC;
  float pull_QC;
  float tana_QC;

  int nChg_ptCut;
  int nChg_QC;
  int nChg_ptCut_QC;
  int nNeutral_ptCut;
  int nPFCand_QC_ptCut;

  float Rchg;
  float Rneutral;
  float R;
  float Rchg_QC;
  //

  float eChargedHadrons;
  float ePhotons;
  float eNeutralHadrons;
  float eElectrons;
  float eMuons;
  float eHFHadrons;
  float eHFEM;

  int nChargedHadrons;
  int nPhotons;
  int nNeutralHadrons;
  int nElectrons;
  int nMuons;
  int nHFHadrons;
  int nHFEM;

  float pt_preKinFit;
  float eta_preKinFit;
  float phi_preKinFit;
  float e_preKinFit;

  float ptGen;
  float etaGen;
  float phiGen;
  float eGen;

  float ptPart;
  float etaPart;
  float phiPart;
  float ePart;
  int pdgIdPart;
  int pdgIdPartMom;
  int pdgIdPartMomMom;

  float ptPart_status2;
  float etaPart_status2;
  float phiPart_status2;
  float ePart_status2;
  int pdgIdPart_status2;

  bool btag_loose() const;
  bool btag_medium() const;

  //btags:
  float trackCountingHighEffBJetTag;
  float trackCountingHighPurBJetTag;
  float combinedSecondaryVertexBJetTag;
  float simpleSecondaryVertexHighEffBJetTag;
  float simpleSecondaryVertexHighPurBJetTag;
  float jetBProbabilityBJetTag;
  float jetProbabilityBJetTag;

  float beta;
  float betaStar;


private:
  int nConstituents() const {return nChargedHadrons+nPhotons+nNeutralHadrons+nElectrons+nMuons+nHFHadrons+nHFEM ;}
};


#endif