FastJet  3.0.6
ClusterSequenceArea.cc
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
the longitudinally invariant kt algorithm