#ifndef ALIEMCALQADATAMAKERREC_H
#define ALIEMCALQADATAMAKERREC_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/*
  Produces the data needed to calculate the quality assurance. 
  All data must be mergeable objects.

  Based on PHOS code written by
  Y. Schutz CERN July 2007
 
 Created one histogram for QA shifter;-- Yaxian Mao: 11/2009
 The idea:average counts for all the towers should be flat 
 Change all existing histograms as experts
 
 Change histograms for DQM shifter: --  Yaxian Mao 04/2010
 Calcuate the amplitude ratio from current run and the LED reference, for QAChecker use
 Also calculate the ratio of amplitude from LED Monitor system (current/Reference), to check LED system  
 
*/

// --- ROOT system ---
class TH1F ; 
class TH1I ;
class TH2F ;
class TH2 ; 
class TLine ;
class TText ;
class TProfile ;
class TObjArray ; 

// --- AliRoot header files ---
#include "AliQADataMakerRec.h"
class AliCaloRawAnalyzer;
class AliEMCALGeometry;
#include "AliCaloConstants.h"

class AliEMCALQADataMakerRec: public AliQADataMakerRec {

public:
  
  //Histograms for Raw data control
  enum HRawType_t { 
    // first normal Low Gain and High Gain info
    kNsmodLG,kNsmodHG,kTimeLG,kTimeHG,
    kNtotLG,kNtotHG,kSigHG,kSigLG,
    kPedLG,kPedHG,
    k2DRatioAmp,kRatioDist, kLEDMonRatio, kLEDMonRatioDist,
    // then TRU info
    kNsmodTRU,
    kSigTRU,kNtotTRU,
    kNL0TRU, kTimeL0TRU,
		kNL0FirstTRU, kTimeL0FirstTRU,
    // and also LED Mon info
    kNsmodLGLEDMon,kNsmodHGLEDMon,kTimeLGLEDMon,kTimeHGLEDMon,
    kSigLGLEDMon,kSigHGLEDMon,kNtotLGLEDMon,kNtotHGLEDMon,
    kPedLGLEDMon,kPedHGLEDMon,
		//and STU info
		kAmpL1, kGL1, kJL1,
		kGL1V0, kJL1V0, kSTUTRU  
  } ;

  //Histograms for RecPoints  control
  enum HRPType_t {kRecPE,kRecPM,kRecPDigM};

  //Histograms for ESDs  control
  enum HESDType_t {kESDCaloClusE,kESDCaloClusM,kESDCaloCellA,kESDCaloCellM} ;
  
public:
  
  AliEMCALQADataMakerRec(Int_t fitAlgo = Algo::kStandard) ;  // ctor
  // Default fitter is kStandard=0 (see AliCaloConstants)
  
  AliEMCALQADataMakerRec(const AliEMCALQADataMakerRec& qadm) ;   
  AliEMCALQADataMakerRec& operator = (const AliEMCALQADataMakerRec& qadm) ;
  virtual ~AliEMCALQADataMakerRec() {;} // dtor

  Int_t GetFittingAlgorithm() const {return fFittingAlgorithm;}
  void  SetFittingAlgorithm(Int_t val);
  
  AliCaloRawAnalyzer *GetRawAnalyzer()    const {return fRawAnalyzer   ;}
  AliCaloRawAnalyzer *GetRawAnalyzerTRU() const {return fRawAnalyzerTRU;}

  void SetSuperModules(int i) {fSuperModules = i;}; //The number of SuperModules
  int GetSuperModules() const {return fSuperModules;}; //The number of SuperModules

  // for pedestal calculation with raw data
  void SetFirstPedestalSample(int i)  {fFirstPedestalSample = i;}          // first sample
  int  GetFirstPedestalSample() const {return fFirstPedestalSample;}       // first sample
  void SetLastPedestalSample(int i)   {fLastPedestalSample = i;}           // last sample
  int  GetLastPedestalSample() const  {return fLastPedestalSample;}        // last sample
  
  void SetFirstPedestalSampleTRU(int i)  {fFirstPedestalSampleTRU = i;}    // first sample, TRU
  int  GetFirstPedestalSampleTRU() const {return fFirstPedestalSampleTRU;} // first sample, TRU
  void SetLastPedestalSampleTRU(int i)   {fLastPedestalSampleTRU = i;}     // last sample, TRU
  int  GetLastPedestalSampleTRU() const  {return fLastPedestalSampleTRU;}  // last sample, TRU
  
  // for selection of interesting signal (max-min) range 
  // Low Gain channels
  void SetMinSignalLG(int i)  {fMinSignalLG = i;}
  int  GetMinSignalLG() const {return fMinSignalLG;}
  void SetMaxSignalLG(int i)  {fMaxSignalLG = i;}
  int  GetMaxSignalLG() const {return fMaxSignalLG;}
  
  // High Gain channels
  void SetMinSignalHG(int i)  {fMinSignalHG = i;}
  int  GetMinSignalHG() const {return fMinSignalHG;}
  void SetMaxSignalHG(int i)  {fMaxSignalHG = i;}
  int  GetMaxSignalHG() const {return fMaxSignalHG;}
  
  // TRU channels
  void SetMinSignalTRU(int i) {fMinSignalTRU = i;}
  int GetMinSignalTRU() const {return fMinSignalTRU;}
  void SetMaxSignalTRU(int i) {fMaxSignalTRU = i;}
  int GetMaxSignalTRU() const {return fMaxSignalTRU;}
  
  // LEDMon channels
  void SetMinSignalLGLEDMon(int i) {fMinSignalLGLEDMon = i;}
  int GetMinSignalLGLEDMon() const {return fMinSignalLGLEDMon;}
  void SetMaxSignalLGLEDMon(int i) {fMaxSignalLGLEDMon = i;}
  int GetMaxSignalLGLEDMon() const {return fMaxSignalLGLEDMon;}
  void SetMinSignalHGLEDMon(int i) {fMinSignalHGLEDMon = i;}
  int GetMinSignalHGLEDMon() const {return fMinSignalHGLEDMon;}
  void SetMaxSignalHGLEDMon(int i) {fMaxSignalHGLEDMon = i;}
  int GetMaxSignalHGLEDMon() const {return fMaxSignalHGLEDMon;}

  virtual void   EndOfDetectorCycle(AliQAv1::TASKINDEX_t, TObjArray ** list) ;
  void           GetCalibRefFromOCDB() ;
  void		       GetTruChannelPosition( Int_t &globRow, Int_t &globColumn, Int_t module, Int_t ddl, Int_t branch, Int_t column ) const;
  virtual void   InitESDs() ; 
  virtual void   InitDigits() ; 
  virtual void   InitRecPoints() ; 
  virtual void   InitRaws() ; 
  virtual void   MakeESDs(AliESDEvent * esd) ;
  virtual void   MakeDigits() ;
  virtual void   MakeDigits(TTree * digTree) ; 
  virtual void   MakeRecPoints(TTree * recpoTree) ; 
  virtual void   MakeRaws(AliRawReader* rawReader) ;
  virtual void	 MakeRawsSTU(AliRawReader* rawReader);
  virtual void   StartOfDetectorCycle() ; 

private:
  
  void ConvertProfile2H(TProfile * p, TH2 * histo) ; //change the profile plot to a 2D histogram
  
  Int_t fFittingAlgorithm;              // select the fitting algorithm, only used in copy ctor ...

  AliCaloRawAnalyzer *fRawAnalyzer;     // for signal fitting
  AliCaloRawAnalyzer *fRawAnalyzerTRU;  // for TRU signal fitting
  AliEMCALGeometry   *fGeom;            // EMCAL geometry, needed for STU decoding

  int fSuperModules;           // The number of SuperModules activated
  int fFirstPedestalSample;    // first sample for pedestal calculation, in bunch
  int fLastPedestalSample;     // last sample for pedestal calculation, in bunch
  int fFirstPedestalSampleTRU; // first sample for pedestal calculation, in bunch
  int fLastPedestalSampleTRU;  // last sample for pedestal calculation, in bunch
  
  int fMinSignalLG;       // minimum signal, for Low Gain channels
  int fMaxSignalLG;       // maximum signal, for Low Gain channels
  int fMinSignalHG;       // minimum signal, for High Gain channels
  int fMaxSignalHG;       // maximum signal, for High Gain channels
  int fMinSignalTRU;      // minimum signal, for TRU channels
  int fMaxSignalTRU;      // maximum signal, for TRU channels
  int fMinSignalLGLEDMon; // minimum signal, for LEDMon channels, low gain
  int fMaxSignalLGLEDMon; // maximum signal, for LEDMon channels, low gain
  int fMinSignalHGLEDMon; // minimum signal, for LEDMon channels, high gain
  int fMaxSignalHGLEDMon; // maximum signal, for LEDMon channels, high gain

  TProfile * fCalibRefHistoPro ; // Profile reference histogram from LED run
  TH2F     * fCalibRefHistoH2F ; // H2F reference histogram from LED run
  TProfile * fLEDMonRefHistoPro; // Profile reference histogram from LED monitor
  TH2F     * fHighEmcHistoH2F  ; // H2F reference histogram from LED run
  
//  TText **    fTextSM        ; //! Text info for each SM  
//  TLine *     fLineCol       ; //! line to distinguish the different SM side: A side and C side
//  TLine *     fLineRow       ; //! line to distinguish the different SM sector 0 and 1 

  ClassDef(AliEMCALQADataMakerRec,6)  // description

};

#endif // AliEMCALQADATAMAKERREC_H