29 #ifndef __FASTJET_CLUSTERSEQUENCEAREABASE_HH__
30 #define __FASTJET_CLUSTERSEQUENCEAREABASE_HH__
32 #include "fastjet/ClusterSequence.hh"
33 #include "fastjet/LimitedWarning.hh"
34 #include "fastjet/Selector.hh"
36 FASTJET_BEGIN_NAMESPACE
51 (
const std::vector<L> & pseudojets,
53 const bool & writeout_combinations =
false) :
111 virtual double empty_area(
const Selector & selector)
const;
119 double empty_area_from_jets(
const std::vector<PseudoJet> & all_jets,
134 double R = jet_def().R();
135 return empty_area(selector)/(0.55*pi*R*R);
144 double median_pt_per_unit_area(
const Selector & selector)
const;
151 double median_pt_per_unit_area_4vector(
const Selector & selector)
const;
157 double median_pt_per_unit_something(
158 const Selector & selector,
bool use_area_4vector)
const;
184 virtual void get_median_rho_and_sigma(
const Selector & selector,
185 bool use_area_4vector,
186 double & median,
double & sigma,
187 double & mean_area)
const;
208 virtual void get_median_rho_and_sigma(
const std::vector<PseudoJet> & all_jets,
210 bool use_area_4vector,
211 double & median,
double & sigma,
213 bool all_are_inclusive =
false)
const;
221 virtual void get_median_rho_and_sigma(
const Selector & selector,
222 bool use_area_4vector,
223 double & median,
double & sigma)
const {
225 get_median_rho_and_sigma(selector, use_area_4vector,
226 median, sigma, mean_area);
238 virtual void parabolic_pt_per_unit_area(
double & a,
double & b,
240 double exclude_above=-1.0,
241 bool use_area_4vector=
false)
const;
247 std::vector<PseudoJet> subtracted_jets(
const double rho,
248 const double ptmin=0.0)
const;
258 std::vector<PseudoJet> subtracted_jets(
const Selector & selector,
259 const double ptmin=0.0)
const;
263 const double rho)
const;
277 double subtracted_pt(
const PseudoJet & jet,
279 bool use_area_4vector=
false)
const;
288 double subtracted_pt(
const PseudoJet & jet,
290 bool use_area_4vector=
false)
const;
294 void _check_selector_good_for_median(
const Selector &selector)
const;
304 void _check_jet_alg_good_for_median()
const;
310 FASTJET_END_NAMESPACE
312 #endif // __FASTJET_CLUSTERSEQUENCEAREABASE_HH__