Thrust Class Reference

Get the e+ e- thrust basis and the thrust, thrust major and thrust minor scalars. More...

#include <Thrust.hh>

Inheritance diagram for Thrust:

Collaboration diagram for Thrust:

List of all members.

Public Member Functions
	Thrust ()
	Constructor.
	Thrust (const FinalState &fsp)
virtual const Projection *	clone () const
	Clone on the heap.
bool	before (const Projection &p) const
virtual const std::set< PdgIdPair >	beamPairs () const
virtual std::string	name () const
	Get the name of the projection.
Projection &	addPdgIdPair (PdgId beam1, PdgId beam2)
	Add a colliding beam pair.
Log &	getLog () 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.
double	thrust () const
double	thrustMajor () const
	The thrust major scalar, , (thrust along thrust major axis).
double	thrustMinor () const
	The thrust minor scalar, , (thrust along thrust minor axis).
double	oblateness () const
	The oblateness, .
const Vector3 &	thrustAxis () const
const Vector3 &	thrustMajorAxis () const
	The thrust major axis (axis of max thrust perpendicular to thrust axis).
const Vector3 &	thrustMinorAxis () const
	The thrust minor axis (axis perpendicular to thrust and thrust major).
const Vector3 &	axis1 () const
	AxesDefinition axis accessors.
const Vector3 &	axis2 () const
	The 2nd most significant ("major") axis.
const Vector3 &	axis3 () const
	The least significant ("minor") axis.
Direct methods
Ways to do the calculation directly, without engaging the caching system
void	calc (const FinalState &fs)
	Manually calculate the thrust, without engaging the caching system.
void	calc (const vector< Particle > &fsparticles)
	Manually calculate the thrust, without engaging the caching system.
void	calc (const vector< FourMomentum > &fsmomenta)
	Manually calculate the thrust, without engaging the caching system.
void	calc (const vector< Vector3 > &threeMomenta)
	Manually calculate the thrust, without engaging the caching system.
Projection "getting" functions
std::set< ConstProjectionPtr >	getProjections () const
	Get the contained projections, including recursion.
template<typename PROJ >
const PROJ &	getProjection (const std::string &name) const
	Get the named projection, specifying return type via a template argument.
const Projection &	getProjection (const std::string &name) const
Projection applying functions
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const PROJ &proj) const
	Apply the supplied projection on event.
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const Projection &proj) const
	Apply the supplied projection on event.
template<typename PROJ >
const PROJ &	applyProjection (const Event &evt, const std::string &name) const
	Apply the named projection on event.
Protected Member Functions
void	project (const Event &e)
	Perform the projection on the Event.
int	compare (const Projection &p) const
	Compare projections.
Cmp< Projection >	mkNamedPCmp (const Projection &otherparent, const std::string &pname) const
Cmp< Projection >	mkPCmp (const Projection &otherparent, const std::string &pname) const
ProjectionHandler &	getProjHandler () const
	Get a reference to the ProjectionHandler for this thread.
Projection registration functions
template<typename PROJ >
const PROJ &	addProjection (const PROJ &proj, const std::string &name)
const Projection &	_addProjection (const Projection &proj, const std::string &name)
	Untemplated function to do the work...
Protected Attributes
bool	_allowProjReg
	Flag to forbid projection registration in analyses until the init phase.
Private Member Functions
void	_calcThrust (const vector< Vector3 > &fsmomenta)
	Explicitly calculate the thrust values.
Private Attributes
vector< double >	_thrusts
	The thrust scalars.
vector< Vector3 >	_thrustAxes
	The thrust axes.
Friends
class	Event
	Event is a friend.
class	Cmp< Projection >
	The Cmp specialization for Projection is a friend.

---

Detailed Description

Get the e+ e- thrust basis and the thrust, thrust major and thrust minor scalars.

Author:

Andy Buckley

The scalar (maximum) thrust is defined as

, with the direction of the unit vector which maximises being identified as the thrust axis. The unit vector which maximises the thrust scalar in the plane perpendicular to is the "thrust major" direction, and the vector perpendicular to both the thrust and thrust major directions is the thrust minor. Both the major and minor directions have associated thrust scalars.

Thrust calculations have particularly simple forms for less than 4 particles, and in those cases this projection is computationally minimal. For 4 or more particles, a more general calculation must be carried out, based on the Brandt/Dahmen method from Z. Phys. C1 (1978). While a polynomial improvement on the exponential scaling of the naive method, this algorithm scales asymptotically as . Be aware that the thrust may easily be the most computationally demanding projection in Rivet for large events!

The Rivet implementation of thrust is based heavily on Stefan Gieseke's Herwig++ re-coding of the 'tasso' code from HERWIG.

NB. special case with >= 4 coplanar particles will still fail. NB. Thrust assumes all momenta are in the CoM system: no explicit boost is performed. This can be dealt with by appropriate choice of the supplied FinalState.

Definition at line 44 of file Thrust.hh.

---

Constructor & Destructor Documentation

Thrust

(

)

[inline]

Constructor.

Definition at line 48 of file Thrust.hh.

Referenced by Thrust::clone().

{}

Thrust

(

const FinalState &

fsp

)

[inline]

Definition at line 50 of file Thrust.hh.

References ProjectionApplier::addProjection(), and Projection::setName().

                                  {
      setName("Thrust");
      addProjection(fsp, "FS");
    }

---

Member Function Documentation

const Projection & _addProjection	(	const Projection &	proj,
		const std::string &	name
	)		[protected, inherited]

Untemplated function to do the work...

Definition at line 33 of file ProjectionApplier.cc.

References ProjectionApplier::_allowProjReg, ProjectionApplier::getProjHandler(), ProjectionApplier::name(), Projection::name(), and ProjectionHandler::registerProjection().

Referenced by ProjectionApplier::addProjection().

                                                                             {
    if (!_allowProjReg) {
      cerr << "Trying to register projection '"
           << proj.name() << "' before init phase in '" << this->name() << "'." << endl;
      exit(2);
    }
    const Projection& reg = getProjHandler().registerProjection(*this, proj, name);
    return reg;
  }

void _calcThrust

(

const vector< Vector3 > &

fsmomenta

)

[private]

Explicitly calculate the thrust values.

Todo:

Improve this --- special directions bad... (a,b,c) _|_ 1/(a^2+b^2) (b,-a,0) etc., but which combination minimises error?

Definition at line 101 of file Thrust.cc.

References Rivet::_calcT(), Thrust::_thrustAxes, Thrust::_thrusts, Vector3::cross(), Rivet::cross(), Rivet::dot(), Vector< N >::mod(), MSG_DEBUG, Vector3::unit(), Vector3::x(), and Vector3::z().

Referenced by Thrust::calc().

                                                           {
    // Make a vector of the three-momenta in the final state
    double momentumSum(0.0);
    foreach (const Vector3& p3, fsmomenta) {
      momentumSum += p3.mod();
    }
    MSG_DEBUG("Number of particles = " << fsmomenta.size());


    // Clear the caches
    _thrusts.clear();
    _thrustAxes.clear();


    // If there are fewer than 2 visible particles, we can't do much
    if (fsmomenta.size() < 2) {
      for (int i = 0; i < 3; ++i) {
        _thrusts.push_back(-1);
        _thrustAxes.push_back(Vector3(0,0,0));
      }
      return;
    }


    // Handle special case of thrust = 1 if there are only 2 particles
    if (fsmomenta.size() == 2) {
      Vector3 axis(0,0,0);
      _thrusts.push_back(1.0);
      _thrusts.push_back(0.0);
      _thrusts.push_back(0.0);
      axis = fsmomenta[0].unit();
      if (axis.z() < 0) axis = -axis;
      _thrustAxes.push_back(axis);
      /// @todo Improve this --- special directions bad...
      /// (a,b,c) _|_ 1/(a^2+b^2) (b,-a,0) etc., but which combination minimises error?
      if (axis.z() < 0.75)
        _thrustAxes.push_back( (axis.cross(Vector3(0,0,1))).unit() );
      else
        _thrustAxes.push_back( (axis.cross(Vector3(0,1,0))).unit() );
      _thrustAxes.push_back( _thrustAxes[0].cross(_thrustAxes[1]) );
      return;
    }



    // Temporary variables for calcs
    Vector3 axis(0,0,0);
    double val = 0.;

    // Get thrust
    _calcT(fsmomenta, val, axis);
    MSG_DEBUG("Mom sum = " << momentumSum);
    _thrusts.push_back(val / momentumSum);
    // Make sure that thrust always points along the +ve z-axis.
    if (axis.z() < 0) axis = -axis;
    axis = axis.unit();
    MSG_DEBUG("Axis = " << axis);
    _thrustAxes.push_back(axis);

    // Get thrust major
    vector<Vector3> threeMomenta;
    foreach (const Vector3& v, fsmomenta) {
      // Get the part of each 3-momentum which is perpendicular to the thrust axis
      const Vector3 vpar = dot(v, axis.unit()) * axis.unit();
      threeMomenta.push_back(v - vpar);
    }
    _calcT(threeMomenta, val, axis);
    _thrusts.push_back(val / momentumSum);
    if (axis.x() < 0) axis = -axis;
    axis = axis.unit();
    _thrustAxes.push_back(axis);

    // Get thrust minor
    if (_thrustAxes[0].dot(_thrustAxes[1]) < 1e-10) {
      axis = _thrustAxes[0].cross(_thrustAxes[1]);
      _thrustAxes.push_back(axis);
      val = 0.0;
      foreach (const Vector3& v, fsmomenta) {
        val += fabs(dot(axis, v));
      }
      _thrusts.push_back(val / momentumSum);
    } else {
      _thrusts.push_back(-1.0);
      _thrustAxes.push_back(Vector3(0,0,0));
    }

  }

