Corrected UInt_t <-> Int_t conversion
[u/mrichter/AliRoot.git] / TPC / AliTPCRecoParam.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16
17 ///////////////////////////////////////////////////////////////////////////////
18 //                                                                           //
19 // Class with TPC reconstruction parameters                                  //
20 //                                                                           //  
21 //
22 /*
23   The reconstruction parameters are used in the AliTPCclustererMI and AliTPCtrackerMI
24   
25   They are retrieved:
26   0. User speciefied it in reconstruction macro
27   1. if (not 0) from OCDB  - AliTPCcalibDB::GetRecoParam(eventtype)
28   2. if (not 0 or 1) default parameter - High flux enevironment used  
29
30   FIXME:
31   In the future  reconstruction parameters should be changed on event basis
32   But for the moment, event types are still not defined 
33
34
35   // Setting for systematic errors addition
36   [0] - systematic RMSY
37   [1] - systematic RMSZ
38   [2] - systematic RMSSNP
39   [3] - systematic RMSTheta
40   [4] - systematic RMSCuravture -  systematic error in 1/cm not in 1/pt
41   //
42   //  How to add it example - 3 mm systematic error y, 3 cm systematic error z (drift)
43   Double_t sysError[5]={0.3,3, 0.3/150., 3./150.,0.3/(150*150.)}
44   param->SetSystematicError(sysError);
45
46 */
47                                                                            //
48 ///////////////////////////////////////////////////////////////////////////////
49
50
51 #include "AliTPCRecoParam.h"
52
53 ClassImp(AliTPCRecoParam)
54
55
56 Bool_t AliTPCRecoParam::fgUseTimeCalibration=kTRUE; // flag usage the time dependent calibration
57                                       // to be switched off for pass 0 reconstruction
58                                       // Use static function, other option will be to use 
59                                       // additional specific storage ?
60
61 //_____________________________________________________________________________
62 AliTPCRecoParam::AliTPCRecoParam():
63   AliDetectorRecoParam(),
64   fBClusterSharing(kTRUE),
65   fCtgRange(1.05),       
66   fMaxSnpTracker(0.95),
67   fMaxSnpTrack(0.999),
68   fUseOuterDetectors(kFALSE),
69   fDumpSignal(kFALSE),
70   fFirstBin(0),
71   fLastBin(-1),
72   fBCalcPedestal(kFALSE),
73   fBDoUnfold(kTRUE),
74   fDumpAmplitudeMin(100),
75   fMaxNoise(2.),
76   //
77   fUseOnePadCluster(kTRUE),
78   fMinMaxCutAbs(4.),
79   fMinLeftRightCutAbs(6.),
80   fMinUpDownCutAbs(6.),
81   //
82   fMinMaxCutSigma(4.),
83   fMinLeftRightCutSigma(7.),
84   fMinUpDownCutSigma(8.),
85   fMaxC(0.3),
86   fBSpecialSeeding(kFALSE),
87   fBKinkFinder(kTRUE),
88   fLastSeedRowSec(120),
89   fSeedGapPrim(6),
90   fSeedGapSec(6),
91   fUseFieldCorrection(2),      // use field correction
92   fUseComposedCorrection(kFALSE),      // use field correction
93   fUseRPHICorrection(0),      // use rphi correction
94   fUseRadialCorrection(0),    // use radial correction
95   fUseQuadrantAlignment(0),   // use quadrant alignment
96   fUseSectorAlignment(0),     // use sector alignment
97   fUseDriftCorrectionTime(1), // use drift correction time
98   fUseDriftCorrectionGY(1),   // use drif correction global y
99   fUseGainCorrectionTime(0),  // use gain correction time
100   fUseExBCorrection(1),  // use ExB correction
101   fUseMultiplicityCorrectionDedx(kTRUE), // use Dedx multiplicity correction
102   fUseAlignmentTime(kTRUE),              // use time dependent alignment correction
103   //
104   fUseTotCharge(kTRUE),          // switch use total or max charge
105   fMinFraction(0.01),           // truncated mean - lower threshold
106   fMaxFaction(0.7),            // truncated mean - upper threshold
107   fNeighborRowsDedx(2),           // neighbour rows for below threshold dEdx calculation
108   fUseTOFCorrection(kTRUE),
109   fUseSystematicCorrelation(kTRUE)
110 {
111   //
112   // constructor
113   //
114   SetName("TPC");
115   SetTitle("TPC");
116   for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
117   fCutSharedClusters[0]=0.5; // maximal allowed fraction of shared clusters - shorter track
118   fCutSharedClusters[1]=0.25; // maximal allowed fraction of shared clusters - longer  track
119   fClusterMaxRange[0]=1;     // y - pad      range
120   fClusterMaxRange[1]=1;     // z - time bin range
121   fKinkAngleCutChi2[0]=9;    // angular cut for kink finder - to create a kink
122                              // ~ about 5 % rate  for high pt kink finder
123   fKinkAngleCutChi2[1]=12;    // angular cut for kink finder - to use the partial track                             // form kink 
124                              // ~ about 2 % rate  for high pt kink finder
125 }
126
127 //_____________________________________________________________________________
128 AliTPCRecoParam::~AliTPCRecoParam() 
129 {
130   //
131   // destructor
132   //  
133 }
134
135 void AliTPCRecoParam::Print(const Option_t* /*option*/) const{
136   //
137   //
138   //
139   AliTPCRecoParam::Dump();
140   printf("Systematic errors:\n");
141   const char * cherrs[5]={"sy=","sz=","ssnp=","stheta=","s1pt="};
142   for (Int_t i=0; i<5; i++){
143     printf("%s%f\n",cherrs[i],fSystematicErrors[i]);
144   }
145 }
146
147
148 AliTPCRecoParam *AliTPCRecoParam::GetLowFluxParam(){
149   //
150   // make default reconstruction  parameters for low  flux env.
151   //
152   AliTPCRecoParam *param = new AliTPCRecoParam;
153   param->fCtgRange = 10;
154   param->fFirstBin = 0;
155   param->fLastBin  = 1000;
156   param->SetName("Low Flux");
157   param->SetTitle("Low Flux");
158   return param;
159 }
160
161 AliTPCRecoParam *AliTPCRecoParam::GetHighFluxParam(){
162   //
163   // make reco parameters for high flux env.
164   //
165   AliTPCRecoParam *param = new AliTPCRecoParam;
166   param->fCtgRange = 1.05;
167   param->fFirstBin = 0;
168   param->fLastBin  = 1000;  
169   param->fUseTotCharge=kFALSE;
170   param->SetName("High Flux");
171   param->SetTitle("High Flux");
172   return param;
173 }
174
175 AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
176   //
177   // make reco parameters for high flux env.
178   //
179   AliTPCRecoParam *param = new AliTPCRecoParam;
180   param->fCtgRange = 1.05;
181   param->fFirstBin = 80;
182   param->fLastBin  = 1000;  
183   param->fMaxSnpTracker = 0.9; 
184   param->fMaxC          = 0.06; 
185   //
186   param->SetName("Hlt Param");
187   param->SetTitle("Hlt Param"); 
188   param->fBKinkFinder   = kFALSE;
189   return param;
190 }
191
192 AliTPCRecoParam *AliTPCRecoParam::GetLaserTestParam(Bool_t bPedestal){
193   //
194   // special setting for laser
195   //
196   AliTPCRecoParam *param = new AliTPCRecoParam;
197   param->fDumpSignal=kTRUE;
198   param->fCtgRange = 10.05;
199   param->fFirstBin = 0;
200   param->fLastBin  = 1000;
201   param->fBCalcPedestal = bPedestal;
202   param->fBDoUnfold     = kFALSE;
203   param->fDumpAmplitudeMin = 150;
204   param->fBKinkFinder   = kFALSE;
205   param->fMaxSnpTracker = 0.98;
206   param->fMaxC          = 0.02;
207   param->fBSpecialSeeding = kTRUE;
208   param->fUseTOFCorrection=kFALSE;
209   //
210   //
211   param->SetName("Laser Flux");
212   param->SetTitle("Laser Flux");
213   return param;
214 }
215
216 AliTPCRecoParam *AliTPCRecoParam::GetCosmicTestParam(Bool_t bPedestal){
217   //
218   // special setting for cosmic 
219   // 
220   AliTPCRecoParam *param = new AliTPCRecoParam;
221   param->fDumpSignal=kTRUE;
222   param->fCtgRange = 10.05;    // full TPC
223   param->fFirstBin = 60;
224   param->fLastBin  = 1000;
225   param->fBCalcPedestal = bPedestal;
226   param->fBDoUnfold     = kFALSE;
227   param->fBSpecialSeeding = kTRUE;
228   param->fMaxC          = 0.07;
229   param->fBKinkFinder   = kFALSE;
230   param->fUseTOFCorrection =kFALSE;
231   param->SetName("Cosmic Flux");
232   param->SetTitle("Cosmic Flux");
233
234   return param;
235 }
236
237
238 Bool_t  AliTPCRecoParam::GetUseTimeCalibration(){ 
239   //
240   // get
241   //
242   return fgUseTimeCalibration;
243 }
244 void    AliTPCRecoParam::SetUseTimeCalibration(Bool_t useTimeCalibration) {
245   //
246   // set 
247   //
248   fgUseTimeCalibration = useTimeCalibration;
249 }
250