| 1 | #ifndef _PHYSICS_TOWER_HH_
|
|---|
| 2 | #define _PHYSICS_TOWER_HH_
|
|---|
| 3 |
|
|---|
| 4 | #include "../interface/LorentzVector.hh"
|
|---|
| 5 | #include "../interface/CalTower.hh"
|
|---|
| 6 |
|
|---|
| 7 | class PhysicsTower
|
|---|
| 8 | {
|
|---|
| 9 | public:
|
|---|
| 10 |
|
|---|
| 11 | LorentzVector fourVector;
|
|---|
| 12 | CalTower calTower;
|
|---|
| 13 |
|
|---|
| 14 | PhysicsTower(): fourVector(LorentzVector()), calTower(CalTower()), fjindex(-1) {}
|
|---|
| 15 | PhysicsTower(LorentzVector v, CalTower c): fourVector(v), calTower(c), fjindex(-1) {}
|
|---|
| 16 | PhysicsTower(const PhysicsTower& p): fourVector(p.fourVector), calTower(p.calTower), fjindex(p.fjindex) {}
|
|---|
| 17 | PhysicsTower(CalTower c):
|
|---|
| 18 | fourVector(LorentzVector(c.Et*cos(c.phi),c.Et*sin(c.phi),c.Et*sinh(c.eta),c.Et*cosh(c.eta))), calTower(c), fjindex(-1) {}
|
|---|
| 19 | PhysicsTower(LorentzVector v): fourVector(v), calTower(CalTower(v.Et(),v.eta(),v.phi())), fjindex(-1) {}
|
|---|
| 20 | double Et() const {return calTower.Et;}
|
|---|
| 21 | double eta() const {return calTower.eta;}
|
|---|
| 22 | double phi() const {return calTower.phi;}
|
|---|
| 23 | int iEta() const {return calTower.iEta;}
|
|---|
| 24 | int iPhi() const {return calTower.iPhi;}
|
|---|
| 25 | bool isEqual(PhysicsTower p)
|
|---|
| 26 | {
|
|---|
| 27 | return fourVector.isEqual(p.fourVector) && calTower.isEqual(p.calTower);
|
|---|
| 28 | }
|
|---|
| 29 | /// addition by GPS (2008-08-15) for tracking within fastjet
|
|---|
| 30 | int fjindex;
|
|---|
| 31 | };
|
|---|
| 32 |
|
|---|
| 33 | #endif
|
|---|
