TRT_SeededSpacePointFinder_ATL.cxx File Reference

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/ServiceHandle.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "TrkSpacePoint/SpacePointCLASS_DEF.h"
#include "TRT_SeededSpacePointFinderTool/TRT_SeededSpacePointFinder_ATL.h"
#include "SiSpacePointsSeed/SiSpacePointsSeed.h"
#include "InDetPrepRawData/SCT_Cluster.h"
#include "InDetIdentifier/SCT_ID.h"
#include "StoreGate/ReadHandle.h"
#include <ostream>
#include <iomanip>
#include <set>

Go to the source code of this file.

Functions
double	rollrating (double angle)
void	bracket_angle (double angle, double delta, double *min, double *max)
double	rotrating (double y, double x)

Function Documentation

bracket_angle()

void bracket_angle	(	double	angle,
		double	delta,
		double *	min,
		double *	max )

Definition at line 612 of file TRT_SeededSpacePointFinder_ATL.cxx.

                                    {
  double amin = rollrating(angle - delta);
  double amax = rollrating(angle + delta);
  if (delta >= M_PI) {
    amin = 0.0;
    amax = 4.0;
  }
  *min = amin;
  *max = amax;
}

rollrating()

double rollrating

(

double

angle

)

Definition at line 582 of file TRT_SeededSpacePointFinder_ATL.cxx.

                         {
  while (angle < 0) {
    angle += 2*M_PI;
  }
  while (angle > 2*M_PI) {
    angle -= 2*M_PI;
  }
  double rotations = angle/(2*M_PI);
  double rquadrant = rotations*4.0;
  long quadrant = (long)rquadrant & 3;
  double twist;
  if ((quadrant & 1) != 0) {
    twist = cos(angle);
  }
  else {
    twist = sin(angle);
  }
  twist *= twist;
  quadrant -= ((quadrant & 2) << 1);
  return quadrant + twist;
}

rotrating()

double rotrating ( double y, double x )

inline

Definition at line 629 of file TRT_SeededSpacePointFinder_ATL.cxx.

                              {
  long asign_x = (long)(x < 0.0);
  long asign_y = (long)(y < 0.0);
  long quadrant = -(asign_y << 1) + (asign_y ^ asign_x);
  double x2 = x*x;
  double y2 = y*y;
  double denominator = x2 + y2;
  double numerator = ((quadrant & 1) != 0) ? x2 : y2;
  return (double)quadrant + numerator/denominator;
}