TRTTimeCorrection Class Reference

Time correction. More...

#include <TRTTimeCorrection.h>

Inheritance diagram for TRTTimeCorrection:

Collaboration diagram for TRTTimeCorrection:

Public Member Functions
	TRTTimeCorrection (const TRTDigSettings *digset, const InDetDD::TRT_DetectorManager *detmgr, const TRT_ID *, const ITRT_CalDbTool *)
	~TRTTimeCorrection ()
double	TimeShift (const int &strawID, const InDetDD::TRT_DetElementContainer *detElements)
	Returns the time it would take to travel at light-speed from (0,0,0) to the farthest end of the wire (where the electronics are).
void	PropagationTime (const int &strawID, const double &meanZ, double &propagationTime1, double &propagationTime2)
	Calculates the time between the signal reaching the wire and when it reaches the electronics.
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.

Protected Attributes
const TRT_ID *	m_trt_id

Private Member Functions
Identifier	getIdentifier (int hitID, bool &statusok)
void	Initialize ()
double	calculateTimeShift_Barrel (const unsigned int &iPhi, const unsigned int &iRing, const unsigned int &iLayer, const unsigned int &iStraw, const int strawID, const InDetDD::TRT_DetElementContainer *detElements)
	Time shift for barrel straws.
double	calculateTimeShift_EndCap (const unsigned int &iPhi, const unsigned int &iWheel, const unsigned int &iLayer, const int strawID, const InDetDD::TRT_DetElementContainer *detElements)
	Time shift for end cap straws.
double	calculateTimeShiftFromStrawEnds (const Amg::Vector3D &strawend1_globalcoord, const Amg::Vector3D &strawend2_globalcoord, const int strawID)
	Time shift from straw endpoints in global system.
void	calculateSignalDists_Barrel (const unsigned int &iRing, const unsigned int &iLayer, double &direct_dist, double &reflect_dist) const
	Calculate the distance along the wire the signal travels before reaching the electronics.
void	calculateSignalDists_EndCap (const unsigned int &iWheel, double &direct_dist, double &reflect_dist) const
	Calculate the distance along the wire the signal travels before reaching the electronics.
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
const TRTDigSettings *	m_settings
const InDetDD::TRT_DetectorManager *	m_detmgr
std::vector< std::vector< std::vector< double > > >	m_timeShiftForEndCapPlanes
	Cached timeshifts.
std::vector< std::vector< std::vector< std::vector< double > > > >	m_timeShiftForBarrelStraws
	Cached timeshifts.
std::vector< double >	m_directDistsForEndCapWheels
	Cached distances.
std::vector< double >	m_reflectedDistsForEndCapWheels
	Cached distances.
std::vector< std::vector< double > >	m_directDistsForBarrelLayers
	Cached distances.
std::vector< std::vector< double > >	m_reflectedDistsForBarrelLayers
	Cached distances.
const unsigned int	m_subdetectorMask
const unsigned int	m_right5Bits
const unsigned int	m_shift5Bits
const unsigned int	m_shift10Bits
const unsigned int	m_shift15Bits
const double	m_notInitVal
	Value used to denote an uninitialized value.
double	m_signalPropagationSpeed = 0.0
double	m_lengthDeadRegion = 0.0
double	m_maxVertexDisplacement = 0.0
bool	m_timeShiftPhiSectSymmetry = false
bool	m_getT0FromData = false
const ITRT_CalDbTool *	m_trtcaldbtool
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

Time correction.

Definition at line 27 of file TRTTimeCorrection.h.

Constructor & Destructor Documentation

TRTTimeCorrection()

TRTTimeCorrection::TRTTimeCorrection	(	const TRTDigSettings *	digset,
		const InDetDD::TRT_DetectorManager *	detmgr,
		const TRT_ID *	trt_id,
		const ITRT_CalDbTool *	calDbTool )

Definition at line 23 of file TRTTimeCorrection.cxx.

  : AthMessaging("TRTTimeCorrection"),
    m_settings(digset), m_detmgr(detmgr), m_trt_id(trt_id),
    m_subdetectorMask(0x00200000), m_right5Bits(0x0000001F),
    m_shift5Bits(5), m_shift10Bits(10), m_shift15Bits(15), m_notInitVal(-999999.0), m_trtcaldbtool(calDbTool)
{
  Initialize();
}

~TRTTimeCorrection()

TRTTimeCorrection::~TRTTimeCorrection ( )

default

Member Function Documentation

calculateSignalDists_Barrel()

void TRTTimeCorrection::calculateSignalDists_Barrel ( const unsigned int & iRing, const unsigned int & iLayer, double & direct_dist, double & reflect_dist ) const

private

Calculate the distance along the wire the signal travels before reaching the electronics.

Both the direct and reflected signal. The signal starts in the middle of the straw.

Definition at line 337 of file TRTTimeCorrection.cxx.

                                                                                                      {
 
  //We need to calculate the distance along the wire that the signal
  //has to travel before it reaches the electronics. Both if it goes
  //directly to the electronics and if it goes via a reflection in the
  //other end of the wire.
  //
  //The signal starts in the middle of the active region.
  //
  //In addition to the length of the active gas, the signal also has
  //to go through the little dead region at the end.
 
  const InDetDD::TRT_BarrelElement * barrel_element(m_detmgr->getBarrelElement(0,//positive,
                                                                               iRing,//moduleIndex,
                                                                               0,//int phiIndex,
                                                                               iLayer));//strawLayerIndex
 
  direct_dist  = 0.5*barrel_element->strawLength() + m_lengthDeadRegion;
  reflect_dist = 1.5*barrel_element->strawLength() + 3*m_lengthDeadRegion;
}

calculateSignalDists_EndCap()

void TRTTimeCorrection::calculateSignalDists_EndCap ( const unsigned int & iWheel, double & direct_dist, double & reflect_dist ) const

private

Calculate the distance along the wire the signal travels before reaching the electronics.

Both the direct and reflected signal. The signal starts in the middle of the straw.

Definition at line 360 of file TRTTimeCorrection.cxx.

                                                                                 {
 
  //For an explanation, please read the comment in calculateSignalDists_Barrel
 
  const InDetDD::TRT_EndcapElement * ec_element(m_detmgr->getEndcapElement(0,//positive,
                                                                           iWheel,//wheelIndex,
                                                                           0,//strawLayerIndex,
                                                                           0));//phiIndex
 
  direct_dist  = 0.5*ec_element->strawLength() + m_lengthDeadRegion;
  reflect_dist = 1.5*ec_element->strawLength() + 3*m_lengthDeadRegion;
}

calculateTimeShift_Barrel()

double TRTTimeCorrection::calculateTimeShift_Barrel ( const unsigned int & iPhi, const unsigned int & iRing, const unsigned int & iLayer, const unsigned int & iStraw, const int strawID, const InDetDD::TRT_DetElementContainer * detElements )

private

Time shift for barrel straws.

Definition at line 173 of file TRTTimeCorrection.cxx.

                                                                                 {
 
  const InDetDD::TRT_BarrelElement * barrel_element(detElements->getBarrelDetElement(0/*positive*/,
                                             iRing, iPhi, iLayer ));
 
  //Sanity checks:
  if (!barrel_element) {
    ATH_MSG_ERROR("calculateTimeShift_Barrel: Could not get element for iRing = "
                  << iRing <<" and iLayer = "<<iLayer<<". Timeshift becomes 0.");
    return 0.0;
  }
 
  if (iStraw >= barrel_element->nStraws()) {
    ATH_MSG_ERROR("calculateTimeShift_Barrel: Trying to access iStraw "
                  << iStraw <<" in an element with "<<barrel_element->nStraws()<<" straws (iRing="
                  << iRing <<",iLayer="<<iLayer<<"). Timeshift becomes 0.");
    return 0.0;
  }
 
  //Straw endpoints for the straw lying along the z-axis:
  Amg::Vector3D strawend1(0,0,   barrel_element->strawLength()*0.5  );
  Amg::Vector3D strawend2(0,0, -(barrel_element->strawLength()*0.5) );
 
  //And here put into their proper global place:
  strawend1 = barrel_element->strawTransform(iStraw) * strawend1;
  strawend2 = barrel_element->strawTransform(iStraw) * strawend2;
 
  return calculateTimeShiftFromStrawEnds(strawend1,strawend2,strawID);
 
}

calculateTimeShift_EndCap()

double TRTTimeCorrection::calculateTimeShift_EndCap ( const unsigned int & iPhi, const unsigned int & iWheel, const unsigned int & iLayer, const int strawID, const InDetDD::TRT_DetElementContainer * detElements )

private

Time shift for end cap straws.

Definition at line 210 of file TRTTimeCorrection.cxx.

                                                                                {
 
  const InDetDD::TRT_EndcapElement * ec_element(detElements->getEndcapDetElement(0/*positive*/,
                                         iWheel, iLayer, iPhi ));
 
  //Sanity check:
  if (!ec_element) {
    ATH_MSG_ERROR("calculateTimeShift_EndCap: Could not get element for iWheel = "
                  << iWheel <<" and iLayer = "<<iLayer<<". Timeshift becomes 0.");
    return 0.0;
  }
 
  //Straw endpoints for the straw lying along the z-axis:
  Amg::Vector3D strawend1(0,0, ec_element->strawLength() * 0.5 );
  Amg::Vector3D strawend2(0,0, ec_element->strawLength() * (-0.5) );
 
  //And here put into their proper global place:
  strawend1 = ec_element->strawTransform(0) * strawend1;
  strawend2 = ec_element->strawTransform(0) * strawend2;
 
  return calculateTimeShiftFromStrawEnds(strawend1,strawend2,strawID); //,m_lvl
 
}

calculateTimeShiftFromStrawEnds()

double TRTTimeCorrection::calculateTimeShiftFromStrawEnds ( const Amg::Vector3D & strawend1_globalcoord, const Amg::Vector3D & strawend2_globalcoord, const int strawID )

private

Time shift from straw endpoints in global system.

Definition at line 239 of file TRTTimeCorrection.cxx.

                                                                                {
 
  //The two (hopefully relevant) extreme points of the vertex region:
  Amg::Vector3D vertexExtension1( m_settings->timeOffsetCalcVertexX(),
                                  m_settings->timeOffsetCalcVertexY(),
                                  m_settings->timeOffsetCalcVertexZ() + m_maxVertexDisplacement);
  Amg::Vector3D vertexExtension2( m_settings->timeOffsetCalcVertexX(),
                                  m_settings->timeOffsetCalcVertexY(),
                                  m_settings->timeOffsetCalcVertexZ() - m_maxVertexDisplacement);
 
  //Minimum distance between vertex region and the straw ends:
  // const double mindisttoend1(std::min(strawend1_globalcoord.distance(vertexExtension1),
  //                                  strawend1_globalcoord.distance(vertexExtension2)));
  const double mindisttoend1(std::min((strawend1_globalcoord-vertexExtension1).mag(),
                                      (strawend1_globalcoord-vertexExtension2).mag()));
  // const double mindisttoend2(std::min(strawend2_globalcoord.distance(vertexExtension1),
  //                                  strawend2_globalcoord.distance(vertexExtension2)));
  const double mindisttoend2(std::min((strawend2_globalcoord-vertexExtension1).mag(),
                                      (strawend2_globalcoord-vertexExtension2).mag()));
 
  //Just a sanity check here:
  if ( (mindisttoend2<mindisttoend1) == m_settings->electronicsAreAtFarEnd() ) {
    ATH_MSG_WARNING("It would seem that the local z-coordinate of a test straw grows TOWARDS"
                    <<" the electronics ends. This will give trouble elsewhere!!");
  }
 
  double shift = 1.0; // 1 ns (negative) overall shift for the whole TRT detector. Now set in stone.
                      // Used to be set with overallT0Shift() & overallT0ShiftShortBarrel()
                      // Note: if you change this then you need to set ToolSvc.InDetTRT_DriftFunctionTool.MCTuningShift
 
  if (m_settings->getT0FromData()) {
    bool identifierOK;
    const Identifier idStraw(getIdentifier(strawID, identifierOK));
    if (identifierOK) {
      shift = m_trtcaldbtool->getT0(idStraw);
    } else {
      ATH_MSG_ERROR("Attempt to use t0 from data failed: TRTCalDbSvc was not able to supply t0 for straw with identifier: "
                    << idStraw << ". Please set getT0FromData=false in jobOptions and run again");
    }
  }
 
  if (m_settings->electronicsAreAtFarEnd())
    return std::max(mindisttoend1,mindisttoend2) / (m_settings->distanceToTimeFactor() * CLHEP::c_light) - shift;
  else
    return std::min(mindisttoend1,mindisttoend2) / (m_settings->distanceToTimeFactor() * CLHEP::c_light) - shift;
 
}

getIdentifier()

Identifier TRTTimeCorrection::getIdentifier ( int hitID, bool & statusok )

private

Definition at line 376 of file TRTTimeCorrection.cxx.

{
  statusok = true;
 
  Identifier IdStraw;
  Identifier IdLayer;
 
  const int mask(0x0000001F);
  const int word_shift(5);
  int trtID, ringID, moduleID, layerID, strawID;
  int wheelID, planeID, sectorID;
 
  const InDetDD::TRT_BarrelElement *barrelElement;
  const InDetDD::TRT_EndcapElement *endcapElement;
 
  if ( !(hitID & 0x00200000) ) {      // barrel
    strawID   = hitID & mask;
    hitID   >>= word_shift;
    layerID   = hitID & mask;
    hitID   >>= word_shift;
    moduleID  = hitID & mask;
    hitID   >>= word_shift;
    ringID    = hitID & mask;
    trtID     = hitID >> word_shift;
 
    barrelElement =
      m_detmgr->getBarrelElement(trtID, ringID, moduleID, layerID);
    if ( barrelElement ) {
 
      IdLayer = barrelElement->identify();
      IdStraw = m_trt_id->straw_id(IdLayer, strawID);
    } else {
      ATH_MSG_ERROR("Could not find detector element for barrel identifier with "
                    << "(ipos,iring,imod,ilayer,istraw) = ("
                    << trtID << ", " << ringID << ", " << moduleID << ", "
                    << layerID << ", " << strawID << ")");
      statusok = false;
    }
  } else {                           // endcap
    strawID   = hitID & mask;
    hitID   >>= word_shift;
    planeID   = hitID & mask;
    hitID   >>= word_shift;
    sectorID  = hitID & mask;
    hitID   >>= word_shift;
    wheelID   = hitID & mask;
    trtID     = hitID >> word_shift;
 
    // change trtID (which is 2/3 for endcaps) to use 0/1 in getEndcapElement
    if (trtID == 3) trtID = 0;
    else            trtID = 1;
 
    endcapElement =
      m_detmgr->getEndcapElement(trtID, wheelID, planeID, sectorID);
 
    if ( endcapElement ) {
      IdLayer = endcapElement->identify();
      IdStraw = m_trt_id->straw_id(IdLayer, strawID);
    } else {
      ATH_MSG_ERROR("Could not find detector element for endcap identifier with "
                    << "(ipos,iwheel,isector,iplane,istraw) = ("
                    << trtID << ", " << wheelID << ", " << sectorID << ", "
                    << planeID << ", " << strawID << ")");
      ATH_MSG_ERROR("If this happens very rarely, don't be alarmed (it is a Geant4 'feature')");
      ATH_MSG_ERROR("If it happens a lot, you probably have misconfigured geometry in the sim. job.");
      statusok = false;
    }
 
  }
 
  return IdStraw;
}

Initialize()

void TRTTimeCorrection::Initialize ( )

private

