#include <MdtCalibrationTool.h>

Inheritance diagram for MdtCalibrationTool:

Collaboration diagram for MdtCalibrationTool:

Public Types
using	CorrectionPtr = MuonCalib::MdtFullCalibData::CorrectionPtr

using	RtRelationPtr = MuonCalib::MdtFullCalibData::RtRelationPtr

using	TubeContainerPtr = MuonCalib::MdtFullCalibData::TubeContainerPtr

Public Member Functions
	MdtCalibrationTool (const std::string &type, const std::string &name, const IInterface *parent)
	constructor More...

virtual	~MdtCalibrationTool ()=default
	destructor More...

virtual StatusCode	initialize () override final
	initialization More...

virtual MdtCalibOutput	calibrate (const EventContext &ctx, const MdtCalibInput &hit, bool resolFromRtrack=false) const override final
	Convert the raw MDT time (+charge) into a drift radius + error. More...

virtual MdtCalibTwinOutput	calibrateTwinTubes (const EventContext &ctx, MdtCalibInput &&hit, MdtCalibInput &&twinHit) const override final
	Convert the raw MDT times of two twin hits into a Twin position (coordinate along tube) It returns whether the conversion was successful. More...

virtual double	getResolutionFromRt (const EventContext &ctx, const Identifier &module, const double time) const override final

virtual ToolSettings	getSettings () const override final

virtual const MuonCalib::MdtFullCalibData *	getCalibConstants (const EventContext &ctx, const Identifier &channelId) const override final

Private Member Functions
Muon::MdtDriftCircleStatus	driftTimeStatus (double driftTime, const MuonCalib::MdtRtRelation &rtRelation) const

Private Attributes
Gaudi::Property< int >	m_windowSetting {this, "TimeWindowSetting", timeWindowMode::Default}

Gaudi::Property< double >	m_timeWindowLowerBound {this, "TimeWindowLowerBound", 0.}

Gaudi::Property< double >	m_timeWindowUpperBound {this, "TimeWindowUpperBound", 0.}

Gaudi::Property< bool >	m_doTof {this, "DoTofCorrection", true}

Gaudi::Property< bool >	m_doProp {this, "DoPropagationCorrection", true}

Gaudi::Property< bool >	m_doTemp {this, "DoTemperatureCorrection", false}

Gaudi::Property< bool >	m_doField {this,"DoMagneticFieldCorrection", false}

Gaudi::Property< bool >	m_doSlew {this, "DoSlewingCorrection", false}

Gaudi::Property< bool >	m_doBkg {this, "DoBackgroundCorrection", false}

Gaudi::Property< bool >	m_doPropUncert {this, "DoPropagationTimeUncert", false}

ToolHandle< MuonCalib::IShiftMapTools >	m_t0ShiftTool {this, "T0ShiftTool", ""}

ToolHandle< MuonCalib::IShiftMapTools >	m_tMaxShiftTool {this, "TShiftMaxTool", ""}

Gaudi::Property< bool >	m_doT0Shift {this, "DoT0Shift", false}

Gaudi::Property< bool >	m_doTMaxShift {this, "DoTMaxShift", false}

Gaudi::Property< double >	m_unphysicalHitRadiusUpperBound {this, "UpperBoundHitRadius", 20.}

Gaudi::Property< double >	m_unphysicalHitRadiusLowerBound {this, "LowerBoundHitRadius" , 0.}

Gaudi::Property< double >	m_resTwin {this, "ResolutionTwinTube" , 1.05, "Twin tube resolution"}

SG::ReadCondHandleKey< MuonCalib::MdtCalibDataContainer >	m_calibDbKey

SG::ReadCondHandleKey< AtlasFieldCacheCondObj >	m_fieldCacheCondObjInputKey

ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

Detailed Description

the Mdt Calib Service provides, on request, the drift radius and its error, computed applying the proper calibration, for any hit in Mdt chambers

Author: Martin Woudstra, Niels van Eldik

Definition at line 35 of file MdtCalibrationTool.h.

Member Typedef Documentation

◆ CorrectionPtr

using MdtCalibrationTool::CorrectionPtr = MuonCalib::MdtFullCalibData::CorrectionPtr

Definition at line 39 of file MdtCalibrationTool.h.

◆ RtRelationPtr

using MdtCalibrationTool::RtRelationPtr = MuonCalib::MdtFullCalibData::RtRelationPtr

Definition at line 40 of file MdtCalibrationTool.h.

◆ TubeContainerPtr

using MdtCalibrationTool::TubeContainerPtr = MuonCalib::MdtFullCalibData::TubeContainerPtr

Definition at line 41 of file MdtCalibrationTool.h.

Constructor & Destructor Documentation

◆ MdtCalibrationTool()

MdtCalibrationTool::MdtCalibrationTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

constructor

Definition at line 33 of file MdtCalibrationTool.cxx.

33 :

34 base_class(type, name, parent) {}

◆ ~MdtCalibrationTool()

virtual MdtCalibrationTool::~MdtCalibrationTool ( )

virtualdefault

destructor

Member Function Documentation

◆ calibrate()

MdtCalibOutput MdtCalibrationTool::calibrate	(	const EventContext &	ctx,
		const MdtCalibInput &	hit,
		bool	resolFromRtrack = `false`
	)		const

finaloverridevirtual

Convert the raw MDT time (+charge) into a drift radius + error.

It returns whether the conversion was successful.

Parameters

[in,out]	hit	Hit must have pointer set to the MdtDigit, as well as the global hit position (including the position along the tube!)
[in]	signedTracklength	the track length from the 'triggerpoint' to the hit. It is used for the time-of-flight correction. This triggerpoint is the I.P. for ATLAS p-p collisions, typically scintillators in test-beam and cosmic teststands, and not clear yet what is it is for cosmics in ATLAS. The sign is for determining the sign of the time-of-flight correction. If a muon first passes the triggerpoint, and then the MDT tube, the sign should be positive (the case for ATLAS p-p and H8 test-beam). If a muon first passes the MDT tube, and then de triggerpoint, the sign should be negative (typically the case for cosmic-ray teststands).
[in]	triggerTime	the time of the 'triggerpoint' in ns. This is the time (measured with the same clock as the MDT TDC's) when the muon passed a known point in space: the 'triggerpoint'. For ATLAS this is 0.0 since the TDC's are synchonised w.r.t. the bunch-crossings. For H8 test-beam it is the trigger time, which is time when the muon passed the trigger scintillators. For cosmic-ray teststands it is usually also the time when the muon passed the trigger scintillators. For cosmics in ATLAS it is not clear yet.
[in]	resolFromRtrack	indicates the method to provide the resolution as a function of the distance of the reconstructed track from the anode wire instead of the drift radius

Retrieve the constants for the specific tube

calculate drift time

slewing corrections

now check whether we are outside the time window

Definition at line 107 of file MdtCalibrationTool.cxx.

                                                                          {
   
   const Identifier& id{calibIn.identify()};
   const MuonCalib::MdtFullCalibData* calibConstants = getCalibConstants(ctx, id);
   if (!calibConstants) {
      ATH_MSG_WARNING("Could not find calibration data for channel "<<m_idHelperSvc->toString(id));
      return MdtCalibOutput{};
   }
  
   // require at least the MdtRtRelation to be available
   const RtRelationPtr& rtRelation{calibConstants->rtRelation};
   // Hardcoded MDT tube radius 14.6mm here - not correct for sMDT
   // on the other hand it should be rare that a tube does not have an RT
   if(!rtRelation) {
     ATH_MSG_WARNING("No rtRelation found, cannot calibrate tube "<<m_idHelperSvc->toString(id));
     return MdtCalibOutput{};
   }
   if (!calibConstants->tubeCalib) {
      ATH_MSG_WARNING("Cannot extract the single tube calibrations for tube "<<m_idHelperSvc->toString(id));
      return MdtCalibOutput{};
   }
   const SingleTubeCalib* singleTubeData = calibConstants->tubeCalib->getCalib(id);
   if (!singleTubeData) {
      ATH_MSG_WARNING("Failed to access tubedata for " << m_idHelperSvc->toString(id));
      return MdtCalibOutput{};
   }
   
   MdtCalibOutput calibResult{};
   // correct for global t0 of rt-region
   
   calibResult.setTubeT0(singleTubeData->t0 + rtRelation->t0Global());
   calibResult.setMeanAdc(singleTubeData->adcCal);
  
   // set propagation delay
   if (m_doProp) {   
     const double propagationDistance = calibIn.signalPropagationDistance(); 
     ATH_MSG_VERBOSE("Calibration of "<<m_idHelperSvc->toString(id)<<", propagation distance: "<<propagationDistance<<" -> "
                   <<(singleTubeData->inversePropSpeed * propagationDistance));
     calibResult.setPropagationTime(singleTubeData->inversePropSpeed * propagationDistance);
   }
   
   const double driftTime = calibIn.tdc() * tdcBinSize 
                          - (m_doTof ? calibIn.timeOfFlight() : 0.)
                          - calibIn.triggerTime()
                          - calibResult.tubeT0() 
                          - calibResult.signalPropagationTime();
  
   calibResult.setDriftTime(driftTime);
   // apply corrections
   double corrTime{0.};
   const bool doCorrections = m_doField || m_doTemp || m_doBkg;
   if (doCorrections) {
     const CorrectionPtr& corrections{calibConstants->corrections};
     const RtRelationPtr& rtRelation{calibConstants->rtRelation};
     const MuonCalib::IRtRelation* rt = rtRelation->rt();
     ATH_MSG_VERBOSE("There are correction functions.");
       if (m_doSlew && corrections->slewing()) {
         double slewTime=corrections->slewing()->correction(calibResult.driftTime(), calibIn.adc());
         corrTime -= slewTime;
         calibResult.setSlewingTime(slewTime);
       }
  
    
       if (m_doField && corrections->bField()) {
         MagField::AtlasFieldCache fieldCache{};
  
         SG::ReadCondHandle readHandle{m_fieldCacheCondObjInputKey, ctx};
         if (!readHandle.isValid()) {
           THROW_EXCEPTION("calibrate: Failed to retrieve AtlasFieldCacheCondObj with key " << m_fieldCacheCondObjInputKey.key());
         }
         readHandle->getInitializedCache(fieldCache);
  
         Amg::Vector3D  globalB{Amg::Vector3D::Zero()};
         fieldCache.getField(calibIn.closestApproach().data(), globalB.data());
         const Amg::Vector2D locBField = calibIn.projectMagneticField(globalB);
         using BFieldComp = MdtCalibInput::BFieldComp;
         calibResult.setLorentzTime(corrections->bField()->correction(calibResult.driftTime(), 
                                                                      locBField[static_cast<int>(BFieldComp::alongWire)], 
                                                                      locBField[static_cast<int>(BFieldComp::alongTrack)]));
         corrTime -= calibResult.lorentzTime();
       }
       if(m_doTemp && rt && rt->hasTmaxDiff()) {
         const MdtIdHelper& id_helper{m_idHelperSvc->mdtIdHelper()};
         const int mL = id_helper.multilayer(id);
         const double tempTime = MuonCalib::RtScaleFunction(calibResult.driftTime(), mL == 2, *rt);
         calibResult.setTemperatureTime(tempTime);
         corrTime-=calibResult.temperatureTime();
       }
       // background corrections (I guess this is never active)
       if (m_doBkg && corrections->background()) {
         double bgLevel{0.};
         calibResult.setBackgroundTime(corrections->background()->correction(calibResult.driftTime(), bgLevel ));
         corrTime += calibResult.backgroundTime();
       }
   }
  
   calibResult.setDriftTime(calibResult.driftTime() + corrTime);
  
   // calculate drift radius + error
   double r{0.}, reso{0.};
   double t = calibResult.driftTime();
   double t_inrange = t;
   Muon::MdtDriftCircleStatus timeStatus = driftTimeStatus(t, *rtRelation);
   
   assert(rtRelation->rt() != nullptr);
   r = rtRelation->rt()->radius(t);
   // apply tUpper gshift
   if (m_doTMaxShift) {
     float tShift = m_tMaxShiftTool->getValue(id);
     r = rtRelation->rt()->radius( t * (1 + tShift) );
   }
   // check whether drift times are within range, if not fix them to the min/max range
   if ( t < rtRelation->rt()->tLower() ) {
     t_inrange = rtRelation->rt()->tLower();
     double rmin = rtRelation->rt()->radius( t_inrange );
     double drdt = (rtRelation->rt()->radius( t_inrange + 30. ) - rmin)/30.;
     if (timeStatus == Muon::MdtStatusBeforeSpectrum) {
       t = rtRelation->rt()->tLower() - m_timeWindowLowerBound;
     }
     // if we get here we are outside the rt range but inside the window.
     r = std::max(rmin + drdt*(t-t_inrange), m_unphysicalHitRadiusLowerBound.value());
   } else if( t > rtRelation->rt()->tUpper() ) {
     t_inrange = rtRelation->rt()->tUpper();
     double rmax = rtRelation->rt()->radius( t_inrange );
     double drdt = (rmax - rtRelation->rt()->radius( t_inrange - 30. ))/30.;
     // now check whether we are outside the time window
     if ( timeStatus == Muon::MdtStatusAfterSpectrum ) {
       t = rtRelation->rt()->tUpper() + m_timeWindowUpperBound;
     }
     // if we get here we are outside the rt range but inside the window.
     r = rmax + drdt*(t-t_inrange);
   }
  
   assert(rtRelation->rtRes() != nullptr);
   if (!resolFromRtrack) {
     reso = rtRelation->rtRes()->resolution( t_inrange );
   } else {
     const std::optional<double> tFromR = rtRelation->tr()->driftTime(std::abs(calibIn.distanceToTrack()));
     reso = rtRelation->rtRes()->resolution(tFromR.value_or(0.));
   }
   
  
   if (m_doPropUncert) {
       assert(rtRelation->rt() != nullptr);
  
       const double driftTimeUp = std::min(rtRelation->rt()->tUpper(),
                                           calibIn.tdc() * tdcBinSize 
                                         - (m_doTof ? calibIn.timeOfFlight() : 0.)
                                         - calibIn.triggerTime()
                                         - calibResult.tubeT0());
  
       const double driftTimeDn = std::max(rtRelation->rt()->tLower(),
                                           calibIn.tdc() * tdcBinSize 
                                         - (m_doTof ? calibIn.timeOfFlight() : 0.)
                                         - calibIn.triggerTime()
                                         - calibResult.tubeT0()
                                         - calibIn.tubeLength() * singleTubeData->inversePropSpeed);
  
       const double radiusUp = rtRelation->rt()->radius(driftTimeUp);
       const double radiusDn = rtRelation->rt()->radius(driftTimeDn);
       ATH_MSG_VERBOSE("Measurement "<<m_idHelperSvc->toString(calibIn.identify())
           <<" nominal drift time "<<driftTime<<", down: "<<driftTimeDn<<", up: "<<driftTimeUp
           <<" --> driftRadius: "<<r<<" pm "<<reso<<", prop-up: "<<radiusUp<<", prop-dn: "<<radiusDn
           <<" delta: "<<(radiusUp-radiusDn));
       calibResult.setDriftUncertSigProp(0.5*std::abs(radiusUp - radiusDn));
   }
  
   calibResult.setDriftRadius(r, reso);
   calibResult.setStatus(timeStatus);
   // summary
   ATH_MSG_VERBOSE( "Calibration for tube " << m_idHelperSvc->toString(id)
                    <<" passed. "<<std::endl<<"Input: "<<calibIn<<std::endl<<"Extracted calib constants: "<<calibResult<<std::endl);
   return calibResult;
 }  //end MdtCalibrationTool::calibrate

◆ calibrateTwinTubes()

MdtCalibTwinOutput MdtCalibrationTool::calibrateTwinTubes	(	const EventContext &	ctx,
		MdtCalibInput &&	hit,
		MdtCalibInput &&	twinHit
	)		const

finaloverridevirtual

Convert the raw MDT times of two twin hits into a Twin position (coordinate along tube) It returns whether the conversion was successful.

The primary hit has a smaller tdc...

Propagation time difference inside the primary tube

HVPropTime - ROPropTime = invPropSpeed*(twinZ - HVPos) - (ROPos - twinZ) = 2*invPropSpeed*twinZ - (HVPos + ROPos)

Definition at line 288 of file MdtCalibrationTool.cxx.

                                                                                          {
   
   MdtCalibOutput primResult = calibrate(ctx, primHit);
   MdtCalibOutput twinResult = calibrate(ctx, twinHit);
  
   // get Identifier and MdtReadOutElement for twin tubes
   // get 'raw' drifttimes of twin pair; we don't use timeofFlight or propagationTime cause they are irrelevant for twin coordinate
   double primdriftTime = primHit.tdc()*tdcBinSize - primResult.tubeT0();
   double twinDriftTime = twinHit.tdc()*tdcBinSize - twinResult.tubeT0();
   if (primdriftTime >= twinDriftTime) {
      ATH_MSG_VERBOSE("Swap "<<m_idHelperSvc->toString(primHit.identify())<<" & "<<m_idHelperSvc->toString(twinHit.identify())
                     <<" primDriftTime: "<<primdriftTime<<", secondTime: "<<twinDriftTime);
      std::swap(primdriftTime, twinDriftTime);
      std::swap(primResult, twinResult);
      std::swap(primHit, twinHit);
   }
   const Identifier& primId = primHit.identify();
   const Identifier& twinId = twinHit.identify();
  
   // get calibration constants from DbTool
   const MuonCalib::MdtFullCalibData* data1st = getCalibConstants(ctx, primId);
   const MuonCalib::MdtFullCalibData* data2nd = getCalibConstants(ctx, twinId);
   if (!data1st || !data2nd) {
     ATH_MSG_WARNING(__FILE__<<":"<<__LINE__<<" Failed to access calibration constants for tubes "<<
                     m_idHelperSvc->toString(primId)<<" & "<<m_idHelperSvc->toString(twinId));
     return MdtCalibTwinOutput{};
   }
   const SingleTubeCalib* calibSingleTube1st = data1st->tubeCalib->getCalib(primId);
   const SingleTubeCalib* calibSingleTube2nd = data2nd->tubeCalib->getCalib(twinId);
   if (!calibSingleTube1st || !calibSingleTube2nd) {
     ATH_MSG_WARNING(__FILE__<<":"<<__LINE__<<" Failed to access calibration constants for tubes "<<
                   m_idHelperSvc->toString(primId)<<" & "<<m_idHelperSvc->toString(twinId));
     return MdtCalibTwinOutput{};
   }
  
   // get tubelength and set HV-delay (~6ns)
   constexpr double HVdelay = 6.;
  
   double twin_timedif =  twinDriftTime - primdriftTime - calibSingleTube2nd->inversePropSpeed * twinHit.tubeLength() - HVdelay;
   
   const double tubeHalfLength = 0.5*primHit.tubeLength();
   const double zTwin = std::clamp(0.5* primHit.readOutSide()* twin_timedif / calibSingleTube1st->inversePropSpeed, 
                                   -tubeHalfLength, tubeHalfLength);
   const double errZTwin = m_resTwin / calibSingleTube1st->inversePropSpeed;
  
   ATH_MSG_VERBOSE( "Twin calibration -  tube: " << m_idHelperSvc->toString(primId)<< " twintube: " << m_idHelperSvc->toString(twinId)<<endmsg
                << " prompthit tdc = " << primHit.tdc() << "  twinhit tdc = " << twinHit.tdc()
                << " tube driftTime = " << primResult<< "  second tube driftTime = " << twinResult<<endmsg
                << " Time difference =" << twin_timedif  << " zTwin=" << zTwin<<", errorZ="<<errZTwin);
   MdtCalibTwinOutput calibResult{primHit,twinHit,primResult, twinResult};     
  
  
   calibResult.setLocZ(zTwin, errZTwin);
   return calibResult;
 } 

◆ driftTimeStatus()

Muon::MdtDriftCircleStatus MdtCalibrationTool::driftTimeStatus	(	double	driftTime,
		const MuonCalib::MdtRtRelation &	rtRelation
	)		const

private

Definition at line 350 of file MdtCalibrationTool.cxx.

                                                                                                                 {  
   if (rtRelation.rt()) {
     if(driftTime < rtRelation.rt()->tLower() - m_timeWindowLowerBound) {
         ATH_MSG_VERBOSE( " drift time outside time window "
                       << driftTime << ". Mininum time = "
                       << rtRelation.rt()->tLower() - m_timeWindowLowerBound );
         return Muon::MdtStatusBeforeSpectrum;
     } else if (driftTime > rtRelation.rt()->tUpper() + m_timeWindowUpperBound) {
         ATH_MSG_VERBOSE( " drift time outside time window "
                       << driftTime << ". Maximum time  = "
                       << rtRelation.rt()->tUpper() + m_timeWindowUpperBound);
         return Muon::MdtStatusAfterSpectrum;
     }
   } else {
     ATH_MSG_WARNING( "No valid rt relation supplied for driftTimeStatus method" );
     return Muon::MdtStatusUnDefined;
   }
   return Muon::MdtStatusDriftTime;
 }

◆ getCalibConstants()

const MuonCalib::MdtFullCalibData * MdtCalibrationTool::getCalibConstants	(	const EventContext &	ctx,
		const Identifier &	channelId
	)		const

finaloverridevirtual

Definition at line 99 of file MdtCalibrationTool.cxx.

                                                                                                             {
     SG::ReadCondHandle constantHandle{m_calibDbKey, ctx};
     if (!constantHandle.isValid()){
          THROW_EXCEPTION("Failed to retrieve the Mdt calibration constants "<<m_calibDbKey.fullKey());
     }
     return constantHandle->getCalibData(channelId, msgStream());
 }

◆ getResolutionFromRt()

double MdtCalibrationTool::getResolutionFromRt	(	const EventContext &	ctx,
		const Identifier &	module,
		const double	time
	)		const

finaloverridevirtual

Definition at line 370 of file MdtCalibrationTool.cxx.

                                                                                                                             {
   
   const MuonCalib::MdtFullCalibData* moduleConstants = getCalibConstants(ctx, moduleID);
   if (!moduleConstants){
       THROW_EXCEPTION("Failed to retrieve set of calibration constants for "<<m_idHelperSvc->toString(moduleID));
   }
   const RtRelationPtr& rtRel{moduleConstants->rtRelation};
   if (!rtRel) {
     THROW_EXCEPTION("No rt-relation found for "<<m_idHelperSvc->toString(moduleID));
   }
   const double t = std::min(std::max(time, rtRel->rt()->tLower()), rtRel->rt()->tUpper());
   return rtRel->rtRes()->resolution(t);
 }

◆ getSettings()

ToolSettings MdtCalibrationTool::getSettings ( ) const

finaloverridevirtual

Definition at line 36 of file MdtCalibrationTool.cxx.

                                                    {
     ToolSettings settings{};
     using Property = ToolSettings::Property;
     settings.setBit(Property::TofCorrection, m_doTof);
     settings.setBit(Property::PropCorrection, m_doProp);
     settings.setBit(Property::TempCorrection, m_doTemp);
     settings.setBit(Property::MagFieldCorrection, m_doField);
     settings.setBit(Property::SlewCorrection, m_doSlew);
     settings.setBit(Property::BackgroundCorrection, m_doBkg);
     settings.window = static_cast<timeWindowMode>(m_windowSetting.value()); 
     return settings;
 }

◆ initialize()

StatusCode MdtCalibrationTool::initialize ( )

finaloverridevirtual

initialization

Ensure that the conditions dependency is properly declared

Shifting tools to evaluate systematic uncertainties on the T0 timing

Definition at line 48 of file MdtCalibrationTool.cxx.

                                           {  
   ATH_MSG_DEBUG( "Initializing" );
  
   switch(m_windowSetting.value()) {
     case timeWindowMode::UserDefined:
        ATH_MSG_DEBUG("Use predefined user values of "<<m_timeWindowLowerBound<<" & "<<m_timeWindowUpperBound);
        break;
     case timeWindowMode::Default:
       ATH_MSG_DEBUG("Use 1000. & 2000. as the lower and upper time window values ");
       m_timeWindowLowerBound = 1000.;
       m_timeWindowUpperBound = 2000.;
       break;
     case timeWindowMode::CollisionG4:
        ATH_MSG_DEBUG("Use Geant4 collision time window of 20-30");
        m_timeWindowLowerBound = 20.;
        m_timeWindowUpperBound = 30.;
        break;
     case timeWindowMode::CollisionData:
         ATH_MSG_DEBUG("Use collision data time window of 10 to 30");
         m_timeWindowLowerBound = 10.;
         m_timeWindowUpperBound = 30.;
         break;
     case timeWindowMode::CollisionFitT0:
         ATH_MSG_DEBUG("Use collision data time window of 50 to 100 to fit T0 in the end");
         m_timeWindowLowerBound = 50.;
         m_timeWindowUpperBound = 100.;
         break;
     default:
        ATH_MSG_FATAL("Unknown time window setting "<<m_windowSetting<<" provided.");
        return StatusCode::FAILURE;
   };
  
   ATH_CHECK(m_idHelperSvc.retrieve());
   ATH_CHECK(m_fieldCacheCondObjInputKey.initialize());
   ATH_CHECK(m_calibDbKey.initialize());
   ATH_CHECK(m_t0ShiftTool.retrieve(EnableTool{m_doT0Shift}));
   ATH_CHECK(m_tMaxShiftTool.retrieve(EnableTool{m_doTMaxShift}));
   ATH_MSG_DEBUG("Initialization finalized "<<std::endl
                 <<"  TimeWindow: ["<< m_timeWindowLowerBound.value()<<";"<<m_timeWindowUpperBound.value()<<"]"<<std::endl
                 <<"   Correct time of flight "<<(m_doTof ? "yay" : "nay")<<std::endl
                 <<"   Correct propagation time "<<(m_doProp ? "si" : "no")<<std::endl
                 <<"   Correct temperature "<<(m_doTemp ? "si" : "no")<<std::endl
                 <<"   Correct magnetic field "<<(m_doField ? "si" : "no")<<std::endl
                 <<"   Correct time slew "<<(m_doSlew ? "si" : "no")<<std::endl
                 <<"   Correct background "<<(m_doBkg ? "si" : "no"));
   return StatusCode::SUCCESS;
 }  //end MdtCalibrationTool::initialize

Member Data Documentation

◆ m_calibDbKey

SG::ReadCondHandleKey<MuonCalib::MdtCalibDataContainer> MdtCalibrationTool::m_calibDbKey

private

Initial value:

{this, "CalibDataKey", "MdtCalibConstants",

"Conditions object containing the calibrations"}

Definition at line 127 of file MdtCalibrationTool.h.

◆ m_doBkg

Gaudi::Property<bool> MdtCalibrationTool::m_doBkg {this, "DoBackgroundCorrection", false}

private

Definition at line 111 of file MdtCalibrationTool.h.

◆ m_doField

Gaudi::Property<bool> MdtCalibrationTool::m_doField {this,"DoMagneticFieldCorrection", false}

private

Definition at line 109 of file MdtCalibrationTool.h.

◆ m_doProp

Gaudi::Property<bool> MdtCalibrationTool::m_doProp {this, "DoPropagationCorrection", true}

private

Definition at line 107 of file MdtCalibrationTool.h.

◆ m_doPropUncert

Gaudi::Property<bool> MdtCalibrationTool::m_doPropUncert {this, "DoPropagationTimeUncert", false}

private

Definition at line 112 of file MdtCalibrationTool.h.

◆ m_doSlew

Gaudi::Property<bool> MdtCalibrationTool::m_doSlew {this, "DoSlewingCorrection", false}

private

Definition at line 110 of file MdtCalibrationTool.h.

◆ m_doT0Shift

Gaudi::Property<bool> MdtCalibrationTool::m_doT0Shift {this, "DoT0Shift", false}

private

Definition at line 119 of file MdtCalibrationTool.h.

◆ m_doTemp

Gaudi::Property<bool> MdtCalibrationTool::m_doTemp {this, "DoTemperatureCorrection", false}

private

Definition at line 108 of file MdtCalibrationTool.h.

◆ m_doTMaxShift

Gaudi::Property<bool> MdtCalibrationTool::m_doTMaxShift {this, "DoTMaxShift", false}

private

Definition at line 120 of file MdtCalibrationTool.h.

◆ m_doTof

Gaudi::Property<bool> MdtCalibrationTool::m_doTof {this, "DoTofCorrection", true}

private

Definition at line 106 of file MdtCalibrationTool.h.

◆ m_fieldCacheCondObjInputKey

SG::ReadCondHandleKey<AtlasFieldCacheCondObj> MdtCalibrationTool::m_fieldCacheCondObjInputKey

private

Initial value:

{this, "AtlasFieldCacheCondObj", "fieldCondObj",

"Name of the Magnetic Field conditions object key"}

Definition at line 134 of file MdtCalibrationTool.h.

◆ m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MdtCalibrationTool::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 136 of file MdtCalibrationTool.h.

◆ m_resTwin

Gaudi::Property<double> MdtCalibrationTool::m_resTwin {this, "ResolutionTwinTube" , 1.05, "Twin tube resolution"}

private

Definition at line 124 of file MdtCalibrationTool.h.

◆ m_t0ShiftTool

ToolHandle<MuonCalib::IShiftMapTools> MdtCalibrationTool::m_t0ShiftTool {this, "T0ShiftTool", ""}

private

Definition at line 114 of file MdtCalibrationTool.h.

◆ m_timeWindowLowerBound

Gaudi::Property<double> MdtCalibrationTool::m_timeWindowLowerBound {this, "TimeWindowLowerBound", 0.}

private

Definition at line 104 of file MdtCalibrationTool.h.

◆ m_timeWindowUpperBound

Gaudi::Property<double> MdtCalibrationTool::m_timeWindowUpperBound {this, "TimeWindowUpperBound", 0.}

private

Definition at line 105 of file MdtCalibrationTool.h.

◆ m_tMaxShiftTool

ToolHandle<MuonCalib::IShiftMapTools> MdtCalibrationTool::m_tMaxShiftTool {this, "TShiftMaxTool", ""}

private

Definition at line 116 of file MdtCalibrationTool.h.

◆ m_unphysicalHitRadiusLowerBound

Gaudi::Property<double> MdtCalibrationTool::m_unphysicalHitRadiusLowerBound {this, "LowerBoundHitRadius" , 0.}

private

Definition at line 123 of file MdtCalibrationTool.h.

◆ m_unphysicalHitRadiusUpperBound

Gaudi::Property<double> MdtCalibrationTool::m_unphysicalHitRadiusUpperBound {this, "UpperBoundHitRadius", 20.}

private

Definition at line 122 of file MdtCalibrationTool.h.

◆ m_windowSetting

Gaudi::Property<int> MdtCalibrationTool::m_windowSetting {this, "TimeWindowSetting", timeWindowMode::Default}

private

Definition at line 103 of file MdtCalibrationTool.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

◆ CorrectionPtr

◆ RtRelationPtr

◆ TubeContainerPtr

Constructor & Destructor Documentation

◆ MdtCalibrationTool()

◆ ~MdtCalibrationTool()

Member Function Documentation

◆ calibrate()

◆ calibrateTwinTubes()

◆ driftTimeStatus()

◆ getCalibConstants()

◆ getResolutionFromRt()

◆ getSettings()

◆ initialize()

Member Data Documentation

◆ m_calibDbKey

◆ m_doBkg

◆ m_doField

◆ m_doProp

◆ m_doPropUncert

◆ m_doSlew

◆ m_doT0Shift

◆ m_doTemp

◆ m_doTMaxShift

◆ m_doTof

◆ m_fieldCacheCondObjInputKey

◆ m_idHelperSvc

◆ m_resTwin

◆ m_t0ShiftTool

◆ m_timeWindowLowerBound

◆ m_timeWindowUpperBound

◆ m_tMaxShiftTool

◆ m_unphysicalHitRadiusLowerBound

◆ m_unphysicalHitRadiusUpperBound

◆ m_windowSetting