rivet is hosted by Hepforge, IPPP Durham
Rivet 3.1.6

Projection for calculating correlators for flow measurements. More...

#include <Correlators.hh>

Inheritance diagram for Rivet::Correlators:
Rivet::Projection Rivet::ProjectionApplier

Public Member Functions

 Correlators (const ParticleFinder &fsp, int nMaxIn=2, int pMaxIn=0, vector< double > pTbinEdgesIn={})
 
 Correlators (const ParticleFinder &fsp, int nMaxIn, int pMaxIn, const YODA::Scatter2D hIn)
 
const pair< double, double > intCorrelator (vector< int > n) const
 Integrated correlator of n harmonic, with the number of powers being the size of n. E.G. n should be:
<<2>>_2 => n = {2, -2}
<<4>>_2 => n = {2, 2, -2, -2}
<<2>>_4 => n = {4, -4}
<<4>>_4 => n = {4, 4, -4, 4}, and so on.
 
const vector< pair< double, double > > pTBinnedCorrelators (vector< int > n, bool overflow=false) const
 pT differential correlator of n harmonic, for number of powers n.size() More...
 
const pair< double, double > intCorrelatorGap (const Correlators &other, vector< int > n1, vector< int > n2) const
 Integrated correlator of n1 harmonic, for number of powers n1.size() More...
 
const vector< pair< double, double > > pTBinnedCorrelatorsGap (const Correlators &other, vector< int > n1, vector< int > n2, bool overflow=false) const
 pT differential correlators of n1 harmonic, for number n1.size() More...
 
 DEFAULT_RIVET_PROJ_CLONE (Correlators)
 
Standard constructors and destructors.
virtual unique_ptr< Projectionclone () const =0
 Clone on the heap. More...
 
virtual std::string name () const
 Get the name of the projection. More...
 
bool valid () const
 Get the state of the projetion.
 
bool failed () const
 Get the state of the projetion.
 
Projection operation and comparison
bool before (const Projection &p) const
 
Projection "getting" functions
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
 
const ProjectiongetProjection (const std::string &name) const
 
template<typename PROJ >
const PROJ & get (const std::string &name) const
 

Static Public Member Functions

static vector< int > hVec (int n, int m)
 Construct a harmonic vectors from n harmonics and m number of particles. More...
 
static pair< int, int > getMaxValues (vector< vector< int > > &hList)
 Return the maximal values for n, p to be used in the constructor of Correlators(xxx, nMax, pMax, xxxx)
 

Protected Member Functions

void project (const Event &e)
 
CmpState compare (const Projection &p) const
 
void fillCorrelators (const Particle &p, const double &weight)
 
const complex< double > getQ (int n, int p) const
 
const complex< double > getP (int n, int p, double pT=0.) const
 
Projection registration functions
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...
 

Beam configuration

Todo:
Does it really make sense to restrict Projections to particular beam configs? Do we use this in practice?
virtual const std::set< PdgIdPair > beamPairs () const
 
ProjectionaddPdgIdPair (PdgId beam1, PdgId beam2)
 
LoggetLog () const
 Get a Log object based on the getName() property of the calling projection object.
 
void setName (const std::string &name)
 Used by derived classes to set their name.
 
void fail ()
 Set the projection in an unvalid state.
 
Cmp< ProjectionmkNamedPCmp (const Projection &otherparent, const std::string &pname) const
 
Cmp< ProjectionmkPCmp (const Projection &otherparent, const std::string &pname) const
 

Projection applying functions

template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > applyProjection (const Event &evt, const Projection &proj) const
 
template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > applyProjection (const Event &evt, const PROJ &proj) const
 
template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > applyProjection (const Event &evt, const std::string &name) const
 
template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > apply (const Event &evt, const Projection &proj) const
 Apply the supplied projection on event evt (user-facing alias).
 
template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > apply (const Event &evt, const PROJ &proj) const
 Apply the supplied projection on event evt (user-facing alias).
 
template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > apply (const Event &evt, const std::string &name) const
 Apply the supplied projection on event evt (user-facing alias).
 
template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > apply (const std::string &name, const Event &evt) const
 Apply the supplied projection on event evt (convenience arg-reordering alias).
 
void markAsOwned () const
 Mark this object as owned by a proj-handler.
 
