00001 //STARTHEADER 00002 // $Id: RangeDefinition.cc 760 2007-06-25 16:07:05Z salam $ 00003 // 00004 // Copyright (c) 2005-2007, Matteo Cacciari and Gavin Salam 00005 // 00006 //---------------------------------------------------------------------- 00007 // This file is part of FastJet. 00008 // 00009 // FastJet is free software; you can redistribute it and/or modify 00010 // it under the terms of the GNU General Public License as published by 00011 // the Free Software Foundation; either version 2 of the License, or 00012 // (at your option) any later version. 00013 // 00014 // The algorithms that underlie FastJet have required considerable 00015 // development and are described in hep-ph/0512210. If you use 00016 // FastJet as part of work towards a scientific publication, please 00017 // include a citation to the FastJet paper. 00018 // 00019 // FastJet is distributed in the hope that it will be useful, 00020 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00021 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00022 // GNU General Public License for more details. 00023 // 00024 // You should have received a copy of the GNU General Public License 00025 // along with FastJet; if not, write to the Free Software 00026 // Foundation, Inc.: 00027 // 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00028 //---------------------------------------------------------------------- 00029 //ENDHEADER 00030 00031 #include "fastjet/RangeDefinition.hh" 00032 00033 FASTJET_BEGIN_NAMESPACE // defined in fastjet/internal/base.hh 00034 00035 using namespace std; 00036 00039 void RangeDefinition::_numerical_total_area(double rapmax, int npoints) { 00040 00041 int count = 0; 00042 double deltaphi = twopi/double(npoints); 00043 double deltarap = 2.0*rapmax/double(npoints); 00044 double phi = 0.0; 00045 for(int i = 0; i < npoints; i++) { 00046 double rap = -rapmax; 00047 for (int j = 0; j < npoints; j++) { 00048 if ( is_in_range(rap,phi) ) { count++; } 00049 rap += deltarap; 00050 } 00051 phi += deltaphi; 00052 } 00053 00054 _total_area = double(count)/double(npoints*npoints)*2.0*twopi*rapmax; 00055 } 00056 00057 FASTJET_END_NAMESPACE