#include <AddTRTMomConstr.h>

Inheritance diagram for AddTRTMomConstr:

Collaboration diagram for AddTRTMomConstr:

Public Member Functions
	AddTRTMomConstr (const std::string &name, ISvcLocator *pSvcLocator)
	The AddTRTMomConstr is an implementation to add the so-called TRT momentum constraint on a track.
	~AddTRTMomConstr ()
StatusCode	initialize ()
	initialize method of this algorithm.
StatusCode	execute ()
	execute method of this algorithm that is called for each event
StatusCode	finalize ()
	finalize method of this algorithm.
MsgStream &	dump (MsgStream &outst) const
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
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.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
const Trk::TrackStateOnSurface *	findouterscthit (const Trk::Track *track)
	Returns the outermost SCT TSOS on the given Track.
const Trk::TrackStateOnSurface *	findinnertrthit (const Trk::Track *track)
	Returns the innermost TRT TSOS on the given Track.
bool	accept (const Trk::Track &track)
	Verifies if the given track passes the track selection criteria specified via the jobOptions.
const Trk::MeasurementSet	addPM (Trk::MeasurementSet &ms, const Trk::PseudoMeasurementOnTrack *pm)
	adds a PseudoMeasurement to a MeasurementSet
const Trk::PseudoMeasurementOnTrack *	createPMfromSi (const Trk::Perigee *mp)
	creates a PseudoMeasurement with (z0, theta) from extended track perigee parameters
const Trk::PseudoMeasurementOnTrack *	createPMfromTRT (const Trk::Perigee mpSi, const Trk::Perigee mpTRT)
	creates a PseudoMeasurement with (z0, theta, q/p) from TRT track perigee parameters
Trk::Track *	addTRTMomentumConstraint (const Trk::Track *track)
	Strips of all TRT hits from the track and replaces them with a TRT momentum constraint from a the TRT segment belonging to the (extended) track (returns NULL in case of failure)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
const AtlasDetectorID *	m_idHelper
	Detector ID helper.
const TRT_ID *	m_trtid
	Identifier helper class for the TRT.
std::string	m_trackListInput
	Name of the TrackCollection (input)
std::string	m_trackListOutput
	Name of the TrackCollection (Output)
ToolHandle< Trk::ITrackSummaryTool >	m_trackSummaryTool
	The summary tool to make the track selection.
ToolHandle< Trk::ITrackFitter >	m_trackFitter
	The TrackFitter to refit the tracks (segment, momentum constraint)
bool	m_useThetaCorrection
	correct theta after TRT+PM(z0,theta) refit, small diff.
double	m_thetaCorr
	theta can be different after TRT+PM fit by < o(10e-4).
bool	m_applyTrkSel
	apply a selection on tracks or not
double	m_selPtMin
	minimal pT cut value for the TrackSelection
double	m_selEtaMin
	minimal eta cut value for the TrackSelection
double	m_selEtaMax
	maximal eta cut value for the TrackSelection
int	m_selNHitPIXMin
	minimal number of PIX hits cut value for the TrackSelection
int	m_selNHitSCTMin
	minimal number of SCT hits cut value for the TrackSelection
int	m_selNHitTRTMin
	minimal number of TRT hits cut value for the TrackSelection
int	m_selNHitTRTMinCrack
	minimal number of TRT hits cut value in trans.
double	m_selEtaCrackMin
	minimum \|eta\| to define the transition region
double	m_selEtaCrackMax
	maximum \|eta\| to define the transition region
double	m_selChiSqPerDOFMin
	maximal Chisquare per degree of freedom cut value for the TrackSelection
size_t	m_nTracksProcessed
	Counter for number of tracks processed.
size_t	m_nTracksPresel
	Counter for number of tracks passing the preselection.
size_t	m_nTracksAccepted
	Counter for number of tracks passing the preselection and with PM.
size_t	m_nRejectPIX
	Counter for number of tracks failing the min number of PIX hits req.
size_t	m_nRejectSCT
	Counter for number of tracks failing the min number of SCT hits req.
size_t	m_nRejectTRT
	Counter for number of tracks failing the min number of TRT hits req.
size_t	m_nRejectPM
	Counter for number of tracks failing the addition of a pseudo-measurement (PM)
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 36 of file AddTRTMomConstr.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ AddTRTMomConstr()

AddTRTMomConstr::AddTRTMomConstr	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

The AddTRTMomConstr is an implementation to add the so-called TRT momentum constraint on a track.

It was developed in the context of alignment, and its main idea is that if the TRT is free from "weak modes", its curvature measurement can be added as a PseudoMeasurement (PM) to a Si-only track. This extra information, taken into account while refitting the Si-only part of the track or read in by the alignment algorithms directly, can serve to eliminate the alignment weak modes from the Si tracker of ATLAS. More details can be found at: https://twiki.cern.ch/twiki/bin/view/Atlas/TRTMomConstraint

