Added method for filling TProfile
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecParam.h
1 #ifndef ALIEMCALRECPARAM_H
2 #define ALIEMCALRECPARAM_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice                               */
5
6 /* $Id$ */
7
8 //-----------------------------------------------------------------------------
9 // Container of EMCAL reconstruction parameters
10 // The purpose of this object is to store it to OCDB
11 // and retrieve it in AliEMCALClusterizerv1, AliEMCALPID,
12 // AliEMCALTracker and use it to configure AliEMCALRawUtils
13 // 
14 // 
15 // Author: Yuri Kharlov
16 //-----------------------------------------------------------------------------
17
18 // --- ROOT system ---
19
20 #include "AliDetectorRecoParam.h" 
21 #include "AliLog.h"
22
23 class AliEMCALRecParam : public AliDetectorRecoParam
24 {
25  public:
26
27   enum AliEMCALClusterizerFlag
28   {
29     kClusterizerv1  = 0,
30     kClusterizerNxN = 1   
31   };
32   
33   AliEMCALRecParam() ;
34   AliEMCALRecParam(const AliEMCALRecParam& recParam);
35   AliEMCALRecParam& operator = (const AliEMCALRecParam& recParam);
36   virtual ~AliEMCALRecParam() {}
37   
38   //Clustering (Unfolding : Cynthia)
39   Float_t GetClusteringThreshold() const     {return fClusteringThreshold ;}
40   Float_t GetW0                 () const     {return fW0                  ;}
41   Float_t GetMinECut            () const     {return fMinECut             ;}
42   Float_t GetLocMaxCut          () const     {return fLocMaxCut           ;}
43   Float_t GetTimeCut            () const     {return fTimeCut             ;}
44   Float_t GetTimeMin            () const     {return fTimeMin             ;}
45   Float_t GetTimeMax            () const     {return fTimeMax             ;}
46   Bool_t  GetUnfold             () const     {return fUnfold              ;}
47         
48   void SetClusteringThreshold(Float_t thrsh)     {fClusteringThreshold = thrsh;}
49   void SetW0                 (Float_t w0)        {fW0        = w0         ;}
50   void SetMinECut            (Float_t ecut)      {fMinECut   = ecut       ;}
51   void SetLocMaxCut          (Float_t locMaxCut) {fLocMaxCut = locMaxCut  ;}
52   void SetTimeCut            (Float_t t)         {fTimeCut   = t          ;}
53   void SetTimeMin            (Float_t t)         {fTimeMin   = t          ;}
54   void SetTimeMax            (Float_t t)         {fTimeMax   = t          ;}
55   void SetUnfold             (Bool_t unfold)     {fUnfold = unfold ;}
56   
57   //PID (Guenole)
58   Double_t GetGamma(Int_t i, Int_t j) const       {return fGamma[i][j];} 
59   Double_t GetGammaEnergyProb(Int_t i) const      {return fGammaEnergyProb[i];} 
60   Double_t GetGamma1to10(Int_t i, Int_t j) const  {return fGamma1to10[i][j];}   // not used
61   Double_t GetHadron(Int_t i, Int_t j) const      {return fHadron[i][j];}
62   Double_t GetHadron1to10(Int_t i, Int_t j) const {return fHadron1to10[i][j];}   // not used
63   Double_t GetHadronEnergyProb(Int_t i) const     {return fHadronEnergyProb[i];}
64   Double_t GetPiZero(Int_t i, Int_t j) const      {return fPiZero[i][j];}
65   Double_t GetPiZeroEnergyProb(Int_t i) const     {return fPiZeroEnergyProb[i];}
66   
67   void SetGamma(Int_t i, Int_t j,Double_t param )       {fGamma[i][j]=param;}
68   void SetGammaEnergyProb(Int_t i, Double_t param )     {fGammaEnergyProb[i]=param;}
69   void SetGamma1to10(Int_t i, Int_t j,Double_t param )  {fGamma1to10[i][j]=param;}
70   void SetHadron(Int_t i, Int_t j,Double_t param )      {fHadron[i][j]=param;}
71   void SetHadron1to10(Int_t i, Int_t j,Double_t param ) {fHadron1to10[i][j]=param;}
72   void SetHadronEnergyProb(Int_t i,Double_t param )     {fHadronEnergyProb[i]=param;}
73   void SetPiZero(Int_t i, Int_t j,Double_t param)       {fPiZero[i][j]=param;}
74   void SetPiZeroEnergyProb(Int_t i,Double_t param)      {fPiZeroEnergyProb[i]=param;}
75   
76   //Track Matching (Alberto)
77   /* track matching cut setters */
78   void SetTrkCutX(Double_t value)        {fTrkCutX = value;}
79   void SetTrkCutY(Double_t value)        {fTrkCutY = value;}
80   void SetTrkCutZ(Double_t value)        {fTrkCutZ = value;}
81   void SetTrkCutR(Double_t value)        {fTrkCutR = value;}
82   void SetTrkCutAlphaMin(Double_t value) {fTrkCutAlphaMin = value;}
83   void SetTrkCutAlphaMax(Double_t value) {fTrkCutAlphaMax = value;}
84   void SetTrkCutAngle(Double_t value)    {fTrkCutAngle = value;}
85   void SetTrkCutNITS(Double_t value)     {fTrkCutNITS = value;}
86   void SetTrkCutNTPC(Double_t value)     {fTrkCutNTPC = value;}
87   /* track matching cut getters */
88   Double_t GetTrkCutX() const        {return fTrkCutX;}
89   Double_t GetTrkCutY() const        {return fTrkCutY;}
90   Double_t GetTrkCutZ() const        {return fTrkCutZ;}
91   Double_t GetTrkCutR() const        {return fTrkCutR;}
92   Double_t GetTrkCutAlphaMin() const {return fTrkCutAlphaMin;}
93   Double_t GetTrkCutAlphaMax() const {return fTrkCutAlphaMax;}
94   Double_t GetTrkCutAngle() const    {return fTrkCutAngle;}
95   Double_t GetTrkCutNITS() const     {return fTrkCutNITS;}
96   Double_t GetTrkCutNTPC() const     {return fTrkCutNTPC;}
97   
98   //Raw signal fitting (Jenn)
99   /* raw signal setters */
100   void SetHighLowGainFactor(Double_t value) {fHighLowGainFactor = value;}
101   void SetOrderParameter(Int_t value)       {fOrderParameter = value;}
102   void SetTau(Double_t value)               {fTau = value;}
103   void SetNoiseThreshold(Int_t value)       {fNoiseThreshold = value;}
104   void SetNPedSamples(Int_t value)          {fNPedSamples = value;} 
105   void SetRemoveBadChannels(Bool_t val)     {fRemoveBadChannels=val; }
106   void SetFittingAlgorithm(Int_t val)       {fFittingAlgorithm=val; }
107   void SetFALTROUsage(Bool_t val)           {fUseFALTRO=val; }
108   void SetLEDFit(Bool_t val)                {fFitLEDEvents=val; }
109
110         
111   /* raw signal getters */
112   Double_t GetHighLowGainFactor() const {return fHighLowGainFactor;}
113   Int_t    GetOrderParameter()    const {return fOrderParameter;}
114   Double_t GetTau()               const {return fTau;}
115   Int_t    GetNoiseThreshold()    const {return fNoiseThreshold;}
116   Int_t    GetNPedSamples()       const {return fNPedSamples;}
117   Bool_t   GetRemoveBadChannels() const {return fRemoveBadChannels;}
118   Int_t    GetFittingAlgorithm()  const {return fFittingAlgorithm; }
119   Bool_t   UseFALTRO()            const {return fUseFALTRO; }
120   Bool_t   FitLEDEvents()         const {return fFitLEDEvents; }
121
122   //Unfolding (Adam)
123   Double_t GetSSPars(Int_t i) const   {return fSSPars[i];}
124   Double_t GetPar5(Int_t i) const     {return fPar5[i];}
125   Double_t GetPar6(Int_t i) const     {return fPar6[i];}
126   void SetSSPars(Int_t i, Double_t param )     {fSSPars[i]=param;}
127   void SetPar5(Int_t i, Double_t param )       {fPar5[i]=param;}
128   void SetPar6(Int_t i, Double_t param )       {fPar6[i]=param;}
129
130
131   virtual void Print(Option_t * option="") const ;
132   
133   static AliEMCALRecParam* GetDefaultParameters();
134   static AliEMCALRecParam* GetLowFluxParam();
135   static AliEMCALRecParam* GetHighFluxParam();
136   static AliEMCALRecParam* GetCalibParam();
137   static AliEMCALRecParam* GetCosmicParam();
138
139   static const  TObjArray* GetMappings();
140   
141   void    SetClusterizerFlag(Short_t val) { fClusterizerFlag = val;  }
142   Short_t GetClusterizerFlag() const      { return fClusterizerFlag; }
143   
144  private:
145   //Clustering
146   Float_t fClusteringThreshold ; // Minimum energy to seed a EC digit in a cluster
147   Float_t fW0 ;                  // Logarithmic weight for the cluster center of gravity calculation
148   Float_t fMinECut;              // Minimum energy for a digit to be a member of a cluster
149   Bool_t  fUnfold;               // Flag to perform cluster unfolding
150   Float_t fLocMaxCut;            // Minimum energy difference to consider local maxima in a cluster
151   Float_t fTimeCut ;             // Maximum time of digits with respect to EMC cluster max.
152   Float_t fTimeMin ;             // Minimum time of digits
153   Float_t fTimeMax ;             // Maximum time of digits
154   Short_t fClusterizerFlag ;     // Choice of the clusterizer; Default selection (v1) is zero
155
156   //PID (Guenole)
157   Double_t fGamma[6][6];         // Parameter to Compute PID for photons     
158   Double_t fGamma1to10[6][6];    // Parameter to Compute PID not used
159   Double_t fHadron[6][6];            // Parameter to Compute PID for hadrons     
160   Double_t fHadron1to10[6][6];   // Parameter to Compute PID for hadrons between 1 and 10 GeV    
161   Double_t fHadronEnergyProb[6]; // Parameter to Compute PID for energy ponderation for hadrons          
162   Double_t fPiZeroEnergyProb[6]; // Parameter to Compute PID for energy ponderation for Pi0      
163   Double_t fGammaEnergyProb[6];  // Parameter to Compute PID for energy ponderation for gamma    
164   Double_t fPiZero[6][6];        // Parameter to Compute PID for pi0     
165   
166   
167   //Track-Matching (Alberto)
168   Double_t  fTrkCutX;              // X-difference cut for track matching
169   Double_t  fTrkCutY;              // Y-difference cut for track matching
170   Double_t  fTrkCutZ;              // Z-difference cut for track matching
171   Double_t  fTrkCutR;              // cut on allowed track-cluster distance
172   Double_t  fTrkCutAlphaMin;       // cut on 'alpha' parameter for track matching (min)
173   Double_t  fTrkCutAlphaMax;       // cut on 'alpha' parameter for track matching (min)
174   Double_t  fTrkCutAngle;          // cut on relative angle between different track points for track matching
175   Double_t  fTrkCutNITS;           // Number of ITS hits for track matching
176   Double_t  fTrkCutNTPC;           // Number of TPC hits for track matching
177   
178   //Raw signal fitting parameters (Jenn)
179   Double_t fHighLowGainFactor;     // gain factor to convert between high and low gain
180   Int_t    fOrderParameter;        // order parameter for raw signal fit
181   Double_t fTau;                   // decay constant for raw signal fit
182   Int_t    fNoiseThreshold;        // threshold to consider signal or noise
183   Int_t    fNPedSamples;           // number of time samples to use in pedestal calculation
184   Bool_t   fRemoveBadChannels;     // select if bad channels are removed before fitting
185   Int_t    fFittingAlgorithm;      // select the fitting algorithm
186   Bool_t   fUseFALTRO;             // get FALTRO (trigger) and put it on trigger digits.
187   Bool_t   fFitLEDEvents;          // fit LED events or not
188         
189   //Shower shape parameters (Adam)
190   Double_t fSSPars[8]; // Unfolding shower shape parameters
191   Double_t fPar5[3];   // UF SSPar nr 5
192   Double_t fPar6[3];   // UF SSPar nr 6
193
194   static TObjArray* fgkMaps;       // ALTRO mappings for RCU0..RCUX
195   
196   ClassDef(AliEMCALRecParam,14)     // Reconstruction parameters
197     
198     } ;
199
200 #endif //  ALIEMCALRECPARAM_H
201