Possibility of selecting larger windows in SA tracker, in order to use the tracker...
[u/mrichter/AliRoot.git] / ITS / AliITSRecoParam.cxx
CommitLineData
44347160 1/**************************************************************************
572f41f9 2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
44347160 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
572f41f9 16#include "AliITSRecoParam.h"
17
18/* $Id$ */
44347160 19
20///////////////////////////////////////////////////////////////////////////////
21// //
22// Class with ITS reconstruction parameters //
23// Origin: andrea.dainese@lnl.infn.it //
24// //
25///////////////////////////////////////////////////////////////////////////////
26
27
44347160 28
29ClassImp(AliITSRecoParam)
30
e50912db 31const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
32const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
33const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
34const Double_t AliITSRecoParam::fgkriw=80.0;
35const Double_t AliITSRecoParam::fgkdiw=0.0053;
36const Double_t AliITSRecoParam::fgkX0iw=30.0;
37const Double_t AliITSRecoParam::fgkrcd=61.0;
38const Double_t AliITSRecoParam::fgkdcd=0.0053;
39const Double_t AliITSRecoParam::fgkX0cd=30.0;
40const Double_t AliITSRecoParam::fgkyr=12.8;
41const Double_t AliITSRecoParam::fgkdr=0.03;
42const Double_t AliITSRecoParam::fgkzm=0.2;
43const Double_t AliITSRecoParam::fgkdm=0.40;
44const Double_t AliITSRecoParam::fgkrs=50.0;
45const Double_t AliITSRecoParam::fgkds=0.001;
46const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
47const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
48const Double_t AliITSRecoParam::fgkrPipe=3.;
49const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
50const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
51const Double_t AliITSRecoParam::fgkdPipe=0.0028;
52const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
53const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
54const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
55const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
56const Double_t AliITSRecoParam::fgkX0Air=21.82;
57const Double_t AliITSRecoParam::fgkX0Be=65.19;
58const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
59const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
60const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
61const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
44347160 62
63//_____________________________________________________________________________
572f41f9 64AliITSRecoParam::AliITSRecoParam() : TObject(),
e50912db 65fMaxSnp(1.),
572f41f9 66fNSigmaYLayerForRoadY(0),
67fNSigmaRoadY(0),
68fNSigmaZLayerForRoadZ(0),
69fNSigmaRoadZ(0),
70fNSigma2RoadZC(0),
71fNSigma2RoadYC(0),
72fNSigma2RoadZNonC(0),
73fNSigma2RoadYNonC(0),
74fMaxNormChi2NonCForHypothesis(0),
75fMaxChi2(0),
76fMaxRoad(0),
77fMaxChi2In(0),
78fChi2PerCluster(0),
79fXV(0),
80fYV(0),
81fZV(0),
82fSigmaXV(0),
83fSigmaYV(0),
84fSigmaZV(0),
85fVertexCut(0),
86fMaxDZforPrimTrk(0),
87fMaxDZToUseConstraint(0),
88fMaxDforV0dghtrForProlongation(0),
89fMaxDForProlongation(0),
90fMaxDZForProlongation(0),
91fMinPtForProlongation(0),
92fAddVirtualClustersInDeadZone(kTRUE),
93fZWindowDeadZone(0),
94fSigmaXDeadZoneHit2(0),
95fSigmaZDeadZoneHit2(0),
96fXPassDeadZoneHits(0),
e50912db 97fUseTGeoInTracker(0),
572f41f9 98fAllowSharedClusters(kTRUE),
e50912db 99fClusterErrorsParam(1),
e340bb86 100fFindV0s(kTRUE),
101fFactorSAWindowSizes(1.)
44347160 102{
103 //
104 // constructor
105 //
106 SetLayersParameters();
e50912db 107 SetUseTGeoInTracker(0);
afd25725 108 SetAllowSharedClusters(kTRUE);
109 SetFindV0s(kTRUE);
110 SetAddVirtualClustersInDeadZone(kTRUE);
572f41f9 111 SetUseAmplitudeInfo(kTRUE);
e50912db 112 SetClusterErrorsParam(1);
44347160 113}
114//_____________________________________________________________________________
115AliITSRecoParam::~AliITSRecoParam()
116{
117 //
118 // destructor
119 //
120}
121//_____________________________________________________________________________
122AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
123{
124 //
125 // make default reconstruction parameters for hig flux env.
126 //
127 AliITSRecoParam *param = new AliITSRecoParam();
128
129 param->fMaxSnp = 0.95;
130
131 param->fNSigmaYLayerForRoadY = 4.;
132 param->fNSigmaRoadY = 7.5;
133 param->fNSigmaZLayerForRoadZ = 4.;
134 param->fNSigmaRoadZ = 7.5;
135
136 param->fNSigma2RoadZC = 60.; //7.75^2
137 param->fNSigma2RoadYC = 60.; //7.75^2
138 param->fNSigma2RoadZNonC = 50.; //7.07^2
139 param->fNSigma2RoadYNonC = 50.; //7.07^2
140
141 param->fMaxChi2PerCluster[0] = 11.; //7
142 param->fMaxChi2PerCluster[1] = 12.; //5
143 param->fMaxChi2PerCluster[2] = 12.; //8
144 param->fMaxChi2PerCluster[3] = 5.; //8
145 param->fMaxChi2PerCluster[4] = 12.; //6.5
146
147 param->fMaxNormChi2NonC[0] = 7.;
148 param->fMaxNormChi2NonC[1] = 8.;
149 param->fMaxNormChi2NonC[2] = 8.;
150 param->fMaxNormChi2NonC[3] = 11.;
151 param->fMaxNormChi2NonC[4] = 14.;
152 param->fMaxNormChi2NonC[5] = 25.;
153
154 param->fMaxNormChi2C[0] = 11.;
155 param->fMaxNormChi2C[1] = 13.;
156 param->fMaxNormChi2C[2] = 15.;
157 param->fMaxNormChi2C[3] = 18.;
158 param->fMaxNormChi2C[4] = 30.;
159 param->fMaxNormChi2C[5] = 35.;
afd25725 160
161 param->fMaxNormChi2NonCForHypothesis = 7.;
44347160 162
163 param->fMaxChi2 = 35.;
164
165 param->fMaxChi2s[0] = 25.; //40
166 param->fMaxChi2s[1] = 25.; //40
167 param->fMaxChi2s[2] = 25.; //40
168 param->fMaxChi2s[3] = 25.; //40
169 param->fMaxChi2s[4] = 40.; //40
170 param->fMaxChi2s[5] = 50.; //40
171
172 param->fMaxRoad = 6.;
173
174 // not used
175 param->fMaxChi2In = 16.;
176
177 param->fMaxChi2sR[0] = 10.;
178 param->fMaxChi2sR[1] = 10.;
179 param->fMaxChi2sR[2] = 10.;
180 param->fMaxChi2sR[3] = 10.;
181 param->fMaxChi2sR[4] = 30.;
182 param->fMaxChi2sR[5] = 40.;
183
184 param->fChi2PerCluster = 9.;
185 // not used
186
187 param->fXV = 0.;
188 param->fYV = 0.;
189 param->fZV = 0.;
190 param->fSigmaXV = 0.0050;
191 param->fSigmaYV = 0.0050;
192 param->fSigmaZV = 0.0100;
193
194 param->fVertexCut = 25.;
afd25725 195
196 param->fMaxDZforPrimTrk = 0.4;
197 param->fMaxDZToUseConstraint = 3.;
198
199 param->fMaxDforV0dghtrForProlongation = 30.;
200 param->fMaxDForProlongation = 10.;
201 param->fMaxDZForProlongation = 20.;
202 param->fMinPtForProlongation = 0.120;
203
204 param->fZWindowDeadZone = 2.0;
205 param->fSigmaXDeadZoneHit2 = 0.004/12.;
206 param->fSigmaZDeadZoneHit2 = 0.001/12.;
207 param->fXPassDeadZoneHits = 0.018;
208
44347160 209
210 return param;
211}
212//_____________________________________________________________________________
213AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
214{
215 //
216 // make default reconstruction parameters for low flux env.
217 //
218 return GetHighFluxParam();
219}
220//_____________________________________________________________________________
221AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
222{
223 //
224 // make default reconstruction parameters for cosmics
225 //
e340bb86 226 AliITSRecoParam *param = new AliITSRecoParam();
227
228 param->SetFactorSAWindowSizes(3.);
229
230 param->fMaxSnp = 0.95;
231
232 param->fNSigmaYLayerForRoadY = 4.;
233 param->fNSigmaRoadY = 7.5;
234 param->fNSigmaZLayerForRoadZ = 4.;
235 param->fNSigmaRoadZ = 7.5;
236
237 param->fNSigma2RoadZC = 60.; //7.75^2
238 param->fNSigma2RoadYC = 60.; //7.75^2
239 param->fNSigma2RoadZNonC = 50.; //7.07^2
240 param->fNSigma2RoadYNonC = 50.; //7.07^2
241
242 param->fMaxChi2PerCluster[0] = 11.; //7
243 param->fMaxChi2PerCluster[1] = 12.; //5
244 param->fMaxChi2PerCluster[2] = 12.; //8
245 param->fMaxChi2PerCluster[3] = 5.; //8
246 param->fMaxChi2PerCluster[4] = 12.; //6.5
247
248 param->fMaxNormChi2NonC[0] = 7.;
249 param->fMaxNormChi2NonC[1] = 8.;
250 param->fMaxNormChi2NonC[2] = 8.;
251 param->fMaxNormChi2NonC[3] = 11.;
252 param->fMaxNormChi2NonC[4] = 14.;
253 param->fMaxNormChi2NonC[5] = 25.;
254
255 param->fMaxNormChi2C[0] = 11.;
256 param->fMaxNormChi2C[1] = 13.;
257 param->fMaxNormChi2C[2] = 15.;
258 param->fMaxNormChi2C[3] = 18.;
259 param->fMaxNormChi2C[4] = 30.;
260 param->fMaxNormChi2C[5] = 35.;
261
262 param->fMaxNormChi2NonCForHypothesis = 7.;
263
264 param->fMaxChi2 = 35.;
265
266 param->fMaxChi2s[0] = 25.; //40
267 param->fMaxChi2s[1] = 25.; //40
268 param->fMaxChi2s[2] = 25.; //40
269 param->fMaxChi2s[3] = 25.; //40
270 param->fMaxChi2s[4] = 40.; //40
271 param->fMaxChi2s[5] = 50.; //40
272
273 param->fMaxRoad = 6.;
274
275 // not used
276 param->fMaxChi2In = 16.;
277
278 param->fMaxChi2sR[0] = 10.;
279 param->fMaxChi2sR[1] = 10.;
280 param->fMaxChi2sR[2] = 10.;
281 param->fMaxChi2sR[3] = 10.;
282 param->fMaxChi2sR[4] = 30.;
283 param->fMaxChi2sR[5] = 40.;
284
285 param->fChi2PerCluster = 9.;
286 // not used
287
288 param->fXV = 0.;
289 param->fYV = 0.;
290 param->fZV = 0.;
291 param->fSigmaXV = 0.0050;
292 param->fSigmaYV = 0.0050;
293 param->fSigmaZV = 0.0100;
294
295 param->fVertexCut = 25.;
296
297 param->fMaxDZforPrimTrk = 0.4;
298 param->fMaxDZToUseConstraint = 3.;
299
300 param->fMaxDforV0dghtrForProlongation = 30.;
301 param->fMaxDForProlongation = 10.;
302 param->fMaxDZForProlongation = 20.;
303 param->fMinPtForProlongation = 0.120;
304
305 param->fZWindowDeadZone = 2.0;
306 param->fSigmaXDeadZoneHit2 = 0.004/12.;
307 param->fSigmaZDeadZoneHit2 = 0.001/12.;
308 param->fXPassDeadZoneHits = 0.018;
309
310
311 return param;
44347160 312}
313//_____________________________________________________________________________
314void AliITSRecoParam::SetLayersParameters()
315{
316 //
317 // number of layers and layers spatial resolutions
318 //
319
320 // spatial resolutions of the detectors
321 // y: 12 12 38 38 20 20 micron
322 fSigmaY2[0]=1.44e-6;
323 fSigmaY2[1]=1.44e-6;
324 fSigmaY2[2]=1.444e-5;
325 fSigmaY2[3]=1.444e-5;
326 fSigmaY2[4]=4.0e-6;
327 fSigmaY2[5]=4.0e-6;
328 // z: 120 120 28 28 830 830 micron
329 fSigmaZ2[0]=1.44e-4;
330 fSigmaZ2[1]=1.44e-4;
331 fSigmaZ2[2]=7.84e-6;
332 fSigmaZ2[3]=7.84e-6;
333 fSigmaZ2[4]=6.889e-3;
334 fSigmaZ2[5]=6.889e-3;
335
336 return;
337}
338