99#include "TRandom3.h"
1010#include "TSystem.h"
1111#include "fairlogger/Logger.h"
12- #include "fastjet/ClusterSequence.hh"
12+ #include <fastjet/AreaDefinition.hh>
13+ #include <fastjet/ClusterSequence.hh>
14+ #include <fastjet/ClusterSequenceArea.hh>
15+ #include <fastjet/GhostedAreaSpec.hh>
16+ #include <fastjet/PseudoJet.hh>
17+ #include <fastjet/Selector.hh>
18+ #include <fastjet/tools/JetMedianBackgroundEstimator.hh>
19+ #include <fastjet/tools/Subtractor.hh>
1320#include <cmath>
1421#include <fstream>
1522#include <string>
@@ -113,7 +120,7 @@ protected:
113120 const std ::vector < int > pdgXiOmega = {3312 , -3312 , 3334 , -3334 };
114121 const double mpi = 0.1395704 ;
115122
116- std ::vector < PseudoJet > fjParticles ;
123+ std ::vector < fastjet :: PseudoJet > fjParticles ;
117124
118125 for (int i = 0 ; i < event .size (); ++ i ) {
119126 const auto& p = event [i ];
@@ -129,16 +136,16 @@ protected:
129136
130137 double energy = std ::sqrt (p .p () * p .p () + mpi * mpi );
131138
132- PseudoJet pj (p .px (), p .py (), p .pz (), energy );
139+ fastjet :: PseudoJet pj (p .px (), p .py (), p .pz (), energy );
133140 pj .set_user_index (i ); // map back to Pythia index
134141 fjParticles .push_back (pj );
135142 }
136143
137144 if (fjParticles .empty ()) return false;
138145
139- JetDefinition jetDef (antikt_algorithm , mJetR );
140- ClusterSequence cs (fjParticles , jetDef );
141- auto jets = sorted_by_pt (cs .inclusive_jets (mPtJetThreshold ));
146+ fastjet :: JetDefinition jetDef (fastjet :: antikt_algorithm , mJetR );
147+ fastjet :: ClusterSequence cs (fjParticles , jetDef );
148+ auto jets = fastjet :: sorted_by_pt (cs .inclusive_jets (mPtJetThreshold ));
142149
143150 for (const auto& jet : jets ) {
144151 for (const auto& c : jet .constituents ()) {
0 commit comments