X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSRecoParam.cxx;h=d9db7884c8a92494434eb6ce65f445347d808e98;hb=d934d2578e7cb4b77f8144a7153aa6e478f64809;hp=05e20f6b2bc73e02bb00132fe40b30a38725a8dc;hpb=afd25725ca399dedf6e4866e4bbd8589c16ff257;p=u%2Fmrichter%2FAliRoot.git

diff --git a/ITS/AliITSRecoParam.cxx b/ITS/AliITSRecoParam.cxx
index 05e20f6b2bc..d9db7884c8a 100644
--- a/ITS/AliITSRecoParam.cxx
+++ b/ITS/AliITSRecoParam.cxx
@@ -1,5 +1,5 @@
 /**************************************************************************
- * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
+ * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
  * Author: The ALICE Off-line Project.                                    *
  * Contributors are mentioned in the code where appropriate.              *
@@ -13,6 +13,9 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
+#include "AliITSRecoParam.h"
+
+/* $Id$ */
 
 ///////////////////////////////////////////////////////////////////////////////
 //                                                                           //
@@ -22,22 +25,145 @@
 ///////////////////////////////////////////////////////////////////////////////
 
 
-#include "AliITSRecoParam.h"
 
 ClassImp(AliITSRecoParam)
 
+const Int_t AliITSRecoParam::fgkLayersNotToSkip[AliITSgeomTGeo::kNLayers]={0,0,0,0,0,0};
+const Int_t AliITSRecoParam::fgkLastLayerToTrackTo=0;
+const Int_t AliITSRecoParam::fgkMaxDetectorPerLayer=1000;
+const Double_t AliITSRecoParam::fgkriw=80.0;
+const Double_t AliITSRecoParam::fgkdiw=0.0053;
+const Double_t AliITSRecoParam::fgkX0iw=30.0;
+const Double_t AliITSRecoParam::fgkrcd=61.0;
+const Double_t AliITSRecoParam::fgkdcd=0.0053;
+const Double_t AliITSRecoParam::fgkX0cd=30.0;
+const Double_t AliITSRecoParam::fgkyr=12.8;
+const Double_t AliITSRecoParam::fgkdr=0.03;
+const Double_t AliITSRecoParam::fgkzm=0.2;
+const Double_t AliITSRecoParam::fgkdm=0.40;
+const Double_t AliITSRecoParam::fgkrs=50.0;
+const Double_t AliITSRecoParam::fgkds=0.001;
+const Double_t AliITSRecoParam::fgkrInsideITSscreen=49.0;
+const Double_t AliITSRecoParam::fgkrInsideSPD1=3.7;
+const Double_t AliITSRecoParam::fgkrPipe=3.;
+const Double_t AliITSRecoParam::fgkrInsidePipe=2.7;
+const Double_t AliITSRecoParam::fgkrOutsidePipe=3.3;
+const Double_t AliITSRecoParam::fgkdPipe=0.0028;
+const Double_t AliITSRecoParam::fgkrInsideShield[2]={7.5,25.0};
+const Double_t AliITSRecoParam::fgkrOutsideShield[2]={10.5,30.0};
+const Double_t AliITSRecoParam::fgkdshield[2]={0.0097,0.0034};
+const Double_t AliITSRecoParam::fgkX0shield[2]={38.6,42.0};
+const Double_t AliITSRecoParam::fgkX0Air=21.82;
+const Double_t AliITSRecoParam::fgkX0Be=65.19;
+const Double_t AliITSRecoParam::fgkBoundaryWidth=0.2;
+const Double_t AliITSRecoParam::fgkDeltaXNeighbDets=0.5;
+const Double_t AliITSRecoParam::fgkSPDdetzlength=6.960; // 7.072-2*0.056
+const Double_t AliITSRecoParam::fgkSPDdetxlength=1.298; // 1.410-2*0.056
 
 //_____________________________________________________________________________
