30 #include "fastjet/ClusterSequence.hh"
31 #include "fastjet/EECambridgePlugin.hh"
32 #include "fastjet/NNH.hh"
38 FASTJET_BEGIN_NAMESPACE
46 void init(
const PseudoJet & jet) {
47 double norm = 1.0/sqrt(jet.modp2());
53 double distance(
const EECamBriefJet * jet)
const {
54 double dij = 1 - nx*jet->nx
60 double beam_distance()
const {
61 return numeric_limits<double>::max();
69 string EECambridgePlugin::description ()
const {
71 desc <<
"EECambridge plugin with ycut = " << ycut() ;
76 int njets = cs.
jets().size();
93 double scale = min(cs.
jets()[i].E(), cs.
jets()[j].E());
94 dij = 2 * vij * scale * scale;
95 if (dij > Q2 * ycut()) {
97 if (cs.
jets()[i].E() > cs.
jets()[j].E()) std::swap(i,j);
117 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
double Q2() const
return Q()^2
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...