Definition at line 34 of file AddTRTMomConstr.cxx.

  : AthAlgorithm   (name, pSvcLocator)
  , m_idHelper(nullptr)
  , m_trtid(nullptr)
  , m_trackListInput("ExtendedTracks")
  , m_trackListOutput("SiTracksWithTRTMomConstr")
  , m_trackSummaryTool("Trk::ITrackSummaryTool/ITrackSummaryTool")
  , m_trackFitter("Trk::KalmanFitter/TrkKalmanFitter")
  , m_useThetaCorrection(false)
  , m_thetaCorr(1.)
  , m_applyTrkSel(false)
  , m_selPtMin(0)
  , m_selEtaMin(-2.5)
  , m_selEtaMax(2.5)
  , m_selNHitPIXMin(0)
  , m_selNHitSCTMin(0)
  , m_selNHitTRTMin(20)
  , m_selNHitTRTMinCrack(15)
  , m_selEtaCrackMin(0.7)
  , m_selEtaCrackMax(1.)
  , m_selChiSqPerDOFMin(100)
  , m_nTracksProcessed(0)
  , m_nTracksPresel(0)
  , m_nTracksAccepted(0)
  , m_nRejectPIX(0)
  , m_nRejectSCT(0)
  , m_nRejectTRT(0)
  , m_nRejectPM(0)
{
  // declare algorithm parameters
  declareProperty("TrackListInput"     , m_trackListInput) ;
  declareProperty("TrackListOutput"    , m_trackListOutput) ;
  declareProperty("TrackFitter"        , m_trackFitter) ;
  declareProperty("TrackSummaryTool"   , m_trackSummaryTool) ;
  declareProperty("UseThetaCorrection" , m_useThetaCorrection) ;
  declareProperty("ApplyTrkSel"        , m_applyTrkSel) ;
  declareProperty("SelPtMin"           , m_selPtMin) ;
  declareProperty("SelEtaMin"          , m_selEtaMin) ;
  declareProperty("SelEtaMax"          , m_selEtaMax) ;
  declareProperty("SelNHitPIXMin"      , m_selNHitPIXMin) ;
  declareProperty("SelNHitSCTMin"      , m_selNHitSCTMin) ;
  declareProperty("SelNHitTRTMin"      , m_selNHitTRTMin) ;
  declareProperty("SelNHitTRTMinCrack" , m_selNHitTRTMinCrack) ;
  declareProperty("SelEtaCrackMin"     , m_selEtaCrackMin) ;
  declareProperty("SelEtaCrackMax"     , m_selEtaCrackMax) ;
  declareProperty("SelChiSqPerDOFMin"  , m_selChiSqPerDOFMin) ;
}

◆ ~AddTRTMomConstr()

AddTRTMomConstr::~AddTRTMomConstr ( )

Definition at line 82 of file AddTRTMomConstr.cxx.

82{}

Member Function Documentation

◆ accept()

bool AddTRTMomConstr::accept ( const Trk::Track & track )

private

Verifies if the given track passes the track selection criteria specified via the jobOptions.

