#include <Rtypes.h>
#ifndef __eeSinglePeak_hh__
#define __eeSinglePeak_hh__

#include <vector>
struct eeSinglePeakParam_t {
  //-- Fit parameters --
  Double_t yield;      // Integrated yield
  Double_t mean;       // Mean
  Double_t sigma;      // Sigma
  Double_t ryield;     // Relative yield for 2nd gaussian
  Double_t rsigma;     // Relative sigma for 2nd gaussian
  //-- Fit errors --
  Double_t eyield;     // Integrated yield
  Double_t emean;      // Mean
  Double_t esigma;     // Sigma
  Double_t eryield;    // Relative yield for 2nd gaussian
  Double_t ersigma;    // Relative sigma for 2nd gaussian
  //-- Limits --
  Double_t min_yield;
  Double_t min;
  Double_t min_sigma;
  Double_t min_ryield;
  Double_t min_rsigma;
  Double_t max_yield;
  Double_t max;
  Double_t max_sigma;
  Double_t max_ryield;
  Double_t max_rsigma;  
  //-- Fit statistics --
  Double_t chiSquared; // Chi^2 of the fit
  Int_t    ndf;        // Degrees of freedom
  //-- Seed tower of cluster + plane of fit --
  Int_t sector;        // Sector [0..11]
  Int_t subsector;     // Subsector [0..4]
  Int_t etabin;        // Eta bin (aka tile number) [0..11]
  Int_t plane;         // Plane (U,V) [0..1]
};

//-- a vector of single peak parameters --
typedef std::vector<eeSinglePeakParam_t> eeSinglePeakParamVec_t;

Double_t eeSinglePeak ( Double_t *x, Double_t *p );
Double_t eeSinglePeak ( Double_t *X, eeSinglePeakParam_t *p);

#endif