FitInterface.h

Go to the documentation of this file.

#ifndef _SdHorizontalReconstruction_FitInterface_h_
#define _SdHorizontalReconstruction_FitInterface_h_

#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/CovarianceMatrix.h>
#include <vector>

namespace SdHorizontalReconstructionNS {

  // (corex,corey) in barycenter CS
  // (theta, phi) wrt barycenter
  // (distance, ctime) of origin wrt barycenter
  // (coreXExt,coreYExt) in Barycenter CS
  // (thetaExt, phiExt) in Barycenter CS

  enum SizeDataParameter {
    eN19 = 0, eCoreX = 1, eCoreY = 2
  };

  enum AxisDataParameter {
    eTheta = 0, ePhi = 1, eDistance = 2, eCTime = 3
  };

  enum ExternalGeometryParameter {
    eCoreXExt = 0, eCoreYExt = 1,
    eThetaExt = 2, ePhiExt = 3
  };

  struct SizeData {

    static const std::size_t size = 3;

    SizeData()
      : fCov(size),
        fRecStage(0),
        fNFree(0),
        fNLogLike(0)
    {}

    double fPar[size];
    utl::CovarianceMatrix fCov;
    double fRecStage;
    int fNFree;
    double fNLogLike;
  };

  struct AxisData {

    static const std::size_t size = 4;

    AxisData()
      : fCov(size),
        fNFree(0),
        fChi2(0)
    {}

    double fPar[size];
    utl::CovarianceMatrix fCov;
    int fNFree;
    double fChi2;
  };

  struct SilentStationData {
    int fId;
    utl::Point fPos;
  };

  struct StationData : public SilentStationData {
    bool fRejected;

    double fSignal;
    double fT50;
    double fCTime; // relative to barycenter time
    bool fSaturated;
    double fRecoveryErr;
  };

  typedef std::vector<SilentStationData> SilentStationList;
  typedef std::vector<StationData> StationList;

  struct ExternalGeometryData {

    static const std::size_t size = 4;

    ExternalGeometryData()
    {
      for (size_t i = 0; i < 4; ++i)
      {
        fPar[i] = 0;
        for (size_t j = 0; j < 4; ++j)
          fInvCov[i][j] = 0;
      }
    }

    double fPar[size];
    double fInvCov[size][size];
  };

} // SdHorizontalReconstructionNS

#endif