fastjet::ClusterSequenceAreaBase Class Reference
base class that sets interface for extensions of ClusterSequence that provide information about the area of each jet; More...
#include <ClusterSequenceAreaBase.hh>
Public Member Functions | |
| template<class L > | |
| ClusterSequenceAreaBase (const std::vector< L > &pseudojets, const JetDefinition &jet_def, const bool &writeout_combinations=false) | |
| a constructor which just carries out the construction of the parent class | |
| ClusterSequenceAreaBase () | |
| default constructor | |
| virtual | ~ClusterSequenceAreaBase () |
| destructor | |
| virtual double | area (const PseudoJet &jet) const |
| return the area associated with the given jet; this base class returns 0. | |
| virtual double | area_error (const PseudoJet &jet) const |
| return the error (uncertainty) associated with the determination of the area of this jet; this base class returns 0. | |
| virtual PseudoJet | area_4vector (const PseudoJet &jet) const |
| return a PseudoJet whose 4-vector is defined by the following integral | |
| virtual bool | is_pure_ghost (const PseudoJet &jet) const |
| true if a jet is made exclusively of ghosts | |
| virtual bool | has_explicit_ghosts () const |
| returns true if ghosts are explicitly included within jets for this ClusterSequence; | |
| virtual double | empty_area (const RangeDefinition &range) const |
| return the total area, within range, that is free of jets, in general based on the inclusive jets | |
| double | empty_area_from_jets (const std::vector< PseudoJet > &all_jets, const RangeDefinition &range) const |
| return the total area, within range, that is free of jets, based on the supplied all_jets | |
| virtual double | n_empty_jets (const RangeDefinition &range) const |
| return something similar to the number of pure ghost jets in the given range in an active area case. | |
| double | median_pt_per_unit_area (const RangeDefinition &range) const |
| the median of (pt/area) for jets contained within range, making use also of the info on n_empty_jets | |
| double | median_pt_per_unit_area_4vector (const RangeDefinition &range) const |
| the median of (pt/area_4vector) for jets contained within making use also of the info on n_empty_jets | |
| double | median_pt_per_unit_something (const RangeDefinition &range, bool use_area_4vector) const |
the function that does the work for median_pt_per_unit_area and median_pt_per_unit_area_4vector:
| |
| virtual void | get_median_rho_and_sigma (const RangeDefinition &range, bool use_area_4vector, double &median, double &sigma, double &mean_area) const |
| using jets withing range (and with 4-vector areas if use_area_4vector), calculate the median pt/area, as well as an "error" (uncertainty), which is defined as the 1-sigma half-width of the distribution of pt/A, obtained by looking for the point below which we have (1-0.6827)/2 of the jets (including empty jets). | |
| virtual void | get_median_rho_and_sigma (const std::vector< PseudoJet > &all_jets, const RangeDefinition &range, bool use_area_4vector, double &median, double &sigma, double &mean_area, bool all_are_inclusive=false) const |
| a more advanced version of get_median_rho_and_sigma, which allows one to use any "view" of the event containing all jets (so that, e.g. | |
| virtual void | get_median_rho_and_sigma (const RangeDefinition &range, 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 | |
| virtual void | parabolic_pt_per_unit_area (double &a, double &b, const RangeDefinition &range, double exclude_above=-1.0, bool use_area_4vector=false) const |
| fits a form pt_per_unit_area(y) = a + b*y^2 in the range "range". | |
| std::vector< PseudoJet > | subtracted_jets (const double rho, const double ptmin=0.0) const |
| return a vector of all subtracted jets, using area_4vector, given rho. | |
| std::vector< PseudoJet > | subtracted_jets (const RangeDefinition &range, const double ptmin=0.0) const |
| return a vector of subtracted jets, using area_4vector. | |
| PseudoJet | subtracted_jet (const PseudoJet &jet, const double rho) const |
| return a subtracted jet, using area_4vector, given rho | |
| PseudoJet | subtracted_jet (const PseudoJet &jet, const RangeDefinition &range) const |
| return a subtracted jet, using area_4vector; note that this is potentially inefficient if repeatedly used for many different jets, because rho will be recalculated each time around. | |
| double | subtracted_pt (const PseudoJet &jet, const double rho, bool use_area_4vector=false) const |
| return the subtracted pt, given rho | |
| double | subtracted_pt (const PseudoJet &jet, const RangeDefinition &range, bool use_area_4vector=false) const |
| return the subtracted pt; note that this is potentially inefficient if repeatedly used for many different jets, because rho will be recalculated each time around. | |
Private Member Functions | |
| void | _check_jet_alg_good_for_median () const |
| check the jet algorithm is suitable (and if not issue a warning) | |
Static Private Attributes | |
| static LimitedWarning | _warnings |
| handle warning messages | |
Detailed Description
base class that sets interface for extensions of ClusterSequence that provide information about the area of each jet;the virtual functions here all return 0, since no area determination is implemented.
Definition at line 45 of file ClusterSequenceAreaBase.hh.
Constructor & Destructor Documentation
| fastjet::ClusterSequenceAreaBase::ClusterSequenceAreaBase | ( | const std::vector< L > & | pseudojets, | |
| const JetDefinition & | jet_def, | |||
| const bool & | writeout_combinations = false | |||
| ) | [inline] |
a constructor which just carries out the construction of the parent class
Definition at line 51 of file ClusterSequenceAreaBase.hh.
00053 : 00054 ClusterSequence(pseudojets, jet_def, writeout_combinations) {}
| fastjet::ClusterSequenceAreaBase::ClusterSequenceAreaBase | ( | ) | [inline] |
| virtual fastjet::ClusterSequenceAreaBase::~ClusterSequenceAreaBase | ( | ) | [inline, virtual] |
Member Function Documentation
| void fastjet::ClusterSequenceAreaBase::_check_jet_alg_good_for_median | ( | ) | const [private] |
check the jet algorithm is suitable (and if not issue a warning)
Definition at line 346 of file ClusterSequenceAreaBase.cc.
References _warnings, fastjet::cambridge_algorithm, fastjet::cambridge_for_passive_algorithm, fastjet::ClusterSequence::jet_def(), fastjet::kt_algorithm, and LimitedWarning::warn().
00346 { 00347 if (jet_def().jet_algorithm() != kt_algorithm 00348 && jet_def().jet_algorithm() != cambridge_algorithm 00349 && jet_def().jet_algorithm() != cambridge_for_passive_algorithm) { 00350 _warnings.warn("ClusterSequenceAreaBase: jet_def being used may not be suitable for estimating diffuse backgrounds (good options are kt, cam)"); 00351 } 00352 }
| virtual double fastjet::ClusterSequenceAreaBase::area | ( | const PseudoJet & | jet | ) | const [inline, virtual] |
return the area associated with the given jet; this base class returns 0.
Reimplemented in fastjet::ClusterSequenceActiveArea, fastjet::ClusterSequenceActiveAreaExplicitGhosts, fastjet::ClusterSequenceArea, and fastjet::ClusterSequenceVoronoiArea.
Definition at line 67 of file ClusterSequenceAreaBase.hh.
Referenced by empty_area_from_jets(), parabolic_pt_per_unit_area(), print_jets(), and subtracted_pt().
| virtual PseudoJet fastjet::ClusterSequenceAreaBase::area_4vector | ( | const PseudoJet & | jet | ) | const [inline, virtual] |
return a PseudoJet whose 4-vector is defined by the following integral
drap d PseudoJet("rap,phi,pt=one") * Theta("rap,phi inside jet boundary")
where PseudoJet("rap,phi,pt=one") is a 4-vector with the given rapidity (rap), azimuth (phi) and pt=1, while Theta("rap,phi inside jet boundary") is a function that is 1 when rap,phi define a direction inside the jet boundary and 0 otherwise.
This base class returns a null 4-vector.
Reimplemented in fastjet::ClusterSequenceActiveArea, fastjet::ClusterSequenceActiveAreaExplicitGhosts, fastjet::ClusterSequenceArea, and fastjet::ClusterSequenceVoronoiArea.
Definition at line 84 of file ClusterSequenceAreaBase.hh.
Referenced by parabolic_pt_per_unit_area(), and subtracted_jet().
| virtual double fastjet::ClusterSequenceAreaBase::area_error | ( | const PseudoJet & | jet | ) | const [inline, virtual] |
return the error (uncertainty) associated with the determination of the area of this jet; this base class returns 0.
Reimplemented in fastjet::ClusterSequenceActiveArea, fastjet::ClusterSequenceArea, and fastjet::ClusterSequenceVoronoiArea.
Definition at line 71 of file ClusterSequenceAreaBase.hh.
Referenced by print_jets().
| double fastjet::ClusterSequenceAreaBase::empty_area | ( | const RangeDefinition & | range | ) | const [virtual] |
return the total area, within range, that is free of jets, in general based on the inclusive jets
return the total area, within range, that is free of jets.
Calculate this as (range area) - {i in range} A_i
for ClusterSequences with explicit ghosts, assume that there will never be any empty area, i.e. it is always filled in by pure ghosts jets. This holds for seq.rec. algorithms
Reimplemented in fastjet::ClusterSequenceActiveArea, fastjet::ClusterSequenceActiveAreaExplicitGhosts, fastjet::ClusterSequenceArea, and fastjet::ClusterSequencePassiveArea.
Definition at line 54 of file ClusterSequenceAreaBase.cc.
References empty_area_from_jets(), has_explicit_ghosts(), and fastjet::ClusterSequence::inclusive_jets().
00054 { 00055 00056 if (has_explicit_ghosts()) {return 0.0;} 00057 else { return empty_area_from_jets(inclusive_jets(0.0), range);} 00058 00059 }
| double fastjet::ClusterSequenceAreaBase::empty_area_from_jets | ( | const std::vector< PseudoJet > & | all_jets, | |
| const RangeDefinition & | range | |||
| ) | const |
return the total area, within range, that is free of jets, based on the supplied all_jets
return the total area, within range, that is free of jets.
Calculate this as (range area) - {i in range} A_i
Definition at line 66 of file ClusterSequenceAreaBase.cc.
References area(), fastjet::RangeDefinition::area(), and fastjet::RangeDefinition::is_in_range().
Referenced by empty_area().
00068 { 00069 00070 double empty = range.area(); 00071 for (unsigned i = 0; i < all_jets.size(); i++) { 00072 if (range.is_in_range(all_jets[i])) empty -= area(all_jets[i]); 00073 } 00074 return empty; 00075 }
| virtual void fastjet::ClusterSequenceAreaBase::get_median_rho_and_sigma | ( | const RangeDefinition & | range, | |
| bool | use_area_4vector, | |||
| double & | median, | |||
| double & | sigma | |||
| ) | const [inline, virtual] |
same as the full version of get_median_rho_and_error, but without access to the mean_area
Reimplemented in fastjet::ClusterSequenceArea.
Definition at line 188 of file ClusterSequenceAreaBase.hh.
00190 { 00191 double mean_area; 00192 get_median_rho_and_sigma(range, use_area_4vector, 00193 median, sigma, mean_area); 00194 }
| virtual void fastjet::ClusterSequenceAreaBase::get_median_rho_and_sigma | ( | const std::vector< PseudoJet > & | all_jets, | |
| const RangeDefinition & | range, | |||
| bool | use_area_4vector, | |||
| double & | median, | |||
| double & | sigma, | |||
| double & | mean_area, | |||
| bool | all_are_inclusive = false | |||
| ) | const [virtual] |
a more advanced version of get_median_rho_and_sigma, which allows one to use any "view" of the event containing all jets (so that, e.g.
one might use Cam on a different resolution scale without have to rerun the algorithm).
By default it will assume that "all" are not inclusive jets, so that in dealing with empty area it has to calculate the number of empty jets based on the empty area and the the observed <area> of jets rather than a surmised area
Note that for small effective radii, this can cause problems because the harder jets get an area >> <ghost-jet-area> and so the estimate comes out all wrong. In these situations it is highly advisable to use an area with explicit ghosts, since then the "empty" jets are actually visible.
Reimplemented in fastjet::ClusterSequenceArea.
| void fastjet::ClusterSequenceAreaBase::get_median_rho_and_sigma | ( | const RangeDefinition & | range, | |
| bool | use_area_4vector, | |||
| double & | median, | |||
| double & | sigma, | |||
| double & | mean_area | |||
| ) | const [virtual] |
using jets withing range (and with 4-vector areas if use_area_4vector), calculate the median pt/area, as well as an "error" (uncertainty), which is defined as the 1-sigma half-width of the distribution of pt/A, obtained by looking for the point below which we have (1-0.6827)/2 of the jets (including empty jets).
The subtraction for a jet with uncorrected pt pt^U and area A is
pt^S = pt^U - median*A +- sigma*sqrt(A)
where the error is only that associated with the fluctuations in the noise and not that associated with the noise having caused changes in the hard-particle content of the jet.
NB: subtraction may also be done with 4-vector area of course, and this is recommended for jets with larger values of R, as long as rho has also been determined with a 4-vector area; using a scalar area causes one to neglect terms of relative order $R^2/8$ in the jet $p_t$.
Reimplemented in fastjet::ClusterSequenceArea.
Definition at line 156 of file ClusterSequenceAreaBase.cc.
References fastjet::ClusterSequence::inclusive_jets().
Referenced by median_pt_per_unit_something().
00158 { 00159 00160 vector<PseudoJet> incl_jets = inclusive_jets(); 00161 get_median_rho_and_sigma(incl_jets, range, use_area_4vector, 00162 median, sigma, mean_area, true); 00163 }
| virtual bool fastjet::ClusterSequenceAreaBase::has_explicit_ghosts | ( | ) | const [inline, virtual] |
returns true if ghosts are explicitly included within jets for this ClusterSequence;
Derived classes that do include explicit ghosts should provide an alternative version of this routine and set it properly.
Reimplemented in fastjet::ClusterSequenceActiveAreaExplicitGhosts, and fastjet::ClusterSequenceArea.
Definition at line 101 of file ClusterSequenceAreaBase.hh.
Referenced by empty_area().
| virtual bool fastjet::ClusterSequenceAreaBase::is_pure_ghost | ( | const PseudoJet & | jet | ) | const [inline, virtual] |
true if a jet is made exclusively of ghosts
NB: most area classes do not give any explicit ghost jets, but some do, and they should replace this function with their own version.
Reimplemented in fastjet::ClusterSequenceActiveAreaExplicitGhosts, and fastjet::ClusterSequenceArea.
Definition at line 92 of file ClusterSequenceAreaBase.hh.
| double fastjet::ClusterSequenceAreaBase::median_pt_per_unit_area | ( | const RangeDefinition & | range | ) | const |
the median of (pt/area) for jets contained within range, making use also of the info on n_empty_jets
Definition at line 77 of file ClusterSequenceAreaBase.cc.
References median_pt_per_unit_something().
Referenced by subtracted_pt().
00077 { 00078 return median_pt_per_unit_something(range,false); 00079 }
| double fastjet::ClusterSequenceAreaBase::median_pt_per_unit_area_4vector | ( | const RangeDefinition & | range | ) | const |
the median of (pt/area_4vector) for jets contained within making use also of the info on n_empty_jets
Definition at line 81 of file ClusterSequenceAreaBase.cc.
References median_pt_per_unit_something().
Referenced by subtracted_jet(), and subtracted_jets().
00081 { 00082 return median_pt_per_unit_something(range,true); 00083 }
| double fastjet::ClusterSequenceAreaBase::median_pt_per_unit_something | ( | const RangeDefinition & | range, | |
| bool | use_area_4vector | |||
| ) | const |
the function that does the work for median_pt_per_unit_area and median_pt_per_unit_area_4vector:
- something_is_area_4vect = false -> use plain area
- something_is_area_4vect = true -> use 4-vector area
the median of (pt/area) for jets contained within range, counting the empty area as if it were made up of a collection of empty jets each of area (0.55 * pi R^2).
Definition at line 90 of file ClusterSequenceAreaBase.cc.
References get_median_rho_and_sigma().
Referenced by median_pt_per_unit_area(), and median_pt_per_unit_area_4vector().
00091 { 00092 00093 double median, sigma, mean_area; 00094 get_median_rho_and_sigma(range, use_area_4vector, median, sigma, mean_area); 00095 return median; 00096 00097 }
| virtual double fastjet::ClusterSequenceAreaBase::n_empty_jets | ( | const RangeDefinition & | range | ) | const [inline, virtual] |
return something similar to the number of pure ghost jets in the given range in an active area case.
For the local implementation we return empty_area/(0.55 pi R^2), based on measured properties of ghost jets with kt and cam. Note that the number returned is a double.
Reimplemented in fastjet::ClusterSequence1GhostPassiveArea, fastjet::ClusterSequenceActiveArea, and fastjet::ClusterSequenceArea.
Definition at line 119 of file ClusterSequenceAreaBase.hh.
References fastjet::pi, and fastjet::JetDefinition::R().
00119 { 00120 double R = jet_def().R(); 00121 return empty_area(range)/(0.55*pi*R*R); 00122 }
| void fastjet::ClusterSequenceAreaBase::parabolic_pt_per_unit_area | ( | double & | a, | |
| double & | b, | |||
| const RangeDefinition & | range, | |||
| double | exclude_above = -1.0, |
|||
| bool | use_area_4vector = false | |||
| ) | const [virtual] |
fits a form pt_per_unit_area(y) = a + b*y^2 in the range "range".
fits a form pt_per_unit_area(y) = a + b*y^2 for jets in range.
exclude_above allows one to exclude large values of pt/area from fit. (if negative, the cut is discarded) use_area_4vector = true uses the 4vector areas.
exclude_above allows one to exclude large values of pt/area from fit. use_area_4vector = true uses the 4vector areas.
Reimplemented in fastjet::ClusterSequenceArea.
Definition at line 104 of file ClusterSequenceAreaBase.cc.
References area(), area_4vector(), fastjet::ClusterSequence::inclusive_jets(), fastjet::RangeDefinition::is_in_range(), and fastjet::PseudoJet::perp().
00106 { 00107 00108 int n=0; 00109 int n_excluded = 0; 00110 double mean_f=0, mean_x2=0, mean_x4=0, mean_fx2=0; 00111 00112 vector<PseudoJet> incl_jets = inclusive_jets(); 00113 00114 for (unsigned i = 0; i < incl_jets.size(); i++) { 00115 if (range.is_in_range(incl_jets[i])) { 00116 double this_area; 00117 if ( use_area_4vector ) { 00118 this_area = area_4vector(incl_jets[i]).perp(); 00119 } else { 00120 this_area = area(incl_jets[i]); 00121 } 00122 double f = incl_jets[i].perp()/this_area; 00123 if (exclude_above <= 0.0 || f < exclude_above) { 00124 double x = incl_jets[i].rap(); double x2 = x*x; 00125 mean_f += f; 00126 mean_x2 += x2; 00127 mean_x4 += x2*x2; 00128 mean_fx2 += f*x2; 00129 n++; 00130 } else { 00131 n_excluded++; 00132 } 00133 } 00134 } 00135 00136 if (n <= 1) { 00137 // meaningful results require at least two jets inside the 00138 // area -- mind you if there are empty jets we should be in 00139 // any case doing something special... 00140 a = 0.0; 00141 b = 0.0; 00142 } else { 00143 mean_f /= n; 00144 mean_x2 /= n; 00145 mean_x4 /= n; 00146 mean_fx2 /= n; 00147 00148 b = (mean_f*mean_x2 - mean_fx2)/(mean_x2*mean_x2 - mean_x4); 00149 a = mean_f - b*mean_x2; 00150 } 00151 //cerr << "n_excluded = "<< n_excluded << endl; 00152 }
| PseudoJet fastjet::ClusterSequenceAreaBase::subtracted_jet | ( | const PseudoJet & | jet, | |
| const RangeDefinition & | range | |||
| ) | const |
return a subtracted jet, using area_4vector; note that this is potentially inefficient if repeatedly used for many different jets, because rho will be recalculated each time around.
return a subtracted jet, using area_4vector; note that this is potentially inefficient if repeatedly used for many different jets, because rho will be recalculated each time around.
Definition at line 308 of file ClusterSequenceAreaBase.cc.
References median_pt_per_unit_area_4vector(), and subtracted_jet().
00309 { 00310 double rho = median_pt_per_unit_area_4vector(range); 00311 PseudoJet sub_jet = subtracted_jet(jet, rho); 00312 return sub_jet; 00313 }
| PseudoJet fastjet::ClusterSequenceAreaBase::subtracted_jet | ( | const PseudoJet & | jet, | |
| const double | rho | |||
| ) | const |
return a subtracted jet, using area_4vector, given rho
Definition at line 287 of file ClusterSequenceAreaBase.cc.
References area_4vector(), fastjet::PseudoJet::cluster_hist_index(), fastjet::PseudoJet::perp(), fastjet::PseudoJet::set_cluster_hist_index(), fastjet::PseudoJet::set_user_index(), and fastjet::PseudoJet::user_index().
Referenced by subtracted_jet(), subtracted_jets(), and subtracted_pt().
00288 { 00289 PseudoJet area4vect = area_4vector(jet); 00290 PseudoJet sub_jet; 00291 // sanity check 00292 if (rho*area4vect.perp() < jet.perp() ) { 00293 sub_jet = jet - rho*area4vect; 00294 } else { sub_jet = PseudoJet(0.0,0.0,0.0,0.0); } 00295 00296 // make sure the subtracted jet has the same index (cluster and user) 00297 // (i.e. "looks like") the original jet 00298 sub_jet.set_cluster_hist_index(jet.cluster_hist_index()); 00299 sub_jet.set_user_index(jet.user_index()); 00300 00301 return sub_jet; 00302 }
| vector< PseudoJet > fastjet::ClusterSequenceAreaBase::subtracted_jets | ( | const RangeDefinition & | range, | |
| const double | ptmin = 0.0 | |||
| ) | const |
return a vector of subtracted jets, using area_4vector.
Only inclusive_jets above ptmin are subtracted and returned. the ordering is the same as that of sorted_by_pt(cs.inclusive_jets()), i.e. not necessarily ordered in pt once subtracted
Only inclusive_jets above ptmin are subtracted and returned. the ordering is the same as that of sorted_by_pt(cs.inclusive_jets()), i.e. not necessarily ordered in pt once subtracted
Definition at line 277 of file ClusterSequenceAreaBase.cc.
References median_pt_per_unit_area_4vector(), and subtracted_jets().
00280 { 00281 double rho = median_pt_per_unit_area_4vector(range); 00282 return subtracted_jets(rho,ptmin); 00283 }
| vector< PseudoJet > fastjet::ClusterSequenceAreaBase::subtracted_jets | ( | const double | rho, | |
| const double | ptmin = 0.0 | |||
| ) | const |
return a vector of all subtracted jets, using area_4vector, given rho.
Only inclusive_jets above ptmin are subtracted and returned. the ordering is the same as that of sorted_by_pt(cs.inclusive_jets()), i.e. not necessarily ordered in pt once subtracted
Only inclusive_jets above ptmin are subtracted and returned. the ordering is the same as that of sorted_by_pt(cs.inclusive_jets()), i.e. not necessarily ordered in pt once subtracted
Definition at line 261 of file ClusterSequenceAreaBase.cc.
References fastjet::ClusterSequence::inclusive_jets(), fastjet::ClusterSequence::jets(), fastjet::sorted_by_pt(), and subtracted_jet().
Referenced by subtracted_jets().
00263 { 00264 vector<PseudoJet> sub_jets; 00265 vector<PseudoJet> jets = sorted_by_pt(inclusive_jets(ptmin)); 00266 for (unsigned i=0; i<jets.size(); i++) { 00267 PseudoJet sub_jet = subtracted_jet(jets[i],rho); 00268 sub_jets.push_back(sub_jet); 00269 } 00270 return sub_jets; 00271 }
| double fastjet::ClusterSequenceAreaBase::subtracted_pt | ( | const PseudoJet & | jet, | |
| const RangeDefinition & | range, | |||
| bool | use_area_4vector = false | |||
| ) | const |
return the subtracted pt; note that this is potentially inefficient if repeatedly used for many different jets, because rho will be recalculated each time around.
return the subtracted pt; note that this is potentially inefficient if repeatedly used for many different jets, because rho will be recalculated each time around.
Definition at line 332 of file ClusterSequenceAreaBase.cc.
References median_pt_per_unit_area(), fastjet::PseudoJet::perp(), subtracted_jet(), and subtracted_pt().
00334 { 00335 if ( use_area_4vector ) { 00336 PseudoJet sub_jet = subtracted_jet(jet,range); 00337 return sub_jet.perp(); 00338 } else { 00339 double rho = median_pt_per_unit_area(range); 00340 return subtracted_pt(jet,rho,false); 00341 } 00342 }
| double fastjet::ClusterSequenceAreaBase::subtracted_pt | ( | const PseudoJet & | jet, | |
| const double | rho, | |||
| bool | use_area_4vector = false | |||
| ) | const |
return the subtracted pt, given rho
Definition at line 317 of file ClusterSequenceAreaBase.cc.
References area(), fastjet::PseudoJet::perp(), and subtracted_jet().
Referenced by subtracted_pt().
00319 { 00320 if ( use_area_4vector ) { 00321 PseudoJet sub_jet = subtracted_jet(jet,rho); 00322 return sub_jet.perp(); 00323 } else { 00324 return jet.perp() - rho*area(jet); 00325 } 00326 }
Member Data Documentation
LimitedWarning fastjet::ClusterSequenceAreaBase::_warnings [static, private] |
handle warning messages
allow for warnings
Reimplemented in fastjet::ClusterSequenceActiveAreaExplicitGhosts.
Definition at line 246 of file ClusterSequenceAreaBase.hh.
Referenced by _check_jet_alg_good_for_median().
The documentation for this class was generated from the following files:
- include/fastjet/ClusterSequenceAreaBase.hh
- src/ClusterSequenceAreaBase.cc
