RdHyperbolicWavefrontFit.h

Go to the documentation of this file.

#ifndef _RdHyperbolicWavefrontFit_h_
#define _RdHyperbolicWavefrontFit_h_

#include <fwk/VModule.h>
#include <revt/Station.h>
#include <revt/Channel.h>
#include <string>
#include <vector>
#include <utl/CoordinateSystem.h>
#include <utl/Vector.h>
#include <utl/Point.h>


namespace evt {
  class Event;
}

namespace revt {
  class Channel;
}

namespace RdHyperbolicWavefrontFit {

  struct FitParameters {

    // !!! all angles in rad !!!
    double u = 0;
    double uError = 0;
    double v = 0;
    double vError = 0;

    double gamma = 0;
    double gammaError = 0;

    double t0 = 0;
    double t0Error = 0;

    int status = 0;
    int nCalls = 0;
    double minChi = 0;
    int NDF = 0;
    double recStage = 0; // RecStage defined in trunk/Framework/Shower/ShowerRRecData.cc

  };


  class RdHyperbolicWavefrontFit : public fwk::VModule {

  public:
    RdHyperbolicWavefrontFit() { }
    virtual ~RdHyperbolicWavefrontFit() { }

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

  private:

    //the hyperbolic wavefront function
    static void RdHyperbolicWavefrontFnc(int& nPar, double* const grad, double& value, double* const par,
                                         const int flag);

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

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

    // minuit globals
    static revt::REvent* fgCurrentREvent;
    static evt::Event fgCurrentEvent;
    static utl::CoordinateSystemPtr fgLocalCS;
    static utl::Vector gSAxis;
    static utl::Point gCore;
    static double gTimeOffset;
    static unsigned int nPar;
    static double gTime0;
    static unsigned int gStation0Id;
    static double gMaxAmplitude;

    // xml settings
    std::string fUsedCore;
    bool fDoCoreSearch = false;

    // Calculation of the time residual of the individual station.
    // Also fills them in the StationRecData
    void CalculateTimeResidual(const utl::Vector& axis, const double gamma, const double t0) const;

    REGISTER_MODULE("RdHyperbolicWavefrontFit", RdHyperbolicWavefrontFit);

  };

}


#endif