MuonNavigableLayer Class Reference

#include <MuonNavigableLayer.h>

Inheritance diagram for MuonNavigableLayer:

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...
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
void	setSchool (NavigationSchool const *sh)
virtual	~NavigableLayer ()

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

Detailed Description

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

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.

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

Member Function Documentation

compatibleLayers() [1/2]

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

overridepure virtual

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

compatibleLayers() [2/2]

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

overridepure virtual

NavigableLayer interface.

Reimplemented from NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

detLayer()

const DetLayer* MuonNavigableLayer::detLayer ( ) const

overridepure virtual

return DetLayer

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

isInsideOut()

bool MuonNavigableLayer::isInsideOut

(

const FreeTrajectoryState &

fts

)

const

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

nextLayers() [1/2]

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

overridepure virtual

NavigableLayer interface.

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

nextLayers() [2/2]

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

overridepure virtual

NavigableLayer interface.

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

setDetLayer()

void MuonNavigableLayer::setDetLayer ( const DetLayer *  )

overridepure virtual

set DetLayer

Implements NavigableLayer.

Implemented in MuonBarrelNavigableLayer, and MuonForwardNavigableLayer.

trackingRange()

MuonEtaRange MuonNavigableLayer::trackingRange

(

const FreeTrajectoryState &

fts

)

const

Definition at line 36 of file MuonNavigableLayer.cc.

References calculateEta(), FreeTrajectoryState::curvilinearError(), PVValHelper::eta, egammaIdentification::eta_max, dqm-mbProfile::log, CurvilinearTrajectoryError::matrix(), SiStripPI::min, PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), FastTimerService_cff::range, 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));

  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;
}