#include <TOBRod.h>

Inheritance diagram for TOBRod:

Public Types
typedef PeriodicBinFinderInZ < float >	BinFinderType

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 std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &ts, const Propagator &, const MeasurementEstimator &) 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

	TOBRod (std::vector< const GeomDet * > &innerDets, std::vector< const GeomDet * > &outerDets)

	~TOBRod ()

Public Member Functions inherited from DetRod
virtual const BoundPlane &	specificSurface () const
	Return the rod surface as a BoundPlane. More...

virtual const BoundSurface &	surface () const
	The surface of the GeometricSearchDet. More...

virtual	~DetRod ()

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

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

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

bool	overlap (const GlobalPoint &gpos, const GeomDet &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 GeomDet * > &	subRod (int ind) const

Private Attributes
std::vector< const GeomDet * >	theDets

BinFinderType	theInnerBinFinder

std::vector< const GeomDet * >	theInnerDets

ReferenceCountingPointer < BoundPlane >	theInnerPlane

BinFinderType	theOuterBinFinder

std::vector< const GeomDet * >	theOuterDets

ReferenceCountingPointer < BoundPlane >	theOuterPlane

Additional Inherited Members
Protected Member Functions inherited from DetRod
void	setPlane (BoundPlane *plane)
	Set the rod's plane. More...

Protected Attributes inherited from GeometricSearchDet
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A concrete implementation for TOB Rod

Definition at line 16 of file TOBRod.h.

Member Typedef Documentation

typedef PeriodicBinFinderInZ<float> TOBRod::BinFinderType

Definition at line 18 of file TOBRod.h.

Constructor & Destructor Documentation

TOBRod::TOBRod	(	std::vector< const GeomDet * > &	innerDets,
		std::vector< const GeomDet * > &	outerDets
	)

Definition at line 28 of file TOBRod.cc.

References i, LogDebug, DetRod::setPlane(), python.multivaluedict::sort(), theDets, theInnerBinFinder, theInnerDets, theInnerPlane, theOuterBinFinder, theOuterDets, and theOuterPlane.

                                              :
   theInnerDets(innerDets),theOuterDets(outerDets)
 {
   theDets.assign(theInnerDets.begin(),theInnerDets.end());
   theDets.insert(theDets.end(),theOuterDets.begin(),theOuterDets.end());
 
 
   RodPlaneBuilderFromDet planeBuilder;
   setPlane( planeBuilder( theDets));
   theInnerPlane = planeBuilder( theInnerDets);
   theOuterPlane = planeBuilder( theOuterDets);
 
 
   sort(theDets.begin(),theDets.end(),DetZLess());
   sort(theInnerDets.begin(),theInnerDets.end(),DetZLess());
   sort(theOuterDets.begin(),theOuterDets.end(),DetZLess());
   theInnerBinFinder = BinFinderType(theInnerDets.begin(), theInnerDets.end());
   theOuterBinFinder = BinFinderType(theOuterDets.begin(), theOuterDets.end());
 
 
  
   LogDebug("TkDetLayers") << "==== DEBUG TOBRod =====" ; 
   for (vector<const GeomDet*>::const_iterator i=theInnerDets.begin();
        i != theInnerDets.end(); i++){
     LogDebug("TkDetLayers") << "inner TOBRod's Det pos z,perp,eta,phi: " 
                 << (**i).position().z() << " , " 
                 << (**i).position().perp() << " , " 
                 << (**i).position().eta() << " , " 
                 << (**i).position().phi() ;
   }
   
   for (vector<const GeomDet*>::const_iterator i=theOuterDets.begin();
        i != theOuterDets.end(); i++){
     LogDebug("TkDetLayers") << "outer TOBRod's Det pos z,perp,eta,phi: " 
                 << (**i).position().z() << " , " 
                 << (**i).position().perp() << " , " 
                 << (**i).position().eta() << " , " 
                 << (**i).position().phi() ;
   }
   LogDebug("TkDetLayers") << "==== end DEBUG TOBRod =====" ; 
   
 
 
 }

TOBRod::~TOBRod ( )

Definition at line 74 of file TOBRod.cc.

                {
   
 } 

Member Function Documentation

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

private

Definition at line 173 of file TOBRod.cc.

References CompatibleDetToGroupAdder::add(), SubLayerCrossing::closestDetIndex(), query::result, SubLayerCrossing::subLayerIndex(), and subRod().

Referenced by groupedCompatibleDetsV().

 {
 
   const vector<const GeomDet*>& sRod( subRod( crossing.subLayerIndex()));
   return CompatibleDetToGroupAdder::add( *sRod[crossing.closestDetIndex()], 
                      tsos, prop, est, result);
 }

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

inlinevirtual

Implements GeometricSearchDet.

Definition at line 26 of file TOBRod.h.

References theDets.

26 {return theDets;}

TOBRod::theDets

std::vector< const GeomDet * > theDets

Definition: TOBRod.h:74

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

                                               {
   edm::LogError("TkDetLayers") << "temporary dummy implementation of TOBRod::compatible()!!" ;
   return pair<bool,TrajectoryStateOnSurface>();
 }

const vector< const GeometricSearchDet * > & TOBRod::components ( ) const

virtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 80 of file TOBRod.cc.

                         {
   throw DetLayerException("TOBRod doesn't have GeometricSearchDet components");
 }

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

private

Definition at line 136 of file TOBRod.cc.

References PeriodicBinFinderInZ< T >::binIndex(), PeriodicBinFinderInZ< T >::binPosition(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), rho, theInnerBinFinder, theInnerPlane, theOuterBinFinder, theOuterPlane, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by groupedCompatibleDetsV().

 {
   GlobalPoint startPos( startingState.globalPosition());
   GlobalVector startDir( startingState.globalMomentum());
   double rho( startingState.transverseCurvature());
 
   HelixBarrelPlaneCrossingByCircle crossing( startPos, startDir, rho, propDir);
 
   pair<bool,double> innerPath = crossing.pathLength( *theInnerPlane);
   if (!innerPath.first) return SubLayerCrossings();
 
   GlobalPoint gInnerPoint( crossing.position(innerPath.second));
   int innerIndex = theInnerBinFinder.binIndex(gInnerPoint.z());
   float innerDist = fabs( theInnerBinFinder.binPosition(innerIndex) - gInnerPoint.z());
   SubLayerCrossing innerSLC( 0, innerIndex, gInnerPoint);
 
   pair<bool,double> outerPath = crossing.pathLength( *theOuterPlane);
   if (!outerPath.first) return SubLayerCrossings();
 
   GlobalPoint gOuterPoint( crossing.position(outerPath.second));
   int outerIndex = theOuterBinFinder.binIndex(gOuterPoint.z());
   float outerDist = fabs( theOuterBinFinder.binPosition(outerIndex) - gOuterPoint.z());
   SubLayerCrossing outerSLC( 1, outerIndex, gOuterPoint);
 
   if (innerDist < outerDist) {
     return SubLayerCrossings( innerSLC, outerSLC, 0);
   }
   else {
     return SubLayerCrossings( outerSLC, innerSLC, 1);
   } 
 }

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

private

Definition at line 186 of file TOBRod.cc.

References MeasurementEstimator::maximalLocalDisplacement(), GeomDet::surface(), and detailsBasic3DVector::y.

Referenced by groupedCompatibleDetsV().

 {
   return
     est.maximalLocalDisplacement(tsos, det->surface()).y();
 }

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

virtual

Reimplemented from GeometricSearchDet.

Definition at line 96 of file TOBRod.cc.

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

                                                       {
   
   SubLayerCrossings  crossings; 
   crossings = computeCrossings( tsos, prop.propagationDirection());
   if(! crossings.isValid()) return;
 
   vector<DetGroup> closestResult;
   addClosest( tsos, prop, est, crossings.closest(), closestResult);
   if (closestResult.empty()){
     vector<DetGroup> nextResult;
     addClosest( tsos, prop, est, crossings.other(), nextResult);
     if(nextResult.empty())    return;
 
     DetGroupElement nextGel( nextResult.front().front());  
     int crossingSide = LayerCrossingSide().barrelSide( nextGel.trajectoryState(), prop);
     DetGroupMerger::orderAndMergeTwoLevels( closestResult, nextResult, result, 
                        crossings.closestIndex(), crossingSide);   
   } else {
   
     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 TOBRod::overlap	(	const GlobalPoint &	gpos,
		const GeomDet &	rod,
		float	phiWin
	)		const

private

Definition at line 235 of file TOBRod.cc.

References BoundSurface::bounds(), Bounds::length(), GeomDet::surface(), GloballyPositioned< T >::toLocal(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by searchNeighbors().

 {
   // check if the z window around TSOS overlaps with the detector theDet (with a 1% margin added)
   
   //   const float tolerance = 0.1;
   const float relativeMargin = 1.01;
 
   LocalPoint localCrossPoint( det.surface().toLocal(crossPoint));
   //   if (fabs(localCrossPoint.z()) > tolerance) {
   //     edm::LogInfo(TkDetLayers) << "TOBRod::overlap calculation assumes point on surface, but it is off by "
   //     << localCrossPoint.z() ;
   //   }
 
   float localY = localCrossPoint.y();
   float detHalfLength = det.surface().bounds().length()/2.;
 
   //   edm::LogInfo(TkDetLayers) << "TOBRod::overlap: Det at " << det.position() << " hit at " << localY 
   //        << " Window " << window << " halflength "  << detHalfLength ;
   
   if ( ( fabs(localY)-window) < relativeMargin*detHalfLength ) { // FIXME: margin hard-wired!
     return true;
   } else {
     return false;
   }
 }

void TOBRod::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 197 of file TOBRod.cc.

References Clusterizer1DCommons::add(), SubLayerCrossing::closestDetIndex(), overlap(), SubLayerCrossing::position(), SubLayerCrossing::subLayerIndex(), subRod(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by groupedCompatibleDetsV().

 {
   GlobalPoint gCrossingPos = crossing.position();
 
   const vector<const GeomDet*>& sRod( subRod( 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 (gCrossingPos.z() < sRod[closestIndex]->surface().position().z()) {
       posStartIndex = closestIndex;
     }
     else {
       negStartIndex = closestIndex;
     }
   }
 
   typedef CompatibleDetToGroupAdder Adder;
   for (int idet=negStartIndex; idet >= 0; idet--) {
     if (!overlap( gCrossingPos, *sRod[idet], window)) break;
     if (!Adder::add( *sRod[idet], tsos, prop, est, result)) break;
   }
   for (int idet=posStartIndex; idet < static_cast<int>(sRod.size()); idet++) {
     if (!overlap( gCrossingPos, *sRod[idet], window)) break;
     if (!Adder::add( *sRod[idet], tsos, prop, est, result)) break;
   }
 }

const std::vector<const GeomDet*>& TOBRod::subRod ( int ind ) const

inlineprivate

Definition at line 68 of file TOBRod.h.

References theInnerDets, and theOuterDets.

Referenced by addClosest(), and searchNeighbors().

                                                           {
     return (ind==0 ? theInnerDets : theOuterDets);
   }

Member Data Documentation

std::vector<const GeomDet*> TOBRod::theDets

private

Definition at line 74 of file TOBRod.h.

Referenced by basicComponents(), and TOBRod().

BinFinderType TOBRod::theInnerBinFinder

private

Definition at line 81 of file TOBRod.h.

Referenced by computeCrossings(), and TOBRod().

std::vector<const GeomDet*> TOBRod::theInnerDets

private

Definition at line 75 of file TOBRod.h.

Referenced by subRod(), and TOBRod().

ReferenceCountingPointer<BoundPlane> TOBRod::theInnerPlane

private

Definition at line 78 of file TOBRod.h.

Referenced by computeCrossings(), and TOBRod().

BinFinderType TOBRod::theOuterBinFinder

private

Definition at line 82 of file TOBRod.h.

Referenced by computeCrossings(), and TOBRod().

std::vector<const GeomDet*> TOBRod::theOuterDets

private

Definition at line 76 of file TOBRod.h.

Referenced by subRod(), and TOBRod().

ReferenceCountingPointer<BoundPlane> TOBRod::theOuterPlane

private

Definition at line 79 of file TOBRod.h.

Referenced by computeCrossings(), and TOBRod().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation