#include <CompositeTECPetal.h>

Inheritance diagram for CompositeTECPetal:

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...

	CompositeTECPetal (std::vector< const TECWedge * > &innerWedges, std::vector< const TECWedge * > &outerWedges)

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

	~CompositeTECPetal ()

Public Member Functions inherited from TECPetal
virtual const BoundDiskSector &	specificSurface () const

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

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

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 ()

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

int	findBin (float R, int layer) const

GlobalPoint	findPosition (int index, int diskSectorIndex) 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

Static Private Member Functions
static float	computeWindowSize (const GeomDet *det, const TrajectoryStateOnSurface &tsos, const MeasurementEstimator &est)

static bool	overlap (const GlobalPoint &gpos, const GeometricSearchDet &rod, float window)

Private Attributes
std::vector< float >	theBackBoundaries

std::vector< const GeometricSearchDet * >	theBackComps

ReferenceCountingPointer < BoundDiskSector >	theBackSector

std::vector< const GeomDet * >	theBasicComps

std::vector< const GeometricSearchDet * >	theComps

std::vector< float >	theFrontBoundaries

std::vector< const GeometricSearchDet * >	theFrontComps

ReferenceCountingPointer < BoundDiskSector >	theFrontSector

Additional Inherited Members
Public Types inherited from GeometricSearchDetWithGroups
typedef GeometricSearchDet::DetWithState	DetWithState

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 TECPetal
ReferenceCountingPointer < BoundDiskSector >	theDiskSector

Protected Attributes inherited from GeometricSearchDet
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A concrete implementation for TEC petals

Definition at line 17 of file CompositeTECPetal.h.

Constructor & Destructor Documentation

CompositeTECPetal::CompositeTECPetal	(	std::vector< const TECWedge * > &	innerWedges,
		std::vector< const TECWedge * > &	outerWedges
	)

Definition at line 58 of file CompositeTECPetal.cc.

References LogDebug, theBackBoundaries, theBackComps, theBackSector, theBasicComps, theComps, TECPetal::theDiskSector, theFrontBoundaries, theFrontComps, and theFrontSector.

                                                            : 
   theFrontComps(innerWedges.begin(),innerWedges.end()), 
   theBackComps(outerWedges.begin(),outerWedges.end())
 {
   theComps.assign(theFrontComps.begin(),theFrontComps.end());
   theComps.insert(theComps.end(),theBackComps.begin(),theBackComps.end());
 
   details::fillBoundaries( theFrontComps, theFrontBoundaries);
   details::fillBoundaries( theBackComps, theBackBoundaries);
 
 
   for(vector<const GeometricSearchDet*>::const_iterator it=theComps.begin();
       it!=theComps.end();it++){  
     theBasicComps.insert(theBasicComps.end(),   
              (**it).basicComponents().begin(),
              (**it).basicComponents().end());
   }
 
 
   //the Wedge are already R ordered
   //sort( theWedges.begin(), theWedges.end(), DetLessR());
   //sort( theFrontWedges.begin(), theFrontWedges.end(), DetLessR() );
   //sort( theBackWedges.begin(), theBackWedges.end(), DetLessR() );
   vector<const TECWedge*> allWedges;
   allWedges.assign(innerWedges.begin(),innerWedges.end());
   allWedges.insert(allWedges.end(),outerWedges.begin(),outerWedges.end());
 
   theDiskSector  = ForwardDiskSectorBuilderFromWedges()( allWedges );
   theFrontSector = ForwardDiskSectorBuilderFromWedges()( innerWedges);
   theBackSector  = ForwardDiskSectorBuilderFromWedges()( outerWedges);
 
   //--------- DEBUG INFO --------------
   LogDebug("TkDetLayers") << "DEBUG INFO for CompositeTECPetal" ;
 
   for(vector<const GeometricSearchDet*>::const_iterator it=theFrontComps.begin(); 
       it!=theFrontComps.end(); it++){
     LogDebug("TkDetLayers") << "frontWedge phi,z,r: " 
                 << (*it)->surface().position().phi() << " , "
                 << (*it)->surface().position().z() <<   " , "
                 << (*it)->surface().position().perp() ;
   }
 
   for(vector<const GeometricSearchDet*>::const_iterator it=theBackComps.begin(); 
       it!=theBackComps.end(); it++){
     LogDebug("TkDetLayers") << "backWedge phi,z,r: " 
                 << (*it)->surface().position().phi() << " , "
                 << (*it)->surface().position().z() <<   " , "
                 << (*it)->surface().position().perp() ;
   }
   //----------------------------------- 
 
 
 }

