[u/mrichter/AliRoot.git] / EMCAL / AliEMCALQADataMakerRec.h

#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 ; 

// --- Standard library ---

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

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:
  enum fitAlgorithm {kFastFit=1, kNeuralNet = 2, kLMS = 4, kPeakFinder = 5, kCrude = 6};
  AliEMCALQADataMakerRec(fitAlgorithm fitAlgo = kNeuralNet) ;          // ctor
 
  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;}
  AliCaloRawAnalyzerKStandard *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

 AliCaloRawAnalyzer *fRawAnalyzer;    // for signal fitting
 AliCaloRawAnalyzerKStandard *fRawAnalyzerTRU;    // for signal fitting, for TRU
 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,5)  // description 

};

#endif // AliEMCALQADATAMAKERREC_H
Commit	Line	Data
	1	#ifndef ALIEMCALQADATAMAKERREC_H
	2	#define ALIEMCALQADATAMAKERREC_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/*
	7	Produces the data needed to calculate the quality assurance.
	8	All data must be mergeable objects.
	9
	10	Based on PHOS code written by
	11	Y. Schutz CERN July 2007
	12
	13	Created one histogram for QA shifter;-- Yaxian Mao: 11/2009
	14	The idea:average counts for all the towers should be flat
	15	Change all existing histograms as experts
	16
	17	Change histograms for DQM shifter: -- Yaxian Mao 04/2010
	18	Calcuate the amplitude ratio from current run and the LED reference, for QAChecker use
	19	Also calculate the ratio of amplitude from LED Monitor system (current/Reference), to check LED system
	20
	21	*/
	22
	23
	24	// --- ROOT system ---
	25	class TH1F ;
	26	class TH1I ;
	27	class TH2F ;
	28	class TH2 ;
	29	class TLine ;
	30	class TText ;
	31	class TProfile ;
	32	class TObjArray ;
	33
	34	// --- Standard library ---
	35
	36	// --- AliRoot header files ---
	37	#include "AliQADataMakerRec.h"
	38	class AliCaloRawAnalyzer;
	39	class AliCaloRawAnalyzerKStandard;
	40	//class AliEMCALGeoParams;
	41	class AliEMCALGeometry;
	42
	43	class AliEMCALQADataMakerRec: public AliQADataMakerRec {
	44
	45	public:
	46	//Histograms for Raw data control
	47	enum HRawType_t {
	48	// first normal Low Gain and High Gain info
	49	kNsmodLG,kNsmodHG,kTimeLG,kTimeHG,
	50	kNtotLG,kNtotHG,kSigHG,kSigLG,
	51	kPedLG,kPedHG,
	52	k2DRatioAmp,kRatioDist, kLEDMonRatio, kLEDMonRatioDist,
	53	// then TRU info
	54	kNsmodTRU,
	55	kSigTRU,kNtotTRU,
	56	kNL0TRU, kTimeL0TRU,
	57	kNL0FirstTRU, kTimeL0FirstTRU,
	58	// and also LED Mon info
	59	kNsmodLGLEDMon,kNsmodHGLEDMon,kTimeLGLEDMon,kTimeHGLEDMon,
	60	kSigLGLEDMon,kSigHGLEDMon,kNtotLGLEDMon,kNtotHGLEDMon,
	61	kPedLGLEDMon,kPedHGLEDMon,
	62	//and STU info
	63	kAmpL1, kGL1, kJL1,
	64	kGL1V0, kJL1V0, kSTUTRU
	65	} ;
	66
	67	//Histograms for RecPoints control
	68	enum HRPType_t {kRecPE,kRecPM,kRecPDigM};
	69
	70	//Histograms for ESDs control
	71	enum HESDType_t {kESDCaloClusE,kESDCaloClusM,kESDCaloCellA,kESDCaloCellM} ;
	72
	73
	74	public:
	75	enum fitAlgorithm {kFastFit=1, kNeuralNet = 2, kLMS = 4, kPeakFinder = 5, kCrude = 6};
	76	AliEMCALQADataMakerRec(fitAlgorithm fitAlgo = kNeuralNet) ; // ctor
	77
	78	AliEMCALQADataMakerRec(const AliEMCALQADataMakerRec& qadm) ;
	79	AliEMCALQADataMakerRec& operator = (const AliEMCALQADataMakerRec& qadm) ;
	80	virtual ~AliEMCALQADataMakerRec() {;} // dtor
	81
	82	Int_t GetFittingAlgorithm() const {return fFittingAlgorithm; }
	83	void SetFittingAlgorithm(Int_t val);
	84	AliCaloRawAnalyzer *GetRawAnalyzer() const { return fRawAnalyzer;}
	85	AliCaloRawAnalyzerKStandard *GetRawAnalyzerTRU() const { return fRawAnalyzerTRU;}
	86
	87	void SetSuperModules(int i) {fSuperModules = i;}; //The number of SuperModules
	88	int GetSuperModules() const {return fSuperModules;}; //The number of SuperModules
	89
	90	// for pedestal calculation with raw data
	91	void SetFirstPedestalSample(int i) {fFirstPedestalSample = i;}; // first sample
	92	int GetFirstPedestalSample() const {return fFirstPedestalSample;}; // first sample
	93	void SetLastPedestalSample(int i) {fLastPedestalSample = i;}; // last sample
	94	int GetLastPedestalSample() const {return fLastPedestalSample;}; // last sample
	95	void SetFirstPedestalSampleTRU(int i) {fFirstPedestalSampleTRU = i;}; // first sample, TRU
	96	int GetFirstPedestalSampleTRU() const {return fFirstPedestalSampleTRU;}; // first sample, TRU
	97	void SetLastPedestalSampleTRU(int i) {fLastPedestalSampleTRU = i;}; // last sample, TRU
	98	int GetLastPedestalSampleTRU() const {return fLastPedestalSampleTRU;}; // last sample, TRU
	99
	100	// for selection of interesting signal (max-min) range
	101	// Low Gain channels
	102	void SetMinSignalLG(int i) {fMinSignalLG = i;};
	103	int GetMinSignalLG() const {return fMinSignalLG;};
	104	void SetMaxSignalLG(int i) {fMaxSignalLG = i;};
	105	int GetMaxSignalLG() const {return fMaxSignalLG;};
	106	// High Gain channels
	107	void SetMinSignalHG(int i) {fMinSignalHG = i;};
	108	int GetMinSignalHG() const {return fMinSignalHG;};
	109	void SetMaxSignalHG(int i) {fMaxSignalHG = i;};
	110	int GetMaxSignalHG() const {return fMaxSignalHG;};
	111	// TRU channels
	112	void SetMinSignalTRU(int i) {fMinSignalTRU = i;};
	113	int GetMinSignalTRU() const {return fMinSignalTRU;};
	114	void SetMaxSignalTRU(int i) {fMaxSignalTRU = i;};
	115	int GetMaxSignalTRU() const {return fMaxSignalTRU;};
	116	// LEDMon channels
	117	void SetMinSignalLGLEDMon(int i) {fMinSignalLGLEDMon = i;};
	118	int GetMinSignalLGLEDMon() const {return fMinSignalLGLEDMon;};
	119	void SetMaxSignalLGLEDMon(int i) {fMaxSignalLGLEDMon = i;};
	120	int GetMaxSignalLGLEDMon() const {return fMaxSignalLGLEDMon;};
	121	void SetMinSignalHGLEDMon(int i) {fMinSignalHGLEDMon = i;};
	122	int GetMinSignalHGLEDMon() const {return fMinSignalHGLEDMon;};
	123	void SetMaxSignalHGLEDMon(int i) {fMaxSignalHGLEDMon = i;};
	124	int GetMaxSignalHGLEDMon() const {return fMaxSignalHGLEDMon;};
	125
	126	virtual void EndOfDetectorCycle(AliQAv1::TASKINDEX_t, TObjArray ** list) ;
	127	void GetCalibRefFromOCDB() ;
	128	void GetTruChannelPosition( Int_t &globRow, Int_t &globColumn, Int_t module, Int_t ddl, Int_t branch, Int_t column ) const;
	129	virtual void InitESDs() ;
	130	virtual void InitDigits() ;
	131	virtual void InitRecPoints() ;
	132	virtual void InitRaws() ;
	133	virtual void MakeESDs(AliESDEvent * esd) ;
	134	virtual void MakeDigits() ;
	135	virtual void MakeDigits(TTree * digTree) ;
	136	virtual void MakeRecPoints(TTree * recpoTree) ;
	137	virtual void MakeRaws(AliRawReader* rawReader) ;
	138	virtual void MakeRawsSTU(AliRawReader* rawReader);
	139	virtual void StartOfDetectorCycle() ;
	140
	141	private:
	142	void ConvertProfile2H(TProfile * p, TH2 * histo) ; //change the profile plot to a 2D histogram
	143
	144	Int_t fFittingAlgorithm; // select the fitting algorithm
	145
	146	AliCaloRawAnalyzer *fRawAnalyzer; // for signal fitting
	147	AliCaloRawAnalyzerKStandard *fRawAnalyzerTRU; // for signal fitting, for TRU
	148	AliEMCALGeometry *fGeom; //EMCAL geometry, needed for STU decoding
	149
	150	int fSuperModules; //The number of SuperModules activated
	151	int fFirstPedestalSample; // first sample for pedestal calculation, in bunch
	152	int fLastPedestalSample; // last sample for pedestal calculation, in bunch
	153	int fFirstPedestalSampleTRU; // first sample for pedestal calculation, in bunch
	154	int fLastPedestalSampleTRU; // last sample for pedestal calculation, in bunch
	155	int fMinSignalLG; // minimum signal, for Low Gain channels
	156	int fMaxSignalLG; // maximum signal, for Low Gain channels
	157	int fMinSignalHG; // minimum signal, for High Gain channels
	158	int fMaxSignalHG; // maximum signal, for High Gain channels
	159	int fMinSignalTRU; // minimum signal, for TRU channels
	160	int fMaxSignalTRU; // maximum signal, for TRU channels
	161	int fMinSignalLGLEDMon; // minimum signal, for LEDMon channels, low gain
	162	int fMaxSignalLGLEDMon; // maximum signal, for LEDMon channels, low gain
	163	int fMinSignalHGLEDMon; // minimum signal, for LEDMon channels, high gain
	164	int fMaxSignalHGLEDMon; // maximum signal, for LEDMon channels, high gain
	165
	166	TProfile * fCalibRefHistoPro ; // Profile reference histogram from LED run
	167	TH2F * fCalibRefHistoH2F ; // H2F reference histogram from LED run
	168	TProfile * fLEDMonRefHistoPro ; // Profile reference histogram from LED monitor
	169	TH2F * fHighEmcHistoH2F ; // H2F reference histogram from LED run
	170	// TText ** fTextSM ; //! Text info for each SM
	171	// TLine * fLineCol ; //! line to distinguish the different SM side: A side and C side
	172	// TLine * fLineRow ; //! line to distinguish the different SM sector 0 and 1
	173
	174	ClassDef(AliEMCALQADataMakerRec,5) // description
	175
	176	};
	177
	178	#endif // AliEMCALQADATAMAKERREC_H