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();
100 nnh.merge_jets(i, j, cs.
jets()[k], k);
102 double diB = cs.
jets()[i].E()*cs.
jets()[i].E();
110 FASTJET_END_NAMESPACE
double dij_min(int &iA, int &iB)
return the dij_min and indices iA, iB, for the corresponding jets.
double norm(const VPoint p)
norm of a vector
Help solve closest pair problems with generic interparticle and beam distance.
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], with the specified dij, and return the index (newjet_k) allocated to the new jet, whose momentum is assumed to be the 4-vector sum of that of jet_i and jet_j
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...
const std::vector< PseudoJet > & jets() const
allow the user to access the internally stored _jets() array, which contains both the initial particl...