StripTopologyBuilder.cc

Go to the documentation of this file.

#include "Geometry/TrackerGeometryBuilder/interface/StripTopologyBuilder.h"
#include "Geometry/CommonTopologies/interface/RectangularStripTopology.h"
#include "Geometry/CommonTopologies/interface/RadialStripTopology.h"
#include "DataFormats/GeometrySurface/interface/Bounds.h"
#include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"



StripTopologyBuilder::StripTopologyBuilder(){}

StripTopology* StripTopologyBuilder::build(const Bounds* bs,double apvnumb,std::string part)
{
  theAPVNumb = apvnumb;

  StripTopology* result;
  if (part == "barrel") {
    result = constructBarrel( bs->length(), bs->width());
  }
  else {
    int yAx = (dynamic_cast<const TrapezoidalPlaneBounds*>(bs))->yAxisOrientation();
    result = constructForward( bs->length(), bs->width(),bs->widthAtHalfLength(),yAx);
  }
  return result;
}

StripTopology* StripTopologyBuilder::constructBarrel( float length, float width)
{
  int nstrip = int(128*theAPVNumb);
  float pitch = width/nstrip;
  
  return new RectangularStripTopology(nstrip,pitch,length);
}

StripTopology* StripTopologyBuilder::constructForward( float length, float width, float widthAtHalf, int yAxOr)
{
  int nstrip = int(128*theAPVNumb);
  float rCross = widthAtHalf*length/(2*(width-widthAtHalf));
  float aw = atan2(widthAtHalf/2., static_cast<double>(rCross))/(nstrip/2);
  return new RadialStripTopology(nstrip, aw, length, rCross, yAxOr);
}