#include <DirectTrackerNavigation.h>

Public Member Functions
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection) const
	find compatible layers for a given trajectory state More...

	DirectTrackerNavigation (const edm::ESHandle< GeometricSearchTracker > &, bool outOnly=true)
	constructor More...

Private Member Functions
float	calculateEta (float r, float z) const

bool	checkCompatible (const FreeTrajectoryState &, const BarrelDetLayer *) const

bool	checkCompatible (const FreeTrajectoryState &, const ForwardDetLayer *) const

void	inOutBPx (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutBTEC (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutBTID (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutFPx (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutFTEC (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutFTID (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutPx (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutTIB (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

void	inOutTOB (const FreeTrajectoryState &, std::vector< const DetLayer * > &) const

bool	outward (const FreeTrajectoryState &) const

Private Attributes
float	theEpsilon

edm::ESHandle< GeometricSearchTracker >	theGeometricSearchTracker

bool	theOutLayerOnlyFlag

Detailed Description

Definition at line 26 of file DirectTrackerNavigation.h.

Constructor & Destructor Documentation

DirectTrackerNavigation::DirectTrackerNavigation	(	const edm::ESHandle< GeometricSearchTracker > &	tkLayout,
		bool	outOnly = `true`
	)

constructor

Definition at line 32 of file DirectTrackerNavigation.cc.

                                                                : 
  theGeometricSearchTracker(tkLayout), theOutLayerOnlyFlag(outOnly), theEpsilon(-0.01) {
 
 }

Member Function Documentation

float DirectTrackerNavigation::calculateEta	(	float	r,
		float	z
	)		const

private

Definition at line 244 of file DirectTrackerNavigation.cc.

References cmsBatch::log, and funct::tan().

Referenced by checkCompatible().

                                                                   {
 
   if ( z > 0 ) return -log((tan(atan(r/z)/2.)));
   return log(-(tan(atan(r/z)/2.)));
 
 }

bool DirectTrackerNavigation::checkCompatible	(	const FreeTrajectoryState &	fts,
		const BarrelDetLayer *	dl
	)		const

private

Definition at line 194 of file DirectTrackerNavigation.cc.

References calculateEta(), stringResolutionProvider_cfi::eta, PV3DBase< T, PVType, FrameType >::eta(), FreeTrajectoryState::position(), TCMET_cfi::radius, BarrelDetLayer::specificSurface(), and theEpsilon.

Referenced by inOutBPx(), inOutBTEC(), inOutBTID(), inOutFPx(), inOutFTEC(), inOutFTID(), inOutPx(), inOutTIB(), and inOutTOB().

                                                                                                            {
 
   float eta0 = fts.position().eta();
 
   const BoundCylinder& bc = dl->specificSurface();
   float radius = bc.radius();
   float length = bc.bounds().length()/2.;
 
   float eta = calculateEta(radius, length);
 
   return ( fabs(eta0) <= (fabs(eta) + theEpsilon) );
 
 }

bool DirectTrackerNavigation::checkCompatible	(	const FreeTrajectoryState &	fts,
		const ForwardDetLayer *	dl
	)		const

private

Definition at line 212 of file DirectTrackerNavigation.cc.

References calculateEta(), PV3DBase< T, PVType, FrameType >::eta(), FreeTrajectoryState::position(), ForwardDetLayer::specificSurface(), theEpsilon, and z.

                                                                                                             {
 
   float eta0 = fts.position().eta();
 
   const BoundDisk& bd = dl->specificSurface();
 
   float outRadius = bd.outerRadius();
   float inRadius = bd.innerRadius();
   float z = bd.position().z();
 
   float etaOut = calculateEta(outRadius, z);
   float etaIn = calculateEta(inRadius, z);
  
   if ( eta0 > 0 ) return ( eta0 > ( etaOut - theEpsilon) && eta0 < (etaIn + theEpsilon) );
   else return ( eta0 < (etaOut + theEpsilon ) && eta0 > ( etaIn - theEpsilon ) );
 
 }

vector< const DetLayer * > DirectTrackerNavigation::compatibleLayers	(	const FreeTrajectoryState &	fts,
		PropagationDirection	timeDirection
	)		const

find compatible layers for a given trajectory state

Definition at line 43 of file DirectTrackerNavigation.cc.

References PV3DBase< T, PVType, FrameType >::eta(), inOutBPx(), inOutBTEC(), inOutBTID(), inOutFPx(), inOutFTEC(), inOutFTID(), inOutPx(), inOutTIB(), inOutTOB(), LogTrace, oppositeToMomentum, convertSQLitetoXML_cfg::output, outward(), FreeTrajectoryState::position(), groupFilesInBlocks::reverse, and theOutLayerOnlyFlag.

Referenced by TSGFromPropagation::trackerSeeds().

                                                                           {
 
   bool inOut = outward(fts);
   double eta0 = fts.position().eta();
 
   vector<const DetLayer*> output;
 
   // check eta of DetLayers for compatibility
 
   if (inOut) { 
 
       if ( !theOutLayerOnlyFlag ) {
          inOutPx(fts,output);
 
          if ( eta0 > 1.55) inOutFPx(fts,output);
          else if ( eta0 < -1.55 ) inOutBPx(fts,output);
 
          if ( fabs(eta0) < 1.67 ) inOutTIB(fts,output); 
 
          if ( eta0 > 1.17 ) inOutFTID(fts,output);
          else if ( eta0 < -1.17 ) inOutBTID(fts,output);
       }
 
       if ( fabs(eta0) < 1.35 ) inOutTOB(fts,output);
 
       if ( eta0 > 0.97 ) inOutFTEC(fts,output);
       else if ( eta0 < -0.97 )  inOutBTEC(fts,output);
 
    } else {
       LogTrace("Muon|RecoMuon|DirectionTrackerNavigation")<<"No implementation for inward state at this moment. ";
    }
 
   if ( dir == oppositeToMomentum ) std::reverse(output.begin(),output.end());
 
   return output;
 
 }

void DirectTrackerNavigation::inOutBPx	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 158 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::negPixelForwardLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                             {
 
   for(const auto i: theGeometricSearchTracker->negPixelForwardLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutBTEC	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 182 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::negTecLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                              {
 
   for(const auto i: theGeometricSearchTracker->negTecLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutBTID	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 170 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::negTidLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                              {
 
   for(const auto i: theGeometricSearchTracker->negTidLayers()) {
       if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutFPx	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 122 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::posPixelForwardLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                             {
 
   for(const auto i : theGeometricSearchTracker->posPixelForwardLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutFTEC	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 146 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::posTecLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                              {
 
   for(const auto i : theGeometricSearchTracker->posTecLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutFTID	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 134 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::posTidLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                              {
 
   for(const auto i: theGeometricSearchTracker->posTidLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutPx	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 86 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, GeometricSearchTracker::pixelBarrelLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

                                                                                                            {
 
   for (const auto i : theGeometricSearchTracker->pixelBarrelLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutTIB	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 98 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, theGeometricSearchTracker, and GeometricSearchTracker::tibLayers().

Referenced by compatibleLayers().

                                                                                                             {
 
   for(const auto i : theGeometricSearchTracker->tibLayers() ) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

void DirectTrackerNavigation::inOutTOB	(	const FreeTrajectoryState &	fts,
		std::vector< const DetLayer * > &	output
	)		const

private

Definition at line 110 of file DirectTrackerNavigation.cc.

References checkCompatible(), mps_fire::i, theGeometricSearchTracker, and GeometricSearchTracker::tobLayers().

Referenced by compatibleLayers().

                                                                                                             {
 
   for( const auto i : theGeometricSearchTracker->tobLayers()) {
     if ( checkCompatible(fts,i) ) output.push_back(i);
   }
 
 }

bool DirectTrackerNavigation::outward ( const FreeTrajectoryState & fts ) const

private

Definition at line 234 of file DirectTrackerNavigation.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), Basic3DVector< T >::dot(), FreeTrajectoryState::momentum(), and FreeTrajectoryState::position().

Referenced by compatibleLayers().

                                                                           {
  
   return (fts.position().basicVector().dot(fts.momentum().basicVector()) > 0);
 
 }