PixelBarrelLayer Class Reference

#include <PixelBarrelLayer.h>

Inheritance diagram for PixelBarrelLayer:

Public Types
typedef PeriodicBinFinderInPhi < double >	BinFinderType
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
Public Types inherited from GeometricSearchDetWithGroups
typedef GeometricSearchDet::DetWithState	DetWithState

Public Member Functions
virtual const std::vector < const GeomDet * > &	basicComponents () const
virtual const std::vector < const GeometricSearchDet * > &	components () const
	Returns basic components, if any. More...
void	groupedCompatibleDetsV (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetGroup > &result) const
	PixelBarrelLayer (std::vector< const PixelRod * > &innerRods, std::vector< const PixelRod * > &outerRods)
virtual SubDetector	subDetector () const
	The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap) More...
	~PixelBarrelLayer ()
Public Member Functions inherited from RodBarrelLayer
	RodBarrelLayer ()
virtual	~RodBarrelLayer ()
Public Member Functions inherited from BarrelDetLayer
	BarrelDetLayer ()
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
std::vector< const DetLayer * >	compatibleLayers (NavigationDirection direction) const
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection, int &counter) const
	DetLayer ()
NavigableLayer *	navigableLayer () const
	Return the NavigableLayer associated with this DetLayer. More...
virtual std::vector< const DetLayer * >	nextLayers (NavigationDirection direction) const
virtual std::vector< const DetLayer * >	nextLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const
virtual void	setNavigableLayer (NavigableLayer *nlp)
	Set the NavigableLayer associated with this DetLayer. More...
virtual	~DetLayer ()
Public Member Functions inherited from GeometricSearchDet
virtual std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
	GeometricSearchDet ()
virtual std::vector< DetGroup >	groupedCompatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
virtual const Surface::PositionType &	position () const
	Returns position of the surface. More...
virtual	~GeometricSearchDet ()
Public Member Functions inherited from GeometricSearchDetWithGroups
void	compatibleDetsV (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetWithState > &result) const
bool	hasGroups () const

Private Member Functions
bool	addClosest (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, std::vector< DetGroup > &result) const
double	calculatePhiWindow (double Xmax, const GeomDet &det, const TrajectoryStateOnSurface &state) const
SubLayerCrossings	computeCrossings (const TrajectoryStateOnSurface &tsos, PropagationDirection propDir) const
float	computeWindowSize (const GeomDet *det, const TrajectoryStateOnSurface &tsos, const MeasurementEstimator &est) const
BoundCylinder *	cylinder (const std::vector< const GeometricSearchDet * > &rods) const
bool	overlap (const GlobalPoint &gpos, const GeometricSearchDet &rod, float phiWin) const
void	searchNeighbors (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, float window, std::vector< DetGroup > &result, bool checkClosest) const
const std::vector< const GeometricSearchDet * > &	subLayer (int ind) const

Private Attributes
std::vector< const GeomDet * >	theBasicComps
std::vector< const GeometricSearchDet * >	theComps
BinFinderType	theInnerBinFinder
std::vector< const GeometricSearchDet * >	theInnerComps
ReferenceCountingPointer < BoundCylinder >	theInnerCylinder
BinFinderType	theOuterBinFinder
std::vector< const GeometricSearchDet * >	theOuterComps
ReferenceCountingPointer < BoundCylinder >	theOuterCylinder

Additional Inherited Members
Protected Member Functions inherited from BarrelDetLayer
virtual BoundCylinder *	computeSurface ()
virtual void	initialize ()
void	setSurface (BoundCylinder *cp)
Protected Attributes inherited from GeometricSearchDet
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A concrete implementation for PixelBarrel layer built out of BarrelPixelRod

Definition at line 17 of file PixelBarrelLayer.h.

Member Typedef Documentation

typedef PeriodicBinFinderInPhi<double> PixelBarrelLayer::BinFinderType

Definition at line 19 of file PixelBarrelLayer.h.

Constructor & Destructor Documentation

PixelBarrelLayer::PixelBarrelLayer	(	std::vector< const PixelRod * > &	innerRods,
		std::vector< const PixelRod * > &	outerRods
	)

Definition at line 20 of file PixelBarrelLayer.cc.

