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 #include "AliESDV0Params.h"
  18 #include "AliLog.h"
  19
  20 /* $Id$ */
  21
  22 ///////////////////////////////////////////////////////////////////////////////
  23 //                                                                           //
  24 // Class with ITS reconstruction parameters                                  //
  25 // Origin: andrea.dainese@lnl.infn.it                                        //
  26 //                                                                           //
  27 ///////////////////////////////////////////////////////////////////////////////
  28
  29
  30
  31 ClassImp(AliITSRecoParam)
  32
  33 const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
  34 const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
  35 const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
  36 const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
  37 const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
  38
  39 //_____________________________________________________________________________
  40 AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
  41 fTracker(0),
  42 fITSonly(kFALSE),
  43 fVertexer(0),
  44 fClusterFinder(0),
  45 fPID(0),
  46 fVtxr3DZCutWide(0.),
  47 fVtxr3DRCutWide(0.),
  48 fVtxr3DZCutNarrow(0.),
  49 fVtxr3DRCutNarrow(0.),
  50 fVtxr3DPhiCutLoose(0.),
  51 fVtxr3DPhiCutTight(0.),
  52 fVtxr3DDCACut(0.),
  53 fVtxr3DPileupAlgo(1),
  54 fMaxSnp(1.),
  55 fNSigmaYLayerForRoadY(0),
  56 fNSigmaRoadY(0),
  57 fNSigmaZLayerForRoadZ(0),
  58 fNSigmaRoadZ(0),
  59 fNSigma2RoadZC(0),
  60 fNSigma2RoadYC(0),
  61 fNSigma2RoadZNonC(0),
  62 fNSigma2RoadYNonC(0),
  63 fRoadMisal(0),
  64 fMaxNormChi2NonCForHypothesis(0),
  65 fMaxChi2(0),
  66 fMaxRoad(0),
  67 fMaxChi2In(0),
  68 fChi2PerCluster(0),
  69 fSearchForExtras(kTRUE),
  70 fXV(0),
  71 fYV(0),
  72 fZV(0),
  73 fSigmaXV(0),
  74 fSigmaYV(0),
  75 fSigmaZV(0),
  76 fVertexCut(0),
  77 fMaxDZforPrimTrk(0),
  78 fMaxDZToUseConstraint(0),
  79 fMaxDforV0dghtrForProlongation(0),
  80 fMaxDForProlongation(0),
  81 fMaxDZForProlongation(0),
  82 fMinPtForProlongation(0),
  83 fAddVirtualClustersInDeadZone(kFALSE),
  84 fZWindowDeadZone(0),
  85 fSigmaXDeadZoneHit2(0),
  86 fSigmaZDeadZoneHit2(0),
  87 fXPassDeadZoneHits(0),
  88 fSkipSubdetsNotInTriggerCluster(kTRUE),
  89 fUseTGeoInTracker(3),
  90 fStepSizeTGeo(5.),
  91 fAllowSharedClusters(kTRUE),
  92 fClusterErrorsParam(1),
  93 fComputePlaneEff(kFALSE),
  94 fHistoPlaneEff(kFALSE),
  95 fUseTrackletsPlaneEff(kFALSE),
  96 fMCTrackletsPlaneEff(kFALSE),
  97 fBkgTrackletsPlaneEff(kFALSE),
  98 fTrackleterPhiWindowL1(0.10),
  99 fTrackleterPhiWindowL2(0.07),
 100 fTrackleterZetaWindowL1(0.6),
 101 fTrackleterZetaWindowL2(0.4),
 102 fUpdateOncePerEventPlaneEff(kTRUE),
 103 fMinContVtxPlaneEff(3),
 104 fIPlanePlaneEff(0),
 105 fReadPlaneEffFromOCDB(kFALSE),
 106 fMinPtPlaneEff(0),
 107 fMaxMissingClustersPlaneEff(5),
 108 fMaxMissingClustersOutPlaneEff(5),
 109 fRequireClusterInOuterLayerPlaneEff(kFALSE),
 110 fRequireClusterInInnerLayerPlaneEff(kFALSE),
 111 fOnlyConstraintPlaneEff(kFALSE),
 112 fNSigXFromBoundaryPlaneEff(1.),
 113 fNSigZFromBoundaryPlaneEff(1.),
 114 fImproveWithVertex(kFALSE),
 115 fExtendedEtaAcceptance(kFALSE),
 116 fUseBadZonesFromOCDB(kTRUE),
 117 fUseSingleBadChannelsFromOCDB(kFALSE),
 118 fMinFractionOfBadInRoad(0),
 119 fAllowProlongationWithEmptyRoad(kFALSE),
 120 fInwardFlagSA(kFALSE),
 121 fOuterStartLayerSA(2),
 122 fInnerStartLayerSA(3),
 123 fMinNPointsSA(3),
 124 fFactorSAWindowSizes(1.),
 125 fNLoopsSA(32),
 126 fMinPhiSA(0.002),
 127 fMaxPhiSA(0.0145),
 128 fMinLambdaSA(0.003),
 129 fMaxLambdaSA(0.008),
 130 fMinClusterChargeSA(0.),
 131 fSAOnePointTracks(kFALSE),
 132 fSAUseAllClusters(kFALSE),
 133 fMaxSPDcontrForSAToUseAllClusters(1000000),
 134 fSAUsedEdxInfo(kFALSE),
 135 fSelectBestMIP03(kFALSE),
 136 fFlagFakes(kFALSE),
 137 fUseImproveKalman(kFALSE),
 138 fFindV0s(kTRUE),
 139 fStoreLikeSignV0s(kFALSE),
 140 fUseUnfoldingInClusterFinderSPD(kFALSE),
 141 fUseUnfoldingInClusterFinderSDD(kTRUE),
 142 fUseUnfoldingInClusterFinderSSD(kFALSE),
 143 fUseBadChannelsInClusterFinderSSD(kFALSE),
 144 fUseSDDCorrectionMaps(kTRUE),
 145 fUseSDDClusterSizeSelection(kFALSE),
 146 fMinClusterChargeSDD(0.),
 147 fUseChargeMatchingInClusterFinderSSD(kTRUE),
 148 fTrackleterPhiWindow(0.08),
 149 fTrackleterThetaWindow(0.025),
 150 fTrackleterPhiShift(0.0045),
 151 fTrackleterRemoveClustersFromOverlaps(kFALSE),
 152 fTrackleterPhiOverlapCut(0.005),
 153 fTrackleterZetaOverlapCut(0.05),
 154 fTrackleterPhiRotationAngle(0.0),
 155 fTrackleterNStdDev(1.5),
 156 fScaleDTBySin2T(kFALSE),
 157 fUseCosmicRunShiftsSSD(kFALSE),
 158 fSPDRemoveNoisyFlag(kTRUE),
 159 fSPDRemoveDeadFlag(kTRUE),
 160 fVertexerFastSmearX(0.005),
 161 fVertexerFastSmearY(0.005),
 162 fVertexerFastSmearZ(0.01),
 163 fAlignFilterCosmics(kFALSE),
 164 fAlignFilterCosmicMergeTracks(kTRUE),
 165 fAlignFilterMinITSPoints(4),
 166 fAlignFilterMinITSPointsMerged(4),
 167 fAlignFilterOnlyITSSATracks(kTRUE),
 168 fAlignFilterOnlyITSTPCTracks(kFALSE),
 169 fAlignFilterSkipExtra(kFALSE),
 170 fAlignFilterMaxMatchingAngle(0.085),
 171 fAlignFilterMinAngleWrtModulePlanes(0.52),
 172 fAlignFilterMinPt(0.),
 173 fAlignFilterMaxPt(1.e10),
 174 fAlignFilterFillQANtuples(kTRUE),
 175 //
 176 fMultCutPxDrSPDin(0.1),
 177 fMultCutPxDrSPDout(0.15),
 178 fMultCutPxDz(0.2),
 179 fMultCutDCArz(0.5),
 180 fMultCutMinElectronProbTPC(0.5),
 181 fMultCutMinElectronProbESD(0.1),
 182 fMultCutMinP(0.05),
 183 fMultCutMinRGamma(2.),
 184 fMultCutMinRK0(1.),
 185 fMultCutMinPointAngle(0.98),
 186 fMultCutMaxDCADauther(0.5),
 187 fMultCutMassGamma(0.03),
 188 fMultCutMassGammaNSigma(5.),
 189 fMultCutMassK0(0.03),
 190 fMultCutMassK0NSigma(5.),
 191 fMultCutChi2cGamma(2.),
 192 fMultCutChi2cK0(2.),
 193 fMultCutGammaSFromDecay(-10.),
 194 fMultCutK0SFromDecay(-10.),
 195 fMultCutMaxDCA(1.),
 196 //
 197 fCorrectLorentzAngleSPD(kTRUE),
 198 fTanLorentzAngleHolesSPD(0.017455), // tan(1 degree)
 199 fCorrectLorentzAngleSSD(kTRUE),
 200 fTanLorentzAngleHolesSSD(0.016),  // tan(0.94 degrees)
 201 fTanLorentzAngleElectronsSSD(0.068), // tan(3.98 degrees)
 202 //
 203 fESDV0Params(NULL)
 204 {
 205   //
 206   // constructor
 207   //
 208   SetName("ITS");
 209   SetTitle("ITS");
 210
 211   SetLayersParameters();
 212   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
 213   SetUseTGeoInTracker(3);
 214   SetStepSizeTGeo(5.);
 215   SetAllowSharedClusters(kTRUE);
 216   SetFindV0s(kTRUE);
 217   SetAddVirtualClustersInDeadZone(kFALSE);
 218   SetUseAmplitudeInfo(kTRUE);
 219   SetClusterErrorsParam(1);
 220   SetClusterMisalError(0.);
 221   SetClusterMisalErrorBOn(0.);
 222   SetVertexer3DDefaults();
 223
 224   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fAlignFilterUseLayer[i]=kTRUE;
 225   fESDV0Params = new AliESDV0Params();
 226
 227   fESDV0Params->SetMaxDist0(0.1);
 228   fESDV0Params->SetMaxDist1(0.1);
 229   fESDV0Params->SetMaxDist(1.);
 230   fESDV0Params->SetMinPointAngle(0.85);
 231   fESDV0Params->SetMinPointAngle2(0.99);
 232   fESDV0Params->SetMinR(0.5);
 233   fESDV0Params->SetMaxR(220.);
 234   fESDV0Params->SetMinPABestConst(0.9999);
 235   fESDV0Params->SetMaxRBestConst(10.);
 236   fESDV0Params->SetCausality0Cut(0.19);
 237   fESDV0Params->SetLikelihood01Cut(0.45);
 238   fESDV0Params->SetLikelihood1Cut(0.5);
 239   fESDV0Params->SetCombinedCut(0.55);
 240   fESDV0Params->SetMinClFullTrk(5.0);
 241   fESDV0Params->SetMinTgl0(1.05);
 242   fESDV0Params->SetMinRTgl0(40.0);
 243
 244   fESDV0Params->SetMinNormDistForbTgl0(3.0);
 245   fESDV0Params->SetMinClForb0(4.5);
 246   fESDV0Params->SetMinNormDistForb1(3.0);
 247   fESDV0Params->SetMinNormDistForb2(2.0);
 248   fESDV0Params->SetMinNormDistForb3(1.0);
 249   fESDV0Params->SetMinNormDistForb4(4.0);
 250   fESDV0Params->SetMinNormDistForb5(5.0);
 251   fESDV0Params->SetMinNormDistForbProt(2.0);
 252   fESDV0Params->SetMaxPidProbPionForb(0.5);
 253
 254   fESDV0Params->SetMinRTPCdensity(40.);
 255   fESDV0Params->SetMaxRTPCdensity0(110.);
 256   fESDV0Params->SetMaxRTPCdensity10(120.);
 257   fESDV0Params->SetMaxRTPCdensity20(130.);
 258   fESDV0Params->SetMaxRTPCdensity30(140.);
 259
 260   fESDV0Params->SetMinTPCdensity(0.6);
 261   fESDV0Params->SetMinTgl1(1.1);
 262   fESDV0Params->SetMinTgl2(1.0);
 263   fESDV0Params->SetMinchi2before0(16.);
 264   fESDV0Params->SetMinchi2before1(16.);
 265   fESDV0Params->SetMinchi2after0(16.);
 266   fESDV0Params->SetMinchi2after1(16.);
 267   fESDV0Params->SetAddchi2SharedCl(18.);
 268   fESDV0Params->SetAddchi2NegCl0(25.);
 269   fESDV0Params->SetAddchi2NegCl1(30.);
 270   fESDV0Params->SetSigp0Par0(0.0001);
 271   fESDV0Params->SetSigp0Par1(0.001);
 272   fESDV0Params->SetSigp0Par2(0.1);
 273   fESDV0Params->SetSigpPar0(0.5);
 274   fESDV0Params->SetSigpPar1(0.6);
 275   fESDV0Params->SetSigpPar2(0.4);
 276   fESDV0Params->SetMaxDcaLh0(0.5);
 277   fESDV0Params->SetStreamLevel(0);
 278   fESDV0Params->SetChi2KF(100);
 279   fESDV0Params->SetRobustChi2KF(100);
 280
 281 }
 282 //_____________________________________________________________________________
 283 AliITSRecoParam::~AliITSRecoParam()
 284 {
 285   //
 286   // destructor
 287   //
 288   if(fESDV0Params){
 289     delete fESDV0Params;
 290     fESDV0Params=NULL;
 291   }
 292 }
 293 //_____________________________________________________________________________
 294 AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
 295 {
 296   //
 297   // make default reconstruction  parameters for hig  flux env.
 298   //
 299   AliITSRecoParam *param = new AliITSRecoParam();
 300   param->SetVertexer3DDefaults();
 301   param->SetSPDVertexerPileupAlgoOff();
 302   // use of bads from OCDB
 303   param->SetUseBadZonesFromOCDB(kTRUE);
 304   param->SetUseSingleBadChannelsFromOCDB(kFALSE);
 305   // use pointing to vertex during prolongation
 306   param->SetImproveWithVertex(kTRUE);
 307   // extended eta acceptance
 308   param->SetExtendedEtaAcceptance(kFALSE);
 309   // allow to skip layer if no cluster and no bad
 310   param->SetAllowProlongationWithEmptyRoad(kFALSE);
 311   // set event specie
 312   param->SetEventSpecie(AliRecoParam::kHighMult);
 313
 314   param->fMaxSnp = 0.95;
 315
 316   param->fNSigmaYLayerForRoadY = 4.;
 317   param->fNSigmaRoadY = 7.5;
 318   param->fNSigmaZLayerForRoadZ = 4.;
 319   param->fNSigmaRoadZ = 7.5;
 320
 321   param->fNSigma2RoadZC = 60.; //7.75^2
 322   param->fNSigma2RoadYC = 60.; //7.75^2
 323   param->fNSigma2RoadZNonC = 50.; //7.07^2
 324   param->fNSigma2RoadYNonC = 50.; //7.07^2
 325
 326   param->fMaxChi2PerCluster[0] = 11.; //7
 327   param->fMaxChi2PerCluster[1] = 12.; //5
 328   param->fMaxChi2PerCluster[2] = 12.; //8
 329   param->fMaxChi2PerCluster[3] = 5.;  //8
 330   param->fMaxChi2PerCluster[4] = 12.; //6.5
 331
 332   param->fMaxNormChi2NonC[0] = 7.;
 333   param->fMaxNormChi2NonC[1] = 8.;
 334   param->fMaxNormChi2NonC[2] = 8.;
 335   param->fMaxNormChi2NonC[3] = 11.;
 336   param->fMaxNormChi2NonC[4] = 14.;
 337   param->fMaxNormChi2NonC[5] = 25.;
 338
 339   param->fMaxNormChi2C[0] = 11.;
 340   param->fMaxNormChi2C[1] = 13.;
 341   param->fMaxNormChi2C[2] = 15.;
 342   param->fMaxNormChi2C[3] = 18.;
 343   param->fMaxNormChi2C[4] = 30.;
 344   param->fMaxNormChi2C[5] = 35.;
 345
 346   param->fMaxNormChi2NonCForHypothesis = 7.;
 347
 348   param->fMaxChi2 = 35.;
 349
 350   param->fMaxChi2s[0] = 25.; //40
 351   param->fMaxChi2s[1] = 25.; //40
 352   param->fMaxChi2s[2] = 25.; //40
 353   param->fMaxChi2s[3] = 25.; //40
 354   param->fMaxChi2s[4] = 40.; //40
 355   param->fMaxChi2s[5] = 50.; //40
 356
 357   param->fMaxRoad = 6.;
 358
 359   // not used
 360   param->fMaxChi2In = 16.;
 361
 362   param->fMaxChi2sR[0] = 10.;
 363   param->fMaxChi2sR[1] = 10.;
 364   param->fMaxChi2sR[2] = 10.;
 365   param->fMaxChi2sR[3] = 10.;
 366   param->fMaxChi2sR[4] = 30.;
 367   param->fMaxChi2sR[5] = 40.;
 368
 369
 370   param->fChi2PerCluster = 9.;
 371   // not used
 372   param->fSearchForExtras = kFALSE;
 373
 374   param->fXV = 0.;
 375   param->fYV = 0.;
 376   param->fZV = 0.;
 377   param->fSigmaXV = 0.0050;
 378   param->fSigmaYV = 0.0050;
 379   param->fSigmaZV = 0.0100;
 380
 381   param->fVertexCut = 25.;
 382
 383   param->fMaxDZforPrimTrk = 0.4;
 384   param->fMaxDZToUseConstraint = 3.;
 385
 386   param->fMaxDforV0dghtrForProlongation = 30.;
 387   param->fMaxDForProlongation = 40.;//10.;
 388   param->fMaxDZForProlongation = 60.;//20.;
 389   param->fMinPtForProlongation = 0.120;
 390
 391   param->fZWindowDeadZone = 2.0;
 392   param->fSigmaXDeadZoneHit2 = 0.004/12.;
 393   param->fSigmaZDeadZoneHit2 = 0.001/12.;
 394   param->fXPassDeadZoneHits = 0.018;
 395   //
 396   // mult reco
 397   param->fMultCutPxDrSPDin = 0.1;
 398   param->fMultCutPxDrSPDout = 0.15;
 399   param->fMultCutPxDz = 0.2;
 400   param->fMultCutDCArz = 0.5;
 401   param->fMultCutMinElectronProbTPC = 0.5;
 402   param->fMultCutMinElectronProbESD = 0.1;
 403   param->fMultCutMinP = 0.05;
 404   param->fMultCutMinRGamma = 2.;
 405   param->fMultCutMinRK0 = 1.;
 406   param->fMultCutMinPointAngle = 0.98;
 407   param->fMultCutMaxDCADauther = 0.5;
 408   param->fMultCutMassGamma = 0.03;
 409   param->fMultCutMassGammaNSigma = 5.;
 410   param->fMultCutMassK0 = 0.03;
 411   param->fMultCutMassK0NSigma = 5.;
 412   param->fMultCutChi2cGamma = 2.;
 413   param->fMultCutChi2cK0 = 2.;
 414   param->fMultCutGammaSFromDecay = -10.;
 415   param->fMultCutK0SFromDecay = -10.;
 416   param->fMultCutMaxDCA = 1.;
 417   //
 418   // trackleter
 419   param->fTrackleterPhiWindow = 0.06;
 420   param->fScaleDTBySin2T = kTRUE;
 421   //
 422   param->fSelectBestMIP03 = kFALSE;//kTRUE;
 423   param->fFlagFakes       = kTRUE;
 424   param->fUseImproveKalman= kFALSE;
 425   //
 426   return param;
 427 }
 428 //_____________________________________________________________________________
 429 AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
 430 {
 431   //
 432   // make default reconstruction  parameters for low  flux env.
 433   //
 434   AliITSRecoParam *param = new AliITSRecoParam();
 435   param->SetVertexer3DDefaults();
 436
 437   // full use of bads from OCDB
 438   param->SetUseBadZonesFromOCDB(kTRUE);
 439   param->SetUseSingleBadChannelsFromOCDB(kTRUE);
 440   // extended eta acceptance
 441   param->SetExtendedEtaAcceptance(kTRUE);
 442   // allow to skip layer if no cluster and no bad
 443   param->SetAllowProlongationWithEmptyRoad(kTRUE);
 444   // set event specie
 445   param->SetEventSpecie(AliRecoParam::kLowMult);
 446
 447   param->fMaxSnp = 0.95;
 448
 449   param->fNSigmaYLayerForRoadY = 4.;
 450   param->fNSigmaRoadY = 7.5;
 451   param->fNSigmaZLayerForRoadZ = 4.;
 452   param->fNSigmaRoadZ = 7.5;
 453
 454   param->fNSigma2RoadZC = 60.; //7.75^2
 455   param->fNSigma2RoadYC = 60.; //7.75^2
 456   param->fNSigma2RoadZNonC = 50.; //7.07^2
 457   param->fNSigma2RoadYNonC = 50.; //7.07^2
 458
 459   param->fMaxChi2PerCluster[0] = 11.; //7
 460   param->fMaxChi2PerCluster[1] = 12.; //5
 461   param->fMaxChi2PerCluster[2] = 12.; //8
 462   param->fMaxChi2PerCluster[3] = 5.;  //8
 463   param->fMaxChi2PerCluster[4] = 12.; //6.5
 464
 465   param->fMaxNormChi2NonC[0] = 7.;
 466   param->fMaxNormChi2NonC[1] = 8.;
 467   param->fMaxNormChi2NonC[2] = 8.;
 468   param->fMaxNormChi2NonC[3] = 11.;
 469   param->fMaxNormChi2NonC[4] = 14.;
 470   param->fMaxNormChi2NonC[5] = 25.;
 471
 472   param->fMaxNormChi2C[0] = 11.;
 473   param->fMaxNormChi2C[1] = 13.;
 474   param->fMaxNormChi2C[2] = 15.;
 475   param->fMaxNormChi2C[3] = 18.;
 476   param->fMaxNormChi2C[4] = 30.;
 477   param->fMaxNormChi2C[5] = 35.;
 478
 479   param->fMaxNormChi2NonCForHypothesis = 7.;
 480
 481   param->fMaxChi2 = 35.;
 482
 483   param->fMaxChi2s[0] = 25.; //40
 484   param->fMaxChi2s[1] = 25.; //40
 485   param->fMaxChi2s[2] = 25.; //40
 486   param->fMaxChi2s[3] = 25.; //40
 487   param->fMaxChi2s[4] = 40.; //40
 488   param->fMaxChi2s[5] = 50.; //40
 489
 490   param->fMaxRoad = 6.;
 491
 492   // not used
 493   param->fMaxChi2In = 16.;
 494
 495   param->fMaxChi2sR[0] = 10.;
 496   param->fMaxChi2sR[1] = 10.;
 497   param->fMaxChi2sR[2] = 10.;
 498   param->fMaxChi2sR[3] = 10.;
 499   param->fMaxChi2sR[4] = 30.;
 500   param->fMaxChi2sR[5] = 40.;
 501
 502   param->fChi2PerCluster = 9.;
 503   // not used
 504
 505   param->fXV = 0.;
 506   param->fYV = 0.;
 507   param->fZV = 0.;
 508   param->fSigmaXV = 0.0050;
 509   param->fSigmaYV = 0.0050;
 510   param->fSigmaZV = 0.0100;
 511
 512   param->fVertexCut = 25.;
 513
 514   param->fMaxDZforPrimTrk = 0.4;
 515   param->fMaxDZToUseConstraint = 3.;
 516
 517   param->fMaxDforV0dghtrForProlongation = 30.;
 518   param->fMaxDForProlongation = 40.;//10.;
 519   param->fMaxDZForProlongation = 60.;//20.;
 520   param->fMinPtForProlongation = 0.120;
 521
 522   param->fZWindowDeadZone = 2.0;
 523   param->fSigmaXDeadZoneHit2 = 0.004/12.;
 524   param->fSigmaZDeadZoneHit2 = 0.001/12.;
 525   param->fXPassDeadZoneHits = 0.018;
 526   param->SetNLoopsSA(20);
 527   param->fMaxPhiSA = 0.07;
 528   param->fMaxLambdaSA = 0.04;
 529
 530
 531   param->GetESDV0Params()->SetMinRTgl0(0.5);
 532   param->GetESDV0Params()->SetMaxDist(1.5);
 533   param->GetESDV0Params()->SetMaxDcaLh0(1.5);
 534   param->GetESDV0Params()->SetMaxRBestConst(80);
 535   param->GetESDV0Params()->SetMinPABestConst(0.99);
 536   param->GetESDV0Params()->SetMinNormDistForbTgl0(1.);
 537   param->GetESDV0Params()->SetMinNormDistForb1(2.);
 538   param->GetESDV0Params()->SetMinNormDistForbProt(1.);
 539   param->GetESDV0Params()->SetMaxPidProbPionForb(0.7);
 540   param->GetESDV0Params()->SetLikelihood01Cut(0.3);
 541   param->GetESDV0Params()->SetLikelihood1Cut(0.35);
 542   param->GetESDV0Params()->SetCombinedCut(0.4);
 543
 544   // mult reco
 545   param->fMultCutPxDrSPDin = 0.1;
 546   param->fMultCutPxDrSPDout = 0.15;
 547   param->fMultCutPxDz = 0.2;
 548   param->fMultCutDCArz = 0.5;
 549   param->fMultCutMinElectronProbTPC = 0.5;
 550   param->fMultCutMinElectronProbESD = 0.1;
 551   param->fMultCutMinP = 0.05;
 552   param->fMultCutMinRGamma = 2.;
 553   param->fMultCutMinRK0 = 1.;
 554   param->fMultCutMinPointAngle = 0.98;
 555   param->fMultCutMaxDCADauther = 0.5;
 556   param->fMultCutMassGamma = 0.03;
 557   param->fMultCutMassGammaNSigma = 5.;
 558   param->fMultCutMassK0 = 0.03;
 559   param->fMultCutMassK0NSigma = 5.;
 560   param->fMultCutChi2cGamma = 2.;
 561   param->fMultCutChi2cK0 = 2.;
 562   param->fMultCutGammaSFromDecay = -10.;
 563   param->fMultCutK0SFromDecay = -10.;
 564   param->fMultCutMaxDCA = 1.;
 565   //
 566
 567   return param;
 568 }
 569 //_____________________________________________________________________________
 570 AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
 571 {
 572   //
 573   // make default reconstruction  parameters for cosmics
 574   //
 575   AliITSRecoParam *param = new AliITSRecoParam();
 576
 577   // vertexer for cosmics
 578   param->SetVertexer(2);
 579
 580   param->SetClusterErrorsParam(2);
 581   param->SetFindV0s(kFALSE);
 582   param->SetAddVirtualClustersInDeadZone(kFALSE);
 583   param->SetUseAmplitudeInfo(kFALSE);
 584
 585   // set event specie
 586   param->SetEventSpecie(AliRecoParam::kCosmic);
 587
 588   // full use of bads from OCDB
 589   param->SetUseBadZonesFromOCDB(kTRUE);
 590   param->SetUseSingleBadChannelsFromOCDB(kTRUE);
 591
 592   // find independently ITS SA tracks
 593   param->SetSAUseAllClusters();
 594   param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
 595
 596   //****** COSMICS 2009 (same as COSMICS 2008) *********************
 597
 598   // to maximize efficiency
 599   param->SetAllowProlongationWithEmptyRoad();
 600   param->SetMinNPointsSA(2);
 601
 602   // larger seach windows for SA (in case of large misalignments)
 603   param->SetNLoopsSA(32);
 604   param->SetFactorSAWindowSizes(20);
 605
 606   // additional error due to misal (B off)
 607   param->SetClusterMisalErrorY(1.0,1.0,1.0,1.0,1.0,1.0); // [cm]
 608   param->SetClusterMisalErrorZ(1.0,1.0,1.0,1.0,1.0,1.0); // [cm]
 609   // additional error due to misal (B on)
 610   param->SetClusterMisalErrorYBOn(0.0,0.0,0.1,0.1,0.1,0.1); // [cm]
 611   param->SetClusterMisalErrorZBOn(0.1,0.1,0.1,0.1,0.1,0.1); // [cm]
 612
 613
 614   // SDD configuration
 615   param->fUseSDDCorrectionMaps = kFALSE;
 616   param->fUseSDDClusterSizeSelection=kTRUE;
 617   param->fMinClusterChargeSDD=30.;
 618
 619
 620   // alignment data filter
 621   param->SetAlignFilterCosmics(kTRUE);
 622   param->SetAlignFilterCosmicMergeTracks(kTRUE);
 623   param->SetAlignFilterMinITSPoints(1);
 624   param->SetAlignFilterMinITSPointsMerged(3);
 625   param->SetAlignFilterOnlyITSSATracks(kTRUE);
 626   param->SetAlignFilterOnlyITSTPCTracks(kFALSE);
 627   param->SetAlignFilterSkipExtra(kFALSE);
 628   param->SetAlignFilterMaxMatchingAngle(0.085/*5deg*/);
 629   param->SetAlignFilterMinPt(0.2);
 630   param->SetAlignFilterMaxPt(1.e10);
 631   param->SetAlignFilterFillQANtuples(kTRUE);
 632
 633   //******************************************************************
 634
 635   param->fMaxSnp = 0.95;
 636
 637   param->fNSigmaYLayerForRoadY = 4.;
 638   param->fNSigmaRoadY = 7.5;
 639   param->fNSigmaZLayerForRoadZ = 4.;
 640   param->fNSigmaRoadZ = 7.5;
 641
 642   param->fNSigma2RoadZC = 60.; //7.75^2
 643   param->fNSigma2RoadYC = 60.; //7.75^2
 644   param->fNSigma2RoadZNonC = 50.; //7.07^2
 645   param->fNSigma2RoadYNonC = 50.; //7.07^2
 646
 647   param->fMaxChi2PerCluster[0] = 11.; //7
 648   param->fMaxChi2PerCluster[1] = 12.; //5
 649   param->fMaxChi2PerCluster[2] = 12.; //8
 650   param->fMaxChi2PerCluster[3] = 5.;  //8
 651   param->fMaxChi2PerCluster[4] = 12.; //6.5
 652
 653   param->fMaxNormChi2NonC[0] = 7.;
 654   param->fMaxNormChi2NonC[1] = 8.;
 655   param->fMaxNormChi2NonC[2] = 8.;
 656   param->fMaxNormChi2NonC[3] = 11.;
 657   param->fMaxNormChi2NonC[4] = 14.;
 658   param->fMaxNormChi2NonC[5] = 25.;
 659
 660   param->fMaxNormChi2C[0] = 11.;
 661   param->fMaxNormChi2C[1] = 13.;
 662   param->fMaxNormChi2C[2] = 15.;
 663   param->fMaxNormChi2C[3] = 18.;
 664   param->fMaxNormChi2C[4] = 30.;
 665   param->fMaxNormChi2C[5] = 35.;
 666
 667   param->fMaxNormChi2NonCForHypothesis = 7.;
 668
 669   param->fMaxChi2 = 35.;
 670
 671   param->fMaxChi2s[0] = 25.; //40
 672   param->fMaxChi2s[1] = 25.; //40
 673   param->fMaxChi2s[2] = 25.; //40
 674   param->fMaxChi2s[3] = 25.; //40
 675   param->fMaxChi2s[4] = 40.; //40
 676   param->fMaxChi2s[5] = 50.; //40
 677
 678   param->fMaxRoad = 6.;
 679
 680   // not used
 681   param->fMaxChi2In = 16.;
 682
 683   param->fMaxChi2sR[0] = 10.;
 684   param->fMaxChi2sR[1] = 10.;
 685   param->fMaxChi2sR[2] = 10.;
 686   param->fMaxChi2sR[3] = 10.;
 687   param->fMaxChi2sR[4] = 30.;
 688   param->fMaxChi2sR[5] = 40.;
 689
 690   param->fChi2PerCluster = 9.;
 691   // not used
 692
 693   param->fXV = 0.;
 694   param->fYV = 0.;
 695   param->fZV = 0.;
 696   param->fSigmaXV = 0.0050;
 697   param->fSigmaYV = 0.0050;
 698   param->fSigmaZV = 0.0100;
 699
 700   param->fVertexCut = 25.;
 701
 702   param->fMaxDZforPrimTrk = 0.4;
 703   param->fMaxDZToUseConstraint = 3.;
 704
 705   param->fMaxDforV0dghtrForProlongation = 30.;
 706   param->fMaxDForProlongation = 10.;
 707   param->fMaxDZForProlongation = 20.;
 708   param->fMinPtForProlongation = 0.120;
 709
 710   param->fZWindowDeadZone = 2.0;
 711   param->fSigmaXDeadZoneHit2 = 0.004/12.;
 712   param->fSigmaZDeadZoneHit2 = 0.001/12.;
 713   param->fXPassDeadZoneHits = 0.018;
 714
 715   return param;
 716 }
 717 //_____________________________________________________________________________
 718 AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
 719 {
 720   //
 721   // make special reconstruction parameters for Plane Efficiency study on layer i (0,5)
 722   //
 723   // if i=-1, then the evaluation for both pixel layers is tried with the tracklet method
 724   //
 725  if (i<-1 || i>=AliITSgeomTGeo::kNLayers) {
 726     printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane nor in the range [0,5] neither =-1\n");
 727     printf("returning null pointer");
 728     return NULL;
 729   }
 730   if(i>=0) {  // Method using tracks (remove given plane from tracking)
 731     AliITSRecoParam *param;
 732     param = GetLowFluxParam();
 733     param->SetClusterErrorsParam(2);
 734     // find independently ITS SA tracks
 735     param->SetSAUseAllClusters();
 736     param->SetOuterStartLayerSA(2);
 737     param->SetAllowProlongationWithEmptyRoad(kTRUE);
 738     // larger seach windows for SA (in case of large misalignments)
 739     param->SetFactorSAWindowSizes(2);
 740
 741     // Misalignment syst errors decided at ITS meeting 25.03.2010
 742     // additional error due to misal (B off)
 743     param->SetClusterMisalErrorY(0.0010,0.0010,0.0300,0.0300,0.0020,0.0020); // [cm]
 744     param->SetClusterMisalErrorZ(0.0100,0.0100,0.0100,0.0100,0.0500,0.0500); // [cm]
 745     // additional error due to misal (B on)
 746     param->SetClusterMisalErrorYBOn(0.0010,0.0030,0.0500,0.0500,0.0020,0.0020); // [cm]
 747     param->SetClusterMisalErrorZBOn(0.0100,0.0100,0.0100,0.0100,0.0500,0.0500); // [cm]
 748     //----
 749
 750     // SDD configuration
 751     param->SetUseSDDCorrectionMaps(kTRUE); // changed 30.04.2010
 752     param->SetUseSDDClusterSizeSelection(kTRUE);
 753     param->SetMinClusterChargeSDD(30.);
 754     param->SetUseUnfoldingInClusterFinderSDD(kFALSE);
 755     //
 756     param->SetComputePlaneEff();
 757     param->SetLayerToSkip(i);
 758     param->SetIPlanePlaneEff(i);
 759     param->fNSigXFromBoundaryPlaneEff= 1.;
 760     param->fNSigZFromBoundaryPlaneEff= 1.;
 761     // optimized setting for SPD0 (i==0)
 762     if (i==0) {
 763       param->fMinPtPlaneEff = 0.200; // high pt particles
 764       param->fMaxMissingClustersPlaneEff = 2; // at most 2 layers out of 5 without cluster
 765       param->fMaxMissingClustersOutPlaneEff = 2; // at most 2 layers out of 5 external ones without cluster
 766       param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
 767       //param->fOnlyConstraintPlaneEff = kTRUE;
 768     }
 769     if (i==1 ) {
 770       param->fMinPtPlaneEff = 0.200; // high pt particles
 771       param->fMaxMissingClustersPlaneEff = 2; // at most 2 layer out of 5 without cluster
 772       param->fMaxMissingClustersOutPlaneEff = 2; // at most 2 layer out of 4 external ones without cluster
 773       //param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SSD1 must be
 774     }
 775     if (i==2) {
 776       param->fMinPtPlaneEff = 0.200; // high pt particles
 777       param->fMaxMissingClustersPlaneEff = 2; // at most 2 layer out of 5 without cluster
 778       param->fMaxMissingClustersOutPlaneEff = 2; // at most 2 layer out of 3 external ones without cluster
 779       //param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
 780       //param->fOnlyConstraintPlaneEff = kTRUE;
 781     }
 782     if (i==3) {
 783       param->fMinPtPlaneEff = 0.200; // high pt particles
 784       param->fMaxMissingClustersPlaneEff = 2; // at most 2 layer out of 5 without cluster
 785       param->fMaxMissingClustersOutPlaneEff = 1; // at most 1 layer out of 2 external ones without cluster
 786     }
 787     if (i==4) {
 788       param->fMinPtPlaneEff = 0.200; // high pt particles
 789       param->fMaxMissingClustersPlaneEff = 2; // at most 2 layer out of 5 without cluster
 790       param->fMaxMissingClustersOutPlaneEff = 1; // at most 1 layer out of 1 external ones without cluster
 791       // param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
 792       //param->fOnlyConstraintPlaneEff = kTRUE;
 793     }
 794     if (i==5) {
 795       param->fMinPtPlaneEff = 0.200; // high pt particles
 796       param->fMaxMissingClustersPlaneEff = 2; // at most 2 layer out of 5 without cluster
 797     }
 798     //
 799     return param;
 800   }
 801   else if (i==-1) { // Method using tracklets
 802     AliITSRecoParam *param;
 803     param = GetLowFluxParam();
 804     param->SetIPlanePlaneEff(i);
 805     param->SetComputePlaneEff(kTRUE,kFALSE);
 806     param->SetUseTrackletsPlaneEff(kTRUE);
 807     param->SetTrackleterPhiWindowL2(0.07);
 808     param->SetTrackleterZetaWindowL2(0.4);
 809     param->SetTrackleterPhiWindowL1(0.10);
 810     param->SetTrackleterZetaWindowL1(0.6);
 811     param->SetUpdateOncePerEventPlaneEff(kTRUE);
 812     param->SetMinContVtxPlaneEff(3);
 813     return param;
 814   }
 815   else {
 816     AliErrorGeneral("AliITSRecoParam",Form("Unrecognised value of i %d\n",i));
 817     return 0;
 818   }
 819 }
 820
 821 //_____________________________________________________________________________
 822 void AliITSRecoParam::SetLayersParameters()
 823 {
 824   //
 825   // number of layers and layers spatial resolutions
 826   //
 827
 828   // spatial resolutions of the detectors
 829   // y: 12 12 38 38 20 20 micron
 830   fSigmaY2[0]=1.44e-6;
 831   fSigmaY2[1]=1.44e-6;
 832   fSigmaY2[2]=1.444e-5;
 833   fSigmaY2[3]=1.444e-5;
 834   fSigmaY2[4]=4.0e-6;
 835   fSigmaY2[5]=4.0e-6;
 836   // z: 120 120 28 28 830 830 micron
 837   fSigmaZ2[0]=1.44e-4;
 838   fSigmaZ2[1]=1.44e-4;
 839   fSigmaZ2[2]=7.84e-6;
 840   fSigmaZ2[3]=7.84e-6;
 841   fSigmaZ2[4]=6.889e-3;
 842   fSigmaZ2[5]=6.889e-3;
 843
 844   return;
 845 }
 846 //_____________________________________________________________________________
 847 void AliITSRecoParam::PrintParameters() const
 848 {
 849   //
 850   // print parameters
 851   //
 852
 853   printf("=============================  AliITSRecoParam::PrintParameters ");
 854   printf("============================= \n\n");
 855   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
 856     if(!fLayersToSkip[i]) {
 857       printf("ITS Traking: using layer %d\n",i);
 858     } else {
 859       printf("ITS Traking: skipping layer %d\n",i);
 860     }
 861   }
 862   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
 863     if(fUseAmplitudeInfo[i]) {
 864       printf("ITS Traking: use amplitude info for layer %d\n",i);
 865     } else {
 866       printf("ITS Traking: don't use amplitude info for layer %d\n",i);
 867     }
 868   }
 869   for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
 870     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]);
 871
 872
 873   Dump();
 874
 875   return;
 876 }
 877
 878