Definition at line 40 of file TRTTimeCorrection.cxx.

                                   {
 
  m_signalPropagationSpeed = m_settings->signalPropagationSpeed() ;
  m_lengthDeadRegion = m_settings->lengthOfDeadRegion();
  m_maxVertexDisplacement = m_settings->maxVertexDisplacement();
  m_timeShiftPhiSectSymmetry = m_settings->timeshiftsSymmetricForPhiSectors();
  m_getT0FromData = m_settings->getT0FromData();
 
  if ( m_getT0FromData && !m_trtcaldbtool ) {
    ATH_MSG_ERROR("Could not find TRT_CalDbTool => cannot use t0 of data.");
    m_getT0FromData=false;
  }
 
  const InDetDD::TRT_Numerology *numerology(m_detmgr->getNumerology());
 
  const unsigned int nbarrelphi(m_timeShiftPhiSectSymmetry ? 1 : numerology->getNBarrelPhi());
  const unsigned int nendcapphi(m_timeShiftPhiSectSymmetry ? 1 : numerology->getNEndcapPhi());
 
  //Initialize barrel max timeshift arrays:
  m_timeShiftForBarrelStraws.resize(nbarrelphi);
  for (auto & timeShiftForBarrelStraw : m_timeShiftForBarrelStraws) {
    timeShiftForBarrelStraw.resize(numerology->getNBarrelRings());
    for (unsigned int iRing = 0; iRing < timeShiftForBarrelStraw.size(); ++iRing) {
      timeShiftForBarrelStraw[iRing].resize(numerology->getNBarrelLayers(iRing));
      for (unsigned int iLayer = 0; iLayer < timeShiftForBarrelStraw[iRing].size(); ++iLayer) {
        if (m_detmgr->getBarrelElement(0,iRing,0,iLayer)) {
          unsigned int nstraws_in_layer = m_detmgr->getBarrelElement(0,iRing,0,iLayer)->nStraws();
          timeShiftForBarrelStraw[iRing][iLayer].assign(nstraws_in_layer,m_notInitVal);
        }
      }
    }
  }
 
  //Initialize endcap max timeshift arrays:
  m_timeShiftForEndCapPlanes.resize(nendcapphi);
  for (auto & timeShiftForEndCapPlane : m_timeShiftForEndCapPlanes) {
    timeShiftForEndCapPlane.resize(numerology->getNEndcapWheels());
    for (unsigned int iWheel = 0; iWheel < timeShiftForEndCapPlane.size(); ++iWheel) {
      timeShiftForEndCapPlane[iWheel].assign(numerology->getNEndcapLayers(iWheel),m_notInitVal);
    }
  }
 
  //Initialize barrel direct/reflected distances array
  m_directDistsForBarrelLayers.resize(numerology->getNBarrelRings());
  m_reflectedDistsForBarrelLayers.resize(numerology->getNBarrelRings());
  for (unsigned int iRing = 0; iRing < m_directDistsForBarrelLayers.size(); ++iRing) {
    m_directDistsForBarrelLayers[iRing].assign(numerology->getNBarrelLayers(iRing),m_notInitVal);
    m_reflectedDistsForBarrelLayers[iRing].assign(numerology->getNBarrelLayers(iRing),m_notInitVal);
  };
 
  //Initialize endcap direct/reflected distances array
  m_directDistsForEndCapWheels.assign(numerology->getNEndcapWheels(),m_notInitVal);
  m_reflectedDistsForEndCapWheels.assign(numerology->getNEndcapWheels(),m_notInitVal);
 
}

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.

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

PropagationTime()

void TRTTimeCorrection::PropagationTime	(	const int &	strawID,
		const double &	meanZ,
		double &	propagationTime1,
		double &	propagationTime2 )

Calculates the time between the signal reaching the wire and when it reaches the electronics.

Result are two times: direct and reflected signal.

Parameters

strawID	straw ID
meanZ	z coordinate of signal in straw local coordinates (z counted positive away from electronics)
propagationTime1	propagation time - direct signal
propagationTime2	propagation time - reflected signal

Definition at line 290 of file TRTTimeCorrection.cxx.

                                                                                            {
 
  double direct_distance, reflect_distance;
 
  if (strawID & m_subdetectorMask) {
 
    //===// EndCap //===//
    const unsigned int iWheel((strawID >> m_shift15Bits) & m_right5Bits);
    direct_distance  = m_directDistsForEndCapWheels[iWheel];
    reflect_distance = m_reflectedDistsForEndCapWheels[iWheel];
 
    if (direct_distance==m_notInitVal) {
      //We need to initialize
      calculateSignalDists_EndCap(iWheel,direct_distance,reflect_distance);
      m_directDistsForEndCapWheels[iWheel] = direct_distance;
      m_reflectedDistsForEndCapWheels[iWheel] = reflect_distance;
    };
 
    //Z is counted positive AWAY from the electronics:
    propagationTime1 = (direct_distance  + meanZ) / m_signalPropagationSpeed;
    propagationTime2 = (reflect_distance - meanZ) / m_signalPropagationSpeed;
 
  } else {
 
    //===// Barrel //===//
    const unsigned int iRing((strawID >> m_shift15Bits) & m_right5Bits);
    const unsigned int iLayer((strawID >> m_shift5Bits) & m_right5Bits);
    direct_distance  = m_directDistsForBarrelLayers[iRing][iLayer];
    reflect_distance = m_reflectedDistsForBarrelLayers[iRing][iLayer];
 
    if (direct_distance==m_notInitVal) {
      //We need to initialize
      calculateSignalDists_Barrel(iRing,iLayer,direct_distance,reflect_distance);
      m_directDistsForBarrelLayers[iRing][iLayer] = direct_distance;
      m_reflectedDistsForBarrelLayers[iRing][iLayer] = reflect_distance;
    };
 
    //Z is counted positive AWAY from the electronics:
    propagationTime1 = (direct_distance  + meanZ) / m_signalPropagationSpeed;
    propagationTime2 = (reflect_distance - meanZ) / m_signalPropagationSpeed;
 
  }
 
  }

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

