#include <AlignmentErrorTool.h>

Inheritance diagram for MuonAlign::AlignmentErrorTool:

Collaboration diagram for MuonAlign::AlignmentErrorTool:

Classes
struct	deviationSummary_t

Public Member Functions
	AlignmentErrorTool (const std::string &, const std::string &, const IInterface *)

	~AlignmentErrorTool () override=default

StatusCode	initialize () override

void	makeAlignmentDeviations (const Trk::Track &track, std::vector< Trk::AlignmentDeviation * > &deviations) const override
	Compute alignment deviations, given a track as input. More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

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
std::string	hardwareName (MuonCalib::MuonFixedLongId calibId) const

std::string	sectorString (MuonCalib::MuonFixedLongId calibId) const

int	sector (MuonCalib::MuonFixedLongId calibId) const

int	hardwareEta (MuonCalib::MuonFixedLongId calibId) 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...

Static Private Member Functions
static std::string_view	side (MuonCalib::MuonFixedLongId calibId)

static bool	isSmallSector (MuonCalib::MuonFixedLongId calibId)

Private Attributes
ToolHandle< MuonCalib::IIdToFixedIdTool >	m_idTool {this, "idTool", "MuonCalib::IdToFixedIdTool"}

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

SG::ReadCondHandleKey< MuonAlignmentErrorData >	m_readKey {this, "ReadKey", "MuonAlignmentErrorData", "Key of MuonAlignmentErrorData"}

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 29 of file AlignmentErrorTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ AlignmentErrorTool()

MuonAlign::AlignmentErrorTool::AlignmentErrorTool	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 27 of file AlignmentErrorTool.cxx.

                                                                                                   : AthAlgTool(t, n, p) {
     declareInterface<Trk::ITrkAlignmentDeviationTool>(this);
 }

◆ ~AlignmentErrorTool()

MuonAlign::AlignmentErrorTool::~AlignmentErrorTool ( )

overridedefault

Member Function Documentation

◆ declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

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

inlineprivateinherited

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

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation 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.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation 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.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

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

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation 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.

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

◆ declareProperty() [5/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation 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.

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

◆ declareProperty() [6/6]

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

◆ hardwareEta()

int MuonAlign::AlignmentErrorTool::hardwareEta ( MuonCalib::MuonFixedLongId calibId ) const

inlineprivate

Definition at line 361 of file AlignmentErrorTool.cxx.

                                                                                  {
     using StationName = MuonCalib::MuonFixedLongId::StationName;
     switch (calibId.stationName()) {
         case StationName::BML:
             {
                 if (sector(calibId)==13) {
                     switch (calibId.eta()) {
                         case 4: return 5;
                         case 5: return 6;
                         case 6: return 7;
                         case -4: return -5;
                         case -5: return -6;
                         case -6: return -7;
                     }
                 }
                 return calibId.eta();
             }
         case StationName::BOL:
             {
                 if (sector(calibId)==13) {
                     if (calibId.eta()== 7) return 1; // BOE1A13 not BOL7A13
                     if (calibId.eta()==-7) return -1; // BOE1C13 not BOL7C13
                 }
                 return calibId.eta();
             }
         case StationName::BOF:
             return calibId.eta()>0 ? calibId.eta()*2-1 : calibId.eta()*2+1;
         case StationName::BOG:
             return calibId.eta()*2;
         case StationName::EIL:
             {
                 if ((sector(calibId) == 1) || (sector(calibId) == 9)) {
                     switch (calibId.eta()) {
                         case 4: return 5;
                         case 5: return 4;
                         case -4: return -5;
                         case -5: return -4;
                     }
                 }
                 return calibId.eta();
             }
         case StationName::EEL:
             {
                 if ((sector(calibId) == 5) && (calibId.eta() == 1)) return 2;
                 if ((sector(calibId) == 5) && (calibId.eta() == -1)) return -2;
                 return calibId.eta();
             }
         default: return calibId.eta();
     }
 }

◆ hardwareName()

std::string MuonAlign::AlignmentErrorTool::hardwareName ( MuonCalib::MuonFixedLongId calibId ) const

inlineprivate

Definition at line 295 of file AlignmentErrorTool.cxx.

                                                                                         {
     using StationName = MuonCalib::MuonFixedLongId::StationName;
  
     // The only exception that cannot be caught by hardwareEta() above
     if (sector(calibId)==13) {
         if (calibId.eta()== 7 && calibId.stationName()==StationName::BOL) return "BOE1A13"; // BOE1A13 not BOL7A13
         if (calibId.eta()==-7 && calibId.stationName()==StationName::BOL) return "BOE1C13"; // BOE1C13 not BOL7C13
         if (calibId.eta()== 8 && calibId.stationName()==StationName::BOL) return "BOE2A13"; // BOE2A13 not BOL8A13
         if (calibId.eta()==-8 && calibId.stationName()==StationName::BOL) return "BOE2C13"; // BOE2C13 not BOL8C13
     }
  
     std::string ret { calibId.stationNameString() };
     ret.push_back(static_cast<char>('0'+std::abs(hardwareEta(calibId))));
     ret.append(side(calibId)).append(sectorString(calibId));
  
     return ret;
 }

◆ initialize()

StatusCode MuonAlign::AlignmentErrorTool::initialize ( )

override

Definition at line 31 of file AlignmentErrorTool.cxx.

                                           {
     ATH_MSG_INFO("*****************************************");
     ATH_MSG_INFO("AlignmentErrorTool::initialize()");
  
     ATH_CHECK(m_idTool.retrieve());
     ATH_CHECK(m_idHelperSvc.retrieve());
     ATH_CHECK(m_readKey.initialize());
  
     return StatusCode::SUCCESS;
 }

◆ inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

◆ interfaceID()

const InterfaceID & Trk::ITrkAlignmentDeviationTool::interfaceID ( )

inlinestaticinherited

Definition at line 37 of file ITrkAlignmentDeviationTool.h.

 {
   return IID_ITrkAlignmentDeviationTool;
 }

◆ isSmallSector()

bool MuonAlign::AlignmentErrorTool::isSmallSector ( MuonCalib::MuonFixedLongId calibId )

inlinestaticprivate

Definition at line 337 of file AlignmentErrorTool.cxx.

                                                                               {
     using StationName = MuonCalib::MuonFixedLongId::StationName;
     switch (calibId.stationName()) {
         case StationName::BIS:
         case StationName::BMS:
         case StationName::BOS:
         case StationName::BEE:
         case StationName::BMF:
         case StationName::BOF:
         case StationName::BOG:
         case StationName::EES:
         case StationName::EMS:
         case StationName::EOS:
         case StationName::EIS:
         case StationName::CSS:
         case StationName::BMG:
         case StationName::MMS:
         case StationName::STS:
             return true;
         default:
             return false;
     }
 }

◆ makeAlignmentDeviations()

void MuonAlign::AlignmentErrorTool::makeAlignmentDeviations	(	const Trk::Track &	track,
		std::vector< Trk::AlignmentDeviation * > &	deviations
	)		const

overridevirtual

Compute alignment deviations, given a track as input.

Implements Trk::ITrkAlignmentDeviationTool.

Definition at line 42 of file AlignmentErrorTool.cxx.

                                                                                                                            {
     ATH_MSG_DEBUG("AlignmentErrorTool::makeAlignmentDeviations()");
  
     SG::ReadCondHandle<MuonAlignmentErrorData> readHandle{m_readKey};
     const MuonAlignmentErrorData* readCdo{*readHandle};
     if (!readCdo) {
         ATH_MSG_ERROR("nullptr to the read conditions object");
         return;
     }
     const auto& deviationVec = readCdo->getDeviations();
     std::vector<deviationSummary_t> devSumVec;
     devSumVec.reserve(deviationVec.size());
     deviationSummary_t aDevSumm;
     for (const auto & i : deviationVec) {
         aDevSumm.translation = i.translation;
         aDevSumm.rotation = i.rotation;
         aDevSumm.stationName = i.stationName;
         aDevSumm.multilayer = i.multilayer;
         devSumVec.emplace_back(std::move(aDevSumm));
     }
  
     typedef Trk::TrackStates tsosc_t;
     const tsosc_t* tsosc = track.trackStateOnSurfaces();
  
     // LOOP ON HITS ON TRACK //
     unsigned int nPrecisionHits = 0;
     for (const auto *tsos : *tsosc) {
         if (!tsos->type(Trk::TrackStateOnSurface::Measurement)) {
             continue;
         }
         const Trk::MeasurementBase* meas = tsos->measurementOnTrack();
         const auto* rot = dynamic_cast<const Trk::RIO_OnTrack*>(meas);
  
         if (!rot) {
             const auto* crot = dynamic_cast<const Trk::CompetingRIOsOnTrack*>(meas);
             if (crot) {
                 unsigned int index = crot->indexOfMaxAssignProb();
                 rot = &(crot->rioOnTrack(index));
             }
         }
         if (!rot) continue;
  
         Identifier channelId = rot->identify();
         if (!m_idHelperSvc->isMuon(channelId)) {
             // the RIO_OnTrack Identifiers could also come from ID or Calo, but this tool is only interested in MS hits
             ATH_MSG_VERBOSE("Given Identifier " << channelId.get_compact() << " is not a muon Identifier, continuing");
             continue;
         }
  
         // Keep only the precision coordinate hits
         if (m_idHelperSvc->isRpc(channelId)
             || m_idHelperSvc->isTgc(channelId)
             || (m_idHelperSvc->isCsc(channelId) && m_idHelperSvc->cscIdHelper().measuresPhi(channelId) == 1)
             || (m_idHelperSvc->issTgc(channelId) && m_idHelperSvc->stgcIdHelper().channelType(channelId) != sTgcIdHelper::sTgcChannelTypes::Strip)) {
             continue;
         }
  
         // To maintain backward compatibility with the old error CLOBs, activate
         // the NSW hits only if specified in the CLOB, else disccard them
         if (!readCdo->hasNswHits()
             && (m_idHelperSvc->issTgc(channelId) || m_idHelperSvc->isMM(channelId))) {
             continue;
         }
  
         MuonCalib::MuonFixedLongId calibId = m_idTool->idToFixedLongId(channelId);
         if (!calibId.isValid()) continue;
  
         // GATHERING INFORMATION TO PUT TOGETHER THE STATION NAME //
         std::string alignStationName = hardwareName(calibId);
         int multilayer = 1;
         if (calibId.is_mdt()) {
             multilayer = calibId.mdtMultilayer();
         } else if (calibId.is_mmg()) {
             multilayer = calibId.mmgMultilayer();
         } else if (calibId.is_stg()) {
             multilayer = calibId.stgMultilayer();
         }
         std::string multilayerName = std::to_string(multilayer);
  
         ATH_MSG_DEBUG("Hit is in station " << alignStationName << " multilayer " << multilayerName);
         ++nPrecisionHits;
  
         // Compute deviationSummary_t building blocks
         const Trk::PrepRawData* prd = rot->prepRawData();
         const Trk::Surface& sur = prd->detectorElement()->surface(prd->identify());
  
         double w2 = 1.0 / (rot->localCovariance()(Trk::loc1, Trk::loc1));
         Amg::Vector3D hitP = tsos->trackParameters()->position();
         Amg::Vector3D hitU = tsos->trackParameters()->momentum().unit();
  
         // Wire direction for MDT, strip direction for MM or sTGC
         int icol = (calibId.is_mdt()||calibId.is_csc()) ? 2 : 1;
         Amg::Vector3D hitV = sur.transform().rotation().col(icol);
  
         // Enforce orientation of the V vectors
         static const Amg::Vector3D zATLAS(0., 0., 1.);
         if (hitP.cross(zATLAS).dot(hitV)<0.0) {
             hitV *= -1.0;
         }
  
         // FOR CROSS-CHECK
         bool is_matched = false;
  
         // LOOP ON STATION DEVIATIONS EXTRACTED FROM INPUT FILE //
         for (auto & iDev : devSumVec) {
             // try to regexp-match the station name and the multilayer name
             if (!boost::regex_match(alignStationName, iDev.stationName)) {
                 continue;
             }
             if (!boost::regex_match(multilayerName, iDev.multilayer)) {
                 continue;
             }
  
             // ASSOCIATE EACH NUISANCE TO A LIST OF HITS
             iDev.hits.push_back(rot);
  
             iDev.sumW2 += w2;
             iDev.sumP += w2 * hitP;
             iDev.sumU += w2 * hitU;
             iDev.sumV += w2 * hitV;
  
             // FOR CROSS-CHECK
             is_matched = true;
  
  
         }  // LOOP ON DEVIATIONS
  
         if (!is_matched) {
             ATH_MSG_WARNING("The hits in the station " << alignStationName << ", multilayer " << multilayerName
                             << " couldn't be matched to any deviation regexp in the list.");
         }
  
     }  // LOOP ON TSOS
  
     // Nuisance parameters covering the complete track are not wanted. (MS/ID
     // error treated differently for now). Removing the deviations covering the
     // full track in further processing.
     for (auto& dev : devSumVec) {
         if (dev.hits.size() == nPrecisionHits) {
             dev.hits.clear();
         }
     }
  
     // GET RID OF PERFECT OVERLAPS BY COMBINING ERRORS
     std::vector<const Trk::RIO_OnTrack*> v1, v2;
     for (auto iti = devSumVec.begin(); iti != devSumVec.end(); ++iti) {
  
         v1 = iti->hits;
         if (v1.empty()) {
             continue;
         }
  
         std::stable_sort(v1.begin(), v1.end());
  
         for (auto itj = iti+1; itj != devSumVec.end(); ++itj) {
  
             if (iti->hits.size() != itj->hits.size()) {
                 continue;
             }
  
             v2 = itj->hits;
             std::stable_sort(v2.begin(), v2.end());
  
             if (v1 == v2) {
                 auto iDev = std::distance(devSumVec.begin(), iti);
                 auto jDev = std::distance(devSumVec.begin(), itj);
                 ATH_MSG_DEBUG("Found deviations " << iDev << " and " << jDev << " related to the same list of hits. Merging...");
                 ATH_MSG_DEBUG("old (translation, rotation) systematic uncertainties for "
                               << iDev << ": " << iti->translation << ", " << iti->rotation);
                 ATH_MSG_DEBUG("old (translation, rotation) systematic uncertainties for "
                               << jDev << ": " << itj->translation << ", " << itj->rotation);
  
                 // MERGE THE TWO DEVIATIONS ASSOCIATED TO THE SAME LIST OF HITS //
                 double new_translation = std::hypot(iti->translation, itj->translation);
                 double new_rotation = std::hypot(iti->rotation, itj->rotation);
  
                 // NOW PREPARE TO ERASE ONE OF THE TWO COPIES //
                 itj->hits.clear();
  
                 // ASSIGN NEW TRASLATION/ROTATION TO THE REMAINING COPY //
                 iti->translation = new_translation;
                 iti->rotation = new_rotation;
                 ATH_MSG_DEBUG("New combined (translation, rotation) systematic uncertainties: " << new_translation << ", " << new_rotation);
  
             }  // FIND AN OVERLAP IN THE HITS LISTS
         }      // SECOND LOOP ON DEVIATIONS
     }          // FIRST LOOP ON DEVIATIONS
  
  
     // NOW BUILD THE DEVIATIONS
     deviations.clear();
     ATH_MSG_DEBUG("************************************");
     ATH_MSG_DEBUG("FINAL LIST OF DEVIATIONS");
     for (const auto & iDev : devSumVec) {
         if (iDev.hits.empty()) {
             continue;
         }
  
         double rotation = iDev.rotation;
         double translation = iDev.translation;
  
         Amg::Vector3D sumP = iDev.sumP;
         Amg::Vector3D sumU = iDev.sumU;
         Amg::Vector3D sumV = iDev.sumV;
         double sumW2 = iDev.sumW2;
  
         sumP *= (1. / sumW2);
         sumU *= (1. / sumW2);
         sumV *= (1. / sumW2);
  
         if (translation >= 0.001 * Gaudi::Units::mm) {
             std::size_t hitshash = 0;
             for (const auto *it : iDev.hits) {
                 boost::hash_combine(hitshash, (it->identify()).get_compact());
             }
             deviations.push_back(
                 new AlignmentTranslationDeviation(sumU.cross(sumV), translation * Gaudi::Units::mm, iDev.hits));
             deviations.back()->setHashOfHits(hitshash);
  
             ATH_MSG_DEBUG("A translation along ("
                           << sumU.x() << ", " << sumU.y() << ", " << sumU.z() << ") with sigma=" << translation * Gaudi::Units::mm
                           << " mm was applied to " << iDev.hits.size()
                           << " hits matching the station: " << iDev.stationName.str() << " and the multilayer "
                           << iDev.multilayer.str());
         }
         if (rotation >= 0.000001 * Gaudi::Units::rad) {
             std::size_t hitshash = 0;
             for (const auto *it : iDev.hits) {
                 boost::hash_combine(hitshash, (it->identify()).get_compact());
             }
             deviations.push_back(new AlignmentRotationDeviation(sumP, sumV, rotation * Gaudi::Units::rad, iDev.hits));
             deviations.back()->setHashOfHits(hitshash);
  
             ATH_MSG_DEBUG("A rotation around the center = (" << sumP.x() << ", " << sumP.y() << ", " << sumP.z() << ") and axis = ("
                                                              << sumV.x() << ", " << sumV.y() << ", " << sumV.z()
                                                              << ") with sigma=" << rotation / Gaudi::Units::mrad << " mrad was applied to "
                                                              << iDev.hits.size() << " hits matching the station "
                                                              << iDev.stationName.str() << " and the multilayer "
                                                              << iDev.multilayer.str());
         }
  
     }  // LOOP ON NUISANCES
  
     ATH_MSG_DEBUG("******************************");
     ATH_MSG_DEBUG("FINAL CHECKUP");
     ATH_MSG_DEBUG("Found " << deviations.size() << " nuisances after duplicates merging");
     ATH_MSG_DEBUG("******************************");
 }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

◆ renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sector()

int MuonAlign::AlignmentErrorTool::sector ( MuonCalib::MuonFixedLongId calibId ) const

inlineprivate

Definition at line 328 of file AlignmentErrorTool.cxx.

                                                                             {
     if (calibId.is_tgc()) {
         // TGC sector convention is special
         return calibId.phi();
     } else {
         return isSmallSector(calibId) ? 2*calibId.phi() : 2*calibId.phi()-1;
     }
 }

◆ sectorString()

std::string MuonAlign::AlignmentErrorTool::sectorString ( MuonCalib::MuonFixedLongId calibId ) const

inlineprivate

Definition at line 317 of file AlignmentErrorTool.cxx.

                                                                                         {
     int sec = sector(calibId);
     if (sec<0 || sec > 99) {
         throw std::runtime_error("Unhandled sector number");
     }
     std::string ret = "00";
     ret[0] += (sec/10);
     ret[1] += (sec%10);
     return ret;
 }

◆ side()

std::string_view MuonAlign::AlignmentErrorTool::side ( MuonCalib::MuonFixedLongId calibId )

inlinestaticprivate

Definition at line 313 of file AlignmentErrorTool.cxx.

                                                                                {
     return calibId.eta()>0 ? "A" : calibId.eta()<0 ? "C" : "B";
 }

◆ sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

overridevirtualinherited

Handle START transition.

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

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonAlign::AlignmentErrorTool::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 43 of file AlignmentErrorTool.h.

◆ m_idTool

ToolHandle<MuonCalib::IIdToFixedIdTool> MuonAlign::AlignmentErrorTool::m_idTool {this, "idTool", "MuonCalib::IdToFixedIdTool"}

private

Definition at line 42 of file AlignmentErrorTool.h.

◆ m_readKey

SG::ReadCondHandleKey<MuonAlignmentErrorData> MuonAlign::AlignmentErrorTool::m_readKey {this, "ReadKey", "MuonAlignmentErrorData", "Key of MuonAlignmentErrorData"}

private

Definition at line 72 of file AlignmentErrorTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

The documentation for this class was generated from the following files:

Classes

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ AlignmentErrorTool()

◆ ~AlignmentErrorTool()

Member Function Documentation

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ hardwareEta()

◆ hardwareName()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ isSmallSector()

◆ makeAlignmentDeviations()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sector()

◆ sectorString()

◆ side()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_idHelperSvc

◆ m_idTool

◆ m_readKey

◆ m_varHandleArraysDeclared

◆ m_vhka