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;
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__ virtual ~ClusterSequenceAreaBase()
destructor
virtual double n_empty_jets(const Selector &selector) const
return something similar to the number of pure ghost jets in the given selector's range in an active ...
virtual double area(const PseudoJet &) const
return the area associated with the given jet; this base class returns 0.
class to provide facilities for giving warnings up to some maximum number of times and to provide glo...
virtual PseudoJet area_4vector(const PseudoJet &) const
return a PseudoJet whose 4-vector is defined by the following integral
ClusterSequenceAreaBase()
default constructor
base class that sets interface for extensions of ClusterSequence that provide information about the a...
virtual bool has_explicit_ghosts() const
returns true if ghosts are explicitly included within jets for this ClusterSequence; ...
Class that encodes information about cuts and other selection criteria that can be applied to PseudoJ...
virtual double area_error(const PseudoJet &) const
return the error (uncertainty) associated with the determination of the area of this jet; this base c...
virtual void get_median_rho_and_sigma(const Selector &selector, bool use_area_4vector, double &median, double &sigma) const
same as the full version of get_median_rho_and_error, but without access to the mean_area ...
Class to contain pseudojets, including minimal information of use to jet-clustering routines...
class that is intended to hold a full definition of the jet clusterer
virtual bool is_pure_ghost(const PseudoJet &) const
true if a jet is made exclusively of ghosts