-AliITSRecoParam::AliITSRecoParam()
+AliITSRecoParam::AliITSRecoParam() : AliDetectorRecoParam(),
+fTracker(0),
+fITSonly(kFALSE),
+fVertexer(0),
+fClusterFinder(0),
+fPID(0),
+fVtxr3DZCutWide(0.),
+fVtxr3DRCutWide(0.),
+fVtxr3DZCutNarrow(0.),
+fVtxr3DRCutNarrow(0.),
+fVtxr3DPhiCutLoose(0.),
+fVtxr3DPhiCutTight(0.),
+fVtxr3DDCACut(0.),
+fMaxSnp(1.),
+fNSigmaYLayerForRoadY(0),
+fNSigmaRoadY(0),
+fNSigmaZLayerForRoadZ(0),
+fNSigmaRoadZ(0),
+fNSigma2RoadZC(0),
+fNSigma2RoadYC(0),
+fNSigma2RoadZNonC(0),
+fNSigma2RoadYNonC(0),
+fRoadMisal(0),
+fMaxNormChi2NonCForHypothesis(0),
+fMaxChi2(0),
+fMaxRoad(0),
+fMaxChi2In(0),
+fChi2PerCluster(0),
+fXV(0), 
+fYV(0),
+fZV(0),
+fSigmaXV(0),
+fSigmaYV(0),
+fSigmaZV(0),
+fVertexCut(0),
+fMaxDZforPrimTrk(0),
+fMaxDZToUseConstraint(0), 
+fMaxDforV0dghtrForProlongation(0),
+fMaxDForProlongation(0),
+fMaxDZForProlongation(0),
+fMinPtForProlongation(0),
+fAddVirtualClustersInDeadZone(kFALSE),
+fZWindowDeadZone(0),
+fSigmaXDeadZoneHit2(0),
+fSigmaZDeadZoneHit2(0),
+fXPassDeadZoneHits(0),
+fUseTGeoInTracker(3),
+fStepSizeTGeo(0.1),
+fAllowSharedClusters(kTRUE),
+fClusterErrorsParam(1),
+fComputePlaneEff(kFALSE),
+fHistoPlaneEff(kFALSE),
+fIPlanePlaneEff(0),
+fReadPlaneEffFromOCDB(kFALSE),
+fMinPtPlaneEff(0),
+fMaxMissingClustersPlaneEff(0),
+fRequireClusterInOuterLayerPlaneEff(kFALSE),
+fRequireClusterInInnerLayerPlaneEff(kFALSE),
+fOnlyConstraintPlaneEff(kFALSE),
+fExtendedEtaAcceptance(kFALSE),
+fUseBadZonesFromOCDB(kTRUE),
+fUseSingleBadChannelsFromOCDB(kFALSE),
+fMinFractionOfBadInRoad(0),
+fAllowProlongationWithEmptyRoad(kFALSE),
+fOuterStartLayerSA(0),
+fFactorSAWindowSizes(1.),
+fNLoopsSA(33),
+fMinPhiSA(0.002),
+fMaxPhiSA(0.0145),
+fMinLambdaSA(0.003),
+fMaxLambdaSA(0.008),
+fSAOnePointTracks(kFALSE),
+fSAUseAllClusters(kFALSE),
+fFindV0s(kTRUE),
+fUseUnfoldingInClusterFinderSPD(kFALSE),
+fUseUnfoldingInClusterFinderSDD(kTRUE),
+fUseUnfoldingInClusterFinderSSD(kTRUE),
+fUseChargeMatchingInClusterFinderSSD(kTRUE),
+fTrackleterOnlyOneTrackletPerC2(kTRUE),
+fTrackleterPhiWindow(0.08),
+fTrackleterZetaWindow(1.00),
+fTrackleterRemoveClustersFromOverlaps(kFALSE),
+fTrackleterPhiOverlapCut(0.005),
+fTrackleterZetaOverlapCut(0.05),
+fUseCosmicRunShiftsSSD(kFALSE)
 {
   //
   // constructor
   //
+  SetName("ITS");
+  SetTitle("ITS");
+
   SetLayersParameters();
-  SetUseTGeoInTracker(kFALSE);
+  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) fLayersToSkip[i]=0;
+  SetUseTGeoInTracker(3);
+  SetStepSizeTGeo(0.1);
   SetAllowSharedClusters(kTRUE);
   SetFindV0s(kTRUE);
-  SetAddVirtualClustersInDeadZone(kTRUE);
+  SetAddVirtualClustersInDeadZone(kFALSE);
+  SetUseAmplitudeInfo(kTRUE);
+  SetClusterErrorsParam(1);
+  SetClusterMisalError(0.);
+  SetVertexer3DDefaults();
 }
 //_____________________________________________________________________________
 AliITSRecoParam::~AliITSRecoParam() 
