| 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
|
|---|
