ForwardDiskSectorBuilderFromWedges Class Reference

#include <ForwardDiskSectorBuilderFromWedges.h>

Public Member Functions
BoundDiskSector *	operator() (const std::vector< const TECWedge * > &wedges) const

Private Member Functions
std::pair< DiskSectorBounds *, GlobalVector >	computeBounds (const std::vector< const TECWedge * > &wedges) const
Surface::RotationType	computeRotation (const std::vector< const TECWedge * > &wedges, Surface::PositionType pos) const

Detailed Description

As it's name indicates, it's a builder of a BoundDiskSector from a collection of Wedges (of one petal). The disk sector has the minimal size fully containing all wedges.

Definition at line 15 of file ForwardDiskSectorBuilderFromWedges.h.

Member Function Documentation

pair< DiskSectorBounds *, GlobalVector > ForwardDiskSectorBuilderFromWedges::computeBounds ( const std::vector< const TECWedge * > &  wedges ) const

private

Definition at line 42 of file ForwardDiskSectorBuilderFromWedges.cc.

References funct::cos(), SiStripPI::max, min(), Geom::phiLess(), phimax, phimin, Geom::pi(), and funct::sin().

                                                 {
  // compute maximum and minimum radius and phi
  float rmin((**(wedges.begin())).specificSurface().innerRadius());
  float rmax(rmin);
  float zmin((**(wedges.begin())).surface().position().z());
  float zmax(zmin);
  float phimin((**(wedges.begin())).surface().position().phi());
  float phimax(phimin);

  for (vector<const TECWedge*>::const_iterator iw = wedges.begin(); iw != wedges.end(); iw++) {
    // edm::LogInfo(TkDetLayers) << "---------------------------------------------" ;
    // edm::LogInfo(TkDetLayers) <<   " Builder: Position of wedge     :" << (**iw).position() ;
    float ri = (**iw).specificSurface().innerRadius();
    float ro = (**iw).specificSurface().outerRadius();
    float zmi = (**iw).surface().position().z() - (**iw).specificSurface().bounds().thickness() / 2.;
    float zma = (**iw).surface().position().z() + (**iw).specificSurface().bounds().thickness() / 2.;
    float phi1 = (**iw).surface().position().phi() - (**iw).specificSurface().phiHalfExtension();
    float phi2 = (**iw).surface().position().phi() + (**iw).specificSurface().phiHalfExtension();
    rmin = min(rmin, ri);
    rmax = max(rmax, ro);
    zmin = min(zmin, zmi);
    zmax = max(zmax, zma);
    if (Geom::phiLess(phi1, phimin))
      phimin = phi1;
    if (Geom::phiLess(phimax, phi2))
      phimax = phi2;
  }

  if (!Geom::phiLess(phimin, phimax))
    edm::LogError("TkDetLayers") << " ForwardDiskSectorBuilderFromWedges : "
                                 << "Something went wrong with Phi Sorting !";
  float zPos = (zmax + zmin) / 2.;
  float phiWin = phimax - phimin;
  float phiPos = (phimax + phimin) / 2.;
  float rmed = (rmin + rmax) / 2.;
  if (phiWin < 0.) {
    if ((phimin < Geom::pi() / 2.) || (phimax > -Geom::pi() / 2.)) {
      edm::LogError("TkDetLayers") << " Debug: something strange going on, please check ";
    }
    // edm::LogInfo(TkDetLayers) << " Petal at pi: phi " << phimin << " " << phimax << " " << phiWin
    //   << " " << 2.*Geom::pi()+phiWin << " " ;
    phiWin += 2. * Geom::pi();
    phiPos += Geom::pi();
  }

  GlobalVector pos(rmed * cos(phiPos), rmed * sin(phiPos), zPos);
  return make_pair(new DiskSectorBounds(rmin, rmax, zmin - zPos, zmax - zPos, phiWin), pos);
}

Surface::RotationType ForwardDiskSectorBuilderFromWedges::computeRotation ( const std::vector< const TECWedge * > &  wedges, Surface::PositionType  pos  ) const

private

Definition at line 92 of file ForwardDiskSectorBuilderFromWedges.cc.

References Vector3DBase< T, FrameTag >::cross(), unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and MetAnalyzer::zAxis.

                                                                                                         {
  GlobalVector yAxis = (GlobalVector(pos.x(), pos.y(), 0.)).unit();

  GlobalVector zAxis(0., 0., 1.);
  GlobalVector xAxis = yAxis.cross(zAxis);

  return Surface::RotationType(xAxis, yAxis);
}

BoundDiskSector * ForwardDiskSectorBuilderFromWedges::operator()

(

const std::vector< const TECWedge * > &

wedges

)

const

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

Definition at line 9 of file ForwardDiskSectorBuilderFromWedges.cc.

References funct::abs(), mps_fire::i, PV3DBase< T, PVType, FrameType >::phi(), hitfit::phidiff(), Geom::phiLess(), Geom::pi(), GloballyPositioned< T >::position(), and makeMuonMisalignmentScenario::rot.

                                                                                                           {
  // check first that all wedges are about at the same phi and z !
  float phiStart = wedges.front()->position().phi();
  float zStart = wedges.front()->position().z();
  float wphimin, wphimax;
  for (vector<const TECWedge*>::const_iterator i = wedges.begin(); i != wedges.end(); i++) {
    float zdiff = (**i).surface().position().z() - zStart;
    if (std::abs(zdiff) > 5.)
      edm::LogError("TkDetLayers") << " ForwardDiskSectorBuilderFromWedges: Trying to build "
                                   << "Petal from Wedges at different z ! Delta Z = " << zdiff;
    float wphi = (**i).surface().position().phi();
    if (Geom::phiLess(phiStart, wphi)) {
      wphimin = phiStart;
      wphimax = wphi;
    } else {
      wphimin = wphi;
      wphimax = phiStart;
    }
    float phidiff = wphimax - wphimin;
    if (phidiff < 0.)
      phidiff += 2. * Geom::pi();
    if (phidiff > 0.3)
      edm::LogError("TkDetLayers") << " ForwardDiskSectorBuilderFromWedges: Trying to build "
                                   << "Petal from Wedges at different phi ! Delta phi = " << phidiff;
  }

  auto bo = computeBounds(wedges);

  Surface::PositionType pos(bo.second.x(), bo.second.y(), bo.second.z());
  Surface::RotationType rot = computeRotation(wedges, pos);
  return new BoundDiskSector(pos, rot, bo.first);
}