#ifndef  __JRECONSTRUCTION__JSHOWERFITPARAMETERS__
#define  __JRECONSTRUCTION__JSHOWERFITPARAMETERS__

#include <TROOT.h>
#include <TObject.h>
#include "JFit/JMEstimator.hh"

/**
 * \author adomi, gmaggi
 */
namespace JRECONSTRUCTION {

  /**
   * Data structure for fit parameters.
   */
  struct JShowerFitParameters_t:
    public TObject
  {
    /**
     * Default constructor.
     */
    JShowerFitParameters_t()
    {
      this->reset();
    }

    /**
     * Virtual destructor.
     */
    virtual ~JShowerFitParameters_t(){}
    
    /**
     * Reset fit parameters.
     */
    void reset()
    {
      numberOfPrefits             =   4;
      TMax_ns                     =  30;
      TMin_ns                     = -30;
      roadWidth_m                 =  80;
      R_Hz                        =  10.0e3;
      mestimator                  =  JFIT::EM_NULL;
      fit_step                    =  0.05;
    }

    /**
     * Equality
     *
     * \param  parameters       fit parameters
     * \return                  true if equals; else false
     */
    bool equals(const JShowerFitParameters_t& parameters) const
    {
      return (this->TMax_ns                  == parameters.TMax_ns                  &&         
	      this->TMin_ns                  == parameters.TMin_ns                  &&         
	      this->numberOfPrefits          == parameters.numberOfPrefits          &&         
	      this->roadWidth_m              == parameters.roadWidth_m              &&
	      this->R_Hz                     == parameters.R_Hz                     && 
	      this->mestimator               == parameters.mestimator               &&
	      this->fit_step                 == parameters.fit_step);
    }

    ClassDef(JShowerFitParameters_t, 1);


    size_t  numberOfPrefits;     ///< number of prefits                                                                    
    double  TMax_ns;             ///< maximum time for local coincidences [ns]                                               
    double  TMin_ns;             ///< minimum time for local coincidences [ns]                                             
    double  roadWidth_m;         ///< road width [m]     
    double  R_Hz;                ///< default rate [Hz]                                                                       
    int     mestimator;          ///< M-estimator     
    double  fit_step;            ///< simplex initial step  
  };  
}

/**
  * Read fit parameters from input.
  *
  * \param  in           input stream
  * \param  parameters   fit parameters
  * \return              input stream
  */
 std::istream& operator>>(std::istream& in, JRECONSTRUCTION::JShowerFitParameters_t& parameters);
   
 /**
  * Write fit parameters to output.
  *
  * \param  out          output stream
  * \param  parameters   fit parameters
  * \return              output stream
  */
 std::ostream& operator<<(std::ostream& out,  const JRECONSTRUCTION::JShowerFitParameters_t& parameters);

#endif