RdScintPlaneFit.h

Go to the documentation of this file.

#ifndef _RdScintPlaneFit_h_
#define _RdScintPlaneFit_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 RdScintPlaneFit {

  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 RdScintPlaneFit: public fwk::VModule {

  public:
    // If you define a constructor, you MUST provide an implementation
    // for it, otherwise the REGISTER_MODULE macro will not register
    // your module. (Here, implementation is in RdPlaneFit.cc)
    // RdPlaneFit();
    //virtual ~RdPlaneFit();

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

  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;

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

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

    // 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

    // The REGISTER_MODULE macro makes the framework aware of your module.
    // You need to put this macro at the very end of your module class definition.
    // Give your module a name (first argument) and specify its 
    // class name (second argument).  In the ModuleSequenceExample.xml file,
    // you can see this module being refered to by the name we give it here.

  REGISTER_MODULE("RdScintPlaneFit", RdScintPlaneFit);

  };

}

#endif

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