#include <GlobalChiSquareFitterTool.h>

Collaboration diagram for ActsTrk::PRDSourceLinkCalibratorGX2F:

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
	The following is the implementation of the PRDSourceLinkCalibratorGX2F. More...

Public Attributes
const Trk::IRIO_OnTrackCreator *	rotCreator {nullptr}

const Trk::IRIO_OnTrackCreator *	broadRotCreator {nullptr}

const ActsTrk::IActsToTrkConverterTool *	converterTool {nullptr}

Detailed Description

Definition at line 59 of file GlobalChiSquareFitterTool.h.

Member Function Documentation

◆ calibrate()

template<typename trajectory_t >

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

The following is the implementation of the PRDSourceLinkCalibratorGX2F.

To use this implementation, the flag for a refitting from PRD measurement has to be turned on in the appropriate Acts config.

Definition at line 57 of file GlobalChiSquareFitterTool.cxx.

                                                            {
  
   const Trk::PrepRawData* prd = sl.template get<PRDSourceLinkGX2F>().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 5
   // 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::PRDSourceLinkCalibratorGX2F' 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);
  
   // Writing the calibrated ROT measurement 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.setProjectorSubspaceIndices(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.setProjectorSubspaceIndices(subspaceIndices);
   } else {
     throw std::runtime_error("Dim " + std::to_string(dim) +
                              " currently not supported.");
   }
   delete rot;  // Delete rot as a precaution
 }