TrackerBounds.cc

Go to the documentation of this file.

#include "TrackingTools/GeomPropagators/interface/TrackerBounds.h"

#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"
#include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"
#include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"

//Ported from ORCA

static const float epsilon = 0.001; // should not matter at all

static Cylinder * initCylinder()
{
  Surface::RotationType rot; // unit rotation matrix
  auto cb = new SimpleCylinderBounds( TrackerBounds::radius()-epsilon, 
                      TrackerBounds::radius()+epsilon, 
                      -TrackerBounds::halfLength(),  TrackerBounds::halfLength()
                                     );
 return new Cylinder(Cylinder::computeRadius(*cb), Surface::PositionType(0,0,0), rot, cb);
}

static Disk* initNegative()
{
  Surface::RotationType rot; // unit rotation matrix

  return new Disk( Surface::PositionType( 0, 0, -TrackerBounds::halfLength()), rot, 
           new SimpleDiskBounds( 0, TrackerBounds::radius(), -epsilon, epsilon));
}

static Disk* initPositive()
{
  Surface::RotationType rot; // unit rotation matrix


  return new Disk( Surface::PositionType( 0, 0, TrackerBounds::halfLength()), rot, 
           new SimpleDiskBounds( 0, TrackerBounds::radius(), -epsilon, epsilon));
}

bool TrackerBounds::isInside(const GlobalPoint &point){
  return (point.perp() <= TrackerBounds::radius() &&
      fabs(point.z()) <= TrackerBounds::halfLength());
}


// static initializers

const ReferenceCountingPointer<BoundCylinder>  TrackerBounds::theCylinder = initCylinder();
const ReferenceCountingPointer<BoundDisk>      TrackerBounds::theNegativeDisk = initNegative();
const ReferenceCountingPointer<BoundDisk>      TrackerBounds::thePositiveDisk = initPositive();