#ifndef  __JRECONSTRUCTION__JSHOWERDIRECTIONPREFITPARAMETERS__
#define  __JRECONSTRUCTION__JSHOWERDIRECTIONPREFITPARAMETERS__

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

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

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

    /**
     * Virtual destructor.
     */
    virtual ~JShowerDirectionPrefitParameters_t(){}
    
    /**
     * Reset fit parameters.
     */
    void reset()
    {
      numberOfPrefits             =    1;
      TMax_ns                     =   30;
      TMin_ns                     =  -30;
      DMax_m                      =   80;
      scanAngle_deg               =   10;
      R_Hz                        = 10.0e3;
      VMax_npe                    =   20.0;
      Emin_GeV                    = 1;
      Emax_GeV                    = 1000;
      En                          = 50;
    }

    /**
     * Equality
     *
     * \param  parameters       fit parameters
     * \return                  true if equals; else false
     */
    bool equals(const JShowerDirectionPrefitParameters_t& parameters) const
    {
      return (this->TMax_ns            == parameters.TMax_ns                  &&         
	      this->TMin_ns                  == parameters.TMin_ns                  &&         
 	      this->numberOfPrefits          == parameters.numberOfPrefits          &&                 
	      this->DMax_m                   == parameters.DMax_m                   &&
	      this->R_Hz                     == parameters.R_Hz                     && 
 	      this->scanAngle_deg            == parameters.scanAngle_deg            &&
	      this->VMax_npe                 == parameters.VMax_npe                 &&
              this->Emin_GeV                 == parameters.Emin_GeV                 &&
              this->Emax_GeV                 == parameters.Emax_GeV                 &&
              this->En                       == parameters.En                         ); 
    }

    ClassDef(JShowerDirectionPrefitParameters_t, 2);

    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 DMax_m;              ///< maximal distance to optical module [m]
    double R_Hz;                ///< default rate [Hz]                                                                                                                                        
    double VMax_npe;            ///< maximum number of of photo-electrons                                                   
    double scanAngle_deg;       ///< scanning angle step in [deg]
    double Emin_GeV;            ///< minimum energy to scan
    double Emax_GeV;            ///< maximum energy to scan
    int    En;                  ///< number of points to scan in energy range
  };   
}

/**
  * Read fit parameters from input.
  *
  * \param  in           input stream
  * \param  parameters   fit parameters
  * \return              input stream
  */
 std::istream& operator>>(std::istream& in, JRECONSTRUCTION::JShowerDirectionPrefitParameters_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::JShowerDirectionPrefitParameters_t& parameters);

#endif