@@ -53,6 +179,16 @@ AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
   // make default reconstruction  parameters for hig  flux env.
   //
   AliITSRecoParam *param = new AliITSRecoParam();
+  param->SetVertexerZ();
+
+  // use of bads from OCDB
+  param->SetUseBadZonesFromOCDB(kTRUE);
+  param->SetUseSingleBadChannelsFromOCDB(kFALSE);
+  // extended eta acceptance
+  param->SetExtendedEtaAcceptance(kFALSE);
+  // allow to skip layer if no cluster and no bad
+  param->SetAllowProlongationWithEmptyRoad(kFALSE);
+
 
   param->fMaxSnp = 0.95;
 
@@ -125,15 +261,14 @@ AliITSRecoParam *AliITSRecoParam::GetHighFluxParam()
   param->fMaxDZToUseConstraint = 3.;
 
   param->fMaxDforV0dghtrForProlongation = 30.;
-  param->fMaxDForProlongation = 10.;
-  param->fMaxDZForProlongation = 20.;
+  param->fMaxDForProlongation = 40.;//10.;
+  param->fMaxDZForProlongation = 60.;//20.;
   param->fMinPtForProlongation = 0.120;
 
   param->fZWindowDeadZone = 2.0;
   param->fSigmaXDeadZoneHit2 = 0.004/12.;
   param->fSigmaZDeadZoneHit2 = 0.001/12.;
   param->fXPassDeadZoneHits = 0.018;
-
   
   return param;
 }
@@ -143,7 +278,100 @@ AliITSRecoParam *AliITSRecoParam::GetLowFluxParam()
   //
   // make default reconstruction  parameters for low  flux env.
   //
-  return GetHighFluxParam();
+  AliITSRecoParam *param = new AliITSRecoParam();
+  param->SetVertexer3DDefaults();
+
+  // full use of bads from OCDB
+  param->SetUseBadZonesFromOCDB(kTRUE);
+  param->SetUseSingleBadChannelsFromOCDB(kTRUE);
+  // extended eta acceptance
+  param->SetExtendedEtaAcceptance(kTRUE);
+  // allow to skip layer if no cluster and no bad
+  param->SetAllowProlongationWithEmptyRoad(kTRUE);
+
+
+  param->fMaxSnp = 0.95;
+
+  param->fNSigmaYLayerForRoadY = 4.;
+  param->fNSigmaRoadY = 7.5;
+  param->fNSigmaZLayerForRoadZ = 4.;
+  param->fNSigmaRoadZ = 7.5;
+
+  param->fNSigma2RoadZC = 60.; //7.75^2
+  param->fNSigma2RoadYC = 60.; //7.75^2
+  param->fNSigma2RoadZNonC = 50.; //7.07^2
+  param->fNSigma2RoadYNonC = 50.; //7.07^2
+
+  param->fMaxChi2PerCluster[0] = 11.; //7
+  param->fMaxChi2PerCluster[1] = 12.; //5
+  param->fMaxChi2PerCluster[2] = 12.; //8
+  param->fMaxChi2PerCluster[3] = 5.;  //8
+  param->fMaxChi2PerCluster[4] = 12.; //6.5
+
+  param->fMaxNormChi2NonC[0] = 7.;
+  param->fMaxNormChi2NonC[1] = 8.;
+  param->fMaxNormChi2NonC[2] = 8.;
+  param->fMaxNormChi2NonC[3] = 11.;
+  param->fMaxNormChi2NonC[4] = 14.;
+  param->fMaxNormChi2NonC[5] = 25.;
+
+  param->fMaxNormChi2C[0] = 11.;
+  param->fMaxNormChi2C[1] = 13.;
+  param->fMaxNormChi2C[2] = 15.;
+  param->fMaxNormChi2C[3] = 18.;
+  param->fMaxNormChi2C[4] = 30.;
+  param->fMaxNormChi2C[5] = 35.;
+
+  param->fMaxNormChi2NonCForHypothesis = 7.;
+  
+  param->fMaxChi2 = 35.;
+
+  param->fMaxChi2s[0] = 25.; //40   
+  param->fMaxChi2s[1] = 25.; //40   
+  param->fMaxChi2s[2] = 25.; //40   
+  param->fMaxChi2s[3] = 25.; //40   
+  param->fMaxChi2s[4] = 40.; //40   
+  param->fMaxChi2s[5] = 50.; //40
+
+  param->fMaxRoad = 6.;
+
+  // not used
+  param->fMaxChi2In = 16.;
+   
+  param->fMaxChi2sR[0] = 10.;   
+  param->fMaxChi2sR[1] = 10.;   
+  param->fMaxChi2sR[2] = 10.;   
+  param->fMaxChi2sR[3] = 10.;   
+  param->fMaxChi2sR[4] = 30.;   
+  param->fMaxChi2sR[5] = 40.;   
+
+  param->fChi2PerCluster = 9.;
+  // not used
+
+  param->fXV = 0.;
+  param->fYV = 0.;
+  param->fZV = 0.;
+  param->fSigmaXV = 0.0050;
+  param->fSigmaYV = 0.0050;
+  param->fSigmaZV = 0.0100;
+
+  param->fVertexCut = 25.;
+
+  param->fMaxDZforPrimTrk = 0.4;
+  param->fMaxDZToUseConstraint = 3.;
+
+  param->fMaxDforV0dghtrForProlongation = 30.;
+  param->fMaxDForProlongation = 40.;//10.;
+  param->fMaxDZForProlongation = 60.;//20.;
+  param->fMinPtForProlongation = 0.120;
+
+  param->fZWindowDeadZone = 2.0;
+  param->fSigmaXDeadZoneHit2 = 0.004/12.;
+  param->fSigmaZDeadZoneHit2 = 0.001/12.;
+  param->fXPassDeadZoneHits = 0.018;
+  param->SetNLoopsSA(10);
+
+  return param;
 }
 //_____________________________________________________________________________
 AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam() 
