ForwardRingDiskBuilderFromDet Class Reference

#include <ForwardRingDiskBuilderFromDet.h>

Public Member Functions
std::pair< SimpleDiskBounds *, float >	computeBounds (const std::vector< const GeomDet * > &dets) const
BoundDisk *	operator() (const std::vector< const GeomDet * > &dets) const

Detailed Description

As it's name indicates, it's a builder of BoundDisk from a collection of Dets. The disk has the minimal size fully containing all Dets.

Definition at line 19 of file ForwardRingDiskBuilderFromDet.h.

Member Function Documentation

pair< SimpleDiskBounds *, float > ForwardRingDiskBuilderFromDet::computeBounds

(

const std::vector< const GeomDet * > &

dets

)

const

Definition at line 31 of file ForwardRingDiskBuilderFromDet.cc.

References BoundingBox::corners(), Vector3DBase< T, FrameTag >::dot(), F(), mps_fire::i, SiStripPI::max, min(), alignCSCRings::r, and Vector3DBase< T, FrameTag >::unit().

                                                                                                                    {
  // go over all corners and compute maximum deviations from mean pos.
  float rmin((**(dets.begin())).surface().position().perp());
  float rmax(rmin);
  float zmin((**(dets.begin())).surface().position().z());
  float zmax(zmin);
  for (vector<const GeomDet*>::const_iterator idet = dets.begin(); idet != dets.end(); idet++) {
    /* ---- original implementation. Is it obsolete?
    vector<DetUnit*> detUnits = (**idet).detUnits();
    for (vector<DetUnit*>::const_iterator detu=detUnits.begin();
    detu!=detUnits.end(); detu++) {
    vector<GlobalPoint> corners = BoundingBox().corners(
    dynamic_cast<const Plane&>((**detu).surface()));
    }
    ----- */
    vector<GlobalPoint> corners = BoundingBox().corners((**idet).specificSurface());
    for (vector<GlobalPoint>::const_iterator i = corners.begin(); i != corners.end(); i++) {
      float r = i->perp();
      float z = i->z();
      rmin = min(rmin, r);
      rmax = max(rmax, r);
      zmin = min(zmin, z);
      zmax = max(zmax, z);
    }
    // in addition to the corners we have to check the middle of the
    // det +/- length/2, since the min (max) radius for typical fw
    // dets is reached there

    float rdet = (**idet).position().perp();
    float len = (**idet).surface().bounds().length();
    float width = (**idet).surface().bounds().width();

    GlobalVector xAxis = (**idet).toGlobal(LocalVector(1, 0, 0));
    GlobalVector yAxis = (**idet).toGlobal(LocalVector(0, 1, 0));
    GlobalVector perpDir = GlobalVector((**idet).position() - GlobalPoint(0, 0, (**idet).position().z()));

    double xAxisCos = xAxis.unit().dot(perpDir.unit());
    double yAxisCos = yAxis.unit().dot(perpDir.unit());

    if (fabs(xAxisCos) > fabs(yAxisCos)) {
      rmin = min(rmin, rdet - width / 2.F);
      rmax = max(rmax, rdet + width / 2.F);
    } else {
      rmin = min(rmin, rdet - len / 2.F);
      rmax = max(rmax, rdet + len / 2.F);
    }
  }

  float zPos = (zmax + zmin) / 2.;
  return make_pair(new SimpleDiskBounds(rmin, rmax, zmin - zPos, zmax - zPos), zPos);
}

BoundDisk * ForwardRingDiskBuilderFromDet::operator()

(

const std::vector< const GeomDet * > &

dets

)

const

Warning, remember to assign this pointer to a ReferenceCountingPointer! Should be changed to return a ReferenceCountingPointer<BoundDisk>

Definition at line 8 of file ForwardRingDiskBuilderFromDet.cc.

References makeMuonMisalignmentScenario::rot.

                                                                                             {
  auto bo = computeBounds(dets);

  //   LogDebug("DetLayers") << "Creating disk at Z: " << bo.second << "\n"
  //        << "Bounds are (rmin/rmax/thick) " << bo.first.innerRadius()
  //        << " / " <<  bo.first.outerRadius()
  //        << " / " <<  bo.first.thickness()  ;

  //   typedef Det::PositionType::BasicVectorType Vector;
  //   Vector posSum(0,0,0);
  //   for (vector<Det*>::const_iterator i=dets.begin(); i!=dets.end(); i++) {
  //     Vector pp = (**i).position().basicVector();
  //     //    LogDebug("DetLayers") << "  "<< (int) ( i-dets.begin()) << " at " << pp ;
  //     posSum += pp;
  //   }
  //   Det::PositionType meanPos( posSum/float(dets.size()));
  //   LogDebug("DetLayers") << "  meanPos "<< meanPos ;

  Surface::PositionType pos(0., 0., bo.second);
  Surface::RotationType rot;
  return new BoundDisk(pos, rot, bo.first);
}