Definition at line 149 of file AddTRTMomConstr.cxx.

                                                    {
  // selection for alignment quality tracks
  ATH_MSG_DEBUG( "Inside accept(.)") ;
 
  bool rc=true ;
  const Trk::Track* pTrack = &track ;
  const Trk::Perigee* mesp = pTrack->perigeeParameters();
  if( !mesp ) {
    ATH_MSG_WARNING("MeasuredPerigee pointer 0, track not accepted") ;
    return false ;
  }
 
  float theta = mesp->parameters()[Trk::theta] ;
  float pt = 0. ;
  if( theta != 0 ){
    float ptInv = std::abs( mesp->parameters()[Trk::qOverP] ) / std::sin( theta ) ;
    if( ptInv != 0 ) {
      pt = 1./ptInv ;
      if( m_selPtMin > 0 && pt < m_selPtMin ) rc = false ;
    }
  } else {
    ATH_MSG_WARNING( "theta = 0, q/p conversion to p_T failed!" );
    return false ;
  }
 
  float eta = -std::log( std::tan( theta/2. ) ) ;
  if( m_selEtaMin < m_selEtaMax && (eta < m_selEtaMin || eta > m_selEtaMax ) ) {
    rc = false ;
  }
 
  float chisqpdof = track.fitQuality()->chiSquared() / track.fitQuality()->numberDoF() ;
  if( m_selChiSqPerDOFMin > 0 && chisqpdof > m_selChiSqPerDOFMin ) rc = false ;
  std::unique_ptr<const Trk::TrackSummary> summary
    = m_trackSummaryTool->summary(Gaudi::Hive::currentContext(), *pTrack ) ;
  if( summary->get(Trk::numberOfPixelHits) < m_selNHitPIXMin ) {
    ++m_nRejectPIX ;
    rc = false;
  }
  if( summary->get(Trk::numberOfSCTHits) < m_selNHitSCTMin ) {
    ++m_nRejectSCT ;
    rc = false ;
  }
  // different treatment for TRT in transition region
  bool isInTransitionRegion = std::abs(eta) > m_selEtaCrackMin && std::abs(eta) < m_selEtaCrackMax ;
  if( isInTransitionRegion ) {
    if( int(summary->get(Trk::numberOfTRTHits)) < m_selNHitTRTMinCrack ) {
      ++m_nRejectTRT ;
      rc = false ;
    }
  } else {
    if( int(summary->get(Trk::numberOfTRTHits)) < m_selNHitTRTMin ) {
      ++m_nRejectTRT ;
      rc = false ;
    }
  }
 
  // All cuts passed...
  ATH_MSG_DEBUG( "Track selection:"
        << "\n\t passed : " << rc
        << "\n\t Pt     = " << pt
        << "\n\t phi    = " << mesp->parameters()[Trk::phi]
        << "\n\t eta    = " << eta
        << "\n\t nPIX   = " << summary->get(Trk::numberOfPixelHits)
        << "\n\t nSCT   = " << summary->get(Trk::numberOfSCTHits)
        << "\n\t nTRT   = " << summary->get(Trk::numberOfTRTHits)
               ) ;
  ATH_MSG_DEBUG( "leaving accept(.)" ) ;
  return rc ;
} // accept(.)

◆ addPM()

const Trk::MeasurementSet AddTRTMomConstr::addPM	(	Trk::MeasurementSet &	ms,
		const Trk::PseudoMeasurementOnTrack *	pm )

private

adds a PseudoMeasurement to a MeasurementSet

Definition at line 329 of file AddTRTMomConstr.cxx.

                                                                                                             {
  Trk::MeasurementSet sortedMS ;
  sortedMS.push_back( pm ) ;
  for( int i=0, i_max=ms.size() ; i!=i_max ; ++i ) {
    sortedMS.push_back( ms[i] ) ;
  }
  return sortedMS ;
}

◆ addTRTMomentumConstraint()

Trk::Track * AddTRTMomConstr::addTRTMomentumConstraint ( const Trk::Track * track )

private

Strips of all TRT hits from the track and replaces them with a TRT momentum constraint from a the TRT segment belonging to the (extended) track (returns NULL in case of failure)

