BarrelDetLayer Class Reference

#include <BarrelDetLayer.h>

Inheritance diagram for BarrelDetLayer:

Public Member Functions
	BarrelDetLayer (bool doHaveGroup)
virtual std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &ts, const Propagator &, const MeasurementEstimator &) const
bool	contains (const Local3DPoint &p) const
virtual Location	location () const
	DetLayer interface. More...
virtual const BoundCylinder &	specificSurface () const
	Extension of the interface. More...
virtual const BoundSurface &	surface () const
	GeometricSearchDet interface. More...
virtual	~BarrelDetLayer ()
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
SimpleCylinderBounds const &	bounds () const
virtual BoundCylinder *	computeSurface ()
virtual void	initialize ()
void	setSurface (BoundCylinder *cp)

Private Attributes
ReferenceCountingPointer < BoundCylinder >	theCylinder

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 barrel layers. Barrel layers are cylinders with their axes parallel to the global Z axis. The methods that have a common implementation for all BarrelDetLayers are implemented in this class, but some methods are left abstract.

Definition at line 23 of file BarrelDetLayer.h.

Constructor & Destructor Documentation

BarrelDetLayer::BarrelDetLayer ( bool  doHaveGroup )

inline

Definition at line 26 of file BarrelDetLayer.h.

                                   : DetLayer(doHaveGroup,true),
    theCylinder(0){}

BarrelDetLayer::~BarrelDetLayer ( )

virtual

Definition at line 10 of file BarrelDetLayer.cc.

{}

Member Function Documentation

SimpleCylinderBounds const& BarrelDetLayer::bounds ( ) const

inlineprotected

Definition at line 56 of file BarrelDetLayer.h.

References theCylinder.

{ return static_cast<SimpleCylinderBounds const &>(theCylinder->bounds());} 

pair< bool, TrajectoryStateOnSurface > BarrelDetLayer::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 74 of file BarrelDetLayer.cc.

