30 #include "fastjet/ClusterSequence.hh"
31 #include "fastjet/JadePlugin.hh"
34 #include "fastjet/NNH.hh"
46 FASTJET_BEGIN_NAMESPACE
53 void init(
const PseudoJet & jet) {
54 double norm = 1.0/sqrt(jet.modp2());
58 rt2E = sqrt(2.0)*jet.E();
61 double distance(
const JadeBriefJet * jet)
const {
62 double dij = 1 - nx*jet->nx
65 dij *= rt2E*jet->rt2E;
69 double beam_distance()
const {
70 return numeric_limits<double>::max();
74 double rt2E, nx, ny, nz;
79 string JadePlugin::description ()
const {
81 desc <<
"e+e- JADE algorithm plugin";
87 int njets = cs.
jets().size();
102 double diB = cs.
jets()[i].E()*cs.
jets()[i].E();
110 FASTJET_END_NAMESPACE
double norm(const VPoint p)
norm of a vector
void plugin_record_iB_recombination(int jet_i, double diB)
record the fact that there has been a recombination between jets()[jet_i] and the beam,...
double dij_min(int &iA, int &iB)
return the dij_min and indices iA, iB, for the corresponding jets.
const std::vector< PseudoJet > & jets() const
allow the user to access the internally stored _jets() array, which contains both the initial particl...
void remove_jet(int iA)
remove the jet pointed to by index iA
void plugin_record_ij_recombination(int jet_i, int jet_j, double dij, int &newjet_k)
record the fact that there has been a recombination between jets()[jet_i] and jets()[jet_k],...
void merge_jets(int iA, int iB, const PseudoJet &jet, int jet_index)
merge the jets pointed to by indices A and B and replace them with jet, assigning it an index jet_ind...