[u/mrichter/AliRoot.git] / PHOS / AliPHOSCalibrator.h

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


//_________________________________________________________________________
//  Class for performing calibration in PHOS     
//                  
//*-- Author: D.Peressounko (RRC KI & SUBATECH)


// --- ROOT system ---
#include "TTask.h"
#include "TObjArray.h"
class TH1F ; 

// --- Standard library ---

// --- AliRoot header files ---
class AliPHOSConTableDB ;
class AliPHOSDigit ;

class AliPHOSCalibrator: public TTask{

public:
  AliPHOSCalibrator() ;          // ctor
  AliPHOSCalibrator(const char* run, const char * title = "Default",Bool_t toSplit = kTRUE) ;
  virtual ~AliPHOSCalibrator() ; // dtor

  void AddRun(const char * filename) ; //Add one more file to handle

  virtual void Exec(Option_t * option) ; //Steering method 

  void ScanPedestals(Option_t * option = "append") ;
  void CalculatePedestals(void) ; //calulates pedestals
  void ScanGains(Option_t * opt = "append") ; //calculates gains
  void CalculateGains(void) ; //calculates gains

  void PlotPedestal(Int_t channel) ; //plots distribution of pedestals for given channel
  void PlotPedestals(void) ;
  void PlotGain(Int_t channel) ; //Plot histo with gains for a channel
  void PlotGains() ;             //Plot all gains

  void Print(const Option_t * option) const ;

  TH1F * PedestalHisto(Int_t channel)
    {return dynamic_cast<TH1F* >(fPedHistos->At(channel)) ;} ;
  TH1F * GainHisto(Int_t channel)
    {return dynamic_cast<TH1F* >(fGainHistos->At(channel)) ;} ;

  
  TH1F * Pedestals(void){return fhPedestals ;}
  TH1F * Gains(void){return fhGains ;}

  void SetPedestalPattern(UShort_t pattern = 257)
    {fPedPat = pattern ;} ;   //Sets trigger pattern for PEDESTAL events
  void SetPulserPattern(UShort_t pattern = 33)
    {fPulPat = pattern ;} ;   //Sets trigger pattern for PULSER events
  void SetLEDPattern(UShort_t pattern = 129)
    {fLEDPat = pattern ;} ;   //Sets trigger pattern for LED events
  void SetWideBeamPattern(UShort_t pattern = 1027)
    {fWBPat = pattern ;} ;    //Sets trigger pattern for WIDE BEAM events
  void SetNarrowBeamPattern(UShort_t pattern = 1029)
    {fNBPat = pattern ;} ;    //Sets trigger pattern for NARROW BEAM events

  void SetConTableDB(const char * filename, const char * title = "Default") ;
       //Connection table to convert RawId to AbsId

  void SetNChan(UShort_t nch = 100)
    {fNChan = nch ; }         //Sets number of channels in pedestal histos

  void SetNGainBins(Int_t nbin = 100)
    {fNGainBins = nbin ;}  //Set number of bins in gain histograms

  void SetGainMax(Float_t hmax = 0.01)
    {fGainMax = hmax ;}    //Set range of gain histograms

  void WritePedestals(const char * version="v1",
		      Int_t begValidRange = 0,
		      Int_t endValidRange = 0) ;

  void ReadPedestals(const char * version="v1",
		     Int_t ValidRange = 0) ;
		      
  void WriteGains(const char * version="v1",
		      Int_t begValidRange = 0,
		      Int_t endValidRange = 0) ;

  void ReadGains(const char * version="v1",
		     Int_t ValidRange = 0) ;


private:
  void Init() ;
 
private:
  TList  * fRunList ;          //list of runs to be handled
  TObjArray * fPedHistos ;     //Resulting histograms of pedestals
  TObjArray * fGainHistos;     //Results of Calibration 

  TH1F *   fhPedestals ;      //Mean values of pedestals for different channels
  TH1F *   fhPedestalsWid ;   //Widths of pedestal distributions for different channels
  TH1F *   fhGains ;          //Final Gains from fitting procedure
  TH1F *   fhGainsWid ;       //Width of final gains from fit

  AliPHOSConTableDB * fctdb ; //!Connection map
  TString  fConTableDB ;      //Name of ConTableDB
  TString  fConTableDBFile ;  //File where ConTableDB is stored

  Float_t  fGainAcceptCorr;   //Maximal deviation from mean Gain (factor)
  Float_t  fAcceptCorr ;      //Maximal deviation of Pedestal from mean for good channel

  Float_t  fGainMax ;         //Range used in Gain histos
  Int_t    fNGainBins ;       //Number of bins in Gain histos

  Int_t    fNch ;             //Number of channels to calibrate
  UShort_t fNChan ;           //Number of bins in pedestal histos

  UShort_t fPedPat ;     //trigger pattern for PEDESTAL events
  UShort_t fPulPat ;     //trigger pattern for PULSER events
  UShort_t fLEDPat ;     //trigger pattern for LED events
  UShort_t fWBPat ;      //trigger pattern for WIDE BEAM events
  UShort_t fNBPat ;      //trigger pattern for NARROW BEAM events
  Bool_t   fToSplit ;    //To split digits from header

  ClassDef(AliPHOSCalibrator,1)  // description 

};

#endif // AliPHOSCALIBRATOR_H
Commit	Line	Data
a6edf90e	1	#ifndef ALIPHOSCALIBRATOR_H
	2	#define ALIPHOSCALIBRATOR_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6
	7	//_________________________________________________________________________
	8	// Class for performing calibration in PHOS
	9	//
	10	//*-- Author: D.Peressounko (RRC KI & SUBATECH)
	11
	12
	13	// --- ROOT system ---
	14	#include "TTask.h"
	15	#include "TObjArray.h"
	16	class TH1F ;
	17
	18	// --- Standard library ---
	19
	20	// --- AliRoot header files ---
	21	class AliPHOSConTableDB ;
	22	class AliPHOSDigit ;
	23
	24	class AliPHOSCalibrator: public TTask{
	25
	26	public:
	27	AliPHOSCalibrator() ; // ctor
	28	AliPHOSCalibrator(const char* run, const char * title = "Default",Bool_t toSplit = kTRUE) ;
	29	virtual ~AliPHOSCalibrator() ; // dtor
	30
	31	void AddRun(const char * filename) ; //Add one more file to handle
	32
	33	virtual void Exec(Option_t * option) ; //Steering method
	34
	35	void ScanPedestals(Option_t * option = "append") ;
	36	void CalculatePedestals(void) ; //calulates pedestals
	37	void ScanGains(Option_t * opt = "append") ; //calculates gains
	38	void CalculateGains(void) ; //calculates gains
	39
	40	void PlotPedestal(Int_t channel) ; //plots distribution of pedestals for given channel
	41	void PlotPedestals(void) ;
	42	void PlotGain(Int_t channel) ; //Plot histo with gains for a channel
	43	void PlotGains() ; //Plot all gains
	44
	45	void Print(const Option_t * option) const ;
	46
	47	TH1F * PedestalHisto(Int_t channel)
	48	{return dynamic_cast<TH1F* >(fPedHistos->At(channel)) ;} ;
	49	TH1F * GainHisto(Int_t channel)
	50	{return dynamic_cast<TH1F* >(fGainHistos->At(channel)) ;} ;
	51
	52
	53	TH1F * Pedestals(void){return fhPedestals ;}
	54	TH1F * Gains(void){return fhGains ;}
	55
	56	void SetPedestalPattern(UShort_t pattern = 257)
	57	{fPedPat = pattern ;} ; //Sets trigger pattern for PEDESTAL events
	58	void SetPulserPattern(UShort_t pattern = 33)
	59	{fPulPat = pattern ;} ; //Sets trigger pattern for PULSER events
	60	void SetLEDPattern(UShort_t pattern = 129)
	61	{fLEDPat = pattern ;} ; //Sets trigger pattern for LED events
	62	void SetWideBeamPattern(UShort_t pattern = 1027)
	63	{fWBPat = pattern ;} ; //Sets trigger pattern for WIDE BEAM events
	64	void SetNarrowBeamPattern(UShort_t pattern = 1029)
65	{fNBPat = pattern ;} ; //Sets trigger pattern for NARROW BEAM events
66
67	void SetConTableDB(const char * filename, const char * title = "Default") ;
68	//Connection table to convert RawId to AbsId
69
70	void SetNChan(UShort_t nch = 100)
71	{fNChan = nch ; } //Sets number of channels in pedestal histos
72
73	void SetNGainBins(Int_t nbin = 100)
74	{fNGainBins = nbin ;} //Set number of bins in gain histograms
75
76	void SetGainMax(Float_t hmax = 0.01)
77	{fGainMax = hmax ;} //Set range of gain histograms
78
79	void WritePedestals(const char * version="v1",
80	Int_t begValidRange = 0,
81	Int_t endValidRange = 0) ;
82
83	void ReadPedestals(const char * version="v1",
84	Int_t ValidRange = 0) ;
85
86	void WriteGains(const char * version="v1",
87	Int_t begValidRange = 0,
88	Int_t endValidRange = 0) ;
89
90	void ReadGains(const char * version="v1",
91	Int_t ValidRange = 0) ;
92
93
94
95	private:
96	void Init() ;
97
98	private:
99	TList * fRunList ; //list of runs to be handled
100	TObjArray * fPedHistos ; //Resulting histograms of pedestals
101	TObjArray * fGainHistos; //Results of Calibration
102
103	TH1F * fhPedestals ; //Mean values of pedestals for different channels
104	TH1F * fhPedestalsWid ; //Widths of pedestal distributions for different channels
105	TH1F * fhGains ; //Final Gains from fitting procedure
106	TH1F * fhGainsWid ; //Width of final gains from fit
107
108	AliPHOSConTableDB * fctdb ; //!Connection map
109	TString fConTableDB ; //Name of ConTableDB
110	TString fConTableDBFile ; //File where ConTableDB is stored
111
112	Float_t fGainAcceptCorr; //Maximal deviation from mean Gain (factor)
113	Float_t fAcceptCorr ; //Maximal deviation of Pedestal from mean for good channel
114
115	Float_t fGainMax ; //Range used in Gain histos
116	Int_t fNGainBins ; //Number of bins in Gain histos
117
118	Int_t fNch ; //Number of channels to calibrate
119	UShort_t fNChan ; //Number of bins in pedestal histos
120
121	UShort_t fPedPat ; //trigger pattern for PEDESTAL events
122	UShort_t fPulPat ; //trigger pattern for PULSER events
123	UShort_t fLEDPat ; //trigger pattern for LED events
124	UShort_t fWBPat ; //trigger pattern for WIDE BEAM events
125	UShort_t fNBPat ; //trigger pattern for NARROW BEAM events
126	Bool_t fToSplit ; //To split digits from header
127
128	ClassDef(AliPHOSCalibrator,1) // description
129
130	};
131
132	#endif // AliPHOSCALIBRATOR_H