#include <TECLayer.h>

Inheritance diagram for TECLayer:

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

virtual SubDetector	subDetector () const
	The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap) More...

	TECLayer (std::vector< const TECPetal * > &innerPetals, std::vector< const TECPetal * > &outerPetals)

	~TECLayer ()

Public Member Functions inherited from ForwardDetLayer
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	~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

Protected Types
typedef PeriodicBinFinderInPhi < double >	BinFinderPhi

Protected Member Functions
virtual BoundDisk *	computeDisk (std::vector< const GeometricSearchDet * > &petals) const

Protected Member Functions inherited from ForwardDetLayer
virtual BoundDisk *	computeSurface ()

virtual void	initialize ()

float	rmax () const

float	rmin () const

void	setSurface (BoundDisk *cp)

float	zmax () const

float	zmin () const

Protected Attributes
BinFinderPhi	theBackBinFinder

std::vector< const GeometricSearchDet * >	theBackComps

ReferenceCountingPointer < BoundDisk >	theBackDisk

std::vector< const GeomDet * >	theBasicComps

std::vector< const GeometricSearchDet * >	theComps

BinFinderPhi	theFrontBinFinder

std::vector< const GeometricSearchDet * >	theFrontComps

ReferenceCountingPointer < BoundDisk >	theFrontDisk

Protected Attributes inherited from GeometricSearchDet
GeomDetCompatibilityChecker	theCompatibilityChecker

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 &startingState, PropagationDirection propDir) const

bool	overlap (const GlobalPoint &gpos, const GeometricSearchDet &petal, float window) 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

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

Public Types inherited from GeometricSearchDetWithGroups
typedef GeometricSearchDet::DetWithState	DetWithState

Detailed Description

A concrete implementation for TEC layer built out of TECPetals

Definition at line 16 of file TECLayer.h.

Member Typedef Documentation

typedef PeriodicBinFinderInPhi<double> TECLayer::BinFinderPhi

protected

Definition at line 77 of file TECLayer.h.

Constructor & Destructor Documentation

TECLayer::TECLayer	(	std::vector< const TECPetal * > &	innerPetals,
		std::vector< const TECPetal * > &	outerPetals
	)

Definition at line 21 of file TECLayer.cc.

References computeDisk(), LogDebug, PV3DBase< T, PVType, FrameType >::perp(), GeometricSearchDet::position(), ForwardDetLayer::setSurface(), ForwardDetLayer::specificSurface(), theBackBinFinder, theBackComps, theBackDisk, theBasicComps, theComps, theFrontBinFinder, theFrontComps, theFrontDisk, and PV3DBase< T, PVType, FrameType >::z().

                                                  : 
   theFrontComps(innerPetals.begin(),innerPetals.end()), 
   theBackComps(outerPetals.begin(),outerPetals.end())
 {
   theComps.assign(theFrontComps.begin(),theFrontComps.end());
   theComps.insert(theComps.end(),theBackComps.begin(),theBackComps.end());
 
   for(vector<const GeometricSearchDet*>::const_iterator it=theComps.begin();
       it!=theComps.end();it++){  
     theBasicComps.insert(theBasicComps.end(),   
              (**it).basicComponents().begin(),
              (**it).basicComponents().end());
   }
 
 
   //This should be no necessary. TO BE CHECKED
   //sort(theFrontPetals.begin(), theFrontPetals.end(), PetalLessPhi());
   //sort(theBackPetals.begin(), theBackPetals.end(), PetalLessPhi());
 
   // building disk for front and back petals
   setSurface( computeDisk( theComps ) );
   theFrontDisk = computeDisk( theFrontComps);
   theBackDisk  = computeDisk( theBackComps);
 
   // set up the bin finders
   theFrontBinFinder = BinFinderPhi(theFrontComps.front()->position().phi(),
                    theFrontComps.size());
   theBackBinFinder  = BinFinderPhi(theBackComps.front()->position().phi(),
                    theBackComps.size());  
 
   //--------- DEBUG INFO --------------
   LogDebug("TkDetLayers") << "DEBUG INFO for TECLayer" << "\n"
               << "TECLayer z,perp, innerRadius, outerR: " 
               << this->position().z()    << " , "
               << this->position().perp() << " , "
               << this->specificSurface().innerRadius() << " , "
               << this->specificSurface().outerRadius() ;
   
 
   for(vector<const GeometricSearchDet*>::const_iterator it=theFrontComps.begin(); 
       it!=theFrontComps.end(); it++){
     LogDebug("TkDetLayers") << "frontPetal 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") << "backPetal phi,z,r: " 
      << (*it)->surface().position().phi() << " , "
      << (*it)->surface().position().z() <<   " , "
      << (*it)->surface().position().perp() ;
   }
   //----------------------------------- 
 
 
 }

