/usr/include/fastjet/CircularRange.hh is in libfastjet-dev 3.0.6+dfsg-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 | //STARTHEADER
// $Id: CircularRange.hh 2577 2011-09-13 15:11:38Z salam $
//
// Copyright (c) 2005-2011, Matteo Cacciari, Gavin P. Salam and Gregory Soyez
//
//----------------------------------------------------------------------
// This file is part of FastJet.
//
// FastJet is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// The algorithms that underlie FastJet have required considerable
// development and are described in hep-ph/0512210. If you use
// FastJet as part of work towards a scientific publication, please
// include a citation to the FastJet paper.
//
// FastJet is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with FastJet. If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------
//ENDHEADER
#ifndef __FASTJET_CIRCULARRANGE_HH__
#define __FASTJET_CIRCULARRANGE_HH__
#include "fastjet/RangeDefinition.hh"
#include "fastjet/Error.hh"
// for backwards compatibility: one should now use SelectorCircle,
// defined in fastjet/Selector.hh, instead CircularRange
#warning This file includes fastjet/CircularRange.hh, \
a deprecated FastJet header provided only for backward compatibility. \
This is not guaranteed to work in future releases of FastJet. \
From FastJet 3.0 onwards, please consider using Selector, defined in \
fastjet/Selector.hh, instead of RangeDefinition and, in particular, \
SelectorCircle instead of CircularRange.
FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh
class CircularRange : public fastjet::RangeDefinition {
public:
/// constructor
CircularRange() {_set_invalid_rapphi();}
/// initialise CircularRange with a jet
CircularRange(const fastjet::PseudoJet & jet, double distance) {
_distance = distance;
_rapjet = jet.rap();
_phijet = jet.phi();
_total_area = fastjet::pi*_distance*_distance; }
/// initialise CircularRange with a (rap,phi) point
CircularRange(double rap, double phi, double distance) {
_distance = distance;
_rapjet = rap;
_phijet = phi;
_total_area = fastjet::pi*_distance*_distance; }
/// initialise CircularRange with just the radius parameter
CircularRange(double distance) {
_set_invalid_rapphi();
_distance = distance;
_total_area = fastjet::pi*_distance*_distance; }
/// destructor
virtual ~CircularRange() {}
/// return description of range
virtual inline std::string description() const {
std::ostringstream ostr;
ostr << "CircularRange: within distance "<< _distance << " of given jet or point." ;
return ostr.str(); }
/// returns true since this range is localizable (i.e. set_position
/// does something meaningful)
virtual inline bool is_localizable() const { return true; }
/// return bool according to whether (rap,phi) is in range
virtual inline bool is_in_range(double rap, double phi) const {
if (! _rapphi_are_valid()) {
throw Error("Circular range used without a center having being defined (use set_position())");
}
double deltaphi = _phijet - phi;
if ( deltaphi > pi) { deltaphi -= twopi; }
else if ( deltaphi < -pi) { deltaphi += twopi; }
bool inrange = ( (rap-_rapjet)*(rap-_rapjet) +
deltaphi*deltaphi <= _distance*_distance );
return inrange; }
/// return the minimal and maximal rapidity of this range
virtual inline void get_rap_limits(double & rapmin, double & rapmax) const {
rapmin = _rapjet - _distance;
rapmax = _rapjet + _distance; }
private:
double _distance;
/// value for phi that marks it as invalid
const static double _invalid_phi = -1000.0;
/// set internal phi so as to mark things as invalid
void _set_invalid_rapphi() {_phijet = _invalid_phi;}
/// true if rap,phi are valid (tests only phi)
bool _rapphi_are_valid() const {return _phijet != _invalid_phi;}
};
FASTJET_END_NAMESPACE
#endif // __FASTJET_CIRCULARRANGE_HH__
|