area.h

// -*- 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 : public Cjet{
 public:
  Cjet_area();
 
  Cjet_area(Cjet &j);
 
  ~Cjet_area();
 
  // area information
  double passive_area;   
  double active_area;    
};
 
class Carea : public Csiscone{
 public:
  Carea();
 
  ~Carea();
 
  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);
 
  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);
 
  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;        
  double grid_eta_max;  
  double grid_shift;    
 
  double pt_soft;       
  double pt_shift;      
  double pt_soft_min;   
 
  std::vector<Cjet_area> jet_areas;
};
 
}
#endif