rivet is hosted by Hepforge, IPPP Durham
Classes | Public Member Functions | Protected Member Functions | Protected Attributes | Friends
FourMomentum Class Reference

Specialized version of the FourVector with momentum/energy functionality. More...

#include <Vector4.hh>

Inheritance diagram for FourMomentum:
Inheritance graph
[legend]
Collaboration diagram for FourMomentum:
Collaboration graph
[legend]

List of all members.

Classes

struct  byEAscending
 Struct for sorting by increasing energy. More...
struct  byEDescending
 Struct for sorting by decreasing energy. More...

Public Member Functions

 FourMomentum ()
template<typename V4 >
 FourMomentum (const V4 &other)
 FourMomentum (const Vector< 4 > &other)
 FourMomentum (const double E, const double px, const double py, const double pz)
 ~FourMomentum ()
double E () const
 Get energy $ E $ (time component of momentum).
Vector3 p3 () const
 Get 3-momentum part, $ p $.
double px () const
 Get 3-momentum part, $ p $ (alias).
double py () const
 Get y-component of momentum $ p_y $.
double pz () const
 Get z-component of momentum $ p_z $.
FourMomentum & setE (double E)
 Set energy $ E $ (time component of momentum).
FourMomentum & setPx (double px)
 Set x-component of momentum $ p_x $.
FourMomentum & setPy (double py)
 Set y-component of momentum $ p_y $.
FourMomentum & setPz (double pz)
 Set z-component of momentum $ p_z $.
double mass () const
 Get the mass $ m = \sqrt{E^2 - p^2} $ (the Lorentz self-invariant).
double mass2 () const
 Get the squared mass $ m^2 = E^2 - p^2 $ (the Lorentz self-invariant).
double rapidity () const
 Calculate the rapidity.
double rap () const
 Alias for rapidity.
double absrapidity () const
 Absolute rapidity.
double absrap () const
 Absolute rapidity.
double pT2 () const
 Calculate the squared transverse momentum $ p_T^2 $.
double pt2 () const
 Calculate the squared transverse momentum $ p_T^2 $.
double pT () const
 Calculate the transverse momentum $ p_T $.
double pt () const
 Calculate the transverse momentum $ p_T $.
double Et2 () const
 Calculate the transverse energy $ E_T^2 = E^2 \sin^2{\theta} $.
double Et () const
 Calculate the transverse energy $ E_T = E \sin{\theta} $.
Vector3 boostVector () const
 Calculate the boost vector (in units of $ \beta $).
FourMomentum & operator*= (double a)
 Multiply by a scalar.
FourMomentum & operator/= (double a)
 Divide by a scalar.
FourMomentum & operator+= (const FourMomentum &v)
 Add to this 4-vector. NB time as well as space components are added.
FourMomentum & operator-= (const FourMomentum &v)
 Subtract from this 4-vector. NB time as well as space components are subtracted.
FourMomentum operator- () const
 Multiply all components (time and space) by -1.
double t () const
double x () const
double y () const
double z () const
FourVector & setT (const double t)
FourVector & setX (const double x)
FourVector & setY (const double y)
FourVector & setZ (const double z)
double invariant () const
double angle (const FourVector &v) const
 Angle between this vector and another.
double angle (const Vector3 &v3) const
 Angle between this vector and another (3-vector)
double polarRadius2 () const
 Square of the projection of the 3-vector on to the $ x-y $ plane This is a more efficient function than polarRadius, as it avoids the square root. Use it if you only need the squared value, or e.g. an ordering by magnitude.
double perp2 () const
 Synonym for polarRadius2.
double rho2 () const
 Synonym for polarRadius2.
double polarRadius () const
 Projection of 3-vector on to the $ x-y $ plane.
double perp () const
 Synonym for polarRadius.
double rho () const
 Synonym for polarRadius.
double azimuthalAngle (const PhiMapping mapping=ZERO_2PI) const
 Angle subtended by the 3-vector's projection in x-y and the x-axis.
double phi (const PhiMapping mapping=ZERO_2PI) const
 Synonym for azimuthalAngle.
double polarAngle () const
 Angle subtended by the 3-vector and the z-axis.
double theta () const
 Synonym for polarAngle.
double pseudorapidity () const
 Pseudorapidity (defined purely by the 3-vector components)
double eta () const
 Synonym for pseudorapidity.
double abspseudorapidity () const
 Get the $ |\eta| $ directly.
double abseta () const
 Get the $ |\eta| $ directly (alias).
Vector3 vector3 () const
 Get the spatial part of the 4-vector as a 3-vector.
double contract (const FourVector &v) const
 Contract two 4-vectors, with metric signature (+ - - -).
double dot (const FourVector &v) const
 Contract two 4-vectors, with metric signature (+ - - -).
double operator* (const FourVector &v) const
 Contract two 4-vectors, with metric signature (+ - - -).
FourVector & operator+= (const FourVector &v)
 Add to this 4-vector.
FourVector & operator-= (const FourVector &v)
 Subtract from this 4-vector. NB time as well as space components are subtracted.
const double & get (const size_t index) const
const double & operator[] (const size_t index) const
 Direct access to vector elements by index.
double & operator[] (const size_t index)
 Direct access to vector elements by index.
Vector< N > & set (const size_t index, const double value)
 Set indexed value.
size_t size () const
 Vector dimensionality.
bool isZero (double tolerance=1E-5) const
 Check for nullness, allowing for numerical precision.
double mod2 () const
 Calculate the modulus-squared of a vector. $ \sum_{i=1}^N x_i^2 $.
double mod () const
 Calculate the modulus of a vector. $ \sqrt{\sum_{i=1}^N x_i^2} $.
bool operator== (const Vector< N > &a) const
bool operator!= (const Vector< N > &a) const
bool operator< (const Vector< N > &a) const
bool operator<= (const Vector< N > &a) const
bool operator> (const Vector< N > &a) const
bool operator>= (const Vector< N > &a) const

Protected Member Functions

double & get (const size_t index)

Protected Attributes

Eigen::Vector< double, N > _vec
 Vector.

Friends

FourMomentum multiply (const double a, const FourMomentum &v)
FourMomentum multiply (const FourMomentum &v, const double a)
FourMomentum add (const FourMomentum &a, const FourMomentum &b)
FourMomentum transform (const LorentzTransform &lt, const FourMomentum &v4)

