PCGFitter.h

Go to the documentation of this file.

#ifndef _FdProfileConstrainedGeometryFit_PCGFitter_h_
#define _FdProfileConstrainedGeometryFit_PCGFitter_h_

#include <utl/Minou.h>

#include "ChiZeroRegression.h"
#include "ProfileChi2.h"
#include "ApertureLight.h"

#include <evt/Event.h>

#include <fevt/Eye.h>
#include <fevt/EyeRecData.h>

#include "../FdApertureLightFinderKG/FdApertureLightFinder.h"


namespace FdProfileConstrainedGeometryFit {

  class PCGFitter : public utl::Minou::Base {
  public:
    PCGFitter() : Base({{"chi0"}}) { }
    void Init();

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

    bool Run(fevt::Eye &eye);

  public:
    enum EApLightMethod {
      eExternal,
      eInternal,
      eCallForEach
    };

  protected:
    fevt::Eye* fEye = nullptr;
    ChiZeroRegression* fChi0Regression = nullptr;
    ProfileChi2 fProfileFit;

    FdApertureLightFinderKG::FdApertureLightFinder fApertureLightFinder;
    ApertureLight fApLight;
    EApLightMethod fApLightMethod;
    bool fCheckUnderground = false;
    bool fPrescan = false;
    bool fScanOnly = false;
    double fScanStep = 0;
    double fScanStart = 0;
    double fScanStop = 0;
    bool fUseLightFlux = false;
    bool fDelZeroFlux = false;
    bool fRealAtm = false;
    bool fDeex = false;

    fevt::Eye& CopyEye(const fevt::Eye& eye, evt::Event& eventCopy) const;

    bool ScanChi0(fevt::Eye &eye,
                  double step = 1*utl::degree,
                  double from = 10*utl::degree,
                  double to = 170*utl::degree);


    double CombinedChi2(const double chi0, double &rp, double &t0,
                        const fevt::PCGFData::EPCGFStatus status,
                        const bool verbose = false);

    bool Prescan(const double step, double &from, double &to);

    bool Underground(const double chi0, const fevt::Eye &eye);

    void FillParams(fevt::Eye& eyeCopy,
                    const double chi0, const double chi0Err, const double rp, const double rpErr,
                    const double t0, const double t0Err) const;

    void FillParams(fevt::Eye& eyeCopy, const double chi0, const double rp, const double t0) const;

    bool AdjustGeometry(fevt::Eye& eyeCopy) const;

  private:
    std::vector<fevt::PCGFData> fPCGFData;
  };

}


#endif