MuonNavigableLayer Class Reference

#include <MuonNavigableLayer.h>

Inheritance diagram for MuonNavigableLayer:

Public Member Functions
virtual std::vector< const DetLayer * >	compatibleLayers (NavigationDirection dir) const =0
virtual std::vector< const DetLayer * >	compatibleLayers (const FreeTrajectoryState &fts, PropagationDirection dir) const =0
	NavigableLayer interface. More...
virtual DetLayer *	detLayer () const =0
	return DetLayer More...
bool	isInsideOut (const FreeTrajectoryState &fts) const
virtual std::vector< const DetLayer * >	nextLayers (NavigationDirection dir) const =0
	NavigableLayer interface. More...
virtual std::vector< const DetLayer * >	nextLayers (const FreeTrajectoryState &fts, PropagationDirection dir) const =0
	NavigableLayer interface. More...
virtual void	setDetLayer (DetLayer *)=0
	set DetLayer More...
MuonEtaRange	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	~NavigableLayer ()

Detailed Description

base class for MuonBarrelNavigableLayer and MuonForwardNavigable. trackingRange defines an MuonEtaRange for an FTS, which is used for search compatible DetLayers.

Date:

2007/01/29 16:24:17

Revision:

1.8

Author

: Chang Liu - Purdue University Chang.nosp@m..Liu.nosp@m.@cern.nosp@m..ch

Modification:

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

Date:

2007/01/29 16:25:21

Revision:

1.9

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 28 of file MuonNavigableLayer.h.

Member Function Documentation

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

pure virtual

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

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

pure virtual

NavigableLayer interface.

Reimplemented from NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

virtual DetLayer* MuonNavigableLayer::detLayer ( ) const

pure virtual

return DetLayer

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

bool MuonNavigableLayer::isInsideOut

(

const FreeTrajectoryState &

fts

)

const

Definition at line 74 of file MuonNavigableLayer.cc.

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

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

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

pure virtual

NavigableLayer interface.

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

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

pure virtual

NavigableLayer interface.

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

virtual void MuonNavigableLayer::setDetLayer ( DetLayer *  )

pure virtual

set DetLayer

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

MuonEtaRange MuonNavigableLayer::trackingRange

(

const FreeTrajectoryState &

fts

)

const

Definition at line 39 of file MuonNavigableLayer.cc.

References calculateEta(), FreeTrajectoryState::curvilinearError(), eta(), create_public_lumi_plots::log, CurvilinearTrajectoryError::matrix(), min, PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), alignCSCRings::r, funct::sin(), mathSSE::sqrt(), funct::tan(), theta(), detailsBasic3DVector::z, 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));

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

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

  return range;
}