Projection& addPdgIdPair	(	PdgId	beam1,
		PdgId	beam2
	)		[inline, inherited]

Add a colliding beam pair.

Definition at line 107 of file Projection.hh.

References Projection::_beamPairs.

Referenced by Projection::Projection().

                                                       {
      _beamPairs.insert(PdgIdPair(beam1, beam2));
      return *this;
    }

const PROJ& addProjection	(	const PROJ &	proj,
		const std::string &	name
	)		[inline, protected, inherited]

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.

Definition at line 116 of file ProjectionApplier.hh.

References ProjectionApplier::_addProjection().

Referenced by FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), VetoedFinalState::addVetoOnThisFinalState(), BeamThrust::BeamThrust(), CDF_2009_S8057893::CDF_2009_S8057893::init(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), CMS_2010_S8547297::init(), CMS_2010_S8656010::init(), ATLAS_2010_S8894728::init(), ALICE_2012_I1181770::init(), ATLAS_2011_S8994773::init(), CMS_2011_S8950903::init(), CMS_2012_PAS_QCD_11_010::init(), CMS_2011_S8968497::init(), CMS_2012_I1193338::init(), CMS_2012_I1184941::init(), CMSTOTEM_2014_I1294140::init(), ATLAS_2014_I1298811::init(), ATLAS_2011_I894867::init(), ALICE_2011_S8909580::init(), LHCB_2013_I1208105::init(), CMS_2011_S8941262::init(), CMS_2011_S8973270::init(), CMS_2012_I1090423::init(), ATLAS_2010_CONF_2010_049::init(), CDF_2007_S7057202::init(), CMS_2011_S9086218::init(), TOTEM_2012_I1115294::init(), CMS_2011_S9215166::init(), CMS_2012_I1087342::init(), ATLAS_2014_I1282441::init(), ALICE_2011_S8945144::init(), ATLAS_2012_I1091481::init(), MC_MUONS::init(), MC_TAUS::init(), ATLAS_2010_S8591806::init(), ATLAS_2011_S9002537::init(), CMS_2011_I954992::init(), CMS_2011_S8957746::init(), CMS_2011_S8978280::init(), TOTEM_2012_002::init(), CMS_2011_S9088458::init(), CMS_2013_I1265659::init(), D0_2011_I895662::init(), ATLAS_2011_I930220::init(), LHCF_2012_I1115479::init(), MC_ELECTRONS::init(), CDF_2012_NOTE10874::init(), CMS_2013_I1273574::init(), D0_2010_S8570965::init(), ATLAS_2011_I925932::init(), MC_JETTAGS::init(), ATLAS_2012_I1124167::init(), STAR_2006_S6500200::init(), STAR_2008_S7993412::init(), CMS_2013_I1256943::init(), CMS_QCD_10_024::init(), UA5_1987_S1640666::init(), CDF_1997_S3541940::init(), MC_HINC::init(), MC_JETS::init(), MC_KTSPLITTINGS::init(), CDF_2005_S6080774::init(), CDF_2006_S6450792::init(), SFM_1984_S1178091::init(), ARGUS_1993_S2789213::init(), ATLAS_2011_S9128077::init(), CMS_2013_I1261026::init(), CMS_2013_I1258128::init(), CMS_2014_I1298810::init(), D0_2000_S4480767::init(), BABAR_2007_S7266081::init(), BELLE_2008_I786560::init(), CDF_1993_S2742446::init(), CDF_2000_S4155203::init(), MC_WWINC::init(), CDF_2008_S7782535::init(), MC_ZJETS::init(), CDF_2008_S8093652::init(), MC_ZZINC::init(), STAR_2006_S6870392::init(), UA5_1982_S875503::init(), UA5_1989_S1926373::init(), ATLAS_2013_I1190187::init(), CMS_2013_I1209721::init(), D0_2008_S6879055::init(), ATLAS_2012_I1204447::init(), ARGUS_1993_S2669951::init(), LHCB_2012_I1208102::init(), ATLAS_2011_I944826::init(), MC_ZINC::init(), ALEPH_1991_S2435284::init(), MC_WINC::init(), MC_WJETS::init(), MC_WWKTSPLITTINGS::init(), MC_ZZJETS::init(), ATLAS_2012_I1183818::init(), ATLAS_2011_S9035664::init(), ATLAS_2012_I1204784::init(), CMS_2012_I1102908::init(), ATLAS_2013_I1219109::init(), CMS_2013_I1272853::init(), ATLAS_2013_I1243871::init(), ATLAS_2014_I1268975::init(), D0_2000_I499943::init(), ATLAS_2014_I1325553::init(), E735_1998_S3905616::init(), MC_ZZKTSPLITTINGS::init(), MC_DIJET::init(), MC_DIPHOTON::init(), MC_ZKTSPLITTINGS::init(), MC_HJETS::init(), MC_HKTSPLITTINGS::init(), CDF_2005_S6217184::init(), MC_PHOTONINC::init(), MC_WKTSPLITTINGS::init(), ATLAS_2011_S8924791::init(), OPAL_1995_S3198391::init(), OPAL_1996_S3257789::init(), OPAL_1997_S3608263::init(), OPAL_1998_S3702294::init(), OPAL_1998_S3749908::init(), OPAL_2000_S4418603::init(), ATLAS_2010_S8817804::init(), PDG_TAUS::init(), D0_1996_S3324664::init(), D0_2007_S7075677::init(), D0_2009_S8202443::init(), D0_2010_S8821313::init(), ALEPH_2002_S4823664::init(), DELPHI_1999_S3960137::init(), ALEPH_1996_S3196992::init(), UA5_1986_S1583476::init(), EXAMPLE_CUTS::init(), CDF_1990_S2089246::init(), CDF_1988_S1865951::init(), ATLAS_2011_I945498::init(), MC_PHOTONKTSPLITTINGS::init(), MC_PHOTONJETS::init(), MC_WWJETS::init(), CMS_2012_I1107658::init(), D0_2008_S7554427::init(), D0_2008_S7863608::init(), D0_2010_S8671338::init(), MC_IDENTIFIED::init(), H1_1995_S3167097::init(), CDF_1994_S2952106::init(), ATLAS_2011_I954993::init(), MC_GENERIC::init(), MC_LEADJETUE::init(), MC_TTBAR::init(), CDF_2008_S7828950::init(), ATLAS_2012_I1180197::init(), MC_PHOTONJETUE::init(), ATLAS_2012_I1188891::init(), D0_2008_S7837160::init(), CDF_2008_S7540469::init(), ATLAS_2011_S9131140::init(), D0_2001_S4674421::init(), ATLAS_2011_S9212353::init(), ATLAS_2014_I1300647::init(), ATLAS_2014_I1306615::init(), ATLAS_2014_I1307756::init(), ATLAS_2012_CONF_2012_104::init(), ATLAS_2012_CONF_2012_105::init(), UA1_1990_S2044935::init(), ZEUS_2001_S4815815::init(), LHCB_2011_I919315::init(), ATLAS_2012_I1119557::init(), ATLAS_2012_I1095236::init(), JADE_OPAL_2000_S4300807::init(), D0_1996_S3214044::init(), ATLAS_2011_CONF_2011_098::init(), ALICE_2010_S8624100::init(), ATLAS_2012_CONF_2012_103::init(), CDF_1996_S3108457::init(), EXAMPLE::init(), ATLAS_2012_CONF_2012_109::init(), ATLAS_2012_I1126136::init(), ATLAS_2011_I926145::init(), OPAL_1997_S3396100::init(), CDF_2008_S8095620::init(), MC_QCD_PARTONS::init(), ATLAS_2011_S8983313::init(), ATLAS_2013_I1217863_Z::init(), ATLAS_2012_I1186556::init(), STAR_2006_S6860818::init(), ATLAS_2013_I1217863_W::init(), CMS_2012_I941555::init(), CMS_2014_I1303894::init(), DELPHI_1995_S3137023::init(), STAR_2009_UE_HELEN::init(), JADE_1998_S3612880::init(), ALEPH_2004_S5765862::init(), ATLAS_2012_I1112263::init(), ATLAS_2012_I1125961::init(), ATLAS_2012_I1190891::init(), LHCB_2013_I1218996::init(), CMS_2013_I1218372::init(), ALICE_2010_S8706239::init(), ATLAS_2012_I1125575::init(), ATLAS_2014_I1288706::init(), ATLAS_2011_S9212183::init(), ATLAS_2012_CONF_2012_001::init(), MC_HFJETS::init(), ALICE_2010_S8625980::init(), CDF_2009_NOTE_9936::init(), CDF_1996_S3349578::init(), D0_2009_S8320160::init(), CDF_2000_S4266730::init(), CDF_2001_S4563131::init(), MC_SUSY::init(), ATLAS_2011_S8971293::init(), CDF_2009_S8383952::init(), CDF_2009_S8436959::init(), CDF_1998_S3618439::init(), D0_2009_S8349509::init(), ATLAS_2013_I1263495::init(), ALEPH_2001_S4656318::init(), D0_2006_S6438750::init(), CDF_2001_S4517016::init(), SLD_2002_S4869273::init(), ATLAS_2010_S8919674::init(), ATLAS_2012_I1082936::init(), D0_2010_S8566488::init(), DELPHI_2002_069_CONF_603::init(), LHCB_2011_I917009::init(), SLD_1996_S3398250::init(), ATLAS_2011_I921594::init(), UA5_1988_S1867512::init(), CDF_1996_S3418421::init(), MC_PHOTONS::init(), LHCB_2014_I1281685::init(), LHCB_2012_I1119400::init(), D0_2008_S7662670::init(), ATLAS_2014_I1306294::init(), CDF_2006_S6653332::init(), ATLAS_2014_I1304688::init(), MC_WPOL::init(), CDF_2008_S7541902::init(), ATLAS_2011_S9108483::init(), ATLAS_2013_I1230812::init(), ATLAS_2012_I1199269::init(), ATLAS_2014_I1307243::init(), ATLAS_2012_I1117704::init(), ATLAS_2011_S9225137::init(), ATLAS_2011_S9019561::init(), ATLAS_2012_I1082009::init(), ATLAS_2011_CONF_2011_090::init(), ATLAS_2012_I943401::init(), ATLAS_2012_I946427::init(), OPAL_2002_S5361494::init(), DELPHI_2000_S4328825::init(), TASSO_1990_S2148048::init(), ATLAS_2012_I1084540::init(), CMS_2013_I1224539_WJET::init(), CMS_2013_I1224539_DIJET::init(), ATLAS_2012_I1083318::init(), ATLAS_2012_CONF_2012_153::init(), CDF_2009_S8233977::init(), ATLAS_2013_I1217867::init(), CMS_2013_I1224539_ZJET::init(), CDF_2010_S8591881_DY::init(), CDF_2010_S8591881_QCD::init(), D0_2004_S5992206::init(), ATLAS_2010_S8914702::init(), ATLAS_2014_I1312627::init(), ATLAS_2011_S9120807::init(), CDF_2001_S4751469::init(), D0_2008_S7719523::init(), ATLAS_2010_S8918562::init(), ATLAS_2011_S9041966::init(), ATLAS_2012_I1094568::init(), DELPHI_1996_S3430090::init(), ALEPH_1996_S3486095::init(), ALEPH_1999_S4193598::init(), ATLAS_2012_I1093738::init(), OPAL_2004_S6132243::init(), OPAL_1994_S2927284::init(), BELLE_2013_I1216515::init(), STAR_2008_S7869363::init(), BABAR_2007_S6895344::init(), ATLAS_2011_S9126244::init(), MC_VH2BB::init(), BELLE_2001_S4598261::init(), BABAR_2005_S6181155::init(), ATLAS_2011_I919017::init(), OPAL_2001_S4553896::init(), BABAR_2013_I1238276::init(), DELPHI_2003_WUD_03_11::init(), ATLAS_2012_I1203852::init(), CDF_2004_S5839831::init(), BABAR_2003_I593379::init(), ATLAS_2012_I1094061::init(), ARGUS_1993_S2653028::init(), CLEO_2004_S5809304::init(), OPAL_1998_S3780481::init(), ATLAS_2012_I1093734::init(), OPAL_1993_S2692198::init(), ATLAS_2014_I1279489::init(), ATLAS_2012_I1094564::init(), H1_1994_S2919893::init(), H1_2000_S4129130::init(), SLD_2004_S5693039::init(), SLD_1999_S3743934::init(), PDG_HADRON_MULTIPLICITIES::init(), PDG_HADRON_MULTIPLICITIES_RATIOS::init(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), PrimaryHadrons::PrimaryHadrons(), PromptFinalState::PromptFinalState(), Sphericity::Sphericity(), Spherocity::Spherocity(), TauFinder::TauFinder(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder().

                                                                       {
      const Projection& reg = _addProjection(proj, name);
      const PROJ& rtn = dynamic_cast<const PROJ&>(reg);
      return rtn;
    }

const PROJ& applyProjection	(	const Event &	evt,
		const PROJ &	proj
	)		const [inline, inherited]

Apply the supplied projection on event.

Definition at line 70 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

Referenced by DISFinalState::project().

                                                                          {
      return pcast<PROJ>(_applyProjection(evt, proj));
    }

const PROJ& applyProjection	(	const Event &	evt,
		const Projection &	proj
	)		const [inline, inherited]

Apply the supplied projection on event.

Definition at line 77 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

                                                                                {
      return pcast<PROJ>(_applyProjection(evt, proj));
    }

const PROJ& applyProjection	(	const Event &	evt,
		const std::string &	name
	)		const [inline, inherited]

Apply the named projection on event.

Definition at line 84 of file ProjectionApplier.hh.

References ProjectionApplier::_applyProjection().

                                                                               {
      return pcast<PROJ>(_applyProjection(evt, name));
    }

const Vector3& axis1

(

)

const [inline, virtual]

AxesDefinition axis accessors.

Implements AxesDefinition.

Definition at line 98 of file Thrust.hh.

References Thrust::thrustAxis().

{ return thrustAxis(); }

const Vector3& axis2

(

)

const [inline, virtual]

The 2nd most significant ("major") axis.

Implements AxesDefinition.

Definition at line 99 of file Thrust.hh.

References Thrust::thrustMajorAxis().

{ return thrustMajorAxis(); }

const Vector3& axis3

(

)

const [inline, virtual]

The least significant ("minor") axis.

Implements AxesDefinition.

Definition at line 100 of file Thrust.hh.

References Thrust::thrustMinorAxis().

{ return thrustMinorAxis(); }

const set< PdgIdPair > beamPairs

(

)

const [virtual, inherited]

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.

Definition at line 33 of file Projection.cc.

References Projection::_beamPairs, Projection::beamPairs(), Projection::getLog(), ProjectionApplier::getProjections(), Rivet::intersection(), and Log::TRACE.

Referenced by Projection::beamPairs().

                                                   {
    set<PdgIdPair> ret = _beamPairs;
    set<ConstProjectionPtr> projs = getProjections();
    for (set<ConstProjectionPtr>::const_iterator ip = projs.begin(); ip != projs.end(); ++ip) {
      ConstProjectionPtr p = *ip;
      getLog() << Log::TRACE << "Proj addr = " << p << endl;
      if (p) ret = intersection(ret, p->beamPairs());
    }
    return ret;
  }

bool 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.

Definition at line 22 of file Projection.cc.

References Projection::compare().

Referenced by less< const Rivet::Projection * >::operator()().

                                                   {
    const std::type_info& thisid = typeid(*this);
    const std::type_info& otherid = typeid(p);
    if (thisid == otherid) {
      return compare(p) < 0;
    } else {
      return thisid.before(otherid);
    }
  }

void calc

(

const FinalState &

fs

)

Manually calculate the thrust, without engaging the caching system.

Definition at line 10 of file Thrust.cc.

References ParticleFinder::particles().

Referenced by CMS_2011_S8957746::analyze(), CMS_2013_I1209721::analyze(), ATLAS_2012_I1124167::analyze(), and Thrust::project().

                                        {
    calc(fs.particles());
  }

void calc

(

const vector< Particle > &

fsparticles

)

Manually calculate the thrust, without engaging the caching system.

Definition at line 14 of file Thrust.cc.

References Thrust::_calcThrust(), Particle::momentum(), and FourVector::vector3().

                                                       {
    vector<Vector3> threeMomenta;
    threeMomenta.reserve(fsparticles.size());
    foreach (const Particle& p, fsparticles) {
      threeMomenta.push_back( p.momentum().vector3() );
    }
    _calcThrust(threeMomenta);
  }

void calc

(

const vector< FourMomentum > &

fsmomenta

)

Manually calculate the thrust, without engaging the caching system.

Definition at line 23 of file Thrust.cc.

References Thrust::_calcThrust(), and FourVector::vector3().

                                                         {
    vector<Vector3> threeMomenta;
    threeMomenta.reserve(fsmomenta.size());
    foreach (const FourMomentum& v, fsmomenta) {
      threeMomenta.push_back(v.vector3());
    }
    _calcThrust(threeMomenta);
  }

void calc

(

const vector< Vector3 > &

threeMomenta

)

Manually calculate the thrust, without engaging the caching system.

Definition at line 32 of file Thrust.cc.

References Thrust::_calcThrust().

                                                    {
    _calcThrust(fsmomenta);
  }

virtual const Projection* clone

(

)

const [inline, virtual]

Clone on the heap.

Implements AxesDefinition.

Definition at line 56 of file Thrust.hh.

References Thrust::Thrust().

                                            {
      return new Thrust(*this);
    }

int compare

(

const Projection &

p

)

const [inline, protected, virtual]

Compare projections.

Implements Projection.

Definition at line 70 of file Thrust.hh.

References Projection::mkNamedPCmp().

                                           {
      return mkNamedPCmp(p, "FS");
    }

Log& getLog

(

)

const [inline, inherited]

Get a Log object based on the getName() property of the calling projection object.

Reimplemented from ProjectionApplier.

Definition at line 114 of file Projection.hh.

References Projection::name().

Referenced by Projection::beamPairs(), InvMassFinalState::calc(), ChargedFinalState::project(), InitialQuarks::project(), PromptFinalState::project(), LossyFinalState< ConstRandomFilter >::project(), UnstableFinalState::project(), and VetoedFinalState::project().

                        {
      string logname = "Rivet.Projection." + name();
      return Log::getLog(logname);
    }

const PROJ& getProjection

(

const std::string &

name

)

const [inline, inherited]

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

Definition at line 52 of file ProjectionApplier.hh.

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

Referenced by ProjectionApplier::_applyProjection(), Rivet::pcmp(), and Hemispheres::project().

                                                           {
      const Projection& p = getProjHandler().getProjection(*this, name);
      return pcast<PROJ>(p);
    }

const Projection& getProjection

(

const std::string &

name

)

const [inline, inherited]

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

Definition at line 60 of file ProjectionApplier.hh.

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

                                                                 {
      return getProjHandler().getProjection(*this, name);
    }

std::set<ConstProjectionPtr> getProjections

(

)

const [inline, inherited]

Get the contained projections, including recursion.

Definition at line 45 of file ProjectionApplier.hh.

References ProjectionHandler::DEEP, ProjectionHandler::getChildProjections(), and ProjectionApplier::getProjHandler().

Referenced by Projection::beamPairs().

                                                      {
      return getProjHandler().getChildProjections(*this, ProjectionHandler::DEEP);
    }

ProjectionHandler& getProjHandler

(

)

const [inline, protected, inherited]

Get a reference to the ProjectionHandler for this thread.

Definition at line 97 of file ProjectionApplier.hh.

References ProjectionApplier::_projhandler.

Referenced by ProjectionApplier::_addProjection(), ProjectionApplier::getProjection(), ProjectionApplier::getProjections(), and ProjectionApplier::~ProjectionApplier().

                                              {
      return _projhandler;
    }

Cmp< Projection > mkNamedPCmp	(	const Projection &	otherparent,
		const std::string &	pname
	)		const [protected, inherited]

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

Definition at line 45 of file Projection.cc.

References Rivet::pcmp().

Referenced by BeamThrust::compare(), ChargedLeptons::compare(), FParameter::compare(), CentralEtHCM::compare(), MergedFinalState::compare(), ChargedFinalState::compare(), DISLepton::compare(), DISKinematics::compare(), DISFinalState::compare(), VisibleFinalState::compare(), TauFinder::compare(), PromptFinalState::compare(), NeutralFinalState::compare(), Spherocity::compare(), LeadingParticlesFinalState::compare(), ParisiTensor::compare(), FoxWolframMoments::compare(), MissingMomentum::compare(), LossyFinalState< ConstRandomFilter >::compare(), Thrust::compare(), InvMassFinalState::compare(), Sphericity::compare(), HeavyHadrons::compare(), ZFinder::compare(), Hemispheres::compare(), DressedLeptons::compare(), WFinder::compare(), IdentifiedFinalState::compare(), VetoedFinalState::compare(), JetShape::compare(), and FastJets::compare().

                                                                     {
    return pcmp(*this, otherparent, pname);
  }

Cmp< Projection > mkPCmp	(	const Projection &	otherparent,
		const std::string &	pname
	)		const [protected, inherited]

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

Definition at line 51 of file Projection.cc.

References Rivet::pcmp().

                                                                {
    return pcmp(*this, otherparent, pname);
  }

virtual std::string name

(

)

const [inline, virtual, inherited]

Get the name of the projection.

Implements ProjectionApplier.

Definition at line 101 of file Projection.hh.

References Projection::_name.

Referenced by ProjectionApplier::_addProjection(), ProjectionHandler::_checkDuplicate(), ProjectionHandler::_clone(), ProjectionHandler::_getEquiv(), VetoedFinalState::addVetoOnThisFinalState(), Projection::getLog(), ProjectionHandler::registerProjection(), and Projection::setName().

                                   {
      return _name;
    }

double oblateness

(

)

const [inline]

The oblateness, .

Definition at line 85 of file Thrust.hh.

References Thrust::_thrusts.

Referenced by OPAL_2004_S6132243::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _thrusts[1] - _thrusts[2]; }

