ForwardDetLayer Class Reference

#include <ForwardDetLayer.h>

Inheritance diagram for ForwardDetLayer:

Public Member Functions
virtual std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &, const Propagator &, const MeasurementEstimator &) const
bool	contains (const Local3DPoint &p) const
	ForwardDetLayer (bool doHaveGroups)
virtual Location	location () const
	Which part of the detector (barrel, endcap) More...
virtual const BoundDisk &	specificSurface () const
virtual const BoundSurface &	surface () const
	The surface of the GeometricSearchDet. More...
virtual	~ForwardDetLayer ()
Public Member Functions inherited from DetLayer
	DetLayer (bool doHaveGroup, bool ibar)
bool	isBarrel () const
bool	isForward () const
int	seqNum () const
void	setSeqNum (int sq)
virtual SubDetector	subDetector () const =0
	The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap) More...
virtual	~DetLayer ()
Public Member Functions inherited from GeometricSearchDet
virtual const std::vector < const GeomDet * > &	basicComponents () const =0
virtual std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
virtual void	compatibleDetsV (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetWithState > &result) const
virtual const std::vector < const GeometricSearchDet * > &	components () const =0
	Returns basic components, if any. More...
	GeometricSearchDet (bool doHaveGroups)
virtual std::vector< DetGroup >	groupedCompatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
virtual void	groupedCompatibleDetsV (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetGroup > &result) const
bool	hasGroups () const
virtual const Surface::PositionType &	position () const
	Returns position of the surface. More...
virtual	~GeometricSearchDet ()

Protected Member Functions
SimpleDiskBounds const &	bounds () const
virtual BoundDisk *	computeSurface ()
virtual void	initialize ()
float	rmax () const
float	rmin () const
void	setSurface (BoundDisk *cp)
float	zmax () const
float	zmin () const

Private Attributes
ReferenceCountingPointer < BoundDisk >	theDisk

Additional Inherited Members
Public Types inherited from DetLayer
typedef GeomDetEnumerators::Location	Location
typedef GeomDetEnumerators::SubDetector	SubDetector
Public Types inherited from GeometricSearchDet
typedef std::pair< const GeomDet *, TrajectoryStateOnSurface >	DetWithState
typedef BoundSurface::PositionType	PositionType
typedef BoundSurface::RotationType	RotationType
typedef TrajectoryStateOnSurface	TrajectoryState
Protected Attributes inherited from GeometricSearchDet
bool	haveGroups
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A specialization of the DetLayer interface for forward layers. Forward layers are disks with their axes parallel to the global Z axis. The methods that have a common implementation for all ForwardDetLayers are implemented in this class, but some methods are left abstract.

Definition at line 23 of file ForwardDetLayer.h.

Constructor & Destructor Documentation

ForwardDetLayer::ForwardDetLayer ( bool  doHaveGroups )

inline

Definition at line 26 of file ForwardDetLayer.h.

: DetLayer(doHaveGroups,false) {}

ForwardDetLayer::~ForwardDetLayer ( )

virtual

Definition at line 10 of file ForwardDetLayer.cc.

                                  {
}

Member Function Documentation

SimpleDiskBounds const& ForwardDetLayer::bounds ( ) const

inlineprotected

Definition at line 58 of file ForwardDetLayer.h.

References theDisk.

Referenced by zmax(), and zmin().

{ return static_cast<SimpleDiskBounds const &>(theDisk->bounds());} 

pair< bool, TrajectoryStateOnSurface > ForwardDetLayer::compatible ( const TrajectoryStateOnSurface &  ts, const Propagator &  , const MeasurementEstimator &    ) const

virtual

tests the geometrical compatibility of the Det with the predicted state. The FreeTrajectoryState argument is propagated to the Det surface using the Propagator argument. The resulting TrajectoryStateOnSurface is tested for compatibility with the surface bounds. If compatible, a std::pair< true, propagatedState> is returned. If the propagation fails, or if the state is not compatible, a std::pair< false, propagatedState> is returned.

Implements GeometricSearchDet.

Definition at line 104 of file ForwardDetLayer.cc.

References funct::abs(), deltaR(), f, TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), bookConverter::max, HLT_25ns14e33_v1_cff::nSigma, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::perp2(), LocalTrajectoryError::positionError(), Propagator::propagate(), diffTwoXMLs::r2, mathSSE::sqrt(), unlikely, LocalError::xx(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

{
  if unlikely(theDisk == 0)  edm::LogError("DetLayers") 
    << "ERROR: BarrelDetLayer::compatible() is used before the layer surface is initialized" ;
  // throw an exception? which one?

  TrajectoryStateOnSurface myState = prop.propagate( ts, specificSurface());
  if unlikely( !myState.isValid()) return make_pair( false, myState);


  // check z;  (are we sure?????)
  // auto z = myState.localPosition().z();
  // if ( z<bounds().minZ | z>bounds().maxZ) return make_pair( false, myState);

  // check r
  auto r2 =  myState.localPosition().perp2();
  if ( (r2 > rmin()*rmin()) & (r2< rmax()*rmax()) ) return make_pair( true, myState);


  // take into account the thickness of the layer
  float deltaR = 0.5f*bounds().thickness() *
    myState.localDirection().perp()/std::abs(myState.localDirection().z());

  // and take into account the error on the predicted state
  const float nSigma = 3.;
  if (myState.hasError()) {
    LocalError err = myState.localError().positionError();
    // ignore correlation for the moment...
    deltaR += nSigma * sqrt(err.xx() + err.yy());
  }

  // check r again;
  auto ri2 = std::max(rmin()-deltaR,0.f); ri2*=ri2;
  auto ro2 = rmax()+deltaR; ro2*=ro2;
  return make_pair( (r2>ri2) & (r2<ro2), myState);
}

BoundDisk * ForwardDetLayer::computeSurface ( )

protectedvirtual

Definition at line 29 of file ForwardDetLayer.cc.

References makeMuonMisalignmentScenario::components, BoundingBox::corners(), Vector3DBase< T, FrameTag >::dot(), F(), LogDebug, bookConverter::max, min(), alignCSCRings::r, makeMuonMisalignmentScenario::rot, Vector3DBase< T, FrameTag >::unit(), create_public_lumi_plots::width, detailsBasic3DVector::z, and zPos.

Referenced by PixelForwardLayer::PixelForwardLayer(), and PixelForwardLayerPhase1::PixelForwardLayerPhase1().

                                           {
  LogDebug("DetLayers") << "ForwaLayer::computeSurface callded" ;
  vector<const GeomDet*> comps= basicComponents();

  vector<const GeomDet*>::const_iterator ifirst = comps.begin();
  vector<const GeomDet*>::const_iterator ilast  = comps.end();

  // Find extension in R
  float theRmin = components().front()->position().perp(); 
  float theRmax = theRmin;
  float theZmin = components().back()->position().z();
  float theZmax = theZmin;
  for ( vector<const GeomDet*>::const_iterator deti = ifirst;
    deti != ilast; deti++) {
    vector<GlobalPoint> corners = 
      BoundingBox().corners( dynamic_cast<const Plane&>((**deti).surface()));
    for (vector<GlobalPoint>::const_iterator ic = corners.begin();
     ic != corners.end(); ic++) {
      float r = ic->perp();
      LogDebug("DetLayers") << "corner.perp(): " << r ;
      float z = ic->z();
      theRmin = min( theRmin, r);
      theRmax = max( theRmax, r);
      theZmin = min( theZmin, z);
      theZmax = max( theZmax, 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  = (**deti).position().perp();
    float len   = (**deti).surface().bounds().length();
    float width = (**deti).surface().bounds().width();

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

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

    LogDebug("DetLayers") << "in ForwardDetLayer::computeSurface(),xAxisCos,yAxisCos: " << xAxisCos << " , " << yAxisCos ;
    LogDebug("DetLayers") << "det pos.perp,length,width: " 
              << rdet << " , " 
              << len  << " , "
              << width ;

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


  LogDebug("DetLayers") << "creating SimpleDiskBounds with r range" << theRmin << " " 
            << theRmax << " and z range " << theZmin << " " << theZmax ;

  // By default the forward layers are positioned around the z axis of the
  // global frame, and the axes of their local frame coincide with 
  // those of the global grame (z along the disk axis)
  float zPos = (theZmax+theZmin)/2.;
  PositionType pos(0.,0.,zPos);
  RotationType rot;

  return new BoundDisk( pos, rot, 
            new SimpleDiskBounds( theRmin, theRmax, 
                              theZmin-zPos, theZmax-zPos));
}  

bool ForwardDetLayer::contains

(

const Local3DPoint &

p

)

const

Definition at line 19 of file ForwardDetLayer.cc.

                                                          {
  return surface().bounds().inside(p);
}

void ForwardDetLayer::initialize ( )

protectedvirtual

Definition at line 24 of file ForwardDetLayer.cc.

                                 {
  setSurface( computeSurface());
}

virtual Location ForwardDetLayer::location ( ) const

inlinevirtual

Which part of the detector (barrel, endcap)

Implements DetLayer.

Definition at line 38 of file ForwardDetLayer.h.

References GeomDetEnumerators::endcap.

{return GeomDetEnumerators::endcap;}

float ForwardDetLayer::rmax ( ) const

inlineprotected

Definition at line 51 of file ForwardDetLayer.h.

References theDisk.

{ return theDisk->outerRadius();}

float ForwardDetLayer::rmin ( ) const

inlineprotected

Definition at line 50 of file ForwardDetLayer.h.

References theDisk.

{ return theDisk->innerRadius();}

void ForwardDetLayer::setSurface ( BoundDisk *  cp )

protected

Definition at line 15 of file ForwardDetLayer.cc.

References CommonMethods::cp().

Referenced by PixelForwardLayer::PixelForwardLayer(), PixelForwardLayerPhase1::PixelForwardLayerPhase1(), and TECLayer::TECLayer().

                                               { 
  theDisk = cp;
}

virtual const BoundDisk& ForwardDetLayer::specificSurface ( ) const

inlinevirtual

Definition at line 41 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuonNavigationSchool::addEndcapLayer(), DirectMuonNavigation::checkCompatible(), DirectTrackerNavigation::checkCompatible(), PixelForwardLayer::computeCrossings(), PixelForwardLayerPhase1::computeCrossings(), MuonShowerInformationFiller::crossingPoint(), MultipleScatteringGeometry::detLayers(), RectangularEtaPhiTrackingRegion::estimator(), InnerDeltaPhi::initForwardLayer(), SimpleNavigationSchool::linkNextBarrelLayer(), OutInConversionSeedFinder::makeEstimator(), SymmetricLayerFinder::mirrorPartner(), MSLayer::MSLayer(), DiskLessInnerRadius::operator()(), PixelForwardLayer::PixelForwardLayer(), PixelForwardLayerPhase1::PixelForwardLayerPhase1(), TECLayer::TECLayer(), and SimpleNavigableLayer::wellInside().

{ return *theDisk;}

virtual const BoundSurface& ForwardDetLayer::surface ( ) const

inlinevirtual

The surface of the GeometricSearchDet.

Implements GeometricSearchDet.

Definition at line 31 of file ForwardDetLayer.h.

References theDisk.

Referenced by PixelForwardLayer::computeCrossings(), PixelForwardLayerPhase1::computeCrossings(), RectangularEtaPhiTrackingRegion::estimator(), PixelTripletNoTipGenerator::hitTriplets(), ThirdHitRZPredictionBase::initLayer(), ThirdHitPrediction::initLayer(), OutInConversionSeedFinder::makeEstimator(), ValidHitPairFilter::operator()(), MuonDetLayerComp::operator()(), PixelForwardLayer::PixelForwardLayer(), PixelForwardLayerPhase1::PixelForwardLayerPhase1(), MuonNavigationPrinter::printLayers(), TECLayer::searchNeighbors(), and SimpleNavigableLayer::wellInside().

{ return *theDisk;}

float ForwardDetLayer::zmax ( ) const

inlineprotected

Definition at line 53 of file ForwardDetLayer.h.

References bounds(), theDisk, and SimpleDiskBounds::thickness().

{ return (theDisk->position().z() + bounds().thickness()*0.5f);}

float ForwardDetLayer::zmin ( ) const

inlineprotected

Definition at line 52 of file ForwardDetLayer.h.

References bounds(), theDisk, and SimpleDiskBounds::thickness().

{ return (theDisk->position().z() - bounds().thickness()*0.5f);}

Member Data Documentation

ReferenceCountingPointer<BoundDisk> ForwardDetLayer::theDisk

private

Definition at line 61 of file ForwardDetLayer.h.

Referenced by bounds(), rmax(), rmin(), specificSurface(), surface(), zmax(), and zmin().