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 fInwardFlagSA(kFALSE),
 129 fOuterStartLayerSA(2),
 130 fInnerStartLayerSA(3),
 131 fMinNPointsSA(3),
 132 fFactorSAWindowSizes(1.),
 133 fNLoopsSA(33),
 134 fMinPhiSA(0.002),
 135 fMaxPhiSA(0.0145),
 136 fMinLambdaSA(0.003),
 137 fMaxLambdaSA(0.008),
 138 fMinClusterChargeSA(0.),
 139 fSAOnePointTracks(kFALSE),
 140 fSAUseAllClusters(kFALSE),
 141 fFindV0s(kTRUE),
 142 fStoreLikeSignV0s(kFALSE),
 143 fUseUnfoldingInClusterFinderSPD(kFALSE),
 144 fUseUnfoldingInClusterFinderSDD(kTRUE),
 145 fUseUnfoldingInClusterFinderSSD(kTRUE),
 146 fUseSDDClusterSizeSelection(kFALSE),
 147 fMinClusterChargeSDD(0.),
 148 fUseChargeMatchingInClusterFinderSSD(kTRUE),
 149 fTrackleterOnlyOneTrackletPerC2(kTRUE),
 150 fTrackleterPhiWindow(0.08),
 151 fTrackleterZetaWindow(1.00),
 152 fTrackleterRemoveClustersFromOverlaps(kFALSE),
 153 fTrackleterPhiOverlapCut(0.005),
 154 fTrackleterZetaOverlapCut(0.05),
 155 fUseCosmicRunShiftsSSD(kFALSE)
 156 {
 157   //
 158   // constructor
 159   //
 160   SetName("ITS");
 161   SetTitle("ITS");
 162
 163   SetLayersParameters();
 164   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
 165   SetUseTGeoInTracker(3);
 166   SetStepSizeTGeo(0.1);
 167   SetAllowSharedClusters(kTRUE);
 168   SetFindV0s(kTRUE);
 169   SetAddVirtualClustersInDeadZone(kFALSE);
 170   SetUseAmplitudeInfo(kTRUE);
 171   SetClusterErrorsParam(1);
 172   SetClusterMisalError(0.);
 173   SetVertexer3DDefaults();
 174 }
 175 //_____________________________________________________________________________
 176 AliITSRecoParam::~AliITSRecoParam()
 177 {
 178   //
 179   // destructor
 180   //
 181 }
 182 //_____________________________________________________________________________
 183 AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
 184 {
 185   //
 186   // make default reconstruction  parameters for hig  flux env.
 187   //
 188   AliITSRecoParam *param = new AliITSRecoParam();
 189   param->SetVertexerZ();
 190
 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
 199
 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.;
 231
 232   param->fMaxNormChi2NonCForHypothesis = 7.;
 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.;
 266
 267   param->fMaxDZforPrimTrk = 0.4;
 268   param->fMaxDZToUseConstraint = 3.;
 269
 270   param->fMaxDforV0dghtrForProlongation = 30.;
 271   param->fMaxDForProlongation = 40.;//10.;
 272   param->fMaxDZForProlongation = 60.;//20.;
 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;
 279
 280   return param;
 281 }
 282 //_____________________________________________________________________________
 283 AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
 284 {
 285   //
 286   // make default reconstruction  parameters for low  flux env.
 287   //
 288   AliITSRecoParam *param = new AliITSRecoParam();
 289   param->SetVertexer3DDefaults();
 290
 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
 299
 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;
 379   param->SetNLoopsSA(20);
 380   param->fMaxPhiSA = 0.07;
 381   param->fMaxLambdaSA = 0.04;
 382   return param;
 383 }
 384 //_____________________________________________________________________________
 385 AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
 386 {
 387   //
 388   // make default reconstruction  parameters for cosmics
 389   //
 390   AliITSRecoParam *param = new AliITSRecoParam();
 391
 392   // vertexer for cosmics
 393   param->SetVertexer(2);
 394
 395   // find independently ITS SA tracks
 396   param->SetSAUseAllClusters();
 397   param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
 398
 399   // larger seach windows for SA (in case of large misalignments)
 400   param->SetFactorSAWindowSizes(3.);
 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;
 484 }
 485 //_____________________________________________________________________________
 486 AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
 487 {
 488   //
 489   // make special reconstruction parameters for Plane Efficiency study on layer i
 490   //
 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   }
 496   AliITSRecoParam *param;
 497   param = GetHighFluxParam();
 498   param->SetComputePlaneEff();
 499   param->SetLayerToSkip(i);
 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   //
 524   return param;
 525 }
 526 //_____________________________________________________________________________
 527 void 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 }
 551 //_____________________________________________________________________________
 552 void AliITSRecoParam::PrintParameters() const
 553 {
 554   //
 555   // print parameters
 556   //
 557
 558   printf("=============================  AliITSRecoParam::PrintParameters ");
 559   printf("============================= \n\n");
 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++)
 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]);
 576
 577
 578   Dump();
 579
 580   return;
 581 }