References funct::abs(), TrajectoryStateOnSurface::cartesianError(), GlobalErrorBase< T, ErrorWeightType >::czz(), TrajectoryStateOnSurface::globalDirection(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::isValid(), PV3DBase< T, PVType, FrameType >::perp(), CartesianTrajectoryError::position(), Propagator::propagate(), mathSSE::sqrt(), unlikely, and PV3DBase< T, PVType, FrameType >::z().

{
  if unlikely(theCylinder == 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 assuming symmetric bounds around position().z()
  auto z0 = std::abs(myState.globalPosition().z()-specificSurface().position().z());
  auto deltaZ = 0.5f*bounds().length();
  if (z0<deltaZ) return make_pair( true, myState);

  // take into account the thickness of the layer
  deltaZ += 0.5f*bounds().thickness() *
    std::abs(myState.globalDirection().z())/myState.globalDirection().perp();

  // take also into account the error on the predicted state
  const float nSigma = 3.;
  if (myState.hasError())
    deltaZ += nSigma * sqrt( myState.cartesianError().position().czz());
  //
  //  check z again
  return make_pair(z0<deltaZ, myState);
}

BoundCylinder * BarrelDetLayer::computeSurface ( )

protectedvirtual

Definition at line 31 of file BarrelDetLayer.cc.

References Cylinder::computeRadius(), BoundingBox::corners(), F(), bookConverter::max, min(), alignCSCRings::r, makeMuonMisalignmentScenario::rot, and detailsBasic3DVector::z.

                                              {
  vector< const GeomDet*> comps = basicComponents();

  // Find extension in Z
  float theRmin = comps.front()->position().perp();   
  float theRmax = theRmin;
  float theZmin = comps.front()->position().z(); 
  float theZmax = theZmin;
  for ( vector< const GeomDet*>::const_iterator deti = comps.begin(); 
    deti != comps.end(); 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();
      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 +/- thickness/2
    // , since the min  radius for some barrel dets is reached there
    float rdet = (**deti).position().perp();
    float thick = (**deti).surface().bounds().thickness();
    theRmin = min( theRmin, rdet-thick/2.F);
    theRmax = max( theRmax, rdet+thick/2.F);
  }
  
  // By default the barrel layers are positioned at the center of the 
  // global frame, and the axes of their local frame coincide with 
  // those of the global grame (z along the cylinder axis)
  PositionType pos(0.,0.,0.);
  RotationType rot;

  auto scp = new SimpleCylinderBounds( theRmin, theRmax,
                                       theZmin, theZmax);
  return new Cylinder(Cylinder::computeRadius(*scp), pos, rot, scp);
}  

bool BarrelDetLayer::contains

(

const Local3DPoint &

p

)

const

Definition at line 19 of file BarrelDetLayer.cc.

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

void BarrelDetLayer::initialize ( )

protectedvirtual

Definition at line 23 of file BarrelDetLayer.cc.

Referenced by TBPLayer::construct(), MuRodBarrelLayer::MuRodBarrelLayer(), and TIBLayer::TIBLayer().

{
  setSurface( computeSurface());
}

virtual Location BarrelDetLayer::location ( ) const

inlinevirtual

DetLayer interface.

Implements DetLayer.

Definition at line 39 of file BarrelDetLayer.h.

References GeomDetEnumerators::barrel.

{return GeomDetEnumerators::barrel;}

void BarrelDetLayer::setSurface ( BoundCylinder *  cp )

protected

Definition at line 15 of file BarrelDetLayer.cc.

References CommonMethods::cp().

                                                  { 
  theCylinder = cp;
}

virtual const BoundCylinder& BarrelDetLayer::specificSurface ( ) const

inlinevirtual

Extension of the interface.

Definition at line 43 of file BarrelDetLayer.h.

References theCylinder.

Referenced by MuonNavigationSchool::addBarrelLayer(), GroupedCkfTrajectoryBuilder::advanceOneLayer(), CkfDebugger::analyseRecHitNotFound(), DirectMuonNavigation::checkCompatible(), DirectTrackerNavigation::checkCompatible(), PixelHitMatcher::compatibleHits(), MuonShowerInformationFiller::crossingPoint(), MultipleScatteringGeometry::detLayers(), RectangularEtaPhiTrackingRegion::estimator(), HitPairGeneratorFromLayerPairForPhotonConversion::getLayerRadius(), PixelTripletNoTipGenerator::hitTriplets(), InnerDeltaPhi::initBarrelLayer(), ThirdHitRZPredictionBase::initLayer(), ThirdHitPrediction::initLayer(), TkLayerLess::insideOutLess(), TkLayerLess::insideOutLessSigned(), SimpleNavigationSchool::linkNextForwardLayer(), OutInConversionSeedFinder::makeEstimator(), MSLayer::MSLayer(), MuonDetLayerComp::operator()(), ExtractBarrelDetLayerR::operator()(), MuonNavigationPrinter::printLayers(), TIBLayer::TIBLayer(), and SimpleNavigableLayer::wellInside().

{ return *theCylinder;}

virtual const BoundSurface& BarrelDetLayer::surface ( ) const

inlinevirtual

GeometricSearchDet interface.

Implements GeometricSearchDet.

Definition at line 32 of file BarrelDetLayer.h.

References theCylinder.

Referenced by TkLayerLess::barrelForwardLess(), RectangularEtaPhiTrackingRegion::estimator(), PixelTripletNoTipGenerator::hitTriplets(), ThirdHitRZPredictionBase::initLayer(), ThirdHitPrediction::initLayer(), SimpleNavigationSchool::linkNextForwardLayer(), MSLayer::MSLayer(), and ValidHitPairFilter::operator()().

{ return *theCylinder;}

Member Data Documentation

ReferenceCountingPointer<BoundCylinder> BarrelDetLayer::theCylinder

private

Definition at line 61 of file BarrelDetLayer.h.

Referenced by bounds(), specificSurface(), and surface().