plugin = new fastjet::SISConePlugin(coneRadius, overlapThreshold, nPassMax, ptProtoJetMin, caching);
vector<fastjet::PseudoJet> inputParticles;
+
//******************************** READING OF INPUT PARTICLES
// Here we look for px, py pz and energy of each particle that we gather in a PseudoJet object, and we put all these PseudoJet in a vector of PseudoJets : input_particles.
+ cout<<"Only the 1st (hardest) jet of the event hasn't been taken into account for rho estimation"<<endl;
+ }
+ if (inclusiveJetsForRhoEstimation.size() > (UInt_t)NHardJets) inclusiveJetsForRhoEstimation.erase(inclusiveJetsForRhoEstimation.begin(), inclusiveJetsForRhoEstimation.begin()+NHardJets);
+
+ // Estimation of rho and fluctuations sigma
+
+ //1st method
+ Double_t rho = 0;
+ Double_t sigma = 0;
+ Double_t meanarea = 0;
+ // 3rd argument : 1 = use area 4 vector rather than simple area
+ // last argument : 0 = in case of explicit ghosts use
+ if (inclusiveJetsForRhoEstimation.size() != 0) csForRho.get_median_rho_and_sigma(inclusiveJetsForRhoEstimation, range, 1, rho, sigma, meanarea, 1); // this gives the fluctuations also