Detailed Description

Specialized version of the FourVector with momentum/energy functionality.

Definition at line 274 of file Vector4.hh.

Constructor & Destructor Documentation

FourMomentum ( ) [inline]

Definition at line 281 of file Vector4.hh.

{ }
FourMomentum ( const V4 &  other) [inline]

Definition at line 284 of file Vector4.hh.

References FourMomentum::setE(), FourMomentum::setPx(), FourMomentum::setPy(), and FourMomentum::setPz().

                                  {
      this->setE(other.t());
      this->setPx(other.x());
      this->setPy(other.y());
      this->setPz(other.z());
    }
FourMomentum ( const Vector< 4 > &  other) [inline]

Definition at line 291 of file Vector4.hh.

      : FourVector(other) { }
FourMomentum ( const double  E,
const double  px,
const double  py,
const double  pz 
) [inline]

Definition at line 294 of file Vector4.hh.

References FourMomentum::setE(), FourMomentum::setPx(), FourMomentum::setPy(), and FourMomentum::setPz().

                                                                                    {
      this->setE(E);
      this->setPx(px);
      this->setPy(py);
      this->setPz(pz);
    }
~FourMomentum ( ) [inline]

Definition at line 301 of file Vector4.hh.

{}

Member Function Documentation

double abseta ( ) const [inline, inherited]

Get the $ |\eta| $ directly (alias).

Definition at line 142 of file Vector4.hh.

References FourVector::eta().

Referenced by ParticleBase::abseta(), CMS_2013_I1265659::analyze(), CDF_2001_S4517016::analyze(), CMS_2014_I1303894::analyze(), and ATLAS_2012_I1204447::analyze().

{ return fabs(eta()); }
double abspseudorapidity ( ) const [inline, inherited]

Get the $ |\eta| $ directly.

Definition at line 139 of file Vector4.hh.

References FourVector::eta().

Referenced by ParticleBase::abspseudorapidity().

{ return fabs(eta()); }
double absrap ( ) const [inline]

Absolute rapidity.

Definition at line 367 of file Vector4.hh.

References FourMomentum::rap().

Referenced by ParticleBase::absrap(), CDF_2008_S8093652::analyze(), D0_2010_S8566488::analyze(), D0_2008_S7863608::analyze(), ATLAS_2013_I1230812::analyze(), and ATLAS_2011_I945498::selectJets().

                          {
      return fabs(rap());
    }
double absrapidity ( ) const [inline]

Absolute rapidity.

Definition at line 362 of file Vector4.hh.

References FourMomentum::rapidity().

Referenced by ParticleBase::absrapidity().

                               {
      return fabs(rapidity());
    }
double angle ( const FourVector &  v) const [inline, inherited]

Angle between this vector and another.

Definition at line 67 of file Vector4.hh.

References Vector3::angle(), and FourVector::vector3().

Referenced by D0_1996_S3214044::_fourJetAnalysis(), H1_2000_S4129130::analyze(), H1_1994_S2919893::analyze(), ParticleBase::angle(), Rivet::angle(), and MC_VH2BB::boostAngles().

                                            {
      return vector3().angle( v.vector3() );
    }
double angle ( const Vector3 &  v3) const [inline, inherited]

Angle between this vector and another (3-vector)

Definition at line 72 of file Vector4.hh.

References Vector3::angle(), and FourVector::vector3().

                                          {
      return vector3().angle(v3);
    }
double azimuthalAngle ( const PhiMapping  mapping = ZERO_2PI) const [inline, inherited]

Angle subtended by the 3-vector's projection in x-y and the x-axis.

Definition at line 109 of file Vector4.hh.

References Vector3::azimuthalAngle(), and FourVector::vector3().

Referenced by CDF_2010_S8591881_DY::analyze(), ParticleBase::azimuthalAngle(), Rivet::deltaR(), and DISKinematics::project().

                                                                   {
      return vector3().azimuthalAngle(mapping);
    }
Vector3 boostVector ( ) const [inline]

Calculate the boost vector (in units of $ \beta $).

Todo:
Be careful about c=1 convention...

Definition at line 402 of file Vector4.hh.

References FourMomentum::E(), FourMomentum::px(), FourMomentum::py(), and FourMomentum::pz().

Referenced by CDF_1997_S3541940::_avg_beam_in_lab(), CDF_1996_S3349578::_avg_beam_in_lab(), CDF_1996_S3349578::_fiveJetAnalysis(), D0_1996_S3214044::_fourJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), CDF_1996_S3349578::_threeJetAnalysis(), D0_1996_S3214044::_threeJetAnalysis(), CLEO_2004_S5809304::analyze(), BELLE_2001_S4598261::analyze(), BABAR_2005_S6181155::analyze(), BABAR_2007_S6895344::analyze(), BABAR_2003_I593379::analyze(), BABAR_2013_I1238276::analyze(), ARGUS_1993_S2653028::analyze(), BELLE_2013_I1216515::analyze(), BELLE_2008_I786560::analyze(), ARGUS_1993_S2669951::analyze(), BABAR_2007_S7266081::analyze(), CDF_1996_S3108457::analyze(), CDF_1997_S3541940::analyze(), ARGUS_1993_S2789213::analyze(), PDG_TAUS::analyzeRadiativeDecay(), MC_VH2BB::boostAngles(), and DISKinematics::project().

                                {
      // const Vector3 p3 = vector3();
      // const double m2 = mass2();
      // if (Rivet::isZero(m2)) return p3.unit();
      // else {
      //   // Could also do this via beta = tanh(rapidity), but that's
      //   // probably more messy from a numerical hygiene point of view.
      //   const double p2 = p3.mod2();
      //   const double beta = sqrt( p2 / (m2 + p2) );
      //   return beta * p3.unit();
      // }
      /// @todo Be careful about c=1 convention...
      return Vector3(px()/E(), py()/E(), pz()/E());
    }
double contract ( const FourVector &  v) const [inline, inherited]

Contract two 4-vectors, with metric signature (+ - - -).

Definition at line 153 of file Vector4.hh.

References FourVector::t(), FourVector::x(), FourVector::y(), and FourVector::z().

Referenced by Rivet::contract(), FourVector::dot(), and FourVector::operator*().

                                               {
      const double result = t()*v.t() - x()*v.x() - y()*v.y() - z()*v.z();
      return result;
    }
double dot ( const FourVector &  v) const [inline, inherited]

Contract two 4-vectors, with metric signature (+ - - -).

Definition at line 159 of file Vector4.hh.

References FourVector::contract().

                                          {
      return contract(v);
    }
double E ( ) const [inline]

Get energy $ E $ (time component of momentum).

Definition at line 305 of file Vector4.hh.

References FourVector::t().

Referenced by CDF_1997_S3541940::_avg_beam_in_lab(), CDF_1996_S3349578::_avg_beam_in_lab(), BeamThrust::_calcBeamThrust(), CDF_1996_S3349578::_fiveJetAnalysis(), D0_1996_S3214044::_fourJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), CDF_1996_S3349578::_threeJetAnalysis(), D0_1996_S3214044::_threeJetAnalysis(), OPAL_1993_S2692198::analyze(), H1_2000_S4129130::analyze(), ARGUS_1993_S2669951::analyze(), ALICE_2012_I1181770::analyze(), H1_1994_S2919893::analyze(), D0_2006_S6438750::analyze(), MC_PHOTONINC::analyze(), MC_PHOTONKTSPLITTINGS::analyze(), CDF_1997_S3541940::analyze(), MC_HFJETS::analyze(), MC_PHOTONJETS::analyze(), LHCB_2010_S8758301::analyze(), LHCB_2011_I917009::analyze(), ARGUS_1993_S2789213::analyze(), D0_2008_S7719523::analyze(), ATLAS_2011_S9041966::analyze(), ALEPH_1996_S3486095::analyze(), FourMomentum::boostVector(), Hemispheres::calc(), FastJets::calc(), Rivet::cmpMomByAscE(), Rivet::cmpMomByE(), ATLAS_2014_I1306615::cosTS_CS(), ParticleBase::E(), ParticleBase::energy(), FourMomentum::Et(), FourMomentum::byEAscending::operator()(), FourMomentum::rapidity(), Jet::setState(), Rivet::sqrtS(), and Jet::totalEnergy().

{ return t(); }
double Et ( ) const [inline]

Calculate the transverse energy $ E_T = E \sin{\theta} $.

Definition at line 397 of file Vector4.hh.

References FourMomentum::E(), and FourVector::polarAngle().

Referenced by H1_2000_S4129130::analyze(), H1_1994_S2919893::analyze(), CDF_2005_S6080774::analyze(), D0_1996_S3324664::analyze(), MC_DIPHOTON::analyze(), CDF_2009_S8436959::analyze(), CDF_2001_S4517016::analyze(), D0_2008_S7837160::analyze(), MC_WINC::analyze(), ATLAS_2012_I946427::analyze(), ATLAS_2010_S8919674::analyze(), ATLAS_2014_I1307756::analyze(), CDF_2008_S7541902::analyze(), ATLAS_2012_I1199269::analyze(), ATLAS_2011_I921594::analyze(), ATLAS_2013_I1263495::analyze(), ATLAS_2011_S9120807::analyze(), ATLAS_2010_S8914702::analyze(), ATLAS_2012_I1083318::analyze(), ATLAS_2014_I1306615::analyze(), ATLAS_2012_I1093738::analyze(), Rivet::cmpMomByAscEt(), Rivet::cmpMomByEt(), ParticleBase::Et(), FourMomentum::Et2(), ATLAS_2011_I928289_W::FillPlots1d(), InvMassFinalState::massT(), and MissingMomentum::project().

                      {
      return E() * sin(polarAngle());
    }
double Et2 ( ) const [inline]

Calculate the transverse energy $ E_T^2 = E^2 \sin^2{\theta} $.

Definition at line 392 of file Vector4.hh.

References FourMomentum::Et().

                       {
      return Et() * Et();
    }
double eta ( ) const [inline, inherited]

Synonym for pseudorapidity.

Definition at line 134 of file Vector4.hh.

References Vector3::eta(), and FourVector::vector3().

Referenced by CDF_2004_S5839831::_calcTransCones(), FourVector::abseta(), FourVector::abspseudorapidity(), H1_2000_S4129130::analyze(), CMS_2013_I1265659::analyze(), CMS_2011_S8941262::analyze(), H1_1994_S2919893::analyze(), ATLAS_2011_S9002537::analyze(), CDF_2000_S4266730::analyze(), D0_1996_S3324664::analyze(), MC_ZJETS::analyze(), D0_2008_S6879055::analyze(), CDF_1993_S2742446::analyze(), CDF_2009_S8436959::analyze(), CDF_2008_S7540469::analyze(), D0_2010_S8821313::analyze(), D0_2006_S6438750::analyze(), MC_WINC::analyze(), CDF_1996_S3418421::analyze(), D0_2008_S7837160::analyze(), CDF_1994_S2952106::analyze(), ATLAS_2014_I1298811::analyze(), MC_PHOTONJETS::analyze(), MC_ZZJETS::analyze(), MC_WWJETS::analyze(), ATLAS_2013_I1243871::analyze(), LHCB_2014_I1281685::analyze(), CDF_2008_S7541902::analyze(), ATLAS_2012_I1199269::analyze(), MC_PHOTONJETUE::analyze(), MC_ZZINC::analyze(), MC_WWINC::analyze(), CDF_2010_S8591881_DY::analyze(), LHCB_2012_I1119400::analyze(), MC_TTBAR::analyze(), ATLAS_2010_S8914702::analyze(), D0_2008_S7719523::analyze(), MC_SUSY::analyze(), ATLAS_2014_I1306615::analyze(), ATLAS_2012_I1093738::analyze(), MC_VH2BB::analyze(), CDF_2004_S5839831::analyze(), JetShape::calc(), Rivet::cmpMomByAbsEta(), Rivet::cmpMomByDescAbsEta(), Rivet::cmpMomByEta(), ParticleBase::eta(), ATLAS_2012_I1204784::fillPlots(), and ParticleBase::pseudorapidity().

                       {
      return vector3().eta();
    }
const double& get ( const size_t  index) const [inline, inherited]

