1 #include "fastjet/ClusterSequeneArea.hh"
2 #include "fastjet/tools/Selector.hh"
3 #include "fastjet/tools/BackgroundEstimator"
7 using namespace fastjet;
19 vector<PseudoJet> input_particles;
20 double px, py , pz, E;
21 while (cin >> px >> py >> pz >> E) {
25 if (jets.size() > 2 && jets[0].perp() > 100) {
26 double A_J = (jets[0].perp()-jets[1].perp())/(jets[0].perp()+jets[1].perp());