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 46 of file ForwardDiskSectorBuilderFromWedges.cc.

References funct::cos(), hpstanc_transforms::max, min(), 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 ( PhiLess()( phi1, phimin)) phimin = phi1;
    if ( PhiLess()( phimax, phi2)) phimax = phi2;
  }

  if (!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 97 of file ForwardDiskSectorBuilderFromWedges.cc.

References Vector3DBase< T, FrameTag >::cross(), csvLumiCalc::unit, PV3DBase< T, PVType, FrameType >::x(), MuonErrorMatrixValues_cff::xAxis, PV3DBase< T, PVType, FrameType >::y(), MuonErrorMatrixValues_cff::yAxis, 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 11 of file ForwardDiskSectorBuilderFromWedges.cc.

References funct::abs(), mps_fire::i, PV3DBase< T, PVType, FrameType >::phi(), hitfit::phidiff(), 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 ( 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);
}