| 1 | #ifndef D0RunIconeJets_HepEntityPre96_class
|
|---|
| 2 | #define D0RunIconeJets_HepEntityPre96_class
|
|---|
| 3 |
|
|---|
| 4 | #include "inline_maths.h"
|
|---|
| 5 | #include "HepEntityI.h"
|
|---|
| 6 | #include <fastjet/internal/base.hh>
|
|---|
| 7 |
|
|---|
| 8 | FASTJET_BEGIN_NAMESPACE
|
|---|
| 9 |
|
|---|
| 10 | namespace d0runi{
|
|---|
| 11 |
|
|---|
| 12 | //Author: Lars Sonnenschein 25/Feb/2010
|
|---|
| 13 | //This is an example class fulfilling the minimal requirements needed by the
|
|---|
| 14 | //D0 RunI cone jet algorithm implementation prior to 1996, which is an inlined template class
|
|---|
| 15 | //See FERMILAB-Pub-97-242-E for details
|
|---|
| 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 | // 2011-11-14 Gregory Soyez <soyez@fastjet.fr>
|
|---|
| 30 | //
|
|---|
| 31 | // * removed some harmless warnings coming with the -Wshadow gcc option
|
|---|
| 32 | //
|
|---|
| 33 | // 2011-10-06 Gregory Soyez <soyez@fastjet.fr>
|
|---|
| 34 | //
|
|---|
| 35 | // * put the code in the fastjet::d0runi namespace
|
|---|
| 36 |
|
|---|
| 37 | class HepEntityIpre96 : public HepEntityI {
|
|---|
| 38 |
|
|---|
| 39 | public:
|
|---|
| 40 |
|
|---|
| 41 | HepEntityIpre96() {
|
|---|
| 42 | Et=0.;
|
|---|
| 43 | eta=0.;
|
|---|
| 44 | phi=0.;
|
|---|
| 45 | Ex=0.;
|
|---|
| 46 | Ey=0.;
|
|---|
| 47 | Ez=0.;
|
|---|
| 48 | index = -1;
|
|---|
| 49 | phi_pre96=0.;
|
|---|
| 50 | eta_pre96=0.;
|
|---|
| 51 |
|
|---|
| 52 | return;
|
|---|
| 53 | }
|
|---|
| 54 |
|
|---|
| 55 |
|
|---|
| 56 | HepEntityIpre96(double E_in, double px_in, double py_in, double pz_in,
|
|---|
| 57 | int index_in = -1) : index(index_in) {
|
|---|
| 58 | //Snowmass Et scheme
|
|---|
| 59 | double pt = sqrt(px_in*px_in+py_in*py_in);
|
|---|
| 60 | double p = sqrt(pt*pt+pz_in*pz_in);
|
|---|
| 61 | phi = inline_maths::phi(px_in,py_in);
|
|---|
| 62 | double theta = asin(pt/p);
|
|---|
| 63 | eta = inline_maths::eta(theta);
|
|---|
| 64 |
|
|---|
| 65 | Et = E_in*sin(theta);
|
|---|
| 66 |
|
|---|
| 67 | phi_pre96 = phi;
|
|---|
| 68 | eta_pre96 = eta;
|
|---|
| 69 |
|
|---|
| 70 | Ex = Et*cos(phi_pre96);
|
|---|
| 71 | Ey = Et*sin(phi_pre96);
|
|---|
| 72 | Ez = Et*sinh(eta_pre96);
|
|---|
| 73 |
|
|---|
| 74 | return;
|
|---|
| 75 | }
|
|---|
| 76 |
|
|---|
| 77 |
|
|---|
| 78 | inline double px() const {
|
|---|
| 79 | return Et*cos(phi_pre96);
|
|---|
| 80 | }
|
|---|
| 81 |
|
|---|
| 82 | inline double py() const {
|
|---|
| 83 | return Et*sin(phi_pre96);
|
|---|
| 84 | }
|
|---|
| 85 |
|
|---|
| 86 | inline double pz() const {
|
|---|
| 87 | return Et*sinh(eta_pre96);
|
|---|
| 88 | }
|
|---|
| 89 |
|
|---|
| 90 | inline double E() const {
|
|---|
| 91 | return Et*cosh(eta_pre96);
|
|---|
| 92 | }
|
|---|
| 93 |
|
|---|
| 94 |
|
|---|
| 95 | inline void Add(const HepEntityIpre96 el) {
|
|---|
| 96 | //assumes Et, eta and phi stored accurately
|
|---|
| 97 | double w2 = el.Et;
|
|---|
| 98 | Et += el.Et;
|
|---|
| 99 | w2 /= Et;
|
|---|
| 100 |
|
|---|
| 101 | eta += w2*(el.eta - eta);
|
|---|
| 102 | phi += w2*inline_maths::delta_phi(el.phi, phi);
|
|---|
| 103 |
|
|---|
| 104 |
|
|---|
| 105 | Ex += el.Ex;
|
|---|
| 106 | Ey += el.Ey;
|
|---|
| 107 | Ez += el.Ez;
|
|---|
| 108 | phi_pre96 = atan2(Ey, Ex);
|
|---|
| 109 | double theta_pre96 = atan2(sqrt(Ex*Ex+Ey*Ey),Ez);
|
|---|
| 110 | eta_pre96 = -log(tan(theta_pre96/2.));
|
|---|
| 111 |
|
|---|
| 112 | return;
|
|---|
| 113 | }
|
|---|
| 114 |
|
|---|
| 115 |
|
|---|
| 116 | inline void Fill(double E_in, double px_in, double py_in, double pz_in, int index_in) {
|
|---|
| 117 | double pt = sqrt(px_in*px_in+py_in*py_in);
|
|---|
| 118 | double p = sqrt(pt*pt+pz_in*pz_in);
|
|---|
| 119 | phi = inline_maths::phi(px_in,py_in);
|
|---|
| 120 | double theta = asin(pt/p);
|
|---|
| 121 | eta = inline_maths::eta(theta);
|
|---|
| 122 |
|
|---|
| 123 | Et = E_in*sin(theta);
|
|---|
| 124 |
|
|---|
| 125 |
|
|---|
| 126 | phi_pre96 = phi;
|
|---|
| 127 | eta_pre96 = eta;
|
|---|
| 128 |
|
|---|
| 129 | Ex = Et*cos(phi_pre96);
|
|---|
| 130 | Ey = Et*sin(phi_pre96);
|
|---|
| 131 | Ez = Et*sinh(eta_pre96);
|
|---|
| 132 |
|
|---|
| 133 | index = index_in;
|
|---|
| 134 |
|
|---|
| 135 | return;
|
|---|
| 136 | }
|
|---|
| 137 |
|
|---|
| 138 |
|
|---|
| 139 | double Ex;
|
|---|
| 140 | double Ey;
|
|---|
| 141 | double Ez;
|
|---|
| 142 | int index;
|
|---|
| 143 | double phi_pre96;
|
|---|
| 144 | double eta_pre96;
|
|---|
| 145 |
|
|---|
| 146 | private:
|
|---|
| 147 |
|
|---|
| 148 |
|
|---|
| 149 |
|
|---|
| 150 | };
|
|---|
| 151 | //end of class HepEntityIpre96;
|
|---|
| 152 |
|
|---|
| 153 | } // end of namespace d0runi
|
|---|
| 154 |
|
|---|
| 155 | FASTJET_END_NAMESPACE
|
|---|
| 156 |
|
|---|
| 157 | #endif
|
|---|
