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 */
7 //_________________________________________________________________________
8 // Class for performing calibration in PHOS
10 //*-- Author: D.Peressounko (RRC KI & SUBATECH)
13 // --- ROOT system ---
15 #include "TObjArray.h"
18 // --- Standard library ---
20 // --- AliRoot header files ---
21 class AliPHOSConTableDB ;
24 class AliPHOSCalibrator: public TTask{
27 AliPHOSCalibrator() ; // ctor
28 AliPHOSCalibrator(const char* run, const char * title = "Default") ;
29 AliPHOSCalibrator(const AliPHOSCalibrator & ctor) : TTask(ctor) {
30 // cpy ctor: no implementation yet
31 // requested by the Coding Convention
32 Fatal("cpy ctor", "not implemented") ;
35 virtual ~AliPHOSCalibrator() ; // dtor
37 void AddRun(const char * filename) ; //Add one more file to handle
39 virtual void Exec(Option_t * option) ; //Steering method
41 void ScanPedestals(Option_t * option = "append") ;
42 void CalculatePedestals(void) ; //calulates pedestals
43 void ScanGains(Option_t * opt = "append") ; //calculates gains
44 void CalculateGains(void) ; //calculates gains
46 void PlotPedestal(Int_t channel) ; //plots distribution of pedestals for given channel
47 void PlotPedestals(void) ;
48 void PlotGain(Int_t channel) ; //Plot histo with gains for a channel
49 void PlotGains() ; //Plot all gains
51 virtual void Print() const ;
53 TH1F * PedestalHisto(Int_t channel)
54 {return dynamic_cast<TH1F* >(fPedHistos->At(channel)) ;} ;
55 TH1F * GainHisto(Int_t channel)
56 {return dynamic_cast<TH1F* >(fGainHistos->At(channel)) ;} ;
59 TH1F * Pedestals(void){return fhPedestals ;}
60 TH1F * Gains(void){return fhGains ;}
62 //Clean resulting histograms
63 void Reset(void){if(fhPedestals)fhPedestals->Reset("ICE");
64 if(fhGains)fhGains->Reset("ICE");}
66 //Set energy of beam used to calibrate
67 void SetBeamEnergy(Float_t e){fBeamEnergy = e;}
69 void SetConTableDB(const char * filename, const char * title = "Default") ;
70 //Connection table to convert RawId to AbsId
72 void SetNChan(UShort_t nch = 100)
73 {fNChan = nch ; } //Sets number of channels in pedestal histos
75 void SetNGainBins(Int_t nbin = 100)
76 {fNGainBins = nbin ;} //Set number of bins in gain histograms
78 void SetGainMax(Float_t hmax = 0.01)
79 {fGainMax = hmax ;} //Set range of gain histograms
81 void ReadFromASCII(const char * filename) ; //Read gains and pedestals from ascii file
83 void WritePedestals(const char * version="v1") ;
85 void ReadPedestals(const char * version="v1") ;
87 void WriteGains(const char * version="v1") ;
89 void ReadGains(const char * version="v1") ;
91 AliPHOSCalibrator & operator = (const AliPHOSCalibrator & /*rvalue*/){
92 Fatal("operator =","assigment operator is not implemented") ;
102 TList * fRunList ; //list of runs to be handled
103 TObjArray * fPedHistos ; //Resulting histograms of pedestals
104 TObjArray * fGainHistos; //Results of Calibration
106 TH1F * fhPedestals ; //Mean values of pedestals for different channels
107 TH1F * fhPedestalsWid ; //Widths of pedestal distributions for different channels
108 TH1F * fhGains ; //Final Gains from fitting procedure
109 TH1F * fhGainsWid ; //Width of final gains from fit
111 AliPHOSConTableDB * fctdb ; //!Connection map
112 TString fConTableDB ; //Name of ConTableDB
113 TString fConTableDBFile ; //File where ConTableDB is stored
115 Float_t fBeamEnergy ; //Calibration beam energy
116 Float_t fGainAcceptCorr; //Maximal deviation from mean Gain (factor)
117 Float_t fAcceptCorr ; //Maximal deviation of Pedestal from mean for good channel
119 Float_t fGainMax ; //Range used in Gain histos
120 Int_t fNGainBins ; //Number of bins in Gain histos
122 Int_t fNch ; //Number of channels to calibrate
123 UShort_t fNChan ; //Number of bins in pedestal histos
125 ClassDef(AliPHOSCalibrator,1) // description
129 #endif // AliPHOSCALIBRATOR_H