CylindricalState.h

Go to the documentation of this file.

#ifndef CylindricalState_H
#define CylindricalState_H

#include "TrackPropagation/RungeKutta/interface/RKSmallVector.h"

#include <iostream>

class CylindricalState {
public:

  typedef double                                   Scalar;
  typedef RKSmallVector<Scalar,5>                  Vector;

  CylindricalState() {}

  CylindricalState( const LocalPoint& pos, const LocalVector& mom, Scalar ch) {
      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( Scalar rho, const Vector& par, Scalar prSign) :
      par_(par), rho_(rho), prSign_(prSign) {}


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

  const LocalVector momentum() const {
      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& parameters() const { return par_;}

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

  Scalar rho() const {return rho_;}

  double prSign() const {return prSign_;}

private:

  Vector par_;
  Scalar rho_;
  Scalar prSign_; 

};

#endif