Definition at line 33 of file VectorN.hh.

References Vector< N >::_vec.

                                                {
      if (index >= N) {
        throw std::runtime_error("Tried to access an invalid vector index.");
      } else {
        return _vec(index);
      }
    }
double& get ( const size_t  index) [inline, protected, inherited]

Definition at line 126 of file VectorN.hh.

References Vector< N >::_vec.

                                    {
      if (index >= N) {
        throw std::runtime_error("Tried to access an invalid vector index.");
      } else {
        return _vec(index);
      }
    }
double invariant ( ) const [inline, inherited]

Definition at line 61 of file Vector4.hh.

References FourVector::t(), FourVector::x(), FourVector::y(), and FourVector::z().

Referenced by Rivet::invariant(), and FourMomentum::mass2().

                             {
      // Done this way for numerical precision
      return (t() + z())*(t() - z()) - x()*x() - y()*y();
    }
bool isZero ( double  tolerance = 1E-5) const [inline, inherited]

Check for nullness, allowing for numerical precision.

Definition at line 67 of file VectorN.hh.

References Vector< N >::_vec, and Rivet::isZero().

                                             {
      for (size_t i=0; i < N; ++i) {
        if (! Rivet::isZero(_vec[i], tolerance) ) return false;
      }
      return true;
    }
double mass ( ) const [inline]

Get the mass $ m = \sqrt{E^2 - p^2} $ (the Lorentz self-invariant).

For spacelike momenta, the mass will be -sqrt(|mass2|).

Definition at line 336 of file Vector4.hh.

References FourMomentum::mass2(), and Rivet::sign().

Referenced by OPAL_1993_S2692198::analyze(), CMS_2013_I1265659::analyze(), CDF_2008_S8093652::analyze(), CMS_2011_S8968497::analyze(), CMS_2012_I1090423::analyze(), CDF_2000_S4155203::analyze(), ATLAS_2014_I1282441::analyze(), CDF_2000_S4266730::analyze(), CDF_2005_S6080774::analyze(), EXAMPLE_CUTS::analyze(), MC_DIPHOTON::analyze(), ATLAS_2011_I930220::analyze(), CDF_2008_S7540469::analyze(), D0_2010_S8570965::analyze(), CDF_1996_S3108457::analyze(), D0_2010_S8566488::analyze(), CDF_1996_S3418421::analyze(), D0_2009_S8320160::analyze(), ATLAS_2012_I1126136::analyze(), ATLAS_2014_I1307756::analyze(), ATLAS_2012_I1199269::analyze(), MC_ZZINC::analyze(), CDF_2010_S8591881_DY::analyze(), ATLAS_2011_S9120807::analyze(), MC_TTBAR::analyze(), ATLAS_2011_S9019561::analyze(), ATLAS_2012_CONF_2012_153::analyze(), MC_SUSY::analyze(), MC_VH2BB::analyze(), ATLAS_2012_I1203852::analyze(), ATLAS_2014_I1279489::analyze(), InvMassFinalState::calc(), Rivet::identifyZstates(), ParticleBase::mass(), and ATLAS_2014_I1306615::tau_jet().

                        {
      // assert(Rivet::isZero(mass2()) || mass2() > 0);
      // if (Rivet::isZero(mass2())) {
      //   return 0.0;
      // } else {
      //   return sqrt(mass2());
      // }
      return sign(mass2()) * sqrt(fabs(mass2()));
    }
double mass2 ( ) const [inline]

Get the squared mass $ m^2 = E^2 - p^2 $ (the Lorentz self-invariant).

Definition at line 347 of file Vector4.hh.

References FourVector::invariant().

Referenced by CDF_1997_S3541940::_reduce(), CDF_1996_S3349578::_reduce(), D0_1996_S3214044::_safeMass(), CDF_1997_S3541940::_safeMass(), CDF_1996_S3349578::_safeMass(), ATLAS_2011_I894867::analyze(), CMS_2012_I1193338::analyze(), CDF_2000_S4155203::analyze(), CMS_2013_I1218372::analyze(), D0_2001_S4674421::analyze(), MC_ZZINC::analyze(), MC_WWINC::analyze(), ATLAS_2013_I1217863_W::analyze(), ATLAS_2012_I1083318::analyze(), InvMassFinalState::calc(), FourMomentum::mass(), ParticleBase::mass2(), and DISKinematics::project().

                         {
      return invariant();
    }
double mod ( ) const [inline, inherited]

Calculate the modulus of a vector. $ \sqrt{\sum_{i=1}^N x_i^2} $.

Definition at line 87 of file VectorN.hh.

References Vector< N >::mod2().

                       {
      const double norm = mod2();
      assert(norm >= 0);
      return sqrt(norm);
    }
double mod2 ( ) const [inline, inherited]

Calculate the modulus-squared of a vector. $ \sum_{i=1}^N x_i^2 $.

Definition at line 76 of file VectorN.hh.

References Vector< N >::mod2(), and Vector< N >::size().

                        {
      double mod2 = 0.0;
      for (size_t i = 0; i < size(); ++i) {
        const double element = get(i);
        mod2 += element*element;
      }
      return mod2;
    }
bool operator!= ( const Vector< N > &  a) const [inline, inherited]

Definition at line 104 of file VectorN.hh.

References Vector< N >::_vec.

                                              {
      return _vec != a._vec;
    }
double operator* ( const FourVector &  v) const [inline, inherited]

Contract two 4-vectors, with metric signature (+ - - -).

Definition at line 164 of file Vector4.hh.

References FourVector::contract().

                                                {
      return contract(v);
    }
FourMomentum& operator*= ( double  a) [inline]

Multiply by a scalar.

Reimplemented from FourVector.

Definition at line 443 of file Vector4.hh.

References Vector< 4 >::_vec, Vector< N >::_vec, and FourMomentum::multiply.

                                       {
      _vec = multiply(a, *this)._vec;
      return *this;
    }
FourVector& operator+= ( const FourVector &  v) [inline, inherited]

Add to this 4-vector.

Definition at line 181 of file Vector4.hh.

References Vector< 4 >::_vec, Vector< N >::_vec, and FourVector::add.

                                                {
      _vec = add(*this, v)._vec;
      return *this;
    }
FourMomentum& operator+= ( const FourMomentum &  v) [inline]

