G/PCGFitter.h

Go to the documentation of this file.

#ifndef _FdProfileConstrainedGeometryFitPG_PCGFitter_h_
#define _FdProfileConstrainedGeometryFitPG_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 FdProfileConstrainedGeometryFitPG {

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

    void Init();

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

    bool Run(fevt::Eye &eye);

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

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

    FdApertureLightFinderKG::FdApertureLightFinder fApertureLightFinder;
    ApertureLight fApLight;
    EApLightMethod fApLightMethod;
    bool fCheckUnderground;
    bool fPrescan;
    bool fScanOnly;
    double fScanStep;
    double fScanStart;
    double fScanStop;
    bool fUseLightFlux;
    bool fDelZeroFlux;
    bool fRealAtm;
    bool fDeex;
    bool fEmissionPointCorrection;
    unsigned int fTelId;
    bool fPrecise;
    bool fSkipNegativeT0;

    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,
                  int level = 0);


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

    bool Prescan(const double step, double &from, double &to);
    
    bool Underground(const double chi0, const double rp, 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