This file is indexed.

/usr/include/siscone/area.h is in libsiscone-dev 2.0.6-1.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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
// -*- C++ -*-
///////////////////////////////////////////////////////////////////////////////
// File: area.h                                                              //
// Description: header file for the computation of jet area                  //
// This file is part of the SISCone project.                                 //
// For more details, see http://projects.hepforge.org/siscone                //
//                                                                           //
// Copyright (c) 2006 Gavin Salam and Gregory Soyez                          //
//                                                                           //
// This program 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.                                       //
//                                                                           //
// This program 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 this program; if not, write to the Free Software               //
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA //
//                                                                           //
// $Revision:: 149                                                          $//
// $Date:: 2007-03-15 00:13:58 +0100 (Thu, 15 Mar 2007)                     $//
///////////////////////////////////////////////////////////////////////////////

#ifndef __SISCONE_AREA_H__
#define __SISCONE_AREA_H__

#include "momentum.h"
#include "siscone.h"

namespace siscone{

/**
 * \class Cjet_area
 * real Jet information, including its area(s)
 *
 * This class contains information for one single jet. 
 * That is, first, its momentum carrying information
 * about its centre and pT, and second, its particle
 * contents.
 * Compared to the Cjet class, it also includes the 
 * passive and active areas of the jet computed using 
 * the Carea class.
 */
class Cjet_area : public Cjet{
 public:
  /// default ctor
  Cjet_area();

  /// jet-initialised ctor
  Cjet_area(Cjet &j);

  /// default dtor
  ~Cjet_area();

  // area information
  double passive_area;   ///< passive area
  double active_area;    ///< active area
};

/**
 * \class Carea
 * class for the computation of jet areas.
 * 
 * This is the class user should use whenever you want to compute
 * the jet area (passive and active). .
 * It uses the SISCone algorithm to perform the jet analysis.
 */
class Carea : public Csiscone{
 public:
  /// default ctor
  Carea();

  /// default dtor
  ~Carea();

  /**
   * compute the jet areas from a given particle set.
   * The parameters of this method are the ones which control the jet clustering alghorithn.
   * Note that the pt_min is not allowed here soince the jet-area determination involves soft 
   * particles/jets and thus is used internally.
   * \param _particles   list of particles
   * \param _radius      cone radius
   * \param _f           shared energy threshold for splitting&merging
   * \param _n_pass_max  maximum number of passes (0=full search, the default)
   * \param _split_merge_scale    the scale choice for the split-merge procedure
   *        NOTE: SM_pt leads to IR unsafety for some events with momentum conservation. 
   *              SM_Et is IR safe but not boost invariant and not implemented(!)
   *              SM_mt is IR safe for hadronic events, but not for decays of two 
   *                    back-to-back particles of identical mass
   *              SM_pttilde  
   *                    is always IR safe, and also boost invariant (default)
   * \param _hard_only  when this is set on, only hard jets are computed
   *                    and not the purely ghosted jets (default: false)
   * \return the number of jets (including pure-ghost ones if they are included)
   */
  int compute_areas(std::vector<Cmomentum> &_particles, double _radius, double _f, 
		    int _n_pass_max=0, Esplit_merge_scale _split_merge_scale=SM_pttilde,
		    bool _hard_only=false);

  /**
   * compute the jet active areas from a given particle set.
   * The parameters of this method are the ones which control the jet clustering alghorithn.
   * Note that the pt_min is not allowed here soince the jet-area determination involves soft 
   * particles/jets and thus is used internally.
   * See compute_areas for paramters definition.
   * \return the number of jets (including pure-ghost ones if they are included)
   */
  int compute_active_areas(std::vector<Cmomentum> &_particles, double _radius, double _f, 
			   int _n_pass_max=0, Esplit_merge_scale _split_merge_scale=SM_pttilde,
			   bool _hard_only=false);

  /**
   * compute the jet passive areas from a given particle set.
   * The parameters of this method are the ones which control the jet clustering alghorithn.
   * Note that the pt_min is not allowed here soince the jet-area determination involves soft 
   * particles/jets and thus is used internally.
   * See compute_areas for paramters definition.
   */
  int compute_passive_areas(std::vector<Cmomentum> &_particles, double _radius, double _f, 
			    int _n_pass_max=0, Esplit_merge_scale _split_merge_scale=SM_pttilde);

  int grid_size;        ///< size of the grid we add soft particles on (N_soft=(grid_size^2))
  double grid_eta_max;  ///< maximal value of eta we add soft particles on
  double grid_shift;    ///< fractional (random) displacement of the points om the grid

  double pt_soft;       ///< pt of the soft particles added
  double pt_shift;      ///< amplitude of the pt random shift
  double pt_soft_min;   ///< pt_min used in SM to compute passive areas

  /// jets with their areas
  std::vector<Cjet_area> jet_areas;
};

}
#endif