1 | #ifndef D0RunIIconeJets_HepEntity_class |
---|
2 | #define D0RunIIconeJets_HepEntity_class |
---|
3 | |
---|
4 | #include "inline_maths.h" |
---|
5 | |
---|
6 | #include <fastjet/internal/base.hh> |
---|
7 | |
---|
8 | FASTJET_BEGIN_NAMESPACE |
---|
9 | |
---|
10 | namespace d0{ |
---|
11 | |
---|
12 | //Author: Lars Sonnenschein 28/Mar/2007 |
---|
13 | //This is an example class fulfilling the minimal requirements needed by the |
---|
14 | //D0 RunII cone jet algorithm implementation, which is an inlined template class |
---|
15 | // |
---|
16 | // |
---|
17 | // This file is distributed with FastJet under the terms of the GNU |
---|
18 | // General Public License (v2). Permission to do so has been granted |
---|
19 | // by Lars Sonnenschein and the D0 collaboration (see COPYING for |
---|
20 | // details) |
---|
21 | // |
---|
22 | // History of changes in FastJet compared tothe original version of |
---|
23 | // HepEntity.h |
---|
24 | // |
---|
25 | // 2011-12-13 Gregory Soyez <soyez@fastjet.fr> |
---|
26 | // |
---|
27 | // * added license information |
---|
28 | // |
---|
29 | // 2009-01-17 Gregory Soyez <soyez@fastjet.fr> |
---|
30 | // |
---|
31 | // * put the code in the fastjet::d0 namespace |
---|
32 | // |
---|
33 | // 2007-12-14 Gavin Salam <salam@lpthe.jussieu.fr> |
---|
34 | // |
---|
35 | // * added an index member |
---|
36 | |
---|
37 | class HepEntity { |
---|
38 | |
---|
39 | public: |
---|
40 | |
---|
41 | HepEntity() { |
---|
42 | E=0.; |
---|
43 | px=0.; |
---|
44 | py=0.; |
---|
45 | pz=0.; |
---|
46 | index = -1; |
---|
47 | return; |
---|
48 | } |
---|
49 | |
---|
50 | |
---|
51 | HepEntity(double E_in, double px_in, double py_in, double pz_in, int index_in = -1) : |
---|
52 | E(E_in), px(px_in), py(py_in), pz(pz_in), index(index_in) { |
---|
53 | return; |
---|
54 | } |
---|
55 | |
---|
56 | |
---|
57 | HepEntity(const HepEntity& in) : E(in.E), px(in.px), py(in.py), pz(in.pz), index(in.index) { |
---|
58 | return; |
---|
59 | } |
---|
60 | |
---|
61 | |
---|
62 | inline double y() const { |
---|
63 | return inline_maths::y(E,pz); |
---|
64 | } |
---|
65 | |
---|
66 | |
---|
67 | inline double phi() const { |
---|
68 | return inline_maths::phi(px,py); |
---|
69 | } |
---|
70 | |
---|
71 | |
---|
72 | inline double pT() const { |
---|
73 | return sqrt(inline_maths::sqr(px)+inline_maths::sqr(py)); |
---|
74 | } |
---|
75 | |
---|
76 | |
---|
77 | inline void p4vec(float* p) const { |
---|
78 | p[0] = px; |
---|
79 | p[1] = py; |
---|
80 | p[2] = pz; |
---|
81 | p[3] = E; |
---|
82 | return; |
---|
83 | } |
---|
84 | |
---|
85 | inline void Add(const HepEntity el) { |
---|
86 | E += el.E; |
---|
87 | px += el.px; |
---|
88 | py += el.py; |
---|
89 | pz += el.pz; |
---|
90 | return; |
---|
91 | } |
---|
92 | |
---|
93 | inline void Fill(double E_in, double px_in, double py_in, double pz_in, int index_in = -1) { |
---|
94 | E = E_in; |
---|
95 | px = px_in; |
---|
96 | py = py_in; |
---|
97 | pz = pz_in; |
---|
98 | index = index_in; |
---|
99 | return; |
---|
100 | } |
---|
101 | |
---|
102 | |
---|
103 | double E; |
---|
104 | double px; |
---|
105 | double py; |
---|
106 | double pz; |
---|
107 | int index; |
---|
108 | |
---|
109 | private: |
---|
110 | |
---|
111 | |
---|
112 | |
---|
113 | }; |
---|
114 | //end of class HepEntity; |
---|
115 | |
---|
116 | |
---|
117 | } // namespace d0 |
---|
118 | |
---|
119 | FASTJET_END_NAMESPACE |
---|
120 | |
---|
121 | #endif |
---|