References cylinder(), i, BarrelDetLayer::initialize(), LogDebug, theBasicComps, theComps, theInnerBinFinder, theInnerComps, theInnerCylinder, theOuterBinFinder, theOuterComps, and theOuterCylinder.

                                                       : 
  theInnerComps(innerRods.begin(),innerRods.end()), 
  theOuterComps(outerRods.begin(),outerRods.end())
{
  theComps.assign(theInnerComps.begin(),theInnerComps.end());
  theComps.insert(theComps.end(),theOuterComps.begin(),theOuterComps.end());

  for(vector<const GeometricSearchDet*>::const_iterator it=theComps.begin();
      it!=theComps.end();it++){  
    theBasicComps.insert(theBasicComps.end(),   
             (**it).basicComponents().begin(),
             (**it).basicComponents().end());
  }

  theInnerCylinder = cylinder( theInnerComps);
  theOuterCylinder = cylinder( theOuterComps);

  theInnerBinFinder = BinFinderType(theInnerComps.front()->position().phi(),
                    theInnerComps.size());
  theOuterBinFinder = BinFinderType(theOuterComps.front()->position().phi(),
                    theOuterComps.size());

  
  BarrelDetLayer::initialize();


  //--------- DEBUG INFO --------------
  LogDebug("TkDetLayers") << "==== DEBUG PixelBarrelLayer =====" ; 
  LogDebug("TkDetLayers") << "PixelBarrelLayer innerCyl r,lenght: "
                          << theInnerCylinder->radius() << " , "
                          << theInnerCylinder->bounds().length();

  LogDebug("TkDetLayers") << "PixelBarrelLayer outerCyl r,lenght: " 
              << theOuterCylinder->radius() << " , "
              << theOuterCylinder->bounds().length();


    for (vector<const GeometricSearchDet*>::const_iterator i=theInnerComps.begin();
       i != theInnerComps.end(); i++){
    LogDebug("TkDetLayers") << "inner PixelRod pos z,perp,eta,phi: " 
                << (**i).position().z()    << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta()  << " , " 
                << (**i).position().phi()  ;
  }
  
  for (vector<const GeometricSearchDet*>::const_iterator i=theOuterComps.begin();
       i != theOuterComps.end(); i++){
    LogDebug("TkDetLayers") << "outer PixelRod pos z,perp,eta,phi: " 
                << (**i).position().z()    << " , " 
                << (**i).position().perp() << " , " 
                << (**i).position().eta()  << " , " 
                << (**i).position().phi()  ;
  }
  LogDebug("TkDetLayers") << "==== end DEBUG PixelBarrelLayer =====" ; 
  //----------------------------------- 
}

PixelBarrelLayer::~PixelBarrelLayer

(

)

Definition at line 79 of file PixelBarrelLayer.cc.

References i, and theComps.

                                   {
  vector<const GeometricSearchDet*>::const_iterator i;
  for (i=theComps.begin(); i!=theComps.end(); i++) {
    delete *i;
  }
} 

Member Function Documentation

bool PixelBarrelLayer::addClosest ( const TrajectoryStateOnSurface &  tsos, const Propagator &  prop, const MeasurementEstimator &  est, const SubLayerCrossing &  crossing, std::vector< DetGroup > &  result  ) const

private

Definition at line 169 of file PixelBarrelLayer.cc.

References CompatibleDetToGroupAdder::add(), SubLayerCrossing::closestDetIndex(), subLayer(), and SubLayerCrossing::subLayerIndex().

Referenced by groupedCompatibleDetsV().

{
  const vector<const GeometricSearchDet*>& sub( subLayer( crossing.subLayerIndex()));
  const GeometricSearchDet* det(sub[crossing.closestDetIndex()]);
  return CompatibleDetToGroupAdder::add( *det, tsos, prop, est, result);
}

virtual const std::vector<const GeomDet*>& PixelBarrelLayer::basicComponents ( ) const

inlinevirtual

Implements GeometricSearchDet.

Definition at line 29 of file PixelBarrelLayer.h.

References theBasicComps.

{return theBasicComps;}

double PixelBarrelLayer::calculatePhiWindow ( double  Xmax, const GeomDet &  det, const TrajectoryStateOnSurface &  state  ) const

private

Definition at line 190 of file PixelBarrelLayer.cc.

References TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::localPosition(), min, PV3DBase< T, PVType, FrameType >::phi(), edm::shift, GeomDet::surface(), and Surface::toGlobal().

Referenced by computeWindowSize().

{

  LocalPoint startPoint = state.localPosition();
  LocalVector shift( Xmax , 0. , 0.);
  LocalPoint shift1 = startPoint + shift;
  LocalPoint shift2 = startPoint + (-shift); 
  //LocalPoint shift2( startPoint); //original code;
  //shift2 -= shift;

  double phi1 = det.surface().toGlobal(shift1).phi();
  double phi2 = det.surface().toGlobal(shift2).phi();
  double phiStart = state.globalPosition().phi();
  double phiWin = min(fabs(phiStart-phi1),fabs(phiStart-phi2));

  return phiWin;
}

virtual const std::vector<const GeometricSearchDet*>& PixelBarrelLayer::components ( ) const

inlinevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 31 of file PixelBarrelLayer.h.

References theComps.

{return theComps;}

SubLayerCrossings PixelBarrelLayer::computeCrossings ( const TrajectoryStateOnSurface &  tsos, PropagationDirection  propDir  ) const

private

Definition at line 121 of file PixelBarrelLayer.cc.

References PeriodicBinFinderInPhi< T >::binIndex(), PeriodicBinFinderInPhi< T >::binPosition(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), HelixBarrelCylinderCrossing::hasSolution(), Geom::pi(), HelixBarrelCylinderCrossing::position(), rho, theInnerBinFinder, theInnerCylinder, theOuterBinFinder, theOuterCylinder, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by groupedCompatibleDetsV().

{
  GlobalPoint startPos( startingState.globalPosition());
  GlobalVector startDir( startingState.globalMomentum());
  double rho( startingState.transverseCurvature());
  
  HelixBarrelCylinderCrossing innerCrossing( startPos, startDir, rho,
                         propDir,*theInnerCylinder);

  if (!innerCrossing.hasSolution()) return SubLayerCrossings();
  //{
  //edm::LogInfo(TkDetLayers) << "ERROR in PixelBarrelLayer: inner cylinder not crossed by track" ;
  //throw DetLayerException("TkRodBarrelLayer: inner subRod not crossed by track");
  //}

  GlobalPoint gInnerPoint( innerCrossing.position());
  int innerIndex = theInnerBinFinder.binIndex(gInnerPoint.phi());
  float innerDist = theInnerBinFinder.binPosition(innerIndex) - gInnerPoint.phi();
  SubLayerCrossing innerSLC( 0, innerIndex, gInnerPoint);

  HelixBarrelCylinderCrossing outerCrossing( startPos, startDir, rho,
                         propDir,*theOuterCylinder);

  if (!outerCrossing.hasSolution()) return SubLayerCrossings();
  //if (!outerCrossing.hasSolution()) {
  //  throw DetLayerException("PixelBarrelLayer: inner cylinder not crossed by track");
  //}

  GlobalPoint gOuterPoint( outerCrossing.position());
  int outerIndex = theOuterBinFinder.binIndex(gOuterPoint.phi());
  float outerDist = theOuterBinFinder.binPosition(outerIndex) - gOuterPoint.phi() ;
  SubLayerCrossing outerSLC( 1, outerIndex, gOuterPoint);
  
  innerDist *= PhiLess()( theInnerBinFinder.binPosition(innerIndex),gInnerPoint.phi()) ? -1. : 1.; 
  outerDist *= PhiLess()( theOuterBinFinder.binPosition(outerIndex),gOuterPoint.phi()) ? -1. : 1.; 
  if (innerDist < 0.) { innerDist += 2.*Geom::pi();}
  if (outerDist < 0.) { outerDist += 2.*Geom::pi();}
  

  if (innerDist < outerDist) {
    return SubLayerCrossings( innerSLC, outerSLC, 0);
  }
  else {
    return SubLayerCrossings( outerSLC, innerSLC, 1);
  } 
}

float PixelBarrelLayer::computeWindowSize ( const GeomDet *  det, const TrajectoryStateOnSurface &  tsos, const MeasurementEstimator &  est  ) const

private

Definition at line 180 of file PixelBarrelLayer.cc.

