ATLAS Offline Software
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
TRT_ToT_dEdx Class Referencefinal

#include <TRT_ToT_dEdx.h>

Inheritance diagram for TRT_ToT_dEdx:
Collaboration diagram for TRT_ToT_dEdx:

Public Types

enum  EDataBaseType { kOldDB, kNewDB }
 
enum  EstCalc { kAlgStandard, kAlgReweight, kAlgReweightTrunkOne }
 
enum  EOccupancyCorrection { kRSOnly, kHitBased, kTrackBased, kGlobal }
 
enum  EGasType { kXenon, kArgon, kKrypton, kUnset }
 

Public Member Functions

 TRT_ToT_dEdx (const std::string &, const std::string &, const IInterface *)
 AlgTool like constructor. More...
 
virtual ~TRT_ToT_dEdx ()
 Virtual destructor. More...
 
virtual StatusCode initialize () override
 AlgTool initailize method. More...
 
virtual double dEdx (const EventContext &ctx, const Trk::Track *track, bool useHThits, std::optional< float > localOccupancy=std::nullopt) const override final
 function to calculate sum ToT normalised to number of used hits More...
 
virtual double usedHits (const EventContext &ctx, const Trk::Track *track, bool useHThits=true) const override final
 function to calculate number of used hits More...
 
virtual double getTest (const EventContext &ctx, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, Trk::ParticleHypothesis antihypothesis, int nUsedHits) const override final
 function to calculate likelihood ratio test More...
 
void setStatusCorrection (bool value)
 
bool getStatusCorrection () const
 
void setCorrectionType (EOccupancyCorrection value)
 
bool getCorrectionType () const
 
void setMinRtrack (float minRtrack)
 
float getMinRtrack () const
 
void setMaxRtrack (float maxRtrack)
 
float getMaxRtrack () const
 
void setStatusUseZeroRHitCut (bool value)
 
bool getStatusUseZeroRHitCut () const
 
void setGasTypeFordEdXCalculation (EGasType gasType)
 
int getGasTypeFordEdXCalculation () const
 
void setAlgorithm (EstCalc alg)
 
int getAlgorithm () const
 
virtual double dEdx (const EventContext &ctx, const Trk::Track *track, bool useHitsHT, std::optional< float > localOccupancy=std::nullopt) const=0
 function to calculate sum ToT normalised to number of used hits More...
 
virtual double dEdx (const Trk::Track *track, bool useHitsHT, std::optional< float > localOccupancy=std::nullopt) const
 function to calculate sum ToT normalised to number of used hits More...
 
virtual double usedHits (const EventContext &ctx, const Trk::Track *track, bool useHitsHT=true) const=0
 function to calculate number of used hits More...
 
virtual double usedHits (const Trk::Track *track, bool useHitsHT=true) const
 
virtual double getTest (const EventContext &ctx, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, Trk::ParticleHypothesis antihypothesis, int nUsedHits) const=0
 function to calculate likelihood ratio test More...
 
virtual double getTest (const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, Trk::ParticleHypothesis antihypothesis, int nUsedHits) const
 
virtual double dEdx (const Trk::Track *track, bool useHitsHT, std::optional< float > localOccupancy=std::nullopt) const
 
virtual double usedHits (const Trk::Track *track, bool useHitsHT=true) const
 
virtual double getTest (const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, Trk::ParticleHypothesis antihypothesis, int nUsedHits) const
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Static Public Member Functions

static double calculateTrackLengthInStraw (const Trk::TrackStateOnSurface *trackState, const TRT_ID *identifier)
 
static const InterfaceID & interfaceID ()
 AlgTool interface methods. More...
 

Protected Member Functions

bool isGoodHit (const EventContext &ctx, const Trk::TrackStateOnSurface *trackState, bool useHitsHT, double &length) const
 function to define what is a good hit to be used for dEdx calculation cuts on track level can be made latekAlgStandardr by the user. More...
 
double correctNormalization (const EventContext &ctx, double nVtx=-1) const
 correct overall dEdx normalization on track level More...
 
double getProb (const EventContext &ctx, const Trk::TrackStateOnSurface *itr, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, int nUsedHits) const
 function to calculate likelihood from prediction and resolution More...
 
double getProb (const EventContext &ctx, EGasType gasType, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, int nUsedHits) const
 
double predictdEdx (const EventContext &ctx, const Trk::TrackStateOnSurface *itr, const double pTrk, Trk::ParticleHypothesis hypothesis) const
 function to calculate expectation value for dEdx using BB fit More...
 
double predictdEdx (const EventContext &ctx, EGasType gasType, const double pTrk, Trk::ParticleHypothesis hypothesis) const
 
double mass (const EventContext &ctx, const Trk::TrackStateOnSurface *itr, const double pTrk, double dEdx) const
 function to extract most likely mass in bg [0,3] More...
 
double correctToT_corrRZ (const EventContext &ctx, const Trk::TrackStateOnSurface *itr, double length) const
 main function to correct ToT values on hit level as a function of track radius and z-position More...
 
EGasType gasTypeInStraw (const EventContext &ctx, const Trk::TrackStateOnSurface *itr) const
 return gas type for that hit More...
 
EGasType gasTypeInStraw (const EventContext &ctx, const InDet::TRT_DriftCircleOnTrack *driftcircle) const
 
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

double fitFuncEndcap_corrRZ (const EventContext &ctx, EGasType gas, double driftRadius, double rPosition, int Layer, int sign) const
 function to compute correction factor in endcap region More...
 
double fitFuncBarrel_corrRZ (const EventContext &ctx, EGasType gas, double driftRadius, double zPosition, int Layer, int StrawLayer) const
 function to compute correction factor in barrel region More...
 
double fitFuncBarrelLong_corrRZ (const EventContext &ctx, EGasType gasType, double driftRadius, double zPosition, int Layer, int StrawLayer) const
 function called by fitFuncBarrel_corrRZ for long straws More...
 
double fitFuncBarrelShort_corrRZ (const EventContext &ctx, EGasType gasType, double driftRadius, double zPosition, int StrawLayer) const
 function called by fitFuncBarrel_corrRZ for short straws More...
 
double fitFuncPol_corrRZ (const EventContext &ctx, EGasType gasType, int parameter, double driftRadius, int Layer, int Strawlayer, int sign, int set) const
 function called by fitFuncBarrel_corrRZ and fitFuncEndcap_corrRZ More...
 
double fitFuncEndcap_corrRZL (const EventContext &ctx, EGasType gasType, double driftRadius, double radialPosition, int Layer, int sign) const
 function to compute correction factor in endcap region More...
 
double fitFuncBarrel_corrRZL (const EventContext &ctx, EGasType gasType, double driftRadius, double zPosition, int Layer, int StrawLayer) const
 function to compute correction factor in barrel region
More...
 
double hitOccupancyCorrection (const EventContext &ctx, const Trk::TrackStateOnSurface *itr) const
 
double trackOccupancyCorrection (const EventContext &ctx, const Trk::Track *track, bool useHThits, std::optional< float > localOccupancy) const
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

SG::ReadDecorHandleKey< xAOD::EventInfo > m_rdhkEvtInfo
 
const TRT_ID * m_trtId {nullptr}
 
ToolHandle< ITRT_StrawStatusSummaryTool > m_TRTStrawSummaryTool {this, "TRTStrawSummaryTool", "TRT_StrawStatusSummaryTool"}
 
ToolHandle< Trk::IPRD_AssociationTool > m_assoTool {this, "AssociationTool","InDet::InDetPRD_AssociationToolGangedPixels","PRD_AssociationTool"}
 
ToolHandle< InDet::ITRT_LocalOccupancy > m_localOccTool {this, "TRT_LocalOccupancyTool","", "Local occupancy tool"}
 
BooleanProperty m_corrected
 
BooleanProperty m_divideByL
 
IntegerProperty m_useTrackPartWithGasType
 
IntegerProperty m_toolScenario
 
IntegerProperty m_correctionType {this, "TRT_dEdx_correctionType", kTrackBased, "Type of dEdx correction"}
 
BooleanProperty m_isData {this, "TRT_dEdx_isData", true}
 
FloatProperty m_trackConfig_maxRtrack {this, "TRT_dEdx_trackConfig_maxRtrack", 1.85, "maximum track radius"}
 
FloatProperty m_trackConfig_minRtrack {this, "TRT_dEdx_trackConfig_minRtrack", 0.15, "maximum track radius"}
 
BooleanProperty m_useZeroRHitCut
 
unsigned int m_nTrunkateHits = 1
 
SG::ReadCondHandleKey< TRTDedxcorrection > m_ReadKey {this,"Dedxcorrection","Dedxcorrection","Dedx constants in-key"}
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 48 of file TRT_ToT_dEdx.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ EDataBaseType

enum TRT_ToT_dEdx::EDataBaseType
Enumerator
kOldDB 
kNewDB 

Definition at line 54 of file TRT_ToT_dEdx.h.

54 {kOldDB,kNewDB};

◆ EGasType

enum TRT_ToT_dEdx::EGasType
Enumerator
kXenon 
kArgon 
kKrypton 
kUnset 

Definition at line 57 of file TRT_ToT_dEdx.h.

57 {kXenon,kArgon,kKrypton,kUnset};

◆ EOccupancyCorrection

enum TRT_ToT_dEdx::EOccupancyCorrection
Enumerator
kRSOnly 
kHitBased 
kTrackBased 
kGlobal 

Definition at line 56 of file TRT_ToT_dEdx.h.

56 {kRSOnly, kHitBased, kTrackBased, kGlobal};

◆ EstCalc

enum TRT_ToT_dEdx::EstCalc
Enumerator
kAlgStandard 
kAlgReweight 
kAlgReweightTrunkOne 

Definition at line 55 of file TRT_ToT_dEdx.h.

55 {kAlgStandard,kAlgReweight,kAlgReweightTrunkOne};

Constructor & Destructor Documentation

◆ TRT_ToT_dEdx()

TRT_ToT_dEdx::TRT_ToT_dEdx ( const std::string &  t,
const std::string &  n,
const IInterface *  p 
)

AlgTool like constructor.

Definition at line 35 of file TRT_ToT_dEdx.cxx.

35  : AthAlgTool(t,n,p)
36 {
37  declareInterface<ITRT_ToT_dEdx>(this);
38 }

◆ ~TRT_ToT_dEdx()

TRT_ToT_dEdx::~TRT_ToT_dEdx ( )
virtualdefault

Virtual destructor.

Member Function Documentation

◆ calculateTrackLengthInStraw()

double TRT_ToT_dEdx::calculateTrackLengthInStraw ( const Trk::TrackStateOnSurface *  trackState,
const TRT_ID *  identifier 
)
static

Definition at line 1185 of file TRT_ToT_dEdx.cxx.

1185  {
1186  if (trackState->type(Trk::TrackStateOnSurface::Outlier)) return 0.; //Outliers
1187 
1188  const Trk::MeasurementBase* trkM = trackState->measurementOnTrack();
1189  if (!trkM) {
1190  return 0.;
1191  }
1192 
1193  // Check if this is RIO on track
1194  // and if yes check if is TRT Drift Circle
1195  // then set the ptr
1196  const InDet::TRT_DriftCircleOnTrack* driftcircle = nullptr;
1197  if (trkM->type(Trk::MeasurementBaseType::RIO_OnTrack)) {
1198  const Trk::RIO_OnTrack* tmpRio = static_cast<const Trk::RIO_OnTrack*>(trkM);
1199  if (tmpRio->rioType(Trk::RIO_OnTrackType::TRT_DriftCircle)) {
1200  driftcircle = static_cast<const InDet::TRT_DriftCircleOnTrack*>(tmpRio);
1201  }
1202  }
1203 
1204  if (!driftcircle) {
1205  return 0.;
1206  }
1207 
1208  const Trk::TrackParameters* trkP = trackState->trackParameters();
1209  if(trkP==nullptr) return 0.;
1210 
1211  double Trt_Rtrack = std::abs(trkP->parameters()[Trk::locR]);
1212  double Trt_HitTheta = trkP->parameters()[Trk::theta];
1213  double Trt_HitPhi = trkP->parameters()[Trk::phi];
1214  Identifier DCId = driftcircle->identify();
1215  int HitPart = std::abs(identifier->barrel_ec(DCId));
1216  const InDetDD::TRT_BaseElement* element = driftcircle->detectorElement();
1217  double strawphi = element->center(DCId).phi();
1218 
1219  // check if track is an outlier
1220  if (Trt_Rtrack >= 2.0) {
1221  return 0.;
1222  }
1223 
1224  double length=0;
1225  if (HitPart == 1) { //Barrel
1226  length = 2*std::sqrt(4-Trt_Rtrack*Trt_Rtrack)*1./std::abs(std::sin(Trt_HitTheta));
1227  } else if (HitPart == 2) { //EndCap
1228  length = 2*std::sqrt(4-Trt_Rtrack*Trt_Rtrack)*1./std::sqrt(1-std::sin(Trt_HitTheta)*std::sin(Trt_HitTheta)*std::cos(Trt_HitPhi-strawphi)*std::cos(Trt_HitPhi-strawphi));
1229  } else {
1230  // This should never happen
1231  throw std::runtime_error("Unknown barrel/endcap identifier: " + std::to_string(HitPart) + ". Must be 1(Barrel) or 2(Endcap)");
1232  }
1233 
1234  return length;
1235 }

◆ correctNormalization()

double TRT_ToT_dEdx::correctNormalization ( const EventContext &  ctx,
double  nVtx = -1 
) const
protected

correct overall dEdx normalization on track level

Parameters
numberof primary vertices per event
Returns
scaling variable

Definition at line 684 of file TRT_ToT_dEdx.cxx.

685 {
686  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
687  const TRTDedxcorrection* dEdxCorrection{*readHandle};
688  if(dEdxCorrection==nullptr) {
689  ATH_MSG_ERROR(" correctNormalization: Could not find any dEdxCorrection in CondStore. Return zero.");
690  return 0;
691  }
692 
693  EGasType gasType = static_cast<EGasType> (m_useTrackPartWithGasType.value());
694  if (m_useTrackPartWithGasType==kUnset) gasType=kXenon;
695  if (nVtx<=0) nVtx=dEdxCorrection->normNzero[gasType];
696  double slope = dEdxCorrection->normSlopeTot[gasType];
697  double offset = dEdxCorrection->normOffsetTot[gasType];
698  if (m_divideByL){
699  slope = dEdxCorrection->normSlopeTotDivideByLength[gasType];
700  offset = dEdxCorrection->normOffsetTotDivideByLength[gasType];
701  }
702  double shift = dEdxCorrection->normOffsetData[gasType];
703  if(!m_isData)shift = 0;
704  return (slope*dEdxCorrection->normNzero[gasType]+offset)/(slope*nVtx+offset+shift);
705 }

◆ correctToT_corrRZ()

double TRT_ToT_dEdx::correctToT_corrRZ ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  itr,
double  length 
) const
protected

main function to correct ToT values on hit level as a function of track radius and z-position

Parameters
trackon surface object
boolvariable to specify whether ToT or ToT/L correction
boolvariable to specify whether data or MC correction
boolvariable whether correction should actually be applied
boolvariable whether mimic ToT to other gas hits shoule be used
Returns
corrected value for ToT
Todo:
implement possiblity to set the scaling factor run-by-run from database, should probably be done later on track- level

Definition at line 708 of file TRT_ToT_dEdx.cxx.

711 {
712  const Trk::MeasurementBase* trkM = itr->measurementOnTrack();
713  const Trk::TrackParameters* trkP = itr->trackParameters();
714 
715  // Check if this is RIO on track
716  // annd if yes check if is TRT Drift Circle
717  // then set the ptr
718  const InDet::TRT_DriftCircleOnTrack* driftcircle = nullptr;
719  if (trkM && trkM->type(Trk::MeasurementBaseType::RIO_OnTrack)) {
720  const Trk::RIO_OnTrack* tmpRio = static_cast<const Trk::RIO_OnTrack*>(trkM);
721  if (tmpRio->rioType(Trk::RIO_OnTrackType::TRT_DriftCircle)) {
722  driftcircle = static_cast<const InDet::TRT_DriftCircleOnTrack*>(tmpRio);
723  }
724  }
725 
726  if (!driftcircle) {
727  return 0;
728  }
729  if (driftcircle->prepRawData()==nullptr) {
730  return 0;
731  }
732 
733  Identifier DCId = driftcircle->identify();
734  double timeOverThreshold = driftcircle->prepRawData()->timeOverThreshold();
735  if(timeOverThreshold==0) {
736  ATH_MSG_WARNING("correctToT_corrRZ(const Trk::TrackStateOnSurface *itr):: ToT="<<timeOverThreshold<<". We must cut that hit in isGoodHit() !");
737  return 0;
738  }
739  int hitPart = m_trtId->barrel_ec(DCId);
740  int StrawLayer = m_trtId->straw_layer(DCId);
741  int Layer = m_trtId->layer_or_wheel(DCId);
742  double hitRtrack = std::abs(trkP->parameters()[Trk::locR]);
743  EGasType gasType = gasTypeInStraw(ctx,itr);
744  if(gasType==kUnset) {
745  ATH_MSG_ERROR("correctToT_corrRZ(const Trk::TrackStateOnSurface *itr):: Gas type in straw is kUnset! Return ToT = 0");
746  return 0;
747  }
748 
749  if(m_divideByL && length>0) timeOverThreshold = timeOverThreshold/length;
750  if(!m_corrected) return timeOverThreshold;
751  /* else correct */
752 
753  double hitZ = driftcircle->globalPosition().z();
754  double trackx = driftcircle->globalPosition().x();
755  double tracky = driftcircle->globalPosition().y();
756  double hitPosR = std::sqrt(trackx*trackx+tracky*tracky);
757 
760  double ToTmip = 1;
761  double valToT = 0;
762  if(m_divideByL){
763  if (abs(hitPart)==1) // Barrel
764  valToT = fitFuncBarrel_corrRZL(ctx,gasType, hitRtrack, hitZ, Layer, StrawLayer);
765  else // End-cap
766  valToT = fitFuncEndcap_corrRZL(ctx,gasType, hitRtrack, hitPosR, Layer, hitZ>0?1:(hitZ<0?-1:0));
767  }else{
768  if (abs(hitPart)==1) // Barrel
769  valToT = fitFuncBarrel_corrRZ(ctx,gasType, hitRtrack, hitZ, Layer, StrawLayer);
770  else // End-cap
771  valToT = fitFuncEndcap_corrRZ(ctx,gasType, hitRtrack, hitPosR, Layer, hitZ>0?1:(hitZ<0?-1:0));
772  }
773  if (std::isinf(valToT)) return 0.;
774  if (valToT!=0) return ToTmip*timeOverThreshold/valToT;
775  return 0.;
776 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType &   
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase &  hndl,
const std::string &  doc,
const SG::VarHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey &  hndl,
const std::string &  doc,
const SG::VarHandleKeyType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray &  hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType &   
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ dEdx() [1/4]

virtual double ITRT_ToT_dEdx::dEdx

function to calculate sum ToT normalised to number of used hits

Parameters
trackpointer to track
useHitsHTdecide if HT hits should be used in the estimate
Returns
dEdx value

◆ dEdx() [2/4]

double TRT_ToT_dEdx::dEdx ( const EventContext &  ctx,
const Trk::Track *  track,
bool  useHitsHT,
std::optional< float >  localOccupancy = std::nullopt 
) const
finaloverridevirtual

function to calculate sum ToT normalised to number of used hits

Parameters
Eventcontext ctx
trackpointer to track
useHitsHTdecide if HT hits should be used in the estimate
Returns
dEdx value

Implements ITRT_ToT_dEdx.

Definition at line 133 of file TRT_ToT_dEdx.cxx.

137 {
138  ATH_MSG_DEBUG("dEdx()");
139 
140  double nVtx=-1.;
141  // Event information
142  SG::ReadDecorHandle<xAOD::EventInfo,float> eventInfoDecor(m_rdhkEvtInfo,ctx);
143  if(!eventInfoDecor.isPresent()) {
144  REPORT_MESSAGE(MSG::FATAL) << "EventInfo decoration not available!";
145  return 0;
146  }
147 
148  // Average interactions per crossing for the current BCID
149  // TODO: we should really not hard-code a mu to nVtx conversion
150  double mu = eventInfoDecor(0);
151  if(m_isData) {
152  nVtx = 1.3129 + 0.716194*mu + (-0.00475074)*mu*mu;
153  }
154  else {
155  nVtx = 1.0897 + 0.748287*mu + (-0.00421788)*mu*mu;
156  }
157 
158  if (!track) {
159  return 0;
160  }
161  const Trk::TrackStates* vtsos = track->trackStateOnSurfaces();
162  if (!vtsos) {
163  return 0;
164  }
165 
166  EGasType gasType;
167  Trk::TrackStates::const_iterator itr = vtsos->begin();
168  Trk::TrackStates::const_iterator itre = vtsos->end();
169  size_t vtos_size = vtsos->size();
170  double correctionFactor = 1.;
171 
172  if (m_toolScenario==kAlgStandard) {
173  std::vector<double> vecToT;
174  vecToT.reserve(vtos_size);
175  double ToTsum = 0;
176 
177  for ( ; itr!=itre ; ++itr) {
178  double length = 0;
179  if ( isGoodHit(ctx,(*itr), useHitsHT, length)) {
180  double ToT_correct = correctToT_corrRZ(ctx,*itr, length);
181  if (m_correctionType == kHitBased){
182  correctionFactor = hitOccupancyCorrection(ctx,*itr);
183  ToT_correct*=correctionFactor;
184  }
185  vecToT.push_back(ToT_correct);
186  }
187  }
188 
189  sort(vecToT.begin(), vecToT.end());
190  size_t nhits = vecToT.size();
191 
192  if (m_divideByL) {
193  nhits-=m_nTrunkateHits;
194  }
195 
196  // Boost speed
197  if (nhits<1) return 0.0;
198 
199  ToTsum = std::accumulate(vecToT.begin(), vecToT.end(), 0);
200  if (m_correctionType == kTrackBased) {
201  correctionFactor=trackOccupancyCorrection(ctx,track, useHitsHT,localOccupancy);
202  } else {
203  correctionFactor=correctNormalization(ctx,nVtx);
204  }
205  ToTsum*=correctionFactor;
206  return ToTsum/nhits;
207  }
208 
209  if(m_toolScenario==kAlgReweight || m_toolScenario==kAlgReweightTrunkOne) {
210  std::vector<double> vecToT_Xe;
211  vecToT_Xe.reserve(vtos_size/2);
212  std::vector<double> vecToT_Ar;
213  vecToT_Ar.reserve(vtos_size/2);
214  std::vector<double> vecToT_Kr;
215 
216  if(m_useTrackPartWithGasType!=kUnset) {
217  ATH_MSG_WARNING(
218  "dEdX_Estimator():: Using m_toolScenario="
219  << m_toolScenario << " scenario m_useTrackPartWithGasType is set to"
220  << m_useTrackPartWithGasType
221  << ", but kUnset is requiered. Check you tool configuration.");
222  }
223 
224  for ( ; itr!=itre ; ++itr) {
225  double length=0;
226  if (isGoodHit(ctx,(*itr), useHitsHT, length)) {
227  gasType=gasTypeInStraw(ctx,*itr);
228  double ToT_correct = correctToT_corrRZ(ctx,*itr, length);
229  if (m_correctionType == kHitBased) {
230  correctionFactor = hitOccupancyCorrection(ctx,*itr);
231  }
232  ToT_correct*=correctionFactor;
233  if(gasType==kXenon) {
234  vecToT_Xe.push_back(ToT_correct);
235  } else if (gasType==kArgon) {
236  vecToT_Ar.push_back(ToT_correct);
237  } else if (gasType==kKrypton) {
238  vecToT_Kr.push_back(ToT_correct);
239  } else {
240  ATH_MSG_ERROR("dEdX_Estimator():: During scenario kAlgReweight variable gasTypeInStraw got value kUnset.");
241  }
242  }
243  }
244 
245  sort(vecToT_Xe.begin(), vecToT_Xe.end());
246  sort(vecToT_Ar.begin(), vecToT_Ar.end());
247  sort(vecToT_Kr.begin(), vecToT_Kr.end());
248 
249  size_t nhitsXe = vecToT_Xe.size();
250  size_t nhitsAr = vecToT_Ar.size();
251  size_t nhitsKr = vecToT_Kr.size();
252 
253  if (m_divideByL) {
254  if(m_toolScenario==kAlgReweight){
255  if (nhitsXe>=m_nTrunkateHits) nhitsXe-=m_nTrunkateHits;
256  if (nhitsAr>=m_nTrunkateHits) nhitsAr-=m_nTrunkateHits;
257  if (nhitsKr>=m_nTrunkateHits) nhitsKr-=m_nTrunkateHits;
258  } else {// kAlgReweightTrunkOne
259  int trunkGas = kUnset;
260  double maxToT = 0.;
261  if (nhitsXe>0 && vecToT_Xe.at(nhitsXe-1)>maxToT) {
262  trunkGas = kXenon;
263  maxToT = vecToT_Xe.at(nhitsXe-1);
264  }
265  if(nhitsAr>0 && vecToT_Ar.at(nhitsAr-1)>maxToT){
266  trunkGas = kArgon;
267  maxToT = vecToT_Ar.at(nhitsAr-1);
268  }
269  if (nhitsKr>0 && vecToT_Kr.at(nhitsKr-1)>maxToT) {
270  trunkGas = kKrypton;
271  }
272 
273  if (trunkGas==kXenon) {
274  nhitsXe-=m_nTrunkateHits;
275  } else if (trunkGas==kArgon) {
276  nhitsAr-=m_nTrunkateHits;
277  } else if (trunkGas==kKrypton) {
278  nhitsKr-=m_nTrunkateHits;
279  }
280  }
281  }
282 
283  size_t nhits = nhitsXe + nhitsAr + nhitsKr;
284  if(nhits<1) return 0.0;
285 
286  double ToTsumXe = 0;
287  double ToTsumAr = 0;
288  double ToTsumKr = 0;
289  for (size_t i = 0; i < nhitsXe;i++) {
290  ToTsumXe+=vecToT_Xe[i];
291  }
292  for (size_t i = 0; i < nhitsAr;i++) {
293  ToTsumAr+=vecToT_Ar[i];
294  }
295  for (size_t i = 0; i < nhitsKr;i++) {
296  ToTsumKr+=vecToT_Kr[i];
297  }
298 
299  ToTsumXe = (nhitsXe>0) ? ToTsumXe/nhitsXe : 0;
300  ToTsumAr = (nhitsAr>0) ? ToTsumAr/nhitsAr : 0;
301  ToTsumKr = (nhitsKr>0) ? ToTsumKr/nhitsKr : 0;
302  double ToTsum = ToTsumXe*nhitsXe + ToTsumAr*nhitsAr + ToTsumKr*nhitsKr;
303 
304  if (m_correctionType == kTrackBased) {
305  correctionFactor =
306  trackOccupancyCorrection(ctx, track, true, localOccupancy);
307  } else {
308  correctionFactor = correctNormalization(ctx, nVtx);
309  }
310  ToTsum *= correctionFactor;
311 
312  return ToTsum/nhits;
313  }
314 
315  ATH_MSG_ERROR("dEdX_Estimator():: m_toolScenario has wrong value "<<m_toolScenario<<"");
316  return 0.;
317 }

◆ dEdx() [3/4]

double ITRT_ToT_dEdx::dEdx ( const Trk::Track *  track,
bool  useHitsHT,
std::optional< float >  localOccupancy = std::nullopt 
) const
inlinevirtualinherited

Definition at line 91 of file ITRT_ToT_dEdx.h.

94 {
95  return dEdx(Gaudi::Hive::currentContext(), track, useHitsHT, localOccupancy);
96 }

◆ dEdx() [4/4]

double ITRT_ToT_dEdx::dEdx
inline

function to calculate sum ToT normalised to number of used hits

Parameters
trackpointer to track
useHitsHTdecide if HT hits should be used in the estimate
Returns
dEdx value

Definition at line 45 of file ITRT_ToT_dEdx.h.

94 {
95  return dEdx(Gaudi::Hive::currentContext(), track, useHitsHT, localOccupancy);
96 }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ fitFuncBarrel_corrRZ()

double TRT_ToT_dEdx::fitFuncBarrel_corrRZ ( const EventContext &  ctx,
EGasType  gas,
double  driftRadius,
double  zPosition,
int  Layer,
int  StrawLayer 
) const
private

function to compute correction factor in barrel region

Parameters
driftradius
zposition
layerindex
strawlayer index
boolvariable to specify whether data or MC correction
Returns
correction

Definition at line 779 of file TRT_ToT_dEdx.cxx.

785 {
786  if (Layer == 0 && StrawLayer < 9) {
787  return fitFuncBarrelShort_corrRZ(ctx,gasType, driftRadius, zPosition, StrawLayer);
788  }
789  return fitFuncBarrelLong_corrRZ(ctx,gasType, driftRadius, zPosition, Layer, StrawLayer);
790 }

◆ fitFuncBarrel_corrRZL()

double TRT_ToT_dEdx::fitFuncBarrel_corrRZL ( const EventContext &  ctx,
EGasType  gasType,
double  driftRadius,
double  zPosition,
int  Layer,
int  StrawLayer 
) const
private

function to compute correction factor in barrel region

Parameters
driftradius
zposition
layerindex
strawlayer index
desiredversion of ToT definition
boolvariable to specify whether data or MC correction
Returns
correction

Definition at line 990 of file TRT_ToT_dEdx.cxx.

996 {
997  /*
998  * T(r,z) = T0(r)+ b(r)*z*z
999  */
1000  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
1001  const TRTDedxcorrection* dEdxCorrection{*readHandle};
1002  if(dEdxCorrection==nullptr) {
1003  ATH_MSG_ERROR(" fitFuncBarrel_corrRZL: Could not find any dEdxCorrection in CondStore. Return zero.");
1004  return 0;
1005  }
1006 
1007  double a,b,c,d,e,f,g;
1008  if (Layer==0 && Strawlayer<9) { // short straws
1009  if (m_isData){
1010  a = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(0)*9+Strawlayer];
1011  b = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(1)*9+Strawlayer];
1012  c = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(2)*9+Strawlayer];
1013  d = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(3)*9+Strawlayer];
1014  e = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(4)*9+Strawlayer];
1015  f = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(5)*9+Strawlayer];
1016  g = dEdxCorrection->paraShortCorrRZDivideByLengthDATA[gasType][(6)*9+Strawlayer];
1017  } else {
1018  a = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(0)*9+Strawlayer];
1019  b = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(1)*9+Strawlayer];
1020  c = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(2)*9+Strawlayer];
1021  d = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(3)*9+Strawlayer];
1022  e = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(4)*9+Strawlayer];
1023  f = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(5)*9+Strawlayer];
1024  g = dEdxCorrection->paraShortCorrRZDivideByLengthMC[gasType][(6)*9+Strawlayer];
1025  }
1026  } else {
1027  if (m_isData) {
1028  a = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(0)*30*3+Layer*30+Strawlayer];
1029  b = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(1)*30*3+Layer*30+Strawlayer];
1030  c = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(2)*30*3+Layer*30+Strawlayer];
1031  d = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(3)*30*3+Layer*30+Strawlayer];
1032  e = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(4)*30*3+Layer*30+Strawlayer];
1033  f = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(5)*30*3+Layer*30+Strawlayer];
1034  g = dEdxCorrection->paraLongCorrRZDivideByLengthDATA[gasType][(6)*30*3+Layer*30+Strawlayer];
1035  } else {
1036  a = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(0)*30*3+Layer*30+Strawlayer];
1037  b = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(1)*30*3+Layer*30+Strawlayer];
1038  c = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(2)*30*3+Layer*30+Strawlayer];
1039  d = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(3)*30*3+Layer*30+Strawlayer];
1040  e = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(4)*30*3+Layer*30+Strawlayer];
1041  f = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(5)*30*3+Layer*30+Strawlayer];
1042  g = dEdxCorrection->paraLongCorrRZDivideByLengthMC[gasType][(6)*30*3+Layer*30+Strawlayer];
1043  }
1044  }
1045  double z = std::abs(zPosition);
1046  double r = std::abs(driftRadius);
1047  double T0neg=a;
1048  double T0pos=b;
1049  double T1 = std::exp(-c*r*r)+d*r;
1050  double slope = e*r+f*r*r+g*r*r*r;
1051  double result;
1052  result = T0neg+T1+slope*z;
1053  if (zPosition>0) result = T0pos+T1+slope*z;
1054 
1055  return result;
1056 }

◆ fitFuncBarrelLong_corrRZ()

double TRT_ToT_dEdx::fitFuncBarrelLong_corrRZ ( const EventContext &  ctx,
EGasType  gasType,
double  driftRadius,
double  zPosition,
int  Layer,
int  StrawLayer 
) const
private

function called by fitFuncBarrel_corrRZ for long straws 

              |z|       /|z| - l  \

T(r,z) = T0(r) + — exp|-------—| v(r) \ s(r) /

Definition at line 809 of file TRT_ToT_dEdx.cxx.

815 {
821  double z = std::abs(zPosition);
822  int sign=1;
823  if(zPosition<0)sign=-1;
824  double l = 704.6;
825  // define set of parameters for negative and positive z positions
826  double T0 = fitFuncPol_corrRZ(ctx,gasType, 0,driftRadius,Layer,StrawLayer,sign,0);
827  double v = fitFuncPol_corrRZ(ctx,gasType, 1,driftRadius,Layer,StrawLayer,sign,0);
828  double s = fitFuncPol_corrRZ(ctx,gasType, 2,driftRadius,Layer,StrawLayer,sign,0);
829  //_in theory_ For IEEE-compatible type double, argument causes exp to overflow if outside [-708.4, 709.8]
830  //however, overflow still seen when argument is 702; so I restrict these to -600, 600
831  const double expArg=(z-l)/s;
832  if (not inRange(expArg, -600.0,600.0)){
833  return expArg>0 ? std::numeric_limits<double>::infinity():0.;
834  }
835  return T0+(z/v)*std::exp(expArg);
836 }

◆ fitFuncBarrelShort_corrRZ()

double TRT_ToT_dEdx::fitFuncBarrelShort_corrRZ ( const EventContext &  ctx,
EGasType  gasType,
double  driftRadius,
double  zPosition,
int  StrawLayer 
) const
private

function called by fitFuncBarrel_corrRZ for short straws

T(r,z) = T0(r)+ b(r)*|z|

Definition at line 839 of file TRT_ToT_dEdx.cxx.

844 {
848  double z = std::abs(zPosition);
849  int sign=1;
850  if(zPosition<0)sign=-1;
851  double T0 = fitFuncPol_corrRZ(ctx,gasType, 0,driftRadius,0,StrawLayer,sign,1);
852  double b = fitFuncPol_corrRZ(ctx,gasType, 1,driftRadius,0,StrawLayer,sign,1);
853  return T0+b*z;
854 }

◆ fitFuncEndcap_corrRZ()

double TRT_ToT_dEdx::fitFuncEndcap_corrRZ ( const EventContext &  ctx,
EGasType  gas,
double  driftRadius,
double  rPosition,
int  Layer,
int  sign 
) const
private

function to compute correction factor in endcap region

Parameters
driftradius
radialposition
strawlayer index
signfor positive or negative side
boolvariable to specify whether data or MC correction
Returns
correction

T(r,R) = T0(r)+ a(r)*R

Definition at line 793 of file TRT_ToT_dEdx.cxx.

799 {
803  double T0 = fitFuncPol_corrRZ(ctx,gasType, 0,driftRadius,Layer,0,sign,2);
804  double a = fitFuncPol_corrRZ(ctx,gasType, 1,driftRadius,Layer,0,sign,2);
805  return T0+a*radialPosition;
806 }

◆ fitFuncEndcap_corrRZL()

double TRT_ToT_dEdx::fitFuncEndcap_corrRZL ( const EventContext &  ctx,
EGasType  gasType,
double  driftRadius,
double  radialPosition,
int  Layer,
int  sign 
) const
private

function to compute correction factor in endcap region

Parameters
driftradius
radialposition
strawlayer index
signfor positive or negative side
boolvariable to specify whether data or MC correction
Returns
correction

Definition at line 939 of file TRT_ToT_dEdx.cxx.

945 {
946  /*
947  * T(r,R) = T0(r)+ a(r)*R
948  */
949 
950  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
951  const TRTDedxcorrection* dEdxCorrection{*readHandle};
952  if(dEdxCorrection==nullptr) {
953  ATH_MSG_ERROR(" fitFuncEndcap_corrRZL: Could not find any dEdxCorrection in CondStore. Return zero.");
954  return 0;
955  }
956 
957  double r = std::abs(driftRadius);
958  double a,b,c,d,e,f,g,h,i;
959  if(sign >0) Layer+=14;
960  if(m_isData){
961  a = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(0)*28+Layer];
962  b = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(1)*28+Layer];
963  c = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(2)*28+Layer];
964  d = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(3)*28+Layer];
965  e = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(4)*28+Layer];
966  f = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(5)*28+Layer];
967  g = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(6)*28+Layer];
968  h = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(7)*28+Layer];
969  i = dEdxCorrection->paraEndCorrRZDivideByLengthDATA[gasType][(8)*28+Layer];
970  }else{
971  a = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(0)*28+Layer];
972  b = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(1)*28+Layer];
973  c = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(2)*28+Layer];
974  d = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(3)*28+Layer];
975  e = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(4)*28+Layer];
976  f = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(5)*28+Layer];
977  g = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(6)*28+Layer];
978  h = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(7)*28+Layer];
979  i = dEdxCorrection->paraEndCorrRZDivideByLengthMC[gasType][(8)*28+Layer];
980  }
981 
982  double T1 = b*r+c*r*r+d*r*r*r+e*r*r*r*r+f*r*r*r*r*r;
983  double slope = g+h*r+i*r*r;
984  double T0 = a;
985 
986  return T0+T1+slope*radialPosition;
987 }

◆ fitFuncPol_corrRZ()

double TRT_ToT_dEdx::fitFuncPol_corrRZ ( const EventContext &  ctx,
EGasType  gasType,
int  parameter,
double  driftRadius,
int  Layer,
int  Strawlayer,
int  sign,
int  set 
) const
private

function called by fitFuncBarrel_corrRZ and fitFuncEndcap_corrRZ

Definition at line 857 of file TRT_ToT_dEdx.cxx.