void project

(

const Event &

e

)

[inline, protected, virtual]

Perform the projection on the Event.

Implements Projection.

Definition at line 63 of file Thrust.hh.

References Thrust::calc(), and Rivet::particles().

                                 {
      const vector<Particle> ps
        = applyProjection<FinalState>(e, "FS").particles();
      calc(ps);
    }

void setName

(

const std::string &

name

)

[inline, inherited]

Used by derived classes to set their name.

Definition at line 120 of file Projection.hh.

References Projection::_name, and Projection::name().

Referenced by FastJets::_init1(), FastJets::_init2(), FastJets::_init3(), Beam::Beam(), BeamThrust::BeamThrust(), CentralEtHCM::CentralEtHCM(), ChargedFinalState::ChargedFinalState(), ChargedLeptons::ChargedLeptons(), ConstLossyFinalState::ConstLossyFinalState(), DISFinalState::DISFinalState(), DISKinematics::DISKinematics(), DISLepton::DISLepton(), DressedLeptons::DressedLeptons(), FinalState::FinalState(), FoxWolframMoments::FoxWolframMoments(), FParameter::FParameter(), HadronicFinalState::HadronicFinalState(), HeavyHadrons::HeavyHadrons(), Hemispheres::Hemispheres(), IdentifiedFinalState::IdentifiedFinalState(), InitialQuarks::InitialQuarks(), JetAlg::JetAlg(), JetShape::JetShape(), LeadingParticlesFinalState::LeadingParticlesFinalState(), LossyFinalState< ConstRandomFilter >::LossyFinalState(), MergedFinalState::MergedFinalState(), MissingMomentum::MissingMomentum(), NeutralFinalState::NeutralFinalState(), NonHadronicFinalState::NonHadronicFinalState(), ParisiTensor::ParisiTensor(), PrimaryHadrons::PrimaryHadrons(), PromptFinalState::PromptFinalState(), Sphericity::Sphericity(), Spherocity::Spherocity(), TauFinder::TauFinder(), Thrust::Thrust(), TriggerCDFRun0Run1::TriggerCDFRun0Run1(), TriggerCDFRun2::TriggerCDFRun2(), TriggerUA5::TriggerUA5(), UnstableFinalState::UnstableFinalState(), VetoedFinalState::VetoedFinalState(), VisibleFinalState::VisibleFinalState(), WFinder::WFinder(), and ZFinder::ZFinder().

                                        {
      _name = name;
    }

