00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #ifndef __FASTJET_CLUSTERSEQUENCEAREABASE_HH__
00032 #define __FASTJET_CLUSTERSEQUENCEAREABASE_HH__
00033
00034 #include "fastjet/ClusterSequence.hh"
00035 #include "fastjet/internal/LimitedWarning.hh"
00036 #include "fastjet/RangeDefinition.hh"
00037
00038 FASTJET_BEGIN_NAMESPACE
00039
00045 class ClusterSequenceAreaBase : public ClusterSequence {
00046 public:
00047
00050 template<class L> ClusterSequenceAreaBase
00051 (const std::vector<L> & pseudojets,
00052 const JetDefinition & jet_def,
00053 const bool & writeout_combinations = false) :
00054 ClusterSequence(pseudojets, jet_def, writeout_combinations) {}
00055
00056
00058 ClusterSequenceAreaBase() {}
00059
00060
00062 virtual ~ClusterSequenceAreaBase() {}
00063
00064
00067 virtual double area (const PseudoJet & jet) const {return 0.0;}
00068
00071 virtual double area_error (const PseudoJet & jet) const {return 0.0;}
00072
00084 virtual PseudoJet area_4vector(const PseudoJet & jet) const {
00085 return PseudoJet(0.0,0.0,0.0,0.0);}
00086
00092 virtual bool is_pure_ghost(const PseudoJet & jet) const {
00093 return false;
00094 }
00095
00097 virtual double empty_area(const RangeDefinition & range) const;
00098
00104 virtual double n_empty_jets(const RangeDefinition & range) const {
00105 double R = jet_def().R();
00106 return empty_area(range)/(0.55*pi*R*R);
00107 }
00108
00111 double median_pt_per_unit_area(const RangeDefinition & range) const;
00112
00115 double median_pt_per_unit_area_4vector(const RangeDefinition & range) const;
00116
00121 double median_pt_per_unit_something(
00122 const RangeDefinition & range, bool use_area_4vector) const;
00123
00140 virtual void get_median_rho_and_sigma(const RangeDefinition & range,
00141 bool use_area_4vector,
00142 double & median, double & sigma,
00143 double & mean_area) const;
00144
00147 void get_median_rho_and_sigma(const RangeDefinition & range,
00148 bool use_area_4vector,
00149 double & median, double & sigma) const {
00150 double mean_area;
00151 get_median_rho_and_sigma(range, use_area_4vector,
00152 median, sigma, mean_area);
00153 }
00154
00155
00159 virtual void parabolic_pt_per_unit_area(double & a, double & b,
00160 const RangeDefinition & range,
00161 double exclude_above=-1.0,
00162 bool use_area_4vector=false) const;
00163
00164
00169 std::vector<PseudoJet> subtracted_jets(const double rho,
00170 const double ptmin=0.0) const;
00171
00176 std::vector<PseudoJet> subtracted_jets(const RangeDefinition & range,
00177 const double ptmin=0.0) const;
00178
00180 PseudoJet subtracted_jet(const PseudoJet & jet,
00181 const double rho) const;
00182
00187 PseudoJet subtracted_jet(const PseudoJet & jet,
00188 const RangeDefinition & range) const;
00189
00191 double subtracted_pt(const PseudoJet & jet,
00192 const double rho,
00193 bool use_area_4vector=false) const;
00194
00198 double subtracted_pt(const PseudoJet & jet,
00199 const RangeDefinition & range,
00200 bool use_area_4vector=false) const;
00201
00202
00203 private:
00205 static LimitedWarning _warnings;
00206
00208 void _check_jet_alg_good_for_median() const;
00209
00210 };
00211
00212
00213
00214 FASTJET_END_NAMESPACE
00215
00216 #endif // __FASTJET_CLUSTERSEQUENCEAREABASE_HH__