#ifndef MARS_MCalibrationChargeCalc
#define MARS_MCalibrationChargeCalc

/////////////////////////////////////////////////////////////////////////////
//                                                                         //
// MCalibrationChargeCalc                                                   //
//                                                                         //
// Integrates the time slices of the all pixels of a calibration event     //
// and substract the pedestal value                                        //
//                                                                         //
/////////////////////////////////////////////////////////////////////////////

#ifndef MARS_MTask
#include "MTask.h"
#endif

#ifndef MARS_MBadPixelsPix
#include "MBadPixelsPix.h"
#endif

class MRawEvtData;
class MRawRunHeader;
class MPedestalCam;
class MPedestalPix;
class MCalibrationChargePINDiode;
class MCalibrationChargeBlindPix;
class MCalibrationChargePix;
class MCalibrationChargeCam;
class MCalibrationQECam;
class MGeomCam;
class MExtractedSignalCam;
class MBadPixelsCam;
class MBadPixelsPix;
class MTime;
class MCalibrationChargeCalc : public MTask
{
private:

  static const Float_t fgChargeLimit;       // The default for fChargeLimit
  static const Float_t fgChargeErrLimit;    // The default for fChargeErrLimit
  static const Float_t fgChargeRelErrLimit; // The default for fChargeRelErrLimit
  static const Float_t fgTimeLowerLimit;    // The default for fTimeLowerLimit
  static const Float_t fgTimeUpperLimit;    // The default for fTimeUpperLimit
  
  Float_t fChargeLimit;          // The limit (in units of PedRMS) for acceptance of the fitted mean charge
  Float_t fChargeErrLimit;       // The limit (in units of PedRMS) for acceptance of the fitted charge sigma square
  Float_t fChargeRelErrLimit;    // The limit (in units of Sigma of fitted charge) for acceptance of the fitted mean  
  Float_t fTimeLowerLimit;       // The limit (in units of FADC slices) for distance to first signal extraction slice
  Float_t fTimeUpperLimit;       // The limit (in units of FADC slices) for distance to last signal extraction slice

  MPedestalCam               *fPedestals;     //! Pedestals of all pixels in the camera
  MCalibrationChargeCam      *fCam;           // Calibration events of all pixels in the camera
  MCalibrationQECam          *fQECam;         // Derived Quantum efficiency of all pixels in the camera  
  MRawEvtData                *fRawEvt;        //! raw event data (time slices)
  MRawRunHeader              *fRunHeader;     //! RunHeader information
  MGeomCam                   *fGeom;          //! Geometry information
  MBadPixelsCam              *fBadPixels;     //! Bad Pixels information
  
  MTime                      *fEvtTime;       //! Time of the event

  MExtractedSignalCam        *fSignals;       // Extracted signal of all pixels in the camera
  MCalibrationChargePINDiode *fPINDiode;      // Calibration results of the PIN Diode
  MCalibrationChargeBlindPix *fBlindPixel;    // Calibration results of the Blind Pixel 

  Float_t  fNumHiGainSamples; 
  Float_t  fNumLoGainSamples; 
  Float_t  fSqrtHiGainSamples;
  Float_t  fSqrtLoGainSamples; 
  
  Float_t fConversionHiLo;
  Int_t   fFlags;                              // Flag for the fits used

  enum  { kUseQualityChecks,
          kHiLoGainCalibration };
  
  Int_t  PreProcess(MParList *pList);
  Bool_t ReInit(MParList *pList); 
  Int_t  Process();
  Int_t  PostProcess();

  void FinalizePedestals(const MPedestalPix &ped, MCalibrationChargePix &cal,
                         Float_t &avped, Float_t &avrms);
  void FinalizeAvPedestals(MCalibrationChargePix &cal, Float_t avped, Float_t avrms, Int_t avnum);
  Bool_t FinalizeCharges(MCalibrationChargePix &cal, MBadPixelsPix &bad);

  void PrintUnsuitable(MBadPixelsPix::UnsuitableType_t typ, const char *text) const;   
  void PrintUncalibrated(MBadPixelsPix::UncalibratedType_t typ, const char *text) const;

public:

  MCalibrationChargeCalc(const char *name=NULL, const char *title=NULL);

  void Clear(const Option_t *o="");
  
  void SetChargeLimit    (   const Float_t f=fgChargeLimit       ) { fChargeLimit       = f; }
  void SetChargeErrLimit (   const Float_t f=fgChargeErrLimit    ) { fChargeErrLimit    = f; }
  void SetChargeRelErrLimit( const Float_t f=fgChargeRelErrLimit ) { fChargeRelErrLimit = f; }

  void SetTimeLowerLimit (   const Float_t f=fgTimeLowerLimit    ) { fTimeLowerLimit  = f;   }
  void SetTimeUpperLimit (   const Float_t f=fgTimeUpperLimit    ) { fTimeUpperLimit  = f;   }

  void SkipQualityChecks(Bool_t b=kTRUE)
      {b ? CLRBIT(fFlags, kUseQualityChecks)    : SETBIT(fFlags, kUseQualityChecks);}
  void SkipHiLoGainCalibration(Bool_t b=kTRUE)
      {b ? CLRBIT(fFlags, kHiLoGainCalibration) : SETBIT(fFlags, kHiLoGainCalibration);}

  ClassDef(MCalibrationChargeCalc, 1)   // Task to fill the Calibration Containers from raw data
};

#endif