Add to this 4-vector. NB time as well as space components are added.

Definition at line 455 of file Vector4.hh.

References Vector< 4 >::_vec, Vector< N >::_vec, and FourMomentum::add.

                                                    {
      _vec = add(*this, v)._vec;
      return *this;
    }
FourMomentum operator- ( ) const [inline]

Multiply all components (time and space) by -1.

Reimplemented from FourVector.

Definition at line 467 of file Vector4.hh.

References Vector< N >::_vec, and Vector< 4 >::_vec.

                                   {
      FourMomentum result;
      result._vec = -_vec;
      return result;
    }
FourVector& operator-= ( const FourVector &  v) [inline, inherited]

Subtract from this 4-vector. NB time as well as space components are subtracted.

Definition at line 187 of file Vector4.hh.

References Vector< 4 >::_vec, Vector< N >::_vec, and FourVector::add.

                                                {
      _vec = add(*this, -v)._vec;
      return *this;
    }
FourMomentum& operator-= ( const FourMomentum &  v) [inline]

Subtract from this 4-vector. NB time as well as space components are subtracted.

Definition at line 461 of file Vector4.hh.

References Vector< 4 >::_vec, Vector< N >::_vec, and FourMomentum::add.

                                                    {
      _vec = add(*this, -v)._vec;
      return *this;
    }
FourMomentum& operator/= ( double  a) [inline]

Divide by a scalar.

Reimplemented from FourVector.

Definition at line 449 of file Vector4.hh.

References Vector< 4 >::_vec, Vector< N >::_vec, and FourMomentum::multiply.

                                       {
      _vec = multiply(1.0/a, *this)._vec;
      return *this;
    }
bool operator< ( const Vector< N > &  a) const [inline, inherited]

Definition at line 108 of file VectorN.hh.

References Vector< N >::_vec.

                                             {
      return _vec < a._vec;
    }
bool operator<= ( const Vector< N > &  a) const [inline, inherited]

Definition at line 112 of file VectorN.hh.

References Vector< N >::_vec.

                                              {
      return _vec <= a._vec;
    }
bool operator== ( const Vector< N > &  a) const [inline, inherited]

Definition at line 100 of file VectorN.hh.

References Vector< N >::_vec.

                                              {
      return _vec == a._vec;
    }
bool operator> ( const Vector< N > &  a) const [inline, inherited]

Definition at line 116 of file VectorN.hh.

References Vector< N >::_vec.

                                             {
      return _vec > a._vec;
    }
bool operator>= ( const Vector< N > &  a) const [inline, inherited]

Definition at line 120 of file VectorN.hh.

References Vector< N >::_vec.

                                              {
      return _vec >= a._vec;
    }
const double& operator[] ( const size_t  index) const [inline, inherited]

Direct access to vector elements by index.

Definition at line 42 of file VectorN.hh.

                                                       {
      return get(index);
    }
double& operator[] ( const size_t  index) [inline, inherited]

Direct access to vector elements by index.

Definition at line 47 of file VectorN.hh.

                                           {
      return get(index);
    }
Vector3 p3 ( ) const [inline]

Get 3-momentum part, $ p $.

Definition at line 308 of file Vector4.hh.

References FourVector::vector3().

Referenced by CDF_1996_S3349578::_fiveJetAnalysis(), D0_1996_S3214044::_fourJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), CDF_1997_S3541940::_psi(), CDF_1996_S3349578::_psi(), CDF_1996_S3349578::_threeJetAnalysis(), D0_1996_S3214044::_threeJetAnalysis(), CLEO_2004_S5809304::analyze(), OPAL_1993_S2692198::analyze(), ARGUS_1993_S2669951::analyze(), ALICE_2012_I1181770::analyze(), CDF_1997_S3541940::analyze(), ARGUS_1993_S2789213::analyze(), LHCB_2014_I1281685::analyze(), MC_ZZINC::analyze(), MC_WWINC::analyze(), LHCB_2012_I1119400::analyze(), and ALEPH_1996_S3486095::analyze().

{ return vector3(); }
double perp ( ) const [inline, inherited]

Synonym for polarRadius.

Definition at line 99 of file Vector4.hh.

References Vector3::perp(), and FourVector::vector3().

                        {
      return vector3().perp();
    }
double perp2 ( ) const [inline, inherited]

Synonym for polarRadius2.

Definition at line 84 of file Vector4.hh.

References Vector3::perp2(), and FourVector::vector3().

                         {
      return vector3().perp2();
    }
double phi ( const PhiMapping  mapping = ZERO_2PI) const [inline, inherited]

Synonym for azimuthalAngle.

Definition at line 114 of file Vector4.hh.

References Vector3::phi(), and FourVector::vector3().

Referenced by CDF_2004_S5839831::_calcTransCones(), CMS_2013_I1265659::analyze(), ATLAS_2011_S9002537::analyze(), D0_1996_S3324664::analyze(), D0_2008_S6879055::analyze(), CDF_1993_S2742446::analyze(), CDF_2009_S8436959::analyze(), CDF_2008_S7540469::analyze(), ATLAS_2011_I930220::analyze(), D0_2010_S8821313::analyze(), D0_2006_S6438750::analyze(), CDF_1994_S2952106::analyze(), ATLAS_2014_I1298811::analyze(), MC_PHOTONJETS::analyze(), CMS_2011_S9120041::analyze(), ATLAS_2012_I1199269::analyze(), MC_PHOTONJETUE::analyze(), ATLAS_2011_S9212183::analyze(), MC_WWINC::analyze(), MC_ZZINC::analyze(), ATLAS_2011_CONF_2011_098::analyze(), ATLAS_2014_I1315949::analyze(), ATLAS_2011_S8983313::analyze(), ATLAS_2012_I1095236::analyze(), ATLAS_2012_I1125961::analyze(), ATLAS_2013_I1217867::analyze(), ATLAS_2011_CONF_2011_090::analyze(), ATLAS_2012_CONF_2012_109::analyze(), ATLAS_2010_S8914702::analyze(), D0_2008_S7719523::analyze(), ATLAS_2012_I1083318::analyze(), ATLAS_2014_I1304688::analyze(), ATLAS_2011_S9212353::analyze(), ATLAS_2011_S9041966::analyze(), MC_SUSY::analyze(), ATLAS_2012_I1093738::analyze(), MC_VH2BB::analyze(), CDF_2004_S5839831::analyze(), and ParticleBase::phi().

                                                        {
      return vector3().phi(mapping);
    }
