ClosestApproachInRPhi.h

Go to the documentation of this file.

#ifndef _ClosestApproachInRPhi_H_
#define _ClosestApproachInRPhi_H_

#include "TrackingTools/PatternTools/interface/ClosestApproachOnHelices.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h" 

class ClosestApproachInRPhi : public ClosestApproachOnHelices {

public:

  ClosestApproachInRPhi() {status_ = false;}

  virtual bool calculate(const TrajectoryStateOnSurface & sta, 
             const TrajectoryStateOnSurface & stb);

  virtual bool calculate(const FreeTrajectoryState & sta,
             const FreeTrajectoryState & stb);

  virtual bool status() const {return status_;}

  virtual std::pair<GlobalPoint, GlobalPoint> points() const;

  std::pair <GlobalTrajectoryParameters, GlobalTrajectoryParameters >
  trajectoryParameters () const;

  virtual GlobalPoint crossingPoint() const;

  virtual float distance() const;

  virtual ClosestApproachInRPhi * clone() const {
    return new ClosestApproachInRPhi(* this);
  }

private:

  bool compute(const TrackCharge & chargeA, 
           const GlobalVector & momentumA, 
           const GlobalPoint & positionA, 
           const TrackCharge & chargeB, 
           const GlobalVector & momentumB, 
           const GlobalPoint & positionB);

  // given the old Parameters, and a new GlobalPoint,
  // we return the full new GlobalTrajectoryParameters at the
  // Point.
  static GlobalTrajectoryParameters
  newTrajectory( const GlobalPoint & newpt,
         const GlobalTrajectoryParameters & oldpar, double bz);

  // Computes center coordinates and unsigned radius of circle;
  static void circleParameters(const TrackCharge& charge, 
                   const GlobalVector& momemtum, 
                   const GlobalPoint& position, 
                   double& xc, double& yc, double& r,
                   double bz);

  // Computes crossing points of 2 circles with centres (cx_i, cy_i) 
  // and unsigned radii r_i. 
  // Two cases: - circles have one or two intersection points; 
  //              return value = 1; 
  //            - circles do not cross; computes point of closest approach 
  //              on each circle; return value = 2;
  // if the calculation fails (e.g. concentric circles), return value = 0;

  static int transverseCoord(double cxa, double cya, double ra, 
                 double cxb, double cyb, double rb, 
                 double & xg1, double & yg1, 
                 double & xg2, double & yg2);

  // Computes z-coordinate on helix at given transverse coordinates
  static double zCoord(const GlobalVector& mom, const GlobalPoint& pos, 
               double r, double xc, double yc, double xg, double yg);


private:
  GlobalPoint posA, posB;
  GlobalTrajectoryParameters paramA, paramB;
  double bz;
  bool status_;
  
};

#endif