MuonSegment.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonSegment/MuonSegment.h"
// TrkEvent
#include "CxxUtils/sincos.h"
#include "EventPrimitives/EventPrimitivesToStringConverter.h"
#include "GeoPrimitives/GeoPrimitivesToStringConverter.h"
#include "TrkEventPrimitives/DefinedParameter.h"
#include "TrkEventPrimitives/FitQuality.h"
#include "TrkEventPrimitives/LocalParameters.h"
#include "TrkEventPrimitives/ParamDefs.h"
#include "TrkRIO_OnTrack/RIO_OnTrack.h"
 
#include <cfloat>
#include <iostream>
namespace Muon {
 
const float MuonSegment::kNoValue = FLT_MAX;
 
MuonSegment::MuonSegment()
  : Segment()
  , SpaceTimePointBase(kNoValue, kNoValue, 1.)
  , Trk::SurfacePtrHolderImplDetEl<Trk::PlaneSurface>(nullptr)
  , m_globalPosition()
  , m_globalDirection()
  , m_localDirection()
{}
 
MuonSegment::MuonSegment(const MuonSegment& seg)
  : Segment(seg)
  , SpaceTimePointBase(seg)
  , Trk::SurfacePtrHolderImplDetEl<Trk::PlaneSurface>(nullptr)
  , m_globalPosition(seg.m_globalPosition)
  , m_globalDirection(seg.m_globalDirection)
  , m_localDirection(seg.m_localDirection)
{
  m_associatedSurface = (seg.m_associatedSurface->associatedDetectorElement())
                          ? seg.m_associatedSurface
                          : seg.m_associatedSurface->clone();
}
 
MuonSegment&
MuonSegment::operator=(const MuonSegment& seg)
{
  if (this != &seg) {
    Trk::Segment::operator=(seg);
    Trk::SpaceTimePointBase::operator=(seg);
    Trk::SurfacePtrHolderImplDetEl<Trk::PlaneSurface>::operator=(seg);
    m_globalPosition = seg.m_globalPosition;
    m_globalDirection = seg.m_globalDirection;
    m_localDirection = seg.m_localDirection;
  }
  return (*this);
}
 
MuonSegment::MuonSegment(Trk::LocalParameters&& locpars,
                         Amg::MatrixX&& locerr,
                         Trk::PlaneSurface* psf,
                         DataVector<const Trk::MeasurementBase>&& cmeas,
                         Trk::FitQuality* fqual,
                         Segment::Author author)
  : Segment(std::move(locpars), std::move(locerr),
            std::move(cmeas), fqual, author)
  , SpaceTimePointBase(kNoValue, kNoValue, 1.)
  , Trk::SurfacePtrHolderImplDetEl<Trk::PlaneSurface>(psf)
  , m_globalPosition()
  , m_globalDirection()
  , m_localDirection()
{
  recalculateCache();
  psf->globalToLocalDirection(m_globalDirection, m_localDirection);
}
 
MuonSegment::MuonSegment(const Amg::Vector2D& locSegPos,
                         const Trk::LocalDirection& locSegDir,
                         Amg::MatrixX&& locErr,
                         Trk::PlaneSurface* psf,
                         DataVector<const Trk::MeasurementBase>&& cmeas,
                         Trk::FitQuality* fqual,
                         Segment::Author author)
  : Segment(Trk::LocalParameters(), std::move(locErr),
            std::move(cmeas), fqual, author)
  , SpaceTimePointBase(kNoValue, kNoValue, 1.)
  , Trk::SurfacePtrHolderImplDetEl<Trk::PlaneSurface>(psf)
  , m_globalPosition()
  , m_globalDirection()
  , m_localDirection(locSegDir)
{
  psf->localToGlobalDirection(locSegDir, m_globalDirection);
 
  Amg::Vector2D lpos(locSegPos[Trk::locX], locSegPos[Trk::locY]);
  m_associatedSurface->localToGlobal(lpos, m_globalDirection, m_globalPosition);
 
  double phi = m_globalDirection.phi();
  double theta = m_globalDirection.theta();
  std::array<Trk::DefinedParameter, 4> pars = {
      {{locSegPos[Trk::locX], Trk::locX},
       {locSegPos[Trk::locY], Trk::locY},
       {phi, Trk::phi},
       {theta, Trk::theta}}};
  m_localParams = Trk::LocalParameters(pars);
}
 
MuonSegment::~MuonSegment() = default;
 
void
MuonSegment::recalculateCache()
{
 
  CxxUtils::sincos scphi(m_localParams[Trk::phi]);
  CxxUtils::sincos sctheta(m_localParams[Trk::theta]);
 
  double sintheta = sctheta.sn;
  double costheta = sctheta.cs;
  double sinphi = scphi.sn;
  double cosphi = scphi.cs;
  m_globalDirection =
    Amg::Vector3D(cosphi * sintheta, sinphi * sintheta, costheta);
  Amg::Vector2D lpos(m_localParams[Trk::locX], m_localParams[Trk::locY]);
  m_associatedSurface->localToGlobal(lpos, m_globalPosition, m_globalPosition);
}
 
MsgStream&
MuonSegment::dump(MsgStream& out) const
{
  out << "Muon::MuonSegment " << std::endl;
  out << "  - global position      : " << std::endl
      << Amg::toString(globalPosition()) << std::endl;
  out << std::endl << localDirection() << std::endl;
  out << "  - local parameters     : " << std::endl
      << localParameters() << std::endl;
  out << "  - local Err Matrix     : ";
  out << std::endl << Amg::toString(localCovariance()) << std::endl;
  if (nullptr != m_fitQuality)
    out << "  - fit Quality          : " << std::endl
        << (*fitQuality()) << std::endl;
  else
    out << "  - fit Quality          : NO " << std::endl;
  out << "  - associated Surface   : " << std::endl
      << associatedSurface() << std::endl;
  out << "  - t0 (error) : " << time() << "(" << errorTime() << ")"
      << std::endl;
  out << "  - it contains   : " << numberOfContainedROTs()
      << " RIO_OnTrack object" << std::endl;
  unsigned int numRoT = 1;
  for (const Trk::MeasurementBase* m : containedMeasurements()) {
    const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(m);
    if (rot) {
      out << "RoT " << numRoT++ << std::endl << *rot << std::endl;
    }
  }
  return out;
}
 
std::ostream&
MuonSegment::dump(std::ostream& out) const
{
  out << "Muon::MuonSegment " << std::endl;
  out << "  - global position      : " << std::endl
      << Amg::toString(globalPosition()) << std::endl;
  out << std::endl << localDirection() << std::endl;
  out << "  - local parameters     : " << std::endl
      << localParameters() << std::endl;
  out << "  - local Err Matrix     : ";
  out << std::endl << Amg::toString(localCovariance()) << std::endl;
  if (nullptr != m_fitQuality)
    out << "  - fit Quality          : " << std::endl
        << (*fitQuality()) << std::endl;
  else
    out << "  - fit Quality          : NO " << std::endl;
  out << "  - associated Surface   : " << std::endl
      << associatedSurface() << std::endl;
  out << "  - t0 (error) : " << time() << "(" << errorTime() << ")"
      << std::endl;
  out << "  - it contains   : " << numberOfContainedROTs()
      << " RIO_OnTrack object" << std::endl;
  unsigned int numRoT = 1;
  for (const Trk::MeasurementBase* m : containedMeasurements()) {
    const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(m);
    if (rot) {
      out << "RoT " << numRoT++ << std::endl << *rot << std::endl;
    }
  }
  return out;
}
 
}