References calculatePhiWindow(), MeasurementEstimator::maximalLocalDisplacement(), GeomDet::surface(), and vdt::x.

Referenced by groupedCompatibleDetsV().

{
  double xmax = 
    est.maximalLocalDisplacement(tsos, det->surface()).x();
  return calculatePhiWindow( xmax, *det, tsos);
}

BoundCylinder * PixelBarrelLayer::cylinder ( const std::vector< const GeometricSearchDet * > &  rods ) const

private

Definition at line 271 of file PixelBarrelLayer.cc.

References tmp.

Referenced by PixelBarrelLayer().

{
  vector<const GeomDet*> tmp;
  for (vector<const GeometricSearchDet*>::const_iterator it=rods.begin(); it!=rods.end(); it++) {    
    tmp.insert(tmp.end(),(*it)->basicComponents().begin(),(*it)->basicComponents().end());
  }
  return CylinderBuilderFromDet()( tmp.begin(), tmp.end());
}

void PixelBarrelLayer::groupedCompatibleDetsV ( const TrajectoryStateOnSurface &  tsos, const Propagator &  prop, const MeasurementEstimator &  est, std::vector< DetGroup > &  result  ) const

virtual

Reimplemented from GeometricSearchDet.

Definition at line 88 of file PixelBarrelLayer.cc.

References addClosest(), LayerCrossingSide::barrelSide(), SubLayerCrossings::closest(), SubLayerCrossings::closestIndex(), computeCrossings(), computeWindowSize(), SubLayerCrossings::isValid(), DetGroupMerger::orderAndMergeTwoLevels(), SubLayerCrossings::other(), Propagator::propagationDirection(), query::result, searchNeighbors(), and svgfig::window().

                                                           {
  vector<DetGroup> closestResult;
  SubLayerCrossings  crossings;
  crossings = computeCrossings( tsos, prop.propagationDirection());
  if(! crossings.isValid()) return;  
  
  addClosest( tsos, prop, est, crossings.closest(), closestResult);
  if (closestResult.empty()) {
    addClosest( tsos, prop, est, crossings.other(), result);
    return;
  }

  DetGroupElement closestGel( closestResult.front().front());
  float window = computeWindowSize( closestGel.det(), closestGel.trajectoryState(), est);
  
  searchNeighbors( tsos, prop, est, crossings.closest(), window,
           closestResult, false);
  
  vector<DetGroup> nextResult;
  searchNeighbors( tsos, prop, est, crossings.other(), window,
           nextResult, true);
  
  int crossingSide = LayerCrossingSide().barrelSide( closestGel.trajectoryState(), prop);
  DetGroupMerger::orderAndMergeTwoLevels( closestResult, nextResult, result, 
                      crossings.closestIndex(), crossingSide);
}

bool PixelBarrelLayer::overlap ( const GlobalPoint &  gpos, const GeometricSearchDet &  rod, float  phiWin  ) const

private

Definition at line 253 of file PixelBarrelLayer.cc.

References PV3DBase< T, PVType, FrameType >::phi(), BoundSurface::phiSpan(), rangesIntersect(), and GeometricSearchDet::surface().

Referenced by searchNeighbors().

{
  GlobalPoint crossPoint(gpos);

  // introduce offset (extrapolated point and true propagated point differ by 0.0003 - 0.00033, 
  // due to thickness of Rod of 1 cm) 
  const float phiOffset = 0.00034;  //...TOBE CHECKED LATER...
  phiWin += phiOffset;

  // detector phi range
  
  pair<float,float> phiRange(crossPoint.phi()-phiWin, crossPoint.phi()+phiWin);

  return rangesIntersect(phiRange, gsdet.surface().phiSpan(), PhiLess());

} 

void PixelBarrelLayer::searchNeighbors ( const TrajectoryStateOnSurface &  tsos, const Propagator &  prop, const MeasurementEstimator &  est, const SubLayerCrossing &  crossing, float  window, std::vector< DetGroup > &  result, bool  checkClosest  ) const

private

Definition at line 210 of file PixelBarrelLayer.cc.

References CompatibleDetToGroupAdder::add(), PeriodicBinFinderInPhi< T >::binIndex(), SubLayerCrossing::closestDetIndex(), overlap(), PV3DBase< T, PVType, FrameType >::phi(), SubLayerCrossing::position(), subLayer(), SubLayerCrossing::subLayerIndex(), theInnerBinFinder, and theOuterBinFinder.