TECLayer::~TECLayer ( )

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

private

Definition at line 189 of file TECLayer.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 TECLayer, adding petal at r,z,phi: (" 
     << det->position().perp() << "," 
     << det->position().z() << "," 
     << det->position().phi() << ")" << endl;
 
   return CompatibleDetToGroupAdder().add( *det, tsos, prop, est, result); 
 }

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

inlinevirtual

Implements GeometricSearchDet.

Definition at line 24 of file TECLayer.h.

References theBasicComps.

24 {return theBasicComps;}

TECLayer::theBasicComps

std::vector< const GeomDet * > theBasicComps

Definition: TECLayer.h:71

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

inlinevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 26 of file TECLayer.h.

References theComps.

26 {return theComps;}

TECLayer::theComps

std::vector< const GeometricSearchDet * > theComps

Definition: TECLayer.h:68

SubLayerCrossings TECLayer::computeCrossings	(	const TrajectoryStateOnSurface &	startingState,
		PropagationDirection	propDir
	)		const

private

Definition at line 132 of file TECLayer.cc.

References funct::abs(), PeriodicBinFinderInPhi< T >::binIndex(), Geom::deltaPhi(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), LogDebug, GloballyPositioned< T >::phi(), rho, ForwardDetLayer::surface(), theBackBinFinder, theBackComps, theBackDisk, theFrontBinFinder, theFrontComps, theFrontDisk, 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( *theFrontDisk);
   if (!frontPath.first) SubLayerCrossings();
 
 
   GlobalPoint gFrontPoint(crossing.position(frontPath.second));
   LogDebug("TkDetLayers") 
     << "in TECLayer,front crossing point: r,z,phi: (" 
     << gFrontPoint.perp() << ","
     << gFrontPoint.z() << "," 
     << gFrontPoint.phi() << ")" << endl;
   
 
   int frontIndex = theFrontBinFinder.binIndex(gFrontPoint.barePhi()); 
   SubLayerCrossing frontSLC( 0, frontIndex, gFrontPoint);
 
 
 
   pair<bool,double> backPath = crossing.pathLength( *theBackDisk);
   if (!backPath.first) SubLayerCrossings();
 
 
   GlobalPoint gBackPoint( crossing.position(backPath.second));
   LogDebug("TkDetLayers") 
     << "in TECLayer,back crossing point: r,z,phi: (" 
     << gBackPoint.perp() << "," 
     << gFrontPoint.z() << "," 
     << gBackPoint.phi() << ")" << endl;
 
 
   int backIndex = theBackBinFinder.binIndex(gBackPoint.barePhi());
   SubLayerCrossing backSLC( 1, backIndex, gBackPoint);
 
   
   // 0ss: frontDisk has index=0, backDisk has index=1
   float frontDist = std::abs(Geom::deltaPhi( double(gFrontPoint.barePhi()), 
                          double(theFrontComps[frontIndex]->surface().phi())));
   float backDist = std::abs(Geom::deltaPhi( double(gBackPoint.barePhi()), 
                         double(theBackComps[backIndex]->surface().phi())));
   
 
   if (frontDist < backDist) {
     return SubLayerCrossings( frontSLC, backSLC, 0);
   }
   else {
     return SubLayerCrossings( backSLC, frontSLC, 1);
   } 
 }

