HybridGeometryFinderWG/HybridGeometryFinder.h

Go to the documentation of this file.

#ifndef _HybridGeometryFinderWG_HybridGeometryFinder_h_
#define _HybridGeometryFinderWG_HybridGeometryFinder_h_

#include <fwk/VModule.h>

#include <utl/AxialVector.h>
#include <utl/NumericalErrorPropagation.h>
#include <utl/Point.h>
#include <utl/Vector.h>

#include <map>
#include <vector>


namespace fevt {
  class Eye;
}

namespace sevt {
  class Station;
}

namespace HybridGeometryFinderWG {

  class HybridGeometryFinder : public fwk::VModule {

    public:

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

    private:

      int GetDataFromPreviousFit();

      void CalculateGeometryForTels();

      static std::vector<double> CalculateAxisForTel(const utl::Point, utl::Vector, const double, const int);

      void SelectPixels();

      bool SelectHybridStation(const int);

      bool FitAxis();

      static void MinuitFitFuncAxis(int&, double*, double&, double*, int);

      void RecalculateChiSquare();

      bool ReadmitPixel();

      bool RemovePixel();

      void StoreData();

    private:

      static int fgVerbosity;

      static double fgCoreHeight;

      static double fgSDPError;

      int fTimeFitModel;

      unsigned int fAxisMinPixel;

      double fMaxSDPDist;

      static double fgVerticalCurvature;

      double fMaxStationToSDPDist;

      double fMaxTimeResidue;

      static int fgLaserShots;

      int fFdOnly;

      int fPassThrough;

      evt::Event* fEvent;

      static fevt::Eye* fgCurEye;

      unsigned int fCurTelId;

      // The SDP stuff
      utl::AxialVector fSDP;
      //double fSDPThetaError; // unused. LN.
      //double fSDPPhiError; // unused. LN.
      //double fSDPCorrThetaPhi; // unused. LN.
      double fChiSqSDP;
      unsigned int fNDofSDP;
      std::map<int, double> fChiSqSDPPerTel;
      std::map<int, unsigned int> fNDofSDPPerTel;

      // The axis stuff
      double fChiZero;
      double fRp;
      double fTZero;
      utl::Vector fAxis;
      utl::Point fCore;
      double fAxisTheta;
      double fAxisThetaError;
      double fAxisPhi;
      double fAxisPhiError;
      double fNorthing;
      double fNorthingError;
      double fEasting;
      double fEastingError;
      double fTCore;
      double fTCoreError;
      double fCovMatrixAxis[5][5];
      double fChiSqTime;
      unsigned int fNDofTime;
      std::map<int, double> fChiSqTimePerTel;
      std::map<int, unsigned int> fNDofTimePerTel;
      double fChiSqAxis;
      static unsigned int fgNDofAxis;
      std::map<int, double> fChiSqAxisPerTel;
      std::map<int, unsigned int> fNDofAxisPerTel;

      // The SD stuff
      static std::vector<unsigned short int> fgStationList;
      static sevt::Station* fgFitStation;
      double fFitStationOffset;
      double fDeltaStationAxis;
      static int fgHybridFitMode;

      std::map<int, std::vector<double> > fTelParameters;

      REGISTER_MODULE("HybridGeometryFinderWG", HybridGeometryFinder);

      friend class AxisParameterCalculator;
  };


  class AxisParameterCalculator : public utl::FCNCalculator {

    public:

      AxisParameterCalculator(const HybridGeometryFinder* hybridFinder) {fHybridFinder = hybridFinder;}

      std::vector<double> operator() (const std::vector<double>&) const;

    private:

      const HybridGeometryFinder* fHybridFinder;

  };

} // HybridGeometryFinderWG


#endif // _HybridGeometryFinderWG_HybridGeometryFinder_h_

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