CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
DirectTrackerNavigation Class Reference

#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< GeometricSearchTrackertheGeometricSearchTracker
 
bool theOutLayerOnlyFlag
 

Detailed Description

Definition at line 26 of file DirectTrackerNavigation.h.

Constructor & Destructor Documentation

◆ DirectTrackerNavigation()

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

constructor

Definition at line 32 of file DirectTrackerNavigation.cc.

33  : theGeometricSearchTracker(tkLayout), theOutLayerOnlyFlag(outOnly), theEpsilon(-0.01) {}
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker

Member Function Documentation

◆ calculateEta()

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

Definition at line 219 of file DirectTrackerNavigation.cc.

References dqm-mbProfile::log, alignCSCRings::r, funct::tan(), and z.

Referenced by checkCompatible().

219  {
220  if (z > 0)
221  return -log((tan(atan(r / z) / 2.)));
222  return log(-(tan(atan(r / z) / 2.)));
223 }
Tan< T >::type tan(const T &t)
Definition: Tan.h:22

◆ checkCompatible() [1/2]

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

Definition at line 176 of file DirectTrackerNavigation.cc.

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

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

176  {
177  float eta0 = fts.position().eta();
178 
179  const BoundCylinder& bc = dl->specificSurface();
180  float radius = bc.radius();
181  float length = bc.bounds().length() / 2.;
182 
183  float eta = calculateEta(radius, length);
184 
185  return (fabs(eta0) <= (fabs(eta) + theEpsilon));
186 }
float calculateEta(float r, float z) const
T eta() const
Definition: PV3DBase.h:73
GlobalPoint position() const
virtual const BoundCylinder & specificSurface() const final
Extension of the interface.

◆ checkCompatible() [2/2]

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

Definition at line 191 of file DirectTrackerNavigation.cc.

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

191  {
192  float eta0 = fts.position().eta();
193 
194  const BoundDisk& bd = dl->specificSurface();
195 
196  float outRadius = bd.outerRadius();
197  float inRadius = bd.innerRadius();
198  float z = bd.position().z();
199 
200  float etaOut = calculateEta(outRadius, z);
201  float etaIn = calculateEta(inRadius, z);
202 
203  if (eta0 > 0)
204  return (eta0 > (etaOut - theEpsilon) && eta0 < (etaIn + theEpsilon));
205  else
206  return (eta0 < (etaOut + theEpsilon) && eta0 > (etaIn - theEpsilon));
207 }
float calculateEta(float r, float z) const
T eta() const
Definition: PV3DBase.h:73
GlobalPoint position() const
virtual const BoundDisk & specificSurface() const final

◆ compatibleLayers()

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

find compatible layers for a given trajectory state

Definition at line 38 of file DirectTrackerNavigation.cc.

References DeadROC_duringRun::dir, 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.

