Loading...
Searching...
No Matches
HepEntityIpre96.h
#ifndef D0RunIconeJets_HepEntityPre96_class
#define D0RunIconeJets_HepEntityPre96_class
#include "inline_maths.h"
#include "HepEntityI.h"
#include <fastjet/internal/base.hh>
FASTJET_BEGIN_NAMESPACE
namespace d0runi{
//Author: Lars Sonnenschein 25/Feb/2010
//This is an example class fulfilling the minimal requirements needed by the
//D0 RunI cone jet algorithm implementation prior to 1996, which is an inlined template class
//See FERMILAB-Pub-97-242-E for details
// This file is distributed with FastJet under the terms of the GNU
// General Public License (v2). Permission to do so has been granted
// by Lars Sonnenschein and the D0 collaboration (see COPYING for
// details)
//
// History of changes in FastJet compared tothe original version of
// HepEntity.h
//
// 2011-12-13 Gregory Soyez <soyez@fastjet.fr>
//
// * added license information
//
// 2011-11-14 Gregory Soyez <soyez@fastjet.fr>
//
// * removed some harmless warnings coming with the -Wshadow gcc option
//
// 2011-10-06 Gregory Soyez <soyez@fastjet.fr>
//
// * put the code in the fastjet::d0runi namespace
class HepEntityIpre96 : public HepEntityI {
public:
HepEntityIpre96() {
Et=0.;
eta=0.;
phi=0.;
Ex=0.;
Ey=0.;
Ez=0.;
index = -1;
phi_pre96=0.;
eta_pre96=0.;
return;
}
HepEntityIpre96(double E_in, double px_in, double py_in, double pz_in,
int index_in = -1) : index(index_in) {
//Snowmass Et scheme
double pt = sqrt(px_in*px_in+py_in*py_in);
double p = sqrt(pt*pt+pz_in*pz_in);
phi = inline_maths::phi(px_in,py_in);
double theta = asin(pt/p);
eta = inline_maths::eta(theta);
Et = E_in*sin(theta);
phi_pre96 = phi;
eta_pre96 = eta;
Ex = Et*cos(phi_pre96);
Ey = Et*sin(phi_pre96);
Ez = Et*sinh(eta_pre96);
return;
}
inline double px() const {
return Et*cos(phi_pre96);
}
inline double py() const {
return Et*sin(phi_pre96);
}
inline double pz() const {
return Et*sinh(eta_pre96);
}
inline double E() const {
return Et*cosh(eta_pre96);
}
inline void Add(const HepEntityIpre96 el) {
//assumes Et, eta and phi stored accurately
double w2 = el.Et;
Et += el.Et;
w2 /= Et;
eta += w2*(el.eta - eta);
phi += w2*inline_maths::delta_phi(el.phi, phi);
Ex += el.Ex;
Ey += el.Ey;
Ez += el.Ez;
phi_pre96 = atan2(Ey, Ex);
double theta_pre96 = atan2(sqrt(Ex*Ex+Ey*Ey),Ez);
eta_pre96 = -log(tan(theta_pre96/2.));
return;
}
inline void Fill(double E_in, double px_in, double py_in, double pz_in, int index_in) {
double pt = sqrt(px_in*px_in+py_in*py_in);
double p = sqrt(pt*pt+pz_in*pz_in);
phi = inline_maths::phi(px_in,py_in);
double theta = asin(pt/p);
eta = inline_maths::eta(theta);
Et = E_in*sin(theta);
phi_pre96 = phi;
eta_pre96 = eta;
Ex = Et*cos(phi_pre96);
Ey = Et*sin(phi_pre96);
Ez = Et*sinh(eta_pre96);
index = index_in;
return;
}
double Ex;
double Ey;
double Ez;
int index;
double phi_pre96;
double eta_pre96;
private:
};
//end of class HepEntityIpre96;
} // end of namespace d0runi
FASTJET_END_NAMESPACE
#endif
fastjet::theta
double theta(const PseudoJet &a, const PseudoJet &b)
returns the angle between a and b
Definition PseudoJet.hh:946