1 |
//--------------------------------------------------------------------------------------------------
|
2 |
// $Id: Jet.h,v 1.1 2008/07/07 16:02:25 loizides Exp $
|
3 |
//
|
4 |
// Jet
|
5 |
//
|
6 |
// Details to be worked out...
|
7 |
//
|
8 |
// Authors: C.Loizides
|
9 |
//--------------------------------------------------------------------------------------------------
|
10 |
|
11 |
#ifndef DATATREE_JET_H
|
12 |
#define DATATREE_JET_H
|
13 |
|
14 |
#include "MitAna/DataTree/interface/Particle.h"
|
15 |
#include "MitAna/DataTree/interface/Track.h"
|
16 |
|
17 |
namespace mithep
|
18 |
{
|
19 |
class Jet : public Particle
|
20 |
{
|
21 |
public:
|
22 |
Jet() : fMaxEInEmTowers(0), fMaxEInHadTowers(0), fEnergyFractionH(0), fEnergyFractionEm(0),
|
23 |
fHadEnergyInHB(0), fHadEnergyInHO(0), fHadEnergyInHE(0), fHadEnergyInHF(0),
|
24 |
fEmEnergyInEB(0), fEmEnergyInEE(0), fEmEnergyInHF(0), fTowersArea(0), fN(0),
|
25 |
fN60(0), fN90(0) {}
|
26 |
Jet(Double_t px, Double_t py, Double_t pz, Double_t e) : fMom(px,py,pz,e),
|
27 |
fMaxEInEmTowers(0), fMaxEInHadTowers(0), fEnergyFractionH(0), fEnergyFractionEm(0),
|
28 |
fHadEnergyInHB(0), fHadEnergyInHO(0), fHadEnergyInHE(0), fHadEnergyInHF(0),
|
29 |
fEmEnergyInEB(0), fEmEnergyInEE(0), fEmEnergyInHF(0), fTowersArea(0), fN(0),
|
30 |
fN60(0), fN90(0) {}
|
31 |
~Jet() {}
|
32 |
|
33 |
Double_t MaxEInEmTowers() const { return fMaxEInEmTowers; }
|
34 |
Double_t MaxEInHadTowers() const { return fMaxEInHadTowers; }
|
35 |
Double_t HadEnergyInHO() const { return fHadEnergyInHO; }
|
36 |
Double_t HadEnergyInHB() const { return fHadEnergyInHB; }
|
37 |
Double_t HadEnergyInHF() const { return fHadEnergyInHF; }
|
38 |
Double_t HadEnergyInHE() const { return fHadEnergyInHE; }
|
39 |
Double_t EmEnergyInEB() const { return fEmEnergyInEB; }
|
40 |
Double_t EmEnergyInEE() const { return fEmEnergyInEE; }
|
41 |
Double_t EmEnergyInHF() const { return fEmEnergyInHF; }
|
42 |
Double_t EnergyFractionH() const { return fEnergyFractionH; }
|
43 |
Double_t EnergyFractionEm() const { return fEnergyFractionEm; }
|
44 |
Double_t TowersArea() const { return fTowersArea; }
|
45 |
UShort_t N() const { return fN; }
|
46 |
UShort_t N60() const { return fN60; }
|
47 |
UShort_t N90() const { return fN90; }
|
48 |
|
49 |
void SetMaxEInEmTowers(Double_t val) { fMaxEInEmTowers = val; }
|
50 |
void SetMaxEInHadTowers(Double_t val) { fMaxEInHadTowers = val; }
|
51 |
void SetHadEnergyInHO(Double_t val) { fHadEnergyInHO = val; }
|
52 |
void SetHadEnergyInHB(Double_t val) { fHadEnergyInHB = val; }
|
53 |
void SetHadEnergyInHF(Double_t val) { fHadEnergyInHF = val; }
|
54 |
void SetHadEnergyInHE(Double_t val) { fHadEnergyInHE = val; }
|
55 |
void SetEmEnergyInEB(Double_t val) { fEmEnergyInEB = val; }
|
56 |
void SetEmEnergyInEE(Double_t val) { fEmEnergyInEE = val; }
|
57 |
void SetEmEnergyInHF(Double_t val) { fEmEnergyInHF = val; }
|
58 |
void SetEnergyFractionH(Double_t val) { fEnergyFractionH = val; }
|
59 |
void SetEnergyFractionEm(Double_t val) { fEnergyFractionEm = val; }
|
60 |
void SetTowersArea(Double_t val) { fTowersArea = val; }
|
61 |
void SetN(UShort_t n) { fN = n; }
|
62 |
void SetN60(UShort_t n) { fN60 = n; }
|
63 |
void SetN90(UShort_t n) { fN90 = n; }
|
64 |
|
65 |
protected:
|
66 |
FourVector fMom; //four momentum of jet
|
67 |
Double32_t fMaxEInEmTowers; //maximum energy in EM towers
|
68 |
Double32_t fMaxEInHadTowers; //maximum energy in HCAL towers
|
69 |
Double32_t fEnergyFractionH; //hadronic energy fraction
|
70 |
Double32_t fEnergyFractionEm; //electromagnetic energy fraction
|
71 |
Double32_t fHadEnergyInHB; //hadronic energy in HB
|
72 |
Double32_t fHadEnergyInHO; //hadronic energy in HO
|
73 |
Double32_t fHadEnergyInHE; //hadronic energy in HE
|
74 |
Double32_t fHadEnergyInHF; //hadronic energy in HF
|
75 |
Double32_t fEmEnergyInEB; //electromagnetic energy in EB
|
76 |
Double32_t fEmEnergyInEE; //electromagnetic energy in EE
|
77 |
Double32_t fEmEnergyInHF; //electromagnetic energy extracted from HF
|
78 |
Double32_t fTowersArea; //area of contributing towers
|
79 |
UShort_t fN; //number of constituents
|
80 |
UShort_t fN60; //number constituents with 60% of total energy
|
81 |
UShort_t fN90; //number constituents with 90% of total energy
|
82 |
|
83 |
ClassDef(Jet, 1) // Jet class
|
84 |
};
|
85 |
}
|
86 |
#endif
|