CompositeTECPetal::~CompositeTECPetal ( )

Definition at line 114 of file CompositeTECPetal.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 CompositeTECPetal::addClosest	(	const TrajectoryStateOnSurface &	tsos,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		const SubLayerCrossing &	crossing,
		std::vector< DetGroup > &	result
	)		const

private

Definition at line 226 of file CompositeTECPetal.cc.

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

Referenced by groupedCompatibleDetsV().

 {
   const vector<const GeometricSearchDet*>& sub( subLayer( crossing.subLayerIndex()));
   const GeometricSearchDet* det(sub[crossing.closestDetIndex()]);
 
   LogDebug("TkDetLayers") 
     << "in TECPetal, adding Wedge at r,z,phi: (" 
     << det->position().perp() << "," 
     << det->position().z() << "," 
     << det->position().phi() << ")" ;
   LogDebug("TkDetLayers") 
     << "wedge comps size: " 
     << det->basicComponents().size();
 
   return CompatibleDetToGroupAdder::add( *det, tsos, prop, est, result);
 }

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

inlinevirtual

Implements GeometricSearchDet.

Definition at line 25 of file CompositeTECPetal.h.

References theBasicComps.

25 {return theBasicComps;}

CompositeTECPetal::theBasicComps

std::vector< const GeomDet * > theBasicComps

Definition: CompositeTECPetal.h:82

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

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

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

inlinevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 27 of file CompositeTECPetal.h.

References theComps.

27 {return theComps;}

CompositeTECPetal::theComps

std::vector< const GeometricSearchDet * > theComps

Definition: CompositeTECPetal.h:79

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

private

Definition at line 173 of file CompositeTECPetal.cc.

References findBin(), findPosition(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), LogDebug, perp(), rho, theBackSector, theFrontSector, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by groupedCompatibleDetsV().

 {
   double rho( startingState.transverseCurvature());
   
   HelixPlaneCrossing::PositionType startPos( startingState.globalPosition() );
   HelixPlaneCrossing::DirectionType startDir( startingState.globalMomentum() );
   HelixForwardPlaneCrossing crossing(startPos,startDir,rho,propDir);
   pair<bool,double> frontPath = crossing.pathLength( *theFrontSector);
 
   if (!frontPath.first) return SubLayerCrossings();
 
   GlobalPoint gFrontPoint(crossing.position(frontPath.second));
   LogDebug("TkDetLayers") 
     << "in TECPetal,front crossing : r,z,phi: (" 
     << gFrontPoint.perp() << ","
     << gFrontPoint.z() << "," 
     << gFrontPoint.phi() << ")";
   
 
   int frontIndex = findBin(gFrontPoint.perp(),0);
   float frontDist = fabs( findPosition(frontIndex,0).perp() - gFrontPoint.perp());
   SubLayerCrossing frontSLC( 0, frontIndex, gFrontPoint);
 
 
 
   pair<bool,double> backPath = crossing.pathLength( *theBackSector);
 
   if (!backPath.first) return SubLayerCrossings();
   
 
   GlobalPoint gBackPoint( crossing.position(backPath.second));
   LogDebug("TkDetLayers") 
     << "in TECPetal,back crossing r,z,phi: (" 
     << gBackPoint.perp() << ","
     << gBackPoint.z() << "," 
     << gBackPoint.phi() << ")" ;
 
   int backIndex = findBin(gBackPoint.perp(),1);
   float backDist = fabs( findPosition(backIndex,1).perp() - gBackPoint.perp());
   
   SubLayerCrossing backSLC( 1, backIndex, gBackPoint);
   
   
   // 0ss: frontDisk has index=0, backDisk has index=1
   if (frontDist < backDist) {
     return SubLayerCrossings( frontSLC, backSLC, 0);
   }
   else {
     return SubLayerCrossings( backSLC, frontSLC, 1);
   } 
 }

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

