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 ()
virtual Location	location () const GCC11_FINAL
	Which part of the detector (barrel, endcap) More...
virtual const BoundDisk &	specificSurface () const GCC11_FINAL
virtual const BoundSurface &	surface () const GCC11_FINAL
	The surface of the GeometricSearchDet. More...
virtual	~ForwardDetLayer ()
Public Member Functions inherited from DetLayer
template<typename... Args>
std::vector< const DetLayer * >	compatibleLayers (Args &&...args) const
	Returns all layers compatible. More...
	DetLayer (bool ibar)
bool	isBarrel () const
bool	isForward () const
NavigableLayer *	navigableLayer () const
	Return the NavigableLayer associated with this DetLayer. More...
template<typename... Args>
std::vector< const DetLayer * >	nextLayers (Args &&...args) const
int	seqNum () const
void	setNavigableLayer (NavigableLayer *nlp)
	Set the NavigableLayer associated with this DetLayer. More...
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 ()
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
virtual bool	hasGroups () const =0
virtual const Surface::PositionType &	position () const
	Returns position of the surface. More...
virtual	~GeometricSearchDet ()

Protected Member Functions
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
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 20 of file ForwardDetLayer.h.

Constructor & Destructor Documentation

ForwardDetLayer::ForwardDetLayer ( )

inline

Definition at line 23 of file ForwardDetLayer.h.

: DetLayer(false),  theDisk(0){}

ForwardDetLayer::~ForwardDetLayer ( )

virtual

Definition at line 10 of file ForwardDetLayer.cc.

                                  {
}

Member Function Documentation

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.

Reimplemented in MuRingForwardDoubleLayer, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, and GCC11_FINAL< T, TOPO >.

Definition at line 104 of file ForwardDetLayer.cc.

References funct::abs(), deltaR(), TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), PV3DBase< T, PVType, FrameType >::perp(), LocalTrajectoryError::positionError(), Propagator::propagate(), mathSSE::sqrt(), tmp, unlikely, LocalError::xx(), LocalError::yy(), PV3DBase< T, PVType, FrameType >::z(), SiStripMonitorClusterAlca_cfi::zmax, SiStripMonitorClusterAlca_cfi::zmin, and zPos.

Referenced by MuRingForwardLayer::compatibleDets().

{
  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);

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

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

  float zPos = 0.5f*(zmax()+zmin());
  SimpleDiskBounds tmp( rmin()-deltaR, rmax()+deltaR, 
            zmin()-zPos, zmax()-zPos);

  return make_pair( tmp.inside(myState.localPosition()), myState);
}

BoundDisk * ForwardDetLayer::computeSurface ( )

protectedvirtual

Reimplemented in MuRingForwardDoubleLayer.

Definition at line 29 of file ForwardDetLayer.cc.

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

                                           {
  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 35 of file ForwardDetLayer.h.

References GeomDetEnumerators::endcap.

{return GeomDetEnumerators::endcap;}

float ForwardDetLayer::rmax ( ) const

inlineprotected

Definition at line 48 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuRingForwardDoubleLayer::compatible(), TECLayer::computeDisk(), and MuRingForwardDoubleLayer::computeSurface().

{ return theDisk->outerRadius();}

float ForwardDetLayer::rmin ( ) const

inlineprotected

Definition at line 47 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuRingForwardDoubleLayer::compatible(), TECLayer::computeDisk(), and MuRingForwardDoubleLayer::computeSurface().

{ return theDisk->innerRadius();}

void ForwardDetLayer::setSurface ( BoundDisk *  cp )

protected

Definition at line 15 of file ForwardDetLayer.cc.

References CommonMethods::cp().

Referenced by MuRingForwardDoubleLayer::MuRingForwardDoubleLayer(), MuRingForwardLayer::MuRingForwardLayer(), and TECLayer::TECLayer().

                                               { 
  theDisk = cp;
}

virtual const BoundDisk& ForwardDetLayer::specificSurface ( ) const

inlinevirtual

Reimplemented in GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, and GCC11_FINAL< T, TOPO >.

Definition at line 38 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuonNavigationSchool::addEndcapLayer(), MuonCSCDetLayerGeometryBuilder::buildLayer(), DirectMuonNavigation::checkCompatible(), DirectTrackerNavigation::checkCompatible(), MuRingForwardDoubleLayer::compatible(), MuRingForwardLayer::compatibleDets(), MuRingForwardDoubleLayer::computeSurface(), MuonShowerInformationFiller::crossingPoint(), MultipleScatteringGeometry::detLayers(), InnerDeltaPhi::initForwardLayer(), SimpleNavigationSchool::linkNextBarrelLayer(), OutInConversionSeedFinder::makeEstimator(), MSLayer::MSLayer(), MuRingForwardDoubleLayer::MuRingForwardDoubleLayer(), MuRingForwardLayer::MuRingForwardLayer(), DiskLessInnerRadius::operator()(), ForwardLayerMirrorImage::operator()(), TECLayer::TECLayer(), and SimpleNavigableLayer::wellInside().

{ return *theDisk;}

virtual const BoundSurface& ForwardDetLayer::surface ( ) const

inlinevirtual

The surface of the GeometricSearchDet.

Implements GeometricSearchDet.

Reimplemented in GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, and GCC11_FINAL< T, TOPO >.

Definition at line 28 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuRingForwardDoubleLayer::compatible(), MuRingForwardLayer::compatibleDets(), TECLayer::computeCrossings(), PixelTripletNoTipGenerator::hitTriplets(), ThirdHitRZPredictionBase::initLayer(), ThirdHitPrediction::initLayer(), OutInConversionSeedFinder::makeEstimator(), ValidHitPairFilter::operator()(), MuonDetLayerComp::operator()(), MuonNavigationPrinter::printLayers(), and SimpleNavigableLayer::wellInside().

{ return *theDisk;}

float ForwardDetLayer::zmax ( ) const

inlineprotected

Definition at line 50 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuRingForwardDoubleLayer::compatible(), TECLayer::computeDisk(), and MuRingForwardDoubleLayer::computeSurface().

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

float ForwardDetLayer::zmin ( ) const

inlineprotected

Definition at line 49 of file ForwardDetLayer.h.

References theDisk.

Referenced by MuRingForwardDoubleLayer::compatible(), TECLayer::computeDisk(), and MuRingForwardDoubleLayer::computeSurface().

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

Member Data Documentation

ReferenceCountingPointer<BoundDisk> ForwardDetLayer::theDisk

private

Definition at line 57 of file ForwardDetLayer.h.

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