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 34 of file ForwardRingDiskBuilderFromDet.cc.

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

{
  // 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 9 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);
}