Definition at line 338 of file AddTRTMomConstr.cxx.

                                                                         {
  Trk::MeasurementSet setSi ;
  Trk::MeasurementSet setTRT ;
  //store all silicon measurements into the measurementset
  for (const auto it:*(track->measurementsOnTrack())){
    const Trk::RIO_OnTrack* rio = dynamic_cast <const Trk::RIO_OnTrack*>(it);
    if (rio) {
      const Identifier& surfaceID = (rio->identify()) ;
      if( m_idHelper->is_sct(surfaceID) || m_idHelper->is_pixel(surfaceID) ) {
        setSi.push_back ( it ) ;
      }
      if( m_idHelper->is_trt(surfaceID) ) {
        setTRT.push_back ( it ) ;
      }
    }
  }
  ATH_MSG_DEBUG("TRTMomConstr() : Found " << setSi.size()  << " Si  measurm's!" ) ;
  ATH_MSG_DEBUG("TRTMomConstr() : Found " << setTRT.size() << " TRT measurm's!") ;
  if( int(setTRT.size()) < m_selNHitTRTMin ) {
    ATH_MSG_DEBUG( "TRTMomConstr() : fewer TRT measurements than required: "
               << setTRT.size() << " < " << m_selNHitTRTMin << ", skip track" ) ;
    return nullptr ;
  }
 
  ATH_MSG_DEBUG( std::setiosflags( std::ios::scientific )) ;
  ATH_MSG_DEBUG (std::setprecision( 7 )) ;
  ATH_MSG_VERBOSE  ( std::setiosflags( std::ios::scientific )) ;
  ATH_MSG_VERBOSE  ( std::setprecision( 7 )) ;
  const Trk::Perigee* perTrk = track->perigeeParameters();
  if( !perTrk ) {
    ATH_MSG_WARNING("TRTMomConstr() : No Perigee parameter on track!");
    return nullptr ;
  }
 
  // now add z_0 and theta from original track as PseudoMeas to the TRT MeasurementSet
  const Trk::PseudoMeasurementOnTrack *pmFromSi = createPMfromSi( perTrk ) ;
  if( !pmFromSi ) {
    ATH_MSG_ERROR( "TRTMomConstr() : PseudoMeasurementOnTrack creation failed! " );
    return nullptr ;
  }
  ATH_MSG_DEBUG( "TRTMomConstr() : pmFromSi " << *pmFromSi) ;
  Trk::MeasurementSet setTRTPM = addPM( setTRT, pmFromSi ) ;
 
   ATH_MSG_VERBOSE ( "TRTMomConstr() : TRT+PM MeasurementSet : " );
  for( int i=0, i_max=setTRTPM.size() ; i!=i_max ; ++i ) {
    ATH_MSG_VERBOSE ("============== i=" << i << " =============");
    ATH_MSG_VERBOSE ( *(setTRTPM[i]));
  }
   ATH_MSG_VERBOSE ("==========================================");
 
  // fit TRT part of the track with PseudoMeas on z_0, theta
   Trk::Track* trkTRT = (m_trackFitter->fit(
     Gaudi::Hive::currentContext(), setTRTPM, *perTrk, true, Trk::pion
     //, Trk::muon
     //, Trk::nonInteracting
   )).release();
   if (!trkTRT) {
     ATH_MSG_DEBUG("TRTMomConstr() : Fit of TRT part of the track failed! ");
     return nullptr;
   }
  const Trk::Perigee* perTRT = trkTRT->perigeeParameters();
  ATH_MSG_DEBUG( "TRTMomConstr() : perTRT " << *perTRT) ;
  // the theta value after TRT+PM fit can be different from the initial one by < o(10e-4). Correct q/p optionally
  if( m_useThetaCorrection ) m_thetaCorr = std::sin( perTrk->parameters()[Trk::theta] )
                                         / std::sin( perTRT->parameters()[Trk::theta] ) ;
  ATH_MSG_DEBUG( "TRTMomConstr() : Scalefactor to correct q/p: " << m_thetaCorr) ;
 
  // define new PM with the momentum constraint from the TRT to pass to Si
  const Trk::PseudoMeasurementOnTrack *pmFromTRT = createPMfromTRT( perTrk, perTRT ) ;
  if( !pmFromTRT ) {
    ATH_MSG_ERROR( "TRTMomConstr() : PseudoMeasurementOnTrack creation failed! " );
    return nullptr ;
  }
  ATH_MSG_DEBUG("TRTMomConstr() : pmFromTRT " << *pmFromTRT ) ;
  Trk::MeasurementSet setSiPM = addPM( setSi, pmFromTRT ) ;
 
  // fit Si part of the track with PM from TRT
  Trk::Track* fittedTrack =
    m_trackFitter
      ->fit(Gaudi::Hive::currentContext(), setSiPM, *perTrk, true, Trk::pion
            //, Trk::muon
            //, Trk::nonInteracting
            ).release();
  if (!fittedTrack) {
    ATH_MSG_DEBUG( "TRTMomConstr() : Si+TRT-p_T Track fit failed !" ) ;
  } else {
    const Trk::Perigee* perSi = fittedTrack->perigeeParameters();
 
    ATH_MSG_DEBUG( "TRTMomConstr() : Si+PM(TRT) track parameteres @ perigee: " << *perSi ) ;
  }
 
  return fittedTrack;
} // TRTMomentumConstraintForExtTracks

◆ createPMfromSi()

const Trk::PseudoMeasurementOnTrack * AddTRTMomConstr::createPMfromSi ( const Trk::Perigee * mp )

private

creates a PseudoMeasurement with (z0, theta) from extended track perigee parameters

Definition at line 269 of file AddTRTMomConstr.cxx.

                                                                                          {
  Trk::DefinedParameter z0   ( mp->parameters()[Trk::z0],    Trk::z0    ) ;
  Trk::DefinedParameter theta( mp->parameters()[Trk::theta], Trk::theta ) ;
  std::array<Trk::DefinedParameter, 2> defPar = {z0, theta};
  if( !mp->covariance() )  return nullptr;
  Trk::LocalParameters  parFromSi( defPar ) ;
  AmgSymMatrix(2) covFromSi;
 
  covFromSi( 0, 0 ) = (*mp->covariance())( Trk::z0,Trk::z0 ) ;
  covFromSi( 1, 1 ) = (*mp->covariance())( Trk::theta,Trk::theta ) ;
  covFromSi( 1, 0 ) = (*mp->covariance())( Trk::z0, Trk::theta ) ;
  covFromSi( 0, 1 ) = (*mp->covariance())( Trk::z0, Trk::theta ) ;
 
  const Trk::Surface& mpSurf = mp->associatedSurface() ;
 
  Trk::PseudoMeasurementOnTrack *pm = new Trk::PseudoMeasurementOnTrack( std::move(parFromSi)
                                                                       , covFromSi
                                                                       , mpSurf
                                                                       ) ;
  return pm ;
}

