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