double thrust

(

)

const [inline]

Thrust scalar accessors The thrust scalar, , (maximum thrust).

Definition at line 79 of file Thrust.hh.

References Thrust::_thrusts.

Referenced by CMS_2011_S8957746::analyze(), CMS_2013_I1209721::analyze(), EXAMPLE::analyze(), JADE_1998_S3612880::analyze(), ATLAS_2012_I1124167::analyze(), TASSO_1990_S2148048::analyze(), OPAL_2004_S6132243::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _thrusts[0]; }

const Vector3& thrustAxis

(

)

const [inline]

Thrust axis accessors The thrust axis.

Definition at line 90 of file Thrust.hh.

References Thrust::_thrustAxes.

Referenced by ALEPH_1996_S3196992::analyze(), JADE_1998_S3612880::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), and Thrust::axis1().

{ return _thrustAxes[0]; }

double thrustMajor

(

)

const [inline]

The thrust major scalar, , (thrust along thrust major axis).

Definition at line 81 of file Thrust.hh.

References Thrust::_thrusts.

Referenced by CMS_2011_S8957746::analyze(), EXAMPLE::analyze(), ATLAS_2012_I1124167::analyze(), OPAL_2004_S6132243::analyze(), ALEPH_2004_S5765862::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _thrusts[1]; }