BoundDisk * TECLayer::computeDisk ( std::vector< const GeometricSearchDet * > & petals ) const

protectedvirtual

Definition at line 265 of file TECLayer.cc.

References i, max(), bookConverter::min, ForwardDetLayer::rmax(), ForwardDetLayer::rmin(), makeMuonMisalignmentScenario::rot, ForwardDetLayer::zmax(), ForwardDetLayer::zmin(), and zPos.

Referenced by TECLayer().

 {
   // Attention: it is assumed that the petals do belong to one layer, and are all
   // of the same rmin/rmax extension !!  
   
   const BoundDiskSector&  diskSector = static_cast<const BoundDiskSector&>(petals.front()->surface());
 
   float rmin = diskSector.innerRadius();
   float rmax = diskSector.outerRadius();
   
   float theZmax(petals.front()->position().z());
   float theZmin(theZmax);
   for ( vector<const GeometricSearchDet*>::const_iterator i = petals.begin(); i != petals.end(); i++ ) {
     float zmin = (**i).position().z() - (**i).surface().bounds().thickness()/2.;
     float zmax = (**i).position().z() + (**i).surface().bounds().thickness()/2.;
     theZmax = max( theZmax, zmax);
     theZmin = min( theZmin, zmin);
   }
 
   float zPos = (theZmax+theZmin)/2.;
   PositionType pos(0.,0.,zPos);
   RotationType rot;
 
   return new BoundDisk( pos, rot, new SimpleDiskBounds(rmin, rmax,    
                                    theZmin-zPos, theZmax-zPos));
 }

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

virtual

Reimplemented from GeometricSearchDet.

Definition at line 92 of file TECLayer.cc.

References addClosest(), SubLayerCrossings::closest(), SubLayerCrossings::closestIndex(), computeCrossings(), tkDetUtil::computeWindowSize(), LayerCrossingSide::endcapSide(), SubLayerCrossings::isValid(), LogDebug, DetGroupMerger::orderAndMergeTwoLevels(), SubLayerCrossings::other(), Propagator::propagationDirection(), and searchNeighbors().

                                                            {  
   SubLayerCrossings  crossings; 
   crossings = computeCrossings( tsos, prop.propagationDirection());
   if(! crossings.isValid()) return;
 
   vector<DetGroup> closestResult;
   addClosest( tsos, prop, est, crossings.closest(), closestResult); 
   LogDebug("TkDetLayers") << "in TECLayer, closestResult.size(): " << closestResult.size();
 
   // this differs from other groupedCompatibleDets logic, which DON'T check next in such cases!!!
   if(closestResult.empty()){
     vector<DetGroup> nextResult;
     addClosest( tsos, prop, est, crossings.other(), nextResult);   
     LogDebug("TkDetLayers") << "in TECLayer, nextResult.size(): " << nextResult.size();
     if(nextResult.empty())       return;
     
 
     DetGroupElement nextGel( nextResult.front().front());  
     int crossingSide = LayerCrossingSide().endcapSide( nextGel.trajectoryState(), prop);
     DetGroupMerger::orderAndMergeTwoLevels( std::move(closestResult), std::move(nextResult), result, 
                         crossings.closestIndex(), crossingSide);   
   }  
   else {
     DetGroupElement closestGel( closestResult.front().front());  
     float phiWindow = tkDetUtil::computeWindowSize( closestGel.det(), closestGel.trajectoryState(), est); 
     searchNeighbors( tsos, prop, est, crossings.closest(), phiWindow,
              closestResult, false); 
     vector<DetGroup> nextResult;  
     searchNeighbors( tsos, prop, est, crossings.other(), phiWindow,
              nextResult, true); 
     
     int crossingSide = LayerCrossingSide().endcapSide( closestGel.trajectoryState(), prop);
     DetGroupMerger::orderAndMergeTwoLevels( std::move(closestResult), std::move(nextResult), result,
                         crossings.closestIndex(), crossingSide);
   }
 }

