ClusterSequenceAreaBase.hh

Go to the documentation of this file.
00001 //STARTHEADER
00002 // $Id: ClusterSequenceAreaBase.hh 958 2007-11-28 17:43:48Z cacciari $
00003 //
00004 // Copyright (c) 2005-2006, 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 #ifndef __FASTJET_CLUSTERSEQUENCEAREABASE_HH__
00032 #define __FASTJET_CLUSTERSEQUENCEAREABASE_HH__
00033 
00034 #include "fastjet/ClusterSequence.hh"
00035 #include "fastjet/internal/LimitedWarning.hh"
00036 #include "fastjet/RangeDefinition.hh"
00037 
00038 FASTJET_BEGIN_NAMESPACE
00039 
00045 class ClusterSequenceAreaBase : public ClusterSequence {
00046 public:
00047   
00050   template<class L> ClusterSequenceAreaBase
00051          (const std::vector<L> & pseudojets, 
00052           const JetDefinition & jet_def,
00053           const bool & writeout_combinations = false) :
00054            ClusterSequence(pseudojets, jet_def, writeout_combinations) {}
00055 
00056 
00058   ClusterSequenceAreaBase() {}
00059 
00060 
00062   virtual ~ClusterSequenceAreaBase() {}
00063 
00064 
00067   virtual double area       (const PseudoJet & jet) const {return 0.0;}
00068 
00071   virtual double area_error (const PseudoJet & jet) const {return 0.0;}
00072 
00084   virtual PseudoJet area_4vector(const PseudoJet & jet) const {
00085     return PseudoJet(0.0,0.0,0.0,0.0);}
00086 
00092   virtual bool is_pure_ghost(const PseudoJet & jet) const {
00093     return false;
00094   }
00095 
00097   virtual double empty_area(const RangeDefinition & range) const;
00098 
00104   virtual double n_empty_jets(const RangeDefinition & range) const {
00105     double R = jet_def().R();
00106     return empty_area(range)/(0.55*pi*R*R);
00107   }
00108 
00111   double median_pt_per_unit_area(const RangeDefinition & range) const;
00112 
00115   double median_pt_per_unit_area_4vector(const RangeDefinition & range) const;
00116   
00121   double median_pt_per_unit_something(
00122                     const RangeDefinition & range, bool use_area_4vector) const;
00123 
00140   virtual void get_median_rho_and_sigma(const RangeDefinition & range, 
00141                                         bool use_area_4vector,
00142                                         double & median, double & sigma,
00143                                         double & mean_area) const;
00144 
00147   void get_median_rho_and_sigma(const RangeDefinition & range, 
00148                                 bool use_area_4vector,
00149                                 double & median, double & sigma) const {
00150     double mean_area;
00151     get_median_rho_and_sigma(range,  use_area_4vector,
00152                              median,  sigma, mean_area);
00153   }
00154   
00155 
00159   virtual void parabolic_pt_per_unit_area(double & a, double & b, 
00160                                           const RangeDefinition & range, 
00161                                           double exclude_above=-1.0, 
00162                                           bool use_area_4vector=false) const;
00163 
00164 
00169   std::vector<PseudoJet> subtracted_jets(const double rho,
00170                                          const double ptmin=0.0) const;
00171 
00176   std::vector<PseudoJet> subtracted_jets(const RangeDefinition & range, 
00177                                          const double ptmin=0.0) const;
00178 
00180   PseudoJet subtracted_jet(const PseudoJet & jet,
00181                            const double rho) const;
00182 
00187   PseudoJet subtracted_jet(const PseudoJet & jet,
00188                            const RangeDefinition & range) const;
00189 
00191   double subtracted_pt(const PseudoJet & jet,
00192                        const double rho,
00193                        bool use_area_4vector=false) const;
00194 
00198   double subtracted_pt(const PseudoJet & jet,
00199                        const RangeDefinition & range,
00200                        bool use_area_4vector=false) const;
00201 
00202 
00203 private:
00205   static LimitedWarning _warnings;
00206 
00208   void _check_jet_alg_good_for_median() const;
00209   
00210 };
00211 
00212 
00213 
00214 FASTJET_END_NAMESPACE
00215 
00216 #endif // __FASTJET_CLUSTERSEQUENCEAREABASE_HH__

Generated on Tue Dec 18 17:05:02 2007 for fastjet by  doxygen 1.5.2