ActsTrk::PRDSourceLinkCalibrator Struct Reference

#include <KalmanFitter.h>

Collaboration diagram for ActsTrk::PRDSourceLinkCalibrator:

Public Member Functions
template<typename trajectory_t >
void	calibrate (const Acts::GeometryContext &gctx, const Acts::CalibrationContext &cctx, const Acts::SourceLink &sl, typename trajectory_t::TrackStateProxy trackState) const

Public Attributes
const Trk::IRIO_OnTrackCreator *	rotCreator {nullptr}
const Trk::IRIO_OnTrackCreator *	broadRotCreator {nullptr}
const ActsTrk::IActsToTrkConverterTool *	converterTool {nullptr}

Detailed Description

Definition at line 69 of file KalmanFitter.h.

Member Function Documentation

calibrate()

template<typename trajectory_t >

void ActsTrk::PRDSourceLinkCalibrator::calibrate	(	const Acts::GeometryContext &	gctx,
		const Acts::CalibrationContext &	cctx,
		const Acts::SourceLink &	sl,
		typename trajectory_t::TrackStateProxy	trackState
	)		const

Definition at line 61 of file KalmanFitter.cxx.

                                                           {
    
    const Trk::PrepRawData* prd = sl.template get<PRDSourceLink>().prd;
    trackState.setUncalibratedSourceLink(Acts::SourceLink{sl});
 
    const Acts::BoundTrackParameters actsParam(trackState.referenceSurface().getSharedPtr(),
             trackState.predicted(),
             trackState.predictedCovariance(),
             Acts::ParticleHypothesis::pion());
    std::unique_ptr<const Trk::TrackParameters> trkParam = converterTool->actsTrackParametersToTrkParameters(actsParam, gctx);
    
    //RIO_OnTrack creation from the PrepRawData
    const Trk::Surface & prdsurf = prd->detectorElement()->surface(prd->identify());
    const Trk::RIO_OnTrack *rot = nullptr;
    const Trk::PlaneSurface *plsurf = nullptr;
    if (prdsurf.type() == Trk::SurfaceType::Plane)
      plsurf = static_cast < const Trk::PlaneSurface *>(&prdsurf);
 
    const Trk::StraightLineSurface *slsurf = nullptr;
        
    if (prdsurf.type() == Trk::SurfaceType::Line)
      slsurf = static_cast < const Trk::StraightLineSurface *>(&prdsurf);
 
    //@TODO, there is no way to put a MSG in this framework yet. So the next 3 lines are commented
    //if ((slsurf == nullptr) && (plsurf == nullptr)) {
      //msg(MSG::ERROR) << "Surface is neither PlaneSurface nor StraightLineSurface!" << endmsg; //ATH_MSG_ERROR doesn't work here because (ActsTrk::PRDSourceLinkCalibrator' has no member named 'msg')
    // }
    if (slsurf != nullptr) {
      Trk::AtaStraightLine atasl(
        slsurf->center(), 
        trackState.predicted()[Trk::phi],
        trackState.predicted()[Trk::theta],
        trackState.predicted()[Trk::qOverP], 
        *slsurf
       );
      if(broadRotCreator){
        rot = broadRotCreator->correct(*prd, atasl,Gaudi::Hive::currentContext()); 
        }
      else{
        rot = rotCreator->correct(*prd, atasl,Gaudi::Hive::currentContext()); 
        }
     } else if (plsurf != nullptr) {
      if ((trkParam.get())->covariance() != nullptr) { 
        Trk::AtaPlane atapl(
          plsurf->center(), 
          trackState.predicted()[Trk::phi],
          trackState.predicted()[Trk::theta],
          trackState.predicted()[Trk::qOverP], 
          *plsurf,
          AmgSymMatrix(5)(*trkParam.get()->covariance()) 
        );
        rot = rotCreator->correct(*prd, atapl,Gaudi::Hive::currentContext());
      } else {
        Trk::AtaPlane atapl(
          plsurf->center(), 
          trackState.predicted()[Trk::phi],
          trackState.predicted()[Trk::theta],
          trackState.predicted()[Trk::qOverP], 
          *plsurf
        );
        rot = rotCreator->correct(*prd, atapl,Gaudi::Hive::currentContext()); 
      }
    } // End of RIO_OnTrack creation from the PrepRawData
 
    int dim = (*rot).localParameters().dimension();
    trackState.allocateCalibrated(dim);
    
    //Writting the calibrated ROT measurment into the trackState
    if (dim == 0)
      {
        throw std::runtime_error("Cannot create dim 0 measurement");
      } else if (dim == 1) {
        trackState.template calibrated<1>() = (*rot).localParameters().template head<1>(); 
        trackState.template calibratedCovariance<1>() = (*rot).localCovariance().template topLeftCorner<1, 1>(); 
        Acts::BoundSubspaceIndices subspaceIndices;
        if ((*rot).associatedSurface().bounds().type() == Trk::SurfaceBounds::Annulus) { 
          subspaceIndices = {Acts::eBoundLoc1}; // y coordinate is l0
        } else {
          subspaceIndices = {Acts::eBoundLoc0}; // x coordinate is l0
        }
        trackState.setBoundSubspaceIndices(subspaceIndices);
      }
      else if (dim == 2)
        {
          trackState.template calibrated<2>() = (*rot).localParameters().template head<2>();
          trackState.template calibratedCovariance<2>() = (*rot).localCovariance().template topLeftCorner<2, 2>();
          Acts::BoundSubspaceIndices subspaceIndices = {Acts::eBoundLoc0, Acts::eBoundLoc1};
          trackState.setBoundSubspaceIndices(subspaceIndices);
        }
      else
      {
        throw std::runtime_error("Dim " + std::to_string(dim) +
                                " currently not supported.");
      }
    delete rot; //Delete rot as a precaution
 
}

Member Data Documentation

broadRotCreator

const Trk::IRIO_OnTrackCreator* ActsTrk::PRDSourceLinkCalibrator::broadRotCreator {nullptr}

Definition at line 77 of file KalmanFitter.h.

converterTool

const ActsTrk::IActsToTrkConverterTool* ActsTrk::PRDSourceLinkCalibrator::converterTool {nullptr}

Definition at line 78 of file KalmanFitter.h.

rotCreator

const Trk::IRIO_OnTrackCreator* ActsTrk::PRDSourceLinkCalibrator::rotCreator {nullptr}

Definition at line 76 of file KalmanFitter.h.

---

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

• KalmanFitter.h
• KalmanFitter.cxx