@@ -151,7 +379,141 @@ AliITSRecoParam *AliITSRecoParam::GetCosmicTestParam()
   //
   // make default reconstruction  parameters for cosmics
   //
-  return GetHighFluxParam();
+  AliITSRecoParam *param = new AliITSRecoParam();
+
+  // vertexer for cosmics
+  param->SetVertexer(2);
+
+  // find independently ITS SA tracks 
+  param->SetSAUseAllClusters();
+  param->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);
+
+  // larger seach windows for SA (in case of large misalignments)
+  param->SetFactorSAWindowSizes(3.); 
+
+  param->fMaxSnp = 0.95;
+
+  param->fNSigmaYLayerForRoadY = 4.;
+  param->fNSigmaRoadY = 7.5;
+  param->fNSigmaZLayerForRoadZ = 4.;
+  param->fNSigmaRoadZ = 7.5;
+
+  param->fNSigma2RoadZC = 60.; //7.75^2
+  param->fNSigma2RoadYC = 60.; //7.75^2
+  param->fNSigma2RoadZNonC = 50.; //7.07^2
+  param->fNSigma2RoadYNonC = 50.; //7.07^2
+
+  param->fMaxChi2PerCluster[0] = 11.; //7
+  param->fMaxChi2PerCluster[1] = 12.; //5
+  param->fMaxChi2PerCluster[2] = 12.; //8
+  param->fMaxChi2PerCluster[3] = 5.;  //8
+  param->fMaxChi2PerCluster[4] = 12.; //6.5
+
+  param->fMaxNormChi2NonC[0] = 7.;
+  param->fMaxNormChi2NonC[1] = 8.;
+  param->fMaxNormChi2NonC[2] = 8.;
+  param->fMaxNormChi2NonC[3] = 11.;
+  param->fMaxNormChi2NonC[4] = 14.;
+  param->fMaxNormChi2NonC[5] = 25.;
+
+  param->fMaxNormChi2C[0] = 11.;
+  param->fMaxNormChi2C[1] = 13.;
+  param->fMaxNormChi2C[2] = 15.;
+  param->fMaxNormChi2C[3] = 18.;
+  param->fMaxNormChi2C[4] = 30.;
+  param->fMaxNormChi2C[5] = 35.;
+
+  param->fMaxNormChi2NonCForHypothesis = 7.;
+  
+  param->fMaxChi2 = 35.;
+
+  param->fMaxChi2s[0] = 25.; //40   
+  param->fMaxChi2s[1] = 25.; //40   
+  param->fMaxChi2s[2] = 25.; //40   
+  param->fMaxChi2s[3] = 25.; //40   
+  param->fMaxChi2s[4] = 40.; //40   
+  param->fMaxChi2s[5] = 50.; //40
+
+  param->fMaxRoad = 6.;
+
+  // not used
+  param->fMaxChi2In = 16.;
+   
+  param->fMaxChi2sR[0] = 10.;   
+  param->fMaxChi2sR[1] = 10.;   
+  param->fMaxChi2sR[2] = 10.;   
+  param->fMaxChi2sR[3] = 10.;   
+  param->fMaxChi2sR[4] = 30.;   
+  param->fMaxChi2sR[5] = 40.;   
+
+  param->fChi2PerCluster = 9.;
+  // not used
+
+  param->fXV = 0.;
+  param->fYV = 0.;
+  param->fZV = 0.;
+  param->fSigmaXV = 0.0050;
+  param->fSigmaYV = 0.0050;
+  param->fSigmaZV = 0.0100;
+
+  param->fVertexCut = 25.;
+
+  param->fMaxDZforPrimTrk = 0.4;
+  param->fMaxDZToUseConstraint = 3.;
+
+  param->fMaxDforV0dghtrForProlongation = 30.;
+  param->fMaxDForProlongation = 10.;
+  param->fMaxDZForProlongation = 20.;
+  param->fMinPtForProlongation = 0.120;
+
+  param->fZWindowDeadZone = 2.0;
+  param->fSigmaXDeadZoneHit2 = 0.004/12.;
+  param->fSigmaZDeadZoneHit2 = 0.001/12.;
+  param->fXPassDeadZoneHits = 0.018;
+
+  
+  return param;
+}
+//_____________________________________________________________________________
+AliITSRecoParam *AliITSRecoParam::GetPlaneEffParam(Int_t i)
+{
+  //
+  // make special reconstruction parameters for Plane Efficiency study on layer i
+  //
+  if (i<0 || i>=AliITSgeomTGeo::kNLayers) {
+    printf("AliITSRecoParam::GetPlaneEffParam: index of ITS Plane not in the range [0,5]\n");
+    printf("returning null pointer");
+    return NULL;
+  }
+  AliITSRecoParam *param;
+  param = GetHighFluxParam();
+  param->SetComputePlaneEff();
+  param->SetLayerToSkip(i);
+  param->SetIPlanePlaneEff(i);
+  // optimized setting for SPD0 (i==0)
+  if (i==0 || i==1) {
+    param->fMinPtPlaneEff = 0.200; // high pt particles
+    param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
+    param->fRequireClusterInOuterLayerPlaneEff = kTRUE; // cluster on SPD1 must be
+    //param->fOnlyConstraintPlaneEff = kTRUE;
+  }
+  if (i==2 || i==3) {
+    param->fMinPtPlaneEff = 0.200; // high pt particles
+    param->fMaxMissingClustersPlaneEff = 1; // at most 1 layer out of 5 without cluster
+    param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
+    //param->fOnlyConstraintPlaneEff = kTRUE;
+  }
+  if (i==4) {
+    param->fMinPtPlaneEff = 0.200; // high pt particles
+    param->fMaxMissingClustersPlaneEff = 0; // at most 1 layer out of 5 without cluster
+    param->fRequireClusterInOuterLayerPlaneEff = kTRUE;
+    //param->fOnlyConstraintPlaneEff = kTRUE;
+  }
+  if (i==5) {
+    param->fMinPtPlaneEff = 0.200; // high pt particles
+  }
+  //
+  return param;
 }
 //_____________________________________________________________________________
 void AliITSRecoParam::SetLayersParameters() 
@@ -178,4 +540,34 @@ void AliITSRecoParam::SetLayersParameters()
 
   return;
 }
+//_____________________________________________________________________________
+void AliITSRecoParam::PrintParameters() const 
+{
+  //
+  // print parameters
+  //
+
+  printf("=============================  AliITSRecoParam::PrintParameters ");
+  printf("============================= \n\n");
+  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
+    if(!fLayersToSkip[i]) {
+      printf("ITS Traking: using layer %d\n",i);
+    } else {
+      printf("ITS Traking: skipping layer %d\n",i);
+    }
+  }
+  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++) {
+    if(fUseAmplitudeInfo[i]) {
+      printf("ITS Traking: use amplitude info for layer %d\n",i);
+    } else {
+      printf("ITS Traking: don't use amplitude info for layer %d\n",i);
+    }
+  }
+  for(Int_t i=0; i<AliITSgeomTGeo::kNLayers; i++)
+    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]);
+
 
+  Dump();
+
+  return;
+}