7 #include "fastjet/ClusterSequence.hh"
20 std::vector<float> s_e;
21 std::vector<float> s_px;
22 std::vector<float> s_py;
23 std::vector<float> s_pz;
31 }
else if (
m_type ==
"Kt") {
33 }
else if (
m_type ==
"CamKt") {
37 fastjet::JetDefinition microjet_def(jetalg,
m_R);
39 std::vector<fastjet::PseudoJet> microjets = fastjet::sorted_by_pt(microjet_cs.inclusive_jets(
m_minPt));
40 for (
size_t z = 0;
z < microjets.size(); ++
z) {
41 s_e.push_back(microjets[
z].
e());
42 s_px.push_back(microjets[
z].
px());
43 s_py.push_back(microjets[
z].
py());
44 s_pz.push_back(microjets[
z].
pz());
50 jet.setAttribute(
ss.str().c_str(), s_e);
53 jet.setAttribute(
ss.str().c_str(), s_px);
56 jet.setAttribute(
ss.str().c_str(), s_py);
59 jet.setAttribute(
ss.str().c_str(), s_pz);