Main Page | Namespace List | Class Hierarchy | Class List | Directories | File List | Namespace Members | Class Members | File Members
fastjet::ActiveAreaSpec

fastjet::ActiveAreaSpec Class Reference

Class that defines the parameters that go into the measurement of active jet areas. More...

#include <ActiveAreaSpec.hh>

Collaboration diagram for fastjet::ActiveAreaSpec:

Collaboration graph
[legend]
List of all members.

Public Member Functions

 ActiveAreaSpec ()
 default constructor
 ActiveAreaSpec (double ghost_maxrap, int repeat=1, double ghost_area=0.01, double grid_scatter=1e-4, double kt_scatter=0.1, double mean_ghost_kt=1e-100)
 explicit constructor
void _initialize ()
 sets the detailed parameters for the ghosts (which may not be quite the same as those requested -- this is in order for things to fit in nicely into 2pi etc.
double ghost_etamax () const
double ghost_maxrap () const
double ghost_area () const
double grid_scatter () const
double kt_scatter () const
double mean_ghost_kt () const
int repeat () const
double actual_ghost_area () const
int n_ghosts () const
void set_ghost_area (double val)
void set_ghost_etamax (double val)
void set_ghost_maxrap (double val)
void set_grid_scatter (double val)
void set_kt_scatter (double val)
void set_mean_ghost_kt (double val)
void set_repeat (int val)
std::string description ()
 for a summary
void add_ghosts (std::vector< PseudoJet > &) const
 adds the ghost 4-momenta to the vector of PseudoJet's

Private Member Functions

double _our_rand () const

Private Attributes

double _ghost_maxrap
int _repeat
double _ghost_area
double _grid_scatter
double _kt_scatter
double _mean_ghost_kt
double _actual_ghost_area
double _dphi
double _drap
int _n_ghosts
int _nphi
int _nrap

Static Private Attributes

static BasicRandom< double > _random_generator

Detailed Description

Class that defines the parameters that go into the measurement of active jet areas.

Definition at line 48 of file ActiveAreaSpec.hh.

Constructor & Destructor Documentation

fastjet::ActiveAreaSpec::ActiveAreaSpec (   )  [inline]
 

default constructor

Definition at line 51 of file ActiveAreaSpec.hh.

00051                   : _ghost_maxrap(6.0), _repeat(1), _ghost_area(0.01), 
00052                     _grid_scatter(1e-4), _kt_scatter(0.1), 
00053                     _mean_ghost_kt(1e-100),
00054                     _actual_ghost_area(-1.0) {_initialize();};

fastjet::ActiveAreaSpec::ActiveAreaSpec (  double  ghost_maxrap,
int  repeat = 1,
double  ghost_area = 0.01,
double  grid_scatter = 1e-4,
double  kt_scatter = 0.1,
double  mean_ghost_kt = 1e-100
)  [inline, explicit]
 

explicit constructor

Definition at line 57 of file ActiveAreaSpec.hh.

00062                                                         : 
00063     _ghost_maxrap(ghost_maxrap), 
00064     _repeat(repeat), 
00065     _ghost_area(ghost_area), 
00066     _grid_scatter(grid_scatter),  
00067     _kt_scatter(kt_scatter), 
00068     _mean_ghost_kt(mean_ghost_kt),
00069     _actual_ghost_area(-1.0) {_initialize();};

Member Function Documentation

void fastjet::ActiveAreaSpec::_initialize (   ) 
 

sets the detailed parameters for the ghosts (which may not be quite the same as those requested -- this is in order for things to fit in nicely into 2pi etc.

..

Definition at line 45 of file ActiveAreaSpec.cc.

References _actual_ghost_area, _dphi, _drap, _ghost_area, _ghost_maxrap, _n_ghosts, _nphi, _nrap, and fastjet::twopi.

00045                                  {
00046   // add on area-measuring dummy particles
00047   _drap = sqrt(_ghost_area);
00048   _dphi = _drap;
00049   _nphi = int(ceil(twopi/_dphi)); _dphi = twopi/_nphi;
00050   _nrap = int(ceil(_ghost_maxrap/_drap)); _drap = _ghost_maxrap / _nrap;
00051   _actual_ghost_area = _dphi * _drap;
00052   _n_ghosts   = (2*_nrap+1)*_nphi;
00053 
00054   //_random_generator.info(cerr);
00055 }

double fastjet::ActiveAreaSpec::_our_rand (   )  const [inline, private]
 

Definition at line 118 of file ActiveAreaSpec.hh.

Referenced by add_ghosts().

00118 {return _random_generator();};

double fastjet::ActiveAreaSpec::actual_ghost_area (   )  const [inline]
 

Definition at line 85 of file ActiveAreaSpec.hh.

Referenced by fastjet::ClusterSequenceActiveAreaExplicitGhosts::_add_ghosts(), and description().

00085 {return _actual_ghost_area;};

void fastjet::ActiveAreaSpec::add_ghosts (  std::vector< PseudoJet > &   )  const
 

adds the ghost 4-momenta to the vector of PseudoJet's

Definition at line 59 of file ActiveAreaSpec.cc.

References _dphi, _drap, _grid_scatter, _kt_scatter, _mean_ghost_kt, _nphi, _nrap, and _our_rand().

Referenced by fastjet::ClusterSequenceActiveAreaExplicitGhosts::_add_ghosts().

00059                                                                {
00060   // add momenta for ghosts
00061   for (int irap = -_nrap; irap <= _nrap; irap++) {
00062     for (int iphi = 0; iphi < _nphi; iphi++) {
00063       // // include random offsets for all quantities
00064       // double phi = (iphi+0.5) * _dphi + _dphi*rand()*_grid_scatter/RAND_MAX;
00065       // double rap = irap * _drap + _drap*rand()*_grid_scatter/RAND_MAX;
00066       // //double phi = (iphi+0.5) * _dphi* + rand()*_grid_scatter/RAND_MAX;
00067       // //double rap = irap * _drap + rand()*_grid_scatter/RAND_MAX;
00068       // double kt = _mean_ghost_kt*(1+rand()*_kt_scatter/RAND_MAX);
00069      
00070       // include random offsets for all quantities
00071       double phi = (iphi+0.5) * _dphi + _dphi*_our_rand()*_grid_scatter;
00072       double rap = irap * _drap + _drap*_our_rand()*_grid_scatter;
00073       //double phi = (iphi+0.5) * _dphi* + _our_rand()*_grid_scatter;
00074       //double rap = irap * _drap + _our_rand()*_grid_scatter;
00075       double kt = _mean_ghost_kt*(1+_our_rand()*_kt_scatter);
00076 
00077 
00078 
00079       double pminus = kt*exp(-rap);
00080       double pplus  = kt*exp(+rap);
00081       double px = kt*sin(phi);
00082       double py = kt*cos(phi);
00083       //cout << kt<<" "<<rap<<" "<<phi<<"\n";
00084       //if (phi>=twopi || phi < 0.0) cout << "Hey: "<< phi-twopi<<"\n";
00085       PseudoJet mom(px,py,0.5*(pplus-pminus),0.5*(pplus+pminus));
00086       //mom.set_user_index(1);  // for ghost particles (user index now lost...)
00087       event.push_back(mom);
00088       //_is_pure_ghost.push_back(true);
00089     }
00090   }
00091 }

string fastjet::ActiveAreaSpec::description (   ) 
 

for a summary

Definition at line 93 of file ActiveAreaSpec.cc.

References actual_ghost_area(), ghost_area(), ghost_maxrap(), grid_scatter(), kt_scatter(), mean_ghost_kt(), and repeat().

00093                                    {
00094 
00095   ostringstream ostr;
00096   ostr << "Active area specification with ghosts of area " << actual_ghost_area() 
00097        << " (had requested " << ghost_area() << ")"
00098        << ", placed up to y = " << ghost_maxrap() 
00099        << ", scattered rel. to perfect grid by " << grid_scatter() 
00100        << ", mean_ghost_kt = " << mean_ghost_kt()
00101        << ", rel kt_scatter =  " << kt_scatter()
00102        << ", n repetitions of ghost distributions =  " << repeat();
00103   return ostr.str();
00104 }

double fastjet::ActiveAreaSpec::ghost_area (   )  const [inline]
 

Definition at line 78 of file ActiveAreaSpec.hh.

Referenced by description().

00078 {return _ghost_area   ;};

double fastjet::ActiveAreaSpec::ghost_etamax (   )  const [inline]
 

Definition at line 76 of file ActiveAreaSpec.hh.

00076 {return _ghost_maxrap;};

double fastjet::ActiveAreaSpec::ghost_maxrap (   )  const [inline]
 

Definition at line 77 of file ActiveAreaSpec.hh.

Referenced by fastjet::ClusterSequenceActiveArea::_initialise_and_run_AA(), and description().

00077 {return _ghost_maxrap;};

double fastjet::ActiveAreaSpec::grid_scatter (   )  const [inline]
 

Definition at line 79 of file ActiveAreaSpec.hh.

Referenced by description().

00079 {return _grid_scatter;};

double fastjet::ActiveAreaSpec::kt_scatter (   )  const [inline]
 

Definition at line 80 of file ActiveAreaSpec.hh.

Referenced by description().

00080 {return _kt_scatter  ;};

double fastjet::ActiveAreaSpec::mean_ghost_kt (   )  const [inline]
 

Definition at line 81 of file ActiveAreaSpec.hh.

Referenced by description().

00081 {return _mean_ghost_kt  ;};

int fastjet::ActiveAreaSpec::n_ghosts (   )  const [inline]
 

Definition at line 86 of file ActiveAreaSpec.hh.

Referenced by fastjet::ClusterSequenceActiveAreaExplicitGhosts::_add_ghosts().

00086 {return _n_ghosts;};

int fastjet::ActiveAreaSpec::repeat (   )  const [inline]
 

Definition at line 82 of file ActiveAreaSpec.hh.

Referenced by fastjet::ClusterSequenceActiveArea::_initialise_and_run_AA(), and description().

00082 {return _repeat      ;};

void fastjet::ActiveAreaSpec::set_ghost_area (  double  val  )  [inline]
 

Definition at line 89 of file ActiveAreaSpec.hh.

00089 {_ghost_area    = val; _initialize();};

void fastjet::ActiveAreaSpec::set_ghost_etamax (  double  val  )  [inline]
 

Definition at line 90 of file ActiveAreaSpec.hh.

00090 {_ghost_maxrap = val; _initialize();};

void fastjet::ActiveAreaSpec::set_ghost_maxrap (  double  val  )  [inline]
 

Definition at line 91 of file ActiveAreaSpec.hh.

00091 {_ghost_maxrap = val; _initialize();};

void fastjet::ActiveAreaSpec::set_grid_scatter (  double  val  )  [inline]
 

Definition at line 92 of file ActiveAreaSpec.hh.

00092 {_grid_scatter   = val; };

void fastjet::ActiveAreaSpec::set_kt_scatter (  double  val  )  [inline]
 

Definition at line 93 of file ActiveAreaSpec.hh.

00093 {_kt_scatter     = val; };

void fastjet::ActiveAreaSpec::set_mean_ghost_kt (  double  val  )  [inline]
 

Definition at line 94 of file ActiveAreaSpec.hh.

00094 {_mean_ghost_kt  = val; };

void fastjet::ActiveAreaSpec::set_repeat (  int  val  )  [inline]
 

Definition at line 95 of file ActiveAreaSpec.hh.

00095 {_repeat         = val; };

Member Data Documentation

double fastjet::ActiveAreaSpec::_actual_ghost_area [private]
 

Definition at line 114 of file ActiveAreaSpec.hh.

Referenced by _initialize().

double fastjet::ActiveAreaSpec::_dphi [private]
 

Definition at line 114 of file ActiveAreaSpec.hh.

Referenced by _initialize(), and add_ghosts().

double fastjet::ActiveAreaSpec::_drap [private]
 

Definition at line 114 of file ActiveAreaSpec.hh.

Referenced by _initialize(), and add_ghosts().

double fastjet::ActiveAreaSpec::_ghost_area [private]
 

Definition at line 108 of file ActiveAreaSpec.hh.

Referenced by _initialize().

double fastjet::ActiveAreaSpec::_ghost_maxrap [private]
 

Definition at line 106 of file ActiveAreaSpec.hh.

Referenced by _initialize().

double fastjet::ActiveAreaSpec::_grid_scatter [private]
 

Definition at line 109 of file ActiveAreaSpec.hh.

Referenced by add_ghosts().

double fastjet::ActiveAreaSpec::_kt_scatter [private]
 

Definition at line 110 of file ActiveAreaSpec.hh.

Referenced by add_ghosts().

double fastjet::ActiveAreaSpec::_mean_ghost_kt [private]
 

Definition at line 111 of file ActiveAreaSpec.hh.

Referenced by add_ghosts().

int fastjet::ActiveAreaSpec::_n_ghosts [private]
 

Definition at line 115 of file ActiveAreaSpec.hh.

Referenced by _initialize().

int fastjet::ActiveAreaSpec::_nphi [private]
 

Definition at line 115 of file ActiveAreaSpec.hh.

Referenced by _initialize(), and add_ghosts().

int fastjet::ActiveAreaSpec::_nrap [private]
 

Definition at line 115 of file ActiveAreaSpec.hh.

Referenced by _initialize(), and add_ghosts().

BasicRandom< double > fastjet::ActiveAreaSpec::_random_generator [static, private]
 

Definition at line 39 of file ActiveAreaSpec.cc.

int fastjet::ActiveAreaSpec::_repeat [private]
 

Definition at line 107 of file ActiveAreaSpec.hh.

The documentation for this class was generated from the following files:
Generated on Mon Apr 2 20:58:16 2007 for fastjet by  doxygen 1.4.2