◆ createPMfromTRT()

const Trk::PseudoMeasurementOnTrack * AddTRTMomConstr::createPMfromTRT	(	const Trk::Perigee *	mpSi,
		const Trk::Perigee *	mpTRT )

private

creates a PseudoMeasurement with (z0, theta, q/p) from TRT track perigee parameters

< from Si

< from TRT

Definition at line 291 of file AddTRTMomConstr.cxx.

                                                                                                 {
 
  if( !mpSi->covariance() || !mpTRT->covariance() )  return nullptr;
 
  // since z0, theta not measured by TRT, take back initial Si values with blown up errors.
  Trk::DefinedParameter z0TRT    ( mpSi->parameters()[Trk::z0],                 Trk::z0     ) ; 
  Trk::DefinedParameter thetaTRT ( mpSi->parameters()[Trk::theta],              Trk::theta  ) ; 
  Trk::DefinedParameter qOvewPTRT( mpTRT->parameters()[Trk::qOverP]*m_thetaCorr, Trk::qOverP ) ; 
  std::array<Trk::DefinedParameter, 3> parFromTRT_vec = {z0TRT, thetaTRT,
                                                         qOvewPTRT};
  Trk::LocalParameters parFromTRT(parFromTRT_vec);
  double A = 1e6 ; // scale up the errorfor z0 and theta to avoid double-counting when including the PM
 
 
  AmgSymMatrix(3) covFromTRT;
 
  // take the z0 theta part from extended track (correl. messed up by TRT fit) and scale it up with A
  covFromTRT( 0, 0 ) = (*mpSi->covariance())( Trk::z0, Trk::z0      ) * A * A ;
  covFromTRT( 1, 1 ) = (*mpSi->covariance())( Trk::theta, Trk::theta  ) * A * A ;
  covFromTRT( 1, 0 ) = (*mpSi->covariance())( Trk::z0   , Trk::theta  ) * A * A ;
  covFromTRT( 0, 1 ) = (*mpSi->covariance())( Trk::z0   , Trk::theta  ) * A * A ;
  // scale up the cov terms dependant on z0 and theta by A, and m_thetaCorr from q/p. Same m_thetaCorr for theta
  covFromTRT( 2, 2 ) = (*mpTRT->covariance())( Trk::qOverP, Trk::qOverP ) * m_thetaCorr * m_thetaCorr ;
  covFromTRT( 2, 0 ) = (*mpTRT->covariance())( Trk::z0   , Trk::qOverP ) * A * m_thetaCorr ;
  covFromTRT( 2, 1 ) = (*mpTRT->covariance())( Trk::theta, Trk::qOverP ) * A * m_thetaCorr * m_thetaCorr ;
  covFromTRT( 0, 2 ) = (*mpTRT->covariance())( Trk::z0   , Trk::qOverP ) * A * m_thetaCorr ;
  covFromTRT( 1, 2 ) = (*mpTRT->covariance())( Trk::theta, Trk::qOverP ) * A * m_thetaCorr * m_thetaCorr ;
 
  const Trk::Surface& mpTRTSurf = mpSi->associatedSurface()  ;
 
  Trk::PseudoMeasurementOnTrack *pm = new Trk::PseudoMeasurementOnTrack( std::move(parFromTRT)
                                                                       , covFromTRT
                                                                       , mpTRTSurf
                                                                       ) ;
  return pm ;
}

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	hndl,
		const SG::VarHandleKeyType &	)

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Algorithm > >::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; }

AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ dump()

MsgStream & AddTRTMomConstr::dump ( MsgStream & outst ) const

Definition at line 436 of file AddTRTMomConstr.cxx.

                                                         {
  outst << "\n|-------------------------------------------------------------------";
  outst << "-----------------------------|\n" ;
  outst << "|  processed                      : "
        << std::setw(7) << m_nTracksProcessed
        << " tracks                                               |\n";
  outst << "|  accepted by track presel.      : "
        << std::setw(7) << m_nTracksPresel
        << " tracks                                               |\n";
  outst << "|  accepted by track presel. + PM : "
        << std::setw(7) << m_nTracksAccepted
        << " tracks                                               |\n";
  outst << "| ------------------------------------------------------------------";
  outst << "---------------------------- |\n" ;
  outst << "|  reject by # PIX hits           : "
        << std::setw(7) << m_nRejectPIX
        << " tracks                                               |\n";
  outst << "|  reject by # SCT hits           : "
        << std::setw(7) << m_nRejectSCT
        << " tracks                                               |\n";
  outst << "|  reject by # TRT hits           : "
        << std::setw(7) << m_nRejectTRT
        << " tracks                                               |\n";
  outst << "| ------------------------------------------------------------------";
  outst << "---------------------------- |\n" ;
  outst << "|  reject by exist. PM(TRT)       : "
        << std::setw(7) << m_nRejectPM
        << " tracks                                               |\n";
 
  outst << "|-------------------------------------------------------------------";
  outst << "-----------------------------|\n" ;
  return outst ;
}

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Algorithm > >::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; }

AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode AddTRTMomConstr::execute ( )

execute method of this algorithm that is called for each event

Definition at line 117 of file AddTRTMomConstr.cxx.

                                    {
  if( !m_trackListOutput.empty() && !m_trackListInput.empty() ) {
    const TrackCollection* inputtracks = nullptr;
    ATH_CHECK(evtStore()->retrieve( inputtracks, m_trackListInput) );
    auto outputtracks = std::make_unique<ConstDataVector<TrackCollection> >( SG::VIEW_ELEMENTS ) ;
    for (const auto it : *inputtracks){
      auto track=it;
      if( m_applyTrkSel && !accept( *track ) ) {
        ATH_MSG_DEBUG( "kinematic requirements not passed, skip track") ;
        continue ;
      }
      ++m_nTracksPresel ;
      track = addTRTMomentumConstraint( track ) ;
      if( track ) {
        ATH_MSG_VERBOSE ( "addTRTMomentumConstraint(.) produced new track" ) ;
        outputtracks->push_back( track ) ;
      } else {
        ATH_MSG_DEBUG( "addTRTMomentumConstraint(.) returned 0 pointer, skip track") ;
        ++m_nRejectPM ;
      }
    }
 
    m_nTracksProcessed += inputtracks->size() ;
    m_nTracksAccepted += outputtracks->size() ;
 
    if( evtStore()->record( std::move(outputtracks), m_trackListOutput ).isFailure() ) {
      ATH_MSG_ERROR( "Failed to record trackcollection with name " << m_trackListOutput );
    }
  }
  return StatusCode::SUCCESS;
} // execute(.)

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return Algorithm::extraOutputDeps();
}

◆ finalize()

StatusCode AddTRTMomConstr::finalize ( )

finalize method of this algorithm.

Prints out a summary of all events

Definition at line 112 of file AddTRTMomConstr.cxx.

                                     {
  ATH_MSG_DEBUG( (*this) ) ;
  return StatusCode::SUCCESS;
}

◆ findinnertrthit()

const Trk::TrackStateOnSurface * AddTRTMomConstr::findinnertrthit ( const Trk::Track * track )

private

Returns the innermost TRT TSOS on the given Track.

Definition at line 244 of file AddTRTMomConstr.cxx.

                                                                                      {
  ATH_MSG_VERBOSE ( "Inside findinnerscthit: " << track->trackStateOnSurfaces()->size() );
  double rmin=9999999. ;
  const Trk::TrackStateOnSurface* rc=nullptr ;
  for (const Trk::TrackStateOnSurface* tsos : *track->trackStateOnSurfaces()) {
    if( !tsos->type(Trk::TrackStateOnSurface::Outlier) &&
         tsos->measurementOnTrack() &&
         tsos->trackParameters() ) {
      ATH_MSG_VERBOSE ( "Found good measuredtrackparameters" );
      const InDet::TRT_DriftCircleOnTrack* trthit = dynamic_cast<const InDet::TRT_DriftCircleOnTrack*>(tsos->measurementOnTrack());
      if(trthit) {
        ATH_MSG_VERBOSE ( "Found TRT_DriftCircleOnTrack" );
        const Amg::Vector3D& pos = trthit->globalPosition() ;
        double r = sqrt(pos.x()*pos.x() + pos.y()*pos.y()) ;
        if(rc==nullptr || r<rmin) {
          rc = tsos ;
          rmin = r ;
        }
      }
    }
  }
  ATH_MSG_VERBOSE ( "Returning surface with adress: " << rc << " at radius " << rmin );
  return rc ;
}

◆ findouterscthit()

const Trk::TrackStateOnSurface * AddTRTMomConstr::findouterscthit ( const Trk::Track * track )

private

Returns the outermost SCT TSOS on the given Track.

Definition at line 219 of file AddTRTMomConstr.cxx.

                                                                                      {
  ATH_MSG_VERBOSE ("Inside findouterscthit: " << track->trackStateOnSurfaces()->size());
  double rmax=0. ;
  const Trk::TrackStateOnSurface* rc=nullptr ;
  for (const Trk::TrackStateOnSurface* tsos : *track->trackStateOnSurfaces()) {
    if( !tsos->type(Trk::TrackStateOnSurface::Outlier) &&
        tsos->measurementOnTrack() &&
        tsos->trackParameters() ) {
      ATH_MSG_VERBOSE ( "Found good measuredtrackparameters" );
      const InDet::SCT_ClusterOnTrack* sctclus = dynamic_cast<const InDet::SCT_ClusterOnTrack*>(tsos->measurementOnTrack());
      if(sctclus) {
        ATH_MSG_VERBOSE ( "Found SCT_ClusterOnTrack");
        const Amg::Vector3D& pos = sctclus->globalPosition() ;
        double r = std::sqrt(pos.x()*pos.x() + pos.y()*pos.y()) ;
        if(rc==nullptr || r>rmax) {
          rc = tsos ;
          rmax = r ;
        }
      }
    }
  }
  ATH_MSG_VERBOSE ( "Returning surface with adress: " << rc << " at radius " << rmax);
  return rc ;
}

◆ initialize()

StatusCode AddTRTMomConstr::initialize ( )

initialize method of this algorithm.

Definition at line 84 of file AddTRTMomConstr.cxx.

                                       {
  //retrieve the DetectorStore service
  ATH_CHECK(detStore().retrieve()) ;
  ATH_CHECK(detStore()->retrieve(m_trtid));
  ATH_CHECK(detStore()->retrieve(m_idHelper, "AtlasID"));
  // get TrackSummaryTool
  ATH_CHECK(m_trackSummaryTool.retrieve());
  // Set up fitter
  ATH_CHECK(m_trackFitter.retrieve());
  // Print input properties
  if( m_applyTrkSel ) {
    ATH_MSG_DEBUG(
             "Track selection will be applied:"
          << "\n " << m_selEtaMin << " < eta < " << m_selEtaMax
          << "\n trans. region: " << m_selEtaCrackMin << " < |eta| < " << m_selEtaCrackMax
          << "\n pT                            >  " << m_selPtMin/Gaudi::Units::GeV << " GeV"
          << "\n number_of_PIXEL_hits          >= " << m_selNHitPIXMin
          << "\n number_of_SCT_hits            >= " << m_selNHitSCTMin
          << "\n number_of_TRT_hits            >= " << m_selNHitTRTMin
          << "\n number_of_TRT_hits in crack   >= " << m_selNHitTRTMinCrack
          << "\n chi2/dof                      <  " << m_selChiSqPerDOFMin
                ) ;
  } else {
    ATH_MSG_DEBUG( "NO selection will be applied to tracks" ) ;
  }
  return StatusCode::SUCCESS;
} // initialize(.)

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ msg()

MsgStream & AthCommonMsg< Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ 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< Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, and PyAthena::Alg.

Definition at line 66 of file AthAlgorithm.cxx.

                                       {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;
}

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

◆ m_applyTrkSel

bool AddTRTMomConstr::m_applyTrkSel

private

apply a selection on tracks or not

Definition at line 82 of file AddTRTMomConstr.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_idHelper

const AtlasDetectorID* AddTRTMomConstr::m_idHelper

private

Detector ID helper.

Definition at line 72 of file AddTRTMomConstr.h.

◆ m_nRejectPIX

size_t AddTRTMomConstr::m_nRejectPIX

private

Counter for number of tracks failing the min number of PIX hits req.

Definition at line 96 of file AddTRTMomConstr.h.

◆ m_nRejectPM

size_t AddTRTMomConstr::m_nRejectPM

private

Counter for number of tracks failing the addition of a pseudo-measurement (PM)

Definition at line 99 of file AddTRTMomConstr.h.

◆ m_nRejectSCT

size_t AddTRTMomConstr::m_nRejectSCT

private

Counter for number of tracks failing the min number of SCT hits req.

Definition at line 97 of file AddTRTMomConstr.h.

◆ m_nRejectTRT

size_t AddTRTMomConstr::m_nRejectTRT

private

Counter for number of tracks failing the min number of TRT hits req.

Definition at line 98 of file AddTRTMomConstr.h.

◆ m_nTracksAccepted

size_t AddTRTMomConstr::m_nTracksAccepted

private

Counter for number of tracks passing the preselection and with PM.

Definition at line 95 of file AddTRTMomConstr.h.

◆ m_nTracksPresel

size_t AddTRTMomConstr::m_nTracksPresel

private

Counter for number of tracks passing the preselection.

Definition at line 94 of file AddTRTMomConstr.h.

◆ m_nTracksProcessed

size_t AddTRTMomConstr::m_nTracksProcessed

private

Counter for number of tracks processed.

Definition at line 93 of file AddTRTMomConstr.h.

◆ m_selChiSqPerDOFMin

double AddTRTMomConstr::m_selChiSqPerDOFMin

private

maximal Chisquare per degree of freedom cut value for the TrackSelection

Definition at line 92 of file AddTRTMomConstr.h.

