PerigeeParameters.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Oliver Frost                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/geometry/UncertainParameters.h>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    namespace NPerigeeParameterIndices {
 
      enum EPerigeeParameter {
        c_Curv = 0,
 
        c_Phi0 = 1,
 
        c_I = 2,
 
        c_N = 3,
      };
    }
 
    using EPerigeeParameter = NPerigeeParameterIndices::EPerigeeParameter;
 
    // Guard to prevent repeated template symbol emission
    struct PerigeeUtil;
    extern template struct UncertainParametersUtil<PerigeeUtil, EPerigeeParameter>;
 
    struct PerigeeUtil : UncertainParametersUtil<PerigeeUtil, EPerigeeParameter> {
 
      static ParameterVector reversalSigns()
      {
        return ParameterVector({ -1.0, 1.0, -1.0});
      }
 
      static double average(const PerigeeUtil::ParameterVector& fromPar,
                            const PerigeeUtil::CovarianceMatrix& fromCov,
                            const PerigeeUtil::ParameterVector& toPar,
                            const PerigeeUtil::CovarianceMatrix& toCov,
                            PerigeeUtil::ParameterVector& avgPar,
                            PerigeeUtil::CovarianceMatrix& avgCov);
    };
 
    using PerigeeParameters = PerigeeUtil::ParameterVector;
 
    using PerigeeCovariance = PerigeeUtil::CovarianceMatrix;
 
    using PerigeePrecision = PerigeeUtil::PrecisionMatrix;
 
    using PerigeeJacobian = PerigeeUtil::JacobianMatrix;
 
  }
 
}