MTDNavigableLayer Class Reference

#include <MTDNavigableLayer.h>

Inheritance diagram for MTDNavigableLayer:

Public Member Functions
std::vector< const DetLayer * >	compatibleLayers (NavigationDirection dir) const override=0
std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection dir) const override=0
	NavigableLayer interface. More...
const DetLayer *	detLayer () const override=0
	return DetLayer More...
bool	isInsideOut (const FreeTrajectoryState &fts) const
std::vector< const DetLayer * >	nextLayers (NavigationDirection dir) const override=0
	NavigableLayer interface. More...
std::vector< const DetLayer * >	nextLayers (const FreeTrajectoryState &fts, PropagationDirection dir) const override=0
	NavigableLayer interface. More...
void	setDetLayer (const DetLayer *) override=0
	set DetLayer More...
MTDEtaRange	trackingRange (const FreeTrajectoryState &fts) const
Public Member Functions inherited from NavigableLayer
virtual std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection timeDirection, int &counter) const
virtual void	setDetLayer (DetLayer const *dl)=0
void	setSchool (NavigationSchool const *sh)
virtual	~NavigableLayer ()

Additional Inherited Members
Protected Attributes inherited from NavigableLayer
NavigationSchool const *	school = 0

Detailed Description

base class for BTLNavigableLayer and ETLNavigableLayer. trackingRange defines an MTDEtaRange for an FTS, which is used for search compatible DetLayers.

Author

: L. Gray - FNAL

Modification:

base class for MTDNavigableLayer and MuonForwardNavigableLayer. trackingRange defines an MTDEtaRange for an FTS, which is used for searching compatible DetLayers.

Author

: Chang Liu - Purdue University Chang.nosp@m..Liu.nosp@m.@cern.nosp@m..ch with contributions from: R. Bellan - INFN Torino

code of trackingRange is from MuonGlobalNavigation in ORCA whose author is Stefano Lacaprara - INFN Padova Modification:

Definition at line 26 of file MTDNavigableLayer.h.

Member Function Documentation

std::vector<const DetLayer*> MTDNavigableLayer::compatibleLayers ( NavigationDirection  dir ) const

overridepure virtual

Implements NavigableLayer.

Implemented in BTLNavigableLayer, and ETLNavigableLayer.

std::vector<const DetLayer*> MTDNavigableLayer::compatibleLayers ( const FreeTrajectoryState &  fts, PropagationDirection  dir  ) const

overridepure virtual

NavigableLayer interface.

Reimplemented from NavigableLayer.

Implemented in BTLNavigableLayer, and ETLNavigableLayer.

const DetLayer* MTDNavigableLayer::detLayer ( ) const

overridepure virtual

return DetLayer

Implements NavigableLayer.

Implemented in BTLNavigableLayer, and ETLNavigableLayer.

bool MTDNavigableLayer::isInsideOut

(

const FreeTrajectoryState &

fts

)

const

Definition at line 72 of file MTDNavigableLayer.cc.

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

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

std::vector<const DetLayer*> MTDNavigableLayer::nextLayers ( NavigationDirection  dir ) const

overridepure virtual

NavigableLayer interface.

Implements NavigableLayer.

Implemented in BTLNavigableLayer, and ETLNavigableLayer.

std::vector<const DetLayer*> MTDNavigableLayer::nextLayers ( const FreeTrajectoryState &  fts, PropagationDirection  dir  ) const

overridepure virtual

NavigableLayer interface.

Implements NavigableLayer.

Implemented in BTLNavigableLayer, and ETLNavigableLayer.

void MTDNavigableLayer::setDetLayer ( const DetLayer *  )

overridepure virtual

set DetLayer

Implemented in BTLNavigableLayer, and ETLNavigableLayer.

MTDEtaRange MTDNavigableLayer::trackingRange

(

const FreeTrajectoryState &

fts

)

const

Definition at line 37 of file MTDNavigableLayer.cc.

References calculateEta(), FreeTrajectoryState::curvilinearError(), PVValHelper::eta, MuonTCMETValueMapProducer_cff::eta_max, cmsBatch::log, CurvilinearTrajectoryError::matrix(), min(), PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), alignCSCRings::r, funct::sin(), mathSSE::sqrt(), funct::tan(), theta(), and PV3DBase< T, PVType, FrameType >::z().

{  
  float z = fts.position().z();
  float r = fts.position().perp();
  float eta;
  if ( z>0 ) eta = -log((tan(atan(r/z)/2.)));
  else eta = log(-(tan(atan(r/z)/2.)));

  double theta = atan(r/z);

  double spread = 5.0*sqrt(fts.curvilinearError().matrix()(2,2))/fabs(sin(theta));  //5*sigma(eta)

  //C.L.: this spread could be too large to use.
  // convert it to a smaller one by assuming a virtual radius
  // that transforms the error on angle to error on z axis.
  // not accurate, but works!

  double eta_max = 0;

  if ( z > 0 ) eta_max = calculateEta(r, z+spread); 
  else eta_max = calculateEta(r, z-spread); 

  spread = std::min(0.07, fabs(eta_max-eta));

  MTDEtaRange range(eta+spread,eta-spread);

  spread = 0.07; 
  // special treatment for special geometry in overlap region
  
  if ( eta > 1.0 && eta < 1.1 )  range = MTDEtaRange(eta+3.0*spread,eta-spread);
  if ( eta < -1.0 && eta > -1.1 ) range = MTDEtaRange(eta+spread,eta-3.0*spread);

  return range;
}