source:
svn/trunk/external/fastjet/ClusterSequenceArea.cc@
1360
Last change on this file since 1360 was 859, checked in by , 12 years ago | |
---|---|
|
|
File size: 1.1 KB |
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.