TimeShift()

double TRTTimeCorrection::TimeShift	(	const int &	strawID,
		const InDetDD::TRT_DetElementContainer *	detElements )

Returns the time it would take to travel at light-speed from (0,0,0) to the farthest end of the wire (where the electronics are).

Parameters

strawID

straw ID

Returns

time correction

Definition at line 97 of file TRTTimeCorrection.cxx.

                                                                                                         {
 
  //TODO: Use hit id helpers (but resolve efficiency issues first).
 
  double timeshift=0.;
 
  //Layer and phi index are needed for both endcap and barrel:
  const unsigned int iLayer((strawID >> m_shift5Bits) & m_right5Bits);
  const unsigned int iPhi(m_timeShiftPhiSectSymmetry ? 0 : ( (strawID >> m_shift10Bits) & m_right5Bits ));
 
  if (strawID & m_subdetectorMask) {
 
    //===// EndCap //===//
 
    const unsigned int iWheel((strawID >> m_shift15Bits) & m_right5Bits);
 
    //Sanity check:
    if (iPhi>=m_timeShiftForEndCapPlanes.size()||
        iWheel>=m_timeShiftForEndCapPlanes[iPhi].size()||
        iLayer>=m_timeShiftForEndCapPlanes[iPhi][iWheel].size()) {
      if (m_timeShiftPhiSectSymmetry) {
        ATH_MSG_ERROR("TimeCorrection::TimeShift: (iWheel,iLayer) = ("
                      << iWheel << ", " << iLayer << ") out of bounds! Returning 0.");
      } else {
        ATH_MSG_ERROR("TimeCorrection::TimeShift: (iPhi,iWheel,iLayer) = ("
                      << iPhi << ", " << iWheel << ", " << iLayer << ") out of bounds! Returning 0.");
      }
      return 0.0;
    }
 
    timeshift = m_timeShiftForEndCapPlanes[iPhi][iWheel][iLayer];
 
    if (timeshift==m_notInitVal) {
      //We need to initialize
      timeshift = calculateTimeShift_EndCap(iPhi,iWheel,iLayer,strawID,detElements);
      m_timeShiftForEndCapPlanes[iPhi][iWheel][iLayer] = timeshift;
    }
 
  } else {
 
    //===// Barrel //===//
 
    const unsigned int iRing((strawID >> m_shift15Bits) & m_right5Bits);
    const unsigned int iStraw(strawID & m_right5Bits);
 
    //Sanity check:
    if (iPhi>=m_timeShiftForBarrelStraws.size()||
        iRing>=m_timeShiftForBarrelStraws[iPhi].size()||
        iLayer>=m_timeShiftForBarrelStraws[iPhi][iRing].size()||
        iStraw>=m_timeShiftForBarrelStraws[iPhi][iRing][iLayer].size()) {
      if (m_timeShiftPhiSectSymmetry) {
        ATH_MSG_ERROR("TimeCorrection::TimeShift: (iRing,iLayer,iStraw) = ("
                      << iRing << ", " << iLayer << ", " << iStraw
                      << ") out of bounds! Returning 0.");
      } else {
        ATH_MSG_ERROR("TimeCorrection::TimeShift: (iPhi,iRing,iLayer,iStraw) = ("
                      << iPhi << ", " << iRing << ", " << iLayer << ", " << iStraw
                      << ") out of bounds! Returning 0.");
      }
      return 0.0;
    }
 
    timeshift = m_timeShiftForBarrelStraws[iPhi][iRing][iLayer][iStraw];
 
    if (timeshift==m_notInitVal) {
      //We need to initialize
      timeshift = calculateTimeShift_Barrel(iPhi,iRing,iLayer,iStraw,strawID,detElements);
      m_timeShiftForBarrelStraws[iPhi][iRing][iLayer][iStraw] = timeshift;
    }
 
  }
 
  return timeshift;
}

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_detmgr