ProjectionHandlergetProjHandler () const
 Get a reference to the ProjectionHandler for this thread.
 

Detailed Description

Projection for calculating correlators for flow measurements.

A projection which calculates Q-vectors and P-vectors, and projects them out as correlators. Implementation follows the description of the ''Generic Framework'':
Phys. Rev. C 83 (2011) 044913, arXiv: 1010.0233
Phys. Rev. C 89 (2014) 064904, arXiv: 1312.3572

Constructor & Destructor Documentation

◆ Correlators()

Rivet::Correlators::Correlators ( const ParticleFinder fsp,
int  nMaxIn = 2,
int  pMaxIn = 0,
vector< double >  pTbinEdgesIn = {} 
)

Constructor

Parameters
fspThe FinalState projection that the correlators should be constructed from
nMaxInThe maximal sum of harmonics, e.g. for
c_2{2} = {2,-2} = 2 + 2 = 4
c_2{4} = {2,2,-2,-2} = 2 + 2 + 2 + 2 = 8
c_4{2} = {4,-4} = 4 + 4 = 8
c_4{4} = {4,4,-4,-4} = 4 + 4 + 4 + 4 = 16.
pMaxInThe maximal number of particles you want to correlate
pTbinEdgesInThe (lower) edges of pT bins, the last one the upper edge of the final bin.

Member Function Documentation

◆ addPdgIdPair()

Projection & Rivet::Projection::addPdgIdPair ( PdgId  beam1,
PdgId  beam2 
)
inlineinherited

Add a colliding beam pair.

Todo:
This deserves a better name!

◆ applyProjection() [1/3]

template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > Rivet::ProjectionApplier::applyProjection ( const Event evt,
const PROJ &  proj 
) const
inlineinherited

Apply the supplied projection on event evt.

Deprecated:
Prefer the simpler apply<> form

◆ applyProjection() [2/3]

template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > Rivet::ProjectionApplier::applyProjection ( const Event evt,
const Projection proj 
) const
inlineinherited

Apply the supplied projection on event evt.

Deprecated:
Prefer the simpler apply<> form

◆ applyProjection() [3/3]

template<typename PROJ = Projection>
std::enable_if_t< std::is_base_of< Projection, PROJ >::value, const PROJ & > Rivet::ProjectionApplier::applyProjection ( const Event evt,
const std::string &  name 
) const
inlineinherited

Apply the named projection on event evt.

Deprecated:
Prefer the simpler apply<> form

References Rivet::ProjectionApplier::name().

◆ beamPairs()

virtual const std::set< PdgIdPair > Rivet::Projection::beamPairs ( ) const
virtualinherited

Return the allowed beam pairs on which this projection can operate, not including recursion. Derived classes should ensure that all contained projections are registered in the _projections set for the beam constraint chaining to work.

Todo:
Remove the beam constraints system from projections.

◆ before()

bool Rivet::Projection::before ( const Projection p) const
inherited

Determine whether this object should be ordered before the object p given as argument. If p is of a different class than this, the before() function of the corresponding type_info objects is used. Otherwise, if the objects are of the same class, the virtual compare(const Projection &) will be returned.

◆ clone()

◆ compare()

CmpState Rivet::Correlators::compare ( const Projection p) const
inlineprotectedvirtual

This function is used to define a unique ordering between different Projection objects of the same class. If this is considered to be equivalent to the Projector object, p, in the argument the function should return 0. If this object should be ordered before p a negative value should be returned, otherwise a positive value should be returned. This function must never be called explicitly, but should only be called from the operator<(const Projection &). When implementing the function in concrete sub-classes, it is then guaranteed that the Projection object p in the argument is of the same class as the sub-class and can be safely dynamically casted to that class.

When implementing this function in a sub-class, the immediate base class version of the function should be called first. If the base class function returns a non-zero value, that value should be returned immediately. Only if zero is returned should this function check the member variables of the sub-class to determine whether this should be ordered before or after p, or if it is equivalent with p.

Implements Rivet::Projection.

References Rivet::Projection::mkPCmp(), and Rivet::Kin::p().

◆ declare() [1/2]

template<typename PROJ >
const PROJ & Rivet::ProjectionApplier::declare ( const PROJ &  proj,
const std::string &  name 
)
inlineprotectedinherited

Register a contained projection (user-facing version)