staticprivate

Definition at line 323 of file CompositeTECPetal.cc.

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

Referenced by groupedCompatibleDetsV().

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

int CompositeTECPetal::findBin	(	float	R,
		int	layer
	)		const

private

Definition at line 331 of file CompositeTECPetal.cc.

References theBackBoundaries, and theFrontBoundaries.

Referenced by computeCrossings().

 {
   return details::findBin(diskSectorType==0 ? theFrontBoundaries : theBackBoundaries,R);
 }

GlobalPoint CompositeTECPetal::findPosition	(	int	index,
		int	diskSectorIndex
	)		const

private

Definition at line 339 of file CompositeTECPetal.cc.

References theBackComps, and theFrontComps.

Referenced by computeCrossings().

 {
   vector<const GeometricSearchDet*> const & diskSector = diskSectorType == 0 ? theFrontComps : theBackComps; 
   return (diskSector[index])->position();
 }

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

virtual

Reimplemented from GeometricSearchDet.

Definition at line 131 of file CompositeTECPetal.cc.

References addClosest(), SubLayerCrossings::closest(), SubLayerCrossings::closestIndex(), computeCrossings(), computeWindowSize(), LayerCrossingSide::endcapSide(), SubLayerCrossings::isValid(), LogDebug, DetGroupMerger::orderAndMergeTwoLevels(), SubLayerCrossings::other(), Propagator::propagationDirection(), 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); 
   LogDebug("TkDetLayers") << "in TECPetal, closestResult.size(): "<< closestResult.size();
 
   if (closestResult.empty()){
     vector<DetGroup> nextResult;
     addClosest( tsos, prop, est, crossings.other(), nextResult); 
     LogDebug("TkDetLayers") << "in TECPetal, nextResult.size(): "<< nextResult.size() ;
     if(nextResult.empty())    return;
     
     DetGroupElement nextGel( nextResult.front().front());  
     int crossingSide = LayerCrossingSide().endcapSide( 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().endcapSide( closestGel.trajectoryState(), prop);
   DetGroupMerger::orderAndMergeTwoLevels( closestResult, nextResult, result, 
                           crossings.closestIndex(), crossingSide);
   }
 }

bool CompositeTECPetal::overlap	(	const GlobalPoint &	gpos,
		const GeometricSearchDet &	rod,
		float	window
	)

staticprivate

Definition at line 297 of file CompositeTECPetal.cc.