865 {
866  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
867  const TRTDedxcorrection* dEdxCorrection{*readHandle};
868  if(dEdxCorrection==nullptr)
869  {
870  ATH_MSG_ERROR(" fitFuncPol_corrRZ: Could not find any dEdxCorrection in CondStore. Return zero.");
871  return 0;
872  }
873 
874  double a = 0;
875  double b = 0;
876  double c = 0;
877  double d = 0;
878  double e = 0;
879  double f = 0;
880  double r = driftRadius;
881  int offset = 0;
882  if(m_isData){
883  if(set==0){ // long straws in barrel
884  a = dEdxCorrection->paraLongCorrRZ[gasType][(6*parameter+0)*30*3+Layer*30+Strawlayer+offset];
885  b = dEdxCorrection->paraLongCorrRZ[gasType][(6*parameter+1)*30*3+Layer*30+Strawlayer+offset];
886  c = dEdxCorrection->paraLongCorrRZ[gasType][(6*parameter+2)*30*3+Layer*30+Strawlayer+offset];
887  d = dEdxCorrection->paraLongCorrRZ[gasType][(6*parameter+3)*30*3+Layer*30+Strawlayer+offset];
888  e = dEdxCorrection->paraLongCorrRZ[gasType][(6*parameter+4)*30*3+Layer*30+Strawlayer+offset];
889  f = dEdxCorrection->paraLongCorrRZ[gasType][(6*parameter+5)*30*3+Layer*30+Strawlayer+offset];
890 
891  }else if (set ==1) { // short straws in barrel
892  if(sign > 0) offset+=108;
893  a = dEdxCorrection->paraShortCorrRZ[gasType][(6*parameter+0)*9+Layer+offset];
894  b = dEdxCorrection->paraShortCorrRZ[gasType][(6*parameter+1)*9+Layer+offset];
895  c = dEdxCorrection->paraShortCorrRZ[gasType][(6*parameter+2)*9+Layer+offset];
896  d = dEdxCorrection->paraShortCorrRZ[gasType][(6*parameter+3)*9+Layer+offset];
897  e = dEdxCorrection->paraShortCorrRZ[gasType][(6*parameter+4)*9+Layer+offset];
898  f = dEdxCorrection->paraShortCorrRZ[gasType][(6*parameter+5)*9+Layer+offset];
899  }else{ // straws in endcap
900  if(sign >0) Layer+=14;
901  a = dEdxCorrection->paraEndCorrRZ[gasType][(6*parameter+0)*28+Layer];
902  b = dEdxCorrection->paraEndCorrRZ[gasType][(6*parameter+1)*28+Layer];
903  c = dEdxCorrection->paraEndCorrRZ[gasType][(6*parameter+2)*28+Layer];
904  d = dEdxCorrection->paraEndCorrRZ[gasType][(6*parameter+3)*28+Layer];
905  e = dEdxCorrection->paraEndCorrRZ[gasType][(6*parameter+4)*28+Layer];
906  f = dEdxCorrection->paraEndCorrRZ[gasType][(6*parameter+5)*28+Layer];
907  }
908  }else{
909  if(set==0){ // long straws in barrel
910  if(sign > 0) offset=1620;
911  a = dEdxCorrection->paraLongCorrRZMC[gasType][(6*parameter+0)*30*3+Layer*30+Strawlayer+offset];
912  b = dEdxCorrection->paraLongCorrRZMC[gasType][(6*parameter+1)*30*3+Layer*30+Strawlayer+offset];
913  c = dEdxCorrection->paraLongCorrRZMC[gasType][(6*parameter+2)*30*3+Layer*30+Strawlayer+offset];
914  d = dEdxCorrection->paraLongCorrRZMC[gasType][(6*parameter+3)*30*3+Layer*30+Strawlayer+offset];
915  e = dEdxCorrection->paraLongCorrRZMC[gasType][(6*parameter+4)*30*3+Layer*30+Strawlayer+offset];
916  f = dEdxCorrection->paraLongCorrRZMC[gasType][(6*parameter+5)*30*3+Layer*30+Strawlayer+offset];
917  }else if (set ==1) { // short straws in barrel
918  if(sign > 0) offset+=108;
919  a = dEdxCorrection->paraShortCorrRZMC[gasType][(6*parameter+0)*9+Layer+offset];
920  b = dEdxCorrection->paraShortCorrRZMC[gasType][(6*parameter+1)*9+Layer+offset];
921  c = dEdxCorrection->paraShortCorrRZMC[gasType][(6*parameter+2)*9+Layer+offset];
922  d = dEdxCorrection->paraShortCorrRZMC[gasType][(6*parameter+3)*9+Layer+offset];
923  e = dEdxCorrection->paraShortCorrRZMC[gasType][(6*parameter+4)*9+Layer+offset];
924  f = dEdxCorrection->paraShortCorrRZMC[gasType][(6*parameter+5)*9+Layer+offset];
925  }else{ // straws in endcap
926  if(sign >0) Layer+=14;
927  a = dEdxCorrection->paraEndCorrRZMC[gasType][(6*parameter+0)*28+Layer];
928  b = dEdxCorrection->paraEndCorrRZMC[gasType][(6*parameter+1)*28+Layer];
929  c = dEdxCorrection->paraEndCorrRZMC[gasType][(6*parameter+2)*28+Layer];
930  d = dEdxCorrection->paraEndCorrRZMC[gasType][(6*parameter+3)*28+Layer];
931  e = dEdxCorrection->paraEndCorrRZMC[gasType][(6*parameter+4)*28+Layer];
932  f = dEdxCorrection->paraEndCorrRZMC[gasType][(6*parameter+5)*28+Layer];
933  }
934  }
935  return a+b*r+c*r*r+d*r*r*r+e*r*r*r*r+f*r*r*r*r*r;
936 }

◆ gasTypeInStraw() [1/2]

TRT_ToT_dEdx::EGasType TRT_ToT_dEdx::gasTypeInStraw ( const EventContext &  ctx,
const InDet::TRT_DriftCircleOnTrack *  driftcircle 
) const
protected

Definition at line 649 of file TRT_ToT_dEdx.cxx.

652 {
653  Identifier DCid = driftcircle->identify();
654 
655  // getStatusHT returns enum {Undefined, Dead, Good, Xenon, Argon, Krypton, EmulatedArgon, EmulatedKrypton}.
656  // Our representation of 'GasType' is 0:Xenon, 1:Argon, 2:Krypton
657  EGasType GasType=kUnset; // kUnset is default
658  if (!m_TRTStrawSummaryTool.empty()) {
659  int stat = m_TRTStrawSummaryTool->getStatusHT(DCid,ctx);
660  if ( stat==2 || stat==3 ) { GasType = kXenon; } // Xe
661  else if ( stat==1 || stat==4 ) { GasType = kArgon; } // Ar
662  else if ( stat==5 ) { GasType = kKrypton; } // Kr
663  else if ( stat==6 ) { GasType = kArgon; } // Emulated Ar
664  else if ( stat==7 ) { GasType = kKrypton;
665  } // Emulated Kr
666  else {
667  ATH_MSG_FATAL(
668  "getStatusHT = "
669  << stat
670  << ", must be 'Good(2)||Xenon(3)' or 'Dead(1)||Argon(4)' or "
671  "'Krypton(5)' or 'EmulatedArgon(6)' or 'EmulatedKr(7)'!");
672  throw std::exception();
673  }
674  }
675 
676  return GasType;
677 }

◆ gasTypeInStraw() [2/2]

TRT_ToT_dEdx::EGasType TRT_ToT_dEdx::gasTypeInStraw ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  itr 
) const
protected

return gas type for that hit

Parameters
trackon surface object
Returns
gasType

Definition at line 622 of file TRT_ToT_dEdx.cxx.

624 {
625  const Trk::MeasurementBase* trkM = itr->measurementOnTrack();
626  if (!trkM) {
627  return kUnset;
628  }
629 
630  // Check if this is RIO on track
631  //annd if yes check if is TRT Drift Circle
632  //then set the ptr
633  const InDet::TRT_DriftCircleOnTrack* driftcircle = nullptr;
634  if (trkM->type(Trk::MeasurementBaseType::RIO_OnTrack)) {
635  const Trk::RIO_OnTrack* tmpRio = static_cast<const Trk::RIO_OnTrack*>(trkM);
636  if (tmpRio->rioType(Trk::RIO_OnTrackType::TRT_DriftCircle)) {
637  driftcircle = static_cast<const InDet::TRT_DriftCircleOnTrack*>(tmpRio);
638  }
639  }
640 
641  if (!driftcircle) {
642  return kUnset;
643  }
644 
645  return gasTypeInStraw(ctx,driftcircle);
646 }

◆ getAlgorithm()

int TRT_ToT_dEdx::getAlgorithm ( ) const
inline

Definition at line 383 of file TRT_ToT_dEdx.h.

383 { return m_toolScenario; }

◆ getCorrectionType()

bool TRT_ToT_dEdx::getCorrectionType ( ) const
inline

Definition at line 368 of file TRT_ToT_dEdx.h.

368 { return m_correctionType; }

◆ getGasTypeFordEdXCalculation()

int TRT_ToT_dEdx::getGasTypeFordEdXCalculation ( ) const
inline

Definition at line 380 of file TRT_ToT_dEdx.h.

380 { return m_useTrackPartWithGasType; }

◆ getMaxRtrack()

float TRT_ToT_dEdx::getMaxRtrack ( ) const
inline

Definition at line 374 of file TRT_ToT_dEdx.h.

374 { return m_trackConfig_maxRtrack; }

◆ getMinRtrack()

float TRT_ToT_dEdx::getMinRtrack ( ) const
inline

Definition at line 371 of file TRT_ToT_dEdx.h.

371 { return m_trackConfig_minRtrack; }

◆ getProb() [1/2]

double TRT_ToT_dEdx::getProb ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  itr,
const double  dEdx_obs,
const double  pTrk,
Trk::ParticleHypothesis  hypothesis,
int  nUsedHits 
) const
protected

function to calculate likelihood from prediction and resolution

Parameters
observeddEdx
trackparameter
particlehypothesis
numberof used hits
Returns
brobability value between 0 and 1

Definition at line 399 of file TRT_ToT_dEdx.cxx.

405 {
406  EGasType gasType = gasTypeInStraw(ctx, itr);
407  return getProb(ctx,gasType, dEdx_obs, pTrk, hypothesis, nUsedHits);
408 }

◆ getProb() [2/2]

double TRT_ToT_dEdx::getProb ( const EventContext &  ctx,
EGasType  gasType,
const double  dEdx_obs,
const double  pTrk,
Trk::ParticleHypothesis  hypothesis,
int  nUsedHits 
) const
protected

Definition at line 411 of file TRT_ToT_dEdx.cxx.

417 {
418  ATH_MSG_DEBUG("getProb():: gasTypeInStraw = "<<gasType<<"");
419 
420  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
421  const TRTDedxcorrection* dEdxCorrection{*readHandle};
422  if (dEdxCorrection==nullptr) {
423  ATH_MSG_ERROR(" getProb: Could not find any dEdxCorrection in CondStore. Return zero.");
424  return 0;
425  }
426 
427  if(gasType==kUnset) {
428  ATH_MSG_DEBUG("getProb():: gasTypeInStraw set kUnset that is not allowed! Use gasTypeInStraw(*itr) to get gas type info for that hit first!");
429  ATH_MSG_DEBUG("getProb():: Now gasTypeInStraw sets to kXenon.");
430  gasType = kXenon;
431  }
432 
433  double dEdx_pred = predictdEdx(ctx,gasType, pTrk, hypothesis);
434  if (dEdx_pred==0) return 0.0;
435  if (hypothesis==Trk::electron) {
436  // correction for pTrk in [MeV]
437  double factor = 1;
438  double correct = 1+factor*(0.045*log10(pTrk)+0.885-1);
439  dEdx_pred= dEdx_pred/correct;
440  }
441 
442  double Resolution =
443  dEdxCorrection->resolution[gasType][0] +
444  dEdxCorrection->resolution[gasType][1] * (nUsedHits + 0.5) +
445  dEdxCorrection->resolution[gasType][2] * (nUsedHits + 0.5) *
446  (nUsedHits + 0.5) +
447  dEdxCorrection->resolution[gasType][3] * (nUsedHits + 0.5) *
448  (nUsedHits + 0.5) * (nUsedHits + 0.5);
449  if (hypothesis == Trk::electron) {
450  Resolution =
451  dEdxCorrection->resolutionElectron[gasType][0] +
452  dEdxCorrection->resolutionElectron[gasType][1] * (nUsedHits + 0.5) +
453  dEdxCorrection->resolutionElectron[gasType][2] * (nUsedHits + 0.5) *
454  (nUsedHits + 0.5) +
455  dEdxCorrection->resolutionElectron[gasType][3] * (nUsedHits + 0.5) *
456  (nUsedHits + 0.5) * (nUsedHits + 0.5);
457  }
458 
459  double prob = std::exp( -0.5 * ( ( ( dEdx_obs - dEdx_pred ) / (Resolution*dEdx_pred) ) *
460  ( ( dEdx_obs - dEdx_pred ) / (Resolution*dEdx_pred) ) )) ;
461 
462  ATH_MSG_DEBUG("getProb():: return "<<prob<<"");
463  return prob;
464 }

◆ getStatusCorrection()

bool TRT_ToT_dEdx::getStatusCorrection ( ) const
inline

Definition at line 365 of file TRT_ToT_dEdx.h.

365 { return m_corrected; }

◆ getStatusUseZeroRHitCut()

bool TRT_ToT_dEdx::getStatusUseZeroRHitCut ( ) const
inline

Definition at line 377 of file TRT_ToT_dEdx.h.

377 { return m_useZeroRHitCut; }

◆ getTest() [1/4]

double ITRT_ToT_dEdx::getTest ( const double  dEdx_obs,
const double  pTrk,
Trk::ParticleHypothesis  hypothesis,
Trk::ParticleHypothesis  antihypothesis,
int  nUsedHits 
) const
inlinevirtualinherited

Definition at line 105 of file ITRT_ToT_dEdx.h.

110 {
111  return getTest(Gaudi::Hive::currentContext(),
112  dEdx_obs,
113  pTrk,
114  hypothesis,
115  antihypothesis,
116  nUsedHits);
117 }

◆ getTest() [2/4]

double ITRT_ToT_dEdx::getTest
inline

Definition at line 81 of file ITRT_ToT_dEdx.h.

110 {
111  return getTest(Gaudi::Hive::currentContext(),
112  dEdx_obs,
113  pTrk,
114  hypothesis,
115  antihypothesis,
116  nUsedHits);
117 }

◆ getTest() [3/4]

double TRT_ToT_dEdx::getTest ( const EventContext &  ctx,
const double  dEdx_obs,
const double  pTrk,
Trk::ParticleHypothesis  hypothesis,
Trk::ParticleHypothesis  antihypothesis,
int  nUsedHits 
) const
finaloverridevirtual

function to calculate likelihood ratio test

Parameters
observeddEdx
trackparameter
particlehypothesis
antihypothesis
numberof used hits
Returns
test value between 0 and 1

Implements ITRT_ToT_dEdx.

Definition at line 467 of file TRT_ToT_dEdx.cxx.

473 {
474  ATH_MSG_DEBUG("getTest()");
475 
476  EGasType gasType = kUnset;
477  if ( dEdx_obs<=0. || pTrk<=0. || nUsedHits<=0 ) return 0.5;
478 
479  double Pone = getProb(ctx,gasType, dEdx_obs,pTrk,hypothesis,nUsedHits);
480  double Ptwo = getProb(ctx,gasType, dEdx_obs,pTrk,antihypothesis,nUsedHits);
481  if ((Pone+Ptwo) != 0) {
482  ATH_MSG_DEBUG("getTest():: return "<<Pone/(Pone+Ptwo)<<"");
483  return Pone/(Pone+Ptwo);
484  } else {
485  return 0.5;
486  }
487 }

◆ getTest() [4/4]

virtual double ITRT_ToT_dEdx::getTest

function to calculate likelihood ratio test

Parameters
Eventcontext ctx
observeddEdx
trackparameter
particlehypothesis
antihypothesis
numberof used hits
Returns
test value between 0 and 1

◆ hitOccupancyCorrection()

double TRT_ToT_dEdx::hitOccupancyCorrection ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  itr 
) const
private

Definition at line 1059 of file TRT_ToT_dEdx.cxx.

1061 {
1062  SG::ReadCondHandle<TRTDedxcorrection> readHandle(m_ReadKey,ctx);
1063  const TRTDedxcorrection* dEdxCorrection{*readHandle};
1064 
1065  const Trk::MeasurementBase* trkM = itr->measurementOnTrack();
1066 
1067  // Check if this is RIO on track
1068  // and if yes check if is TRT Drift Circle
1069  // then set the ptr
1070  const InDet::TRT_DriftCircleOnTrack* driftcircle = nullptr;
1071  if (trkM && trkM->type(Trk::MeasurementBaseType::RIO_OnTrack)) {
1072  const Trk::RIO_OnTrack* tmpRio = static_cast<const Trk::RIO_OnTrack*>(trkM);
1073  if (tmpRio->rioType(Trk::RIO_OnTrackType::TRT_DriftCircle)) {
1074  driftcircle = static_cast<const InDet::TRT_DriftCircleOnTrack*>(tmpRio);
1075  }
1076  }
1077 
1078  if (!driftcircle) {
1079  return 0.;
1080  }
1081 
1082  const Trk::TrackParameters* trkP = itr->trackParameters();
1083  Identifier DCId = driftcircle->identify();
1084  int isHT = driftcircle->highLevel();
1085  int isShared=0;
1086 
1087  const Trk::RIO_OnTrack* hit_trt = nullptr;
1088  if (trkM && trkM->type(Trk::MeasurementBaseType::RIO_OnTrack)) {
1089  hit_trt = static_cast<const Trk::RIO_OnTrack*>(trkM);
1090  }
1091 
1092  if (hit_trt) {
1093  if ( m_assoTool->isShared(*(hit_trt->prepRawData())) ) isShared=1;
1094  }
1095  int layer = m_trtId->layer_or_wheel(DCId);
1096  int phimodule = m_trtId->phi_module(DCId);
1097  int HitPart = m_trtId->barrel_ec(DCId);
1098  double Trt_HitTheta = trkP->parameters()[Trk::theta];
1099  double trackEta = -log(tan(Trt_HitTheta/2.0));
1100 
1101  double localOccupancy = m_localOccTool->LocalOccupancy(ctx, trackEta, phimodule);
1102  double ToTmip = 1;
1103  double valToT = 1.;
1104 
1105  double p0=0., p1=0., p2=0., p0_flat=0.;
1106 
1107  //the calibration array is structured as follows (hence the non intuitive numbers)
1108  //the first 36 parameters are for barrel, the last 168 for the endcap
1109  //each of these subarrays is divided into smaller subarrays of length 12 (barrel has 3 and endcap 14 layers)
1110  int nBarrelLayers = 3, nEndcapLayers = 14;
1111  //each layer has 3 parameters (function 2nd order), so a subarray for each layer has 3 parameters
1112  int nParametersPerLayer = 3;
1113  //this subarray for every layer exits for HT/LT hits, so 6 parameters
1114  int nHTConfigurations = 2;
1115  //this subarray exists for shared/non-shared hits, so 12 parameters
1116  int nSharedConfigurations = 2;
1117  int num =
1118  layer * nParametersPerLayer +
1119  isHT * ((abs(HitPart) - 1) * (nEndcapLayers - nBarrelLayers) *
1120  nParametersPerLayer +
1121  nBarrelLayers * nParametersPerLayer) +
1122  isShared * ((abs(HitPart) - 1) * (nEndcapLayers - nBarrelLayers) *
1123  nParametersPerLayer * nHTConfigurations +
1124  nBarrelLayers * nParametersPerLayer * nHTConfigurations) +
1125  (abs(HitPart) - 1) * nParametersPerLayer * nBarrelLayers *
1126  nHTConfigurations * nSharedConfigurations;
1127  //number for that given hit for non-shared conditions
1128  int num_flat = layer * 3 +
1129  isHT * ((abs(HitPart) - 1) * (nEndcapLayers - nBarrelLayers) *
1130  nParametersPerLayer +
1131  nBarrelLayers * nParametersPerLayer) +
1132  (abs(HitPart) - 1) * nParametersPerLayer * nBarrelLayers *
1133  nHTConfigurations * nSharedConfigurations;
1134 
1135  p0 = dEdxCorrection->hitOccPar[num];
1136  p1 = dEdxCorrection->hitOccPar[num+1];
1137  p2 = dEdxCorrection->hitOccPar[num+2];
1138  p0_flat = dEdxCorrection->hitOccPar[num_flat];
1139 
1140  //fitting function is a polynomial 2nd order f(x)=a+b*x+x^2
1141  //Hence the tot value is divided by the value of the function
1142  //multiplied to the non-shared intercept
1143  valToT = p0_flat/(p0+p1*localOccupancy+p2*localOccupancy*localOccupancy);
1144 
1145  return ToTmip*valToT;
1146 }

◆ initialize()

StatusCode TRT_ToT_dEdx::initialize ( )
overridevirtual

AlgTool initailize method.

Definition at line 47 of file TRT_ToT_dEdx.cxx.

48 {
49 
50  MsgStream log(msgSvc(), name());
51 
52  // retrieve TRT-ID helper
53 
54  StatusCode sc = detStore()->retrieve(m_trtId, "TRT_ID");
55  if (sc.isFailure()){
56  ATH_MSG_ERROR ( "Could not get TRT_ID helper !" );
57  return StatusCode::FAILURE;
58  }
59 
60  // Initialize ReadHandleKey and ReadCondHandleKey
61  ATH_CHECK(m_rdhkEvtInfo.initialize());
62  ATH_CHECK(m_ReadKey.initialize());
63  ATH_CHECK(m_assoTool.retrieve());
64  ATH_CHECK( m_localOccTool.retrieve() );
65  ATH_CHECK(m_TRTStrawSummaryTool.retrieve());
66 
67  if (m_useTrackPartWithGasType > EGasType::kUnset ||
68  m_useTrackPartWithGasType < EGasType::kXenon) {
69  ATH_MSG_ERROR("Property TRT_dEdx_useTrackPartWithGasType has an invalid "
70  << "value: " << m_useTrackPartWithGasType);
71  return StatusCode::FAILURE;
72  }
73 
74  return StatusCode::SUCCESS;
75 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& ITRT_ToT_dEdx::interfaceID ( )
inlinestaticinherited

AlgTool interface methods.

Definition at line 30 of file ITRT_ToT_dEdx.h.

30 { return IID_ITRT_ToT_dEdx; };

◆ isGoodHit()

bool TRT_ToT_dEdx::isGoodHit ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  trackState,
bool  useHitsHT,
double &  length 
) const
protected

function to define what is a good hit to be used for dEdx calculation cuts on track level can be made latekAlgStandardr by the user.

Also returns the length in the straw.

Parameters
trackStatemeasurement on track
useHitsHTdecide if HT hits should be used in the estimate
lengthlength in straw
Returns
decision

Definition at line 78 of file TRT_ToT_dEdx.cxx.

82 {
83  const Trk::MeasurementBase* trkM = trackState->measurementOnTrack();
84  if (!trkM) {
85  return false;
86  }
87 
88  // Check if this is RIO on track
89  // and if yes check if is TRT Drift Circle
90  // then set the ptr
91  const InDet::TRT_DriftCircleOnTrack* driftcircle = nullptr;
92  if (trkM->type(Trk::MeasurementBaseType::RIO_OnTrack)) {
93  const Trk::RIO_OnTrack* tmpRio = static_cast<const Trk::RIO_OnTrack*>(trkM);
94  if (tmpRio->rioType(Trk::RIO_OnTrackType::TRT_DriftCircle)) {
95  driftcircle = static_cast<const InDet::TRT_DriftCircleOnTrack*>(tmpRio);
96  }
97  }
98 
99  if (!driftcircle) {
100  return false;
101  }
102 
103  const Trk::TrackParameters* trkP = trackState->trackParameters();
104  if(trkP==nullptr)return false;
105 
106  double Trt_Rtrack = std::abs(trkP->parameters()[Trk::locR]);
107  double Trt_RHit = std::abs(driftcircle->localParameters()[Trk::driftRadius]);
108  double error = std::sqrt(driftcircle->localCovariance()(Trk::driftRadius,Trk::driftRadius));
109 
110  if (trackState->type(Trk::TrackStateOnSurface::Outlier)) return false; //Outliers
111  if (m_useZeroRHitCut && Trt_RHit==0 && error>1.) return false; //Select precision hits only
112  if ((Trt_Rtrack >= m_trackConfig_maxRtrack) || (Trt_Rtrack <= m_trackConfig_minRtrack)) return false; // drift radius close to wire or wall
113 
114  length = calculateTrackLengthInStraw(trackState, m_trtId);
115 
116  if (m_divideByL and length < 1.7) return false; // Length in the straw
117 
118  if (!useHitsHT) {
119  int TrtHl = driftcircle->highLevel();
120  if (TrtHl==1) return false;
121  }
122 
123  if (m_useTrackPartWithGasType != kUnset) { // don't preselect hits
124  if(m_useTrackPartWithGasType != gasTypeInStraw(ctx,trackState)) return false;
125  }
126 
127  if (driftcircle->prepRawData()->timeOverThreshold()==0.) return false; // If ToT for this hit equal 0, skip it.
128 
129  return true;
130 }

◆ mass()

double TRT_ToT_dEdx::mass ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  itr,
const double  pTrk,
double  dEdx 
) const
protected

function to extract most likely mass in bg [0,3]

Parameters
trackmomentum
measureddEdx
Returns
mass
Todo:
make failsafe

