#include <MmCTPClusterCalibData.h>

Inheritance diagram for Muon::MmCTPClusterCalibData:

Collaboration diagram for Muon::MmCTPClusterCalibData:

Public Types
using	CTPParameters = std::array<double, 2>

Public Member Functions
	MmCTPClusterCalibData (const Muon::IMuonIdHelperSvc *idHelperSvc)
	~MmCTPClusterCalibData ()=default
StatusCode	storeConstants (const Identifier &gasGapId, CTPParameters &&newConstants)
double	getCTPCorrectedDriftVelocity (const Identifier &identifier, const double theta) const
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Private Types
using	parameterMap_t = std::vector<std::unique_ptr<CTPParameters>>

Private Member Functions
std::uint32_t	convertHash (const Identifier &gasGapId) const
	Converts the identifier to a continious hash used to access the stored parameters.
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
const Muon::IMuonIdHelperSvc *	m_idHelperSvc {nullptr}
parameterMap_t	m_database {}
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels).
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging).

Detailed Description

Definition at line 16 of file MmCTPClusterCalibData.h.

Member Typedef Documentation

◆ CTPParameters

using Muon::MmCTPClusterCalibData::CTPParameters = std::array<double, 2>

Definition at line 19 of file MmCTPClusterCalibData.h.

◆ parameterMap_t

using Muon::MmCTPClusterCalibData::parameterMap_t = std::vector<std::unique_ptr<CTPParameters>>

private

Definition at line 33 of file MmCTPClusterCalibData.h.

Constructor & Destructor Documentation

◆ MmCTPClusterCalibData()

Muon::MmCTPClusterCalibData::MmCTPClusterCalibData ( const Muon::IMuonIdHelperSvc * idHelperSvc )

Definition at line 8 of file MmCTPClusterCalibData.cxx.

                                                                                   :
    AthMessaging{"MmCTPClusterCalibData"},
    m_idHelperSvc{idHelperSvc} {
    m_database.resize(m_idHelperSvc->mmIdHelper().detectorElement_hash_max() * 4);
}

◆ ~MmCTPClusterCalibData()

Muon::MmCTPClusterCalibData::~MmCTPClusterCalibData ( )

default

Member Function Documentation

◆ convertHash()

std::uint32_t Muon::MmCTPClusterCalibData::convertHash ( const Identifier & gasGapId ) const

private

Converts the identifier to a continious hash used to access the stored parameters.

Parameters

gasgGapId Identifier to convert

Definition at line 16 of file MmCTPClusterCalibData.cxx.

                                                                               {
    IdentifierHash detHash{};
    const MmIdHelper& idHelper{m_idHelperSvc->mmIdHelper()};
    if(idHelper.get_detectorElement_hash(gasGapId, detHash)){
        ATH_MSG_WARNING(__func__<<"() - "<<m_idHelperSvc->toString(gasGapId)<<" is not a valid Mm identifier");
        return -1;
    }
    return static_cast<unsigned>(detHash)*4 + (idHelper.gasGap(gasGapId) -1); 
 
}

◆ getCTPCorrectedDriftVelocity()

double Muon::MmCTPClusterCalibData::getCTPCorrectedDriftVelocity	(	const Identifier &	identifier,
		const double	theta ) const

Definition at line 39 of file MmCTPClusterCalibData.cxx.

                                                                                                                       {
    //Identifier: There is no pcb segmentation for these corrections, stored with pcb = 1 as default! Therefore expects pcb = 1 
    //Theta: Expected in degrees
    //Drift velocity: will be fully taken from the parametrisation from calibration file 
 
    //If not present in the map return original value
    const auto& calibPars = m_database.at(convertHash(gasGapIdentifier));
    if(!calibPars) {
        ATH_MSG_WARNING("There's no drift velocity calibration available for gasGap " 
                << m_idHelperSvc->toStringGasGap(gasGapIdentifier)<<".");
        return std::numeric_limits<double>::max();
    }
 
    ATH_MSG_VERBOSE( "Retrieving drift velocity for stName" << m_idHelperSvc->toStringGasGap(gasGapIdentifier) );
 
    //Conversions of incident angle
    double trf_theta_in_degrees = (theta > 90.) ? (180. - theta) : theta;
    double tan_theta            = std::tan(trf_theta_in_degrees*Gaudi::Units::deg);
 
    //Parametrisation was derived as delta_residual = x_true - x_charge_weighted = (time_true - time_charge_weighted) * v_drift * tan(theta)
    //The fit used is a linear fit from data obtained for VMM 17 to VMM 95. Data outside this fit range were excluted due to problems with the magnetic field at the inner and outer PCBs, see figure 11.8 on page 135 of https://cds.cern.ch/record/2839930/files/CERN-THESIS-2021-354.pdf. Therefore, it is expected, that x^2 and x^3 are 0.
    //Parametrisation from Stefanie Gotz and Fabian Vogel, for more details see:
    //https://indico.cern.ch/event/1501492/contributions/6321472/attachments/3013528/5314031/SpatialResolution14.pdf#page=5
    //https://indico.cern.ch/event/1458120/contributions/6138857/attachments/2942828/5170893/MuonWeekOctober24FV1.pdf
    double vDrift = (*calibPars)[0] + (*calibPars)[1] * trf_theta_in_degrees;
    //Remove tantheta and always return positive drift velocity!
 
    vDrift =  (tan_theta != 0 ) ? std::abs(vDrift/tan_theta) : std::abs(vDrift);
 
 
    ATH_MSG_VERBOSE( "New drift velocity: " << vDrift  << " for theta: " << trf_theta_in_degrees << " degrees" );
 
    return vDrift;
}

◆ 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();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

◆ 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 167 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 182 of file AthMessaging.h.

183{ return msg() << lvl; }

AthMessaging::msg

MsgStream & msg() const

The standard message stream.

Definition AthMessaging.h:167

◆ 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 user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

◆ 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;
}

◆ storeConstants()

StatusCode Muon::MmCTPClusterCalibData::storeConstants	(	const Identifier &	gasGapId,
		CTPParameters &&	newConstants )

Definition at line 26 of file MmCTPClusterCalibData.cxx.

                                                                               {    
 
    auto& slot = m_database.at(convertHash(gasGapIdentifier));
    if (slot) {
        ATH_MSG_ERROR("The drift velocity calibration constants for gas gap "<<m_idHelperSvc->toStringGasGap(gasGapIdentifier) 
                    <<" already exist. Overwriting is not allowed");
        return StatusCode::FAILURE;
    }
    slot = std::make_unique<CTPParameters>(std::move(newConstants));
    return StatusCode::SUCCESS;
}

Member Data Documentation

◆ ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging).

Definition at line 141 of file AthMessaging.h.

◆ m_database

parameterMap_t Muon::MmCTPClusterCalibData::m_database {}

private

Definition at line 34 of file MmCTPClusterCalibData.h.

34{};

◆ m_idHelperSvc

const Muon::IMuonIdHelperSvc* Muon::MmCTPClusterCalibData::m_idHelperSvc {nullptr}

private

Definition at line 32 of file MmCTPClusterCalibData.h.

32{nullptr};

◆ m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

135{ nullptr };

◆ m_lvl

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

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

138{ MSG::NIL };

◆ 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_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

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

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ CTPParameters

◆ parameterMap_t

Constructor & Destructor Documentation

◆ MmCTPClusterCalibData()

◆ ~MmCTPClusterCalibData()

Member Function Documentation

◆ convertHash()

◆ getCTPCorrectedDriftVelocity()

◆ initMessaging()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ setLevel()

◆ storeConstants()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_database

◆ m_idHelperSvc

◆ m_imsg

◆ m_lvl

◆ m_msg_tls

◆ m_nm