ATLAS Offline Software
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
Muon::MuonTGMeasurementTool Class Reference

#include <MuonTGMeasurementTool.h>

Inheritance diagram for Muon::MuonTGMeasurementTool:
Collaboration diagram for Muon::MuonTGMeasurementTool:

Public Types

typedef Trk::Track Track
 
typedef std::pair< const Trk::Layer *, std::vector< Identifier > * > PairOfLayerID
 
typedef std::pair< const Trk::Layer *, std::vector< const Trk::PrepRawData * > * > PairOfLayerPrd
 

Public Member Functions

 MuonTGMeasurementTool (const std::string &type, const std::string &name, const IInterface *)
 Constructor with AlgTool parameters. More...
 
virtual ~MuonTGMeasurementTool ()=default
 
virtual StatusCode initialize () override
 
const Trk::TrackParameters * layerToDetEl (const Trk::Layer *, const Trk::TrackParameters *, Identifier) const override
 
const Trk::TrackParameters * detElToLayer (const Trk::Layer *, const Trk::TrackParameters *, Identifier) const override
 
const Trk::RIO_OnTrack * measToLayer (const Trk::Layer *, const Trk::TrackParameters *, const Trk::RIO_OnTrack *) const override
 
double residual (const Trk::Layer *, const Trk::TrackParameters *, const Trk::RIO_OnTrack *) const override
 
double residual (const Trk::Layer *, const Trk::TrackParameters *, Identifier) const override
 
double residual (const Trk::TrackParameters *, const Trk::RIO_OnTrack *) const override
 
double residual (const Trk::TrackParameters *, Identifier &) const override
 
const Identifier nearestDetEl (const Trk::Layer *, const Trk::TrackParameters *, bool measPhi, double &pitch) const override
 
const Trk::Layer * associatedLayer (Identifier id, Amg::Vector3D &gp) const override
 
const Trk::Layer * associatedLayer (Identifier id, const Trk::TrackingVolume *vol) const override
 
const Trk::Layer * match (Identifier id, const Trk::Layer *lay) const override
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Static Public Member Functions

static const InterfaceID & interfaceID ()
 AlgTool interface method. More...
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

const Trk::TrackingGeometry * getGeometry () const
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
 
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager > m_DetectorManagerKey
 
ServiceHandle< Trk::ITrackingGeometrySvc > m_trackingGeometrySvc
 
SG::ReadCondHandleKey< Trk::TrackingGeometry > m_trackingGeometryReadKey
 
const MuonGM::MuonDetectorManager * m_muonDetMgr
 
std::unique_ptr< AmgMatrix(5, 5)> m_tgcProjEta
 
std::unique_ptr< AmgMatrix(5, 5)> m_tgcProjPhi
 
std::unique_ptr< AmgMatrix(5, 5)> m_rpcProjEta
 
std::unique_ptr< AmgMatrix(5, 5)> m_rpcProjPhi
 
Gaudi::Property< bool > m_alignedMode {this, "AlignedMode", true}
 
Gaudi::Property< bool > m_useDSManager {this, "UseDSManager", false}
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

contain application tools for (muon) tracking geometry

Author
Sarka.nosp@m..Tod.nosp@m.orova.nosp@m.@cer.nosp@m.n.ch

Definition at line 49 of file MuonTGMeasurementTool.h.

Member Typedef Documentation

◆ PairOfLayerID

typedef std::pair<const Trk::Layer*, std::vector<Identifier>*> Muon::MuonTGMeasurementTool::PairOfLayerID

Definition at line 52 of file MuonTGMeasurementTool.h.

◆ PairOfLayerPrd

typedef std::pair<const Trk::Layer*, std::vector<const Trk::PrepRawData*>*> Muon::MuonTGMeasurementTool::PairOfLayerPrd

Definition at line 53 of file MuonTGMeasurementTool.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ Track

typedef Trk::Track Muon::MuonTGMeasurementTool::Track

Definition at line 51 of file MuonTGMeasurementTool.h.

Constructor & Destructor Documentation

◆ MuonTGMeasurementTool()

Muon::MuonTGMeasurementTool::MuonTGMeasurementTool ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Constructor with AlgTool parameters.

Definition at line 34 of file MuonTGMeasurementTool.cxx.

34  :
35  AthAlgTool(type, name, parent), m_muonDetMgr(nullptr) {
36  declareInterface<Muon::IMuonTGMeasTool>(this);
37 }

◆ ~MuonTGMeasurementTool()

virtual Muon::MuonTGMeasurementTool::~MuonTGMeasurementTool ( )
virtualdefault

Member Function Documentation

◆ associatedLayer() [1/2]

const Trk::Layer * Muon::MuonTGMeasurementTool::associatedLayer ( Identifier  id,
Amg::Vector3D &  gp 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1112 of file MuonTGMeasurementTool.cxx.

1112  {
1113  // Get the messaging service, print where you are
1114  ATH_MSG_DEBUG("MuonTGMeasurementTool::associatedLayer");
1115  const Trk::Layer* lay = nullptr;
1116  // check input
1117  if (!id.get_identifier32().get_compact()) return lay;
1118 
1119  // rely on having misalignment uncertainty covered by span safety marge ( don't loose station from static volume
1120  // when misaligned
1121  const Trk::TrackingVolume* staticVol = getGeometry()->lowestStaticTrackingVolume(gp);
1122  const Trk::DetachedTrackingVolume* station = nullptr;
1123  if (staticVol && !staticVol->confinedDetachedVolumes().empty()) {
1124  Trk::ArraySpan<const Trk::DetachedTrackingVolume* const> detTV = staticVol->confinedDetachedVolumes();
1125  for (const auto *i : detTV) {
1126  if (i->layerRepresentation() && i->layerRepresentation()->layerType() > 0) {
1127  Identifier stId(i->layerRepresentation()->layerType());
1128  if (m_idHelperSvc->mdtIdHelper().stationName(stId) == m_idHelperSvc->mdtIdHelper().stationName(id) &&
1129  m_idHelperSvc->mdtIdHelper().stationEta(stId) == m_idHelperSvc->mdtIdHelper().stationEta(id) &&
1130  m_idHelperSvc->mdtIdHelper().stationPhi(stId) == m_idHelperSvc->mdtIdHelper().stationPhi(id)) {
1131  station = i;
1132  break;
1133  }
1134  }
1135  }
1136  }
1137 
1138  if (station) lay = associatedLayer(id, station->trackingVolume());
1139 
1140  return lay;
1141 }

◆ associatedLayer() [2/2]

const Trk::Layer * Muon::MuonTGMeasurementTool::associatedLayer ( Identifier  id,
const Trk::TrackingVolume *  vol 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1143 of file MuonTGMeasurementTool.cxx.

1143  {
1144  const Trk::Layer* lay = nullptr;
1145  // check input
1146  if (!vol) return lay;
1147 
1148  if (vol->confinedVolumes()) {
1149  Trk::BinnedArraySpan<Trk::TrackingVolume const * const > subVols = vol->confinedVolumes()->arrayObjects();
1150  Trk::BinnedArraySpan<Trk::TrackingVolume const * const >::iterator iter = subVols.begin();
1151  while (!lay && iter != subVols.end()) {
1152  lay = associatedLayer(id, *iter);
1153  if (lay) break;
1154  ++iter;
1155  }
1156  if (lay) return lay;
1157  }
1158 
1159  if (vol->confinedLayers()) {
1160  Trk::BinnedArraySpan<Trk::Layer const * const > ordLay = vol->confinedLayers()->arrayObjects();
1161  Trk::BinnedArraySpan<Trk::Layer const * const >::iterator iter = ordLay.begin();
1162  while (!lay && iter != ordLay.end()) {
1163  lay = match(id, *iter);
1164  if (lay) break;
1165  ++iter;
1166  }
1167  if (lay) return lay;
1168  }
1169 
1170  if (!vol->confinedArbitraryLayers().empty()) {
1171  Trk::ArraySpan<const Trk::Layer* const> unOrdLay = vol->confinedArbitraryLayers();
1172  Trk::ArraySpan<const Trk::Layer* const>::iterator iter = unOrdLay.begin();
1173  while (!lay && iter != unOrdLay.end()) {
1174  lay = match(id, *iter);
1175  if (lay) break;
1176  ++iter;
1177  }
1178  if (lay) return lay;
1179  }
1180 
1181  return lay;
1182 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType &   
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase &  hndl,
const std::string &  doc,
const SG::VarHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey &  hndl,
const std::string &  doc,
const SG::VarHandleKeyType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray &  hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType &   
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ detElToLayer()

const Trk::TrackParameters * Muon::MuonTGMeasurementTool::detElToLayer ( const Trk::Layer *  lay,
const Trk::TrackParameters *  parm,
Identifier  id 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 366 of file MuonTGMeasurementTool.cxx.

367  {
368  // Get the messaging service, print where you are
369  ATH_MSG_DEBUG("MuonTGMeasurementTool::detElToLayer");
370  const Trk::TrackParameters* projPar = nullptr;
371 
372  // check input
373  if (!lay || !parm || !(id.get_identifier32().get_compact() > 0)) return projPar;
374 
375  // check compatibility of layer info and required id ? this was already done when associating !
376  if (!lay->layerType()) return projPar;
377  Identifier layId(lay->layerType());
378 
379  ATH_MSG_DEBUG("MuonTGMeasurementTool::input ok");
380 
381  unsigned int hitType = 0;
382  if (m_idHelperSvc->isMdt(id)) hitType = 1;
383  if (m_idHelperSvc->isRpc(id)) hitType = 2;
384  if (m_idHelperSvc->isCsc(id)) hitType = 3;
385  if (m_idHelperSvc->isTgc(id)) hitType = 4;
386 
387  unsigned int layType = 0;
388  if (m_idHelperSvc->isMdt(layId)) layType = 1;
389  if (m_idHelperSvc->isRpc(layId)) layType = 2;
390  if (m_idHelperSvc->isCsc(layId)) layType = 3;
391  if (m_idHelperSvc->isTgc(layId)) layType = 4;
392 
393  if (layType != hitType) return projPar;
394 
395  if (hitType == 0) {
396  ATH_MSG_DEBUG("unknown hit technology");
397  return projPar;
398  }
399 
400  if (hitType == 1) {
401  // local position of the tube
402  int tube = m_idHelperSvc->mdtIdHelper().tube(id);
403  Amg::Vector2D locWire(0., lay->getRef() + (tube - 1) * 30.035);
404  if (fabs(lay->getRef()) > 10e6) {
405  double sign = (lay->getRef() > 0) ? 1. : -1.;
406  int dec = int(lay->getRef() / 1.e5);
407  double ref0 = dec / 1.e5;
408  double ref1 = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
409  locWire[0] = ref0;
410  locWire[1] = ref1;
411  int tube = m_idHelperSvc->mdtIdHelper().tube(id);
412  int tubeMax = m_idHelperSvc->mdtIdHelper().tubeMax(id);
413  if (tube > 6 && tubeMax - tube > 5) locWire[0] = 0.;
414  }
415  if (sqrt(locWire[0] * locWire[0] + locWire[1] * locWire[1]) > 2000.) {
416  ATH_MSG_WARNING(name() << " wire shift out bounds for MDT tube :" << m_idHelperSvc->mdtIdHelper().stationName(id) << ","
417  << m_idHelperSvc->mdtIdHelper().stationEta(id) << "," << m_idHelperSvc->mdtIdHelper().stationPhi(id)
418  << ": abandon projection");
419  return projPar;
420  }
421  // local position (tube)
422  // const Amg::Vector2D locPos = parm->localPosition();
423  // direction at the layer
424  Amg::Vector3D dir = parm->momentum().unit();
425  Amg::Vector3D locDir = (lay->surfaceRepresentation().transform().inverse()) * dir;
426  Amg::Vector3D normal = lay->surfaceRepresentation().normal();
427  Amg::Vector3D wireDir(1., 0., 0.);
428  const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
429  if (!laySurf) return projPar;
430  //
431  double ND = dir.dot(normal);
432  double DL = locDir.dot(wireDir);
433  double A = sqrt(1. - DL * DL);
434  AmgMatrix(5, 5) mdtProj;
435  mdtProj.setIdentity();
436  //
437  double ori = (laySurf->transform().inverse() * parm->associatedSurface().transform()).rotation()(1, 1) > 0. ? -1. : 1.;
438  if (A > 10e-6) {
439  mdtProj(0, 1) = ND / A;
440  mdtProj(1, 0) = ori;
441  mdtProj(1, 1) = DL / A * sqrt(A * A - ND * ND) / A;
442  mdtProj(0, 0) = 0.;
443  } else {
444  mdtProj(0, 1) = 1.;
445  mdtProj(1, 0) = ori;
446  mdtProj(1, 1) = 0.;
447  mdtProj(0, 0) = 0.;
448  }
449  Amg::VectorX locPar = parm->parameters();
450  AmgMatrix(5, 5) mdtProjInv = mdtProj.inverse();
451  Amg::VectorX pPar = mdtProjInv * locPar;
452  pPar[0] += locWire[0];
453  pPar[1] += locWire[1];
454  ATH_MSG_DEBUG("back projected parameters(MDT->layer):" << pPar);
455 
456  std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(mdtProjInv);
457 
458  ATH_MSG_DEBUG("back projected covariance(MDT->layer):" << (*projEM));
459  projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *laySurf, projEM);
460  // delete projEM;
461  }
462 
463  if (hitType == 2) {
464  // local position at detEl
465  // const Amg::Vector2D locPos = parm->localPosition();
466  const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
467  if (!laySurf) return projPar;
468  //
469  Amg::VectorX locPar = parm->parameters();
470  AmgMatrix(5, 5)* pMx = nullptr;
471  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
472  pMx = m_rpcProjPhi.get();
473  else
474  pMx = m_rpcProjEta.get();
475 
476  double eta = 1.;
477  double sign = (m_idHelperSvc->rpcIdHelper().measuresPhi(id) && m_idHelperSvc->rpcIdHelper().doubletPhi(id) == 2) ? -1. : 1.;
478 
479  double ref = (lay->getRef() > 10000.) ? lay->getRef() - 20000. : lay->getRef();
480  double zswap = (lay->getRef() > 10000.) ? -1. : 1.;
481 
482  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
483  locPar[0] *= eta;
484  else
485  locPar[1] *= -eta;
486  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
487  locPar[1] *= zswap;
488  else
489  locPar[0] *= zswap;
490 
491  Amg::VectorX pPar = (*pMx) * locPar;
492  pPar[0] += sign * ref;
493 
494  //
495  ATH_MSG_DEBUG("back projected parameters(RPC->layer):" << m_idHelperSvc->rpcIdHelper().measuresPhi(id) << "," << pPar);
496 
497  std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarityT(*pMx);
498 
499  ATH_MSG_DEBUG("back projected covariance(RPC->layer):" << (*projEM));
500  projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *laySurf, projEM);
501  // delete projEM;
502 
503  if ((parm->position() - projPar->position()).mag() > 0.001) {
504  ATH_MSG_DEBUG("geometrical RPC projection (detElToLayer) for hit : "
505  << m_idHelperSvc->rpcIdHelper().measuresPhi(id) << "," << m_idHelperSvc->rpcIdHelper().stationName(id) << ","
506  << m_idHelperSvc->rpcIdHelper().stationEta(id) << "," << m_idHelperSvc->rpcIdHelper().stationPhi(id) << ","
507  << m_idHelperSvc->rpcIdHelper().doubletPhi(id) << "," << m_idHelperSvc->rpcIdHelper().doubletR(id) << ","
508  << m_idHelperSvc->rpcIdHelper().doubletZ(id));
509  Amg::Vector2D locPos;
510  const Amg::Vector3D& globPos = parm->position();
511  bool onSurface = laySurf->globalToLocal(globPos, globPos, locPos);
512  if (onSurface) {
513  pPar[0] = locPos[0];
514  pPar[1] = locPos[1];
515  delete projPar;
516  projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *laySurf, projEM);
517  // delete projEM;
518  } else {
519  delete projPar;
520  return nullptr;
521  }
522  }
523  }
524 
525  if (hitType == 3) {
526  // local position of detEl
527  // const Amg::Vector2D locPos = parm->localPosition();
528  const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
529  if (!laySurf) return projPar;
530  // dealing with parallel displaced planes
531  double diff = m_idHelperSvc->cscIdHelper().measuresPhi(id) ? 1.55 : -1.55;
532  //
533  Amg::Vector3D layNormal = lay->surfaceRepresentation().normal();
534  double DN = parm->momentum().dot(layNormal);
535  double t = diff / DN;
536  // displacement of planes
537  const Amg::Vector2D csc_shift(0., lay->getRef());
538  // projection : take into account possible misalignment ;
539  AmgMatrix(5, 5)* pMx = nullptr;
540  if (m_idHelperSvc->cscIdHelper().measuresPhi(id))
541  pMx = m_tgcProjPhi.get();
542  else
543  pMx = m_tgcProjEta.get();
544  AmgMatrix(5, 5) pMxInv = pMx->inverse();
545  Amg::VectorX parProj = pMxInv * parm->parameters();
546  Amg::Vector3D corrLocPos = lay->surfaceRepresentation().center() - t * parm->momentum() + t * DN * layNormal;
547  Amg::Vector2D locCorrLay;
548  bool onSurface = lay->surfaceRepresentation().globalToLocal(corrLocPos, corrLocPos, locCorrLay);
549  if (onSurface && locCorrLay.size() > 0) {
550  parProj[0] += locCorrLay[Trk::locX] + csc_shift[Trk::locX];
551  parProj[1] += locCorrLay[Trk::locY] + csc_shift[Trk::locY];
552  ATH_MSG_DEBUG("back projected parameters(CSC->layer):" << m_idHelperSvc->cscIdHelper().measuresPhi(id) << "," << parProj);
553  }
554  std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(pMxInv);
555 
556  ATH_MSG_DEBUG("back projected covariance(CSC->layer):" << (*projEM));
557  projPar = new Trk::AtaPlane(parProj[0], parProj[1], parProj[2], parProj[3], parProj[4], *laySurf, projEM);
558 
559  ATH_MSG_DEBUG("test CSC projection:detElToLayer:" << parm->position() << "," << projPar->position());
560  }
561 
562  if (hitType == 4) {
563  const Trk::PlaneSurface* laySurf = dynamic_cast<const Trk::PlaneSurface*>(&(lay->surfaceRepresentation()));
564  if (!laySurf) { return projPar; }
565 
566  AmgMatrix(5, 5)* pMx = nullptr;
567  if (m_idHelperSvc->tgcIdHelper().isStrip(id))
568  pMx = m_tgcProjPhi.get();
569  else
570  pMx = m_tgcProjEta.get();
571  AmgMatrix(5, 5) pMxInv = pMx->inverse();
572  Amg::VectorX locPar = pMxInv * parm->parameters();
573  ATH_MSG_DEBUG("back projected parameters(TGC->layer):" << m_idHelperSvc->tgcIdHelper().isStrip(id) << "," << locPar);
574 
575  std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(pMxInv);
576 
577  ATH_MSG_DEBUG("back projected covariance(TGC->layer):" << (*projEM));
578  projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *laySurf, projEM);
579 
580  // verify
581  if (m_alignedMode && (parm->position() - projPar->position()).mag() > 0.001) {
582  ATH_MSG_DEBUG("geometrical TGC projection ( detEl2Layer ):"
583  << m_idHelperSvc->tgcIdHelper().stationName(id) << "," << m_idHelperSvc->tgcIdHelper().stationEta(id) << ","
584  << m_idHelperSvc->tgcIdHelper().stationPhi(id) << "," << m_idHelperSvc->tgcIdHelper().isStrip(id));
585  Amg::Vector2D locPos;
586  const Amg::Vector3D& globPos = parm->position();
587  bool onSurface = laySurf->globalToLocal(globPos, globPos, locPos);
588  if (onSurface) {
589  locPar[0] = locPos[0];
590  locPar[1] = locPos[1];
591  delete projPar;
592  projPar = nullptr;
593  projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *laySurf, projEM);
594  } else {
595  delete projPar;
596  return nullptr;
597  }
598  }
599  }
600  return projPar;
601 }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ getGeometry()

const Trk::TrackingGeometry* Muon::MuonTGMeasurementTool::getGeometry ( ) const
inlineprivate

Good old way of retrieving the volume via the geometry service

Definition at line 98 of file MuonTGMeasurementTool.h.

98  {
100  if (m_trackingGeometryReadKey.empty()) { return m_trackingGeometrySvc->trackingGeometry(); }
101  SG::ReadCondHandle<Trk::TrackingGeometry> handle(m_trackingGeometryReadKey, Gaudi::Hive::currentContext());
102  if (!handle.isValid()) {
103  ATH_MSG_WARNING("Could not retrieve a valid tracking geometry");
104  return nullptr;
105  }
106  return handle.cptr();
107  }

◆ initialize()

StatusCode Muon::MuonTGMeasurementTool::initialize ( )
overridevirtual

Definition at line 40 of file MuonTGMeasurementTool.cxx.

40  {
41  // Get the messaging service, print where you are
42  ATH_MSG_INFO("MuonTGMeasurementTool::initialize()");
43 
44  ATH_CHECK(m_idHelperSvc.retrieve());
45  ATH_CHECK(m_DetectorManagerKey.initialize(!m_useDSManager));
46  if (m_useDSManager) { ATH_CHECK(detStore()->retrieve(m_muonDetMgr)); }
47 
48  // define projection matrices
49  m_tgcProjEta = std::make_unique<AmgMatrix(5, 5)>();
50  m_tgcProjEta->setIdentity();
51  (*m_tgcProjEta)(0, 0) = 0.;
52  (*m_tgcProjEta)(1, 1) = 0.;
53  (*m_tgcProjEta)(0, 1) = 1.;
54  (*m_tgcProjEta)(1, 0) = -1.;
55  m_tgcProjPhi = std::make_unique<AmgMatrix(5, 5)>();
56  m_tgcProjPhi->setIdentity();
57 
58  m_rpcProjEta = std::make_unique<AmgMatrix(5, 5)>();
59  m_rpcProjEta->setIdentity();
60  (*m_rpcProjEta)(0, 0) = 0.;
61  (*m_rpcProjEta)(1, 1) = 0.;
62  (*m_rpcProjEta)(0, 1) = 1.;
63  (*m_rpcProjEta)(1, 0) = 1.;
64  m_rpcProjPhi = std::make_unique<AmgMatrix(5, 5)>();
65  m_rpcProjPhi->setIdentity();
66 
67  if (!m_trackingGeometryReadKey.empty()) {
68  ATH_CHECK(m_trackingGeometryReadKey.initialize());
69  } else {
70  ATH_CHECK(m_trackingGeometrySvc.retrieve());
71  }
72  return StatusCode::SUCCESS;
73 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& Muon::IMuonTGMeasTool::interfaceID ( )
inlinestaticinherited

AlgTool interface method.

Definition at line 29 of file IMuonTGMeasTool.h.

29  {
30  static const InterfaceID IID_IMuonTGMeasTool("IMuonTGMeasTool", 1, 0);
31  return IID_IMuonTGMeasTool;
32  };

◆ layerToDetEl()

const Trk::TrackParameters * Muon::MuonTGMeasurementTool::layerToDetEl ( const Trk::Layer *  lay,
const Trk::TrackParameters *  parm,
Identifier  id 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 74 of file MuonTGMeasurementTool.cxx.

75  {
76  const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
77  if (!m_useDSManager) {
78  SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
79  MuonDetMgr = DetectorManagerHandle.cptr();
80  if (MuonDetMgr == nullptr) {
81  ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object");
82  // return StatusCode::FAILURE;
83  }
84  }
85 
86  // Get the messaging service, print where you are
87  ATH_MSG_DEBUG("MuonTGMeasurementTool::layerToDetEl");
88  const Trk::TrackParameters* projPar = nullptr;
89  // const Amg::Vector2D* locPos = 0;
90  // check input
91  if (!lay || !parm || !id.get_identifier32().get_compact()) return projPar;
92 
93  // get tracking geometry
94 
95  // check compatibility of layer info and required id ? this was already done when associating !
96  if (!lay->layerType()) return projPar;
97  Identifier layId(lay->layerType());
98 
99  unsigned int hitType = 0;
100  if (m_idHelperSvc->isMdt(id)) hitType = 1;
101  if (m_idHelperSvc->isRpc(id)) hitType = 2;
102  if (m_idHelperSvc->isCsc(id)) hitType = 3;
103  if (m_idHelperSvc->isTgc(id)) hitType = 4;
104 
105  unsigned int layType = 0;
106  if (m_idHelperSvc->isMdt(layId)) layType = 1;
107  if (m_idHelperSvc->isRpc(layId)) layType = 2;
108  if (m_idHelperSvc->isCsc(layId)) layType = 3;
109  if (m_idHelperSvc->isTgc(layId)) layType = 4;
110 
111  if (layType != hitType) return projPar;
112 
113  if (hitType == 0) {
114  ATH_MSG_ERROR("unknown hit technology");
115  return projPar;
116  } else {
117  ATH_MSG_DEBUG("hit technology:" << hitType);
118  }
119 
120  ATH_MSG_DEBUG("extrapolated covariance:" << parm->covariance());
121 
122  if (hitType == 1) {
123  const MuonGM::MdtReadoutElement* mdtROE = MuonDetMgr->getMdtReadoutElement(id);
124  if (!mdtROE) {
125  ATH_MSG_WARNING(name() << "MDT readout element not found");
126  return projPar;
127  }
128  // local position of tube
129  int tube = m_idHelperSvc->mdtIdHelper().tube(id);
130  Amg::Vector2D locWire(0., lay->getRef() + (tube - 1) * 30.035);
131  if (fabs(lay->getRef()) > 10e6) {
132  double sign = (lay->getRef() > 0) ? 1. : -1.;
133  int dec = int(lay->getRef() / 1.e5);
134  double ref0 = dec / 1.e5;
135  double ref1 = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
136  locWire[0] = ref0;
137  locWire[1] = ref1;
138  int tube = m_idHelperSvc->mdtIdHelper().tube(id);
139  int tubeMax = m_idHelperSvc->mdtIdHelper().tubeMax(id);
140  if (tube > 6 && tubeMax - tube > 5) locWire[0] = 0.;
141  }
142  if (sqrt(locWire[0] * locWire[0] + locWire[1] * locWire[1]) > 2000.) {
143  ATH_MSG_WARNING(name() << " wire shift out bounds for MDT tube :" << m_idHelperSvc->mdtIdHelper().stationName(id) << ","
144  << m_idHelperSvc->mdtIdHelper().stationEta(id) << "," << m_idHelperSvc->mdtIdHelper().stationPhi(id)
145  << ": abandon projection");
146  return projPar;
147  }
148  // direction at the layer
149  Amg::Vector3D dir = parm->momentum().unit();
150  Amg::Vector3D locDir = (lay->surfaceRepresentation().transform().inverse()) * dir;
151  Amg::Vector3D wireDir(1., 0., 0.);
152  //
153  double ND = dir.dot(lay->surfaceRepresentation().normal());
154  double DL = locDir.dot(wireDir);
155  double A = sqrt(1. - DL * DL);
156  AmgMatrix(5, 5) mdtProj;
157  mdtProj.setIdentity();
158  //
159  const Trk::StraightLineSurface* tubeSurf = dynamic_cast<const Trk::StraightLineSurface*>(&(mdtROE->surface(id)));
160  double ori = (lay->surfaceRepresentation().transform().inverse() * tubeSurf->transform()).rotation()(1, 1) > 0. ? -1. : 1.;
161  if (A > 10e-6) {
162  mdtProj(0, 1) = ND / A;
163  mdtProj(1, 0) = ori;
164  mdtProj(1, 1) = DL / A * sqrt(A * A - ND * ND) / A;
165  mdtProj(0, 0) = 0.;
166  } else {
167  mdtProj(0, 1) = 1.;
168  mdtProj(1, 0) = ori;
169  mdtProj(1, 1) = 0.;
170  mdtProj(0, 0) = 0.;
171  }
172  //
173  Amg::VectorX locPar = parm->parameters();
174  locPar[0] -= locWire[0];
175  locPar[1] -= locWire[1];
176  Amg::VectorX pPar = mdtProj * locPar;
177  ATH_MSG_DEBUG("projected parameters(layer->MDT):" << pPar);
178 
179  if (parm->covariance()) {
180  std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(mdtProj);
181  ATH_MSG_DEBUG("projected covariance(layer->MDT):" << (*projEM));
182  projPar = new Trk::AtaStraightLine(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *tubeSurf, projEM);
183 
184  } else {
185  projPar = new Trk::AtaStraightLine(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *tubeSurf);
186  }
187  }
188 
189  if (hitType == 2) {
190  //
191  const MuonGM::RpcReadoutElement* rpcROE = MuonDetMgr->getRpcReadoutElement(id);
192  if (!rpcROE) return projPar;
193  const Trk::PlaneSurface* stripSurf = dynamic_cast<const Trk::PlaneSurface*>(&(rpcROE->surface(id)));
194  if (!stripSurf) return projPar;
195  // decode ref position
196  Amg::VectorX locPar = parm->parameters();
197  // projection matrix
198  AmgMatrix(5, 5)* pMx = nullptr;
199  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
200  pMx = m_rpcProjPhi.get();
201  else
202  pMx = m_rpcProjEta.get();
203  // projected parameters
204  double eta = 1.;
205  double sign = (m_idHelperSvc->rpcIdHelper().measuresPhi(id) && m_idHelperSvc->rpcIdHelper().doubletPhi(id) == 2) ? -1. : 1.;
206  double zswap = (lay->getRef() > 10000.) ? -1. : 1.;
207  double ref = (zswap < 0.) ? lay->getRef() - 20000. : lay->getRef();
208  locPar[0] -= sign * ref;
209  Amg::VectorX pPar = (*pMx) * locPar;
210  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
211  pPar[0] *= eta;
212  else
213  pPar[1] *= -eta;
214  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
215  pPar[1] *= zswap;
216  else
217  pPar[0] *= zswap;
218 
219  std::optional<AmgMatrix(5, 5)> projEM;
220  projEM.emplace();
221  if (parm->covariance()) {
222  projEM = parm->covariance()->similarity(*pMx);
223  projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf, projEM);
224 
225  } else {
226  projPar = new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf);
227  projEM = std::nullopt;
228  }
229 
230  if (m_alignedMode && (parm->position() - projPar->position()).mag() > 0.001) {
231  ATH_MSG_DEBUG("geometrical RPC projection (layerToDetEl) for hit : "
232  << m_idHelperSvc->rpcIdHelper().measuresPhi(id) << "," << m_idHelperSvc->rpcIdHelper().stationName(id) << ","
233  << m_idHelperSvc->rpcIdHelper().stationEta(id) << "," << m_idHelperSvc->rpcIdHelper().stationPhi(id) << ","
234  << m_idHelperSvc->rpcIdHelper().doubletPhi(id) << "," << m_idHelperSvc->rpcIdHelper().doubletR(id) << ","
235  << m_idHelperSvc->rpcIdHelper().doubletZ(id));
236  Amg::Vector2D locPos;
237  const Amg::Vector3D& globPos = parm->position();
238  bool onSurface = stripSurf->globalToLocal(globPos, globPos, locPos);
239  if (onSurface) {
240  pPar[0] = locPos[0];
241  pPar[1] = locPos[1];
242  delete projPar;
243  projPar = nullptr;
244  projPar = projEM ? new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf, projEM)
245  : new Trk::AtaPlane(pPar[0], pPar[1], pPar[2], pPar[3], pPar[4], *stripSurf);
246  } else {
247  delete projPar;
248  return nullptr;
249  }
250  }
251  }
252 
253  if (hitType == 3) {
254  // local position of detEl
255  const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(id);
256  if (!cscROE) {
257  ATH_MSG_WARNING(name() << "CSC readout element not found");
258  return projPar;
259  }
260 
261  const Trk::PlaneSurface* stripSurf = dynamic_cast<const Trk::PlaneSurface*>(&(cscROE->surface(id)));
262  if (!stripSurf) return projPar;
263  // dealing with displaced planes, possibly sligthly displaced
264  double diff = m_idHelperSvc->cscIdHelper().measuresPhi(id) ? 1.55 : -1.55;
265  // distance between planes (assuming parallel planes)
266  Amg::Vector3D layNormal = lay->surfaceRepresentation().normal();
267  double DN = parm->momentum().dot(layNormal);
268  double t = diff / DN;
269  // displacement ( as projection of strip surface center on the layer )
270  const Amg::Vector2D csc_shift(0., lay->getRef());
271  // correct local position due to plane distance : use method independent on misalignment
272  Amg::Vector3D corrLocPos = parm->position() - t * parm->momentum() + t * DN * layNormal;
273  Amg::Vector2D locCorrLay;
274  bool onSurface = lay->surfaceRepresentation().globalToLocal(corrLocPos, corrLocPos, locCorrLay);
275 
276  if (!onSurface) {
277  ATH_MSG_WARNING(name() << ": misplaced CSC " << id);
278  return projPar;
279  }
280 
281  Amg::VectorX parProj(5);
282  parProj[0] = locCorrLay[Trk::locX] - csc_shift[Trk::locX];
283  parProj[1] = locCorrLay[Trk::locY] - csc_shift[Trk::locY];
284  parProj[2] = parm->parameters()[Trk::phi];
285  parProj[3] = parm->parameters()[Trk::theta];
286  parProj[4] = parm->parameters()[Trk::qOverP];
287  //
288  AmgMatrix(5, 5)* pMx = nullptr;
289  if (m_idHelperSvc->cscIdHelper().measuresPhi(id))
290  pMx = m_tgcProjPhi.get();
291  else
292  pMx = m_tgcProjEta.get();
293  Amg::VectorX locPar = (*pMx) * parProj;
294  ATH_MSG_DEBUG("projected parameters (layer->CSC):" << m_idHelperSvc->cscIdHelper().measuresPhi(id) << "," << locPar);
295 
296  if (parm->covariance()) {
297  std::optional<AmgMatrix(5, 5)> projEM = parm->covariance()->similarity(*pMx);
298 
299  ATH_MSG_DEBUG("projected covariance (layer->CSC):" << *projEM);
300 
301  projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf, projEM);
302 
303  } else
304  projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf);
305 
306  ATH_MSG_DEBUG("test CSC projection:layerToDetEl:" << parm->position() << "," << projPar->position());
307  }
308 
309  if (hitType == 4) {
310  // local position at layer
311  const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(id);
312  if (!tgcROE) {
313  ATH_MSG_WARNING(name() << "TGC readout element not found");
314  return projPar;
315  }
316  const Trk::PlaneSurface* stripSurf = dynamic_cast<const Trk::PlaneSurface*>(&(tgcROE->surface(id)));
317  if (!stripSurf) return projPar;
318  //
319  AmgMatrix(5, 5)* pMx = nullptr;
320  if (m_idHelperSvc->tgcIdHelper().isStrip(id))
321  pMx = m_tgcProjPhi.get();
322  else
323  pMx = m_tgcProjEta.get();
324  Amg::VectorX locPar = (*pMx) * parm->parameters();
325  ATH_MSG_DEBUG("projected parameters (layer->TGC):" << m_idHelperSvc->tgcIdHelper().isStrip(id) << "," << locPar << ","
326  << stripSurf);
327 
328  std::optional<AmgMatrix(5, 5)> projEM = std::nullopt;
329  bool bcov = false;
330 
331  if (parm->covariance()) {
332  projEM = parm->covariance()->similarity(*pMx);
333  bcov = true;
334  ATH_MSG_DEBUG("projected covariance (layer->TGC):" << (*projEM));
335  projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf, projEM);
336 
337  } else {
338  projPar = new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf);
339  }
340 
341  // verify
342  if (m_alignedMode && (parm->position() - projPar->position()).mag() > 0.001) {
343  ATH_MSG_DEBUG("geometrical TGC projection ( layer2detEl ):"
344  << m_idHelperSvc->tgcIdHelper().stationName(id) << "," << m_idHelperSvc->tgcIdHelper().stationEta(id) << ","
345  << m_idHelperSvc->tgcIdHelper().stationPhi(id) << "," << m_idHelperSvc->tgcIdHelper().isStrip(id));
346  Amg::Vector2D locPos;
347  const Amg::Vector3D& globPos = parm->position();
348  bool onSurface = stripSurf->globalToLocal(globPos, globPos, locPos);
349  if (onSurface) {
350  locPar[0] = locPos[0];
351  locPar[1] = locPos[1];
352  delete projPar;
353  projPar = nullptr;
354  projPar = bcov ? new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf, projEM)
355  : new Trk::AtaPlane(locPar[0], locPar[1], locPar[2], locPar[3], locPar[4], *stripSurf);
356  } else {
357  delete projPar;
358  return nullptr;
359  }
360  }
361  }
362 
363  return projPar;
364 }

◆ match()

const Trk::Layer * Muon::MuonTGMeasurementTool::match ( Identifier  id,
const Trk::Layer *  lay 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1184 of file MuonTGMeasurementTool.cxx.

1184  {
1185  const Trk::Layer* mLay = nullptr;
1186  if (!id.get_identifier32().get_compact() || !lay || !lay->layerType()) return mLay;
1187 
1188  Identifier layId(lay->layerType());
1189 
1190  if (m_idHelperSvc->isMdt(id) && m_idHelperSvc->isMdt(layId)) {
1191  if (m_idHelperSvc->mdtIdHelper().multilayer(layId) == m_idHelperSvc->mdtIdHelper().multilayer(id) &&
1192  m_idHelperSvc->mdtIdHelper().tubeLayer(layId) == m_idHelperSvc->mdtIdHelper().tubeLayer(id))
1193  return lay;
1194  }
1195 
1196  if (m_idHelperSvc->isRpc(id) && m_idHelperSvc->isRpc(layId)) {
1197  if (m_idHelperSvc->rpcIdHelper().doubletR(layId) == m_idHelperSvc->rpcIdHelper().doubletR(id) &&
1198  m_idHelperSvc->rpcIdHelper().doubletZ(layId) == m_idHelperSvc->rpcIdHelper().doubletZ(id) &&
1199  m_idHelperSvc->rpcIdHelper().gasGap(layId) == m_idHelperSvc->rpcIdHelper().gasGap(id))
1200  return lay;
1201  }
1202 
1203  if (m_idHelperSvc->isTgc(id) && m_idHelperSvc->isTgc(layId)) {
1204  if (m_idHelperSvc->tgcIdHelper().gasGap(layId) == m_idHelperSvc->tgcIdHelper().gasGap(id)) return lay;
1205  }
1206 
1207  if (m_idHelperSvc->isCsc(id) && m_idHelperSvc->isCsc(layId)) {
1208  if (m_idHelperSvc->cscIdHelper().chamberLayer(layId) == m_idHelperSvc->cscIdHelper().chamberLayer(id) &&
1209  m_idHelperSvc->cscIdHelper().wireLayer(layId) == m_idHelperSvc->cscIdHelper().wireLayer(id))
1210  return lay;
1211  }
1212 
1213  return mLay;
1214 }

◆ measToLayer()

const Trk::RIO_OnTrack * Muon::MuonTGMeasurementTool::measToLayer ( const Trk::Layer *  lay,
const Trk::TrackParameters *  parm,
const Trk::RIO_OnTrack *  rio 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 603 of file MuonTGMeasurementTool.cxx.

604  {
605  const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
606  if (!m_useDSManager) {
607  SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
608  MuonDetMgr = DetectorManagerHandle.cptr();
609  if (MuonDetMgr == nullptr) {
610  ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object");
611  // return StatusCode::FAILURE;
612  }
613  }
614  // Get the messaging service, print where you are
615  ATH_MSG_DEBUG("MuonTGMeasurementTool::measToLayer");
616  const Trk::RIO_OnTrack* projRIO = nullptr;
617 
618  // check input
619  if (!lay || !parm || !rio) return projRIO;
620 
621  // check compatibility of layer info and required id ? this was already done when associating !
622  Identifier id = rio->identify();
623 
624  ATH_MSG_DEBUG("MuonTGMeasurementTool::input ok");
625 
626  unsigned int hitType = 0;
627  if (m_idHelperSvc->isMdt(id)) hitType = 1;
628  if (m_idHelperSvc->isRpc(id)) hitType = 2;
629  if (m_idHelperSvc->isCsc(id)) hitType = 3;
630  if (m_idHelperSvc->isTgc(id)) hitType = 4;
631 
632  if (hitType == 1) {
633  // local position of the tube
634  int tube = m_idHelperSvc->mdtIdHelper().tube(id);
635  Amg::Vector2D locWire(0., lay->getRef() + (tube - 1) * 30.035);
636  if (fabs(lay->getRef()) > 10e6) {
637  double sign = (lay->getRef() > 0) ? 1. : -1.;
638  int dec = int(lay->getRef() / 1.e5);
639  double ref0 = dec / 1.e5;
640  double ref1 = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
641  locWire[0] = ref0;
642  locWire[1] = ref1;
643  int tube = m_idHelperSvc->mdtIdHelper().tube(id);
644  int tubeMax = m_idHelperSvc->mdtIdHelper().tubeMax(id);
645  if (tube > 6 && tubeMax - tube > 5) locWire[0] = 0.;
646  }
647  // direction at the layer
648  Amg::Vector3D dir = parm->momentum().unit();
649  Amg::Vector3D locDir = (lay->surfaceRepresentation().transform().inverse()) * dir;
650  Amg::Vector3D normal = lay->surfaceRepresentation().normal();
651  Amg::Vector3D wireDir(1., 0., 0.);
652  //
653  double ND = dir.dot(normal);
654  double DL = locDir.dot(wireDir);
655  double A = sqrt(1. - DL * DL);
656  double A_ND = A / ND;
657  //
658  double locLay = A_ND * rio->localParameters()[Trk::locR] + locWire[1];
659  // create (fake!) rio ( rio image on TG layer )
660  IdentifierHash idHash(0);
661  const MuonGM::MdtReadoutElement* mdtROE = MuonDetMgr->getMdtReadoutElement(id);
662  auto cov = Amg::MatrixX();
663  cov = A_ND * A_ND * rio->localCovariance();
664  Muon::MdtDriftCircleStatus status = Muon::MdtStatusDriftTime;
665  const Muon::MdtPrepData* mdtPrd =
666  new Muon::MdtPrepData(id,
667  Amg::Vector2D(locLay, 0.),
668  cov,
669  mdtROE,
670  0,
671  0,
672  status);
673  Muon::MuonDriftCircleErrorStrategyInput stratbits(
674  MuonDriftCircleErrorStrategy::UnknownStrategy);
675  Muon::MuonDriftCircleErrorStrategy strat(stratbits);
676  const Muon::MdtDriftCircleOnTrack* mdtRio =
677  new Muon::MdtDriftCircleOnTrack(mdtPrd, Trk::LocalParameters(Trk::DefinedParameter(locLay, Trk::locY)),
678  Amg::MatrixX(mdtPrd->localCovariance()), 0., Trk::DECIDED, 0.0, strat);
679  return mdtRio;
680  }
681 
682  else if (hitType == 2) {
683  //
684  double eta = (m_idHelperSvc->rpcIdHelper().stationEta(id) < 0) ? -1. : 1.;
685  double sign = (m_idHelperSvc->rpcIdHelper().measuresPhi(id) && m_idHelperSvc->rpcIdHelper().doubletPhi(id) == 2) ? -1. : 1.;
686 
687  double locPos = rio->localParameters()[Trk::locX];
688 
689  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id)) locPos *= eta;
690 
691  double ref = (lay->getRef() > 10000.) ? lay->getRef() - 20000. : lay->getRef();
692  double zswap = (lay->getRef() > 10000.) ? -1. : 1.;
693 
694  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
695  locPos += sign * ref;
696  else
697  locPos *= zswap;
698  //
699  const MuonGM::RpcReadoutElement* rpcROE = MuonDetMgr->getRpcReadoutElement(id);
700  const IdentifierHash idHash(0);
701  std::vector<Identifier> rdoList;
702  rdoList.push_back(id);
703  const Muon::RpcPrepData* rpcPrd = new Muon::RpcPrepData(id, idHash, Amg::Vector2D(locPos, 0.), rdoList,
704  Amg::MatrixX(rio->localCovariance()), rpcROE, float(0.), 0, 0);
705  const Muon::RpcClusterOnTrack* rpcRio = nullptr;
706  if (m_idHelperSvc->rpcIdHelper().measuresPhi(id))
707  rpcRio = new Muon::RpcClusterOnTrack(rpcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locX)),
708  Amg::MatrixX(rpcPrd->localCovariance()), parm->localPosition()[Trk::locY]);
709  else
710  rpcRio = new Muon::RpcClusterOnTrack(rpcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locY)),
711  Amg::MatrixX(rpcPrd->localCovariance()), parm->localPosition()[Trk::locX]);
712  return rpcRio;
713  }
714 
715  else if (hitType == 3) {
716  // dealing with parallel displaced planes
717  double diff = m_idHelperSvc->cscIdHelper().measuresPhi(id) ? 1.55 : -1.55;
718  //
719  Amg::Vector3D layNormal = lay->surfaceRepresentation().normal();
720  double DN = parm->momentum().dot(layNormal);
721  double t = diff / DN;
722  // displacement of planes
723  const Amg::Vector2D csc_shift(0., lay->getRef());
724  Amg::Vector3D corrLocPos = lay->surfaceRepresentation().center() - t * parm->momentum() + t * DN * layNormal;
725  Amg::Vector2D locCorrLay;
726  bool onSurface = lay->surfaceRepresentation().globalToLocal(corrLocPos, corrLocPos, locCorrLay);
727  //
728  // if( !locCorrLay.size()==2) { ??
729  if (!onSurface) { return projRIO; }
730  double locPos;
731  if (m_idHelperSvc->cscIdHelper().measuresPhi(id)) {
732  locPos = rio->localParameters()[Trk::locX] + locCorrLay[Trk::locX] + csc_shift[Trk::locX];
733  } else {
734  locPos = rio->localParameters()[Trk::locX] + locCorrLay[Trk::locY] + csc_shift[Trk::locY];
735  }
736  //
737  const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(id);
738  IdentifierHash idHash(0);
739  std::vector<Identifier> rdoList;
740  rdoList.push_back(id);
741  Muon::CscClusterStatus status = Muon::CscStatusSimple;
742  const Muon::CscPrepData* cscPrd =
743  new Muon::CscPrepData(id,
744  idHash,
745  Amg::Vector2D(locPos, 0.),
746  rdoList,
747  Amg::MatrixX(rio->localCovariance()),
748  cscROE,
749  0,
750  0.,
751  status);
752  const Muon::CscClusterOnTrack* cscRio = nullptr;
753  if (m_idHelperSvc->cscIdHelper().measuresPhi(id))
754  cscRio = new Muon::CscClusterOnTrack(cscPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locX)),
755  Amg::MatrixX(cscPrd->localCovariance()), parm->localPosition()[Trk::locY], cscPrd->status());
756  else
757  cscRio = new Muon::CscClusterOnTrack(cscPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locY)),
758  Amg::MatrixX(cscPrd->localCovariance()), parm->localPosition()[Trk::locX], cscPrd->status());
759  return cscRio;
760  }
761 
762  else if (hitType == 4) {
763  //
764  double locPos = rio->localParameters()[Trk::locX];
765  //
766  const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(id);
767  IdentifierHash idHash(0);
768  std::vector<Identifier> rdoList;
769  rdoList.push_back(id);
770  const Muon::TgcPrepData* tgcPrd =
771  new Muon::TgcPrepData(id,
772  idHash,
773  Amg::Vector2D(locPos, 0.),
774  rdoList,
775  Amg::MatrixX(rio->localCovariance()),
776  tgcROE);
777  const Muon::TgcClusterOnTrack* tgcRio = nullptr;
778  if (m_idHelperSvc->tgcIdHelper().isStrip(id)) {
779  Amg::Vector2D loc(locPos, parm->localPosition()[Trk::locY]);
780  tgcRio = new Muon::TgcClusterOnTrack(tgcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locX)),
781  Amg::MatrixX(tgcPrd->localCovariance()), parm->localPosition()[Trk::locY]);
782  } else {
783  tgcRio = new Muon::TgcClusterOnTrack(tgcPrd, Trk::LocalParameters(Trk::DefinedParameter(locPos, Trk::locY)),
784  Amg::MatrixX(tgcPrd->localCovariance()), parm->localPosition()[Trk::locX]);
785  }
786  return tgcRio;
787  }
788 
789  ATH_MSG_DEBUG("unknown hit technology");
790  return projRIO;
791 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ nearestDetEl()

const Identifier Muon::MuonTGMeasurementTool::nearestDetEl ( const Trk::Layer *  lay,
const Trk::TrackParameters *  parm,
bool  measPhi,
double &  pitch 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 793 of file MuonTGMeasurementTool.cxx.

794  {
795  const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
796  if (!m_useDSManager) {
797  SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
798  MuonDetMgr = DetectorManagerHandle.cptr();
799  if (MuonDetMgr == nullptr) { ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object"); }
800  }
801  // Get the messaging service, print where you are
802  ATH_MSG_DEBUG("MuonTGMeasurementTool::nearestDetEl");
803  Identifier nid(0);
804  // check input
805  if (!lay || !parm || !lay->layerType()) return nid;
806 
807  // check compatibility of layer info and required id ? this was already done when associating !
808  Identifier layId(lay->layerType());
809 
810  unsigned int hitType = 0;
811  if (m_idHelperSvc->isMdt(layId)) hitType = 1;
812  if (m_idHelperSvc->isRpc(layId)) hitType = 2;
813  if (m_idHelperSvc->isCsc(layId)) hitType = 3;
814  if (m_idHelperSvc->isTgc(layId)) hitType = 4;
815 
816  if (hitType == 0) {
817  ATH_MSG_DEBUG("unknown hit technology");
818  return nid;
819  }
820 
821  if (hitType == 1) {
822  const MuonGM::MdtReadoutElement* mdtROE = MuonDetMgr->getMdtReadoutElement(layId);
823  if (!mdtROE) return nid;
824  int tMax = mdtROE->getNtubesperlayer();
825  // local position at layer
826  Amg::Vector2D locLay;
827  const Amg::Vector3D& globPos = parm->position();
828  bool onSurface = lay->surfaceRepresentation().globalToLocal(globPos, globPos, locLay);
829  if (!onSurface || measPhi) { return nid; }
830  // nearest tube index
831  double refL = lay->getRef();
832  if (fabs(refL) > 10e6) {
833  double sign = (lay->getRef() > 0) ? 1. : -1.;
834  int dec = int(lay->getRef() / 1.e5);
835  refL = lay->getRef() - dec * 1e5 - 0.5 * (sign + 1) * 1e5;
836  }
837  double dloc = locLay[Trk::locY] - refL + 15.0175;
838  int itube = int(dloc / 30.035) + 1;
839  if (itube < 1 || itube > tMax || dloc < 0.) { return nid; }
840  // tube id
841  Identifier nearId = m_idHelperSvc->mdtIdHelper().channelID(
842  m_idHelperSvc->mdtIdHelper().stationName(layId), m_idHelperSvc->mdtIdHelper().stationEta(layId),
843  m_idHelperSvc->mdtIdHelper().stationPhi(layId), m_idHelperSvc->mdtIdHelper().multilayer(layId),
844  m_idHelperSvc->mdtIdHelper().tubeLayer(layId), itube);
845  // check if position within active volume
846  if (fabs(locLay[Trk::locX]) > 0.5 * mdtROE->getActiveTubeLength(m_idHelperSvc->mdtIdHelper().tubeLayer(layId), itube)) {
847  return nid;
848  }
849  //
850  if (m_idHelperSvc->mdtIdHelper().valid(nearId))
851  return nearId;
852  else
853  return nid;
854  }
855 
856  if (hitType == 2) {
857  // identify phi doublet first
858  Identifier phiId(0);
859  // both phi and eta needs to be within bounds !!!
860  int doubletPhi = 1;
861  bool foundDoubletPhi = false;
862  const MuonGM::RpcReadoutElement* rpcROE = nullptr;
863  while (!foundDoubletPhi && doubletPhi < 3) {
864  Identifier refPhi1 = m_idHelperSvc->rpcIdHelper().channelID(
865  m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
866  m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId),
867  m_idHelperSvc->rpcIdHelper().doubletZ(layId), doubletPhi, m_idHelperSvc->rpcIdHelper().gasGap(layId), 1, 1);
868  rpcROE = MuonDetMgr->getRpcReadoutElement(refPhi1);
869  if (!rpcROE) return nid;
870  if (!m_idHelperSvc->rpcIdHelper().valid(refPhi1)) return nid;
871  int nStripPhi = rpcROE->Nstrips(1);
872 
873  Identifier refPhiN = m_idHelperSvc->rpcIdHelper().channelID(
874  m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
875  m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId),
876  m_idHelperSvc->rpcIdHelper().doubletZ(layId), doubletPhi, m_idHelperSvc->rpcIdHelper().gasGap(layId), 1, nStripPhi);
877  //
878  Amg::Vector2D loc1;
879  const Amg::Vector3D globPos1 = rpcROE->stripPos(refPhi1);
880  bool onSurface1 = rpcROE->surface(refPhi1).globalToLocal(globPos1, globPos1, loc1);
881  Amg::Vector2D locN;
882  const Amg::Vector3D globPosN = rpcROE->stripPos(refPhiN);
883  bool onSurfaceN = rpcROE->surface(refPhiN).globalToLocal(globPosN, globPosN, locN);
884  Amg::Vector2D loc;
885  const Amg::Vector3D& glob = parm->position();
886  bool onSurface = rpcROE->surface(refPhiN).globalToLocal(glob, glob, loc);
887 
888  int strip = 0;
889  if (onSurface1 && onSurfaceN && onSurface) {
890  pitch = (locN[Trk::locX] - loc1[Trk::locX]) / fmax(1, nStripPhi - 1);
891  double dstrip = (loc[Trk::locX] - loc1[Trk::locX]) / pitch + 0.5;
892  strip = dstrip >= 0. ? int(dstrip) + 1 : 0;
893 
894  if (strip > 0 && strip <= nStripPhi) {
895  // correct doublet
896  phiId = m_idHelperSvc->rpcIdHelper().channelID(
897  m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
898  m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId),
899  m_idHelperSvc->rpcIdHelper().doubletZ(layId), doubletPhi, m_idHelperSvc->rpcIdHelper().gasGap(layId), 1, strip);
900  foundDoubletPhi = true;
901 
902  } else {
903  doubletPhi += 1;
904  }
905  } else
906  doubletPhi += 1;
907  }
908  // no valid phi strip - skip search
909  if (!foundDoubletPhi) return nid;
910 
911  // eta strip
912  Identifier etaId(0);
913  int doubletZ = 1;
914  while (doubletZ < 4) {
915  Identifier refEta1 = m_idHelperSvc->rpcIdHelper().channelID(
916  m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
917  m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId), doubletZ, doubletPhi,
918  m_idHelperSvc->rpcIdHelper().gasGap(layId), 0, 1);
919  rpcROE = MuonDetMgr->getRpcReadoutElement(refEta1);
920  if (!rpcROE) return nid;
921  if (!m_idHelperSvc->rpcIdHelper().valid(refEta1)) return nid;
922  int nStrips = rpcROE->Nstrips(0);
923 
924  Identifier refEtaN = m_idHelperSvc->rpcIdHelper().channelID(
925  m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
926  m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId), doubletZ, doubletPhi,
927  m_idHelperSvc->rpcIdHelper().gasGap(layId), 0, nStrips);
928  //
929 
930  Amg::Vector2D loc1;
931  const Amg::Vector3D globPos1 = rpcROE->stripPos(refEta1);
932  bool onSurface1 = rpcROE->surface(refEta1).globalToLocal(globPos1, globPos1, loc1);
933  Amg::Vector2D locN;
934  const Amg::Vector3D globPosN = rpcROE->stripPos(refEtaN);
935  bool onSurfaceN = rpcROE->surface(refEtaN).globalToLocal(globPosN, globPosN, locN);
936  Amg::Vector2D loc;
937  const Amg::Vector3D& glob = parm->position();
938  bool onSurface = rpcROE->surface(refEtaN).globalToLocal(glob, glob, loc);
939 
940  int strip = 0;
941  if (onSurface && onSurface1 && onSurfaceN) {
942  pitch = (locN[Trk::locX] - loc1[Trk::locX]) / fmax(1, nStrips - 1);
943  // strip = int( (refPar->localPosition()[Trk::locX]-loc1[Trk::locX])/pitch+0.5 )+1;
944  double dstrip = (loc[Trk::locX] - loc1[Trk::locX]) / pitch + 0.5;
945  strip = dstrip >= 0. ? int(dstrip) + 1 : 0;
946 
947  if (strip > 0 && strip <= nStrips) {
948  etaId = m_idHelperSvc->rpcIdHelper().channelID(
949  m_idHelperSvc->rpcIdHelper().stationName(layId), m_idHelperSvc->rpcIdHelper().stationEta(layId),
950  m_idHelperSvc->rpcIdHelper().stationPhi(layId), m_idHelperSvc->rpcIdHelper().doubletR(layId), doubletZ, doubletPhi,
951  m_idHelperSvc->rpcIdHelper().gasGap(layId), 0, strip);
952  if (measPhi)
953  return phiId;
954  else
955  return etaId;
956  }
957  }
958  doubletZ += 1;
959  }
960 
961  return nid;
962  }
963 
964  if (hitType == 3) {
965  // ref id
966  Identifier refId = m_idHelperSvc->cscIdHelper().channelID(
967  m_idHelperSvc->cscIdHelper().stationName(layId), m_idHelperSvc->cscIdHelper().stationEta(layId),
968  m_idHelperSvc->cscIdHelper().stationPhi(layId), m_idHelperSvc->cscIdHelper().chamberLayer(layId),
969  m_idHelperSvc->cscIdHelper().wireLayer(layId), measPhi, 1);
970  if (!m_idHelperSvc->cscIdHelper().valid(refId)) return nid;
971  // residual in ref frame
972  const Trk::TrackParameters* refPar = layerToDetEl(lay, parm, refId);
973  if (!refPar) return nid;
974  //
975  const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(refId);
976  if (!cscROE) {
977  delete refPar;
978  return nid;
979  }
980  pitch = cscROE->StripPitch(measPhi);
981  int nStrips = m_idHelperSvc->cscIdHelper().stripMax(refId);
982  if (nStrips < 1) {
983  delete refPar;
984  return nid;
985  }
986 
987  Identifier refIdN = m_idHelperSvc->cscIdHelper().channelID(
988  m_idHelperSvc->cscIdHelper().stationName(refId), m_idHelperSvc->cscIdHelper().stationEta(refId),
989  m_idHelperSvc->cscIdHelper().stationPhi(refId), m_idHelperSvc->cscIdHelper().chamberLayer(refId),
990  m_idHelperSvc->cscIdHelper().wireLayer(refId), measPhi, nStrips);
991  //
992  Amg::Vector3D loc1 = cscROE->surface(refId).transform().inverse() * cscROE->stripPos(refId);
993  Amg::Vector3D locN = cscROE->surface(refIdN).transform().inverse() * cscROE->stripPos(refIdN);
994  int strip = 0;
995 
996  pitch = (locN[0] - loc1[0]) / fmax(1, nStrips - 1);
997  strip = int((refPar->localPosition()[Trk::locX] - loc1[0]) / pitch + 0.5) + 1;
998 
999  delete refPar;
1000  refPar = nullptr;
1001  if (strip > 0 && strip <= nStrips) {
1002  // strip id
1003  Identifier nearId = m_idHelperSvc->cscIdHelper().channelID(
1004  m_idHelperSvc->cscIdHelper().stationName(layId), m_idHelperSvc->cscIdHelper().stationEta(layId),
1005  m_idHelperSvc->cscIdHelper().stationPhi(layId), m_idHelperSvc->cscIdHelper().chamberLayer(layId),
1006  m_idHelperSvc->cscIdHelper().wireLayer(layId), measPhi, strip);
1007  if (fabs(residual(parm, nearId)) > 0.5 * pitch)
1008  ATH_MSG_DEBUG("nearest CSC channel residual too large: " << residual(parm, nearId));
1009  return nearId;
1010  }
1011  }
1012 
1013  if (hitType == 4) { // ref id
1014  if (measPhi && m_idHelperSvc->tgcIdHelper().gasGap(layId) == 2 && m_idHelperSvc->tgcIdHelper().gasGapMax(layId) == 3)
1015  return nid; // no phi strips here
1016  // ref id
1017  Identifier refId = layId;
1018  if (measPhi)
1019  refId = m_idHelperSvc->tgcIdHelper().channelID(
1020  m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
1021  m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, 1);
1022  if (!m_idHelperSvc->tgcIdHelper().valid(refId)) return nid;
1023  // residual in ref frame
1024  const Trk::TrackParameters* refPar = layerToDetEl(lay, parm, refId);
1025  if (!refPar) return nid;
1026  //
1027  const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(layId);
1028  if (!tgcROE) {
1029  delete refPar;
1030  return nid;
1031  }
1032  int nStrips = m_idHelperSvc->tgcIdHelper().channelMax(refId);
1033 
1034  if (nStrips < 1) {
1035  delete refPar;
1036  return nid;
1037  }
1038 
1039  Identifier refIdN = m_idHelperSvc->tgcIdHelper().channelID(
1040  m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
1041  m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, nStrips);
1042  //
1043  Amg::Vector2D loc1;
1044  const Amg::Vector3D glob1 = tgcROE->channelPos(refId);
1045  bool onSurface1 = tgcROE->surface(refId).globalToLocal(glob1, glob1, loc1);
1046  Amg::Vector2D locN;
1047  const Amg::Vector3D globN = tgcROE->channelPos(refIdN);
1048  bool onSurfaceN = tgcROE->surface(refIdN).globalToLocal(globN, globN, locN);
1049 
1050  int strip = 0;
1051  if (onSurface1 && onSurfaceN) {
1052  pitch = (locN[Trk::locX] - loc1[Trk::locX]) / fmax(1, nStrips - 1);
1053  strip = int((refPar->localPosition()[Trk::locX] - loc1[Trk::locX]) / pitch + 0.5) + 1;
1054  } else {
1055  ATH_MSG_DEBUG("local position of boundary elements not retrieved, return 0 ");
1056  delete refPar;
1057  return nid;
1058  }
1059 
1060  if (strip > 0 && strip <= nStrips) {
1061  // check second coordinate for active volume
1062  if (!measPhi && fabs(refPar->localPosition()[Trk::locY]) >
1063  tgcROE->gangCentralWidth(m_idHelperSvc->tgcIdHelper().gasGap(layId), strip)) {
1064  delete refPar;
1065  return nid;
1066  }
1067  Identifier nearId = m_idHelperSvc->tgcIdHelper().channelID(
1068  m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
1069  m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, strip);
1070  Amg::Vector3D stripposition = tgcROE->surface(nearId).transform().inverse() * tgcROE->channelPos(nearId);
1071  Amg::Vector3D localhit = tgcROE->surface(nearId).transform().inverse() * parm->position();
1072 
1073  int plane = m_idHelperSvc->tgcIdHelper().gasGap(nearId);
1074  if (m_idHelperSvc->tgcIdHelper().isStrip(nearId))
1075  pitch = tgcROE->stripPitch(plane, m_idHelperSvc->tgcIdHelper().channel(nearId), localhit[1]);
1076  int last = 0;
1077  while (fabs(stripposition[0] - localhit[0]) > 0.5 * pitch) {
1078  if (stripposition[0] < localhit[0]) {
1079  if (last != -1) {
1080  strip += 1;
1081  last = 1;
1082  } else
1083  break;
1084  } else {
1085  if (last != 1) {
1086  strip -= 1;
1087  last = -1;
1088  } else
1089  break;
1090  }
1091  if (strip < 1 || strip > nStrips) break;
1092  nearId = m_idHelperSvc->tgcIdHelper().channelID(
1093  m_idHelperSvc->tgcIdHelper().stationName(layId), m_idHelperSvc->tgcIdHelper().stationEta(layId),
1094  m_idHelperSvc->tgcIdHelper().stationPhi(layId), m_idHelperSvc->tgcIdHelper().gasGap(layId), measPhi, strip);
1095  stripposition = tgcROE->surface(nearId).transform().inverse() * tgcROE->channelPos(nearId);
1096  localhit = tgcROE->surface(nearId).transform().inverse() * parm->position();
1097  if (m_idHelperSvc->tgcIdHelper().isStrip(nearId))
1098  pitch = tgcROE->stripPitch(plane, m_idHelperSvc->tgcIdHelper().channel(nearId), localhit[1]);
1099  }
1100  delete refPar;
1101  if (strip < 1 || strip > nStrips) return nid;
1102  if (m_idHelperSvc->tgcIdHelper().valid(nearId))
1103  return nearId;
1104  else
1105  return nid;
1106  }
1107  delete refPar;
1108  }
1109  return nid;
1110 }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ residual() [1/4]

double Muon::MuonTGMeasurementTool::residual ( const Trk::Layer *  layer,
const Trk::TrackParameters *  layPar,
const Trk::RIO_OnTrack *  rio 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1238 of file MuonTGMeasurementTool.cxx.

1239  {
1240  double res = 10000.;
1241  if (!layer || !layPar || !rio) return res;
1242 
1243  const Trk::TrackParameters* detElPar = layerToDetEl(layer, layPar, rio->identify());
1244  if (!detElPar) return res;
1245  if (m_idHelperSvc->isMdt(rio->identify())) {
1246  res = fabs(detElPar->localPosition()[Trk::locR] - rio->localParameters()[Trk::locR]);
1247  } else {
1248  res = fabs(detElPar->localPosition()[Trk::locX] - rio->localParameters()[Trk::locX]);
1249  }
1250  delete detElPar;
1251  return res;
1252 }

◆ residual() [2/4]

double Muon::MuonTGMeasurementTool::residual ( const Trk::Layer *  layer,
const Trk::TrackParameters *  layPar,
Identifier  id 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1254 of file MuonTGMeasurementTool.cxx.

1254  {
1255  const MuonGM::MuonDetectorManager* MuonDetMgr = m_muonDetMgr;
1256  if (!m_useDSManager) {
1257  SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
1258  MuonDetMgr = DetectorManagerHandle.cptr();
1259  if (MuonDetMgr == nullptr) {
1260  ATH_MSG_ERROR("Null pointer to the read MuonDetectorManager conditions object");
1261  // return StatusCode::FAILURE;
1262  }
1263  }
1264 
1265  double res = 10000.;
1266  if (!layer || !layPar || !id.get_identifier32().get_compact()) return res;
1267 
1268  const Trk::TrackParameters* detElPar = layerToDetEl(layer, layPar, id);
1269  if (!detElPar) return res;
1270  if (m_idHelperSvc->isMdt(id)) {
1271  res = detElPar->localPosition()[Trk::locR];
1272  } else if (m_idHelperSvc->isRpc(id)) {
1273  const MuonGM::RpcReadoutElement* rpcROE = MuonDetMgr->getRpcReadoutElement(id);
1274  if (rpcROE)
1275  res = detElPar->localPosition()[Trk::locX] -
1276  (detElPar->associatedSurface().transform().inverse() * (rpcROE->stripPos(id)))[Trk::locX];
1277  } else if (m_idHelperSvc->isCsc(id)) {
1278  const MuonGM::CscReadoutElement* cscROE = MuonDetMgr->getCscReadoutElement(id);
1279  if (cscROE)
1280  res = detElPar->localPosition()[Trk::locX] - (detElPar->associatedSurface().transform().inverse() * (cscROE->stripPos(id)))[0];
1281  } else if (m_idHelperSvc->isTgc(id)) {
1282  if (m_idHelperSvc->tgcIdHelper().isStrip(id) && m_idHelperSvc->tgcIdHelper().gasGap(id) == 2 &&
1283  m_idHelperSvc->tgcIdHelper().gasGapMax(id) == 3) {
1284  delete detElPar;
1285  return res; // no phi strips here
1286  }
1287  const MuonGM::TgcReadoutElement* tgcROE = MuonDetMgr->getTgcReadoutElement(id);
1288  if (tgcROE) {
1289  Amg::Vector2D locPos;
1290  const Amg::Vector3D globPos = tgcROE->channelPos(id);
1291  bool onSurface = detElPar->associatedSurface().globalToLocal(globPos, globPos, locPos);
1292  if (onSurface) res = detElPar->localPosition()[Trk::locX] - locPos[Trk::locX];
1293  }
1294  }
1295  delete detElPar;
1296  return res;
1297 }

◆ residual() [3/4]

double Muon::MuonTGMeasurementTool::residual ( const Trk::TrackParameters *  layPar,
const Trk::RIO_OnTrack *  rio 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1216 of file MuonTGMeasurementTool.cxx.

1216  {
1217  double res = 10000.;
1218  if (!layPar || !rio) return res;
1219 
1220  Amg::Vector3D gp = layPar->position();
1221  const Trk::Layer* layer = associatedLayer(rio->identify(), gp);
1222  if (!layer) return res;
1223 
1224  return residual(layer, layPar, rio);
1225 }

◆ residual() [4/4]

double Muon::MuonTGMeasurementTool::residual ( const Trk::TrackParameters *  layPar,
Identifier &  id 
) const
overridevirtual

Implements Muon::IMuonTGMeasTool.

Definition at line 1227 of file MuonTGMeasurementTool.cxx.

1227  {
1228  double res = 10000.;
1229  if (!layPar || !id.get_identifier32().get_compact()) return res;
1230 
1231  Amg::Vector3D gp = layPar->position();
1232  const Trk::Layer* layer = associatedLayer(id, gp);
1233  if (!layer) return res;
1234 
1235  return residual(layer, layPar, id);
1236 }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_alignedMode

Gaudi::Property<bool> Muon::MuonTGMeasurementTool::m_alignedMode {this, "AlignedMode", true}
private

Definition at line 95 of file MuonTGMeasurementTool.h.

◆ m_DetectorManagerKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> Muon::MuonTGMeasurementTool::m_DetectorManagerKey
private
Initial value:
{this, "DetectorManagerKey", "MuonDetectorManager",
"Key of input MuonDetectorManager condition data"}

Definition at line 77 of file MuonTGMeasurementTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> Muon::MuonTGMeasurementTool::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
private

Definition at line 75 of file MuonTGMeasurementTool.h.

◆ m_muonDetMgr

const MuonGM::MuonDetectorManager* Muon::MuonTGMeasurementTool::m_muonDetMgr
private

Definition at line 86 of file MuonTGMeasurementTool.h.

◆ m_rpcProjEta

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_rpcProjEta
private

Definition at line 91 of file MuonTGMeasurementTool.h.

◆ m_rpcProjPhi

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_rpcProjPhi
private

Definition at line 92 of file MuonTGMeasurementTool.h.

◆ m_tgcProjEta

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_tgcProjEta
private

Definition at line 89 of file MuonTGMeasurementTool.h.

◆ m_tgcProjPhi

std::unique_ptr<AmgMatrix(5, 5)> Muon::MuonTGMeasurementTool::m_tgcProjPhi
private

Definition at line 90 of file MuonTGMeasurementTool.h.

◆ m_trackingGeometryReadKey

SG::ReadCondHandleKey<Trk::TrackingGeometry> Muon::MuonTGMeasurementTool::m_trackingGeometryReadKey
private
Initial value:
{this, "TrackingGeometryReadKey", "",
"Key of input TrackingGeometry"}

Definition at line 83 of file MuonTGMeasurementTool.h.

◆ m_trackingGeometrySvc

ServiceHandle<Trk::ITrackingGeometrySvc> Muon::MuonTGMeasurementTool::m_trackingGeometrySvc
private
Initial value:
{this, "TrackingGeometrySvc",
"TrackingGeometrySvc/AtlasTrackingGeometrySvc"}

Definition at line 80 of file MuonTGMeasurementTool.h.

◆ m_useDSManager

Gaudi::Property<bool> Muon::MuonTGMeasurementTool::m_useDSManager {this, "UseDSManager", false}
private

Definition at line 96 of file MuonTGMeasurementTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
MuonGM::MuonDetectorManager::getRpcReadoutElement
const RpcReadoutElement * getRpcReadoutElement(const Identifier &id) const
access via extended identifier (requires unpacking)
Definition: MuonDetDescr/MuonReadoutGeometry/src/MuonDetectorManager.cxx:168
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
Trk::LocalParameters
Definition: LocalParameters.h:98
Trk::PlaneSurface::globalToLocal
virtual bool globalToLocal(const Amg::Vector3D &glob, const Amg::Vector3D &mom, Amg::Vector2D &loc) const override final
Specified for PlaneSurface: GlobalToLocal method without dynamic memory allocation - boolean checks i...
Definition: PlaneSurface.cxx:213
Muon::MuonDriftCircleErrorStrategy
Definition: MuonDriftCircleErrorStrategy.h:15
test_pyathena.eta
eta
Definition: test_pyathena.py:10
Amg::VectorX
Eigen::Matrix< double, Eigen::Dynamic, 1 > VectorX
Dynamic Vector - dynamic allocation.
Definition: EventPrimitives.h:30
Amg::MatrixX
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Definition: EventPrimitives.h:27
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
Trk::locX
@ locX
Definition: ParamDefs.h:37
Trk::locY
@ locY
local cartesian
Definition: ParamDefs.h:38
Muon::MuonTGMeasurementTool::m_useDSManager
Gaudi::Property< bool > m_useDSManager
Definition: MuonTGMeasurementTool.h:96
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
Trk::ParametersBase::position
const Amg::Vector3D & position() const
Access method for the position.
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
Trk::DetachedTrackingVolume::trackingVolume
const TrackingVolume * trackingVolume() const
returns the TrackingVolume
Definition: DetachedTrackingVolume.h:125
Trk::ParametersBase::associatedSurface
virtual const Surface & associatedSurface() const override=0
Access to the Surface associated to the Parameters.
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
Muon::MuonTGMeasurementTool::m_rpcProjEta
std::unique_ptr< AmgMatrix(5, 5)> m_rpcProjEta
Definition: MuonTGMeasurementTool.h:91
Muon::MuonTGMeasurementTool::m_tgcProjPhi
std::unique_ptr< AmgMatrix(5, 5)> m_tgcProjPhi
Definition: MuonTGMeasurementTool.h:90
Muon::MuonTGMeasurementTool::getGeometry
const Trk::TrackingGeometry * getGeometry() const
Definition: MuonTGMeasurementTool.h:98
Trk::Layer::getRef
double getRef() const
get the reference measure
Trk::PrepRawData::localCovariance
const Amg::MatrixX & localCovariance() const
return const ref to the error matrix
CSV_InDetExporter.new
new
Definition: CSV_InDetExporter.py:145
plotBeamSpotVxVal.cov
cov
Definition: plotBeamSpotVxVal.py:201
Muon::MuonTGMeasurementTool::m_muonDetMgr
const MuonGM::MuonDetectorManager * m_muonDetMgr
Definition: MuonTGMeasurementTool.h:86
Muon::TgcClusterOnTrack
Class to represent calibrated clusters formed from TGC strips.
Definition: TgcClusterOnTrack.h:46
Muon::MuonTGMeasurementTool::m_tgcProjEta
std::unique_ptr< AmgMatrix(5, 5)> m_tgcProjEta
Definition: MuonTGMeasurementTool.h:89
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
Muon::MdtDriftCircleStatus
MdtDriftCircleStatus
Enum to represent the 'status' of Mdt measurements e.g.
Definition: MdtDriftCircleStatus.h:25
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
mc.diff
diff
Definition: mc.SFGenPy8_MuMu_DD.py:14
Trk::TrackingVolume::confinedDetachedVolumes
ArraySpan< DetachedTrackingVolume const *const > confinedDetachedVolumes() const
Return detached subVolumes - not the ownership.
Muon::MuonTGMeasurementTool::m_alignedMode
Gaudi::Property< bool > m_alignedMode
Definition: MuonTGMeasurementTool.h:95
MuonGM::RpcReadoutElement
An RpcReadoutElement corresponds to a single RPC module; therefore typicaly a barrel muon station con...
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/RpcReadoutElement.h:55
Trk::RIO_OnTrack
Definition: RIO_OnTrack.h:70
Muon::RpcClusterOnTrack
Class to represent calibrated clusters formed from RPC strips.
Definition: RpcClusterOnTrack.h:35
Trk::locR
@ locR
Definition: ParamDefs.h:44
MuonGM::CscReadoutElement
Definition: CscReadoutElement.h:56
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
const
bool const RAWDATA *ch2 const
Definition: LArRodBlockPhysicsV0.cxx:560
Trk::Surface::center
const Amg::Vector3D & center() const
Returns the center position of the Surface.
SG::VarHandleKey::empty
bool empty() const
Test if the key is blank.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:150
MuonGM::MuonClusterReadoutElement::surface
virtual const Trk::PlaneSurface & surface() const override
access to chamber surface (phi orientation), uses the first gas gap
Definition: MuonClusterReadoutElement.h:123
Muon::MdtStatusDriftTime
@ MdtStatusDriftTime
The tube produced a vaild measurement.
Definition: MdtDriftCircleStatus.h:34
keylayer_zslicemap.strip
strip
Definition: keylayer_zslicemap.py:151
AmgMatrix
#define AmgMatrix(rows, cols)
Definition: EventPrimitives.h:49
Trk::Layer::layerType
int layerType() const
get the Layer coding
similarity
Matrix< Scalar, OtherDerived::RowsAtCompileTime, OtherDerived::RowsAtCompileTime > similarity(const MatrixBase< OtherDerived > &m) const
similarity method : yields ms = m*s*m^T
Definition: AmgMatrixBasePlugin.h:133
Trk::DefinedParameter
std::pair< double, ParamDefs > DefinedParameter
Definition: DefinedParameter.h:27
MuonGM::TgcReadoutElement::channelPos
Amg::Vector3D channelPos(const Identifier &id) const
Returns the position of the active channel (wireGang or strip)
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
Trk::AtaStraightLine
ParametersT< TrackParametersDim, Charged, StraightLineSurface > AtaStraightLine
Definition: Tracking/TrkEvent/TrkParameters/TrkParameters/TrackParameters.h:35
Trk::TrackingVolume::confinedLayers
const LayerArray * confinedLayers() const
Return the subLayer array.
MuonGM::MuonDetectorManager::getTgcReadoutElement
const TgcReadoutElement * getTgcReadoutElement(const Identifier &id) const
access via extended identifier (requires unpacking)
Definition: MuonDetDescr/MuonReadoutGeometry/src/MuonDetectorManager.cxx:247
python.setupRTTAlg.size
int size
Definition: setupRTTAlg.py:39
tubeMax
double tubeMax
Definition: MDT_ResponseTest.cxx:31
A
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
Muon::MuonTGMeasurementTool::match
const Trk::Layer * match(Identifier id, const Trk::Layer *lay) const override
Definition: MuonTGMeasurementTool.cxx:1184
Muon::MuonTGMeasurementTool::associatedLayer
const Trk::Layer * associatedLayer(Identifier id, Amg::Vector3D &gp) const override
Definition: MuonTGMeasurementTool.cxx:1112
AthCommonDataStore
Definition: AthCommonDataStore.h:52
Muon::MuonDriftCircleErrorStrategyInput
std::bitset< 23 > MuonDriftCircleErrorStrategyInput
Definition: MuonDriftCircleErrorStrategy.h:13
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
Trk::Layer::surfaceRepresentation
virtual const Surface & surfaceRepresentation() const =0
Transforms the layer into a Surface representation for extrapolation.
MuonGM::MuonDetectorManager::getMdtReadoutElement
const MdtReadoutElement * getMdtReadoutElement(const Identifier &id) const
access via extended identifier (requires unpacking)
Definition: MuonDetDescr/MuonReadoutGeometry/src/MuonDetectorManager.cxx:204
ParticleGun_EoverP_Config.momentum
momentum
Definition: ParticleGun_EoverP_Config.py:63
MuonGM::MdtReadoutElement
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MdtReadoutElement.h:51
lumiFormat.i
int i
Definition: lumiFormat.py:85
Muon::RpcPrepData
Class to represent RPC measurements.
Definition: RpcPrepData.h:35
TGC
Definition: TgcBase.h:6
Trk::ArraySpan
std::span< T > ArraySpan
Definition: DetachedTrackingVolume.h:34
MuonGM::MuonDetectorManager::getCscReadoutElement
const CscReadoutElement * getCscReadoutElement(const Identifier &id) const
access via extended identifier (requires unpacking)
Definition: MuonDetDescr/MuonReadoutGeometry/src/MuonDetectorManager.cxx:225
Muon::CscPrepData
Class representing clusters from the CSC.
Definition: CscPrepData.h:39
Trk::TrackingVolume::confinedArbitraryLayers
ArraySpan< Layer const *const > confinedArbitraryLayers() const
Return the confined subLayer array.
Trk::theta
@ theta
Definition: ParamDefs.h:66
MuonR4::inverse
CalibratedSpacePoint::Covariance_t inverse(const CalibratedSpacePoint::Covariance_t &mat)
Inverts the parsed matrix.
Definition: MuonSpectrometer/MuonPhaseII/Event/MuonSpacePoint/src/UtilFunctions.cxx:65
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
xAOD::rotation
rotation
Definition: TrackSurface_v1.cxx:15
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
Muon::CscStatusSimple
@ CscStatusSimple
Cluster with non-precision fit.
Definition: CscClusterStatus.h:29
Trk::TrackingGeometry::lowestStaticTrackingVolume
const TrackingVolume * lowestStaticTrackingVolume(const Amg::Vector3D &gp) const
return the lowest static tracking Volume
Definition: TrackingGeometry.cxx:77
MuonGM::MdtReadoutElement::getActiveTubeLength
double getActiveTubeLength(const int tubeLayer, const int tube) const
Muon::CscPrepData::status
CscClusterStatus status() const
Returns the Csc status (position measurement) flag.
Definition: CscPrepData.h:163
res
std::pair< std::vector< unsigned int >, bool > res
Definition: JetGroupProductTest.cxx:14
Amg::transform
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Definition: GeoPrimitivesHelpers.h:156
Muon::MuonTGMeasurementTool::residual
double residual(const Trk::Layer *, const Trk::TrackParameters *, const Trk::RIO_OnTrack *) const override
Definition: MuonTGMeasurementTool.cxx:1238
Muon::MuonDriftCircleErrorStrategy::UnknownStrategy
@ UnknownStrategy
Definition: MuonDriftCircleErrorStrategy.h:17
MuonGM::CscReadoutElement::StripPitch
double StripPitch(int chlayer, int measphi) const
Definition: CscReadoutElement.cxx:265
MuonGM::TgcReadoutElement
A TgcReadoutElement corresponds to a single TGC chamber; therefore typically a TGC station contains s...
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/TgcReadoutElement.h:42
test_pyathena.parent
parent
Definition: test_pyathena.py:15
sign
int sign(int a)
Definition: TRT_StrawNeighbourSvc.h:107
Trk::Surface::normal
virtual const Amg::Vector3D & normal() const
Returns the normal vector of the Surface (i.e.
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
Trk::ParametersBase
Definition: ParametersBase.h:55
MuonGM::TgcReadoutElement::gangCentralWidth
double gangCentralWidth(int gasGap, int gang) const
Returns the length of the central wire in the gang.
Definition: MuonDetDescr/MuonReadoutGeometry/src/TgcReadoutElement.cxx:92
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
Trk::MeasurementBase::localCovariance
const Amg::MatrixX & localCovariance() const
Interface method to get the localError.
Definition: MeasurementBase.h:138
beamspotman.dir
string dir
Definition: beamspotman.py:623
Trk::DECIDED
@ DECIDED
sign of drift radius has been determined
Definition: DriftCircleStatus.h:18
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
MuonGM::nStrips
int nStrips(const MuonGM::TgcReadoutElement &readoutEle, int layer)
Definition: MuonDetDescr/MuonGeoModelTest/src/GeoModelTgcTest.cxx:46
Trk::ParametersCommon::localPosition
Amg::Vector2D localPosition() const
Access method for the local coordinates, local parameter definitions differ for each surface type.
Trk
Ensure that the ATLAS eigen extensions are properly loaded.
Definition: FakeTrackBuilder.h:9
dumpTgcDigiThreshold.isStrip
list isStrip
Definition: dumpTgcDigiThreshold.py:33
Amg
Definition of ATLAS Math & Geometry primitives (Amg)
Definition: AmgStringHelpers.h:19
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
MuonGM::RpcReadoutElement::stripPos
Amg::Vector3D stripPos(const Identifier &id) const
Definition: MuonDetDescr/MuonReadoutGeometry/src/RpcReadoutElement.cxx:177
Trk::BinnedArray::arrayObjects
virtual BinnedArraySpan< T *const > arrayObjects()=0
Return all objects of the Array non-const we can still modify the T.
Muon::MdtDriftCircleOnTrack
This class represents the corrected MDT measurements, where the corrections include the effects of wi...
Definition: MdtDriftCircleOnTrack.h:37
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
Muon::MdtPrepData
Class to represent measurements from the Monitored Drift Tubes.
Definition: MdtPrepData.h:33
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
python.copyTCTOutput.locDir
locDir
Definition: copyTCTOutput.py:113
Trk::MeasurementBase::localParameters
const LocalParameters & localParameters() const
Interface method to get the LocalParameters.
Definition: MeasurementBase.h:132
MuonGM::TgcReadoutElement::stripPitch
double stripPitch(int gasGap, int strip) const
Returns the pitch of the given strip in gasGap i.
Definition: MuonDetDescr/MuonReadoutGeometry/src/TgcReadoutElement.cxx:176
Trk::ParametersBase::momentum
const Amg::Vector3D & momentum() const
Access method for the momentum.
Muon::CscClusterOnTrack
Class to represent the calibrated clusters created from CSC strips.
Definition: CscClusterOnTrack.h:47
a
TList * a
Definition: liststreamerinfos.cxx:10
h
Trk::Surface::globalToLocal
virtual bool globalToLocal(const Amg::Vector3D &glob, const Amg::Vector3D &mom, Amg::Vector2D &loc) const =0
Specified by each surface type: GlobalToLocal method without dynamic memory allocation - boolean chec...
MuonGM::MuonDetectorManager
The MuonDetectorManager stores the transient representation of the Muon Spectrometer geometry and pro...
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonDetectorManager.h:50
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
ref
const boost::regex ref(r_ef)
Trk::PlaneSurface
Definition: PlaneSurface.h:64
Muon::MuonTGMeasurementTool::m_DetectorManagerKey
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager > m_DetectorManagerKey
Definition: MuonTGMeasurementTool.h:77
Muon::MuonTGMeasurementTool::m_rpcProjPhi
std::unique_ptr< AmgMatrix(5, 5)> m_rpcProjPhi
Definition: MuonTGMeasurementTool.h:92
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
Trk::TrackingVolume::confinedVolumes
const TrackingVolumeArray * confinedVolumes() const
Return the subLayer array.
Trk::RIO_OnTrack::identify
Identifier identify() const
return the identifier -extends MeasurementBase
Definition: RIO_OnTrack.h:152
Trk::AtaPlane
ParametersT< TrackParametersDim, Charged, PlaneSurface > AtaPlane
Definition: Tracking/TrkEvent/TrkParameters/TrkParameters/TrackParameters.h:34
Trk::qOverP
@ qOverP
perigee
Definition: ParamDefs.h:67
Muon::TgcPrepData
Class to represent TGC measurements.
Definition: TgcPrepData.h:32
if
if(febId1==febId2)
Definition: LArRodBlockPhysicsV0.cxx:567
Muon::MuonTGMeasurementTool::m_idHelperSvc
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
Definition: MuonTGMeasurementTool.h:75
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
MuonGM::RpcReadoutElement::Nstrips
int Nstrips(bool measphi) const
returns the number of strips for the phi or eta plane
physics_parameters.parameters
parameters
Definition: physics_parameters.py:144
createCablingJSON.doubletPhi
int doubletPhi
Definition: createCablingJSON.py:16
Trk::phi
@ phi
Definition: ParamDefs.h:75
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
merge.status
status
Definition: merge.py:17
MuonGM::CscReadoutElement::stripPos
Amg::Vector3D stripPos(const Identifier &id) const
takes into account internal alignment parameters, hence gives accurate answer
Definition: CscReadoutElement.cxx:227
Muon::CscClusterStatus
CscClusterStatus
Enum to represent the cluster status - see the specific enum values for more details.
Definition: CscClusterStatus.h:23
Muon::MuonTGMeasurementTool::layerToDetEl
const Trk::TrackParameters * layerToDetEl(const Trk::Layer *, const Trk::TrackParameters *, Identifier) const override
Definition: MuonTGMeasurementTool.cxx:74
IdentifierHash
This is a "hash" representation of an Identifier. This encodes a 32 bit index which can be used to lo...
Definition: IdentifierHash.h:25
MuonGM::MdtReadoutElement::getNtubesperlayer
int getNtubesperlayer() const
Returns the number of tubes in each tube layer.
Trk::DetachedTrackingVolume
Definition: DetachedTrackingVolume.h:46
Trk::loc1
@ loc1
Definition: ParamDefs.h:34
Muon::MuonTGMeasurementTool::m_trackingGeometryReadKey
SG::ReadCondHandleKey< Trk::TrackingGeometry > m_trackingGeometryReadKey
Definition: MuonTGMeasurementTool.h:83
Trk::TrackingVolume
Definition: TrackingVolume.h:121
Trk::Surface::transform
const Amg::Transform3D & transform() const
Returns HepGeom::Transform3D by reference.
NSWL1::globalToLocal
Polygon globalToLocal(const Polygon &pol, float z, const Trk::PlaneSurface &surf)
Definition: GeoUtils.cxx:103
Trk::BinnedArraySpan
std::span< T > BinnedArraySpan
Definition: BinnedArray.h:34
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
mag
Scalar mag() const
mag method
Definition: AmgMatrixBasePlugin.h:26
calibdata.tube
tube
Definition: calibdata.py:31
Muon::MuonTGMeasurementTool::m_trackingGeometrySvc
ServiceHandle< Trk::ITrackingGeometrySvc > m_trackingGeometrySvc
Definition: MuonTGMeasurementTool.h:80
MDT
@ MDT
Definition: RegSelEnums.h:31
fitman.k
k
Definition: fitman.py:528
Trk::Layer
Definition: Layer.h:73
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Sun Dec 22 2024 22:03:52 for ATLAS Offline Software by   doxygen 1.8.18