const Vector3& thrustMajorAxis

(

)

const [inline]

The thrust major axis (axis of max thrust perpendicular to thrust axis).

Definition at line 92 of file Thrust.hh.

References Thrust::_thrustAxes.

Referenced by ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), and Thrust::axis2().

{ return _thrustAxes[1]; }

double thrustMinor

(

)

const [inline]

The thrust minor scalar, , (thrust along thrust minor axis).

Definition at line 83 of file Thrust.hh.

References Thrust::_thrusts.

Referenced by OPAL_2004_S6132243::analyze(), ALEPH_1996_S3486095::analyze(), ALEPH_2004_S5765862::analyze(), and DELPHI_1996_S3430090::analyze().

{ return _thrusts[2]; }

const Vector3& thrustMinorAxis

(

)

const [inline]

The thrust minor axis (axis perpendicular to thrust and thrust major).

Definition at line 94 of file Thrust.hh.

References Thrust::_thrustAxes.

Referenced by ALEPH_2004_S5765862::analyze(), DELPHI_1996_S3430090::analyze(), and Thrust::axis3().

{ return _thrustAxes[2]; }

---

Friends And Related Function Documentation

friend class Cmp< Projection > [friend, inherited]

The Cmp specialization for Projection is a friend.

Definition at line 35 of file Projection.hh.

friend class Event [friend, inherited]

Event is a friend.

Definition at line 32 of file Projection.hh.

---

Member Data Documentation

bool _allowProjReg [protected, inherited]

Flag to forbid projection registration in analyses until the init phase.

Definition at line 143 of file ProjectionApplier.hh.

Referenced by ProjectionApplier::_addProjection(), and Analysis::Analysis().

vector<Vector3> _thrustAxes [private]

The thrust axes.

Definition at line 131 of file Thrust.hh.

Referenced by Thrust::_calcThrust(), Thrust::thrustAxis(), Thrust::thrustMajorAxis(), and Thrust::thrustMinorAxis().

vector<double> _thrusts [private]

The thrust scalars.

Definition at line 128 of file Thrust.hh.

Referenced by Thrust::_calcThrust(), Thrust::oblateness(), Thrust::thrust(), Thrust::thrustMajor(), and Thrust::thrustMinor().

---

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

• Thrust.hh
• Thrust.cc