const InDetDD::TRT_DetectorManager* TRTTimeCorrection::m_detmgr

private

Definition at line 109 of file TRTTimeCorrection.h.

m_directDistsForBarrelLayers

std::vector< std::vector<double> > TRTTimeCorrection::m_directDistsForBarrelLayers

private

Cached distances.

Definition at line 125 of file TRTTimeCorrection.h.

m_directDistsForEndCapWheels

std::vector<double> TRTTimeCorrection::m_directDistsForEndCapWheels

private

Cached distances.

Definition at line 121 of file TRTTimeCorrection.h.

m_getT0FromData

bool TRTTimeCorrection::m_getT0FromData = false

private

Definition at line 144 of file TRTTimeCorrection.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lengthDeadRegion

double TRTTimeCorrection::m_lengthDeadRegion = 0.0

private

Definition at line 141 of file TRTTimeCorrection.h.

m_lvl

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

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_maxVertexDisplacement

double TRTTimeCorrection::m_maxVertexDisplacement = 0.0

private

Definition at line 142 of file TRTTimeCorrection.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_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_notInitVal

const double TRTTimeCorrection::m_notInitVal

private

Value used to denote an uninitialized value.

Definition at line 137 of file TRTTimeCorrection.h.

m_reflectedDistsForBarrelLayers

std::vector< std::vector<double> > TRTTimeCorrection::m_reflectedDistsForBarrelLayers

private

Cached distances.

Definition at line 127 of file TRTTimeCorrection.h.

m_reflectedDistsForEndCapWheels

std::vector<double> TRTTimeCorrection::m_reflectedDistsForEndCapWheels

private

Cached distances.

Definition at line 123 of file TRTTimeCorrection.h.

m_right5Bits

const unsigned int TRTTimeCorrection::m_right5Bits

private

Definition at line 131 of file TRTTimeCorrection.h.

m_settings

const TRTDigSettings* TRTTimeCorrection::m_settings

private

Definition at line 108 of file TRTTimeCorrection.h.

m_shift10Bits

const unsigned int TRTTimeCorrection::m_shift10Bits

private

Definition at line 133 of file TRTTimeCorrection.h.

m_shift15Bits

const unsigned int TRTTimeCorrection::m_shift15Bits

private

Definition at line 134 of file TRTTimeCorrection.h.

m_shift5Bits

const unsigned int TRTTimeCorrection::m_shift5Bits

private

Definition at line 132 of file TRTTimeCorrection.h.

m_signalPropagationSpeed

double TRTTimeCorrection::m_signalPropagationSpeed = 0.0

private

Definition at line 140 of file TRTTimeCorrection.h.

m_subdetectorMask

const unsigned int TRTTimeCorrection::m_subdetectorMask

private

Definition at line 130 of file TRTTimeCorrection.h.

m_timeShiftForBarrelStraws

std::vector< std::vector< std::vector< std::vector<double> > > > TRTTimeCorrection::m_timeShiftForBarrelStraws

private

Cached timeshifts.

Definition at line 118 of file TRTTimeCorrection.h.

m_timeShiftForEndCapPlanes

std::vector< std::vector< std::vector<double> > > TRTTimeCorrection::m_timeShiftForEndCapPlanes

private

Cached timeshifts.

Definition at line 115 of file TRTTimeCorrection.h.

m_timeShiftPhiSectSymmetry

bool TRTTimeCorrection::m_timeShiftPhiSectSymmetry = false

private

Definition at line 143 of file TRTTimeCorrection.h.

m_trt_id

const TRT_ID* TRTTimeCorrection::m_trt_id

protected

Definition at line 111 of file TRTTimeCorrection.h.

m_trtcaldbtool

const ITRT_CalDbTool* TRTTimeCorrection::m_trtcaldbtool

private

Definition at line 146 of file TRTTimeCorrection.h.

---

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

• TRTTimeCorrection.h
• TRTTimeCorrection.cxx