Definition at line 570 of file TRT_ToT_dEdx.cxx.

574 {
575  EGasType gasType = gasTypeInStraw(ctx,itr);
576 
577  ATH_MSG_DEBUG("mass(): gasTypeInStraw = "<<gasType<<"");
578 
579  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
580  const TRTDedxcorrection* dEdxCorrection{*readHandle};
581  if(dEdxCorrection==nullptr)
582  {
583  ATH_MSG_ERROR(" mass: Could not find any dEdxCorrection in CondStore. Return zero.");
584  return 0;
585  }
586 
587  if(gasType==kUnset)
588  {
589  ATH_MSG_WARNING("mass():: gasTypeInStraw set kUnset that is not allowed! Use gasTypeInStraw(*itr) to get gas type info for that hit first!");
590  ATH_MSG_WARNING("mass():: Now gasTypeInStraw sets to kXenon.");
591  gasType = kXenon;
592  }
593  // only for testing purposes!!!!
594  // note that dE/dx has to be corrected on track level first before value can be transferred to mass
595  // this has to be tuned on user level
597  static const double bg_min = 0.001;
598  static const double bg_max = 3; // maximal allowed bg
599 
600  static const std::string blumRolandiFunction =
601  "( [0]/sqrt( (x*x/([5]*[5]))/(1.+(x*x/([5]*[5]))) )^[3] ) * ([1] - sqrt( "
602  "(x*x/([5]*[5]))/(1.+(x*x/([5]*[5]))) )^[3] - log([2]+1./((x/[5])^[4]) ) )";
603 
604  TF1 blumRolandi( "BR", blumRolandiFunction.c_str(), 0.7, 100000);
605 
606  blumRolandi.SetParameters(dEdxCorrection->paraDedxP1[gasType],
607  dEdxCorrection->paraDedxP2[gasType],
608  dEdxCorrection->paraDedxP3[gasType],
609  dEdxCorrection->paraDedxP4[gasType],
610  dEdxCorrection->paraDedxP5[gasType],
611  1.);
612  //blumRolandi.SetParameters(&dEdxCorrection->para_dEdx_BB);
613  double betaGamma = blumRolandi.GetX(dEdx, bg_min, bg_max);
614 
615  ATH_MSG_DEBUG("mass():: return "<<pTrk/betaGamma<<"");
616 
617  return pTrk/betaGamma;
618 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ predictdEdx() [1/2]

double TRT_ToT_dEdx::predictdEdx ( const EventContext &  ctx,
const Trk::TrackStateOnSurface *  itr,
const double  pTrk,
Trk::ParticleHypothesis  hypothesis 
) const
protected

function to calculate expectation value for dEdx using BB fit

Parameters
trackmomentum
hypothesis
Returns
dEdx_pred

Definition at line 490 of file TRT_ToT_dEdx.cxx.

494 {
495  EGasType gasType = gasTypeInStraw(ctx, itr);
496  return predictdEdx(ctx,gasType, pTrk, hypothesis);
497 }

◆ predictdEdx() [2/2]

double TRT_ToT_dEdx::predictdEdx ( const EventContext &  ctx,
EGasType  gasType,
const double  pTrk,
Trk::ParticleHypothesis  hypothesis 
) const
protected
Todo:
why is it possible that a 20 MeV particle reaches this point? (see Savannah bug 94644)

Definition at line 500 of file TRT_ToT_dEdx.cxx.

504 {
505  ATH_MSG_DEBUG("predictdEdx(): gasTypeInStraw = " << gasType << "");
506 
507  SG::ReadCondHandle<TRTDedxcorrection> readHandle{ m_ReadKey,ctx };
508  const TRTDedxcorrection* dEdxCorrection{ *readHandle };
509  if (dEdxCorrection == nullptr) {
510  ATH_MSG_ERROR(" predictdEdx: Could not find any dEdxCorrection in "
511  "CondStore. Return zero.");
512  return 0;
513  }
514 
515  if (gasType == kUnset) {
516  ATH_MSG_DEBUG(
517  "predictdEdx():: gasTypeInStraw set kUnset that is not allowed! Use "
518  "gasTypeInStraw(*itr) to get gas type info for that hit first!");
519  ATH_MSG_DEBUG("predictdEdx():: Now gasTypeInStraw sets to kXenon.");
520  gasType = kXenon;
521  }
522 
523  double mass = Trk::ParticleMasses::mass[hypothesis];
524 
525  double betaGamma = pTrk/mass;
527  // low momentum particle can create floating point error
528  // do we need the check in the log parameter in addition? will create CPU increase
529  // do we want to throw an assertion here?
530  if (pTrk < 100)
531  return 0;
532  if (m_divideByL) {
533  if (dEdxCorrection->paraDivideByLengthDedxP3[gasType] +
534  1. / (std::pow(betaGamma,
535  dEdxCorrection->paraDivideByLengthDedxP5[gasType])) <=
536  0)
537  return 0;
538  return dEdxCorrection->paraDivideByLengthDedxP1[gasType] /
539  std::pow(std::sqrt((betaGamma * betaGamma) /
540  (1. + (betaGamma * betaGamma))),
541  dEdxCorrection->paraDivideByLengthDedxP4[gasType]) *
542  (dEdxCorrection->paraDivideByLengthDedxP2[gasType] -
543  std::pow(std::sqrt((betaGamma * betaGamma) /
544  (1. + (betaGamma * betaGamma))),
545  dEdxCorrection->paraDivideByLengthDedxP4[gasType]) -
546  log(dEdxCorrection->paraDivideByLengthDedxP3[gasType] +
547  1. / (std::pow(
548  betaGamma,
549  dEdxCorrection->paraDivideByLengthDedxP5[gasType]))));
550  }
551  if (dEdxCorrection->paraDedxP3[gasType] +
552  1. / (std::pow(betaGamma, dEdxCorrection->paraDedxP5[gasType])) <=
553  0)
554  return 0;
555  return dEdxCorrection->paraDedxP1[gasType] /
556  std::pow(
557  std::sqrt((betaGamma * betaGamma) / (1. + (betaGamma * betaGamma))),
558  dEdxCorrection->paraDedxP4[gasType]) *
559  (dEdxCorrection->paraDedxP2[gasType] -
560  std::pow(
561  std::sqrt((betaGamma * betaGamma) / (1. + (betaGamma * betaGamma))),
562  dEdxCorrection->paraDedxP4[gasType]) -
563  log(dEdxCorrection->paraDedxP3[gasType] +
564  1. / (std::pow(betaGamma, dEdxCorrection->paraDedxP5[gasType]))));
565 
566  //return 0;
567 }

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ setAlgorithm()

void TRT_ToT_dEdx::setAlgorithm ( EstCalc  alg)
inline

Definition at line 382 of file TRT_ToT_dEdx.h.

382 { m_toolScenario=alg; }

◆ setCorrectionType()

void TRT_ToT_dEdx::setCorrectionType ( EOccupancyCorrection  value)
inline

Definition at line 367 of file TRT_ToT_dEdx.h.

367 { m_correctionType = value; }

◆ setGasTypeFordEdXCalculation()

void TRT_ToT_dEdx::setGasTypeFordEdXCalculation ( EGasType  gasType)
inline

Definition at line 379 of file TRT_ToT_dEdx.h.

379 { m_useTrackPartWithGasType = gasType; }

◆ setMaxRtrack()

void TRT_ToT_dEdx::setMaxRtrack ( float  maxRtrack)
inline

Definition at line 373 of file TRT_ToT_dEdx.h.

373 { m_trackConfig_maxRtrack=maxRtrack;}

◆ setMinRtrack()

void TRT_ToT_dEdx::setMinRtrack ( float  minRtrack)
inline

Definition at line 370 of file TRT_ToT_dEdx.h.

370 { m_trackConfig_minRtrack=minRtrack;}

◆ setStatusCorrection()

void TRT_ToT_dEdx::setStatusCorrection ( bool  value)
inline

Definition at line 364 of file TRT_ToT_dEdx.h.

364 { m_corrected = value; }

◆ setStatusUseZeroRHitCut()

void TRT_ToT_dEdx::setStatusUseZeroRHitCut ( bool  value)
inline

Definition at line 376 of file TRT_ToT_dEdx.h.

376 { m_useZeroRHitCut = value; }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ trackOccupancyCorrection()

double TRT_ToT_dEdx::trackOccupancyCorrection ( const EventContext &  ctx,
const Trk::Track *  track,
bool  useHThits,
std::optional< float >  localOccupancy 
) const
private

Definition at line 1149 of file TRT_ToT_dEdx.cxx.

1153 {
1154  SG::ReadCondHandle<TRTDedxcorrection> readHandle{m_ReadKey,ctx};
1155  const TRTDedxcorrection* dEdxCorrection{*readHandle};
1156 
1157  double corr=-999.;
1158 
1159  double trackOcc = localOccupancy != std::nullopt
1160  ? *localOccupancy
1161  : m_localOccTool->LocalOccupancy(ctx, *track);
1162  const Trk::TrackParameters* perigee = track->perigeeParameters();
1163  const AmgVector(Trk::TrackParameters::dim)& parameterVector = perigee->parameters();
1164  double theta = parameterVector[Trk::theta];
1165  double trackEta = -log(tan(theta/2.0));
1166 
1167  // the correction constants were determined in 100 bins of 0.04 in the eta range between -2 and 2
1168  int index = int(25*(trackEta+2.)); // determine the bin of the corrections
1169  if (index < 0) index = 0; // lower bound corresponding to eta = -2
1170  else if (index > 99) index = 99; // upper bound corresponding to eta = +2
1171 
1172  //Function of the from f(x)=a+b*x+c*x^2 was used as a fitting function, separately for tracks with and without excluding HT hits
1173  if (useHitsHT) {
1174  corr=dEdxCorrection->trackOccPar0[index]+dEdxCorrection->trackOccPar1[index]*trackOcc+dEdxCorrection->trackOccPar2[index]*pow(trackOcc,2);
1175  } else {
1176  corr=dEdxCorrection->trackOccPar0NoHt[index]+dEdxCorrection->trackOccPar1NoHt[index]*trackOcc+dEdxCorrection->trackOccPar2NoHt[index]*pow(trackOcc,2);
1177  }
1178 
1179  if (corr != 0) {
1180  return 1./corr;
1181  }
1182  return 0.;
1183 }

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

◆ usedHits() [1/4]

virtual double ITRT_ToT_dEdx::usedHits

function to calculate number of used hits

Parameters
Eventcontext ctx
trackpointer
useHitsHTdecide if HT hits should be used in the estimate
Returns
nHits

◆ usedHits() [2/4]

double TRT_ToT_dEdx::usedHits ( const EventContext &  ctx,
const Trk::Track *  track,
bool  useHThits = true 
) const
finaloverridevirtual

function to calculate number of used hits

Parameters
trackpointer
useHitsHTdecide if HT hits should be used in the estimate
Returns
nHits

Implements ITRT_ToT_dEdx.

Definition at line 320 of file TRT_ToT_dEdx.cxx.

323 {
324  ATH_MSG_DEBUG("usedHits()");
325  EGasType gasType = kUnset;
326 
327  if (!track) {
328  return 0;
329  }
330  const Trk::TrackStates* vtsos = track->trackStateOnSurfaces();
331  if (!vtsos) {
332  return 0;
333  }
334 
335  Trk::TrackStates::const_iterator itr = vtsos->begin();
336  Trk::TrackStates::const_iterator itre = vtsos->end();
337 
338  if (m_toolScenario==kAlgStandard) {
339  int nhits =0;
340 
341  for ( ; itr!=itre ; ++itr) {
342  double length=0;
343  if (isGoodHit(ctx,(*itr), useHitsHT, length)) {
344  nhits++;
345  }
346  }
347  if (m_divideByL) nhits -= m_nTrunkateHits;
348  return nhits;
349  } else if (m_toolScenario==kAlgReweight || m_toolScenario==kAlgReweightTrunkOne) {
350  int nhits = 0;
351  int nhitsXe = 0;
352  int nhitsAr = 0;
353  int nhitsKr = 0;
354 
355  if(m_useTrackPartWithGasType!=kUnset) {
356  ATH_MSG_WARNING(
357  "usedHits_Estimator():: Using m_toolScenario="
358  << m_toolScenario << " scenario m_useTrackPartWithGasType is set to "
359  << m_useTrackPartWithGasType
360  << ", but kUnset is required. Check you tool configuration.");
361  }
362 
363  for ( ; itr!=itre ; ++itr) {
364  double length=0;
365  if ( isGoodHit(ctx,(*itr), useHitsHT, length)) {
366  gasType=gasTypeInStraw(ctx,*itr);
367  if (gasType==kXenon) {
368  nhitsXe++;
369  } else if (gasType==kArgon) {
370  nhitsAr++;
371  } else if (gasType==kKrypton) {
372  nhitsKr++;
373  } else {
374  ATH_MSG_ERROR("usedHits_Estimator():: During scenario kAlgReweight variable gasTypeInStraw got value kUnset.");
375  }
376  }
377  }
378 
379  if (m_divideByL) {
380  if(m_toolScenario==kAlgReweight){
381  if(nhitsXe>0) nhitsXe -= m_nTrunkateHits;
382  if(nhitsAr>0) nhitsAr -= m_nTrunkateHits;
383  if(nhitsKr>0) nhitsKr -= m_nTrunkateHits;
384  } else { // kAlgReweightTrunkOne
385  if(nhitsXe>0 || nhitsAr>0 || nhitsKr>0)
386  nhitsXe -= m_nTrunkateHits;
387  }
388  }
389 
390  nhits = nhitsXe + nhitsAr + nhitsKr;
391  return nhits;
392  }
393 
394  ATH_MSG_ERROR("usedHits_Estimator():: m_toolScenario has wrong value "<<m_toolScenario<<"");
395  return 0;
396 }

◆ usedHits() [3/4]

double ITRT_ToT_dEdx::usedHits ( const Trk::Track *  track,
bool  useHitsHT = true 
) const
inlinevirtualinherited

Definition at line 99 of file ITRT_ToT_dEdx.h.

100 {
101  return usedHits(Gaudi::Hive::currentContext(), track, useHitsHT);
102 }

◆ usedHits() [4/4]

double ITRT_ToT_dEdx::usedHits
inline

Definition at line 60 of file ITRT_ToT_dEdx.h.

100 {
101  return usedHits(Gaudi::Hive::currentContext(), track, useHitsHT);
102 }

Member Data Documentation

◆ m_assoTool

ToolHandle< Trk::IPRD_AssociationTool > TRT_ToT_dEdx::m_assoTool {this, "AssociationTool","InDet::InDetPRD_AssociationToolGangedPixels","PRD_AssociationTool"}
private

Definition at line 70 of file TRT_ToT_dEdx.h.

◆ m_corrected

BooleanProperty TRT_ToT_dEdx::m_corrected
private
Initial value:
{this, "TRT_dEdx_corrected", true,
"If true - make correction using rs-distributions"}

Definition at line 75 of file TRT_ToT_dEdx.h.

◆ m_correctionType

IntegerProperty TRT_ToT_dEdx::m_correctionType {this, "TRT_dEdx_correctionType", kTrackBased, "Type of dEdx correction"}
private

Definition at line 86 of file TRT_ToT_dEdx.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_divideByL

BooleanProperty TRT_ToT_dEdx::m_divideByL
private
Initial value:
{this, "TRT_dEdx_divideByL", true,
"If true - divide ToT to the L of track in straw"}

Definition at line 77 of file TRT_ToT_dEdx.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_isData

BooleanProperty TRT_ToT_dEdx::m_isData {this, "TRT_dEdx_isData", true}
private

Definition at line 99 of file TRT_ToT_dEdx.h.

◆ m_localOccTool

ToolHandle< InDet::ITRT_LocalOccupancy > TRT_ToT_dEdx::m_localOccTool {this, "TRT_LocalOccupancyTool","", "Local occupancy tool"}
private

Definition at line 71 of file TRT_ToT_dEdx.h.

◆ m_nTrunkateHits

unsigned int TRT_ToT_dEdx::m_nTrunkateHits = 1
private

Definition at line 110 of file TRT_ToT_dEdx.h.

◆ m_rdhkEvtInfo

SG::ReadDecorHandleKey<xAOD::EventInfo> TRT_ToT_dEdx::m_rdhkEvtInfo
private
Initial value:
{this
, "averageInteractionsPerCrossingKey"
, "EventInfo.averageInteractionsPerCrossing"
, "Decoration for Average Interaction Per Crossing"}

Definition at line 61 of file TRT_ToT_dEdx.h.

◆ m_ReadKey

SG::ReadCondHandleKey<TRTDedxcorrection> TRT_ToT_dEdx::m_ReadKey {this,"Dedxcorrection","Dedxcorrection","Dedx constants in-key"}
private

Definition at line 256 of file TRT_ToT_dEdx.h.

◆ m_toolScenario

IntegerProperty TRT_ToT_dEdx::m_toolScenario
private
Initial value:
{this, "TRT_dEdx_toolScenario", kAlgReweightTrunkOne,
"Algorithm type for dEdX estimator calculation"}

Definition at line 83 of file TRT_ToT_dEdx.h.

◆ m_trackConfig_maxRtrack

FloatProperty TRT_ToT_dEdx::m_trackConfig_maxRtrack {this, "TRT_dEdx_trackConfig_maxRtrack", 1.85, "maximum track radius"}
private

Definition at line 102 of file TRT_ToT_dEdx.h.

◆ m_trackConfig_minRtrack

FloatProperty TRT_ToT_dEdx::m_trackConfig_minRtrack {this, "TRT_dEdx_trackConfig_minRtrack", 0.15, "maximum track radius"}
private

Definition at line 104 of file TRT_ToT_dEdx.h.

◆ m_trtId

const TRT_ID* TRT_ToT_dEdx::m_trtId {nullptr}
private

Definition at line 65 of file TRT_ToT_dEdx.h.

◆ m_TRTStrawSummaryTool

ToolHandle<ITRT_StrawStatusSummaryTool> TRT_ToT_dEdx::m_TRTStrawSummaryTool {this, "TRTStrawSummaryTool", "TRT_StrawStatusSummaryTool"}
private

Definition at line 67 of file TRT_ToT_dEdx.h.

◆ m_useTrackPartWithGasType

IntegerProperty TRT_ToT_dEdx::m_useTrackPartWithGasType
private
Initial value:
{this, "TRT_dEdx_useTrackPartWithGasType", kUnset,
"If kUnset - use any gas for dEdX calculation"}

Definition at line 80 of file TRT_ToT_dEdx.h.

◆ m_useZeroRHitCut

BooleanProperty TRT_ToT_dEdx::m_useZeroRHitCut
private
Initial value:
{this, "TRT_dEdx_useZeroRHitCut", true,
"skip tracks with where RHit=0"}

Definition at line 107 of file TRT_ToT_dEdx.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

The documentation for this class was generated from the following files:
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
TRT_ToT_dEdx::m_TRTStrawSummaryTool
ToolHandle< ITRT_StrawStatusSummaryTool > m_TRTStrawSummaryTool
Definition: TRT_ToT_dEdx.h:68
beamspotman.r
def r
Definition: beamspotman.py:676
Trk::TrackStateOnSurface::trackParameters
const TrackParameters * trackParameters() const
return ptr to trackparameters const overload
TRT_ToT_dEdx::usedHits
virtual double usedHits(const EventContext &ctx, const Trk::Track *track, bool useHThits=true) const override final
function to calculate number of used hits
Definition: TRT_ToT_dEdx.cxx:320
TRT_ToT_dEdx::isGoodHit
bool isGoodHit(const EventContext &ctx, const Trk::TrackStateOnSurface *trackState, bool useHitsHT, double &length) const
function to define what is a good hit to be used for dEdx calculation cuts on track level can be made...
Definition: TRT_ToT_dEdx.cxx:78
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
DataModel_detail::const_iterator
Const iterator class for DataVector/DataList.
Definition: DVLIterator.h:82
InDet::TRT_DriftCircleOnTrack::prepRawData
virtual const TRT_DriftCircle * prepRawData() const override final
returns the PrepRawData - is a TRT_DriftCircle in this scope
Definition: TRT_DriftCircleOnTrack.h:202
python.SystemOfUnits.s
int s
Definition: SystemOfUnits.py:131
TRT_ToT_dEdx::kTrackBased
@ kTrackBased
Definition: TRT_ToT_dEdx.h:56
SGout2dot.alg
alg
Definition: SGout2dot.py:243
get_generator_info.result
result
Definition: get_generator_info.py:21
TRT_ToT_dEdx::getProb
double getProb(const EventContext &ctx, const Trk::TrackStateOnSurface *itr, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, int nUsedHits) const
function to calculate likelihood from prediction and resolution
Definition: TRT_ToT_dEdx.cxx:399
python.Constants.FATAL
int FATAL
Definition: Control/AthenaCommon/python/Constants.py:19
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
TRT_ToT_dEdx::EGasType
EGasType
Definition: TRT_ToT_dEdx.h:57
index
Definition: index.py:1
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
hist_file_dump.d
d
Definition: hist_file_dump.py:137
TRT_ToT_dEdx::fitFuncBarrelShort_corrRZ
double fitFuncBarrelShort_corrRZ(const EventContext &ctx, EGasType gasType, double driftRadius, double zPosition, int StrawLayer) const
function called by fitFuncBarrel_corrRZ for short straws
Definition: TRT_ToT_dEdx.cxx:839
InDet::TRT_DriftCircle::timeOverThreshold
double timeOverThreshold() const
returns Time over threshold in ns
TRT_ToT_dEdx::kKrypton
@ kKrypton
Definition: TRT_ToT_dEdx.h:57
accumulate
bool accumulate(AccumulateMap &map, std::vector< module_t > const &modules, FPGATrackSimMatrixAccumulator const &acc)
Accumulates an accumulator (e.g.
Definition: FPGATrackSimMatrixAccumulator.cxx:22
theta
Scalar theta() const
theta method
Definition: AmgMatrixBasePlugin.h:75
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
TRTCalib_cfilter.p1
p1
Definition: TRTCalib_cfilter.py:130
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
TRT_ToT_dEdx::getTest
virtual double getTest(const EventContext &ctx, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, Trk::ParticleHypothesis antihypothesis, int nUsedHits) const override final
function to calculate likelihood ratio test
Definition: TRT_ToT_dEdx.cxx:467
TRT_ToT_dEdx::dEdx
virtual double dEdx(const EventContext &ctx, const Trk::Track *track, bool useHThits, std::optional< float > localOccupancy=std::nullopt) const override final
function to calculate sum ToT normalised to number of used hits
Definition: TRT_ToT_dEdx.cxx:133
Trk::TrackStateOnSurface::measurementOnTrack
const MeasurementBase * measurementOnTrack() const
returns MeasurementBase const overload
athena.value
value
Definition: athena.py:124
Trk::RIO_OnTrack::rioType
virtual bool rioType(RIO_OnTrackType::Type type) const =0
Method checking the Rio On Track type.
ITRT_ToT_dEdx::usedHits
virtual double usedHits(const EventContext &ctx, const Trk::Track *track, bool useHitsHT=true) const =0
function to calculate number of used hits
UploadAMITag.l
list l
Definition: UploadAMITag.larcaf.py:158
Trk::RIO_OnTrack
Definition: RIO_OnTrack.h:70
InDet::TRT_DriftCircleOnTrack::highLevel
bool highLevel() const
returns true if the high level threshold was passed
Definition: TRT_DriftCircleOnTrack.h:234
ITRT_ToT_dEdx::dEdx
virtual double dEdx(const EventContext &ctx, const Trk::Track *track, bool useHitsHT, std::optional< float > localOccupancy=std::nullopt) const =0
function to calculate sum ToT normalised to number of used hits
Trk::locR
@ locR
Definition: ParamDefs.h:44
TRT_ToT_dEdx::trackOccupancyCorrection
double trackOccupancyCorrection(const EventContext &ctx, const Trk::Track *track, bool useHThits, std::optional< float > localOccupancy) const
Definition: TRT_ToT_dEdx.cxx:1149
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
xAOD::identifier
identifier
Definition: UncalibratedMeasurement_v1.cxx:15
drawFromPickle.exp
exp
Definition: drawFromPickle.py:36
covarianceTool.prob
prob
Definition: covarianceTool.py:678
InDet::TRT_DriftCircleOnTrack
Definition: TRT_DriftCircleOnTrack.h:53
InDet::TRT_DriftCircleOnTrack::detectorElement
virtual const InDetDD::TRT_BaseElement * detectorElement() const override final
returns the detector element, assoicated with the PRD of this class
Definition: TRT_DriftCircleOnTrack.h:224
Trk::ParametersCommon::dim
static constexpr int dim
Definition: ParametersCommon.h:50
std::sort
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.
Definition: DVL_algorithms.h:554
TRT_ToT_dEdx::kNewDB
@ kNewDB
Definition: TRT_ToT_dEdx.h:54
python.RingerConstants.Layer
Layer
Definition: RingerConstants.py:42
TRT_ToT_dEdx::fitFuncBarrelLong_corrRZ
double fitFuncBarrelLong_corrRZ(const EventContext &ctx, EGasType gasType, double driftRadius, double zPosition, int Layer, int StrawLayer) const
function called by fitFuncBarrel_corrRZ for long straws
Definition: TRT_ToT_dEdx.cxx:809
TRT_ToT_dEdx::kRSOnly
@ kRSOnly
Definition: TRT_ToT_dEdx.h:56
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
Trk::TrackStateOnSurface::Outlier
@ Outlier
This TSoS contains an outlier, that is, it contains a MeasurementBase/RIO_OnTrack which was not used ...
Definition: TrackStateOnSurface.h:122
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:27
TRT_ToT_dEdx::correctToT_corrRZ
double correctToT_corrRZ(const EventContext &ctx, const Trk::TrackStateOnSurface *itr, double length) const
main function to correct ToT values on hit level as a function of track radius and z-position
Definition: TRT_ToT_dEdx.cxx:708
TRT_ToT_dEdx::fitFuncBarrel_corrRZL
double fitFuncBarrel_corrRZL(const EventContext &ctx, EGasType gasType, double driftRadius, double zPosition, int Layer, int StrawLayer) const
function to compute correction factor in barrel region
Definition: TRT_ToT_dEdx.cxx:990
Trk::TrackStateOnSurface::type
bool type(const TrackStateOnSurfaceType type) const
Use this method to find out if the TSoS is of a certain type: i.e.
TRTCalib_cfilter.p2
p2
Definition: TRTCalib_cfilter.py:131
TRT_ToT_dEdx::kArgon
@ kArgon
Definition: TRT_ToT_dEdx.h:57
TRT_ToT_dEdx::kUnset
@ kUnset
Definition: TRT_ToT_dEdx.h:57
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
TRT_ToT_dEdx::fitFuncBarrel_corrRZ
double fitFuncBarrel_corrRZ(const EventContext &ctx, EGasType gas, double driftRadius, double zPosition, int Layer, int StrawLayer) const
function to compute correction factor in barrel region
Definition: TRT_ToT_dEdx.cxx:779
SG::ReadDecorHandle
Handle class for reading a decoration on an object.
Definition: StoreGate/StoreGate/ReadDecorHandle.h:94
TRT_ToT_dEdx::m_assoTool
ToolHandle< Trk::IPRD_AssociationTool > m_assoTool
Definition: TRT_ToT_dEdx.h:70
StdJOSetup.msgSvc
msgSvc
Provide convenience handles for various services.
Definition: StdJOSetup.py:36
lumiFormat.i
int i
Definition: lumiFormat.py:85
z
#define z
python.CaloCondTools.g
g
Definition: CaloCondTools.py:15
beamspotman.n
n
Definition: beamspotman.py:731
Trk::theta
@ theta
Definition: ParamDefs.h:66
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
Trk::electron
@ electron
Definition: ParticleHypothesis.h:27
extractSporadic.h
list h
Definition: extractSporadic.py:97
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AmgVector
AmgVector(4) T2BSTrackFilterTool
Definition: T2BSTrackFilterTool.cxx:114
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
Trk::driftRadius
@ driftRadius
trt, straws
Definition: ParamDefs.h:53
TRT_ToT_dEdx::kHitBased
@ kHitBased
Definition: TRT_ToT_dEdx.h:56
calibdata.exception
exception
Definition: calibdata.py:496
ITRT_ToT_dEdx::getTest
virtual double getTest(const EventContext &ctx, const double dEdx_obs, const double pTrk, Trk::ParticleHypothesis hypothesis, Trk::ParticleHypothesis antihypothesis, int nUsedHits) const =0
function to calculate likelihood ratio test
sign
int sign(int a)
Definition: TRT_StrawNeighbourSvc.h:107
TRTDedxcorrection
Definition: TRTDedxcorrection.h:9
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
hist_file_dump.f
f
Definition: hist_file_dump.py:135
Trk::MeasurementBase::type
virtual bool type(MeasurementBaseType::Type type) const =0
Interface method checking the type.
drawFromPickle.tan
tan
Definition: drawFromPickle.py:36
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
Trk::ParametersBase
Definition: ParametersBase.h:55
TRT_ToT_dEdx::mass
double mass(const EventContext &ctx, const Trk::TrackStateOnSurface *itr, const double pTrk, double dEdx) const
function to extract most likely mass in bg [0,3]
Definition: TRT_ToT_dEdx.cxx:570
TRT_ID::barrel_ec
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
Definition: TRT_ID.h:866
TRT_ID::straw_layer
int straw_layer(const Identifier &id) const
Definition: TRT_ID.h:893
TRT_ToT_dEdx::kAlgStandard
@ kAlgStandard
Definition: TRT_ToT_dEdx.h:55
TRT_ToT_dEdx::kXenon
@ kXenon
Definition: TRT_ToT_dEdx.h:57
TRT_ToT_dEdx::kGlobal
@ kGlobal
Definition: TRT_ToT_dEdx.h:56
DataVector< const Trk::TrackStateOnSurface >
TRT_ToT_dEdx::m_useZeroRHitCut
BooleanProperty m_useZeroRHitCut
Definition: TRT_ToT_dEdx.h:107
TRT_ID::layer_or_wheel
int layer_or_wheel(const Identifier &id) const
Definition: TRT_ID.h:884
beamspotman.stat
stat
Definition: beamspotman.py:266
TRT_ToT_dEdx::m_rdhkEvtInfo
SG::ReadDecorHandleKey< xAOD::EventInfo > m_rdhkEvtInfo
Definition: TRT_ToT_dEdx.h:61
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
Trk::MeasurementBase::localCovariance
const Amg::MatrixX & localCovariance() const
Interface method to get the localError.
Definition: MeasurementBase.h:138
inRange
bool inRange(const double *boundaries, const double value, const double tolerance=0.02)
Definition: LArSCIdVsIdTest.cxx:5
CxxUtils::set
constexpr std::enable_if_t< is_bitmask_v< E >, E & > set(E &lhs, E rhs)
Convenience function to set bits in a class enum bitmask.
Definition: bitmask.h:232
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
Trk::RIO_OnTrackType::TRT_DriftCircle
@ TRT_DriftCircle
Definition: RIO_OnTrack.h:59
trigbs_pickEvents.num
num
Definition: trigbs_pickEvents.py:76
TRT_ToT_dEdx::m_trtId
const TRT_ID * m_trtId
Definition: TRT_ToT_dEdx.h:65
Trk::MeasurementBase
Definition: MeasurementBase.h:58
Trk::ParticleMasses::mass
constexpr double mass[PARTICLEHYPOTHESES]
the array of masses
Definition: ParticleHypothesis.h:53
TRT_ToT_dEdx::m_ReadKey
SG::ReadCondHandleKey< TRTDedxcorrection > m_ReadKey
Definition: TRT_ToT_dEdx.h:256
TRT_ToT_dEdx::fitFuncEndcap_corrRZ
double fitFuncEndcap_corrRZ(const EventContext &ctx, EGasType gas, double driftRadius, double rPosition, int Layer, int sign) const
function to compute correction factor in endcap region
Definition: TRT_ToT_dEdx.cxx:793
TRT_ToT_dEdx::predictdEdx
double predictdEdx(const EventContext &ctx, const Trk::TrackStateOnSurface *itr, const double pTrk, Trk::ParticleHypothesis hypothesis) const
function to calculate expectation value for dEdx using BB fit
Definition: TRT_ToT_dEdx.cxx:490
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
TRT_ToT_dEdx::kOldDB
@ kOldDB
Definition: TRT_ToT_dEdx.h:54
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
TRT_ToT_dEdx::correctNormalization
double correctNormalization(const EventContext &ctx, double nVtx=-1) const
correct overall dEdx normalization on track level
Definition: TRT_ToT_dEdx.cxx:684
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
TRT_ID::phi_module
int phi_module(const Identifier &id) const
Definition: TRT_ID.h:875
TRT_ToT_dEdx::calculateTrackLengthInStraw
static double calculateTrackLengthInStraw(const Trk::TrackStateOnSurface *trackState, const TRT_ID *identifier)
Definition: TRT_ToT_dEdx.cxx:1185
Trk::RIO_OnTrack::prepRawData
virtual const Trk::PrepRawData * prepRawData() const =0
returns the PrepRawData (also known as RIO) object to which this RIO_OnTrack is associated.
Trk::MeasurementBaseType::RIO_OnTrack
@ RIO_OnTrack
Definition: MeasurementBase.h:49
TRT_ToT_dEdx::m_correctionType
IntegerProperty m_correctionType
Definition: TRT_ToT_dEdx.h:87
TRT_ToT_dEdx::m_localOccTool
ToolHandle< InDet::ITRT_LocalOccupancy > m_localOccTool
Definition: TRT_ToT_dEdx.h:71
REPORT_MESSAGE
#define REPORT_MESSAGE(LVL)
Report a message.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:365
Trk::MeasurementBase::localParameters
const LocalParameters & localParameters() const
Interface method to get the LocalParameters.
Definition: MeasurementBase.h:132
DataVector::end
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
python.PyAthena.v
v
Definition: PyAthena.py:154
TRT_ToT_dEdx::fitFuncPol_corrRZ
double fitFuncPol_corrRZ(const EventContext &ctx, EGasType gasType, int parameter, double driftRadius, int Layer, int Strawlayer, int sign, int set) const
function called by fitFuncBarrel_corrRZ and fitFuncEndcap_corrRZ
Definition: TRT_ToT_dEdx.cxx:857
DeMoScan.index
string index
Definition: DeMoScan.py:364
TRT_ToT_dEdx::fitFuncEndcap_corrRZL
double fitFuncEndcap_corrRZL(const EventContext &ctx, EGasType gasType, double driftRadius, double radialPosition, int Layer, int sign) const
function to compute correction factor in endcap region
Definition: TRT_ToT_dEdx.cxx:939
a
TList * a
Definition: liststreamerinfos.cxx:10
h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
TRT_ToT_dEdx::hitOccupancyCorrection
double hitOccupancyCorrection(const EventContext &ctx, const Trk::TrackStateOnSurface *itr) const
Definition: TRT_ToT_dEdx.cxx:1059
Trk::RIO_OnTrack::identify
Identifier identify() const
return the identifier -extends MeasurementBase
Definition: RIO_OnTrack.h:152
InDetDD::TRT_BaseElement::center
virtual const Amg::Vector3D & center() const override final
Element Surface: center of a straw layer.
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
TRT_ToT_dEdx::kAlgReweight
@ kAlgReweight
Definition: TRT_ToT_dEdx.h:55
convertTimingResiduals.offset
offset
Definition: convertTimingResiduals.py:71
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
timeOverThreshold
double timeOverThreshold(unsigned int m_word)
Definition: driftCircle.h:116
Trk::phi
@ phi
Definition: ParamDefs.h:75
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
TRT_ToT_dEdx::m_isData
BooleanProperty m_isData
Definition: TRT_ToT_dEdx.h:99
TRT_ToT_dEdx::m_useTrackPartWithGasType
IntegerProperty m_useTrackPartWithGasType
Definition: TRT_ToT_dEdx.h:81
xAOD::track
@ track
Definition: TrackingPrimitives.h:512
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
InDet::TRT_DriftCircleOnTrack::globalPosition
virtual const Amg::Vector3D & globalPosition() const override final
return the global position of this RIO_OnTrack
Definition: TRT_DriftCircleOnTrack.cxx:160
CaloNoise_fillDB.mu
mu
Definition: CaloNoise_fillDB.py:53
error
Definition: IImpactPoint3dEstimator.h:70
TRT_ToT_dEdx::m_corrected
BooleanProperty m_corrected
Definition: TRT_ToT_dEdx.h:75
TRT_ToT_dEdx::m_toolScenario
IntegerProperty m_toolScenario
Definition: TRT_ToT_dEdx.h:84
python.compressB64.c
def c
Definition: compressB64.py:93
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
TRT_ToT_dEdx::kAlgReweightTrunkOne
@ kAlgReweightTrunkOne
Definition: TRT_ToT_dEdx.h:55
length
double length(const pvec &v)
Definition: FPGATrackSimLLPDoubletHoughTransformTool.cxx:26
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
TRT_ToT_dEdx::gasTypeInStraw
EGasType gasTypeInStraw(const EventContext &ctx, const Trk::TrackStateOnSurface *itr) const
return gas type for that hit
Definition: TRT_ToT_dEdx.cxx:622
TRT_ToT_dEdx::m_divideByL
BooleanProperty m_divideByL
Definition: TRT_ToT_dEdx.h:77
TRT_ToT_dEdx::m_trackConfig_minRtrack
FloatProperty m_trackConfig_minRtrack
Definition: TRT_ToT_dEdx.h:105
TauGNNUtils::Variables::Track::trackEta
bool trackEta(const xAOD::TauJet &, const xAOD::TauTrack &track, double &out)
Definition: TauGNNUtils.cxx:477
TRTCalib_cfilter.p0
p0
Definition: TRTCalib_cfilter.py:129
DataVector::begin
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
fitman.k
k
Definition: fitman.py:528
TRT_ToT_dEdx::m_nTrunkateHits
unsigned int m_nTrunkateHits
Definition: TRT_ToT_dEdx.h:110
InDetDD::TRT_BaseElement
Definition: TRT_BaseElement.h:57
Identifier
Definition: IdentifierFieldParser.cxx:14
TRT_ToT_dEdx::m_trackConfig_maxRtrack
FloatProperty m_trackConfig_maxRtrack
Definition: TRT_ToT_dEdx.h:103
Generated on Thu Nov 7 2024 22:06:59 for ATLAS Offline Software by   doxygen 1.8.18