20 |
|
{ |
21 |
|
public: |
22 |
|
CaloTower() {} |
23 |
< |
virtual ~CaloTower() {} |
23 |
> |
~CaloTower() {} |
24 |
|
|
25 |
< |
void SetEmEnergy(Double_t EmEnergy) { fEmEnergy = EmEnergy; } |
26 |
< |
void SetHadEnergy(Double_t HadEnergy) { fHadEnergy = HadEnergy; } |
27 |
< |
void SetOuterEnergy(Double_t OuterEnergy) { fOuterEnergy = OuterEnergy; } |
28 |
< |
void SetEmLvl1(Int_t EmLvl1) { fEmLvl1 = EmLvl1; } |
29 |
< |
void SetHadLvl1(Int_t HadLvl1) { fHadLvl1 = HadLvl1; } |
30 |
< |
void SetMomentum(Double_t X, Double_t Y, Double_t Z, Double_t E) |
31 |
< |
{ fMomentum.SetPxPyPzE(X,Y,Z,E); } |
25 |
> |
void SetEmEnergy(Double_t EmEnergy) { fEmEnergy = EmEnergy; } |
26 |
> |
void SetEmPosition(Double_t x, Double_t y, Double_t z) |
27 |
> |
{ fEmPosition.SetXYZ(x,y,z); } |
28 |
> |
void SetHadEnergy(Double_t HadEnergy) { fHadEnergy = HadEnergy; } |
29 |
> |
void SetHadPosition(Double_t x, Double_t y, Double_t z) |
30 |
> |
{ fHadPosition.SetXYZ(x,y,z); } |
31 |
> |
void SetOuterEnergy(Double_t OuterEnergy) { fOuterEnergy = OuterEnergy; } |
32 |
> |
void SetEmLvl1(Int_t EmLvl1) { fEmLvl1 = EmLvl1; } |
33 |
> |
void SetHadLvl1(Int_t HadLvl1) { fHadLvl1 = HadLvl1; } |
34 |
> |
void SetMomentum(Double_t x, Double_t y, Double_t z, Double_t e) |
35 |
> |
{ fMomentum.SetPxPyPzE(x,y,z,e); } |
36 |
|
|
37 |
< |
FourVector Mom() const { return fMomentum; } |
38 |
< |
Double_t Eta() const { return fMomentum.Eta(); } |
39 |
< |
Double_t Phi() const { return fMomentum.Phi(); } |
40 |
< |
Double_t Energy() const { return fMomentum.E(); } |
41 |
< |
Double_t Et() const { return fMomentum.E()*sin(fMomentum.Theta()); } |
42 |
< |
Double_t EmEnergy() const { return fEmEnergy; } |
43 |
< |
Double_t HadEnergy() const { return fHadEnergy; } |
44 |
< |
Double_t OuterEnergy() const { return fOuterEnergy; } |
45 |
< |
Double_t EmEt() const { return fEmEnergy*sin(fMomentum.Theta()); } |
46 |
< |
Double_t HadEt() const { return fHadEnergy*sin(fMomentum.Theta()); } |
47 |
< |
Double_t OuterEt() const { return fOuterEnergy*sin(fMomentum.Theta()); } |
48 |
< |
Double_t EmLvl1() const { return fEmLvl1; } |
49 |
< |
Double_t HadLvl1() const { return fHadLvl1; } |
37 |
> |
const FourVector &Mom() const { return fMomentum; } |
38 |
> |
Double_t Eta() const { return fMomentum.Eta(); } |
39 |
> |
Double_t Phi() const { return fMomentum.Phi(); } |
40 |
> |
Double_t Theta() const { return fMomentum.Theta(); } |
41 |
> |
Double_t E() const { return fMomentum.E(); } |
42 |
> |
Double_t Et() const { return E()*sin(Theta()); } |
43 |
> |
Double_t EmEnergy() const { return fEmEnergy; } |
44 |
> |
const ThreeVector &EmPosition() const { return fEmPosition; } |
45 |
> |
Double_t HadEnergy() const { return fHadEnergy; } |
46 |
> |
const ThreeVector &HadPosition() const { return fHadPosition; } |
47 |
> |
Double_t OuterEnergy() const { return fOuterEnergy; } |
48 |
> |
Double_t EmEt() const { return fEmEnergy*sin(Theta()); } |
49 |
> |
Double_t HadEt() const { return fHadEnergy*sin(Theta()); } |
50 |
> |
Double_t OuterEt() const { return fOuterEnergy*sin(Theta()); } |
51 |
> |
Double_t EmLvl1() const { return fEmLvl1; } |
52 |
> |
Double_t HadLvl1() const { return fHadLvl1; } |
53 |
|
|
54 |
|
protected: |
55 |
< |
|
56 |
< |
FourVector fMomentum; |
57 |
< |
Double_t fEmEnergy; |
58 |
< |
Double_t fHadEnergy; |
59 |
< |
Double_t fOuterEnergy; |
60 |
< |
Int_t fEmLvl1; |
61 |
< |
Int_t fHadLvl1; |
55 |
> |
FourVector fMomentum; |
56 |
> |
ThreeVector fEmPosition; //Position of Ecal shower center |
57 |
> |
ThreeVector fHadPosition; //Position of Hcal shower center |
58 |
> |
Double_t fEmEnergy; //tower energy in Ecal |
59 |
> |
Double_t fHadEnergy; //tower energy in Hcal |
60 |
> |
Double_t fOuterEnergy; |
61 |
> |
Int_t fEmLvl1; |
62 |
> |
Int_t fHadLvl1; |
63 |
|
|
64 |
|
ClassDef(CaloTower, 1) // Generic particle class |
65 |
|
}; |