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