FdSDPFinder.h

Go to the documentation of this file.

#ifndef _FdSDPFinderOG_FdSDPFinder_h_
#define _FdSDPFinderOG_FdSDPFinder_h_


namespace utl {class AxialVector;}
namespace fevt{
  class Eye;
  class Pixel;
}


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


namespace FdSDPFinderOG {

  class FdSDPFinder :public fwk::VModule{

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

  private:
    
    utl::AxialVector    FindSDPFirstGuess(fevt::Eye& eye);
    void RemoveNoise(fevt::Eye& eye);
    void ReadmitPixel(fevt::Eye& eye);
    void FindSDP(fevt::Eye& eye);
    void EstimateChi_i(fevt::Eye& eye);
    void RemoveOutliers(fevt::Eye& eye);
    bool IsIsolated(const fevt::Pixel&, const fevt::Eye& eye ) ;

    static void  MinuitFitFunc(int& npar, double *gin, double& f, 
                               double *par, int iflag);

    static fevt::Eye * fCurEye;
    
    utl::Vector      fSDP;
    utl::Vector      fSDPFirstGuess;

    double           fSDPThetaError;
    double           fSDPPhiError;
    double           fSDPCorrThetaPhi;
    static double    fChi2;
    static unsigned int fNDof;

    unsigned int     fMinPixels; 
    double           fMaxDistance;
    bool           fFirstGuessFromCharge;

    double           fIsolationLimit;
    bool             fGoodEvent;
    REGISTER_MODULE("FdSDPFinderOG",FdSDPFinder);
  };


} // FdSDPFinderOG


#endif // _FdSDPFinderOG_FdSDPFinder_h_

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