FdEnergyFinder.h

Go to the documentation of this file.

#ifndef _FdEnergyFinderOG_FdEnergyFinder_h_
#define _FdEnergyFinderOG_FdEnergyFinder_h_

#include <fwk/VModule.h>
#include <evt/GaisserHillas4Parameter.h>
#include <vector>


namespace utl{
  class TabulatedFunctionErrors;
}


namespace FdEnergyFinderOG {

  class FdEnergyFinder :public fwk::VModule {
    
  public: 

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

  private:

    void FitProfile();
    
    void FindEmEnergy();
    
    void FindEmEnergyError();
    
    void FindTotalEnergy();
    
    // void FindTotalEnergy(double Eem);
    
    // double PolinomialInterpolation(double xa[], double ya[], 
    //                               int n, double exp);

    static void MinuitFitFunction(int& npar, double *gin, double& f, 
                                  double *par, int iflag);

    static const utl::TabulatedFunctionErrors* fLongProfile;
    
    evt::GaisserHillas4Parameter* fGHParameters;

    bool  fMinuitFailed;

    double fNmaxGH;
    double fXmaxGH;
    double fX0GH;
    double fEemGH;
    double fEtotalGH;
    double fNmaxGH_error;
    double fXmaxGH_error;
    double fX0GH_error;
    double fEemGH_error;
    static double fChisqGH;
    static double fNdof;
    double fEnergyCutoff;
    unsigned int fFixX0;
    unsigned int fComposition;
    unsigned int fUnseenEnergyModel;


    REGISTER_MODULE("FdEnergyFinderOG",FdEnergyFinder);

  }; // FdEnergyFinder


} // FdEnergyFinderOG


#endif // _FdEnergyFinderOG_FdEnergyFinder_h_

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