IdealHelixParameters.h

Go to the documentation of this file.

#ifndef IDEALHELIXPARAMETERS_H
#define IDEALHELIXPARAMETERS_H

#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "TrackingTools/TrajectoryParametrization/interface/GlobalTrajectoryParameters.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

/* 
Given a track, evaluates the status of the track finding in the XY plane 
the point of tangent of the circle to the reference point(in general the PrimaryVertex or the BeamSpot)

The approach is the following.
Given a track, 
-extract the innermomentum at the innerpoint
-evalaute the radius of curvature and the center of the circle
-solving geometrical equation, evaluate the point of tangence to the circle, starting from the PV
-evaluate the status at this point of tangence T

 */
class IdealHelixParameters{

 public:

 IdealHelixParameters():
  _magnField(nullptr),_track(nullptr),
    _refPoint               (math::XYZVector(0,0,0)),
    _tangentPoint           (math::XYZVector(0,0,0)),
    _MomentumAtTangentPoint(math::XYZVector(0,0,0)){};
  ~IdealHelixParameters(){};

  inline void setMagnField(const MagneticField* magnField){_magnField=magnField;}
  inline void setData(const reco::Track* track, const math::XYZVector& refPoint=math::XYZVector(0,0,0));
  inline void setData(const reco::Track* track, const math::XYZPoint& ref);

  inline bool isTangentPointDistanceLessThan(float rmax, const reco::Track* track, const math::XYZVector& refPoint);

  math::XYZVector   GetCircleCenter() const {return _circleCenter;}
  math::XYZVector   GetTangentPoint() const {return _tangentPoint;}
  math::XYZVector   GetMomentumAtTangentPoint() const {return _MomentumAtTangentPoint;}
  float          GetTransverseIPAtTangent() const {return _transverseIP;}
  float          GetRotationAngle() const {return _rotationAngle;}
  
 private:

  inline void calculate();
  inline void evalCircleCenter();
  inline void evalTangentPoint();
  inline void evalMomentumatTangentPoint();

  const MagneticField *_magnField;
  const reco::Track   *_track;
  float  _radius;
  math::XYZVector  _circleCenter;
  math::XYZVector  _refPoint;
  math::XYZVector  _tangentPoint;
  math::XYZVector _MomentumAtTangentPoint;
  float _transverseIP;
  float _rotationAngle;

  //  std::stringstream ss;

};
#include "IdealHelixParameters.icc"

#endif