| 62 |
|
fastjet::Nsubjettiness Nsubj(N, mode, beta, R0, Rcutoff); |
| 63 |
|
|
| 64 |
|
double sub=Nsubj.result(fjet); |
| 65 |
– |
if (sub < -0.5){ |
| 66 |
– |
cout << "Nsub jettiness " << N << " = " << sub << endl; |
| 67 |
– |
cout << "jet: pt = " << m_jet->pt() << " phi = " << m_jet->phi() << " eta = " << m_jet->eta() << " N = " << m_jet->pfconstituents_indices().size() << endl; |
| 68 |
– |
for (unsigned int i=0; i<FJparticles.size(); ++i){ |
| 69 |
– |
fastjet::PseudoJet p = FJparticles[i]; |
| 70 |
– |
cout << "part " << i << ": pt = " << p.pt() << " phi = " << p.phi() << " eta = " << p.pseudorapidity() << endl; |
| 71 |
– |
} |
| 72 |
– |
} |
| 73 |
– |
|
| 65 |
|
|
| 66 |
|
return sub; |
| 67 |
|
|