Referenced by groupedCompatibleDetsV().

{
  GlobalPoint gCrossingPos = crossing.position();

  const vector<const GeometricSearchDet*>& sLayer( subLayer( crossing.subLayerIndex()));
 
  int closestIndex = crossing.closestDetIndex();
  int negStartIndex = closestIndex-1;
  int posStartIndex = closestIndex+1;

  if (checkClosest) { // must decide if the closest is on the neg or pos side
    if ( PhiLess()( gCrossingPos.phi(), sLayer[closestIndex]->position().phi())) {
      posStartIndex = closestIndex;
    }
    else {
      negStartIndex = closestIndex;
    }
  }

  const BinFinderType& binFinder = (crossing.subLayerIndex()==0 ? theInnerBinFinder : theOuterBinFinder);

  CompatibleDetToGroupAdder adder;
  int quarter = sLayer.size()/4;
  for (int idet=negStartIndex; idet >= negStartIndex - quarter; idet--) {
    const GeometricSearchDet* neighborRod = sLayer[binFinder.binIndex(idet)];
    if (!overlap( gCrossingPos, *neighborRod, window)) break;
    if (!adder.add( *neighborRod, tsos, prop, est, result)) break;
    // maybe also add shallow crossing angle test here???
  }
  for (int idet=posStartIndex; idet < posStartIndex + quarter; idet++) {
    const GeometricSearchDet* neighborRod = sLayer[binFinder.binIndex(idet)];
    if (!overlap( gCrossingPos, *neighborRod, window)) break;
    if (!adder.add( *neighborRod, tsos, prop, est, result)) break;
    // maybe also add shallow crossing angle test here???
  }
}

virtual SubDetector PixelBarrelLayer::subDetector ( ) const

inlinevirtual

The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap)

Implements DetLayer.

Definition at line 40 of file PixelBarrelLayer.h.

References GeomDetEnumerators::PixelBarrel.

{ return GeomDetEnumerators::PixelBarrel;}

const std::vector<const GeometricSearchDet*>& PixelBarrelLayer::subLayer ( int  ind ) const

inlineprivate

Definition at line 75 of file PixelBarrelLayer.h.

References theInnerComps, and theOuterComps.

Referenced by addClosest(), and searchNeighbors().

                                                                       {
    return (ind==0 ? theInnerComps : theOuterComps);}

Member Data Documentation

std::vector<const GeomDet*> PixelBarrelLayer::theBasicComps

private

Definition at line 85 of file PixelBarrelLayer.h.

Referenced by basicComponents(), and PixelBarrelLayer().

std::vector<const GeometricSearchDet*> PixelBarrelLayer::theComps

private

Definition at line 82 of file PixelBarrelLayer.h.

Referenced by components(), PixelBarrelLayer(), and ~PixelBarrelLayer().

BinFinderType PixelBarrelLayer::theInnerBinFinder

private

Definition at line 87 of file PixelBarrelLayer.h.

Referenced by computeCrossings(), PixelBarrelLayer(), and searchNeighbors().

std::vector<const GeometricSearchDet*> PixelBarrelLayer::theInnerComps

private

Definition at line 83 of file PixelBarrelLayer.h.

Referenced by PixelBarrelLayer(), and subLayer().

ReferenceCountingPointer<BoundCylinder> PixelBarrelLayer::theInnerCylinder

private

Definition at line 90 of file PixelBarrelLayer.h.

Referenced by computeCrossings(), and PixelBarrelLayer().

BinFinderType PixelBarrelLayer::theOuterBinFinder

private

Definition at line 88 of file PixelBarrelLayer.h.

Referenced by computeCrossings(), PixelBarrelLayer(), and searchNeighbors().

std::vector<const GeometricSearchDet*> PixelBarrelLayer::theOuterComps

private

Definition at line 84 of file PixelBarrelLayer.h.

Referenced by PixelBarrelLayer(), and subLayer().

ReferenceCountingPointer<BoundCylinder> PixelBarrelLayer::theOuterCylinder

private

Definition at line 91 of file PixelBarrelLayer.h.

Referenced by computeCrossings(), and PixelBarrelLayer().