◆ m_selEtaCrackMax

double AddTRTMomConstr::m_selEtaCrackMax

private

maximum |eta| to define the transition region

Definition at line 91 of file AddTRTMomConstr.h.

◆ m_selEtaCrackMin

double AddTRTMomConstr::m_selEtaCrackMin

private

minimum |eta| to define the transition region

Definition at line 90 of file AddTRTMomConstr.h.

◆ m_selEtaMax

double AddTRTMomConstr::m_selEtaMax

private

maximal eta cut value for the TrackSelection

Definition at line 85 of file AddTRTMomConstr.h.

◆ m_selEtaMin

double AddTRTMomConstr::m_selEtaMin

private

minimal eta cut value for the TrackSelection

Definition at line 84 of file AddTRTMomConstr.h.

◆ m_selNHitPIXMin

int AddTRTMomConstr::m_selNHitPIXMin

private

minimal number of PIX hits cut value for the TrackSelection

Definition at line 86 of file AddTRTMomConstr.h.

◆ m_selNHitSCTMin

int AddTRTMomConstr::m_selNHitSCTMin

private

minimal number of SCT hits cut value for the TrackSelection

Definition at line 87 of file AddTRTMomConstr.h.

◆ m_selNHitTRTMin

int AddTRTMomConstr::m_selNHitTRTMin

private

minimal number of TRT hits cut value for the TrackSelection

Definition at line 88 of file AddTRTMomConstr.h.

◆ m_selNHitTRTMinCrack

int AddTRTMomConstr::m_selNHitTRTMinCrack

private

minimal number of TRT hits cut value in trans.

region for the TrackSelection

Definition at line 89 of file AddTRTMomConstr.h.

◆ m_selPtMin

double AddTRTMomConstr::m_selPtMin

private

minimal pT cut value for the TrackSelection

Definition at line 83 of file AddTRTMomConstr.h.

◆ m_thetaCorr

double AddTRTMomConstr::m_thetaCorr

private

theta can be different after TRT+PM fit by < o(10e-4).

Correct q/p optionally

Definition at line 81 of file AddTRTMomConstr.h.

◆ m_trackFitter

ToolHandle<Trk::ITrackFitter> AddTRTMomConstr::m_trackFitter

private

The TrackFitter to refit the tracks (segment, momentum constraint)

Definition at line 77 of file AddTRTMomConstr.h.

◆ m_trackListInput

std::string AddTRTMomConstr::m_trackListInput

private

Name of the TrackCollection (input)

Definition at line 74 of file AddTRTMomConstr.h.

◆ m_trackListOutput

std::string AddTRTMomConstr::m_trackListOutput

private

Name of the TrackCollection (Output)

Definition at line 75 of file AddTRTMomConstr.h.

◆ m_trackSummaryTool

ToolHandle<Trk::ITrackSummaryTool> AddTRTMomConstr::m_trackSummaryTool

private

The summary tool to make the track selection.

Definition at line 76 of file AddTRTMomConstr.h.

◆ m_trtid

const TRT_ID* AddTRTMomConstr::m_trtid

private

Identifier helper class for the TRT.

Definition at line 73 of file AddTRTMomConstr.h.

◆ m_useThetaCorrection

bool AddTRTMomConstr::m_useThetaCorrection

private

correct theta after TRT+PM(z0,theta) refit, small diff.

typically < o(3x10e-4)

Definition at line 80 of file AddTRTMomConstr.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ AddTRTMomConstr()

◆ ~AddTRTMomConstr()

Member Function Documentation

◆ accept()

◆ addPM()

◆ addTRTMomentumConstraint()

◆ createPMfromSi()

◆ createPMfromTRT()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ dump()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ finalize()

◆ findinnertrthit()

◆ findouterscthit()

◆ initialize()

◆ inputHandles()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_applyTrkSel

◆ m_detStore

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_idHelper

◆ m_nRejectPIX

◆ m_nRejectPM

◆ m_nRejectSCT

◆ m_nRejectTRT

◆ m_nTracksAccepted

◆ m_nTracksPresel

◆ m_nTracksProcessed

◆ m_selChiSqPerDOFMin

◆ m_selEtaCrackMax

◆ m_selEtaCrackMin

◆ m_selEtaMax

◆ m_selEtaMin

◆ m_selNHitPIXMin

◆ m_selNHitSCTMin

◆ m_selNHitTRTMin

◆ m_selNHitTRTMinCrack

◆ m_selPtMin

◆ m_thetaCorr

◆ m_trackFitter

◆ m_trackListInput

◆ m_trackListOutput

◆ m_trackSummaryTool

◆ m_trtid

◆ m_useThetaCorrection

◆ m_varHandleArraysDeclared

◆ m_vhka