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