bool TECLayer::overlap	(	const GlobalPoint &	gpos,
		const GeometricSearchDet &	petal,
		float	window
	)		const

private

Definition at line 250 of file TECLayer.cc.

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

Referenced by searchNeighbors().

 {
   float phi = gpos.barePhi();
   const BoundDiskSector&  diskSector = static_cast<const BoundDiskSector&>(gsdet.surface());
   pair<float,float> phiRange(phi-phiWin,phi+phiWin);
   pair<float,float> petalPhiRange(diskSector.phi() - 0.5*diskSector.phiExtension(),
                   diskSector.phi() + 0.5*diskSector.phiExtension());
 
 
   return rangesIntersect(phiRange, petalPhiRange, PhiLess());
 }

void TECLayer::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 207 of file TECLayer.cc.

References Clusterizer1DCommons::add(), PeriodicBinFinderInPhi< T >::binIndex(), SubLayerCrossing::closestDetIndex(), overlap(), Geom::Phi< T >::phi(), PV3DBase< T, PVType, FrameType >::phi(), SubLayerCrossing::position(), subLayer(), SubLayerCrossing::subLayerIndex(), theBackBinFinder, and theFrontBinFinder.

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 BinFinderPhi& binFinder = (crossing.subLayerIndex()==0 ? theFrontBinFinder : theBackBinFinder);
 
   typedef CompatibleDetToGroupAdder Adder;
   int half = sLayer.size()/2;  // to check if dets are called twice....
   for (int idet=negStartIndex; idet >= negStartIndex - half; idet--) {
     const GeometricSearchDet & neighborPetal = *sLayer[binFinder.binIndex(idet)];
     if (!overlap( gCrossingPos, neighborPetal, window)) break;
     if (!Adder::add( neighborPetal, tsos, prop, est, result)) break;
     // maybe also add shallow crossing angle test here???
   }
   for (int idet=posStartIndex; idet < posStartIndex + half; idet++) {
     const GeometricSearchDet & neighborPetal = *sLayer[binFinder.binIndex(idet)];
     if (!overlap( gCrossingPos, neighborPetal, window)) break;
     if (!Adder::add( neighborPetal, tsos, prop, est, result)) break;
     // maybe also add shallow crossing angle test here???
   }
 }

virtual SubDetector TECLayer::subDetector ( ) const

inlinevirtual

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

Implements DetLayer.

Definition at line 34 of file TECLayer.h.

References GeomDetEnumerators::TEC.

34 {return GeomDetEnumerators::TEC;}

GeomDetEnumerators::TEC

Definition: GeomDetEnumerators.h:11

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

inlineprivate

Definition at line 60 of file TECLayer.h.

References theBackComps, and theFrontComps.

Referenced by addClosest(), and searchNeighbors().

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

Member Data Documentation

BinFinderPhi TECLayer::theBackBinFinder

protected

Definition at line 80 of file TECLayer.h.

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

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

protected

Definition at line 70 of file TECLayer.h.

Referenced by computeCrossings(), subLayer(), and TECLayer().

ReferenceCountingPointer<BoundDisk> TECLayer::theBackDisk

protected

Definition at line 75 of file TECLayer.h.

Referenced by computeCrossings(), and TECLayer().

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

protected

Definition at line 71 of file TECLayer.h.

Referenced by basicComponents(), and TECLayer().

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

protected

Definition at line 68 of file TECLayer.h.

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

BinFinderPhi TECLayer::theFrontBinFinder

protected

Definition at line 79 of file TECLayer.h.

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

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

protected

Definition at line 69 of file TECLayer.h.

Referenced by computeCrossings(), subLayer(), and TECLayer().

ReferenceCountingPointer<BoundDisk> TECLayer::theFrontDisk

protected

Definition at line 74 of file TECLayer.h.

Referenced by computeCrossings(), and TECLayer().

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation