[1] | 1 | #ifndef _PHYSICS_TOWER_HH_ |
---|
| 2 | #define _PHYSICS_TOWER_HH_ |
---|
| 3 | |
---|
| 4 | //---------------------------------------------------------------------- |
---|
| 5 | // This file distributed with FastJet has been obtained from |
---|
| 6 | // http://www.pa.msu.edu/~huston/Les_Houches_2005/JetClu+Midpoint-StandAlone.tgz |
---|
| 7 | // |
---|
| 8 | // Permission to distribute it with FastJet has been granted by Joey |
---|
| 9 | // Huston (see the COPYING file in the main FastJet directory for |
---|
| 10 | // details). |
---|
| 11 | // Changes from the original file are listed below. |
---|
| 12 | //---------------------------------------------------------------------- |
---|
| 13 | |
---|
| 14 | // History of changes compared to the original PhysicsTower.hh file |
---|
| 15 | // |
---|
| 16 | // 2009-01-17 Gregory Soyez <soyez@fastjet.fr> |
---|
| 17 | // |
---|
| 18 | // * put the code in the fastjet::cdf namespace |
---|
| 19 | // |
---|
| 20 | // 2008-08-15 Gavin Salam <salam@lpthe.jussieu.fr> |
---|
| 21 | // |
---|
| 22 | // * switched JetClu plugin over to proper indexed tracking of |
---|
| 23 | // jet contents rather than a (dodgy) map based on particle |
---|
| 24 | // energies |
---|
| 25 | // |
---|
| 26 | // 2006-09-24 Gavin Salam <salam@lpthe.jussieu.fr> |
---|
| 27 | // |
---|
| 28 | // * replaced the private m_index variable by a public fjindex |
---|
| 29 | // one for tracking within FastJet |
---|
| 30 | // |
---|
| 31 | // 2006-09-24 Gavin Salam <salam@lpthe.jussieu.fr> |
---|
| 32 | // |
---|
| 33 | // * added JetClu+MidPoint to FastJet |
---|
| 34 | |
---|
| 35 | #include "LorentzVector.hh" |
---|
| 36 | #include "CalTower.hh" |
---|
| 37 | |
---|
| 38 | #include <fastjet/internal/base.hh> |
---|
| 39 | |
---|
| 40 | FASTJET_BEGIN_NAMESPACE |
---|
| 41 | |
---|
| 42 | namespace cdf{ |
---|
| 43 | |
---|
| 44 | class PhysicsTower |
---|
| 45 | { |
---|
| 46 | public: |
---|
| 47 | |
---|
| 48 | LorentzVector fourVector; |
---|
| 49 | CalTower calTower; |
---|
| 50 | |
---|
| 51 | PhysicsTower(): fourVector(LorentzVector()), calTower(CalTower()), fjindex(-1) {} |
---|
| 52 | PhysicsTower(LorentzVector v, CalTower c): fourVector(v), calTower(c), fjindex(-1) {} |
---|
| 53 | PhysicsTower(const PhysicsTower& p): fourVector(p.fourVector), calTower(p.calTower), fjindex(p.fjindex) {} |
---|
| 54 | PhysicsTower(CalTower c): |
---|
| 55 | 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) {} |
---|
| 56 | PhysicsTower(LorentzVector v): fourVector(v), calTower(CalTower(v.Et(),v.eta(),v.phi())), fjindex(-1) {} |
---|
| 57 | double Et() const {return calTower.Et;} |
---|
| 58 | double eta() const {return calTower.eta;} |
---|
| 59 | double phi() const {return calTower.phi;} |
---|
| 60 | int iEta() const {return calTower.iEta;} |
---|
| 61 | int iPhi() const {return calTower.iPhi;} |
---|
| 62 | bool isEqual(PhysicsTower p) |
---|
| 63 | { |
---|
| 64 | return fourVector.isEqual(p.fourVector) && calTower.isEqual(p.calTower); |
---|
| 65 | } |
---|
| 66 | /// addition by GPS (2008-08-15) for tracking within fastjet |
---|
| 67 | int fjindex; |
---|
| 68 | }; |
---|
| 69 | |
---|
| 70 | } // namespace cdf |
---|
| 71 | |
---|
| 72 | FASTJET_END_NAMESPACE |
---|
| 73 | |
---|
| 74 | #endif |
---|