#ifndef  __JRECONSTRUCTION__JSHOWERPOSITIONFITPARAMETERS__
#define  __JRECONSTRUCTION__JSHOWERPOSITIONFITPARAMETERS__

#include <TROOT.h>
#include <TObject.h>

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

  struct JShowerPositionFitParameters_t:
    public TObject
  {
    /**
     * Default constructor.
     */
    JShowerPositionFitParameters_t()
    {
      this->reset();
    }

    /**
     * Virtual destructor.
     */
    virtual ~JShowerPositionFitParameters_t()
    {}
    
    /**
     * Reset fit parameters.
     */
    void reset()
    {
      numberOfPrefits  =     0;  // process all inputs from PointSimplex, which should give a max of 100 
      numberOfOutfits  =     5;
      TMax_ns          =   100.0;  
      TMin_ns          =  -100.0;   
      roadWidth_m      =    80.0;
      E_GeV            =    10;
      R_Hz             = 10000;   
      TTS_ns           =     2;
      VMax_npe         =    20.0;
      NMax             =  1000;
    }

    /**
     * Equality
     *
     * \param  parameters       fit parameters
     * \return                  true if equals; else false
     */
    bool equals(const JShowerPositionFitParameters_t& parameters) const
    {
      return (this->numberOfPrefits  == parameters.numberOfPrefits  &&         
	      this->numberOfOutfits  == parameters.numberOfOutfits  &&
	      this->TMax_ns          == parameters.TMax_ns          &&
	      this->TMin_ns          == parameters.TMin_ns          &&
	      this->ctMin            == parameters.ctMin            &&
	      this->roadWidth_m      == parameters.roadWidth_m      &&
	      this->E_GeV            == parameters.E_GeV            &&
	      this->R_Hz             == parameters.R_Hz             &&
	      this->TTS_ns           == parameters.TTS_ns           &&
	      this->VMax_npe         == parameters.VMax_npe         &&
	      this->NMax             == parameters.NMax);
    }

    ClassDef(JShowerPositionFitParameters_t, 1);

    size_t     numberOfPrefits;     ///< number of prefits 
    size_t     numberOfOutfits;     ///< number of fits to be saved in output 
    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     ctMin;               ///< minimal cosine space angle between PMT axes      
    double     E_GeV;               ///< energy [GeV]
    double     R_Hz;                ///< default rate [Hz]
    double     TTS_ns;              ///< transition-time spread [ns]
    double     VMax_npe;            ///< maximum number of of photo-electrons
    int        NMax;                ///< maximum number of iterations
  };
}

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

#endif