34 #include "fastjet/internal/ClusterSequence_N2.icc"
38 FASTJET_BEGIN_NAMESPACE
52 template<>
inline void ClusterSequence::_bj_set_jetinfo(
53 EEBriefJet *
const jetA,
const int _jets_index)
const {
55 double E = _jets[_jets_index].E();
57 double p = jet_def().extra_param();
58 switch (_jet_algorithm) {
60 assert(_Rparam > 2.0);
66 if (p <= 0 && scale < 1e-300) scale = 1e-300;
70 throw Error(
"Unrecognised jet algorithm");
74 double norm = _jets[_jets_index].modp2();
76 norm = 1.0/sqrt(norm);
77 jetA->nx =
norm * _jets[_jets_index].px();
78 jetA->ny =
norm * _jets[_jets_index].py();
79 jetA->nz =
norm * _jets[_jets_index].pz();
85 jetA->_jets_index = _jets_index;
93 template<>
double ClusterSequence::_bj_dist(
94 const EEBriefJet *
const jeta,
95 const EEBriefJet *
const jetb)
const {
108 void ClusterSequence::_simple_N2_cluster_BriefJet() {
109 _simple_N2_cluster<BriefJet>();
114 void ClusterSequence::_simple_N2_cluster_EEBriefJet() {
115 _simple_N2_cluster<EEBriefJet>();
127 FASTJET_END_NAMESPACE