References BoundDiskSector::bounds(), BoundDiskSector::innerRadius(), max(), BoundDiskSector::outerRadius(), PV3DBase< T, PVType, FrameType >::perp(), GloballyPositioned< T >::position(), GeometricSearchDet::surface(), DiskSectorBounds::thickness(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by searchNeighbors().

 {
   // this method is just a duplication of overlapInR 
   // adapeted for groupedCompatibleDets() needs
 
   // assume "fixed theta window", i.e. margin in local y = r is changing linearly with z
   float tsRadius = gpos.perp();
   float thetamin = ( max((float)0.,tsRadius-ymax))/(fabs(gpos.z())+10.); // add 10 cm contingency 
   float thetamax = ( tsRadius + ymax)/(fabs(gpos.z())-10.);
 
   const BoundDiskSector& wedgeSector = static_cast<const BoundDiskSector&>( gsdet.surface());                   
   float wedgeMinZ = fabs( wedgeSector.position().z()) - 0.5*wedgeSector.bounds().thickness();
   float wedgeMaxZ = fabs( wedgeSector.position().z()) + 0.5*wedgeSector.bounds().thickness(); 
   float thetaWedgeMin =  wedgeSector.innerRadius()/ wedgeMaxZ;
   float thetaWedgeMax =  wedgeSector.outerRadius()/ wedgeMinZ;
   
   // do the theta regions overlap ?
 
   return  !( thetamin > thetaWedgeMax || thetaWedgeMin > thetamax);
   
 } 

void CompositeTECPetal::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 250 of file CompositeTECPetal.cc.

References Clusterizer1DCommons::add(), SubLayerCrossing::closestDetIndex(), overlap(), PV3DBase< T, PVType, FrameType >::perp(), SubLayerCrossing::position(), subLayer(), SubLayerCrossing::subLayerIndex(), theBackComps, and theFrontComps.

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 ( gCrossingPos.perp() < sLayer[closestIndex]->position().perp() ) {
       posStartIndex = closestIndex;
     }
     else {
       negStartIndex = closestIndex;
     }
   }
 
 
   //const BinFinderType& binFinder = (crossing.subLayerIndex()==0 ? theFrontBinFinder : theBackBinFinder);
   int theSize = crossing.subLayerIndex()==0 ? theFrontComps.size() : theBackComps.size();
   
   typedef CompatibleDetToGroupAdder Adder;
   for (int idet=negStartIndex; idet >= 0; idet--) {
     //if(idet<0 || idet>= theSize) {edm::LogInfo(TkDetLayers) << "===== error! gone out vector bounds.idet: " << idet ;exit;}
     const GeometricSearchDet & neighborWedge = *sLayer[idet];
     if (!overlap( gCrossingPos, neighborWedge, window)) break;  // --- to check
     if (!Adder::add( neighborWedge, tsos, prop, est, result)) break;
     // maybe also add shallow crossing angle test here???
   }
   for (int idet=posStartIndex; idet <theSize; idet++) {
     //if(idet<0 || idet>= theSize) {edm::LogInfo(TkDetLayers) << "===== error! gone out vector bounds.idet: " << idet ;exit;}
     const GeometricSearchDet & neighborWedge = *sLayer[idet];
     if (!overlap( gCrossingPos, neighborWedge, window)) break;  // ---- to check
     if (!Adder::add( neighborWedge, tsos, prop, est, result)) break;
     // maybe also add shallow crossing angle test here???
   }
 }

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

inlineprivate

Definition at line 73 of file CompositeTECPetal.h.

References theBackComps, and theFrontComps.

Referenced by addClosest(), and searchNeighbors().

                                                                                      {
     return (ind==0 ? theFrontComps : theBackComps);
   }

Member Data Documentation

std::vector<float> CompositeTECPetal::theBackBoundaries

private

Definition at line 85 of file CompositeTECPetal.h.

Referenced by CompositeTECPetal(), and findBin().

std::vector<const GeometricSearchDet*> CompositeTECPetal::theBackComps

private

Definition at line 81 of file CompositeTECPetal.h.

Referenced by CompositeTECPetal(), findPosition(), searchNeighbors(), and subLayer().

ReferenceCountingPointer<BoundDiskSector> CompositeTECPetal::theBackSector

private

Definition at line 88 of file CompositeTECPetal.h.

Referenced by CompositeTECPetal(), and computeCrossings().

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

private

Definition at line 82 of file CompositeTECPetal.h.

Referenced by basicComponents(), and CompositeTECPetal().

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

private

Definition at line 79 of file CompositeTECPetal.h.

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

std::vector<float> CompositeTECPetal::theFrontBoundaries

private

Definition at line 84 of file CompositeTECPetal.h.

Referenced by CompositeTECPetal(), and findBin().

std::vector<const GeometricSearchDet*> CompositeTECPetal::theFrontComps

private

Definition at line 80 of file CompositeTECPetal.h.

Referenced by CompositeTECPetal(), findPosition(), searchNeighbors(), and subLayer().

ReferenceCountingPointer<BoundDiskSector> CompositeTECPetal::theFrontSector

private

Definition at line 87 of file CompositeTECPetal.h.

Referenced by CompositeTECPetal(), and computeCrossings().

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation