ITS/AliITSRecoParam.cxx

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