MdGeometryFitter.h

Go to the documentation of this file.

#ifndef _MdGeometryFitterAG_MdGeometryFitter_h_
#define _MdGeometryFitterAG_MdGeometryFitter_h_

#include <string>

#include <TVector3.h>

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

#include <mevt/MEvent.h>
#include <evt/ShowerMRecData.h>

#include "GeometryFitter.h"


namespace MdGeometryFitterAG {

  class MdGeometryFitter : public fwk::VModule {

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

  private:

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

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

    // Parametrization of the radius of curvature of the spherical shower front
    // radius/m = par1 + par2 * (theta/rad)^2
    double fCurvatureRadiusParameter1, fCurvatureRadiusParameter2;

    // Extract from the MD event of the data used in the fit (t50, t50 error and counter positions)
    // The origins of position and time are set to the core position and time reconstructed by the SD respectively  
    void SetTimeData(const mevt::MEvent* mEvent, utl::Point sdCorePosition, utl::TimeStamp fgSdCoreTime, const TVector3& axis, GeometryFitter& fitter); 

    // Fill the MD event with the reconstructed data 
    void FillRecShower(evt::Event& event, const GeometryFitter& geometryFitter, utl::Point sdCorePosition, utl::TimeStamp sdCoreTime, TVector3 sdAxis);

    // Set time residual of counters
    void FillCounter(mevt::MEvent* mEvent, const GeometryFitter& fitter);

    REGISTER_MODULE("MdGeometryFitterAG", MdGeometryFitter);
  };

} // end namespace MdGeometryFitterAG


#endif

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