double polarAngle ( ) const [inline, inherited]

Angle subtended by the 3-vector and the z-axis.

Definition at line 119 of file Vector4.hh.

References Vector3::polarAngle(), and FourVector::vector3().

Referenced by H1_1994_S2919893::beamAngle(), FourMomentum::Et(), ParticleBase::polarAngle(), and DISKinematics::project().

                              {
      return vector3().polarAngle();
    }
double polarRadius ( ) const [inline, inherited]

Projection of 3-vector on to the $ x-y $ plane.

Definition at line 94 of file Vector4.hh.

References Vector3::polarRadius(), and FourVector::vector3().

                               {
      return vector3().polarRadius();
    }
double polarRadius2 ( ) const [inline, inherited]

Square of the projection of the 3-vector on to the $ x-y $ plane This is a more efficient function than polarRadius, as it avoids the square root. Use it if you only need the squared value, or e.g. an ordering by magnitude.

Definition at line 79 of file Vector4.hh.

References Vector3::polarRadius2(), and FourVector::vector3().

                                {
      return vector3().polarRadius2();
    }
double pseudorapidity ( ) const [inline, inherited]

Pseudorapidity (defined purely by the 3-vector components)

Definition at line 129 of file Vector4.hh.

References Vector3::pseudorapidity(), and FourVector::vector3().

Referenced by Rivet::cmpMomByDescEta().

                                  {
      return vector3().pseudorapidity();
    }
double pT ( ) const [inline]

Calculate the transverse momentum $ p_T $.

Definition at line 382 of file Vector4.hh.

References FourMomentum::pT2().

Referenced by CDF_2004_S5839831::_calcTransCones(), CMS_2013_I1265659::analyze(), CMS_2011_S8941262::analyze(), CMS_2012_PAS_QCD_11_010::analyze(), CDF_2000_S4155203::analyze(), STAR_2006_S6870392::analyze(), H1_1994_S2919893::analyze(), ATLAS_2011_S9002537::analyze(), EXAMPLE_CUTS::analyze(), CDF_2005_S6080774::analyze(), CMS_2013_I1273574::analyze(), MC_DIPHOTON::analyze(), ATLAS_2011_I925932::analyze(), D0_2010_S8570965::analyze(), CDF_2008_S7782535::analyze(), D0_2006_S6438750::analyze(), MC_PHOTONINC::analyze(), MC_HFJETS::analyze(), ATLAS_2014_I1298811::analyze(), D0_2008_S7863608::analyze(), ATLAS_2011_I954993::analyze(), ATLAS_2011_S9131140::analyze(), MC_ZZJETS::analyze(), CDF_2008_S8095620::analyze(), CMS_2011_S9120041::analyze(), MC_WWJETS::analyze(), CMS_2013_I1256943::analyze(), ATLAS_2012_I946427::analyze(), D0_2009_S8349509::analyze(), ATLAS_2012_CONF_2012_104::analyze(), ATLAS_2012_I1126136::analyze(), CMS_2012_I941555::analyze(), ATLAS_2013_I1243871::analyze(), LHCB_2014_I1281685::analyze(), ATLAS_2012_I1190891::analyze(), D0_2001_S4674421::analyze(), CDF_2008_S7541902::analyze(), ATLAS_2012_CONF_2012_105::analyze(), ATLAS_2012_I1199269::analyze(), MC_PHOTONJETUE::analyze(), MC_WWINC::analyze(), ATLAS_2011_S9212183::analyze(), MC_ZZINC::analyze(), ATLAS_2011_CONF_2011_098::analyze(), ATLAS_2014_I1315949::analyze(), CDF_2010_S8591881_DY::analyze(), ATLAS_2012_CONF_2012_103::analyze(), ATLAS_2012_I1112263::analyze(), ATLAS_2011_S8983313::analyze(), ATLAS_2013_I1190187::analyze(), LHCB_2012_I1119400::analyze(), ATLAS_2012_I1095236::analyze(), ATLAS_2011_S9120807::analyze(), ATLAS_2012_I1125961::analyze(), ATLAS_2013_I1217867::analyze(), ATLAS_2012_I1180197::analyze(), ATLAS_2011_CONF_2011_090::analyze(), ATLAS_2012_CONF_2012_001::analyze(), ATLAS_2014_I1307243::analyze(), ATLAS_2012_CONF_2012_109::analyze(), ATLAS_2012_I1117704::analyze(), D0_2008_S7719523::analyze(), CMS_2014_I1303894::analyze(), ATLAS_2011_S9019561::analyze(), ATLAS_2012_I1204447::analyze(), ATLAS_2012_CONF_2012_153::analyze(), ATLAS_2012_I1083318::analyze(), ATLAS_2013_I1230812::analyze(), ATLAS_2014_I1304688::analyze(), ATLAS_2011_S9212353::analyze(), ATLAS_2011_S9041966::analyze(), MC_SUSY::analyze(), ATLAS_2012_I943401::analyze(), ATLAS_2014_I1306615::analyze(), ATLAS_2011_S9225137::analyze(), ATLAS_2012_I1093738::analyze(), MC_VH2BB::analyze(), ATLAS_2012_I1094568::analyze(), CDF_2004_S5839831::analyze(), ATLAS_2012_I1203852::analyze(), JetShape::calc(), Zstate::cmppT(), ATLAS_2014_I1288706::FillPlots(), ATLAS_2014_I1300647::FillPlots1d(), ATLAS_2011_I928289_W::FillPlots1d(), ATLAS_2014_I1300647::FillPlots3d(), WFinder::mT(), ATLAS_2014_I1306615::tau_jet(), and Variables::Variables().

                      {
      return sqrt(pT2());
    }
double pt ( ) const [inline]

Calculate the transverse momentum $ p_T $.

Definition at line 387 of file Vector4.hh.

References FourMomentum::pT2().

Referenced by Rivet::cmpMomByAscPt(), Rivet::cmpMomByPt(), and ParticleBase::pt().

                      {
      return sqrt(pT2());
    }
