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