-
- //fastjet::ClusterSequence clust_seq(inputParticles, jetDef);
- fastjet::ClusterSequenceArea clust_seq(inputParticles, jetDef, areaDef);
-
- // save a comment in the header
-
- TString comment = "Running FastJet algorithm with the following setup. ";
- comment+= "Jet definition: ";
- comment+= TString(jetDef.description());
- comment+= ". Strategy adopted by FastJet: ";
- comment+= TString(clust_seq.strategy_string());
- header->SetComment(comment);
- if(debug){
- cout << "--------------------------------------------------------" << endl;
- cout << comment << endl;
- cout << "--------------------------------------------------------" << endl;
- }
- //header->PrintParameters();
-
- // extract the inclusive jets with pt > ptmin, sorted by pt
- double ptmin = header->GetPtMin();
- vector<fastjet::PseudoJet> inclusiveJets = clust_seq.inclusive_jets(ptmin);
-
- //cout << "Number of unclustered particles: " << clust_seq.unclustered_particles().size() << endl;
-
- vector<fastjet::PseudoJet> jets = sorted_by_pt(inclusiveJets); // Added by me
- for (size_t j = 0; j < jets.size(); j++) { // loop for jets // Added by me
-
- if(debug) printf("Jet found %5d %9.5f %8.5f %10.3f \n",(Int_t)j,jets[j].rap(),jets[j].phi(),jets[j].perp());
+ if(debug) printf("Jet found %5d %9.5f %8.5f %10.3f %8.3f +- %6.3f\n", (Int_t)j,jets[j].rap(),jets[j].phi(),jets[j].perp(), area, areaError);