double pT2 ( ) const [inline]

Calculate the squared transverse momentum $ p_T^2 $.

Definition at line 372 of file Vector4.hh.

References Vector3::polarRadius2(), and FourVector::vector3().

Referenced by H1_1994_S2919893::analyze(), ATLAS_2013_I1217863_W::analyze(), FourMomentum::pT(), and FourMomentum::pt().

                       {
      return vector3().polarRadius2();
    }
double pt2 ( ) const [inline]

Calculate the squared transverse momentum $ p_T^2 $.

Definition at line 377 of file Vector4.hh.

References Vector3::polarRadius2(), and FourVector::vector3().

Referenced by ParticleBase::pt2().

                       {
      return vector3().polarRadius2();
    }
double px ( ) const [inline]

Get 3-momentum part, $ p $ (alias).

Get x-component of momentum $ p_x $.

Definition at line 313 of file Vector4.hh.

References FourVector::x().

Referenced by OPAL_1993_S2692198::analyze(), CMS_2013_I1273574::analyze(), LHCB_2010_S8758301::analyze(), LHCB_2011_I917009::analyze(), ATLAS_2012_CONF_2012_104::analyze(), ATLAS_2012_I1126136::analyze(), CDF_2008_S7541902::analyze(), ATLAS_2012_I1095236::analyze(), ATLAS_2012_I1180197::analyze(), ATLAS_2012_I1203852::analyze(), FourMomentum::boostVector(), FastJets::calc(), DISKinematics::project(), ATLAS_2014_I1306615::pTt(), and Jet::setState().

{ return x(); }
double py ( ) const [inline]

Get y-component of momentum $ p_y $.

Definition at line 316 of file Vector4.hh.

References FourVector::y().

Referenced by OPAL_1993_S2692198::analyze(), CMS_2013_I1273574::analyze(), LHCB_2010_S8758301::analyze(), LHCB_2011_I917009::analyze(), ATLAS_2012_CONF_2012_104::analyze(), ATLAS_2012_I1126136::analyze(), CDF_2008_S7541902::analyze(), ATLAS_2012_I1095236::analyze(), ATLAS_2012_I1180197::analyze(), ATLAS_2012_I1203852::analyze(), FourMomentum::boostVector(), FastJets::calc(), ATLAS_2014_I1306615::pTt(), and Jet::setState().

{ return y(); }
double pz ( ) const [inline]

Get z-component of momentum $ p_z $.

Definition at line 319 of file Vector4.hh.

References FourVector::z().

Referenced by OPAL_1993_S2692198::analyze(), LHCB_2010_S8758301::analyze(), LHCB_2011_I917009::analyze(), FourMomentum::boostVector(), FastJets::calc(), ATLAS_2014_I1306615::cosTS_CS(), DISLepton::project(), FourMomentum::rapidity(), Jet::setState(), and Rivet::sqrtS().

{ return z(); }
double rap ( ) const [inline]

Alias for rapidity.

Definition at line 357 of file Vector4.hh.

References FourMomentum::rapidity().

Referenced by FourMomentum::absrap().

                       {
      return rapidity();
    }
double rapidity ( ) const [inline]

Calculate the rapidity.

Definition at line 352 of file Vector4.hh.

References FourMomentum::E(), and FourMomentum::pz().

Referenced by CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), CDF_1996_S3349578::_threeJetAnalysis(), FourMomentum::absrapidity(), LHCB_2012_I1208102::analyze(), CMS_2011_S8968497::analyze(), ATLAS_2011_I930220::analyze(), D0_2010_S8821313::analyze(), MC_PHOTONINC::analyze(), CDF_1997_S3541940::analyze(), D0_2009_S8320160::analyze(), CMS_2012_I941555::analyze(), D0_2009_S8349509::analyze(), CDF_2008_S7541902::analyze(), ATLAS_2014_I1307243::analyze(), D0_2008_S7719523::analyze(), ATLAS_2012_I1083318::analyze(), ATLAS_2014_I1306615::analyze(), ATLAS_2012_I1093738::analyze(), JetShape::calc(), Rivet::cmpMomByAbsRap(), Rivet::cmpMomByDescAbsRap(), Rivet::cmpMomByDescRap(), Rivet::cmpMomByRap(), Rivet::deltaR(), Rivet::deltaRap(), ATLAS_2011_I928289_Z::FillPlots1d(), ATLAS_2014_I1300647::FillPlots3d(), ParticleBase::rap(), FourMomentum::rap(), ParticleBase::rapidity(), ATLAS_2014_I1306615::tau_jet(), and Variables::Variables().

                            {
      return 0.5 * std::log( (E() + pz()) / (E() - pz()) );
    }
double rho ( ) const [inline, inherited]

Synonym for polarRadius.

Definition at line 104 of file Vector4.hh.

References Vector3::rho(), and FourVector::vector3().

                       {
      return vector3().rho();
    }
double rho2 ( ) const [inline, inherited]

Synonym for polarRadius2.

Definition at line 89 of file Vector4.hh.

References Vector3::rho2(), and FourVector::vector3().

                        {
      return vector3().rho2();
    }
Vector<N>& set ( const size_t  index,
const double  value 
) [inline, inherited]

Set indexed value.

Definition at line 52 of file VectorN.hh.

References Vector< N >::_vec.

                                                           {
      if (index >= N) {
        throw std::runtime_error("Tried to access an invalid vector index.");
      } else {
        _vec[index] = value;
      }
      return *this;
    }
FourMomentum& setE ( double  E) [inline]

Set energy $ E $ (time component of momentum).

Definition at line 322 of file Vector4.hh.

References FourVector::setT().

Referenced by FourMomentum::FourMomentum().

{ setT(E); return *this; }
FourMomentum& setPx ( double  px) [inline]

Set x-component of momentum $ p_x $.

Definition at line 325 of file Vector4.hh.

References FourVector::setX().

Referenced by FourMomentum::FourMomentum().

{ setX(px); return *this; }
FourMomentum& setPy ( double  py) [inline]

Set y-component of momentum $ p_y $.

Definition at line 328 of file Vector4.hh.

References FourVector::setY().

Referenced by FourMomentum::FourMomentum().

{ setY(py); return *this; }
FourMomentum& setPz ( double  pz) [inline]

Set z-component of momentum $ p_z $.

Definition at line 331 of file Vector4.hh.

References FourVector::setZ().

Referenced by FourMomentum::FourMomentum().

{ setZ(pz); return *this; }
FourVector& setT ( const double  t) [inline, inherited]

Definition at line 56 of file Vector4.hh.

References FourVector::t().

Referenced by FourVector::FourVector(), and FourMomentum::setE().

{ set(0, t); return *this; }
FourVector& setX ( const double  x) [inline, inherited]

Definition at line 57 of file Vector4.hh.

References FourVector::x().

Referenced by FourVector::FourVector(), and FourMomentum::setPx().

{ set(1, x); return *this; }
FourVector& setY ( const double  y) [inline, inherited]

Definition at line 58 of file Vector4.hh.

References FourVector::y().

Referenced by FourVector::FourVector(), and FourMomentum::setPy().

{ set(2, y); return *this; }
FourVector& setZ ( const double  z) [inline, inherited]

Definition at line 59 of file Vector4.hh.

References FourVector::z().

Referenced by FourVector::FourVector(), and FourMomentum::setPz().

{ set(3, z); return *this; }
size_t size ( ) const [inline, inherited]

Vector dimensionality.

Definition at line 62 of file VectorN.hh.

                        {
      return N;
    }
double t ( ) const [inline, inherited]

Definition at line 52 of file Vector4.hh.

Referenced by FourVector::contract(), FourMomentum::E(), FourVector::invariant(), FourVector::setT(), and Rivet::toString().

{ return get(0); }
double theta ( ) const [inline, inherited]

Synonym for polarAngle.

Definition at line 124 of file Vector4.hh.

References Vector3::theta(), and FourVector::vector3().

Referenced by D0_1996_S3214044::_fourJetAnalysis(), and ParticleBase::theta().

                         {
      return vector3().theta();
    }
Vector3 vector3 ( ) const [inline, inherited]

Get the spatial part of the 4-vector as a 3-vector.

Definition at line 145 of file Vector4.hh.

Referenced by BABAR_2003_I593379::analyze(), BELLE_2001_S4598261::analyze(), ARGUS_1993_S2653028::analyze(), FourVector::angle(), Rivet::angle(), FourVector::azimuthalAngle(), FParameter::calc(), Hemispheres::calc(), Spherocity::calc(), Thrust::calc(), Sphericity::calc(), Rivet::cmpMomByAscP(), Rivet::cmpMomByP(), Rivet::deltaEta(), Rivet::deltaPhi(), Rivet::deltaR(), FourVector::eta(), ParticleBase::p3(), FourMomentum::p3(), FourVector::perp(), FourVector::perp2(), FourVector::phi(), FourVector::polarAngle(), FourVector::polarRadius(), FourVector::polarRadius2(), DISKinematics::project(), FoxWolframMoments::project(), MissingMomentum::project(), FourVector::pseudorapidity(), FourMomentum::pT2(), FourMomentum::pt2(), FourVector::rho(), FourVector::rho2(), and FourVector::theta().

                            {
      return Vector3(get(1), get(2), get(3));
    }
double x ( ) const [inline, inherited]

Definition at line 53 of file Vector4.hh.

Referenced by FourVector::contract(), FourVector::invariant(), Rivet::mT2::mT2(), ParticleBase::px(), FourMomentum::px(), FourVector::setX(), and Rivet::toString().

{ return get(1); }
double y ( ) const [inline, inherited]

Definition at line 54 of file Vector4.hh.

Referenced by FourVector::contract(), FourVector::invariant(), Rivet::mT2::mT2(), ParticleBase::py(), FourMomentum::py(), FourVector::setY(), and Rivet::toString().

{ return get(2); }
double z ( ) const [inline, inherited]

Definition at line 55 of file Vector4.hh.

Referenced by BeamThrust::_calcBeamThrust(), H1_1994_S2919893::analyze(), FourVector::contract(), FourVector::invariant(), ParticleBase::pz(), FourMomentum::pz(), FourVector::setZ(), and Rivet::toString().

{ return get(3); }

Friends And Related Function Documentation

FourMomentum add ( const FourMomentum &  a,
const FourMomentum &  b 
) [friend]

Definition at line 499 of file Vector4.hh.

Referenced by FourMomentum::operator+=(), and FourMomentum::operator-=().

                                                                        {
    FourMomentum result;
    result._vec = a._vec + b._vec;
    return result;
  }
FourMomentum multiply ( const double  a,
const FourMomentum &  v 
) [friend]

Definition at line 477 of file Vector4.hh.

Referenced by FourMomentum::operator*=(), and FourMomentum::operator/=().

                                                                      {
    FourMomentum result;
    result._vec = a * v._vec;
    return result;
  }
FourMomentum multiply ( const FourMomentum &  v,
const double  a 
) [friend]

Definition at line 483 of file Vector4.hh.

                                                                      {
    return multiply(a, v);
  }
FourMomentum transform ( const LorentzTransform &  lt,
const FourMomentum &  v4 
) [friend]

Referenced by CDF_1997_S3541940::_avg_beam_in_lab(), CDF_1996_S3349578::_avg_beam_in_lab(), CDF_1996_S3349578::_fiveJetAnalysis(), CDF_1996_S3349578::_fourJetAnalysis(), CDF_1996_S3349578::_threeJetAnalysis(), ARGUS_1993_S2653028::analyze(), ARGUS_1993_S2669951::analyze(), CDF_1997_S3541940::analyze(), ARGUS_1993_S2789213::analyze(), and MC_VH2BB::boostAngles().

Member Data Documentation

Eigen::Vector<double,N> _vec [protected, inherited]

Vector.

Definition at line 135 of file VectorN.hh.

Referenced by FourVector::operator*=(), FourMomentum::operator*=(), FourVector::operator+=(), FourMomentum::operator+=(), FourVector::operator-(), FourMomentum::operator-(), FourVector::operator-=(), FourMomentum::operator-=(), FourVector::operator/=(), and FourMomentum::operator/=().

The documentation for this class was generated from the following file:
Generated on Tue Mar 24 2015 17:35:41 for The Rivet MC analysis system by   doxygen 1.7.6.1