Todo:
Add SFINAE to require that PROJ inherit from Projection

References Rivet::ProjectionApplier::declareProjection(), and Rivet::ProjectionApplier::name().

Referenced by Rivet::CentralEtHCM::CentralEtHCM(), Rivet::CentralityEstimator::CentralityEstimator(), Rivet::ChargedLeptons::ChargedLeptons(), Rivet::ALICE::CLMultiplicity< INNER >::CLMultiplicity(), Rivet::DISDiffHadron::DISDiffHadron(), Rivet::DISFinalState::DISFinalState(), Rivet::DISKinematics::DISKinematics(), Rivet::DISLepton::DISLepton(), Rivet::GammaGammaFinalState::GammaGammaFinalState(), Rivet::GammaGammaKinematics::GammaGammaKinematics(), Rivet::GammaGammaLeptons::GammaGammaLeptons(), Rivet::GeneratedCentrality::GeneratedCentrality(), Rivet::HadronicFinalState::HadronicFinalState(), Rivet::HeavyHadrons::HeavyHadrons(), Rivet::Hemispheres::Hemispheres(), Rivet::InvisibleFinalState::InvisibleFinalState(), Rivet::LeadingParticlesFinalState::LeadingParticlesFinalState(), Rivet::LossyFinalState< FILTER >::LossyFinalState(), Rivet::MC_pPbMinBiasTrigger::MC_pPbMinBiasTrigger(), Rivet::MC_SumETFwdPbCentrality::MC_SumETFwdPbCentrality(), Rivet::ATLAS::MinBiasTrigger::MinBiasTrigger(), Rivet::MissingMomentum::MissingMomentum(), Rivet::NeutralFinalState::NeutralFinalState(), Rivet::NonHadronicFinalState::NonHadronicFinalState(), Rivet::ParisiTensor::ParisiTensor(), Rivet::PercentileProjection::PercentileProjection(), Rivet::PrimaryHadrons::PrimaryHadrons(), Rivet::SmearedJets::SmearedJets(), Rivet::SmearedMET::SmearedMET(), Rivet::SmearedParticles::SmearedParticles(), Rivet::Spherocity::Spherocity(), Rivet::ATLAS::SumET_PB_Centrality::SumET_PB_Centrality(), Rivet::ATLAS::SumET_PBPB_Centrality::SumET_PBPB_Centrality(), Rivet::TauFinder::TauFinder(), Rivet::TriggerCDFRun0Run1::TriggerCDFRun0Run1(), Rivet::TriggerCDFRun2::TriggerCDFRun2(), Rivet::UndressBeamLeptons::UndressBeamLeptons(), Rivet::ALICE::V0AndTrigger::V0AndTrigger(), Rivet::ALICE::V0Trigger< MODE >::V0Trigger(), Rivet::VetoedFinalState::VetoedFinalState(), Rivet::VisibleFinalState::VisibleFinalState(), Rivet::CentralityProjection::add(), Rivet::VetoedFinalState::addVetoOnThisFinalState(), and Rivet::CentralityBinner< T, MDist >::setProjection().

◆ declare() [2/2]

template<typename PROJ >
const PROJ & Rivet::ProjectionApplier::declare ( const std::string &  name,
const PROJ &  proj 
)
inlineprotectedinherited

Register a contained projection (user-facing, arg-reordered version)

Todo:
Add SFINAE to require that PROJ inherit from Projection

References Rivet::ProjectionApplier::declareProjection(), and Rivet::ProjectionApplier::name().

◆ declareProjection()

template<typename PROJ >
const PROJ & Rivet::ProjectionApplier::declareProjection ( const PROJ &  proj,
const std::string &  name 
)
inlineprotectedinherited

Register a contained projection.

The type of the argument is used to instantiate a new projection internally: this new object is applied to events rather than the argument object. Hence you are advised to only use locally-scoped Projection objects in your Projection and Analysis constructors, and to avoid polymorphism (e.g. handling ConcreteProjection via a pointer or reference to type Projection) since this will screw up the internal type management.

Todo:
Add SFINAE to require that PROJ inherit from Projection

References Rivet::ProjectionApplier::name().

Referenced by Rivet::ProjectionApplier::declare().

◆ get()

template<typename PROJ >
const PROJ & Rivet::ProjectionApplier::get ( const std::string &  name) const
inlineinherited

