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