Rivet::CentralityBinner< T, MDist > Class Template Reference

#include <CentralityBinner.hh>

Inheritance diagram for Rivet::CentralityBinner< T, MDist >:

Public Member Functions
	CentralityBinner (int maxbins=200, double wlim=0.02)
void	setProjection (const CentralityEstimator &p, string pname)
virtual std::string	name () const
	Return the class name.
void	add (T t, double cmin, double cmax, double cestmin=-1.0, double cestmax=-1.0)
T	select (const Event &event, double weight=1.0)
T	select (double cest, double weight=1.0)
void	finalize ()
void	normalizePerEvent ()
map< double, double >	edges () const
const T &	current () const
	Return the current AnalysisObject from the latest call to select().
double	estimator () const
vector< T >	allObjects ()
void	debug ()
	Print out the _flexiBins to cerr.
void	fulldebug ()
Public Member Functions inherited from Rivet::ProjectionApplier
	ProjectionApplier ()
	Constructor.
void	markAsOwned () const
	Mark this object as owned by a proj-handler.
std::set< ConstProjectionPtr >	getProjections () const
	Get the contained projections, including recursion.
bool	hasProjection (const std::string &name) const
	Does this applier have a projection registered under the name name?
template<typename PROJ >
const PROJ &	getProjection (const std::string &name) const
template<typename PROJ >
const PROJ &	get (const std::string &name) const
const Projection &	getProjection (const std::string &name) const
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const Projection &proj) const
	Apply the supplied projection on event evt. More...
template<typename PROJ >
const PROJ &	apply (const Event &evt, const Projection &proj) const
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const PROJ &proj) const
	Apply the supplied projection on event evt. More...
template<typename PROJ >
const PROJ &	apply (const Event &evt, const PROJ &proj) const
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const std::string &name) const
template<typename PROJ >
const PROJ &	apply (const Event &evt, const std::string &name) const
template<typename PROJ >
const PROJ &	apply (const std::string &name, const Event &evt) const

Protected Types
typedef set< FlexiBin >	FlexiBinSet
	Convenient typedefs.

Additional Inherited Members
Protected Member Functions inherited from Rivet::ProjectionApplier
Log &	getLog () const
ProjectionHandler &	getProjHandler () const
	Get a reference to the ProjectionHandler for this thread.
template<typename PROJ >
const PROJ &	declareProjection (const PROJ &proj, const std::string &name)
	Register a contained projection. More...
template<typename PROJ >
const PROJ &	declare (const PROJ &proj, const std::string &name)
	Register a contained projection (user-facing version) More...
template<typename PROJ >
const PROJ &	declare (const std::string &name, const PROJ &proj)
	Register a contained projection (user-facing, arg-reordered version) More...
template<typename PROJ >
const PROJ &	addProjection (const PROJ &proj, const std::string &name)
	Register a contained projection (user-facing version) More...

Detailed Description

template<typename T = Histo1DPtr, typename MDist = MergeDistance>
class Rivet::CentralityBinner< T, MDist >

CentralityBinner contains a series of AnalysisObject of the same quantity each in a different percentiles of another quantity. For example, a CentralityBinner may e.g. contain histograms of the cross section differential in in different centrality regions for heavy ion collisions based on forward energy flow.

Constructor & Destructor Documentation

CentralityBinner()

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

Rivet::CentralityBinner< T, MDist >::CentralityBinner ( int  maxbins = 200, double  wlim = 0.02  )

inline

Create a new empty CentralityBinner. maxbins is the maximum number of bins used by the binner. Default is 1000, which is typically enough. wlim is the mximum allowed error allowed for the centrality limits before a warning is emitted.

Member Function Documentation

add()

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

void Rivet::CentralityBinner< T, MDist >::add ( T  t, double  cmin, double  cmax, double  cestmin = -1.0, double  cestmax = -1.0  )

inline

Add an AnalysisObject in the region between cmin and cmax to this set of CentralityBinners. The range represent percentiles and must be between 0 and 100. No overlaping bins are allowed. Note that (cmin=0, cmax=5), means the five percent MOST central events although the internal notation is reversed for convenience. Optionally supply corresponding limits cestmin and cestmax of the centrality extimator.

References Rivet::max(), and Rivet::min().

edges()

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

map<double,double> Rivet::CentralityBinner< T, MDist >::edges ( ) const

inline

Return a map bin edges of the centrality extimator indexed by the corresponing percentile.

estimator()

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

double Rivet::CentralityBinner< T, MDist >::estimator ( ) const

inline

Return the value of the centrality estimator set in the latest call to select().

References Rivet::CentralityBinTraits< T >::add(), Rivet::inRange(), Rivet::max(), Rivet::min(), and Rivet::CentralityBinTraits< T >::normalize().

finalize()

template<typename T , typename MDist >

void Rivet::CentralityBinner< T, MDist >::finalize

(

)

At the end of the run, calculate the percentiles and fill the AnalysisObjectss provided with the add() function. This is typically called from the finalize method in an Analysis, but can also be called earlier in which case the the select functions can be continued to run as before with the edges between the centrality regions now fixed.

References Rivet::CentralityBinTraits< T >::add(), Rivet::CentralityBinTraits< T >::clone(), Rivet::max(), Rivet::min(), and Rivet::CentralityBinTraits< T >::scale().

normalizePerEvent()

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

void Rivet::CentralityBinner< T, MDist >::normalizePerEvent ( )

inline

Normalize each AnalysisObjects to the sum of event weights in the corresponding centrality bin.

select() [1/2]

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

T Rivet::CentralityBinner< T, MDist >::select ( const Event &  event, double  weight = 1.0  )

inline

Return one of the AnalysisObjects in the CentralityBinner for the given event. This version requires that a CentralityEstimator object has been assigned that can compute the value of the centrality estimator from the event. Optionally the weight of the event is given. This should be the weight that will be used to fill the AnalysisObject. If the centrality estimate is less than zero, the _devnull object will be returned.

References Rivet::CentralityEstimator::estimate().

select() [2/2]

template<typename T , typename MDist >

T Rivet::CentralityBinner< T, MDist >::select	(	double	cest,
		double	weight = 1.0
	)

Return one of the AnalysisObjecsts in the Setup the CentralityBinner depending on the value of the centrality estimator, cest. Optionally the weight of the event is given. This should be the weight that will be used to fill the AnalysisObject. If the centrality estimate is less than zero, the _devnull object will be returned.

References Rivet::CentralityBinTraits< T >::clone().

setProjection()

template<typename T = Histo1DPtr, typename MDist = MergeDistance>

void Rivet::CentralityBinner< T, MDist >::setProjection ( const CentralityEstimator &  p, string  pname  )

inline

Set the centrality projection to be used. Note that this projection must have already been declared to Rivet.

References Rivet::ProjectionApplier::declare().

---

The documentation for this class was generated from the following file:

• include/Rivet/Tools/CentralityBinner.hh