fastjet: fastjet::ClusterSequenceVoronoiArea Class Reference

List of all members.


Public Member Functions
template<class L>
	ClusterSequenceVoronoiArea (const std::vector< L > &pseudojets, const JetDefinition &jet_def, const VoronoiAreaSpec &spec=VoronoiAreaSpec(), const bool &writeout_combinations=false)
	template ctor
	~ClusterSequenceVoronoiArea ()
	default dtor
virtual double	area (const PseudoJet &jet) const
	return the area associated with the given jet
virtual PseudoJet	area_4vector (const PseudoJet &jet) const
	return a 4-vector area associated with the given jet -- stricly this is not the exact 4-vector area, but rather an approximation made of sums of centres of all Voronoi cells in jet, each contributing with a normalisation equal to the area of the cell
virtual double	area_error (const PseudoJet &jet) const
	return the error of the area associated with the given jet (0 by definition for a voronoi area)
Private Member Functions
void	_initializeVA ()
	initialisation of the Voronoi Area
Private Attributes
std::vector< double >	_voronoi_area
	vector containing the result
std::vector< PseudoJet >	_voronoi_area_4vector
	vector containing approx 4-vect areas
VoronoiAreaCalc *	_pa_calc
	area calculator
double	_effective_Rfact
	effective radius
Classes
class	VoronoiAreaCalc
	class for carrying out a voronoi area calculation on a set of initial vectors More...

Detailed Description

Constructor & Destructor Documentation

template<class L>

fastjet::ClusterSequenceVoronoiArea::ClusterSequenceVoronoiArea	(	const std::vector< L > &	pseudojets,
		const JetDefinition &	jet_def,
		const VoronoiAreaSpec &	spec = `VoronoiAreaSpec()`,
		const bool &	writeout_combinations = `false`
	)			`[inline]`

template ctor

template constructor need to be specified in the header!

Parameters:

	pseudojet	list of jets (template type)
	jet_def	jet definition
	effective_Rfact	effective radius
	writeout_combinations	??????

Definition at line 99 of file ClusterSequenceVoronoiArea.hh.

00102                                      :
00103   _effective_Rfact(spec.effective_Rfact()) {
00104 
00105   // transfer the initial jets (type L) into our own array
00106   _transfer_input_jets(pseudojets);
00107 
00108   // run the clustering
00109   _initialise_and_run(jet_def,writeout_combinations);
00110 
00111   // the jet clustering's already been done, now worry about areas...
00112   _initializeVA();
00113 }

fastjet::ClusterSequenceVoronoiArea::~ClusterSequenceVoronoiArea ( )

default dtor

Definition at line 309 of file ClusterSequenceVoronoiArea.cc.

References _pa_calc.

00309                                                         {
00310   delete _pa_calc;
00311 }

Member Function Documentation

virtual double fastjet::ClusterSequenceVoronoiArea::area ( const PseudoJet & jet ) const [inline, virtual]

return the area associated with the given jet

Reimplemented from fastjet::ClusterSequenceAreaBase.

Definition at line 63 of file ClusterSequenceVoronoiArea.hh.

References fastjet::PseudoJet::cluster_hist_index().

Referenced by fastjet::ClusterSequencePassiveArea::_initialise_and_run_PA(), and _initializeVA().

00063                                                           {
00064     return _voronoi_area[jet.cluster_hist_index()];};

virtual PseudoJet fastjet::ClusterSequenceVoronoiArea::area_4vector ( const PseudoJet & jet ) const [inline, virtual]

return a 4-vector area associated with the given jet -- stricly this is not the exact 4-vector area, but rather an approximation made of sums of centres of all Voronoi cells in jet, each contributing with a normalisation equal to the area of the cell

Reimplemented from fastjet::ClusterSequenceAreaBase.

Definition at line 70 of file ClusterSequenceVoronoiArea.hh.

References fastjet::PseudoJet::cluster_hist_index().

Referenced by fastjet::ClusterSequencePassiveArea::_initialise_and_run_PA().

00070                                                                      {
00071     return _voronoi_area_4vector[jet.cluster_hist_index()];};

virtual double fastjet::ClusterSequenceVoronoiArea::area_error ( const PseudoJet & jet ) const [inline, virtual]

return the error of the area associated with the given jet (0 by definition for a voronoi area)

Reimplemented from fastjet::ClusterSequenceAreaBase.

Definition at line 75 of file ClusterSequenceVoronoiArea.hh.

00075                                                                 {
00076     return 0.0;};

void fastjet::ClusterSequenceVoronoiArea::_initializeVA ( ) [private]

initialisation of the Voronoi Area

Definition at line 267 of file ClusterSequenceVoronoiArea.cc.

References _effective_Rfact, fastjet::ClusterSequence::_history, fastjet::ClusterSequence::_jet_def, fastjet::ClusterSequence::_jets, _pa_calc, _voronoi_area, _voronoi_area_4vector, area(), fastjet::ClusterSequenceVoronoiArea::VoronoiAreaCalc::area(), fastjet::ClusterSequence::n_particles(), and fastjet::JetDefinition::R().

00267                                                 {
00268   // run the VAC on our original particles
00269   _pa_calc = new VAC(_jets.begin(), 
00270                      _jets.begin()+n_particles(),
00271                      _effective_Rfact*_jet_def.R());
00272 
00273   // transfer the areas to our local structure
00274   //  -- first the initial ones
00275   _voronoi_area.reserve(2*n_particles());
00276   for (unsigned int i=0; i<n_particles(); i++) {
00277     _voronoi_area.push_back(_pa_calc->area(i));
00278     // make a stab at a 4-vector area
00279     if (_jets[i].perp2() > 0) {
00280       _voronoi_area_4vector.push_back((_pa_calc->area(i)/_jets[i].perp())
00281                                       * _jets[i]);
00282     } else {
00283       // not sure what to do here -- just put zero (it won't be meaningful
00284       // anyway)
00285       _voronoi_area_4vector.push_back(PseudoJet(0.0,0.0,0.0,0.0));
00286     }
00287   }
00288            
00289   //  -- then the combined areas that arise from the clustering
00290   for (unsigned int i = n_particles(); i < _history.size(); i++) {
00291     double area;
00292     PseudoJet area_4vect;
00293     if (_history[i].parent2 >= 0) {
00294       area = _voronoi_area[_history[i].parent1] + 
00295              _voronoi_area[_history[i].parent2];
00296       area_4vect = _voronoi_area_4vector[_history[i].parent1] + 
00297                    _voronoi_area_4vector[_history[i].parent2];
00298     } else {
00299       area = _voronoi_area[_history[i].parent1];
00300       area_4vect = _voronoi_area_4vector[_history[i].parent1];
00301     }
00302     _voronoi_area.push_back(area);
00303     _voronoi_area_4vector.push_back(area_4vect);
00304   }
00305 
00306 }

Member Data Documentation

std::vector<double> fastjet::ClusterSequenceVoronoiArea::_voronoi_area [private]

vector containing the result

Definition at line 87 of file ClusterSequenceVoronoiArea.hh.

Referenced by _initializeVA().

fastjet::ClusterSequenceVoronoiArea Class Reference

Public Member Functions

Private Member Functions

Private Attributes

Classes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation