sTGCAsBuiltData.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#include "MuonAlignmentData/sTGCAsBuiltData.h"
#include "GeoPrimitives/GeoPrimitivesToStringConverter.h"
 
 
std::ostream& operator<<(std::ostream& ostr, const sTGCAsBuiltData::Parameters& par){
   return ostr << "Offset " << par.offset << " rotxy " << par.rotation << " scale " << par.scale  << " nonPara " <<  par.nonPara << " stripBending " << par.stripBending;
}
 
 
sTGCAsBuiltData::sTGCAsBuiltData(const Muon::IMuonIdHelperSvc* idHelperSvc):
    AthMessaging{"sTGCAsBuiltData"},
    m_idHelperSvc{idHelperSvc} {}
 
 
 
Amg::Vector2D sTGCAsBuiltData::correctPosition(const Identifier& channelId, const Amg::Vector2D& pos) const {
    bool issTgc = m_idHelperSvc->issTgc(channelId);
    Identifier asBuiltId = issTgc ? m_idHelperSvc->gasGapId(channelId) : m_idHelperSvc->mmIdHelper().pcbID(channelId);
 
    ParMap::const_iterator par_itr = m_asBuiltData.find(asBuiltId);
    if(par_itr == m_asBuiltData.end()){
        ATH_MSG_WARNING("Missing as built parameters for gas gap " << m_idHelperSvc->toString(channelId));
        return pos;
    }
    Amg::Vector2D correctedPos = pos;
    Parameters pars = par_itr->second;
    constexpr double convScale = 1.e-3; // parameters are stored in um and mrad therefore dividing by 1000
    
    // For the QL3 modules described by the legacy (Run 3) geometry EDM the origin of the local frame is shifted with respect to the actual center of the gap.
    // Therefore a shift is introduced which stabilizes the as built fit which must also be applied here since the parameters are expressed using the actual center of the gap.
    float shift = (std::fabs(m_idHelperSvc->stationEta(channelId)) == 3 && m_idHelperSvc->stationNameString(channelId) == "STL" ?  24.74 : 0.0);
    ATH_MSG_VERBOSE("applying as built parameters for gas gap " << m_idHelperSvc->toString(m_idHelperSvc->gasGapId(channelId)) << " parameters " << pars);
 
    double correction  = pars.offset * convScale  + pars.rotation *convScale * pos.y() + pars.scale*convScale*(pos.x() + shift ) + pars.nonPara * convScale * convScale * (pos.x() + shift  )*pos.y();
    // we only want to correct the banana shape for the MM
    if(!issTgc) correction += pars.stripBending*convScale*convScale*pos.y()*pos.y();
    correctedPos.x() = pos.x() + correction; 
    return correctedPos;
}
 
 
StatusCode sTGCAsBuiltData::setParameters(const Identifier& channelId, const Parameters& pars){
    const Identifier gasGapId = m_idHelperSvc->gasGapId(channelId);
    auto insert_itr = m_asBuiltData.insert(std::make_pair(gasGapId, pars));
    if (!insert_itr.second) {
        ATH_MSG_ERROR("As built parameters for gasGap "<<m_idHelperSvc->toStringGasGap(gasGapId)
            <<" has already been stored ");
        return StatusCode::FAILURE;
    }
    return StatusCode::SUCCESS;
}