MultiTrajectory.icc

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
 
constexpr bool ActsTrk::MutableMultiTrajectory::hasColumn_impl(
    Acts::HashedString key) const {
  using namespace Acts::HashedStringLiteral;
 
  switch (key) {
    case "previous"_hash:
    case "chi2"_hash:
    case "pathLength"_hash:
    case "typeFlags"_hash:
    case "predicted"_hash:
    case "filtered"_hash:
    case "smoothed"_hash:
    case "jacobian"_hash:
    case "projector"_hash:
    case "uncalibratedSourceLink"_hash:
    case "calibrated"_hash:
    case "calibratedCov"_hash:
    case "measdim"_hash:
    case "referenceSurface"_hash:
      return true;
    default:
      for (auto& d : m_decorations) {
        if (d.hash == key) {
          return true;
        }
      }
      return false;
  }
}
 
template <typename T>
void ActsTrk::MutableMultiTrajectory::addColumn_impl(std::string_view name) {
  // It is actually not clear if we would allow decorating RO MTJ, maybe we do
  if constexpr (ActsTrk::detail::accepted_decoration_types<T>::value) {
    m_decorations.emplace_back( ActsTrk::detail::decoration<T>(
        name,
        ActsTrk::detail::constDecorationGetter<T>,
        ActsTrk::detail::decorationCopier<T>,
        ActsTrk::detail::decorationSetter<T>
    )
  );
  INSPECTCALL("added dynamic decoration " << name << " n decorations now " << m_decorations.size());
 
    // it would be useful to force presence of decoration already here
  } else {
    throw std::runtime_error("Can't add decoration of this type to MutableMultiTrajectory");
  }
}
 
template <typename val_t, typename cov_t>
void ActsTrk::MutableMultiTrajectory::allocateCalibrated_impl(IndexType istate,
                              const Eigen::DenseBase<val_t>& val,
                              const Eigen::DenseBase<cov_t>& cov)
 
  requires(Eigen::PlainObjectBase<val_t>::RowsAtCompileTime > 0 &&
            Eigen::PlainObjectBase<val_t>::RowsAtCompileTime <= Acts::eBoundSize &&
            Eigen::PlainObjectBase<val_t>::RowsAtCompileTime ==
                Eigen::PlainObjectBase<cov_t>::RowsAtCompileTime &&
            Eigen::PlainObjectBase<cov_t>::RowsAtCompileTime ==
                Eigen::PlainObjectBase<cov_t>::ColsAtCompileTime)
{
  constexpr std::size_t measdim = val_t::RowsAtCompileTime;
 
  if(m_trackStatesAux->measDim.at(istate) != kInvalid
     && m_trackStatesAux->measDim.at(istate) != measdim) {
    throw std::invalid_argument{
      "Measurement dimension does not match the allocated dimension"};
  }
 
  m_trackStatesAux->measDim.at(istate) = measdim;
  auto idx = m_trackStatesAux->calibrated.at(istate);
  m_trackMeasurementsAux->meas.at(idx).resize(measdim);
  m_trackMeasurementsAux->covMatrix.at(idx).resize(measdim*measdim);
 
  double* measPtr = m_trackMeasurementsAux->meas[idx].data();
  Eigen::Map<Acts::ActsVector<measdim>> valMap(measPtr);
  valMap = val;
 
  double* covPtr = m_trackMeasurementsAux->covMatrix[idx].data();
  Eigen::Map<Acts::ActsSquareMatrix<measdim>> covMap(covPtr);
  covMap = cov;
}