FdAxisFinder.h

Go to the documentation of this file.

#ifndef _FdAxisFinderOG_FdAxisFinder_h_
#define _FdAxisFinderOG_FdAxisFinder_h_

#include <fwk/VModule.h>
#include <utl/Vector.h>
#include <vector>

namespace fevt {
  class Eye;
  class Pixel;
}





namespace FdAxisFinderOG {

  class FdAxisFinder :public fwk::VModule {

  public:
    
    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& event);
    fwk::VModule::ResultFlag Finish();

  private:
    void AxisFirstGuess(fevt::Eye& eye);
    void FindAxis(fevt::Eye& eye);
    void FindAxisLaser(fevt::Eye& eye);
    void FillPoints(fevt::Eye& eye);
    void RemoveOutliers(fevt::Eye& eye);
    bool ReadmitPixel(fevt::Eye& eye);
    bool RejectPixel(fevt::Eye& eye);
    bool CheckForLasers(evt::Event& event);
    
    static void  MinuitFitFunc(int& npar, double *gin, double& f, 
                               double *par, int iflag);
    static void  MinuitFitFuncLaser(int& npar, double *gin, double& f, 
                               double *par, int iflag);

    static utl::Vector fSDP;
    static double fChi2;
    static unsigned int fNDof;

    double fChi0;
    double fChi0Error;
    double fT0;
    double fT0Error;
    static double fRp;  //static for laser reconstruction
    double fRpError;
    double fCovariance[3][3];

    double fT0First;
    double fChi0First;
    double fRpFirst;

    double fTolerance;

    double fSigmasLimit;
    double fSelectAndRecLasers;
    double   fMinimumTimeExtension;
    
    bool     fIsBlob;

    unsigned int fMinPixels;
    
    static int fDebuging;

    static fevt::Eye * fCurEye;
    static bool fIsCeleste;  //for laser reconstruction
    static bool fIsRamiro;   //for laser reconstruction
    bool fMinuitFailed;

    REGISTER_MODULE("FdAxisFinderOG", FdAxisFinder);
  };

} // FdAxisFinderOG


#endif // _FdAxisFinderOG_FdAxisFinder_h_

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: