ALinePar.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonAlignmentData/ALinePar.h"
#include "GeoPrimitives/GeoPrimitivesHelpers.h"
#include "GaudiKernel/SystemOfUnits.h"
 
std::ostream& operator<<(std::ostream& ostr, const ALinePar& par) {
  using Parameter = ALinePar::Parameter;
  ostr<<"ALine AMDB id (name,eta,phi,job)=(";
  ostr<<par.AmdbStation()<<",";
  ostr<<par.AmdbEta()<<",";
  ostr<<par.AmdbPhi()<<",";
  ostr<<par.AmdbJob()<<"), ";
  ostr<<"translation (S/Y,Z,T/X)= (";
  ostr<<par.getParameter(Parameter::transS)<<",";
  ostr<<par.getParameter(Parameter::transZ)<<",";
  ostr<<par.getParameter(Parameter::transT)<<"), ";
  ostr<<"rotation angle";
  ostr<<" around S="<<par.getParameter(Parameter::rotS)<<",";
  ostr<<" around Z="<<par.getParameter(Parameter::rotZ)<<",";
  ostr<<" around T="<<par.getParameter(Parameter::rotT)<<" ";
  return ostr;
}
void ALinePar::setParameters(float s, float z, float t, float rotS, float rotZ, float rotT) {
    m_payload[static_cast<unsigned int>(Parameter::transS)] = s;
    m_payload[static_cast<unsigned int>(Parameter::transZ)] = z;
    m_payload[static_cast<unsigned int>(Parameter::transT)] = t;
    m_payload[static_cast<unsigned int>(Parameter::rotS)  ] = rotS;
    m_payload[static_cast<unsigned int>(Parameter::rotZ)  ] = rotZ;
    m_payload[static_cast<unsigned int>(Parameter::rotT)  ] = rotT;
}
 
HepGeom::Transform3D ALinePar::deltaTransform() const {  // does NOT account for AMDB origin being different from volume centre;
    // for that you would need access to full station Position info...
    // see MuonGeoModel/Station::getDeltaTransform() for details.
    return HepGeom::TranslateY3D(getParameter(Parameter::transS)) * 
           HepGeom::TranslateZ3D(getParameter(Parameter::transZ)) * 
           HepGeom::TranslateX3D(getParameter(Parameter::transT)) * 
           HepGeom::RotateY3D(getParameter(Parameter::rotS)) *
           HepGeom::RotateZ3D(getParameter(Parameter::rotZ)) * 
           HepGeom::RotateX3D(getParameter(Parameter::rotT));
}
Amg::Transform3D  ALinePar::delta () const{
    return Amg::Translation3D(getParameter(Parameter::transS)*Amg::Vector3D::UnitY()) * 
           Amg::Translation3D(getParameter(Parameter::transZ)*Amg::Vector3D::UnitZ()) * 
           Amg::Translation3D(getParameter(Parameter::transT)*Amg::Vector3D::UnitX()) * 
           Amg::getRotateY3D(getParameter(Parameter::rotS)) *
           Amg::getRotateZ3D(getParameter(Parameter::rotZ)) * 
           Amg::getRotateX3D(getParameter(Parameter::rotT));
}