Get the named projection, specifying return type via a template argument (user-facing alias).

Todo:
Add SFINAE to require that PROJ inherit from Projection

References Rivet::ProjectionApplier::name().

◆ getProjection() [1/2]

template<typename PROJ >
const PROJ & Rivet::ProjectionApplier::getProjection ( const std::string &  name) const
inlineinherited

Get the named projection, specifying return type via a template argument.

Todo:
Add SFINAE to require that PROJ inherit from Projection

References Rivet::ProjectionHandler::getProjection(), Rivet::ProjectionApplier::getProjHandler(), Rivet::ProjectionApplier::name(), and Rivet::Kin::p().

Referenced by Rivet::CentralityProjection::compare(), and Rivet::pcmp().

◆ getProjection() [2/2]

const Projection & Rivet::ProjectionApplier::getProjection ( const std::string &  name) const
inlineinherited

Get the named projection (non-templated, so returns as a reference to a Projection base class).

References Rivet::ProjectionHandler::getProjection(), Rivet::ProjectionApplier::getProjHandler(), and Rivet::ProjectionApplier::name().

◆ hVec()

static vector< int > Rivet::Correlators::hVec ( int  n,
int  m 
)
inlinestatic

Construct a harmonic vectors from n harmonics and m number of particles.

Todo:
In C++14 this can be done much nicer with TMP.

Referenced by Rivet::CumulantAnalysis::bookECorrelator(), and Rivet::CumulantAnalysis::bookECorrelatorGap().

◆ intCorrelatorGap()

const pair< double, double > Rivet::Correlators::intCorrelatorGap ( const Correlators other,
vector< int >  n1,
vector< int >  n2 
) const

Integrated correlator of n1 harmonic, for number of powers n1.size()

This method imposes an eta gap, correlating with another phase space, where another Correlators projection other should be defined. The harmonics of the other phase space is given as n2.

To get e.g. integrated <<4>>_2, n1 should be: n1 = {2, 2} and n2 = {-2, -2}

Referenced by Rivet::CumulantAnalysis::ECorrelator::fill().

◆ mkNamedPCmp()

◆ mkPCmp()

Cmp< Projection > Rivet::Projection::mkPCmp ( const Projection otherparent,
const std::string &  pname 
) const
protectedinherited

Shortcut to make a named Cmp<Projection> comparison with the *this object automatically passed as one of the parent projections.

Note
Alias for mkNamedPCmp

Referenced by Rivet::ALICE::PrimaryParticles::compare(), Rivet::SmearedJets::compare(), Rivet::SmearedMET::compare(), Rivet::SmearedParticles::compare(), and compare().

◆ name()

virtual std::string Rivet::Projection::name ( ) const
inlinevirtualinherited

◆ project()

void Rivet::Correlators::project ( const Event e)
protectedvirtual

Take the information available in the Event and make the calculations necessary to obtain the projection. Note that this function must never be called except inside the Event::applyProjection(Projection *) function.

Implements Rivet::Projection.

◆ pTBinnedCorrelators()

const vector< pair< double, double > > Rivet::Correlators::pTBinnedCorrelators ( vector< int >  n,
bool  overflow = false 
) const

pT differential correlator of n harmonic, for number of powers n.size()

The method can include overflow/underflow bins in the beginning/end of the returned vector, by toggling overflow = true.

Referenced by Rivet::CumulantAnalysis::ECorrelator::fill().

◆ pTBinnedCorrelatorsGap()

const vector< pair< double, double > > Rivet::Correlators::pTBinnedCorrelatorsGap ( const Correlators other,
vector< int >  n1,
vector< int >  n2,
bool  overflow = false 
) const

pT differential correlators of n1 harmonic, for number n1.size()

This method imposes an eta gap, correlating with another phase space, where another Correlators projection other should be defined. The harmonics of the other phase space is given as n2.

To get e.g. differential <<4'>_2, n1 should be: n1 = {2, 2} and n2: n2 = {-2, -2}. To get e.g. differential <<2'>>_4, n1 should be: n1 = {4} and n2: n2 = {-4}. The method can include overflow/underflow bins in the beginning/end of the returned vector, by toggling overflow = true.

Referenced by Rivet::CumulantAnalysis::ECorrelator::fill().


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