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