Store number of TRD tracklets on ESD for tracking and PID. Reuse data
[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//_____________________________________________________________________________
6518a6c5 64AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
f9119eb9 65fTracker(0),
66fITSonly(kFALSE),
67fVertexer(0),
876026b6 68fClusterFinder(0),
69fPID(0),
7203e11a 70fVtxr3DZCutWide(0.),
71fVtxr3DRCutWide(0.),
72fVtxr3DZCutNarrow(0.),
73fVtxr3DRCutNarrow(0.),
74fVtxr3DPhiCutLoose(0.),
75fVtxr3DPhiCutTight(0.),
76fVtxr3DDCACut(0.),
e50912db 77fMaxSnp(1.),
572f41f9 78fNSigmaYLayerForRoadY(0),
79fNSigmaRoadY(0),
80fNSigmaZLayerForRoadZ(0),
81fNSigmaRoadZ(0),
82fNSigma2RoadZC(0),
83fNSigma2RoadYC(0),
84fNSigma2RoadZNonC(0),
85fNSigma2RoadYNonC(0),
1c97ce2f 86fRoadMisal(0),
572f41f9 87fMaxNormChi2NonCForHypothesis(0),
88fMaxChi2(0),
89fMaxRoad(0),
90fMaxChi2In(0),
91fChi2PerCluster(0),
92fXV(0),
93fYV(0),
94fZV(0),
95fSigmaXV(0),
96fSigmaYV(0),
97fSigmaZV(0),
98fVertexCut(0),
99fMaxDZforPrimTrk(0),
100fMaxDZToUseConstraint(0),
101fMaxDforV0dghtrForProlongation(0),
102fMaxDForProlongation(0),
103fMaxDZForProlongation(0),
104fMinPtForProlongation(0),
6518a6c5 105fAddVirtualClustersInDeadZone(kFALSE),
572f41f9 106fZWindowDeadZone(0),
107fSigmaXDeadZoneHit2(0),
108fSigmaZDeadZoneHit2(0),
109fXPassDeadZoneHits(0),
6518a6c5 110fUseTGeoInTracker(3),
9be1d1c7 111fStepSizeTGeo(0.1),
572f41f9 112fAllowSharedClusters(kTRUE),
e50912db 113fClusterErrorsParam(1),
4a66240a 114fComputePlaneEff(kFALSE),
5fbd4fd6 115fHistoPlaneEff(kFALSE),
0ed58a47 116fIPlanePlaneEff(0),
275a301c 117fReadPlaneEffFromOCDB(kFALSE),
0ed58a47 118fMinPtPlaneEff(0),
119fMaxMissingClustersPlaneEff(0),
120fRequireClusterInOuterLayerPlaneEff(kFALSE),
121fRequireClusterInInnerLayerPlaneEff(kFALSE),
122fOnlyConstraintPlaneEff(kFALSE),
2755f080 123fExtendedEtaAcceptance(kFALSE),
b406db2e 124fUseBadZonesFromOCDB(kTRUE),
23197852 125fUseSingleBadChannelsFromOCDB(kFALSE),
126fMinFractionOfBadInRoad(0),
ae00569a 127fAllowProlongationWithEmptyRoad(kFALSE),
12b1afb7 128fInwardFlagSA(kFALSE),
129fOuterStartLayerSA(2),
130fInnerStartLayerSA(3),
131fMinNPointsSA(3),
2755f080 132fFactorSAWindowSizes(1.),
c7d6d7b7 133fNLoopsSA(33),
134fMinPhiSA(0.002),
135fMaxPhiSA(0.0145),
136fMinLambdaSA(0.003),
137fMaxLambdaSA(0.008),
3733ccd2 138fMinClusterChargeSA(0.),
5a03f353 139fSAOnePointTracks(kFALSE),
140fSAUseAllClusters(kFALSE),
6518a6c5 141fFindV0s(kTRUE),
b0160ebc 142fStoreLikeSignV0s(kFALSE),
6518a6c5 143fUseUnfoldingInClusterFinderSPD(kFALSE),
144fUseUnfoldingInClusterFinderSDD(kTRUE),
a86176e3 145fUseUnfoldingInClusterFinderSSD(kTRUE),
b779c7dc 146fUseSDDClusterSizeSelection(kFALSE),
147fMinClusterChargeSDD(0.),
7b116aa1 148fUseChargeMatchingInClusterFinderSSD(kTRUE),
149fTrackleterOnlyOneTrackletPerC2(kTRUE),
150fTrackleterPhiWindow(0.08),
151fTrackleterZetaWindow(1.00),
152fTrackleterRemoveClustersFromOverlaps(kFALSE),
153fTrackleterPhiOverlapCut(0.005),
42ed6062 154fTrackleterZetaOverlapCut(0.05),
155fUseCosmicRunShiftsSSD(kFALSE)
44347160 156{
157 //
158 // constructor
159 //
6518a6c5 160 SetName("ITS");
161 SetTitle("ITS");
162
44347160 163 SetLayersParameters();
2755f080 164 for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
6518a6c5 165 SetUseTGeoInTracker(3);
9be1d1c7 166 SetStepSizeTGeo(0.1);
afd25725 167 SetAllowSharedClusters(kTRUE);
168 SetFindV0s(kTRUE);
6518a6c5 169 SetAddVirtualClustersInDeadZone(kFALSE);
572f41f9 170 SetUseAmplitudeInfo(kTRUE);
e50912db 171 SetClusterErrorsParam(1);
f9119eb9 172 SetClusterMisalError(0.);
7203e11a 173 SetVertexer3DDefaults();
44347160 174}
175//_____________________________________________________________________________
176AliITSRecoParam::~AliITSRecoParam()
177{
178 //
179 // destructor
180 //
181}
182//_____________________________________________________________________________
183AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
184{
185 //
186 // make default reconstruction parameters for hig flux env.
187 //
188 AliITSRecoParam *param = new AliITSRecoParam();
14bcd755 189 param->SetVertexerZ();
44347160 190
1aa655ec 191 // use of bads from OCDB
192 param->SetUseBadZonesFromOCDB(kTRUE);
193 param->SetUseSingleBadChannelsFromOCDB(kFALSE);
194 // extended eta acceptance
195 param->SetExtendedEtaAcceptance(kFALSE);
196 // allow to skip layer if no cluster and no bad
197 param->SetAllowProlongationWithEmptyRoad(kFALSE);
198
9be1d1c7 199
44347160 200 param->fMaxSnp = 0.95;
201
202 param->fNSigmaYLayerForRoadY = 4.;
203 param->fNSigmaRoadY = 7.5;
204 param->fNSigmaZLayerForRoadZ = 4.;
205 param->fNSigmaRoadZ = 7.5;
206
207 param->fNSigma2RoadZC = 60.; //7.75^2
208 param->fNSigma2RoadYC = 60.; //7.75^2
209 param->fNSigma2RoadZNonC = 50.; //7.07^2
210 param->fNSigma2RoadYNonC = 50.; //7.07^2
211
212 param->fMaxChi2PerCluster[0] = 11.; //7
213 param->fMaxChi2PerCluster[1] = 12.; //5
214 param->fMaxChi2PerCluster[2] = 12.; //8
215 param->fMaxChi2PerCluster[3] = 5.; //8
216 param->fMaxChi2PerCluster[4] = 12.; //6.5
217
218 param->fMaxNormChi2NonC[0] = 7.;
219 param->fMaxNormChi2NonC[1] = 8.;
220 param->fMaxNormChi2NonC[2] = 8.;
221 param->fMaxNormChi2NonC[3] = 11.;
222 param->fMaxNormChi2NonC[4] = 14.;
223 param->fMaxNormChi2NonC[5] = 25.;
224
225 param->fMaxNormChi2C[0] = 11.;
226 param->fMaxNormChi2C[1] = 13.;
227 param->fMaxNormChi2C[2] = 15.;
228 param->fMaxNormChi2C[3] = 18.;
229 param->fMaxNormChi2C[4] = 30.;
230 param->fMaxNormChi2C[5] = 35.;
afd25725 231
232 param->fMaxNormChi2NonCForHypothesis = 7.;
44347160 233
234 param->fMaxChi2 = 35.;
235
236 param->fMaxChi2s[0] = 25.; //40
237 param->fMaxChi2s[1] = 25.; //40
238 param->fMaxChi2s[2] = 25.; //40
239 param->fMaxChi2s[3] = 25.; //40
240 param->fMaxChi2s[4] = 40.; //40
241 param->fMaxChi2s[5] = 50.; //40
242
243 param->fMaxRoad = 6.;
244
245 // not used
246 param->fMaxChi2In = 16.;
247
248 param->fMaxChi2sR[0] = 10.;
249 param->fMaxChi2sR[1] = 10.;
250 param->fMaxChi2sR[2] = 10.;
251 param->fMaxChi2sR[3] = 10.;
252 param->fMaxChi2sR[4] = 30.;
253 param->fMaxChi2sR[5] = 40.;
254
255 param->fChi2PerCluster = 9.;
256 // not used
257
258 param->fXV = 0.;
259 param->fYV = 0.;
260 param->fZV = 0.;
261 param->fSigmaXV = 0.0050;
262 param->fSigmaYV = 0.0050;
263 param->fSigmaZV = 0.0100;
264
265 param->fVertexCut = 25.;
afd25725 266
267 param->fMaxDZforPrimTrk = 0.4;
268 param->fMaxDZToUseConstraint = 3.;
269
270 param->fMaxDforV0dghtrForProlongation = 30.;
6518a6c5 271 param->fMaxDForProlongation = 40.;//10.;
272 param->fMaxDZForProlongation = 60.;//20.;
afd25725 273 param->fMinPtForProlongation = 0.120;
274
275 param->fZWindowDeadZone = 2.0;
276 param->fSigmaXDeadZoneHit2 = 0.004/12.;
277 param->fSigmaZDeadZoneHit2 = 0.001/12.;
278 param->fXPassDeadZoneHits = 0.018;
44347160 279
280 return param;
281}
282//_____________________________________________________________________________
283AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
284{
285 //
286 // make default reconstruction parameters for low flux env.
287 //
c7d6d7b7 288 AliITSRecoParam *param = new AliITSRecoParam();
7203e11a 289 param->SetVertexer3DDefaults();
c7d6d7b7 290
1aa655ec 291 // full use of bads from OCDB
292 param->SetUseBadZonesFromOCDB(kTRUE);
293 param->SetUseSingleBadChannelsFromOCDB(kTRUE);
294 // extended eta acceptance
295 param->SetExtendedEtaAcceptance(kTRUE);
296 // allow to skip layer if no cluster and no bad
297 param->SetAllowProlongationWithEmptyRoad(kTRUE);
298
9be1d1c7 299
c7d6d7b7 300 param->fMaxSnp = 0.95;
301
302 param->fNSigmaYLayerForRoadY = 4.;
303 param->fNSigmaRoadY = 7.5;
304 param->fNSigmaZLayerForRoadZ = 4.;
305 param->fNSigmaRoadZ = 7.5;
306
307 param->fNSigma2RoadZC = 60.; //7.75^2
308 param->fNSigma2RoadYC = 60.; //7.75^2
309 param->fNSigma2RoadZNonC = 50.; //7.07^2
310 param->fNSigma2RoadYNonC = 50.; //7.07^2
311
312 param->fMaxChi2PerCluster[0] = 11.; //7
313 param->fMaxChi2PerCluster[1] = 12.; //5
314 param->fMaxChi2PerCluster[2] = 12.; //8
315 param->fMaxChi2PerCluster[3] = 5.; //8
316 param->fMaxChi2PerCluster[4] = 12.; //6.5
317
318 param->fMaxNormChi2NonC[0] = 7.;
319 param->fMaxNormChi2NonC[1] = 8.;
320 param->fMaxNormChi2NonC[2] = 8.;
321 param->fMaxNormChi2NonC[3] = 11.;
322 param->fMaxNormChi2NonC[4] = 14.;
323 param->fMaxNormChi2NonC[5] = 25.;
324
325 param->fMaxNormChi2C[0] = 11.;
326 param->fMaxNormChi2C[1] = 13.;
327 param->fMaxNormChi2C[2] = 15.;
328 param->fMaxNormChi2C[3] = 18.;
329 param->fMaxNormChi2C[4] = 30.;
330 param->fMaxNormChi2C[5] = 35.;
331
332 param->fMaxNormChi2NonCForHypothesis = 7.;
333
334 param->fMaxChi2 = 35.;
335
336 param->fMaxChi2s[0] = 25.; //40
337 param->fMaxChi2s[1] = 25.; //40
338 param->fMaxChi2s[2] = 25.; //40
339 param->fMaxChi2s[3] = 25.; //40
340 param->fMaxChi2s[4] = 40.; //40
341 param->fMaxChi2s[5] = 50.; //40
342
343 param->fMaxRoad = 6.;
344
345 // not used
346 param->fMaxChi2In = 16.;
347
348 param->fMaxChi2sR[0] = 10.;
349 param->fMaxChi2sR[1] = 10.;
350 param->fMaxChi2sR[2] = 10.;
351 param->fMaxChi2sR[3] = 10.;
352 param->fMaxChi2sR[4] = 30.;
353 param->fMaxChi2sR[5] = 40.;
354
355 param->fChi2PerCluster = 9.;
356 // not used
357
358 param->fXV = 0.;
359 param->fYV = 0.;
360 param->fZV = 0.;
361 param->fSigmaXV = 0.0050;
362 param->fSigmaYV = 0.0050;
363 param->fSigmaZV = 0.0100;
364
365 param->fVertexCut = 25.;
366
367 param->fMaxDZforPrimTrk = 0.4;
368 param->fMaxDZToUseConstraint = 3.;
369
370 param->fMaxDforV0dghtrForProlongation = 30.;
371 param->fMaxDForProlongation = 40.;//10.;
372 param->fMaxDZForProlongation = 60.;//20.;
373 param->fMinPtForProlongation = 0.120;
374
375 param->fZWindowDeadZone = 2.0;
376 param->fSigmaXDeadZoneHit2 = 0.004/12.;
377 param->fSigmaZDeadZoneHit2 = 0.001/12.;
378 param->fXPassDeadZoneHits = 0.018;
3733ccd2 379 param->SetNLoopsSA(20);
380 param->fMaxPhiSA = 0.07;
381 param->fMaxLambdaSA = 0.04;
c7d6d7b7 382 return param;
44347160 383}
384//_____________________________________________________________________________
385AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
386{
387 //
388 // make default reconstruction parameters for cosmics
389 //
e340bb86 390 AliITSRecoParam *param = new AliITSRecoParam();
391
f9119eb9 392 // vertexer for cosmics
393 param->SetVertexer(2);
394
9f9cae94 395 // find independently ITS SA tracks
396 param->SetSAUseAllClusters();
397 param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
f9119eb9 398
2755f080 399 // larger seach windows for SA (in case of large misalignments)
400 param->SetFactorSAWindowSizes(3.);
e340bb86 401
402 param->fMaxSnp = 0.95;
403
404 param->fNSigmaYLayerForRoadY = 4.;
405 param->fNSigmaRoadY = 7.5;
406 param->fNSigmaZLayerForRoadZ = 4.;
407 param->fNSigmaRoadZ = 7.5;
408
409 param->fNSigma2RoadZC = 60.; //7.75^2
410 param->fNSigma2RoadYC = 60.; //7.75^2
411 param->fNSigma2RoadZNonC = 50.; //7.07^2
412 param->fNSigma2RoadYNonC = 50.; //7.07^2
413
414 param->fMaxChi2PerCluster[0] = 11.; //7
415 param->fMaxChi2PerCluster[1] = 12.; //5
416 param->fMaxChi2PerCluster[2] = 12.; //8
417 param->fMaxChi2PerCluster[3] = 5.; //8
418 param->fMaxChi2PerCluster[4] = 12.; //6.5
419
420 param->fMaxNormChi2NonC[0] = 7.;
421 param->fMaxNormChi2NonC[1] = 8.;
422 param->fMaxNormChi2NonC[2] = 8.;
423 param->fMaxNormChi2NonC[3] = 11.;
424 param->fMaxNormChi2NonC[4] = 14.;
425 param->fMaxNormChi2NonC[5] = 25.;
426
427 param->fMaxNormChi2C[0] = 11.;
428 param->fMaxNormChi2C[1] = 13.;
429 param->fMaxNormChi2C[2] = 15.;
430 param->fMaxNormChi2C[3] = 18.;
431 param->fMaxNormChi2C[4] = 30.;
432 param->fMaxNormChi2C[5] = 35.;
433
434 param->fMaxNormChi2NonCForHypothesis = 7.;
435
436 param->fMaxChi2 = 35.;
437
438 param->fMaxChi2s[0] = 25.; //40
439 param->fMaxChi2s[1] = 25.; //40
440 param->fMaxChi2s[2] = 25.; //40
441 param->fMaxChi2s[3] = 25.; //40
442 param->fMaxChi2s[4] = 40.; //40
443 param->fMaxChi2s[5] = 50.; //40
444
445 param->fMaxRoad = 6.;
446
447 // not used
448 param->fMaxChi2In = 16.;
449
450 param->fMaxChi2sR[0] = 10.;
451 param->fMaxChi2sR[1] = 10.;
452 param->fMaxChi2sR[2] = 10.;
453 param->fMaxChi2sR[3] = 10.;
454 param->fMaxChi2sR[4] = 30.;
455 param->fMaxChi2sR[5] = 40.;
456
457 param->fChi2PerCluster = 9.;
458 // not used
459
460 param->fXV = 0.;
461 param->fYV = 0.;
462 param->fZV = 0.;
463 param->fSigmaXV = 0.0050;
464 param->fSigmaYV = 0.0050;
465 param->fSigmaZV = 0.0100;
466
467 param->fVertexCut = 25.;
468
469 param->fMaxDZforPrimTrk = 0.4;
470 param->fMaxDZToUseConstraint = 3.;
471
472 param->fMaxDforV0dghtrForProlongation = 30.;
473 param->fMaxDForProlongation = 10.;
474 param->fMaxDZForProlongation = 20.;
475 param->fMinPtForProlongation = 0.120;
476
477 param->fZWindowDeadZone = 2.0;
478 param->fSigmaXDeadZoneHit2 = 0.004/12.;
479 param->fSigmaZDeadZoneHit2 = 0.001/12.;
480 param->fXPassDeadZoneHits = 0.018;
481
482
483 return param;
44347160 484}
485//_____________________________________________________________________________
ae00569a 486AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
487{
488 //
489 // make special reconstruction parameters for Plane Efficiency study on layer i
490 //
0ed58a47 491 if (i<0 || i>=AliITSgeomTGeo::kNLayers) {
492 printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane not in the range [0,5]\n");
493 printf("returning null pointer");
494 return NULL;
495 }
ae00569a 496 AliITSRecoParam *param;
497 param = GetHighFluxParam();
498 param->SetComputePlaneEff();
499 param->SetLayerToSkip(i);
0ed58a47 500 param->SetIPlanePlaneEff(i);
501 // optimized setting for SPD0 (i==0)
502 if (i==0 || i==1) {
503 param->fMinPtPlaneEff = 0.200; // high pt particles
504 param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
505 param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
506 //param->fOnlyConstraintPlaneEff = kTRUE;
507 }
508 if (i==2 || i==3) {
509 param->fMinPtPlaneEff = 0.200; // high pt particles
510 param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
511 param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
512 //param->fOnlyConstraintPlaneEff = kTRUE;
513 }
514 if (i==4) {
515 param->fMinPtPlaneEff = 0.200; // high pt particles
516 param->fMaxMissingClustersPlaneEff = 0; // at most 1 layer out of 5 without cluster
517 param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
518 //param->fOnlyConstraintPlaneEff = kTRUE;
519 }
520 if (i==5) {
521 param->fMinPtPlaneEff = 0.200; // high pt particles
522 }
523 //
ae00569a 524 return param;
525}
526//_____________________________________________________________________________
44347160 527void AliITSRecoParam::SetLayersParameters()
528{
529 //
530 // number of layers and layers spatial resolutions
531 //
532
533 // spatial resolutions of the detectors
534 // y: 12 12 38 38 20 20 micron
535 fSigmaY2[0]=1.44e-6;
536 fSigmaY2[1]=1.44e-6;
537 fSigmaY2[2]=1.444e-5;
538 fSigmaY2[3]=1.444e-5;
539 fSigmaY2[4]=4.0e-6;
540 fSigmaY2[5]=4.0e-6;
541 // z: 120 120 28 28 830 830 micron
542 fSigmaZ2[0]=1.44e-4;
543 fSigmaZ2[1]=1.44e-4;
544 fSigmaZ2[2]=7.84e-6;
545 fSigmaZ2[3]=7.84e-6;
546 fSigmaZ2[4]=6.889e-3;
547 fSigmaZ2[5]=6.889e-3;
548
549 return;
550}
ed446fa3 551//_____________________________________________________________________________
552void AliITSRecoParam::PrintParameters() const
553{
554 //
555 // print parameters
556 //
557
558 printf("============================= AliITSRecoParam::PrintParameters ");
f9119eb9 559 printf("============================= \n\n");
ed446fa3 560 for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
561 if(!fLayersToSkip[i]) {
562 printf("ITS Traking: using layer %d\n",i);
563 } else {
564 printf("ITS Traking: skipping layer %d\n",i);
565 }
566 }
567 for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
568 if(fUseAmplitudeInfo[i]) {
569 printf("ITS Traking: use amplitude info for layer %d\n",i);
570 } else {
571 printf("ITS Traking: don't use amplitude info for layer %d\n",i);
572 }
573 }
574 for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
f9119eb9 575 printf("Layer %d:\n sigmaY2 %f, sigmaZ2 %f\n sigmaMisalY %f, sigmaMisalZ %f\n max norm chi2 for non constrained tracks %f\n max norm chi2 for constrained tracks %f\n max predicted chi2 (cluster & track prol.) %f\n",i,fSigmaY2[i],fSigmaZ2[i],fClusterMisalErrorY[i],fClusterMisalErrorZ[i],fMaxNormChi2NonC[i],fMaxNormChi2C[i],fMaxChi2s[i]);
ed446fa3 576
577
578 Dump();
579
580 return;
581}