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 fMaxSnp(1.),
  66 fNSigmaYLayerForRoadY(0),
  67 fNSigmaRoadY(0),
  68 fNSigmaZLayerForRoadZ(0),
  69 fNSigmaRoadZ(0),
  70 fNSigma2RoadZC(0),
  71 fNSigma2RoadYC(0),
  72 fNSigma2RoadZNonC(0),
  73 fNSigma2RoadYNonC(0),
  74 fMaxNormChi2NonCForHypothesis(0),
  75 fMaxChi2(0),
  76 fMaxRoad(0),
  77 fMaxChi2In(0),
  78 fChi2PerCluster(0),
  79 fXV(0),
  80 fYV(0),
  81 fZV(0),
  82 fSigmaXV(0),
  83 fSigmaYV(0),
  84 fSigmaZV(0),
  85 fVertexCut(0),
  86 fMaxDZforPrimTrk(0),
  87 fMaxDZToUseConstraint(0),
  88 fMaxDforV0dghtrForProlongation(0),
  89 fMaxDForProlongation(0),
  90 fMaxDZForProlongation(0),
  91 fMinPtForProlongation(0),
  92 fAddVirtualClustersInDeadZone(kFALSE),
  93 fZWindowDeadZone(0),
  94 fSigmaXDeadZoneHit2(0),
  95 fSigmaZDeadZoneHit2(0),
  96 fXPassDeadZoneHits(0),
  97 fUseTGeoInTracker(3),
  98 fAllowSharedClusters(kTRUE),
  99 fClusterErrorsParam(1),
 100 fComputePlaneEff(kFALSE),
 101 fHistoPlaneEff(kFALSE),
 102 fExtendedEtaAcceptance(kFALSE),
 103 fUseDeadZonesFromOCDB(kFALSE),
 104 fAllowProlongationWithEmptyRoad(kFALSE),
 105 fFactorSAWindowSizes(1.),
 106 fFindV0s(kTRUE),
 107 fUseUnfoldingInClusterFinderSPD(kFALSE),
 108 fUseUnfoldingInClusterFinderSDD(kTRUE),
 109 fUseUnfoldingInClusterFinderSSD(kTRUE)
 110 {
 111   //
 112   // constructor
 113   //
 114   SetName("ITS");
 115   SetTitle("ITS");
 116
 117   SetLayersParameters();
 118   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
 119   SetUseTGeoInTracker(3);
 120   SetAllowSharedClusters(kTRUE);
 121   SetFindV0s(kTRUE);
 122   SetAddVirtualClustersInDeadZone(kFALSE);
 123   SetUseAmplitudeInfo(kTRUE);
 124   SetClusterErrorsParam(1);
 125 }
 126 //_____________________________________________________________________________
 127 AliITSRecoParam::~AliITSRecoParam()
 128 {
 129   //
 130   // destructor
 131   //
 132 }
 133 //_____________________________________________________________________________
 134 AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
 135 {
 136   //
 137   // make default reconstruction  parameters for hig  flux env.
 138   //
 139   AliITSRecoParam *param = new AliITSRecoParam();
 140
 141   param->fMaxSnp = 0.95;
 142
 143   param->fNSigmaYLayerForRoadY = 4.;
 144   param->fNSigmaRoadY = 7.5;
 145   param->fNSigmaZLayerForRoadZ = 4.;
 146   param->fNSigmaRoadZ = 7.5;
 147
 148   param->fNSigma2RoadZC = 60.; //7.75^2
 149   param->fNSigma2RoadYC = 60.; //7.75^2
 150   param->fNSigma2RoadZNonC = 50.; //7.07^2
 151   param->fNSigma2RoadYNonC = 50.; //7.07^2
 152
 153   param->fMaxChi2PerCluster[0] = 11.; //7
 154   param->fMaxChi2PerCluster[1] = 12.; //5
 155   param->fMaxChi2PerCluster[2] = 12.; //8
 156   param->fMaxChi2PerCluster[3] = 5.;  //8
 157   param->fMaxChi2PerCluster[4] = 12.; //6.5
 158
 159   param->fMaxNormChi2NonC[0] = 7.;
 160   param->fMaxNormChi2NonC[1] = 8.;
 161   param->fMaxNormChi2NonC[2] = 8.;
 162   param->fMaxNormChi2NonC[3] = 11.;
 163   param->fMaxNormChi2NonC[4] = 14.;
 164   param->fMaxNormChi2NonC[5] = 25.;
 165
 166   param->fMaxNormChi2C[0] = 11.;
 167   param->fMaxNormChi2C[1] = 13.;
 168   param->fMaxNormChi2C[2] = 15.;
 169   param->fMaxNormChi2C[3] = 18.;
 170   param->fMaxNormChi2C[4] = 30.;
 171   param->fMaxNormChi2C[5] = 35.;
 172
 173   param->fMaxNormChi2NonCForHypothesis = 7.;
 174
 175   param->fMaxChi2 = 35.;
 176
 177   param->fMaxChi2s[0] = 25.; //40
 178   param->fMaxChi2s[1] = 25.; //40
 179   param->fMaxChi2s[2] = 25.; //40
 180   param->fMaxChi2s[3] = 25.; //40
 181   param->fMaxChi2s[4] = 40.; //40
 182   param->fMaxChi2s[5] = 50.; //40
 183
 184   param->fMaxRoad = 6.;
 185
 186   // not used
 187   param->fMaxChi2In = 16.;
 188
 189   param->fMaxChi2sR[0] = 10.;
 190   param->fMaxChi2sR[1] = 10.;
 191   param->fMaxChi2sR[2] = 10.;
 192   param->fMaxChi2sR[3] = 10.;
 193   param->fMaxChi2sR[4] = 30.;
 194   param->fMaxChi2sR[5] = 40.;
 195
 196   param->fChi2PerCluster = 9.;
 197   // not used
 198
 199   param->fXV = 0.;
 200   param->fYV = 0.;
 201   param->fZV = 0.;
 202   param->fSigmaXV = 0.0050;
 203   param->fSigmaYV = 0.0050;
 204   param->fSigmaZV = 0.0100;
 205
 206   param->fVertexCut = 25.;
 207
 208   param->fMaxDZforPrimTrk = 0.4;
 209   param->fMaxDZToUseConstraint = 3.;
 210
 211   param->fMaxDforV0dghtrForProlongation = 30.;
 212   param->fMaxDForProlongation = 40.;//10.;
 213   param->fMaxDZForProlongation = 60.;//20.;
 214   param->fMinPtForProlongation = 0.120;
 215
 216   param->fZWindowDeadZone = 2.0;
 217   param->fSigmaXDeadZoneHit2 = 0.004/12.;
 218   param->fSigmaZDeadZoneHit2 = 0.001/12.;
 219   param->fXPassDeadZoneHits = 0.018;
 220
 221
 222   return param;
 223 }
 224 //_____________________________________________________________________________
 225 AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
 226 {
 227   //
 228   // make default reconstruction  parameters for low  flux env.
 229   //
 230   return GetHighFluxParam();
 231 }
 232 //_____________________________________________________________________________
 233 AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
 234 {
 235   //
 236   // make default reconstruction  parameters for cosmics
 237   //
 238   AliITSRecoParam *param = new AliITSRecoParam();
 239
 240   // larger seach windows for SA (in case of large misalignments)
 241   param->SetFactorSAWindowSizes(3.);
 242
 243   param->fMaxSnp = 0.95;
 244
 245   param->fNSigmaYLayerForRoadY = 4.;
 246   param->fNSigmaRoadY = 7.5;
 247   param->fNSigmaZLayerForRoadZ = 4.;
 248   param->fNSigmaRoadZ = 7.5;
 249
 250   param->fNSigma2RoadZC = 60.; //7.75^2
 251   param->fNSigma2RoadYC = 60.; //7.75^2
 252   param->fNSigma2RoadZNonC = 50.; //7.07^2
 253   param->fNSigma2RoadYNonC = 50.; //7.07^2
 254
 255   param->fMaxChi2PerCluster[0] = 11.; //7
 256   param->fMaxChi2PerCluster[1] = 12.; //5
 257   param->fMaxChi2PerCluster[2] = 12.; //8
 258   param->fMaxChi2PerCluster[3] = 5.;  //8
 259   param->fMaxChi2PerCluster[4] = 12.; //6.5
 260
 261   param->fMaxNormChi2NonC[0] = 7.;
 262   param->fMaxNormChi2NonC[1] = 8.;
 263   param->fMaxNormChi2NonC[2] = 8.;
 264   param->fMaxNormChi2NonC[3] = 11.;
 265   param->fMaxNormChi2NonC[4] = 14.;
 266   param->fMaxNormChi2NonC[5] = 25.;
 267
 268   param->fMaxNormChi2C[0] = 11.;
 269   param->fMaxNormChi2C[1] = 13.;
 270   param->fMaxNormChi2C[2] = 15.;
 271   param->fMaxNormChi2C[3] = 18.;
 272   param->fMaxNormChi2C[4] = 30.;
 273   param->fMaxNormChi2C[5] = 35.;
 274
 275   param->fMaxNormChi2NonCForHypothesis = 7.;
 276
 277   param->fMaxChi2 = 35.;
 278
 279   param->fMaxChi2s[0] = 25.; //40
 280   param->fMaxChi2s[1] = 25.; //40
 281   param->fMaxChi2s[2] = 25.; //40
 282   param->fMaxChi2s[3] = 25.; //40
 283   param->fMaxChi2s[4] = 40.; //40
 284   param->fMaxChi2s[5] = 50.; //40
 285
 286   param->fMaxRoad = 6.;
 287
 288   // not used
 289   param->fMaxChi2In = 16.;
 290
 291   param->fMaxChi2sR[0] = 10.;
 292   param->fMaxChi2sR[1] = 10.;
 293   param->fMaxChi2sR[2] = 10.;
 294   param->fMaxChi2sR[3] = 10.;
 295   param->fMaxChi2sR[4] = 30.;
 296   param->fMaxChi2sR[5] = 40.;
 297
 298   param->fChi2PerCluster = 9.;
 299   // not used
 300
 301   param->fXV = 0.;
 302   param->fYV = 0.;
 303   param->fZV = 0.;
 304   param->fSigmaXV = 0.0050;
 305   param->fSigmaYV = 0.0050;
 306   param->fSigmaZV = 0.0100;
 307
 308   param->fVertexCut = 25.;
 309
 310   param->fMaxDZforPrimTrk = 0.4;
 311   param->fMaxDZToUseConstraint = 3.;
 312
 313   param->fMaxDforV0dghtrForProlongation = 30.;
 314   param->fMaxDForProlongation = 10.;
 315   param->fMaxDZForProlongation = 20.;
 316   param->fMinPtForProlongation = 0.120;
 317
 318   param->fZWindowDeadZone = 2.0;
 319   param->fSigmaXDeadZoneHit2 = 0.004/12.;
 320   param->fSigmaZDeadZoneHit2 = 0.001/12.;
 321   param->fXPassDeadZoneHits = 0.018;
 322
 323
 324   return param;
 325 }
 326 //_____________________________________________________________________________
 327 AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
 328 {
 329   //
 330   // make special reconstruction parameters for Plane Efficiency study on layer i
 331   //
 332   AliITSRecoParam *param;
 333   param = GetHighFluxParam();
 334   param->SetComputePlaneEff();
 335   param->SetLayerToSkip(i);
 336   return param;
 337 }
 338 //_____________________________________________________________________________
 339 void AliITSRecoParam::SetLayersParameters()
 340 {
 341   //
 342   // number of layers and layers spatial resolutions
 343   //
 344
 345   // spatial resolutions of the detectors
 346   // y: 12 12 38 38 20 20 micron
 347   fSigmaY2[0]=1.44e-6;
 348   fSigmaY2[1]=1.44e-6;
 349   fSigmaY2[2]=1.444e-5;
 350   fSigmaY2[3]=1.444e-5;
 351   fSigmaY2[4]=4.0e-6;
 352   fSigmaY2[5]=4.0e-6;
 353   // z: 120 120 28 28 830 830 micron
 354   fSigmaZ2[0]=1.44e-4;
 355   fSigmaZ2[1]=1.44e-4;
 356   fSigmaZ2[2]=7.84e-6;
 357   fSigmaZ2[3]=7.84e-6;
 358   fSigmaZ2[4]=6.889e-3;
 359   fSigmaZ2[5]=6.889e-3;
 360
 361   return;
 362 }
 363