Fork me on GitHub

source: git/external/fastjet/ClusterSequenceArea.cc@ 41376ac

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 41376ac was d7d2da3, checked in by pavel <pavel@…>, 12 years ago

move branches/ModularDelphes to trunk

  • Property mode set to 100644
File size: 1.1 KB
Line 
1#include "fastjet/ClusterSequenceArea.hh"
2
3FASTJET_BEGIN_NAMESPACE
4
5LimitedWarning ClusterSequenceArea::_range_warnings;
6LimitedWarning 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).
11void 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
28FASTJET_END_NAMESPACE
Note: See TracBrowser for help on using the repository browser.