SphericalFitFunction.h

Go to the documentation of this file.

#ifndef _RdSphericalFit_SphericalFitFunction_h_
#define _RdSphericalFit_SphericalFitFunction_h_

#include <fwk/LocalCoordinateSystem.h>

#include <utl/CoordinateSystemPtr.h>
#include <utl/Vector.h>
#include <utl/Minou.h>

#include <vector>

namespace RdSphericalFit {

  struct StationFitData {
    double signalTime;
    double signalTimeErr;
    utl::Vector antennaPos;
  };


  class SphericalFitFunction : public utl::Minou::Base {

  public:
    SphericalFitFunction( const std::vector<StationFitData>& stationData,
                          const std::vector<utl::Minou::ParameterDef>& pars,
                          const utl::CoordinateSystemPtr& localCS,
                          const double tmean) :
      fStationData(stationData),
      fLocalCS(localCS),
      ftmean(tmean)
    {
      GetParameterDefs() = pars;
    }

    double operator()(const std::vector<double>& pars) const;
    double GetChi2(const std::vector<double>& pars);
    int GetNDF();

  private:
    const std::vector<StationFitData>& fStationData;
    const utl::CoordinateSystemPtr& fLocalCS;
    const double ftmean;

    // fit parameters
    mutable double fRadius = 0;
    mutable double fTheta = 0;
    mutable double fPhi = 0;
    mutable double fGamma = 1;

    void
    UpdateParameters(const std::vector<double>& pars)
    const
    {
      fRadius = pars[0];
      fTheta = pars[1];
      fPhi = pars[2];
      fGamma = pars[3];
    }
  };

}


#endif