Line | |
---|
1 | #include "fastjet/ClusterSequenceArea.hh" |
---|
2 | |
---|
3 | FASTJET_BEGIN_NAMESPACE |
---|
4 | |
---|
5 | LimitedWarning ClusterSequenceArea::_range_warnings; |
---|
6 | LimitedWarning ClusterSequenceArea::_explicit_ghosts_repeats_warnings; |
---|
7 | |
---|
8 | /// print a warning if the range is unsuitable for the current |
---|
9 | /// calculation of the area (e.g. because ghosts do not extend |
---|
10 | /// far enough). |
---|
11 | void ClusterSequenceArea::_warn_if_range_unsuitable(const Selector & selector) const { |
---|
12 | _check_selector_good_for_median(selector); |
---|
13 | |
---|
14 | bool no_ghosts = (_area_def.area_type() == voronoi_area) |
---|
15 | || (_area_def.area_type() == passive_area |
---|
16 | && jet_def().jet_algorithm() == kt_algorithm); |
---|
17 | if (! no_ghosts) { |
---|
18 | double rapmin, rapmax; |
---|
19 | selector.get_rapidity_extent(rapmin, rapmax); |
---|
20 | if (rapmin < -_area_def.ghost_spec().ghost_maxrap()+0.95*jet_def().R() || |
---|
21 | rapmax > _area_def.ghost_spec().ghost_maxrap()-0.95*jet_def().R()) { |
---|
22 | _range_warnings.warn("rapidity range for median (rho) extends beyond +-(ghost_maxrap - 0.95*R); this is likely to cause the results to be unreliable; safest option is to increase ghost_maxrap in the area definition"); |
---|
23 | } |
---|
24 | } |
---|
25 | } |
---|
26 | |
---|
27 | |
---|
28 | FASTJET_END_NAMESPACE |
---|
Note: See
TracBrowser
for help on using the repository browser.