RdPlaneFit.h

Go to the documentation of this file.

#ifndef _RdPlaneFit_RdPlaneFit_h_
#define _RdPlaneFit_RdPlaneFit_h_

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

/*namespace evt {
 class Event;
 }*/

namespace revt {
  class REvent;
}

// you should assign your module to a namespace

namespace RdPlaneFit {

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

    FitParameters() :
        u(0), v(0), w(0), ct0(0), sigmaU2(0), sigmaV2(0), sigmaUV(0), sigmact02(0), stage(0)
    {
    }
  };

  class RdPlaneFit: public fwk::VModule {

  public:

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

  private:
    fwk::VModule::ResultFlag LinearFit(FitParameters& fit, const bool reuseFit = false) const;
    
    fwk::VModule::ResultFlag PlaneFit3DDriver(FitParameters& fit) const;
    
    static void PlaneFit3DFnc(int& nPar, double* const grad, double& value, double* const par,
        const int flag);
    
    static void PlaneFit3DHorizonFnc(int& nPar, double* const grad, double& value,
        double* const par, const int flag);
    
    void CalculateTimeResidual(const utl::Vector& axis, const double ct0, double& chi2) const;

    // xml settings
    bool fMinuitOutput;
    int fMinNumberOfStations;
    bool fAllowUnphysicalCosines;
    bool fLockParameterStorage;

    // minuit globals
    static revt::REvent* fgCurrentREvent;
    static utl::Point fgBarycenter;
    static utl::TimeStamp fgBaryTime;
    static utl::Point fgCoordinateOrigin;
    static utl::CoordinateSystemPtr fgLocalCS;
    static double fgMeanEventTime; // this is the weighted mean of the station signal times with respect to the event time

    REGISTER_MODULE("RdPlaneFit", RdPlaneFit);

  };

}

#endif