]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - TPC/AliTPCRecoParam.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
patch for not checking the number of connections when running proof lite (M.Vala)
[u/mrichter/AliRoot.git] / TPC / AliTPCRecoParam.cxx
... / ...
CommitLineData
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
53ClassImp(AliTPCRecoParam)
54
55
56Bool_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//_____________________________________________________________________________
62AliTPCRecoParam::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{
110 //
111 // constructor
112 //
113 SetName("TPC");
114 SetTitle("TPC");
115 for (Int_t i=0;i<5;i++) fSystematicErrors[i]=0;
116 fCutSharedClusters[0]=0.5; // maximal allowed fraction of shared clusters - shorter track
117 fCutSharedClusters[1]=0.25; // maximal allowed fraction of shared clusters - longer track
118 fClusterMaxRange[0]=1; // y - pad range
119 fClusterMaxRange[1]=1; // z - time bin range
120 fKinkAngleCutChi2[0]=9; // angular cut for kink finder - to create a kink
121 // ~ about 5 % rate for high pt kink finder
122 fKinkAngleCutChi2[1]=12; // angular cut for kink finder - to use the partial track // form kink
123 // ~ about 2 % rate for high pt kink finder
124}
125
126//_____________________________________________________________________________
127AliTPCRecoParam::~AliTPCRecoParam()
128{
129 //
130 // destructor
131 //
132}
133
134void AliTPCRecoParam::Print(const Option_t* /*option*/) const{
135 //
136 //
137 //
138 AliTPCRecoParam::Dump();
139 printf("Systematic errors:\n");
140 const char * cherrs[5]={"sy=","sz=","ssnp=","stheta=","s1pt="};
141 for (Int_t i=0; i<5; i++){
142 printf("%s%f\n",cherrs[i],fSystematicErrors[i]);
143 }
144}
145
146
147AliTPCRecoParam *AliTPCRecoParam::GetLowFluxParam(){
148 //
149 // make default reconstruction parameters for low flux env.
150 //
151 AliTPCRecoParam *param = new AliTPCRecoParam;
152 param->fCtgRange = 10;
153 param->fFirstBin = 0;
154 param->fLastBin = 1000;
155 param->SetName("Low Flux");
156 param->SetTitle("Low Flux");
157 return param;
158}
159
160AliTPCRecoParam *AliTPCRecoParam::GetHighFluxParam(){
161 //
162 // make reco parameters for high flux env.
163 //
164 AliTPCRecoParam *param = new AliTPCRecoParam;
165 param->fCtgRange = 1.05;
166 param->fFirstBin = 0;
167 param->fLastBin = 1000;
168 param->fUseTotCharge=kFALSE;
169 param->SetName("High Flux");
170 param->SetTitle("High Flux");
171 return param;
172}
173
174AliTPCRecoParam *AliTPCRecoParam::GetHLTParam(){
175 //
176 // make reco parameters for high flux env.
177 //
178 AliTPCRecoParam *param = new AliTPCRecoParam;
179 param->fCtgRange = 1.05;
180 param->fFirstBin = 80;
181 param->fLastBin = 1000;
182 param->fMaxSnpTracker = 0.9;
183 param->fMaxC = 0.06;
184 //
185 param->SetName("Hlt Param");
186 param->SetTitle("Hlt Param");
187 param->fBKinkFinder = kFALSE;
188 return param;
189}
190
191AliTPCRecoParam *AliTPCRecoParam::GetLaserTestParam(Bool_t bPedestal){
192 //
193 // special setting for laser
194 //
195 AliTPCRecoParam *param = new AliTPCRecoParam;
196 param->fDumpSignal=kTRUE;
197 param->fCtgRange = 10.05;
198 param->fFirstBin = 0;
199 param->fLastBin = 1000;
200 param->fBCalcPedestal = bPedestal;
201 param->fBDoUnfold = kFALSE;
202 param->fDumpAmplitudeMin = 150;
203 param->fBKinkFinder = kFALSE;
204 param->fMaxSnpTracker = 0.98;
205 param->fMaxC = 0.02;
206 param->fBSpecialSeeding = kTRUE;
207 param->fUseTOFCorrection=kFALSE;
208 //
209 //
210 param->SetName("Laser Flux");
211 param->SetTitle("Laser Flux");
212 return param;
213}
214
215AliTPCRecoParam *AliTPCRecoParam::GetCosmicTestParam(Bool_t bPedestal){
216 //
217 // special setting for cosmic
218 //
219 AliTPCRecoParam *param = new AliTPCRecoParam;
220 param->fDumpSignal=kTRUE;
221 param->fCtgRange = 10.05; // full TPC
222 param->fFirstBin = 60;
223 param->fLastBin = 1000;
224 param->fBCalcPedestal = bPedestal;
225 param->fBDoUnfold = kFALSE;
226 param->fBSpecialSeeding = kTRUE;
227 param->fMaxC = 0.07;
228 param->fBKinkFinder = kFALSE;
229 param->fUseTOFCorrection =kFALSE;
230 param->SetName("Cosmic Flux");
231 param->SetTitle("Cosmic Flux");
232
233 return param;
234}
235
236
237Bool_t AliTPCRecoParam::GetUseTimeCalibration(){
238 //
239 // get
240 //
241 return fgUseTimeCalibration;
242}
243void AliTPCRecoParam::SetUseTimeCalibration(Bool_t useTimeCalibration) {
244 //
245 // set
246 //
247 fgUseTimeCalibration = useTimeCalibration;
248}
249