SdPlaneFit.h

Go to the documentation of this file.

#ifndef _SdPlaneFitOG_SdPlaneFit_h_
#define _SdPlaneFitOG_SdPlaneFit_h_

#include <fwk/VModule.h>
#include <utl/Vector.h>
#include <utl/Point.h>
#include <utl/TimeStamp.h>
#include <utl/CoordinateSystem.h>


namespace sevt {
  class SEvent;
}

namespace SdPlaneFitOG {

  struct FitParameters {
    double u = 0;
    double v = 0;
    double w = 0;
    double ct0 = 0;
    double sigmaU2 = 0;
    double sigmaV2 = 0;
    double sigmaUV = 0;
    double sigmact02 = 0;
    double stage = 0;
  };


  class SdPlaneFit : public fwk::VModule {

  public:
    fwk::VModule::ResultFlag Init() override;
    fwk::VModule::ResultFlag Run(evt::Event& event) override;
    fwk::VModule::ResultFlag Finish() override { return eSuccess; }

  private:
    fwk::VModule::ResultFlag LinearFit(FitParameters& fit,
                                       const bool reuseFit = false) const;
    fwk::VModule::ResultFlag PlaneFit3DDriver(FitParameters& fit) const;
    void CalculateTimeResidual3D(const utl::Vector& axis, const double ct0,
                                 double& mean, double& spread, double& chi2) const;

    static void PlaneFit3DFnc(int& nPar, double* const grad, double& value,
                              double* const par, const int flag);

    enum InfoLevel {
      eNone = 0,
      eFinal = 1,
      eIntermediate = 2,
      eObscure = 3,
      eMinuit = 4
    };

    // xml settings
    int fInfoLevel = 0;
    bool fMinuitOutput = false;
    int fMinNumberOfStations = 0;

    // minuit globals
    static sevt::SEvent* fgCurrentSEvent;
    static utl::Point fgBarycenter;
    static utl::TimeStamp fgBaryTime;
    static utl::CoordinateSystemPtr fgBarycenterCS;

    REGISTER_MODULE("SdPlaneFitOG", SdPlaneFit);
  };

}


#endif