Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
Classes | Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
MuonGMR4::ReadoutGeomCnvAlg Class Reference

#include <ReadoutGeomCnvAlg.h>

Inheritance diagram for MuonGMR4::ReadoutGeomCnvAlg:
Collaboration diagram for MuonGMR4::ReadoutGeomCnvAlg:

Classes

struct  ConstructionCache
 

Public Member Functions

 ~ReadoutGeomCnvAlg ()=default
 
StatusCode execute (const EventContext &ctx) const override
 
StatusCode initialize () override
 
bool isReEntrant () const override
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual bool isClonable () const override
 Specify if the algorithm is clonable. More...
 
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
 
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. More...
 
virtual bool filterPassed (const EventContext &ctx) const
 
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 
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 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
 

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

StatusCode buildStation (const ActsGeometryContext &gctx, const Identifier &stationId, ConstructionCache &cacheObj) const
 builds a station object from readout element. More...
 
StatusCode cloneReadoutVolume (const ActsGeometryContext &gctx, const Identifier &stationId, ConstructionCache &cacheObj, GeoIntrusivePtr< GeoVFullPhysVol > &clonedPhysVol, MuonGM::MuonStation *&station) const
 Clones the fullPhysical volume of the readoutElement and embeds it into the associated station. More...
 
GeoIntrusivePtr< GeoVFullPhysVol > cloneNswWedge (const ActsGeometryContext &gctx, const MuonGMR4::MuonReadoutElement *nswRE, ConstructionCache &cacheObj) const
 Clones the fullPhysicalVolume of the
More...
 
StatusCode buildMdt (const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
 
StatusCode buildRpc (const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
 
StatusCode buildSTGC (const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
 
StatusCode buildMM (const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
 
StatusCode buildTgc (const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
 
StatusCode dumpAndCompare (const ActsGeometryContext &gctx, const MuonGMR4::RpcReadoutElement &refEle, const MuonGM::RpcReadoutElement &testEle) const
 
StatusCode dumpAndCompare (const ActsGeometryContext &gctx, const MuonGMR4::MdtReadoutElement &refEle, const MuonGM::MdtReadoutElement &testEle) const
 
StatusCode dumpAndCompare (const ActsGeometryContext &gctx, const MuonGMR4::MmReadoutElement &refEle, const MuonGM::MMReadoutElement &testEle) const
 
StatusCode dumpAndCompare (const ActsGeometryContext &gctx, const MuonGMR4::TgcReadoutElement &refEle, const MuonGM::TgcReadoutElement &testEle) const
 
StatusCode dumpAndCompare (const ActsGeometryContext &gctx, const MuonGMR4::sTgcReadoutElement &refEle, const MuonGM::sTgcReadoutElement &testEle) const
 
StatusCode checkIdCompability (const MuonGMR4::MuonReadoutElement &refEle, const MuonGM::MuonReadoutElement &testEle) 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::WriteCondHandleKey< MuonGM::MuonDetectorManager > m_writeKey {this, "WriteKey", "MuonDetectorManager"}
 
SG::ReadCondHandleKeyArray< ActsTrk::DetectorAlignStore > m_alignStoreKeys {this, "AlignmentKeys", {}, "Alignment key"}
 
Gaudi::Property< bool > m_checkGeo {this, "checkGeo", false, "Checks the positions of the sensors"}
 
const MuonGMR4::MuonDetectorManager * m_detMgr {nullptr}
 
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
 
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

Definition at line 32 of file ReadoutGeomCnvAlg.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ ~ReadoutGeomCnvAlg()

MuonGMR4::ReadoutGeomCnvAlg::~ReadoutGeomCnvAlg ( )
default

Member Function Documentation

◆ buildMdt()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::buildMdt ( const ActsGeometryContext &  gctx,
ConstructionCache &  cacheObj 
) const
private

Access the B-Line and As-built parameters

Build the mother station

1 cm is added as safety margin to the Mdt multilayer envelope

Need to check how to obtain this parameter from the new geometry newElement->m_cutoutShift;

Determine the tube length's

Define the tube staggering

Add the readout element to the detector manager

Definition at line 648 of file ReadoutGeomCnvAlg.cxx.

648  {
650  SubDetAlignment alignItr = gctx.getStore(ActsTrk::DetectorType::Mdt);
651  const MdtAlignmentStore* alignStore = alignItr ?
652  static_cast<const MdtAlignmentStore*>(alignItr->internalAlignment.get()) : nullptr;
653 
654  const std::vector<const MuonGMR4::MdtReadoutElement*> mdtReadOuts{m_detMgr->getAllMdtReadoutElements()};
655  ATH_MSG_INFO("Copy "<<mdtReadOuts.size()<<" Mdt readout elements to the legacy system");
656  for (const MuonGMR4::MdtReadoutElement* copyMe : mdtReadOuts) {
657  const Identifier reId = copyMe->identify();
658  ATH_MSG_DEBUG("Translate "<<m_idHelperSvc->toStringDetEl(reId));
660  GeoIntrusivePtr<GeoVFullPhysVol> physVol{};
661  MuonGM::MuonStation* station{nullptr};
662  ATH_CHECK(cloneReadoutVolume(gctx,reId, cacheObj, physVol, station));
663  if (copyMe->multilayer() == 1) {
664  const MuonGMR4::MdtReadoutElement* otherRE = copyMe->complementaryRE();
665  const double height = std::max(copyMe->moduleHeight(), otherRE->moduleHeight()) -
666  (copyMe->tubePitch() - 2. * copyMe->tubeRadius());
667 
668  const double modHalTHickO{0.5*otherRE->moduleThickness()},
669  modHalfThick{-0.5*copyMe->moduleThickness()};
670 
671  const double thickness = ( (otherRE->asBuiltRefFrame()*(modHalTHickO* Amg::Vector3D::UnitX())) -
672  (copyMe->asBuiltRefFrame()*(modHalfThick* Amg::Vector3D::UnitX()))).z();
673  if (copyMe->isBarrel()) {
674  station->setMdtZsize(height);
675  station->setMdtRsize(thickness);
676  } else {
677  station->setMdtRsize(height);
678  station->setMdtZsize(thickness);
679  }
680  }
681  const MuonGMR4::MdtReadoutElement::parameterBook& pars{copyMe->getParameters()};
682  auto newElement = std::make_unique<MuonGM::MdtReadoutElement>(physVol,
683  m_idHelperSvc->stationNameString(reId),
684  cacheObj.detMgr.get());
685  newElement->setIdentifier(reId);
686  newElement->setMultilayer(copyMe->multilayer());
687  newElement->setNMdtInStation(m_idHelperSvc->mdtIdHelper().multilayerMax(reId));
688  // cppcheck-suppress invalidLifetime; ok: mgr took ownership.
689  newElement->setParentMuonStation(station);
690 
692  newElement->setLongSsize(2*pars.longHalfX - 1.*Gaudi::Units::cm);
693  newElement->setSsize(2*pars.shortHalfX - 1.*Gaudi::Units::cm);
694  newElement->setLongRsize(2*pars.halfY);
695  newElement->setRsize(2*pars.halfY);
696  newElement->setZsize(2*pars.halfHeight);
697  newElement->setLongZsize(2*pars.halfHeight);
698 
699  newElement->m_nlayers = copyMe->numLayers();
700  newElement->m_ntubesperlayer = copyMe->numTubesInLay();
701  newElement->m_deadlength = pars.deadLength;
702  newElement->m_endpluglength = pars.endPlugLength;
703  newElement->m_innerRadius = pars.tubeInnerRad;
704  newElement->m_tubeWallThickness = pars.tubeWall;
705  newElement->m_tubepitch = pars.tubePitch;
708 
710  const MuonGMR4::MdtTubeLayer& tubeLay{*pars.tubeLayers[0]};
711  unsigned int step{1};
712  double lastLength{2.*tubeLay.uncutHalfLength(1)};
713  for (unsigned tube = 0; tube < copyMe->numTubesInLay(); ++tube) {
714  const double currLength = 2.*tubeLay.uncutHalfLength(tube);
715  ATH_MSG_VERBOSE(m_idHelperSvc->toString(copyMe->identify())<< ", tube "<<tube<<", length: "<<currLength);
716  if (std::abs(lastLength - currLength) > std::numeric_limits<float>::epsilon() ||
717  tube == copyMe->numTubesInLay() -1) {
718  newElement->m_tubelength[step-1] = lastLength;
719  newElement->m_tubelength[step] = currLength;
720  if (step == 1) {
721  newElement->m_ntubesinastep = tube;
722  }
723  lastLength = currLength;
724  ++step;
725  }
726  }
727  newElement->m_nsteps = step;
728 
730  double xOffSet{pars.halfY}, yOffSet{pars.halfHeight};
731  if (newElement->barrel()) std::swap(xOffSet, yOffSet);
732  for (unsigned lay = 1; lay <= copyMe->numLayers(); ++lay) {
733  const IdentifierHash tubeHash{copyMe->measurementHash(lay, 1)};
734  const Amg::Vector3D locTube = copyMe->localTubePos(tubeHash);
735  newElement->m_firstwire_x[lay-1] = locTube.z() + xOffSet;
736  newElement->m_firstwire_y[lay-1] = locTube.x() + yOffSet;
737  }
738  MdtAlignmentStore::chamberDistortions distort = alignStore ? alignStore->getDistortion(reId) :
739  MdtAlignmentStore::chamberDistortions{};
740 
741  if (!station->hasMdtAsBuiltParams()){
742  station->setMdtAsBuiltParams(distort.asBuilt);
743  }
744  if (!station->hasBLines()){
745  station->setBline(distort.bLine);
746  }
747  const Amg::Vector3D refPoint = copyMe->bLineReferencePoint();
748  station->setBlineFixedPointInAmdbLRS(refPoint.x(), refPoint.y(), refPoint.z());
749 
750  newElement->geoInitDone();
751  newElement->setBLinePar(distort.bLine);
752  newElement->fillCache();
753  ATH_CHECK(dumpAndCompare(gctx, *copyMe, *newElement));
755  cacheObj.detMgr->addMdtReadoutElement(std::move(newElement));
756  }
757  return StatusCode::SUCCESS;
758 }

◆ buildMM()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::buildMM ( const ActsGeometryContext &  gctx,
ConstructionCache &  cacheObj 
) const
private

Loop over the gas gaps & efine the

Definition at line 479 of file ReadoutGeomCnvAlg.cxx.

479  {
480 
481  SubDetAlignment alignItr = gctx.getStore(ActsTrk::DetectorType::Mm);
482  const auto alignStore = alignItr ?
483  static_cast<const MmAlignmentStore*>(alignItr->internalAlignment.get()) : nullptr;
484 
485  const std::vector<const MuonGMR4::MmReadoutElement*> mmReadouts{m_detMgr->getAllMmReadoutElements()};
486  ATH_MSG_INFO("Copy "<<mmReadouts.size()<<" Mm readout elements to the legacy system");
487 
488  for (const MuonGMR4::MmReadoutElement* copyMe : mmReadouts) {
489  const Identifier reId = copyMe->identify();
490  GeoIntrusivePtr<GeoVFullPhysVol> physVol{cloneNswWedge(gctx, copyMe, cacheObj)};
491  auto newRE = std::make_unique<MuonGM::MMReadoutElement>(physVol,
492  m_idHelperSvc->stationNameString(reId),
493  copyMe->stationEta(),
494  copyMe->stationPhi(),
495  copyMe->multilayer(), cacheObj.detMgr.get(),
496  alignStore ? alignStore->passivation : nullptr);
498  for (unsigned int gasGap = 0; gasGap < copyMe->nGasGaps(); ++gasGap) {
499  const MuonGMR4::StripLayer& stripLayer{copyMe->stripLayer(MuonGMR4::MmReadoutElement::createHash(gasGap +1, 0))};
500  const MuonGMR4::StripDesign& designFrom{stripLayer.design()};
501 
502  newRE->m_Xlg[gasGap] = stripLayer.toOrigin() *
503  Amg::getRotateZ3D(-designFrom.stereoAngle()) *
504  Amg::getRotateY3D(90. * Gaudi::Units::deg);
505  ATH_MSG_VERBOSE("Layer transform "<<gasGap<<" "<<GeoTrf::toString(newRE->m_Xlg[gasGap], true));
506 
507  MuonGM::MuonChannelDesign& designTo{newRE->m_etaDesign[gasGap]};
508  designTo.defineTrapezoid(designFrom.shortHalfHeight(),
509  designFrom.longHalfHeight(),
510  designFrom.halfWidth(),
511  designFrom.stereoAngle());
512  designTo.type = MuonGM::MuonChannelDesign::ChannelType::etaStrip;
513  designTo.detType = MuonGM::MuonChannelDesign::DetType::MM;
514  designTo.inputPitch = designFrom.stripPitch();
515  designTo.inputWidth = designTo.inputPitch * std::cos(designTo.stereoAngle());
516  designTo.nMissedBottomEta = designTo.nMissedBottomStereo = designFrom.firstStripNumber() - 1;
517  designTo.totalStrips = designFrom.numStrips();
518  designTo.nch = designFrom.numStrips();
519 
520  designTo.setFirstPos(designFrom.firstStripPos().x() + 0.5*designTo.inputPitch);
521  }
522 
523  newRE->fillCache();
524  if (alignStore && alignStore->getBLine(reId)) {
525  newRE->setBLinePar(*alignStore->getBLine(reId));
526  }
527  ATH_CHECK(dumpAndCompare(gctx, *copyMe, *newRE));
528  cacheObj.detMgr->addMMReadoutElement(std::move(newRE));
529  }
530  return StatusCode::SUCCESS;
531 }

◆ buildRpc()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::buildRpc ( const ActsGeometryContext &  gctx,
ConstructionCache &  cacheObj 
) const
private

Define the dimensions

Hack to assign the proper strip positions for REs having doubletPhi =2 in their Identifier

Resignment in face of the rpc readout geometry. Overwrite all the parameters

Definition at line 279 of file ReadoutGeomCnvAlg.cxx.

279  {
280 
281  const std::vector<const MuonGMR4::RpcReadoutElement*> readoutEles = m_detMgr->getAllRpcReadoutElements();
282  ATH_MSG_INFO("Copy "<<readoutEles.size()<<" Rpc readout elements to the legacy system");
283  const RpcIdHelper& idHelper{m_idHelperSvc->rpcIdHelper()};
284  for (const MuonGMR4::RpcReadoutElement* copyMe : readoutEles) {
285  const Identifier reId = copyMe->identify();
286  const MuonGMR4::RpcReadoutElement::parameterBook& pars{copyMe->getParameters()};
287  GeoIntrusivePtr<GeoVFullPhysVol> physVol{};
288  MuonGM::MuonStation* station{nullptr};
289  ATH_CHECK(cloneReadoutVolume(gctx, reId, cacheObj, physVol, station));
290  auto newElement = std::make_unique<MuonGM::RpcReadoutElement>(physVol,
291  m_idHelperSvc->stationNameString(reId),
292  1, 1, false, cacheObj.detMgr.get());
293  const bool aSide{copyMe->stationEta() > 0};
294  newElement->setDoubletPhi(copyMe->doubletPhi());
295  newElement->setDoubletR(copyMe->doubletR());
296  newElement->setDoubletZ(copyMe->doubletZ());
297  newElement->setIdentifier(reId);
298  newElement->setParentMuonStation(station);
299  station->addMuonReadoutElementWithAlTransf(newElement.get(), nullptr, station->nMuonReadoutElements());
300 
302  newElement->setLongZsize(2.*pars.halfLength);
303  newElement->setLongSsize(2.*pars.halfWidth);
304  newElement->setLongRsize(2.*pars.halfThickness);
305  newElement->setZsize(2.*pars.halfLength);
306  newElement->setSsize(2.*pars.halfWidth);
307  newElement->setRsize(2.*pars.halfThickness);
308 
309  newElement->m_nlayers = copyMe->nGasGaps();
310  newElement->m_phistripwidth = copyMe->stripPhiWidth();
311  newElement->m_etastripwidth = copyMe->stripEtaWidth();
312  newElement->m_phistrippitch = copyMe->stripPhiPitch();
313  newElement->m_etastrippitch = (aSide > 0 ? 1. : -1.) *copyMe->stripEtaPitch();
314  newElement->m_phistriplength = copyMe->stripPhiLength();
315  newElement->m_etastriplength = copyMe->stripEtaLength();
316 
317  newElement->m_nphistripsperpanel = copyMe->nPhiStrips();
318  newElement->m_netastripsperpanel = copyMe->nEtaStrips();
319  newElement->m_nphistrippanels = copyMe->nPhiPanels();
320  newElement->m_hasDEDontop = true;
321  newElement->m_descratzneg = false;
322 
323  std::vector<Identifier> gapIds{};
324  for (unsigned int gasGap = 1; gasGap <= copyMe->nGasGaps(); ++gasGap) {
325  for (int doubPhi = copyMe->doubletPhiMax(); doubPhi >= copyMe->doubletPhi(); --doubPhi) {
326  for (bool measPhi : {false, true}) {
327  if (measPhi && copyMe->nPhiStrips()==0) continue;
328  const int channel = 1;
329  const Identifier gapId = idHelper.channelID(copyMe->identify(),
330  copyMe->doubletZ(),
331  doubPhi, gasGap, measPhi,
332  channel);
333 
334  gapIds.push_back(gapId);
335  const Amg::Vector3D locStripPos = copyMe->globalToLocalTrans(gctx) * copyMe->stripPosition(gctx, gapId);
336  ATH_MSG_VERBOSE("GasGap "<<m_idHelperSvc->toString(gapId)<<", local strip position: "<<Amg::toString(locStripPos));
337  newElement->m_gasGap_xPos[gasGap -1] = locStripPos.x();
340  const int dbPIdx = copyMe->doubletPhi() == 2 ? 1 : doubPhi;
341  if (measPhi) {
342  newElement->m_first_phistrip_s[dbPIdx -1] = locStripPos.y();
343  newElement->m_phistrip_z = locStripPos.z();
344  } else{
345  newElement->m_first_etastrip_z = locStripPos.z();
346  newElement->m_etastrip_s[dbPIdx-1] = locStripPos.y();
347  }
348  }
349  }
350  }
351  newElement->fillCache();
353  for (const Identifier& gapId : gapIds) {
354  const int surfaceHash = newElement->surfaceHash(gapId);
355  const int layerHash = newElement->layerHash(gapId);
356  const Amg::Transform3D& refTrf{copyMe->localToGlobalTrans(gctx, gapId)};
357  newElement->m_surfaceData->m_layerTransforms[surfaceHash] = refTrf;
358  newElement->m_surfaceData->m_layerCenters[layerHash] = refTrf.translation();
359  newElement->m_surfaceData->m_layerNormals[layerHash] = refTrf.linear() * Amg::Vector3D::UnitZ();
360  }
361  ATH_CHECK(dumpAndCompare(gctx, *copyMe, *newElement));
362  cacheObj.detMgr->addRpcReadoutElement(std::move(newElement));
363  }
364  return StatusCode::SUCCESS;
365 }

◆ buildStation()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::buildStation ( const ActsGeometryContext &  gctx,
const Identifier &  stationId,
ConstructionCache &  cacheObj 
) const
private

builds a station object from readout element.

The parent PhysVol of the readoutElement is interpreted as embedding station volume and all children which are not fullPhysical volumes are attached to the copied clone.

Fetch the readout element to get its parent volume

Retrieve the full phyiscal volume

Copy the full physical volume of the muon station

Make sure to copy all the children from the original tree that're not FullPhysVols -> represent They represent the passive material inside the station and are needed for the TrackinGeometry building

To create the muon station, we need to extract the dimensions --> Recieve the edge points from the shapes

Transformation to reach from the alignable point to the Muon station

Skip the full physical volumes as they represent the readout elements

Fetch the transform of the detector element which is AlignableNode x Station x [relative pos in station]

S / R / Z size

S / R / Z size (long)

Add the physical volume to the world

Definition at line 144 of file ReadoutGeomCnvAlg.cxx.

146  {
147  const std::string stName{m_idHelperSvc->stationNameString(stationId)};
148  const int stEta{m_idHelperSvc->stationEta(stationId)};
149  const int stPhi{m_idHelperSvc->stationPhi(stationId)};
150  MuonGM::MuonStation* station = cacheObj.detMgr->getMuonStation(stName, stEta, stPhi);
151  if (station) {
152  ATH_MSG_DEBUG("Station "<<stName<<" "<<stEta<<" "<<stPhi<<" already exists.");
153  return StatusCode::SUCCESS;
154  }
156  const MuonGMR4::MuonReadoutElement* copyMe = m_detMgr->getReadoutElement(stationId);
157 
159  const GeoVFullPhysVol* readOutVol = copyMe->getMaterialGeom();
160  PVConstLink parentVolume = readOutVol->getParent(); // This is the physical volume that contains the readOutVol as a child. This is a const link to an existing physical volume
161  cacheObj.translatedStations.insert(parentVolume);
163  PVLink parentPhysVol{make_intrusive<GeoFullPhysVol>(parentVolume->getLogVol())}; // This is a mutable pointer. This creates a new physical volume using the same logical structure (shapes and materials but not position).
164 
167  const std::vector<GeoChildNodeWithTrf> children = getChildrenWithRef(parentVolume, false); // we get the list of child nodes attached to the parent volume
168  double minX{1.e9}, maxX{-1.e9}, minY1{1.e9}, maxY1{-1.e9}, minY2{1.e9}, maxY2{-1.e9}, minZ{1.e9}, maxZ{-1.e9};
169  for (const GeoChildNodeWithTrf& child : children) {
170  std::vector<Amg::Vector3D> edges = getPolyShapeEdges(child.volume->getLogVol()->getShape(), //we are getting corner points of the edges
171  readOutVol->getX().inverse() * child.transform); //getX() returns the transformation from readout physical volume ref to parent physical volume ref. We take the inverse and we chain it with child.transform, which is from child ref to parent ref
172  for (const Amg::Vector3D& edge : edges) {
173  minX = std::min(minX, edge.x());
174  maxX = std::max(maxX, edge.x());
175  minZ = std::min(minZ, edge.z());
176  maxZ = std::max(maxZ, edge.z());
177  if (edge.z() < 0) {
178  minY1 = std::min(minY1, edge.y());
179  maxY1 = std::max(maxY1, edge.y());
180  } else {
181  minY2 = std::min(minY2, edge.y());
182  maxY2 = std::max(maxY2, edge.y());
183  }
184  }
185  }
188  const double shortS = (maxY1 - minY1);
189  const double longS = (maxY2 - minY2);
190  const double lengthR = (maxX - minX);
191  const double lengthZ = (maxZ - minZ);
192 
193  const GeoAlignableTransform* alignTrf{copyMe->alignableTransform()};
195  const Amg::Transform3D stationTransform = alignTrf->getDefTransform().inverse()*parentVolume->getX();
196 
197  for (const GeoChildNodeWithTrf& child : children) {
199  const GeoVPhysVol &childVolRef = *child.volume;
200  if (typeid(childVolRef) == typeid(GeoFullPhysVol)) {
201  continue;
202  }
203  // Add the beam lines / foams inside the station volume
204  PVLink childVol = const_pointer_cast<GeoVPhysVol>(child.volume);
205  parentPhysVol->add(cacheObj.newIdTag());
206  parentPhysVol->add(cacheObj.makeTransform(stationTransform*child.transform));
207  parentPhysVol->add(cloneVolume(childVol));
208  }
210  const Amg::Transform3D alignedTransform = copyMe->localToGlobalTrans(gctx) *
211  (stationTransform * readOutVol->getX()).inverse();
212 
213  ATH_MSG_VERBOSE("stName "<<stName<<","<<stEta<<","<<stPhi<<" -- shortS: "<<shortS<<", longS: "<<longS
214  <<", lengthR: "<<lengthR<<", lengthZ "<<lengthZ
215  <<std::endl<<"AlignableNode: "<<GeoTrf::toString(alignedTransform, true)
216  <<std::endl<<"Station transform: "<<GeoTrf::toString(stationTransform, true)
217  <<std::endl<<"Readout transform: "<<GeoTrf::toString(readOutVol->getX(), true));
218  auto newStation = std::make_unique<MuonGM::MuonStation>(stName,
219  shortS, lengthR, lengthZ,
220  longS, lengthR, lengthZ,
221  stEta, stPhi, false);
222  newStation->setPhysVol(parentPhysVol);
224  cacheObj.world->add(cacheObj.newIdTag());
225  GeoIntrusivePtr<GeoAlignableTransform> trf = make_intrusive<GeoAlignableTransform>(alignedTransform);
226  newStation->setTransform(trf);
227 
228  newStation->setNominalAmdbLRSToGlobal( trf->getTransform());
229 
230  cacheObj.detMgr->addMuonStation(std::move(newStation));
231 
232  cacheObj.world->add(trf);
233  cacheObj.world->add(parentPhysVol);
234 
235  return StatusCode::SUCCESS;
236 }

◆ buildSTGC()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::buildSTGC ( const ActsGeometryContext &  gctx,
ConstructionCache &  cacheObj 
) const
private

Initialize the eta design

Initialize the phi design

To circumvent the yCutout calculations in padCorner function in MuonPadDesign

Definition at line 533 of file ReadoutGeomCnvAlg.cxx.

533  {
534  SubDetAlignment alignItr = gctx.getStore(ActsTrk::DetectorType::sTgc);
535  auto alignStore = alignItr ? static_cast<const sTgcAlignmentStore*>(alignItr->internalAlignment.get()) : nullptr;
536 
537  const std::vector<const MuonGMR4::sTgcReadoutElement*> sTgcReadOuts{m_detMgr->getAllsTgcReadoutElements()};
538  ATH_MSG_INFO("Copy "<<sTgcReadOuts.size()<<" sTgc readout elements to the legacy system");
539 
540  for (const MuonGMR4::sTgcReadoutElement* copyMe : sTgcReadOuts) {
541  const Identifier reId = copyMe->identify();
542  ATH_MSG_DEBUG("Translate readout element "<<m_idHelperSvc->toStringDetEl(reId)<<".");
543  GeoIntrusivePtr<GeoVFullPhysVol> physVol{cloneNswWedge(gctx, copyMe, cacheObj)};
544 
545  auto newRE = std::make_unique<MuonGM::sTgcReadoutElement>(physVol,
546  m_idHelperSvc->stationNameString(reId).substr(1),
547  copyMe->stationEta(),
548  copyMe->stationPhi(),
549  copyMe->multilayer(),
550  cacheObj.detMgr.get());
551 
552  if (alignStore && alignStore->getBLine(reId)) {
553  newRE->setBLinePar(*alignStore->getBLine(reId));
554  }
555  for (unsigned int layer = 1; layer <= copyMe->numLayers(); ++layer) {
556  using channelType = MuonGMR4::sTgcReadoutElement::ReadoutChannelType;
557  using ChannelDesign = MuonGM::MuonChannelDesign;
558  const IdentifierHash layerHash = MuonGMR4::sTgcReadoutElement::createHash(layer,channelType::Strip,0);
559 
560  const MuonGMR4::StripLayer& stripLayer{copyMe->stripLayer(layerHash)};
561  newRE->m_Xlg[layer -1] = stripLayer.toOrigin() * Amg::getRotateY3D(90. * Gaudi::Units::deg) * Amg::getTranslateX3D( layer%2 ? - 0.01 : 0.01 );
562 
563  const MuonGMR4::StripDesign& copyEtaDesign{stripLayer.design()};
564  ATH_MSG_VERBOSE("Layer: "<<layer<<" "<<copyEtaDesign);
566  ChannelDesign& etaDesign{newRE->m_etaDesign[layer-1]};
567  etaDesign.type = ChannelDesign::ChannelType::etaStrip;
568  etaDesign.detType = ChannelDesign::DetType::STGC;
569  if (copyEtaDesign.yCutout()) {
570  etaDesign.defineDiamond(copyEtaDesign.shortHalfHeight(),
571  copyEtaDesign.longHalfHeight(),
572  copyEtaDesign.halfWidth(),
573  copyEtaDesign.yCutout());
574  } else {
575  etaDesign.defineTrapezoid(copyEtaDesign.shortHalfHeight(),
576  copyEtaDesign.longHalfHeight(),
577  copyEtaDesign.halfWidth());
578 
579  }
580  etaDesign.firstPitch = copyEtaDesign.firstStripPos().x() + 0.5*copyEtaDesign.stripPitch() + copyEtaDesign.halfWidth();
581  etaDesign.inputPitch = copyEtaDesign.stripPitch();
582  etaDesign.inputWidth = copyEtaDesign.stripWidth();
583  etaDesign.nch = copyEtaDesign.numStrips();
584  ATH_MSG_VERBOSE(m_idHelperSvc->toStringDetEl(copyMe->identify())<<", layer: "<<layer<<", eta-design: "<< copyEtaDesign);
585  etaDesign.setFirstPos(copyEtaDesign.firstStripPos().x() + 0.5*copyEtaDesign.stripPitch());
587 
588  const MuonGMR4::WireGroupDesign& copyPhiDesign{copyMe->wireDesign(layerHash)};
589 
590  ChannelDesign& phiDesign{newRE->m_phiDesign[layer-1]};
591  phiDesign.type = ChannelDesign::ChannelType::phiStrip;
592  phiDesign.detType = ChannelDesign::DetType::STGC;
593  if (copyPhiDesign.yCutout() == 0.) {
594  phiDesign.defineTrapezoid(copyPhiDesign.shortHalfHeight(),
595  copyPhiDesign.longHalfHeight(),
596  copyPhiDesign.halfWidth());
597  } else {
598  phiDesign.defineDiamond(copyPhiDesign.shortHalfHeight(),
599  copyPhiDesign.longHalfHeight(),
600  copyPhiDesign.halfWidth(),
601  copyPhiDesign.yCutout());
602  }
603  phiDesign.inputPitch = copyPhiDesign.stripPitch();
604  phiDesign.inputWidth = copyPhiDesign.stripWidth();
605  ATH_MSG_VERBOSE(m_idHelperSvc->toStringDetEl(copyMe->identify())<<", layer: "<<layer<<", phi-design: "<< copyPhiDesign);
606  phiDesign.setFirstPos(copyPhiDesign.firstStripPos().x()); // Position of 1st wire, accounts for staggering
607  phiDesign.firstPitch = copyPhiDesign.numWiresInGroup(1); // Number of Wires in 1st group, group staggering
608  phiDesign.groupWidth = copyPhiDesign.numWiresInGroup(2); // Number of Wires normal group
609  phiDesign.nGroups = copyPhiDesign.numStrips(); // Number of Wire Groups
610  phiDesign.wireCutout = copyPhiDesign.wireCutout(); // Size of "active" wire region for digits
611  phiDesign.nch = copyPhiDesign.nAllWires();
612  phiDesign.isConvertedFromPhaseII = true;
613 
614  const MuonGMR4::PadDesign& copyPadDesign{copyMe->padDesign(layerHash)};
615  MuonGM::MuonPadDesign& padDesign{newRE->m_padDesign[layer-1]};
616  padDesign.Length = copyMe->chamberHeight();
617  padDesign.sWidth = copyMe->sChamberLength();
618  padDesign.lWidth = copyMe->lChamberLength();
619  padDesign.Size = 2.*copyPadDesign.halfWidth();
620  padDesign.ysFrame = copyMe->sFrameWidth();
621  padDesign.ylFrame = copyMe->lFrameWidth();
622  padDesign.thickness = copyMe->thickness();
624  if (copyPadDesign.yCutout()) {
625  padDesign.yCutout = copyPadDesign.halfWidth();
626  }
627  padDesign.setR(copyPadDesign.beamlineRadius());
628  padDesign.sPadWidth = 2.*copyPadDesign.shortHalfHeight();
629  padDesign.lPadWidth = 2.*copyPadDesign.longHalfHeight();
630  padDesign.nPadColumns = copyPadDesign.numPadPhi();
631  padDesign.firstPhiPos = copyPadDesign.firstPadPhiDiv();
632  padDesign.inputPhiPitch = copyPadDesign.anglePadPhi();
633  padDesign.PadPhiShift = copyPadDesign.padPhiShift();
634  padDesign.nPadH = copyPadDesign.numPadEta();
635  padDesign.padEtaMax = copyPadDesign.maxPadEta();
636  padDesign.firstRowPos = copyPadDesign.firstPadHeight();
637  padDesign.inputRowPitch = copyPadDesign.padHeight();
638  padDesign.sectorOpeningAngle = copyPadDesign.sectorAngle();
639  padDesign.isConvertedFromPhaseII = true;
640  }
641  newRE->fillCache();
642  ATH_CHECK(dumpAndCompare(gctx, *copyMe, *newRE));
643  cacheObj.detMgr->addsTgcReadoutElement(std::move(newRE));
644 
645  }
646  return StatusCode::SUCCESS;
647 }

◆ buildTgc()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::buildTgc ( const ActsGeometryContext &  gctx,
ConstructionCache &  cacheObj 
) const
private

Define the local gasGap positions

In the sector frame, the gasGap is oriented along the x-axis

Definition at line 368 of file ReadoutGeomCnvAlg.cxx.

368  {
369 
370  std::vector<const MuonGMR4::TgcReadoutElement*> tgcReadouts{m_detMgr->getAllTgcReadoutElements()};
371  std::stable_sort(tgcReadouts.begin(), tgcReadouts.end(),
372  [](const MuonGMR4::TgcReadoutElement* a, const MuonGMR4::TgcReadoutElement* b){
373  return a->stationEta() > b->stationEta();
374  });
375  ATH_MSG_INFO("Copy "<<tgcReadouts.size()<<" Tgc readout elements to the legacy system");
376  const TgcIdHelper& idHelper{m_idHelperSvc->tgcIdHelper()};
377 
378  using TgcReadoutParams = MuonGM::TgcReadoutParams;
379  std::map<std::string, std::shared_ptr<TgcReadoutParams>> readoutParMap{};
380 
381  for (const MuonGMR4::TgcReadoutElement* copyMe: tgcReadouts) {
382  const Identifier reId = copyMe->identify();
383  GeoIntrusivePtr<GeoVFullPhysVol> physVol{};
384  MuonGM::MuonStation* station{nullptr};
385  ATH_CHECK(cloneReadoutVolume(gctx,reId, cacheObj, physVol, station));
386 
387  auto newRE = std::make_unique<MuonGM::TgcReadoutElement>(physVol, m_idHelperSvc->stationNameString(reId),
388  cacheObj.detMgr.get());
389  newRE->setIdentifier(reId);
390  newRE->setParentMuonStation(station);
391 
392  std::shared_ptr<TgcReadoutParams>& readOutPars = readoutParMap[copyMe->chamberDesign()];
393  if (!readOutPars) {
394  using WiregangArray = TgcReadoutParams::WiregangArray;
395  using StripArray = TgcReadoutParams::StripArray;
396  using GasGapIntArray = TgcReadoutParams::GasGapIntArray;
397 
398  std::array<WiregangArray, 3> wires{};
399  GasGapIntArray nWireGangs{}, nStrips{};
400  StripArray botMountings{}, topMountings{};
401  bool stripSet{false};
402  double wirePitch{0.};
403 
404  for (unsigned int gasGap =1; gasGap <= copyMe->nGasGaps(); ++gasGap) {
405  nWireGangs[gasGap -1] = copyMe->numWireGangs(gasGap);
406  nStrips[gasGap -1] = copyMe->numStrips(gasGap);
407  if (nWireGangs[gasGap -1]) {
408  const MuonGMR4::WireGroupDesign& design{copyMe->wireGangLayout(gasGap)};
409  wirePitch = design.stripPitch();
410  WiregangArray& fillMe{wires[gasGap-1]};
411  for (int gang = 1; gang <= design.numStrips(); ++gang) {
412  fillMe[gang -1] = design.numWiresInGroup(gang);
413  }
414  }
415  if (nStrips[gasGap -1] && !stripSet) {
416  const MuonGMR4::RadialStripDesign& design {copyMe->stripLayout(gasGap)};
417  const int nCh = nStrips[gasGap -1];
418  for (int strip = 1; strip <= nCh; ++strip) {
419  botMountings[strip-1] = - design.stripLeftBottom(strip).x();
420  topMountings[strip-1] = - design.stripLeftTop(strip).x();
421  }
422  botMountings[nCh] = - design.stripRightBottom(nCh).x();
423  topMountings[nCh] = - design.stripRightTop(nCh).x();
424 
425  stripSet = true;
426  }
427  }
428  readOutPars = std::make_unique<TgcReadoutParams>(copyMe->chamberDesign(),
429  0, 0, wirePitch,
430  idHelper.stationPhiMax(reId),
431  std::move(nWireGangs),
432  std::move(wires[0]),
433  std::move(wires[1]),
434  std::move(wires[2]),
435  0,
436  std::move(botMountings),
437  std::move(topMountings),
438  std::move(nStrips));
439  }
440 
442  for (unsigned int gasGap = 1; gasGap <= copyMe->nGasGaps(); ++gasGap) {
443  const IdentifierHash layHash{ copyMe->constructHash(0, gasGap, false)};
445  const Amg::Vector3D translation{copyMe->globalToLocalTrans(gctx) * copyMe->center(gctx, layHash)};
446  newRE->setPlaneZ(translation.x(), gasGap);
447  }
448  newRE->setRsize(copyMe->moduleHeight());
449  newRE->setSsize(copyMe->moduleWidthS());
450  newRE->setZsize(copyMe->moduleThickness());
451 
452  newRE->setLongRsize(copyMe->moduleHeight());
453  newRE->setLongSsize(copyMe->moduleWidthL());
454  newRE->setLongZsize(copyMe->moduleThickness());
455 
456  newRE->setReadOutParams(readOutPars);
457  newRE->fillCache();
458  ATH_CHECK(dumpAndCompare(gctx, *copyMe, *newRE));
459  cacheObj.detMgr->addTgcReadoutElement(std::move(newRE));
460 
461  }
462  return StatusCode::SUCCESS;
463 }

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality
overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64 {
65  return 0;
66 }

◆ checkIdCompability()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::checkIdCompability ( const MuonGMR4::MuonReadoutElement &  refEle,
const MuonGM::MuonReadoutElement &  testEle 
) const
private

Definition at line 760 of file ReadoutGeomCnvAlg.cxx.

761  {
762 
763  if (refEle.identify() != testEle.identify()) {
764  ATH_MSG_FATAL("Two different elements are compared "
765  <<m_idHelperSvc->toString(refEle.identify())<<" vs. "
766  <<m_idHelperSvc->toString(testEle.identify()));
767  return StatusCode::FAILURE;
768  }
769  if (refEle.identHash() != testEle.detectorElementHash()) {
770  ATH_MSG_FATAL("The hashes of the two detector elements "<<m_idHelperSvc->toString(refEle.identify())
771  <<" are completely different "<<refEle.identHash()<<" vs. "<<testEle.detectorElementHash());
772  return StatusCode::FAILURE;
773  }
774  return StatusCode::SUCCESS;
775 }

◆ cloneNswWedge()

GeoIntrusivePtr< GeoVFullPhysVol > MuonGMR4::ReadoutGeomCnvAlg::cloneNswWedge ( const ActsGeometryContext &  gctx,
const MuonGMR4::MuonReadoutElement *  nswRE,
ConstructionCache &  cacheObj 
) const
private

Clones the fullPhysicalVolume of the

Definition at line 466 of file ReadoutGeomCnvAlg.cxx.

468  {
469  GeoIntrusivePtr<const GeoVFullPhysVol> readOutVol{copyMe->getMaterialGeom()};
470  cacheObj.translatedStations.insert(readOutVol->getParent());
471 
472  PVLink clonedVol{cloneVolume(const_pointer_cast<GeoVFullPhysVol>(readOutVol))};
473  GeoIntrusivePtr<GeoFullPhysVol> physVol{dynamic_pointer_cast<GeoFullPhysVol>(clonedVol)};
474  cacheObj.world->add(cacheObj.newIdTag());
475  cacheObj.world->add(cacheObj.makeTransform(copyMe->localToGlobalTrans(gctx)));
476  cacheObj.world->add(physVol);
477  return physVol;
478 }

◆ cloneReadoutVolume()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::cloneReadoutVolume ( const ActsGeometryContext &  gctx,
const Identifier &  stationId,
ConstructionCache &  cacheObj,
GeoIntrusivePtr< GeoVFullPhysVol > &  clonedPhysVol,
MuonGM::MuonStation *&  station 
) const
private

Clones the fullPhysical volume of the readoutElement and embeds it into the associated station.

If creations of the needed station fails, failure is returned. The references to the clonedPhysVol & to the station are set if the procedure was successful.

This is a hack to include the BIL Rpcs into the translation. Recall that the BI-RPCs break the station paradigm Hence, in the new description they have their own alignable transform. However, the legacy geometry tries to sort them into the corresponding Mdt station with the same stName, stEta, stPhi. So quite a lot of gymnastics is now needed to place them accordingly into that frame. getX(Mdt) * delta(Mdt) * X = getX(Rpc) * delta(Rpc) X = [getX(Mdt) * delta(Mdt)]^{-1} * getX(Rpc) * delta(Rpc) The second term is the transform of the transform of the Muon station. The second one needs to be rewritten Y = getX(Rpc) * delta(Rpc) * getX(alignNode -> readoutFrame) X = Y * getX(alignNode -> readoutFrame)^{-1}

Clone the detector element with all of its subvolumes

Definition at line 239 of file ReadoutGeomCnvAlg.cxx.

243  {
244 
245  ATH_CHECK(buildStation(gctx, reId, cacheObj));
246  const std::string stName{m_idHelperSvc->stationNameString(reId)};
247  station = cacheObj.detMgr->getMuonStation(stName,
248  m_idHelperSvc->stationEta(reId),
249  m_idHelperSvc->stationPhi(reId));
250 
251  PVLink parentPhysVol{station->getPhysVol()};
252  const MuonGMR4::MuonReadoutElement* copyMe = m_detMgr->getReadoutElement(reId);
253  GeoIntrusivePtr<const GeoVFullPhysVol> readOutVol{copyMe->getMaterialGeom()};
254  parentPhysVol->add(cacheObj.newIdTag());
265  const Amg::Transform3D alignNodeToRE{copyMe->alignableTransform()->getDefTransform().inverse() *
266  readOutVol->getParent()->getX() * readOutVol->getX()};
267  const Amg::Transform3D alignedNode{copyMe->localToGlobalTrans(gctx) * alignNodeToRE.inverse()};
268 
269  const Amg::Transform3D stationTrf{station->getTransform().inverse() * alignedNode};
270 
271  parentPhysVol->add(cacheObj.makeTransform(stationTrf*alignNodeToRE));
273  PVLink clonedVol{cloneVolume(const_pointer_cast<GeoVFullPhysVol>(readOutVol))};
274  physVol = dynamic_pointer_cast<GeoVFullPhysVol>(clonedVol);
275  parentPhysVol->add(physVol);
276  return StatusCode::SUCCESS;
277 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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  }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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; }

◆ dumpAndCompare() [1/5]

StatusCode MuonGMR4::ReadoutGeomCnvAlg::dumpAndCompare ( const ActsGeometryContext &  gctx,
const MuonGMR4::MdtReadoutElement &  refEle,
const MuonGM::MdtReadoutElement &  testEle 
) const
private

Definition at line 825 of file ReadoutGeomCnvAlg.cxx.

827  {
828 
829  if (!m_checkGeo) {
830  return StatusCode::SUCCESS;
831  }
832  ATH_CHECK(checkIdCompability(refEle, testEle));
833 
834  ATH_MSG_VERBOSE("Detector element "<<m_idHelperSvc->toString(refEle.identify())
835  <<std::endl<<GeoTrf::toString(refEle.localToGlobalTrans(gctx))
836  <<std::endl<<GeoTrf::toString(testEle.getMaterialGeom()->getAbsoluteTransform())
837  <<std::endl<<"r-size: "<<testEle.getRsize()<<"/"<<testEle.getLongRsize()
838  <<" s-size: "<<testEle.getSsize()<<"/"<<testEle.getLongSsize()
839  <<" z-size: "<<testEle.getZsize()<<"/"<<testEle.getLongZsize());
840  for (unsigned int lay = 1; lay <= refEle.numLayers(); ++lay){
841  for (unsigned int tube = 1; tube <= refEle.numTubesInLay(); ++tube) {
842  const IdentifierHash tubeHash {refEle.measurementHash(lay,tube)};
843  if (!refEle.isValid(tubeHash)) {
844  ATH_MSG_VERBOSE("SKip layer / tube "<<lay <<","<<tube);
845  continue;
846  }
847  const Amg::Transform3D globToLocal = refEle.globalToLocalTrans(gctx, tubeHash);
848 
849  const Amg::Vector3D refPos = refEle.globalTubePos(gctx, tubeHash);
850  const Amg::Vector3D tubePos = testEle.tubePos(lay, tube);
851 
852 
853  if ( (refPos - tubePos).mag() > Gaudi::Units::micrometer &&
854  (globToLocal*refPos - globToLocal * tubePos).perp() > Gaudi::Units::micrometer) {
855  ATH_MSG_ERROR("Tube positions differ for "<<m_idHelperSvc->toString(refEle.measurementId(tubeHash))
856  <<" reference: "<<GeoTrf::toString(globToLocal*refPos)<<" vs. test: "
857  <<GeoTrf::toString(globToLocal*tubePos) <<" delta: "<<(refPos - tubePos).mag()
858  <<" Transforms "<<std::endl
859  <<" **** "<< GeoTrf::toString(globToLocal.inverse())<<std::endl
860  <<" **** "<< GeoTrf::toString(testEle.transform(lay, tube)));
861  return StatusCode::FAILURE;
862  }
863  ATH_MSG_VERBOSE("Tube positions layer: "<<lay<<", tube: "<<tube
864  <<std::endl<<"reference: "<<GeoTrf::toString(refPos)
865  <<std::endl<<"test: "<<GeoTrf::toString(tubePos)
866  <<std::endl<<testEle.tubeLength(lay, tube)<<"/"<<testEle.getActiveTubeLength(lay, tube)<<"/"
867  <<testEle.getWireLength(lay,tube)<<" vs. "<<refEle.tubeLength(tubeHash)
868  <<"/"<<refEle.activeTubeLength(tubeHash)<<"/"<<refEle.wireLength(tubeHash)
869  <<"/"<<refEle.uncutTubeLength(tubeHash));
870  if (std::abs(testEle.getTubeLengthForCaching(lay,tube) - refEle.uncutTubeLength(tubeHash)) >
871  std::numeric_limits<float>::epsilon() ) {
872  ATH_MSG_FATAL("Different tube length's detected for "<<m_idHelperSvc->toStringDetEl(refEle.identify())
873  << " layer: "<<lay<<", tube: "<<tube<<" -- "<<testEle.getTubeLengthForCaching(lay,tube)<<" (new) vs. "
874  <<refEle.uncutTubeLength(tubeHash)<<" (ref)");
875  return StatusCode::FAILURE;
876  }
877  }
878  }
879 
880  return StatusCode::SUCCESS;
881 }

◆ dumpAndCompare() [2/5]

StatusCode MuonGMR4::ReadoutGeomCnvAlg::dumpAndCompare ( const ActsGeometryContext &  gctx,
const MuonGMR4::MmReadoutElement &  refEle,
const MuonGM::MMReadoutElement &  testEle 
) const
private

Definition at line 777 of file ReadoutGeomCnvAlg.cxx.

779  {
780 
781  if (!m_checkGeo) {
782  return StatusCode::SUCCESS;
783  }
784  ATH_CHECK(checkIdCompability(refEle, testEle));
785 
786  ATH_MSG_VERBOSE("Compare basic readout transforms"<<std::endl
787  <<GeoTrf::toString(testEle.absTransform(),true)<<std::endl
788  <<GeoTrf::toString(refEle.localToGlobalTrans(gctx), true));
789  const MmIdHelper& idHelper{m_idHelperSvc->mmIdHelper()};
790  for (unsigned int gasGap = 1; gasGap <= refEle.nGasGaps(); ++ gasGap) {
791  const Identifier gapId = idHelper.channelID(refEle.identify(), refEle.multilayer(), gasGap, 1);
792 
793  const Amg::Transform3D& refTrf{refEle.localToGlobalTrans(gctx, gapId)};
794  const Amg::Transform3D& testTrf{testEle.transform(gapId)};
795  if (!Amg::isIdentity(refTrf.inverse()*testTrf)) {
796  ATH_MSG_FATAL("The layer "<<m_idHelperSvc->toStringGasGap(gapId)<<" does not transform equally"
797  <<GeoTrf::toString(refTrf, true) <<" vs. "<<GeoTrf::toString(testTrf, true));
798  return StatusCode::FAILURE;
799  }
800  const MuonGMR4::StripDesign& stripDesign{refEle.stripLayer(gapId).design()};
801 
802  for (int strip = stripDesign.firstStripNumber(); strip <= stripDesign.numStrips(); ++strip) {
803  const Identifier stripId = idHelper.channelID(refEle.identify(), refEle.multilayer(), gasGap, strip);
804  const Amg::Vector3D refStripPos{refEle.stripPosition(gctx, stripId)};
805  const Amg::Vector3D refStripDir{refEle.localToGlobalTrans(gctx, refEle.layerHash(stripId)).linear() * Amg::Vector3D::UnitX()};
806 
807  Amg::Vector3D testStripPos{Amg::Vector3D::Zero()};
808  if (!testEle.stripGlobalPosition(stripId, testStripPos)) {
809  ATH_MSG_FATAL("Failed to retrieve strip position "<<m_idHelperSvc->toString(stripId));
810  return StatusCode::FAILURE;
811  }
812  const double dist = refStripDir.dot(refStripPos - testStripPos);
813  if (std::abs(dist) > 10. * Gaudi::Units::micrometer) {
814  ATH_MSG_FATAL("The strip "<<Amg::toString(testStripPos)<<" is not describing the same strip as "
815  <<Amg::toString(refStripPos)<<". Channel "<<m_idHelperSvc->toString(stripId)
816  <<" distance: "<<dist<<" "<<(dist / testEle.m_etaDesign[gasGap -1].inputWidth));
817  return StatusCode::FAILURE;
818  }
819  ATH_MSG_VERBOSE("Channel postion "<<m_idHelperSvc->toString(stripId)<<" match between legacy & new");
820  }
821  }
822  return StatusCode::SUCCESS;
823 }

◆ dumpAndCompare() [3/5]

StatusCode MuonGMR4::ReadoutGeomCnvAlg::dumpAndCompare ( const ActsGeometryContext &  gctx,
const MuonGMR4::RpcReadoutElement &  refEle,
const MuonGM::RpcReadoutElement &  testEle 
) const
private

Definition at line 882 of file ReadoutGeomCnvAlg.cxx.

884  {
885 
886  if (!m_checkGeo) {
887  return StatusCode::SUCCESS;
888  }
889  ATH_CHECK(checkIdCompability(refEle, testEle));
890 
891  ATH_MSG_VERBOSE("Compare basic readout transforms"<<std::endl
892  <<" ref: "<<GeoTrf::toString(refEle.localToGlobalTrans(gctx), true)<<std::endl
893  <<" test: "<<GeoTrf::toString(testEle.absTransform(),true)<<std::endl
894  <<"delta: "<<GeoTrf::toString(testEle.absTransform().inverse()*refEle.localToGlobalTrans(gctx), true ));
895  const RpcIdHelper& idHelper{m_idHelperSvc->rpcIdHelper()};
896  for (unsigned int gasGap = 1; gasGap <= refEle.nGasGaps(); ++gasGap) {
897  for (int doubPhi = refEle.doubletPhi(); doubPhi <= refEle.doubletPhiMax(); ++doubPhi) {
898  for (bool measPhi : {false, true}) {
899  if (measPhi && !refEle.nPhiStrips()) continue;
900  for (int strip = 1; strip <= testEle.Nstrips(measPhi); ++strip) {
901  const Identifier stripId = idHelper.channelID(refEle.identify(),
902  refEle.doubletZ(),
903  doubPhi, gasGap, measPhi, strip);
904 
905  const Amg::Transform3D& refTrans{refEle.localToGlobalTrans(gctx, stripId)};
906  const Amg::Transform3D& testTrans{testEle.transform(stripId)};
907  if (strip == 1 && !Amg::isIdentity(refTrans.inverse()*testTrans)) {
908  ATH_MSG_ERROR("Transformation for "<<m_idHelperSvc->toString(stripId)<<" - "<<refEle.identHash()<<std::endl
909  <<" *** ref: "<<GeoTrf::toString(refTrans)<<std::endl
910  <<" *** test: "<<GeoTrf::toString(testTrans)<<std::endl
911  <<" -> delta: "<<GeoTrf::toString(refTrans.inverse()*testTrans));
912  return StatusCode::FAILURE;
913  }
914 
915  const Amg::Vector3D refStripPos = refEle.stripPosition(gctx, stripId);
916  const Amg::Vector3D testStripPos = testEle.stripPos(stripId);
917  // The threshold here used to be epsilon for a float.
918  // But this was then giving a failure in aarch64,
919  // with a difference of almost exactly 1e-4.
920  // It turned out that that was coming from GeoDeDuplicator,
921  // where to consider two transforms equivalent,
922  // the translations must match to 1e-4. But if the
923  // difference is almost exactly 1e-4, then small FP
924  // differences can be magnified to just about 1e-4.
925  if ((refStripPos - testStripPos).mag() > 2e-4){
926  ATH_MSG_ERROR("Mismatch in strip positions "<<m_idHelperSvc->toString(stripId)
927  <<" ref: "<<Amg::toString(refStripPos)<<" test: "<<Amg::toString(testStripPos)
928  <<" local coordinates -- ref: "<<Amg::toString(refTrans.inverse()*refStripPos)
929  <<" test: "<<Amg::toString(refTrans.inverse()*testStripPos));
930  return StatusCode::FAILURE;
931  }
932  ATH_MSG_VERBOSE("Agreement between new and old geometry for channel "<<m_idHelperSvc->toString(stripId)
933  <<" strip position "<<Amg::toString(refStripPos)
934  <<", local: "<<Amg::toString(refTrans.inverse()*refStripPos));
935  }
936  }
937  }
938  }
939  return StatusCode::SUCCESS;
940 }

◆ dumpAndCompare() [4/5]

StatusCode MuonGMR4::ReadoutGeomCnvAlg::dumpAndCompare ( const ActsGeometryContext &  gctx,
const MuonGMR4::sTgcReadoutElement &  refEle,
const MuonGM::sTgcReadoutElement &  testEle 
) const
private

R4 Geometry function accepts the channelID as input identifier, whereas, R3 geometry needs to be provided padID explicitly.

wireGroup center is defined to be half a wire pitch off in R4 compared to R3 convention

Definition at line 1015 of file ReadoutGeomCnvAlg.cxx.

1017  {
1018  if (!m_checkGeo) {
1019  return StatusCode::SUCCESS;
1020  }
1021  ATH_CHECK(checkIdCompability(refEle, testEle));
1022 
1023  ATH_MSG_VERBOSE("Compare basic readout transforms"<<std::endl
1024  <<GeoTrf::toString(testEle.absTransform(),true)<<std::endl
1025  <<GeoTrf::toString(refEle.localToGlobalTrans(gctx), true));
1026  const sTgcIdHelper& idHelper{m_idHelperSvc->stgcIdHelper()};
1027  for (unsigned int gasGap = 1; gasGap <= refEle.numLayers(); ++gasGap) {
1028  for (int chType : {sTgcIdHelper::sTgcChannelTypes::Pad , sTgcIdHelper::sTgcChannelTypes::Strip, sTgcIdHelper::sTgcChannelTypes::Wire}) {
1029  const Identifier layID = idHelper.channelID(refEle.identify(),
1030  refEle.multilayer(),
1031  gasGap, chType, 1);
1032  const unsigned int numChannel = refEle.numChannels(layID);
1033  constexpr unsigned firstCh = 1;
1034  for (unsigned int channel = firstCh; channel < numChannel ; ++channel) {
1035  const Identifier chID = idHelper.channelID(refEle.identify(),
1036  refEle.multilayer(),
1037  gasGap, chType, channel);
1038 
1039  const Amg::Transform3D& refTrans{refEle.localToGlobalTrans(gctx, chID)};
1040  const Amg::Transform3D& testTrans{testEle.transform(chID)};
1041  if (channel == firstCh && (!Amg::doesNotDeform(testTrans.inverse()*refTrans)
1042  || (testTrans.inverse()*refTrans).translation().perp() > std::numeric_limits<float>::epsilon() ) ) {
1043  ATH_MSG_ERROR(__func__<<"() "<<__LINE__<<" - Transformation for "<<m_idHelperSvc->toString(chID)<<std::endl
1044  <<" *** ref: "<<GeoTrf::toString(refTrans, true)<<std::endl
1045  <<" *** test: "<<GeoTrf::toString(testTrans, true));
1046  return StatusCode::FAILURE;
1047  }
1048  if (chType == sTgcIdHelper::sTgcChannelTypes::Pad) {
1051  bool isValid = false;
1052  const int padEta = refEle.padEta(chID);
1053  const int padPhi = refEle.padPhi(chID);
1054  const Identifier padID = idHelper.padID(refEle.identify(),
1055  refEle.multilayer(),
1056  gasGap, chType, padEta, padPhi, isValid);
1057  if(!isValid) {
1058  ATH_MSG_WARNING("The following pad ID is not valid: " << padID);
1059  }
1060  const Amg::Transform3D& testPadTrans{testEle.transform(padID)};
1061  const Amg::Vector3D refChannelPos = refEle.globalChannelPosition(gctx, chID);
1062  Amg::Vector3D testChannelPos(Amg::Vector3D::Zero());
1063  testEle.stripGlobalPosition(padID, testChannelPos);
1064 
1065  const std::array<Amg::Vector3D,4> refPadCorners = refEle.globalPadCorners(gctx, chID);
1066  std::array<Amg::Vector3D,4> testPadCorners{make_array<Amg::Vector3D, 4>(Amg::Vector3D::Zero())};
1067  testEle.padGlobalCorners(padID, testPadCorners);
1068  for (unsigned int cornerIdx = 0; cornerIdx < refPadCorners.size(); ++cornerIdx) {
1069  const double padCornerDiff = (refPadCorners[cornerIdx] - testPadCorners[cornerIdx]).mag();
1070  if (padCornerDiff - 25. > 1. * Gaudi::Units::micrometer){
1071  ATH_MSG_ERROR("Mismatch in pad Corner " << cornerIdx << ": " <<m_idHelperSvc->toString(padID)
1072  <<" ref: "<<Amg::toString(refPadCorners[cornerIdx])<<" test: "<<Amg::toString(testPadCorners[cornerIdx])
1073  <<" difference: " << padCornerDiff
1074  <<" local coordinates -- ref: "<<Amg::toString(refTrans.inverse()*refPadCorners[cornerIdx])
1075  <<" test: "<<Amg::toString(testPadTrans.inverse()*testPadCorners[cornerIdx]));
1076  return StatusCode::FAILURE;
1077  }
1078  }
1079  const double padChannelDiff = (refChannelPos - testChannelPos).mag();
1080  if (padChannelDiff - 25. > 1. * Gaudi::Units::micrometer){
1081  ATH_MSG_ERROR("Mismatch in pad positions "<<m_idHelperSvc->toString(padID)
1082  <<" ref: "<<Amg::toString(refChannelPos)<<" test: "<<Amg::toString(testChannelPos)
1083  <<" difference: " << padChannelDiff
1084  <<" local coordinates -- ref: "<<Amg::toString(refTrans.inverse()*refChannelPos)
1085  <<" test: "<<Amg::toString(testPadTrans.inverse()*testChannelPos));
1086  return StatusCode::FAILURE;
1087  }
1088  ATH_MSG_VERBOSE("Agreement between new and old geometry for channel "<<m_idHelperSvc->toString(padID)
1089  <<" channel position "<<Amg::toString(refChannelPos));
1090  }
1091  else if (chType == sTgcIdHelper::sTgcChannelTypes::Strip){
1092  const Amg::Vector3D refChannelPos = refEle.globalChannelPosition(gctx, chID);
1093  Amg::Vector3D testChannelPos{Amg::Vector3D::Zero()};
1094  testEle.stripGlobalPosition(chID, testChannelPos);
1095  if ((refChannelPos - testChannelPos).mag() > 1. * Gaudi::Units::micrometer){
1096  ATH_MSG_ERROR("Mismatch in strip positions "<<m_idHelperSvc->toString(chID)
1097  <<" ref: "<<Amg::toString(refChannelPos)<<" test: "<<Amg::toString(testChannelPos)
1098  <<" local coordinates -- ref: "<<Amg::toString(testTrans.inverse()*refChannelPos)
1099  <<" test: "<<Amg::toString(testTrans.inverse()*testChannelPos));
1100  return StatusCode::FAILURE;
1101  }
1102  ATH_MSG_VERBOSE("Agreement between new and old geometry for channel "<<m_idHelperSvc->toString(chID)
1103  <<" channel position "<<Amg::toString(refChannelPos));
1104  } else { // wire
1105  const Amg::Vector3D refChannelPos = refEle.globalChannelPosition(gctx, chID);
1106  Amg::Vector3D testChannelPos{Amg::Vector3D::Zero()};
1107  testEle.stripGlobalPosition(chID, testChannelPos);
1108  Amg::Vector3D localRefPos {testTrans.inverse()*refChannelPos};
1109  Amg::Vector3D localTestPos{testTrans.inverse()*testChannelPos};
1111  if((std::abs(localRefPos.x() -localTestPos.x()) > 1.* Gaudi::Units::micrometer)
1112  || (std::abs(refChannelPos.z() -testChannelPos.z()) > 15.* Gaudi::Units::micrometer)){
1113  ATH_MSG_ERROR("Mismatch in wire positions "<<m_idHelperSvc->toString(chID)
1114  <<" ref: "<<Amg::toString(refChannelPos)<<" test: "<<Amg::toString(testChannelPos)
1115  <<" local coordinates -- ref: "<<Amg::toString(testTrans.inverse()*refChannelPos)
1116  <<" test: "<<Amg::toString(testTrans.inverse()*testChannelPos)
1117  <<"delta X: "<<std::abs(localRefPos.x() -localTestPos.x())
1118  <<", delta Z: "<<std::abs(localRefPos.z() -localTestPos.z()));
1119  return StatusCode::FAILURE;
1120  }
1121  }
1122  }
1123  }
1124  }
1125  return StatusCode::SUCCESS;
1126 }

◆ dumpAndCompare() [5/5]

StatusCode MuonGMR4::ReadoutGeomCnvAlg::dumpAndCompare ( const ActsGeometryContext &  gctx,
const MuonGMR4::TgcReadoutElement &  refEle,
const MuonGM::TgcReadoutElement &  testEle 
) const
private

Definition at line 941 of file ReadoutGeomCnvAlg.cxx.

943  {
944 
945  if (!m_checkGeo) {
946  return StatusCode::SUCCESS;
947  }
948  ATH_CHECK(checkIdCompability(refEle, testEle));
949 
950  const TgcIdHelper& idHelper{m_idHelperSvc->tgcIdHelper()};
951 
952  ATH_MSG_VERBOSE("Detector element "<<m_idHelperSvc->toString(refEle.identify())
953  <<std::endl<<GeoTrf::toString(refEle.localToGlobalTrans(gctx), true)
954  <<std::endl<<GeoTrf::toString(testEle.getMaterialGeom()->getAbsoluteTransform(), true)
955  <<std::endl<<"r-size: "<<testEle.getRsize()<<"/"<<testEle.getLongRsize()
956  <<" s-size: "<<testEle.getSsize()<<"/"<<testEle.getLongSsize()
957  <<" z-size: "<<testEle.getZsize()<<"/"<<testEle.getLongZsize());
958 
959  for (unsigned int gasGap = 1; gasGap <= refEle.nGasGaps(); ++gasGap) {
960  for (bool isStrip : {false, true}) {
961  const IdentifierHash layHash = refEle.constructHash(0, gasGap, isStrip);
962  const Identifier layId = idHelper.channelID(refEle.identify(), gasGap, isStrip, 1);
963  ATH_MSG_VERBOSE("Test layer "<<m_idHelperSvc->toString(layId)<<" "<<refEle.numChannels(layHash)<<" "<<layHash);
964  if (!refEle.numChannels(layHash)) continue;
965  const Amg::Transform3D& refLayerTrf = refEle.localToGlobalTrans(gctx, layHash);
966  const Amg::Transform3D& testLayerTrf = testEle.transform(layId);
967  if (!Amg::isIdentity(refLayerTrf.inverse()* testLayerTrf)) {
968  ATH_MSG_FATAL("The transformations in "<<m_idHelperSvc->toString(layId)
969  <<std::endl<<"ref : "<<GeoTrf::toString(refLayerTrf,true)
970  <<std::endl<<"test: "<<GeoTrf::toString(testLayerTrf,true)
971  <<" are not identical. ");
972  return StatusCode::FAILURE;
973  }
974  ATH_MSG_VERBOSE("Transformations in "<<m_idHelperSvc->toString(layId)
975  <<std::endl<<"ref : "<<GeoTrf::toString(refLayerTrf,true)
976  <<std::endl<<"test: "<<GeoTrf::toString(testLayerTrf,true));
977 
978  for (unsigned int ch = 1; ch <= refEle.numChannels(layHash); ++ch) {
979  const IdentifierHash measHash = refEle.constructHash(ch, gasGap, isStrip);
980  const Identifier measId = refEle.measurementId(measHash);
981  const Amg::Vector3D refChannel = refEle.channelPosition(gctx, measHash);
982  const Amg::Vector3D testChannel = testEle.channelPos(measId);
983  if ((refChannel - testChannel).mag() < std::numeric_limits<float>::epsilon()){
984  continue;
985  }
986  std::stringstream msg{};
987  msg<<"The channel "<<m_idHelperSvc->toString(measId)
988  << " is not at the same position "<<Amg::toString(refChannel)
989  <<" vs. "<<Amg::toString(testChannel)<<". Difference: "
990  <<(refChannel - testChannel).mag();
991  if (!isStrip) {
992  msg<<std::endl<<"*** Test *** - wirePitch: "<<testEle.wirePitch()
993  <<", tot wires "<<testEle.nWires(gasGap)
994  <<", wires to reach "<<testEle.nPitchesToGang(gasGap, ch)
995  <<", wires in gang "<<testEle.nWires(gasGap, ch);
996  const MuonGMR4::WireGroupDesign& design{refEle.wireGangLayout(gasGap)};
997  msg<<std::endl<<"*** Ref *** - wirePitch: "<<design.stripPitch()
998  <<", tot wires "<<testEle.nWires(gasGap)
999  <<", wires to reach "<<design.numPitchesToGroup(ch)
1000  <<", wires in gang "<<design.numWiresInGroup(ch);
1001  } else {
1002  const Amg::Vector3D locRefPos{refLayerTrf.inverse() * refChannel};
1003  const Amg::Vector3D locTestPos{refLayerTrf.inverse()* testChannel};
1004  msg<<std::endl<<"*** Ref **** - "<<Amg::toString(locRefPos)<<std::endl;
1005  msg<<std::endl<<"*** Test **** - "<<Amg::toString(locTestPos)<<std::endl;
1006  }
1007  ATH_MSG_FATAL(msg.str());
1008  return StatusCode::FAILURE;
1009  }
1010 
1011  }
1012  }
1013  return StatusCode::SUCCESS;
1014 }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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; }

◆ execute()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::execute ( const EventContext &  ctx) const
override

Prepare the Geometry context

Ensure that the position & tracking alignment caches are split from the conditions object

Check that for every detector technology there's an DetectorAlignStore in the geometry context Otherwise create an empty one.

Create the cache and populate it with the geoWolds

Build the chamber technologies

Finally add the passives

Move the passives also onto the new world tree...

Final check that all elements are cached properly

Definition at line 76 of file ReadoutGeomCnvAlg.cxx.

76  {
77  SG::WriteCondHandle<MuonGM::MuonDetectorManager> writeHandle{m_writeKey, ctx};
78  if (writeHandle.isValid()) {
79  ATH_MSG_DEBUG("The current readout geometry is still valid.");
80  return StatusCode::SUCCESS;
81  }
82  writeHandle.addDependency(IOVInfiniteRange::infiniteRunLB());
84  ActsGeometryContext geoContext{};
85  using TrackingAlignment = ActsTrk::DetectorAlignStore::TrackingAlignStore;
86  for (const SG::ReadCondHandleKey<ActsTrk::DetectorAlignStore>& key : m_alignStoreKeys) {
87  SG::ReadCondHandle<ActsTrk::DetectorAlignStore> readHandle{key, ctx};
88  if (!readHandle.isValid()) {
89  ATH_MSG_FATAL("Failed to retrieve alignment store "<<key.fullKey());
90  return StatusCode::FAILURE;
91  }
92  writeHandle.addDependency(readHandle);
93  auto alignStore = std::make_unique<ActsTrk::DetectorAlignStore>(**readHandle);
95  if (alignStore->geoModelAlignment) {
96  alignStore->geoModelAlignment->clearPosCache();
97  }
98  alignStore->trackingAlignment = std::make_unique<TrackingAlignment>(alignStore->detType);
99  geoContext.setStore(std::move(alignStore));
100  }
103  std::vector<ActsTrk::DetectorType> presentTechs = m_detMgr->getDetectorTypes();
104  for (const ActsTrk::DetectorType detType : presentTechs) {
105  if (geoContext.getStore(detType)) {
106  continue;
107  }
108  ATH_MSG_WARNING("No external detector alignment has been defined for technology "<<ActsTrk::to_string(detType));
109  geoContext.setStore(std::make_unique<ActsTrk::DetectorAlignStore>(detType));
110  }
112  ConstructionCache cacheObj;
113  cacheObj.detMgr = std::make_unique<MuonGM::MuonDetectorManager>();
114  cacheObj.world = createGeoWorld();
115  cacheObj.detMgr->addTreeTop(cacheObj.world);
117  ATH_CHECK(buildMdt(geoContext, cacheObj));
118  ATH_CHECK(buildTgc(geoContext, cacheObj));
119  ATH_CHECK(buildRpc(geoContext, cacheObj));
120  ATH_CHECK(buildSTGC(geoContext, cacheObj));
121  ATH_CHECK(buildMM(geoContext, cacheObj));
123  std::vector<GeoChildNodeWithTrf> treeTops = getChildrenWithRef(m_detMgr->getTreeTop(0), false);
124 
126  for (const GeoChildNodeWithTrf& treeTop : treeTops) {
127  if (hasStationVolume(treeTop.volume, cacheObj.translatedStations)) continue;
128  ATH_MSG_VERBOSE("Detected passive volume "<<treeTop.nodeName);
129 
130  cacheObj.world->add(const_pointer_cast(treeTop.volume));
131  }
132 
134  if (m_checkGeo) {
135  const std::vector<const MuonGMR4::MuonReadoutElement*> refEles{m_detMgr->getAllReadoutElements()};
136  for (const MuonGMR4::MuonReadoutElement* refEle : refEles) {
137  ATH_CHECK(checkIdCompability(*refEle, *cacheObj.detMgr->getReadoutElement(refEle->identify())));
138  }
139  }
140 
141  ATH_CHECK(writeHandle.record(std::move(cacheObj.detMgr)));
142  return StatusCode::SUCCESS;
143 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

90 {
91  // If we didn't find any symlinks to add, just return the collection
92  // from the base class. Otherwise, return the extended collection.
93  if (!m_extendedExtraObjects.empty()) {
94  return m_extendedExtraObjects;
95  }
96  return BaseAlg::extraOutputDeps();
97 }

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext &  ctx) const
inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96  {
97  return execState( ctx ).filterPassed();
98  }

◆ initialize()

StatusCode MuonGMR4::ReadoutGeomCnvAlg::initialize ( )
override

Definition at line 67 of file ReadoutGeomCnvAlg.cxx.

67  {
68  ATH_CHECK(m_idHelperSvc.retrieve());
69  ATH_CHECK(m_writeKey.initialize());
70  ATH_CHECK(m_alignStoreKeys.initialize());
71  ATH_CHECK(detStore()->retrieve(m_detMgr));
72  return StatusCode::SUCCESS;
73 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable
overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

52 {
53  // Reentrant algorithms are clonable.
54  return true;
55 }

◆ isReEntrant()

bool MuonGMR4::ReadoutGeomCnvAlg::isReEntrant ( ) const
inlineoverride

Definition at line 39 of file ReadoutGeomCnvAlg.h.

39 { return false; }

◆ msg() [1/2]

MsgStream& AthCommonMsg< Gaudi::Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

MsgStream& AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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  }

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool  state,
const EventContext &  ctx 
) const
inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100  {
101  execState( ctx ).setFilterPassed( state );
102  }

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext &  ctx)
overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

77 {
78  return BaseAlg::sysExecute (ctx);
79 }

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in InputMakerBase, and HypoBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107  {
108  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
109 
110  if (sc.isFailure()) {
111  return sc;
112  }
113 
114  ServiceHandle<ICondSvc> cs("CondSvc",name());
115  for (auto h : outputHandles()) {
116  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117  // do this inside the loop so we don't create the CondSvc until needed
118  if ( cs.retrieve().isFailure() ) {
119  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120  return StatusCode::SUCCESS;
121  }
122  if (cs->regHandle(this,*h).isFailure()) {
123  sc = StatusCode::FAILURE;
124  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125  << " with CondSvc");
126  }
127  }
128  }
129  return sc;
130 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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_alignStoreKeys

SG::ReadCondHandleKeyArray<ActsTrk::DetectorAlignStore> MuonGMR4::ReadoutGeomCnvAlg::m_alignStoreKeys {this, "AlignmentKeys", {}, "Alignment key"}
private

Definition at line 124 of file ReadoutGeomCnvAlg.h.

◆ m_checkGeo

Gaudi::Property<bool> MuonGMR4::ReadoutGeomCnvAlg::m_checkGeo {this, "checkGeo", false, "Checks the positions of the sensors"}
private

Definition at line 126 of file ReadoutGeomCnvAlg.h.

◆ m_detMgr

const MuonGMR4::MuonDetectorManager* MuonGMR4::ReadoutGeomCnvAlg::m_detMgr {nullptr}
private

Definition at line 127 of file ReadoutGeomCnvAlg.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_idHelperSvc

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

Definition at line 120 of file ReadoutGeomCnvAlg.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_writeKey

SG::WriteCondHandleKey<MuonGM::MuonDetectorManager> MuonGMR4::ReadoutGeomCnvAlg::m_writeKey {this, "WriteKey", "MuonDetectorManager"}
private

Definition at line 122 of file ReadoutGeomCnvAlg.h.

The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
MuonGMR4::ReadoutGeomCnvAlg::buildSTGC
StatusCode buildSTGC(const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
Definition: ReadoutGeomCnvAlg.cxx:533
MuonGMR4::RpcReadoutElement::doubletPhiMax
int doubletPhiMax() const
Returns the maximum phi panel.
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
make_hlt_rep.pars
pars
Definition: make_hlt_rep.py:90
MuonGMR4::ReadoutGeomCnvAlg::buildMdt
StatusCode buildMdt(const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
Definition: ReadoutGeomCnvAlg.cxx:648
MuonGMR4::ReadoutGeomCnvAlg::dumpAndCompare
StatusCode dumpAndCompare(const ActsGeometryContext &gctx, const MuonGMR4::RpcReadoutElement &refEle, const MuonGM::RpcReadoutElement &testEle) const
Definition: ReadoutGeomCnvAlg.cxx:882
MuonGM::MuonPadDesign
Parameters defining the design of the readout sTGC pads.
Definition: MuonPadDesign.h:40
MuonGMR4::sTgcReadoutElement::multilayer
int multilayer() const
Returns the multilayer of the sTgcReadoutElement.
MuonGM::MuonClusterReadoutElement::transform
virtual const Amg::Transform3D & transform() const override
Return local to global transform.
Definition: MuonClusterReadoutElement.h:124
MuonGMR4::MmReadoutElement
Definition: MmReadoutElement.h:19
STGC
@ STGC
Definition: RegSelEnums.h:39
dumpTgcDigiDeadChambers.gasGap
list gasGap
Definition: dumpTgcDigiDeadChambers.py:33
MuonGMR4::TgcReadoutElement::wireGangLayout
const WireGroupDesign & wireGangLayout(unsigned int gasGap) const
Returns access to the wire group design of the given gasGap [1-3] If the gap does not have a wires an...
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
MuonGMR4::MuonDetectorManager::getDetectorTypes
std::vector< ActsTrk::DetectorType > getDetectorTypes() const
Returns a list of all detector types.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MuonDetectorManager.cxx:143
MuonGMR4::StripDesign
Definition: StripDesign.h:30
sendEI_SPB.ch
ch
Definition: sendEI_SPB.py:35
MuonGMR4::MdtReadoutElement::asBuiltRefFrame
Amg::Transform3D asBuiltRefFrame() const
Returns the transformation to go into the reference frame of the as-buit & b-line model starting from...
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MdtReadoutElement.cxx:211
MuonGMR4::MmReadoutElement::createHash
static IdentifierHash createHash(const int gasGap, const int strip)
MuonGM::MuonReadoutElement::absTransform
const Amg::Transform3D & absTransform() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:210
MuonGMR4::ReadoutGeomCnvAlg::cloneReadoutVolume
StatusCode cloneReadoutVolume(const ActsGeometryContext &gctx, const Identifier &stationId, ConstructionCache &cacheObj, GeoIntrusivePtr< GeoVFullPhysVol > &clonedPhysVol, MuonGM::MuonStation *&station) const
Clones the fullPhysical volume of the readoutElement and embeds it into the associated station.
Definition: ReadoutGeomCnvAlg.cxx:239
MuonGMR4::MmReadoutElement::nGasGaps
unsigned int nGasGaps() const
Returns the number of gas gaps.
MuonGMR4::WireGroupDesign
Definition: WireGroupDesign.h:23
MuonGM::MdtReadoutElement::tubeLength
double tubeLength(const int tubeLayer, const int tube) const
MuonGMR4::sTgcReadoutElement::padEta
unsigned int padEta(const Identifier &measId) const
Returns the Eta index of the pad for the given pad identifier.
MuonGMR4::MdtReadoutElement::numTubesInLay
unsigned int numTubesInLay() const
Returns the number of tubes per layer.
plotting.yearwise_efficiency.channel
channel
Definition: yearwise_efficiency.py:24
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
MuonGMR4::ReadoutGeomCnvAlg::buildStation
StatusCode buildStation(const ActsGeometryContext &gctx, const Identifier &stationId, ConstructionCache &cacheObj) const
builds a station object from readout element.
Definition: ReadoutGeomCnvAlg.cxx:144
ActsTrk::DetectorType
DetectorType
Simple enum to Identify the Type of the ACTS sub detector.
Definition: GeometryDefs.h:17
MdtAlignmentStore::chamberDistortions::bLine
const BLinePar * bLine
Definition: MdtAlignmentStore.h:26
MuonGMR4::ReadoutGeomCnvAlg::buildTgc
StatusCode buildTgc(const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
Definition: ReadoutGeomCnvAlg.cxx:368
perp
Scalar perp() const
perp method - perpenticular length
Definition: AmgMatrixBasePlugin.h:44
MuonGM::MuonStation::getPhysVol
PVConstLink getPhysVol() const
Definition: MuonStation.cxx:398
MuonGM::MuonChannelDesign::defineTrapezoid
void defineTrapezoid(double HalfShortY, double HalfLongY, double HalfHeight)
set the trapezoid dimensions
Definition: MuonChannelDesign.cxx:25
MuonGMR4::sTgcReadoutElement::numChannels
unsigned int numChannels(const Identifier &measId) const
Returns the number of strips / wires / pads in a given gasGap.
ActsTrk::DetectorAlignStore::TrackingAlignStore
Store holding the transfomations used by the Acts algorithms.
Definition: DetectorAlignStore.h:33
TgcIdHelper
Definition: TgcIdHelper.h:50
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
MuonGMR4::MdtReadoutElement::tubeLength
double tubeLength(const IdentifierHash &hash) const
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MdtReadoutElement.cxx:185
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
MuonGMR4::MdtReadoutElement::moduleThickness
double moduleThickness() const
Returns the thickness of the chamber.
MuonGM::GasGapIntArray
TgcReadoutParams::GasGapIntArray GasGapIntArray
Definition: RDBReaderAtlas.cxx:20
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
MuonGMR4::MdtReadoutElement::measurementHash
static IdentifierHash measurementHash(unsigned int layerNumber, unsigned int tubeNumber)
Transform the layer and tube number to the measurementHash.
MuonGMR4::TgcReadoutElement::nGasGaps
unsigned int nGasGaps() const
Returns the number of gasgaps described by this ReadOutElement (usally 2 or 3)
MuonGM::MuonStation::getTransform
Amg::Transform3D getTransform() const
Definition: MuonStation.h:180
MuonGMR4::TgcReadoutElement::channelPosition
Amg::Vector3D channelPosition(const ActsGeometryContext &ctx, const Identifier &measId) const
Returns the center of the measurement channel eta measurement: wire gang center phi measurement: stri...
AthCommonDataStore::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 > renounce(T &h)
Definition: AthCommonDataStore.h:380
MuonGMR4::StripLayer::toOrigin
const Amg::Transform3D & toOrigin() const
Returns the transformation to go from the strip layer center to the origin of the Strip chamber.
MuonGMR4::RpcReadoutElement::nGasGaps
unsigned int nGasGaps() const
Returns the number of gasgaps described by this ReadOutElement (usally 2 or 3)
MuonGMR4::MuonReadoutElement
The MuonReadoutElement is an abstract class representing the geometry representing the muon detector.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/MuonReadoutElement.h:38
MuonGMR4::MdtReadoutElement::activeTubeLength
double activeTubeLength(const IdentifierHash &hash) const
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MdtReadoutElement.cxx:173
MuonGM::MuonReadoutElement::getZsize
double getZsize() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:198
MuonGM::MuonReadoutElement::detectorElementHash
IdentifierHash detectorElementHash() const
Returns the IdentifierHash of the detector element.
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:186
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
deg
#define deg
Definition: SbPolyhedron.cxx:17
MuonGMR4::sTgcReadoutElement::ReadoutChannelType
ReadoutChannelType
ReadoutChannelType to distinguish the available readout channels Pad - pad readout channel Strip - et...
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/sTgcReadoutElement.h:31
MuonGMR4::MdtTubeLayer
Helper struct to retrieve the tube lengths and the tube centers directly from the GeoModel tree.
Definition: MdtTubeLayer.h:28
MuonGMR4::MuonDetectorManager::getAllReadoutElements
std::vector< const MuonReadoutElement * > getAllReadoutElements() const
Returns the list of all detector elements.
MuonGM::MdtReadoutElement::getTubeLengthForCaching
double getTubeLengthForCaching(const int tubeLayer, const int tube) const
Definition: MuonDetDescr/MuonReadoutGeometry/src/MdtReadoutElement.cxx:103
MuonGMR4::MuonReadoutElement::globalToLocalTrans
Amg::Transform3D globalToLocalTrans(const ActsGeometryContext &ctx) const
Transformations to translate between local <-> global coordinates.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MuonReadoutElement.cxx:78
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
isValid
bool isValid(const T &p)
Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine.
Definition: AtlasPID.h:812
MuonGMR4::ReadoutGeomCnvAlg::m_alignStoreKeys
SG::ReadCondHandleKeyArray< ActsTrk::DetectorAlignStore > m_alignStoreKeys
Definition: ReadoutGeomCnvAlg.h:124
ActsTrk::DetectorType::sTgc
@ sTgc
Micromegas (NSW)
MuonGMR4::ReadoutGeomCnvAlg::buildMM
StatusCode buildMM(const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
Definition: ReadoutGeomCnvAlg.cxx:479
sTgcIdHelper::sTgcChannelTypes
sTgcChannelTypes
Definition: sTgcIdHelper.h:190
MuonGM::MuonReadoutElement::getLongSsize
double getLongSsize() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:199
MuonGMR4::ReadoutGeomCnvAlg::m_checkGeo
Gaudi::Property< bool > m_checkGeo
Definition: ReadoutGeomCnvAlg.h:126
MuonGMR4::RpcReadoutElement
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/RpcReadoutElement.h:17
MuonGMR4::TgcReadoutElement::constructHash
static IdentifierHash constructHash(unsigned int measCh, unsigned int gasGap, const bool isStrip)
Constructs the Hash out of the Identifier fields (channel, gasGap, isStrip)
MuonGMR4::ReadoutGeomCnvAlg::buildRpc
StatusCode buildRpc(const ActsGeometryContext &gctx, ConstructionCache &cacheObj) const
Definition: ReadoutGeomCnvAlg.cxx:279
RpcIdHelper
Definition: RpcIdHelper.h:51
keylayer_zslicemap.strip
strip
Definition: keylayer_zslicemap.py:151
MuonGM::WiregangArray
TgcReadoutParams::WiregangArray WiregangArray
Definition: RDBReaderAtlas.cxx:22
std::stable_sort
void stable_sort(std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, std::reverse_iterator< DataModel_detail::iterator< DVL > > end, Compare comp)
Specialization of stable_sort for DataVector/List.
Definition: DVL_algorithms.h:711
MuonGM::sTgcReadoutElement::padGlobalCorners
bool padGlobalCorners(const Identifier &id, std::array< Amg::Vector3D, 4 > &gcorners) const
pad global corners
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/sTgcReadoutElement.h:370
MuonGMR4::RpcReadoutElement::parameterBook
Set of parameters to describe a RPC chamber.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/RpcReadoutElement.h:21
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
MuonGM::TgcReadoutElement::channelPos
Amg::Vector3D channelPos(const Identifier &id) const
Returns the position of the active channel (wireGang or strip)
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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
MuonGMR4::ReadoutGeomCnvAlg::checkIdCompability
StatusCode checkIdCompability(const MuonGMR4::MuonReadoutElement &refEle, const MuonGM::MuonReadoutElement &testEle) const
Definition: ReadoutGeomCnvAlg.cxx:760
Amg::getRotateZ3D
Amg::Transform3D getRotateZ3D(double angle)
get a rotation transformation around Z-axis
Definition: GeoPrimitivesHelpers.h:270
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
MuonGM::MuonReadoutElement::getLongRsize
double getLongRsize() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:200
MuonGM::MuonReadoutElement::getSsize
double getSsize() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:196
MuonGMR4::MdtReadoutElement::measurementId
Identifier measurementId(const IdentifierHash &measHash) const override final
Converts the measurement hash back to the full Identifier.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MdtReadoutElement.cxx:47
MdtAlignmentStore
Helper struct to cache simulatenously the As-built and the BLine corrections of the Mdts for fast acc...
Definition: MdtAlignmentStore.h:19
MuonGMR4::MdtReadoutElement::wireLength
double wireLength(const IdentifierHash &hash) const
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MdtReadoutElement.cxx:188
cm
const double cm
Definition: Simulation/ISF/ISF_FastCaloSim/ISF_FastCaloSimParametrization/tools/FCAL_ChannelMap.cxx:25
MuonGM::MdtReadoutElement::getWireLength
double getWireLength(const int tubeLayer, const int tube) const
MdtAlignmentStore::getDistortion
chamberDistortions getDistortion(const Identifier &detElId) const
Returns a chamber distortion that's cached for the corresponding Mdt chamber element.
Definition: MdtAlignmentStore.h:33
ActsGeometryContext::getStore
AlignmentStorePtr & getStore(const DetectorType type)
Returns the mutable alignable store for the ATLAS detector type (Pixel, Mdt, etc.)
Definition: ActsGeometryContext.h:33
MuonGMR4::ReadoutGeomCnvAlg::cloneNswWedge
GeoIntrusivePtr< GeoVFullPhysVol > cloneNswWedge(const ActsGeometryContext &gctx, const MuonGMR4::MuonReadoutElement *nswRE, ConstructionCache &cacheObj) const
Clones the fullPhysicalVolume of the
Definition: ReadoutGeomCnvAlg.cxx:466
AthCommonReentrantAlgorithm::extraOutputDeps
virtual const DataObjIDColl & extraOutputDeps() const override
Return the list of extra output dependencies.
Definition: AthCommonReentrantAlgorithm.cxx:89
MuonGMR4::MdtReadoutElement::uncutTubeLength
double uncutTubeLength(const IdentifierHash &tubeHash) const
Returns the uncut tube length.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MdtReadoutElement.cxx:179
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
MuonGMR4::sTgcReadoutElement::createHash
static IdentifierHash createHash(const unsigned int gasGap, const unsigned int channelType, const unsigned int channel, const unsigned int wireInGrp=0)
Create a measurement hash from the Identifier fields.
MdtAlignmentStore::chamberDistortions
Helper struct to store the pointer to the Mdt distrotion parameters, namely the As-built and the BLin...
Definition: MdtAlignmentStore.h:25
AthCommonDataStore
Definition: AthCommonDataStore.h:52
Amg::toString
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Definition: GeoPrimitivesToStringConverter.h:40
MuonGMR4::MmReadoutElement::multilayer
int multilayer() const
Returns the multi layer of the element [1-2].
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
MuonGM::TgcReadoutElement::nPitchesToGang
double nPitchesToGang(int gasGap, int gang) const
Returns the number of wire pitches that have to be travelled to reach gang i.
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
MuonGMR4::MmReadoutElement::stripPosition
Amg::Vector3D stripPosition(const ActsGeometryContext &ctx, const Identifier &measId) const
Returns the position of the strip center.
ActsTrk::DetectorType::Mm
@ Mm
Maybe not needed in the migration.
z
#define z
MuonGM::MuonChannelDesign::DetType::MM
@ MM
MuonGMR4::MdtReadoutElement::parameterBook
Set of parameters to describe a MDT chamber.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/MdtReadoutElement.h:26
MuonGM::MuonReadoutElement::getLongZsize
double getLongZsize() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:201
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
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
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
MuonGM::TgcReadoutElement::nWires
int nWires(int gasGap) const
Returns the total number of wires in a given gas gap.
MuonGM::MMReadoutElement::stripGlobalPosition
bool stripGlobalPosition(const Identifier &id, Amg::Vector3D &gpos) const
Definition: MMReadoutElement.h:280
MuonGMR4::RpcReadoutElement::doubletPhi
int doubletPhi() const
Returns the doublet Phi field of the MuonReadoutElement identifier.
MuonGM::MdtReadoutElement::getActiveTubeLength
double getActiveTubeLength(const int tubeLayer, const int tube) const
AthCommonReentrantAlgorithm::sysExecute
virtual StatusCode sysExecute(const EventContext &ctx) override
Execute an algorithm.
Definition: AthCommonReentrantAlgorithm.cxx:76
MuonGMR4::MdtReadoutElement::numLayers
unsigned int numLayers() const
Returns the number of tube layer.
Amg::Transform3D
Eigen::Affine3d Transform3D
Definition: GeoPrimitives.h:46
Amg::doesNotDeform
bool doesNotDeform(const Amg::Transform3D &trans)
Checks whether the linear part of the transformation rotates or stetches any of the basis vectors.
Definition: GeoPrimitivesHelpers.h:383
MuonGMR4::RpcReadoutElement::stripPosition
Amg::Vector3D stripPosition(const ActsGeometryContext &ctx, const Identifier &measId) const
Returns the position of the strip center.
MuonGM::MuonStation
Definition: MuonStation.h:51
IOVInfiniteRange::infiniteRunLB
static EventIDRange infiniteRunLB()
Produces an EventIDRange that is infinite in RunLumi and invalid in Time.
Definition: IOVInfiniteRange.h:39
MuonGMR4::TgcReadoutElement::measurementId
Identifier measurementId(const IdentifierHash &measHash) const override final
Converts the measurement hash back to the full Identifier.
MuonGMR4::RpcReadoutElement::doubletZ
int doubletZ() const
Returns the doublet Z field of the MuonReadoutElement identifier.
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
Amg::isIdentity
bool isIdentity(const Amg::Transform3D &trans)
Checks whether the transformation is the Identity transformation.
Definition: GeoPrimitivesHelpers.h:393
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
python.SystemOfUnits.micrometer
int micrometer
Definition: SystemOfUnits.py:71
MuonGMR4::MuonReadoutElement::identHash
IdentifierHash identHash() const
Returns the Identifier has of the Element that is Identical to the detElHash from the id_helper class...
WriteCalibToCool.swap
swap
Definition: WriteCalibToCool.py:94
ActsGeometryContext
Include the GeoPrimitives which need to be put first.
Definition: ActsGeometryContext.h:27
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
MuonGMR4::MmReadoutElement::stripLayer
const StripLayer & stripLayer(const Identifier &measId) const
MuonGM::nStrips
int nStrips(const MuonGM::TgcReadoutElement &readoutEle, int layer)
Definition: MuonDetDescr/MuonGeoModelTest/src/GeoModelTgcTest.cxx:46
MuonGMR4::MdtReadoutElement::moduleHeight
double moduleHeight() const
Returns the height of the chamber (Distance bottom - topWidth)
MuonGMR4::MuonDetectorManager::getTreeTop
PVConstLink getTreeTop(unsigned int i) const override final
Returns the i-th top node of the MuonSystem trees.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MuonDetectorManager.cxx:133
MdtAlignmentStore::chamberDistortions::asBuilt
const MdtAsBuiltPar * asBuilt
Definition: MdtAlignmentStore.h:27
dumpTgcDigiThreshold.isStrip
list isStrip
Definition: dumpTgcDigiThreshold.py:33
ActsTrk::DetectorType::Mdt
@ Mdt
MuonSpectrometer.
MuonGMR4::ReadoutGeomCnvAlg::m_idHelperSvc
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
Definition: ReadoutGeomCnvAlg.h:120
MuonGM::MMReadoutElement::m_etaDesign
std::array< MuonChannelDesign, 4 > m_etaDesign
Definition: MMReadoutElement.h:157
MuonGMR4::sTgcReadoutElement::globalChannelPosition
Amg::Vector3D globalChannelPosition(const ActsGeometryContext &ctx, const Identifier &measId) const
Returns the global pad/strip/wireGroup position.
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
MmAlignmentStore
Definition: MmAlignmentStore.h:15
MuonGMR4::MdtReadoutElement
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/MdtReadoutElement.h:22
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
MuonGMR4::MuonReadoutElement::alignableTransform
const GeoAlignableTransform * alignableTransform() const
Returnsthe alignable transform of the readout element.
MuonGM::RpcReadoutElement::stripPos
Amg::Vector3D stripPos(const Identifier &id) const
Definition: MuonDetDescr/MuonReadoutGeometry/src/RpcReadoutElement.cxx:177
MuonGM::sTgcReadoutElement::stripGlobalPosition
bool stripGlobalPosition(const Identifier &id, Amg::Vector3D &gpos) const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/sTgcReadoutElement.h:332
sTgcIdHelper
Definition: sTgcIdHelper.h:55
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
Amg::getRotateY3D
Amg::Transform3D getRotateY3D(double angle)
get a rotation transformation around Y-axis
Definition: GeoPrimitivesHelpers.h:261
MuonGMR4::MuonReadoutElement::identify
Identifier identify() const override final
Return the athena identifier.
MuonGMR4::RpcReadoutElement::nPhiStrips
unsigned int nPhiStrips() const
Number of strips measuring the phi coordinate.
sTgcAlignmentStore
Definition: sTgcAlignmentStore.h:15
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
MuonGM::TgcReadoutElement::wirePitch
double wirePitch() const
Returns the pitch of the wires.
MuonGMR4::TgcReadoutElement::numChannels
unsigned int numChannels(const IdentifierHash &measHash) const
Returns the number of readout channels.
MuonGMR4::sTgcReadoutElement::numLayers
unsigned int numLayers() const
Returns the number of gas gap layers.
MuonGMR4::sTgcReadoutElement::padPhi
unsigned int padPhi(const Identifier &measId) const
Returns the Phi index of the pad for the given pad identifier.
MuonGMR4::sTgcReadoutElement
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/sTgcReadoutElement.h:21
MuonGMR4::StripLayer::design
const StripDesign & design() const
Returns the underlying strip design.
SG::ReadCondHandleKey< ActsTrk::DetectorAlignStore >
MuonGM::MuonPadDesign::Length
double Length
Definition: MuonPadDesign.h:55
MuonGM::StripArray
TgcReadoutParams::StripArray StripArray
Definition: RDBReaderAtlas.cxx:23
MuonGM::MuonChannelDesign
Definition: MuonChannelDesign.h:24
MuonGMR4::StripLayer
Definition: StripLayer.h:10
MuonGM::MdtReadoutElement::tubePos
Amg::Vector3D tubePos(const Identifier &id) const
Returns the global position of the given tube.
a
TList * a
Definition: liststreamerinfos.cxx:10
h
MmIdHelper
Definition: MmIdHelper.h:54
MuonGM::MuonReadoutElement::getRsize
double getRsize() const
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:197
MuonGMR4::MmReadoutElement::layerHash
IdentifierHash layerHash(const Identifier &measId) const override final
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
MuonGMR4::sTgcReadoutElement::globalPadCorners
globalCornerArray globalPadCorners(const ActsGeometryContext &ctx, const Identifier &measId) const
AthCommonMsg< Gaudi::Algorithm >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
MuonGM::MuonReadoutElement::identify
Identifier identify() const override final
Returns the ATLAS Identifier of the MuonReadOutElement.
Definition: MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/MuonReadoutElement.h:184
MuonGM::TgcReadoutParams
Definition: TgcReadoutParams.h:38
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
MuonGM::RpcReadoutElement::Nstrips
int Nstrips(bool measphi) const
returns the number of strips for the phi or eta plane
python.DecayParser.children
children
Definition: DecayParser.py:32
MuonGMR4::MuonReadoutElement::localToGlobalTrans
const Amg::Transform3D & localToGlobalTrans(const ActsGeometryContext &ctx) const
Returns the local to global transformation into the ATLAS coordinate system.
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/src/MuonReadoutElement.cxx:81
LArCellBinning.step
step
Definition: LArCellBinning.py:158
MuonGMR4::MdtReadoutElement::globalTubePos
Amg::Vector3D globalTubePos(const ActsGeometryContext &ctx, const Identifier &measId) const
Returns the global position of the tube center.
AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthCommonReentrantAlgorithm.cxx:107
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
MuonGMR4::MdtReadoutElement::complementaryRE
const MdtReadoutElement * complementaryRE() const
Returns the pointer to the other readout element inside the muon station.
MuonGMR4::StripDesign::stripPitch
double stripPitch() const
Distance between two adjacent strips.
CaloLCW_tf.trf
trf
Definition: CaloLCW_tf.py:20
MuonGM::MdtReadoutElement::transform
virtual const Amg::Transform3D & transform(const Identifier &id) const override final
Return local to global transform associated with this identifier.
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
MuonGMR4::PadDesign
Definition: PadDesign.h:24
AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Definition: AthCommonReentrantAlgorithm.h:114
MuonGMR4::RadialStripDesign
Definition: RadialStripDesign.h:23
MuonGMR4::ReadoutGeomCnvAlg::m_writeKey
SG::WriteCondHandleKey< MuonGM::MuonDetectorManager > m_writeKey
Definition: ReadoutGeomCnvAlg.h:122
Amg::getTranslateX3D
Amg::Transform3D getTranslateX3D(const double X)
: Returns a shift transformation along the x-axis
Definition: GeoPrimitivesHelpers.h:277
SG::WriteCondHandle
Definition: WriteCondHandle.h:26
MuonGMR4::MdtReadoutElement::isValid
bool isValid(const IdentifierHash &measHash) const
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
MuonGMR4::TgcReadoutElement
Definition: MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/TgcReadoutElement.h:19
calibdata.tube
tube
Definition: calibdata.py:31
MuonGMR4::ReadoutGeomCnvAlg::m_detMgr
const MuonGMR4::MuonDetectorManager * m_detMgr
Definition: ReadoutGeomCnvAlg.h:127
MuonGM::MuonChannelDesign::ChannelType::etaStrip
@ etaStrip
fitman.k
k
Definition: fitman.py:528
generate::Zero
void Zero(TH1D *hin)
Definition: generate.cxx:32
ServiceHandle< ICondSvc >
MuonGMR4::MuonDetectorManager::getReadoutElement
const MuonReadoutElement * getReadoutElement(const Identifier &id) const
Returns a generic Muon readout element.
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Wed Apr 23 2025 22:01:55 for ATLAS Offline Software by   doxygen 1.8.18