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