CylindricalState Class Reference

#include <CylindricalState.h>

Public Types
typedef double	Scalar
typedef RKSmallVector< Scalar, 5 >	Vector

Public Member Functions
Scalar	charge () const
	CylindricalState ()
	CylindricalState (const LocalPoint &pos, const LocalVector &mom, Scalar ch)
	CylindricalState (Scalar rho, const Vector &par, Scalar prSign)
const LocalVector	momentum () const
const Vector &	parameters () const
const LocalPoint	position () const
double	prSign () const
Scalar	rho () const

Private Attributes
Vector	par_
Scalar	prSign_
	sign of local p_r More...
Scalar	rho_

Detailed Description

State for solving the equation of motion with radius (in cylindrical coordinates) as free variable. The dependent variables are phi - azimuthal angle z - z coordinate dphi/dr - derivative of phi versus r dz/dr - derivative of z versus r q/p - charge over momentum magnitude

The coordinate system is externally defined

Definition at line 23 of file CylindricalState.h.

Member Typedef Documentation

typedef double CylindricalState::Scalar

Definition at line 25 of file CylindricalState.h.

typedef RKSmallVector<Scalar, 5> CylindricalState::Vector

Definition at line 26 of file CylindricalState.h.

Constructor & Destructor Documentation

CylindricalState::CylindricalState ( )

inline

Definition at line 28 of file CylindricalState.h.

{}

CylindricalState::CylindricalState ( const LocalPoint &  pos, const LocalVector &  mom, Scalar  ch  )

inline

Definition at line 30 of file CylindricalState.h.

References RecoTauCleanerPlugins::charge, gather_cfg::cout, PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), position, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                             {
    rho_ = pos.perp();
    Scalar cosphi = pos.x() / rho_;
    Scalar sinphi = pos.y() / rho_;
    Scalar p_rho = mom.x() * cosphi + mom.y() * sinphi;
    Scalar p_phi = -mom.x() * sinphi + mom.y() * cosphi;

    par_(0) = pos.phi();
    par_(1) = pos.z();
    par_(2) = p_phi / (p_rho * rho_);
    par_(3) = mom.z() / p_rho;
    par_(4) = ch / mom.mag();

    prSign_ = p_rho > 0 ? 1.0 : -1.0;

    std::cout << "CylindricalState built from pos " << pos << " mom " << mom << " charge " << ch << std::endl;
    std::cout << "p_rho " << p_rho << " p_phi " << p_phi << " dphi_drho " << par_(2) << std::endl;
    std::cout << "Which results in                " << position() << " mom " << momentum() << " charge " << charge()
              << std::endl;
  }

CylindricalState::CylindricalState ( Scalar  rho, const Vector &  par, Scalar  prSign  )

inline

Definition at line 51 of file CylindricalState.h.

: par_(par), rho_(rho), prSign_(prSign) {}

Member Function Documentation

Scalar CylindricalState::charge ( void  ) const

inline

Definition at line 69 of file CylindricalState.h.

{ return par_(4) > 0 ? 1 : -1; }

const LocalVector CylindricalState::momentum ( ) const

inline

Definition at line 55 of file CylindricalState.h.

References funct::abs(), funct::cos(), mps_fire::result, funct::sin(), and mathSSE::sqrt().

Referenced by RKCylindricalDistance< T, N >::operator()().

                                     {
    Scalar cosphi = cos(par_(0));
    Scalar sinphi = sin(par_(0));
    Scalar Q = sqrt(1 + rho_ * rho_ * par_(2) * par_(2) + par_(3) * par_(3));
    Scalar P = std::abs(1. / par_(4));
    Scalar p_rho = prSign_ * P / Q;
    Scalar p_phi = rho_ * par_(2) * p_rho;
    Scalar p_z = par_(3) * p_rho;
    LocalVector result(p_rho * cosphi - p_phi * sinphi, p_rho * sinphi + p_phi * cosphi, p_z);
    return result;
  }

const Vector& CylindricalState::parameters ( void  ) const

inline

Definition at line 67 of file CylindricalState.h.

{ return par_; }

const LocalPoint CylindricalState::position ( ) const

inline

Definition at line 53 of file CylindricalState.h.

Referenced by RKCylindricalDistance< T, N >::operator()().

{ return LocalPoint(LocalPoint::Cylindrical(rho_, par_(0), par_(1))); }

double CylindricalState::prSign ( ) const

inline

Definition at line 73 of file CylindricalState.h.

{ return prSign_; }

Scalar CylindricalState::rho ( ) const

inline

Definition at line 71 of file CylindricalState.h.

Referenced by Lepton.Lepton::absIsoFromEA(), and Muon.Muon::absIsoWithFSR().

{ return rho_; }

Member Data Documentation

Vector CylindricalState::par_

private

Definition at line 76 of file CylindricalState.h.

Scalar CylindricalState::prSign_

private

sign of local p_r

Definition at line 78 of file CylindricalState.h.

Scalar CylindricalState::rho_

private

Definition at line 77 of file CylindricalState.h.