39  {
40  bool inOut = outward(fts);
41  double eta0 = fts.position().eta();
42 
43  vector<const DetLayer*> output;
44 
45  // check eta of DetLayers for compatibility
46 
47  if (inOut) {
48  if (!theOutLayerOnlyFlag) {
49  inOutPx(fts, output);
50 
51  if (eta0 > 1.55)
52  inOutFPx(fts, output);
53  else if (eta0 < -1.55)
54  inOutBPx(fts, output);
55 
56  if (fabs(eta0) < 1.67)
57  inOutTIB(fts, output);
58 
59  if (eta0 > 1.17)
60  inOutFTID(fts, output);
61  else if (eta0 < -1.17)
62  inOutBTID(fts, output);
63  }
64 
65  if (fabs(eta0) < 1.35)
66  inOutTOB(fts, output);
67 
68  if (eta0 > 0.97)
69  inOutFTEC(fts, output);
70  else if (eta0 < -0.97)
71  inOutBTEC(fts, output);
72 
73  } else {
74  LogTrace("Muon|RecoMuon|DirectionTrackerNavigation") << "No implementation for inward state at this moment. ";
75  }
76 
77  if (dir == oppositeToMomentum)
78  std::reverse(output.begin(), output.end());
79 
80  return output;
81 }
T eta() const
Definition: PV3DBase.h:73
void inOutBTEC(const FreeTrajectoryState &, std::vector< const DetLayer *> &) const
void inOutFTID(const FreeTrajectoryState &, std::vector< const DetLayer *> &) const
#define LogTrace(id)
GlobalPoint position() const
void inOutBPx(const FreeTrajectoryState &, std::vector< const DetLayer *> &) const
bool outward(const FreeTrajectoryState &) const
void inOutFTEC(const FreeTrajectoryState &, std::vector< const DetLayer *> &) const
void inOutTOB(const FreeTrajectoryState &, std::vector< const DetLayer *> &) const
void inOutPx(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 inOutTIB(const FreeTrajectoryState &, std::vector< const DetLayer *> &) const

◆ inOutBPx()

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

Definition at line 146 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

146  {
147  for (const auto i : theGeometricSearchTracker->negPixelForwardLayers()) {
148  if (checkCompatible(fts, i))
149  output.push_back(i);
150  }
151 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< ForwardDetLayer const * > const & negPixelForwardLayers() const

◆ inOutBTEC()

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

Definition at line 166 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

166  {
167  for (const auto i : theGeometricSearchTracker->negTecLayers()) {
168  if (checkCompatible(fts, i))
169  output.push_back(i);
170  }
171 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< ForwardDetLayer const * > const & negTecLayers() const

◆ inOutBTID()

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

Definition at line 156 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

156  {
157  for (const auto i : theGeometricSearchTracker->negTidLayers()) {
158  if (checkCompatible(fts, i))
159  output.push_back(i);
160  }
161 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< ForwardDetLayer const * > const & negTidLayers() const

◆ inOutFPx()

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

Definition at line 116 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

116  {
117  for (const auto i : theGeometricSearchTracker->posPixelForwardLayers()) {
118  if (checkCompatible(fts, i))
119  output.push_back(i);
120  }
121 }
std::vector< ForwardDetLayer const * > const & posPixelForwardLayers() const
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker

◆ inOutFTEC()

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

Definition at line 136 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

136  {
137  for (const auto i : theGeometricSearchTracker->posTecLayers()) {
138  if (checkCompatible(fts, i))
139  output.push_back(i);
140  }
141 }
std::vector< ForwardDetLayer const * > const & posTecLayers() const
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker

◆ inOutFTID()

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

Definition at line 126 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

126  {
127  for (const auto i : theGeometricSearchTracker->posTidLayers()) {
128  if (checkCompatible(fts, i))
129  output.push_back(i);
130  }
131 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< ForwardDetLayer const * > const & posTidLayers() const

◆ inOutPx()

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, convertSQLitetoXML_cfg::output, GeometricSearchTracker::pixelBarrelLayers(), and theGeometricSearchTracker.

Referenced by compatibleLayers().

86  {
87  for (const auto i : theGeometricSearchTracker->pixelBarrelLayers()) {
88  if (checkCompatible(fts, i))
89  output.push_back(i);
90  }
91 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< BarrelDetLayer const * > const & pixelBarrelLayers() const

◆ inOutTIB()

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

Definition at line 96 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

96  {
97  for (const auto i : theGeometricSearchTracker->tibLayers()) {
98  if (checkCompatible(fts, i))
99  output.push_back(i);
100  }
101 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< BarrelDetLayer const * > const & tibLayers() const

◆ inOutTOB()

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

Definition at line 106 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

106  {
107  for (const auto i : theGeometricSearchTracker->tobLayers()) {
108  if (checkCompatible(fts, i))
109  output.push_back(i);
110  }
111 }
bool checkCompatible(const FreeTrajectoryState &, const BarrelDetLayer *) const
edm::ESHandle< GeometricSearchTracker > theGeometricSearchTracker
std::vector< BarrelDetLayer const * > const & tobLayers() const

◆ outward()

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

Definition at line 212 of file DirectTrackerNavigation.cc.

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

Referenced by compatibleLayers().

212  {
213  return (fts.position().basicVector().dot(fts.momentum().basicVector()) > 0);
214 }
GlobalPoint position() const
GlobalVector momentum() const
const BasicVectorType & basicVector() const
Definition: PV3DBase.h:53
T dot(const Basic3DVector &rh) const
Scalar product, or "dot" product, with a vector of same type.

Member Data Documentation

◆ theEpsilon

float DirectTrackerNavigation::theEpsilon
private

Definition at line 67 of file DirectTrackerNavigation.h.

Referenced by checkCompatible().

◆ theGeometricSearchTracker

edm::ESHandle<GeometricSearchTracker> DirectTrackerNavigation::theGeometricSearchTracker
private

◆ theOutLayerOnlyFlag

bool DirectTrackerNavigation::theOutLayerOnlyFlag
private

Definition at line 65 of file DirectTrackerNavigation.h.

Referenced by compatibleLayers().