/data/refman/pasoursint/CMSSW_4_2_9_HLT1_bphpatch4/src/RecoTracker/NuclearSeedGenerator/src/TangentHelix.cc

Go to the documentation of this file.

#include "RecoTracker/NuclearSeedGenerator/interface/TangentHelix.h"

TangentHelix::TangentHelix(const TangentHelix& primCircle, const GlobalPoint& outerPoint, const GlobalPoint& innerPoint) :
     theInnerPoint(innerPoint), theOuterPoint(outerPoint), theCircle(primCircle.circle(), outerPoint, innerPoint) {

   // Calculation of vertex.z :
   GlobalPoint inner_T( innerPoint.x() , innerPoint.y() , 0.0 );
   GlobalPoint outer_T( outerPoint.x() , outerPoint.y() , 0.0 );
   GlobalPoint vtx_T( theCircle.vertexPoint().x() , theCircle.vertexPoint().y() , 0.0 );

   double d1 = (inner_T - vtx_T).mag();
   double d = (inner_T - outer_T).mag();

   theVertexPoint = GlobalPoint(vtx_T.x() , vtx_T.y(), innerPoint.z() - (outerPoint.z() - innerPoint.z()) * d1 / d );
   theDirectionAtVertex = GlobalVector(1000, 1000, 1000);
}

GlobalVector TangentHelix::directionAtVertex() {

   if(theDirectionAtVertex.z() > 999) {
      GlobalPoint inner_T( theInnerPoint.x() , theInnerPoint.y() , 0.0 );
      GlobalPoint outer_T( theOuterPoint.x() , theOuterPoint.y() , 0.0 );
      double p_z = (theOuterPoint.z() - theInnerPoint.z()) / (outer_T - inner_T).mag();

      GlobalVector dir_T = theCircle.directionAtVertex();
      GlobalVector dir( dir_T.x(), dir_T.y(), p_z);

      dir/=dir.mag();
      theDirectionAtVertex = dir;
    }

   return theDirectionAtVertex;
}