#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 __attribute__((cold))
	Returns basic components, if any. More...

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

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

	~TECLayer () __attribute__((cold))

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

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
	DetLayer (bool doHaveGroup, bool ibar)

bool	isBarrel () const

bool	isForward () const

int	seqNum () const

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

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

	GeometricSearchDet (bool doHaveGroups)

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

bool	hasGroups () const

virtual const Surface::PositionType &	position () const
	Returns position of the surface. More...

virtual	~GeometricSearchDet ()

Protected Types
typedef PeriodicBinFinderInPhi < float >	BinFinderPhi

Protected Attributes
BinFinderPhi	theBackBinFinder

std::vector< const TECPetal * >	theBackComps

ReferenceCountingPointer < BoundDisk >	theBackDisk

std::vector< const GeomDet * >	theBasicComps

std::vector< const GeometricSearchDet * >	theComps

BinFinderPhi	theFrontBinFinder

std::vector< const TECPetal * >	theFrontComps

ReferenceCountingPointer < BoundDisk >	theFrontDisk

Protected Attributes inherited from GeometricSearchDet
bool	haveGroups

GeomDetCompatibilityChecker	theCompatibilityChecker

Private Member Functions
bool	addClosest (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, std::vector< DetGroup > &result) const __attribute__((hot))

SubLayerCrossings	computeCrossings (const TrajectoryStateOnSurface &startingState, PropagationDirection propDir) const __attribute__((hot))

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

const std::vector< const TECPetal * > &	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

Protected Member Functions inherited from ForwardDetLayer
SimpleDiskBounds const &	bounds () const

virtual BoundDisk *	computeSurface ()

virtual void	initialize ()

float	rmax () const

float	rmin () const

void	setSurface (BoundDisk *cp)

float	zmax () const

float	zmin () const

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<float> TECLayer::BinFinderPhi

protected

Definition at line 78 of file TECLayer.h.

Constructor & Destructor Documentation

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

Definition at line 56 of file TECLayer.cc.

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

                                                  : 
   ForwardDetLayer(true),
   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(auto 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(auto 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 119 of file TECLayer.cc.

References i, and theComps.

                    {
   for (auto 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 225 of file TECLayer.cc.

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

Referenced by groupedCompatibleDetsV().

 {
   const auto & sub( subLayer( crossing.subLayerIndex()));
   const auto  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:69

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 167 of file TECLayer.cc.

References funct::abs(), PeriodicBinFinderInPhi< T >::binIndex(), SiPixelRawToDigiRegional_cfi::deltaPhi, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), LogDebug, rho, 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( gFrontPoint.barePhi(), 
                          theFrontComps[frontIndex]->surface().phi()) );
   float backDist = std::abs(Geom::deltaPhi( gBackPoint.barePhi(), 
                         theBackComps[backIndex]->surface().phi()) );
   
 
   if (frontDist < backDist) {
     return SubLayerCrossings( frontSLC, backSLC, 0);
   }
   else {
     return SubLayerCrossings( backSLC, frontSLC, 1);
   } 
 }

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

virtual

Reimplemented from GeometricSearchDet.

Definition at line 127 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);
   }
 }

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 259 of file TECLayer.cc.

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

Referenced by groupedCompatibleDetsV().

 {
   GlobalPoint gCrossingPos = crossing.position();
   auto gphi = gCrossingPos.barePhi();  
 
   const auto & 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()( gphi, sLayer[closestIndex]->surface().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 auto & neighborPetal = *sLayer[binFinder.binIndex(idet)];
     if (!overlap( gphi, 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 auto & neighborPetal = *sLayer[binFinder.binIndex(idet)];
     if (!overlap( gphi, 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 TECPetal*>& TECLayer::subLayer ( int ind ) const

inlineprivate

Definition at line 61 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 81 of file TECLayer.h.

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

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

protected

Definition at line 72 of file TECLayer.h.

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

ReferenceCountingPointer<BoundDisk> TECLayer::theBackDisk

protected

Definition at line 76 of file TECLayer.h.

Referenced by computeCrossings(), and TECLayer().

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

protected

Definition at line 69 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 80 of file TECLayer.h.

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

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

protected

Definition at line 71 of file TECLayer.h.

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

ReferenceCountingPointer<BoundDisk> TECLayer::theFrontDisk

protected

Definition at line 75 of file TECLayer.h.

Referenced by computeCrossings(), and TECLayer().

Public Member Functions

Protected Types

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation