#include <CalibSegmentChi2Minimizer.h>

Inheritance diagram for MuonR4::CalibSegmentChi2Minimizer:

Collaboration diagram for MuonR4::CalibSegmentChi2Minimizer:

Public Types
using	HitType = std::unique_ptr< CalibratedSpacePoint >

using	HitVec = std::vector< HitType >

Public Member Functions
	CalibSegmentChi2Minimizer (const std::string &name, const EventContext &ctx, const Amg::Transform3D &locToGlobTrf, HitVec &&hits, const ISpacePointCalibrator *calibrator, bool doT0Fit)
	Constructor taking the input hits to fit. More...

double	DoEval (const double *pars) const override final
	Evaluate the chi2 for the given set of parameters. More...

CalibSegmentChi2Minimizer *	Clone () const override final
	Clone of the object. More...

unsigned int	NDim () const override final
	Degrees of freedom. More...

bool	hasPhiMeas () const
	Are phi measurements in the collection. More...

bool	doTimeFit () const
	Does the fit include the time. More...

int	nDoF () const
	Returns the degrees of freedom from the measurements & beamspot constraint. More...

const HitVec &	measurements () const
	Returns the used measurements. More...

std::vector< double >	chi2Terms (const double *pars) const
	Returns the contribution of each measurement to the chi2. More...

HitVec	release (const double *pars)
	Releases all measurements. More...

const Amg::Transform3D &	localToGlobTrans () const
	Returns the local to global transformation. More...

Private Member Functions
double	segmentChi2TermStrip (const CalibratedSpacePoint &hit, const Amg::Vector3D &chamberIsect, const Amg::Vector3D &segmentDir) const

double	segmentChi2TermMdt (const CalibratedSpacePoint &hit, const Amg::Vector3D &chamberIsect, const Amg::Vector3D &segmentDir) const

bool	msgLvl (const MSG::Level lvl) const
	Test the output level. More...

MsgStream &	msg () const
	The standard message stream. More...

MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...

void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
std::string	m_name {}

const EventContext &	m_ctx

Amg::Transform3D	m_locToGlob {Amg::Transform3D::Identity()}

HitVec m_hits	ATLAS_THREAD_SAFE {}

const ISpacePointCalibrator *	m_calibrator {}

bool	m_doT0Fit {false}

bool	m_hasPhi {false}

unsigned int	m_nMeas {0}

std::string	m_nm
	Message source name. More...

boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...

std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...

std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...

Detailed Description

Definition at line 21 of file CalibSegmentChi2Minimizer.h.

Member Typedef Documentation

◆ HitType

using MuonR4::CalibSegmentChi2Minimizer::HitType = std::unique_ptr<CalibratedSpacePoint>

Definition at line 24 of file CalibSegmentChi2Minimizer.h.

◆ HitVec

using MuonR4::CalibSegmentChi2Minimizer::HitVec = std::vector<HitType>

Definition at line 25 of file CalibSegmentChi2Minimizer.h.

Constructor & Destructor Documentation

◆ CalibSegmentChi2Minimizer()

MuonR4::CalibSegmentChi2Minimizer::CalibSegmentChi2Minimizer	(	const std::string &	name,
		const EventContext &	ctx,
		const Amg::Transform3D &	locToGlobTrf,
		HitVec &&	hits,
		const ISpacePointCalibrator *	calibrator,
		bool	doT0Fit
	)

Constructor taking the input hits to fit.

Beam spot constraint

Definition at line 19 of file CalibSegmentChi2Minimizer.cxx.

                                                                       :
             AthMessaging{name},
             m_name{name},
              m_ctx{ctx},
             m_locToGlob{locToGlobTrf},
             m_hits{std::move(hits)},
             m_calibrator{calibrator},
             m_doT0Fit{doT0Fit} {
         
         
         for (const HitType& hit : m_hits) {
             if (hit->fitState() != CalibratedSpacePoint::State::Valid){
                 continue;
             }
             if (hit->type() == xAOD::UncalibMeasType::Other) {
                 m_nMeas +=2;
                 m_hasPhi = true;
                 continue;
             }
             m_nMeas += hit->spacePoint()->measuresEta();
             m_nMeas += hit->spacePoint()->measuresPhi();
             m_hasPhi |= hit->spacePoint()->measuresPhi();
             if (m_doT0Fit && hit->measuresTime() && hit->type() != xAOD::UncalibMeasType::MdtDriftCircleType) {
                 ++m_nMeas;
             }
         }
         if (m_nMeas - 2 - (m_hasPhi ? 2 : 0) <= 1) {
             m_doT0Fit = false;
         }
         std::ranges::sort(m_hits, 
                   [](const HitType& a, const HitType& b) { 
                     return a->positionInChamber().z() < b->positionInChamber().z();
                   });
     }

Member Function Documentation

◆ chi2Terms()

std::vector<double> MuonR4::CalibSegmentChi2Minimizer::chi2Terms ( const double * pars ) const

Returns the contribution of each measurement to the chi2.

◆ Clone()

CalibSegmentChi2Minimizer * MuonR4::CalibSegmentChi2Minimizer::Clone ( ) const

finaloverride

Clone of the object.

Definition at line 60 of file CalibSegmentChi2Minimizer.cxx.

                                                                       {
         HitVec copyHits{};
         for (const HitType& copyMe : m_hits) {
             copyHits.emplace_back(std::make_unique<CalibratedSpacePoint>(*copyMe));
         }
         return new CalibSegmentChi2Minimizer(m_name, m_ctx, m_locToGlob, std::move(copyHits), m_calibrator, m_doT0Fit);
     }

◆ DoEval()

double MuonR4::CalibSegmentChi2Minimizer::DoEval ( const double * pars ) const

finaloverride

Evaluate the chi2 for the given set of parameters.

Definition at line 97 of file CalibSegmentChi2Minimizer.cxx.

                                                                      {
       
         Parameters parVec{Parameters::Zero()};
         for (int p = 0 ; p < toInt(ParamDefs::nPars); ++p) {
             parVec[p] = pars[p];
         }
         const auto[segPos, segDir] = makeLine(parVec);
    
         const double timeDelay = pars[toInt(ParamDefs::time)];
         ATH_MSG_VERBOSE("Starting parameters  position: "<<Amg::toString(segPos)<<", direction: "<<Amg::toString(segDir)
                        <<", timeDelay: "<<timeDelay<<". Do time fit "<<m_doT0Fit);
    
         m_hits = m_calibrator->calibrate(m_ctx, std::move(m_hits), segPos, segDir, timeDelay);
      
         double chi2{0.};
         std::optional<double> arrivalTime = std::nullopt;
         if (m_doT0Fit) {
             arrivalTime = std::make_optional<double>( (m_locToGlob *segPos).mag() * c_inv + timeDelay);
         }
         for (const HitType& hit : m_hits) {            
             chi2 += SegmentFitHelpers::chiSqTerm(segPos,segDir, timeDelay, arrivalTime, *hit, msg());
         }
         ATH_MSG_VERBOSE("Final chi2: "<<chi2);
         return chi2;
     }

◆ doTimeFit()

bool MuonR4::CalibSegmentChi2Minimizer::doTimeFit ( ) const

Does the fit include the time.

Definition at line 73 of file CalibSegmentChi2Minimizer.cxx.

                                                     {
         return m_doT0Fit;
     }

◆ hasPhiMeas()

bool MuonR4::CalibSegmentChi2Minimizer::hasPhiMeas ( ) const

Are phi measurements in the collection.

Definition at line 79 of file CalibSegmentChi2Minimizer.cxx.

                                                      {
         return m_hasPhi;
     }

◆ initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

 {
   m_imsg = Athena::getMessageSvc();
   m_lvl = m_imsg ?
     static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
     MSG::INFO;
 }

◆ localToGlobTrans()

const Amg::Transform3D & MuonR4::CalibSegmentChi2Minimizer::localToGlobTrans ( ) const

Returns the local to global transformation.

Definition at line 67 of file CalibSegmentChi2Minimizer.cxx.

                                                                             {
         return m_locToGlob;
     }

◆ measurements()

const HitVec & MuonR4::CalibSegmentChi2Minimizer::measurements ( ) const

Returns the used measurements.

Definition at line 82 of file CalibSegmentChi2Minimizer.cxx.

                                                                 {
         return m_hits;
     }

◆ msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 164 of file AthMessaging.h.

 {
   MsgStream* ms = m_msg_tls.get();
   if (!ms) {
     if (!m_initialized.test_and_set()) initMessaging();
     ms = new MsgStream(m_imsg,m_nm);
     m_msg_tls.reset( ms );
   }
  
   ms->setLevel (m_lvl);
   return *ms;
 }

◆ msg() [2/2]

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

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 179 of file AthMessaging.h.

180 { return msg() << lvl; }

◆ msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns: boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

 {
   if (!m_initialized.test_and_set()) initMessaging();
   if (m_lvl <= lvl) {
     msg() << lvl;
     return true;
   } else {
     return false;
   }
 }

◆ NDim()

unsigned int MuonR4::CalibSegmentChi2Minimizer::NDim ( ) const

finaloverride

Degrees of freedom.

Definition at line 70 of file CalibSegmentChi2Minimizer.cxx.

                                                        {
         return toInt(ParamDefs::nPars);
     }

◆ nDoF()

int MuonR4::CalibSegmentChi2Minimizer::nDoF ( ) const

Returns the degrees of freedom from the measurements & beamspot constraint.

Definition at line 76 of file CalibSegmentChi2Minimizer.cxx.

                                               {
         return m_nMeas - 2 - (m_hasPhi ? 2 : 0) - m_doT0Fit;
     }

◆ release()

HitVec MuonR4::CalibSegmentChi2Minimizer::release ( const double * pars )

Releases all measurements.

Definition at line 85 of file CalibSegmentChi2Minimizer.cxx.

                                                                 {
         
         Parameters parVec{Parameters::Zero()};
         for (int p = 0 ; p < toInt(ParamDefs::nPars); ++p) {
             parVec[p] = pars[p];
         }
         const auto[segPos, segDir] = makeLine(parVec);
         const double timeDelay = pars[toInt(ParamDefs::time)];
         HitVec released = m_calibrator->calibrate(m_ctx, std::move(m_hits), segPos, segDir, timeDelay);
         m_hits.clear();
         return released;
     }

◆ segmentChi2TermMdt()

double MuonR4::CalibSegmentChi2Minimizer::segmentChi2TermMdt	(	const CalibratedSpacePoint &	hit,
		const Amg::Vector3D &	chamberIsect,
		const Amg::Vector3D &	segmentDir
	)		const

private

◆ segmentChi2TermStrip()

double MuonR4::CalibSegmentChi2Minimizer::segmentChi2TermStrip	(	const CalibratedSpacePoint &	hit,
		const Amg::Vector3D &	chamberIsect,
		const Amg::Vector3D &	segmentDir
	)		const

private

◆ setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

 {
   m_lvl = lvl;
 }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

HitVec m_hits MuonR4::CalibSegmentChi2Minimizer::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 61 of file CalibSegmentChi2Minimizer.h.

◆ m_calibrator

const ISpacePointCalibrator* MuonR4::CalibSegmentChi2Minimizer::m_calibrator {}

private

Definition at line 62 of file CalibSegmentChi2Minimizer.h.

◆ m_ctx

const EventContext& MuonR4::CalibSegmentChi2Minimizer::m_ctx

private

Definition at line 59 of file CalibSegmentChi2Minimizer.h.

◆ m_doT0Fit

bool MuonR4::CalibSegmentChi2Minimizer::m_doT0Fit {false}

private

Definition at line 63 of file CalibSegmentChi2Minimizer.h.

◆ m_hasPhi

bool MuonR4::CalibSegmentChi2Minimizer::m_hasPhi {false}

private

Definition at line 75 of file CalibSegmentChi2Minimizer.h.

◆ m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_locToGlob

Amg::Transform3D MuonR4::CalibSegmentChi2Minimizer::m_locToGlob {Amg::Transform3D::Identity()}

private

Definition at line 60 of file CalibSegmentChi2Minimizer.h.

◆ m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

◆ m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

◆ m_name

std::string MuonR4::CalibSegmentChi2Minimizer::m_name {}

private

Definition at line 58 of file CalibSegmentChi2Minimizer.h.

◆ m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_nMeas

unsigned int MuonR4::CalibSegmentChi2Minimizer::m_nMeas {0}

private

Definition at line 76 of file CalibSegmentChi2Minimizer.h.

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

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ HitType

◆ HitVec

Constructor & Destructor Documentation

◆ CalibSegmentChi2Minimizer()

Member Function Documentation

◆ chi2Terms()

◆ Clone()

◆ DoEval()

◆ doTimeFit()

◆ hasPhiMeas()

◆ initMessaging()

◆ localToGlobTrans()

◆ measurements()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ NDim()

◆ nDoF()

◆ release()

◆ segmentChi2TermMdt()

◆ segmentChi2TermStrip()

◆ setLevel()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_calibrator

◆ m_ctx

◆ m_doT0Fit

◆ m_hasPhi

◆ m_imsg

◆ m_locToGlob

◆ m_lvl

◆ m_msg_tls

◆ m_name

◆ m_nm

◆ m_nMeas