SegmentSeed.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "MuonPatternEvent/SegmentSeed.h"
#include "CxxUtils/sincos.h"
namespace MuonR4{
    using namespace SegmentFit;
 
     SegmentSeed::SegmentSeed(double tanTheta, double interceptY, double tanPhi,
                              double interceptX, double counts,
                              std::vector<HitType>&& hits,
                              const SpacePointBucket* bucket):
        m_parent{bucket},
        m_hits{std::move(hits)},
        m_hasPhiExt{true},
        m_counts{counts}{
            Amg::Vector3D dir{tanPhi, tanTheta, 1.};
            m_pars[toInt(ParamDefs::x0)] = interceptX;
            m_pars[toInt(ParamDefs::y0)] = interceptY;
            m_pars[toInt(ParamDefs::theta)] = dir.theta();
            m_pars[toInt(ParamDefs::phi)] = dir.phi();
    }
    SegmentSeed::SegmentSeed(const HoughMaximum& toCopy):
        m_parent{toCopy.parentBucket()},
        m_hits{toCopy.getHitsInMax()},
        m_counts{toCopy.getCounts()}{
            Amg::Vector3D dir{0., toCopy.tanTheta(), 1.};
            m_pars[toInt(ParamDefs::y0)]    = toCopy.interceptY();
            m_pars[toInt(ParamDefs::theta)] = dir.theta();
            m_pars[toInt(ParamDefs::phi)] = dir.phi();
    }
    double SegmentSeed::tanPhi() const{ return std::tan(m_pars[toInt(ParamDefs::phi)]); }
    double SegmentSeed::tanTheta() const { return std::tan(m_pars[toInt(ParamDefs::theta)]); }
 
    double SegmentSeed::interceptX() const { return m_pars[toInt(ParamDefs::x0)]; }
    double SegmentSeed::interceptY() const {return m_pars[toInt(ParamDefs::y0)]; }
    const Parameters& SegmentSeed::parameters() const{ return m_pars;}
    double SegmentSeed::getCounts() const{ return m_counts;}
    const std::vector<SegmentSeed::HitType>& SegmentSeed::getHitsInMax() const { return m_hits; }
    const SpacePointBucket* SegmentSeed::parentBucket() const{ return m_parent; }
    const MuonGMR4::MuonChamber* SegmentSeed::chamber() const{ return m_parent->chamber(); }
    bool SegmentSeed::hasPhiExtension() const{ return m_hasPhiExt; }
    Amg::Vector3D SegmentSeed::positionInChamber() const{ return Amg::Vector3D{interceptX(), interceptY(),0.}; }
    Amg::Vector3D SegmentSeed::directionInChamber() const{ 
        const CxxUtils::sincos theta{m_pars[toInt(ParamDefs::theta)]},
                               phi{m_pars[toInt(ParamDefs::phi)]};
        return Amg::Vector3D(phi.cs*theta.sn,phi.sn*theta.sn,theta.cs);
    }
}