Protection added.

[u/mrichter/AliRoot.git] / TRD / AliTRDrecoParam.cxx

diff --git a/TRD/AliTRDrecoParam.cxx b/TRD/AliTRDrecoParam.cxx
index 67f5bedff9c184838ec69b34bad74e7275ce12b9..41f87c82eabd70f65b79a9fbe54a068b63f6309e 100644 (file)
--- a/TRD/AliTRDrecoParam.cxx
+++ b/TRD/AliTRDrecoParam.cxx
@@ -33,6 +33,7 @@ ClassImp(AliTRDrecoParam)
 //______________________________________________________________
 AliTRDrecoParam::AliTRDrecoParam()
   :AliDetectorRecoParam()
+  ,fkdNchdy(12.)
   ,fkMaxTheta(1.0)     
   ,fkMaxPhi(2.0) 
   ,fkRoad0y(6.0)
@@ -45,17 +46,20 @@ AliTRDrecoParam::AliTRDrecoParam()
   ,fkFindable(.333)
   ,fkChi2Z(30./*14.*//*12.5*/)
   ,fkChi2Y(.25)
-  ,fkChi2YCut(7.73)
-  ,fkChi2ZCut(0.069)
-  ,fkPhiCut(10.6)
-  ,fkMeanNclusters(72)
-  ,fkSigmaNclusters(5.2632)
+  ,fkChi2YSlope(7.73)
+  ,fkChi2ZSlope(0.069)
+       ,fkChi2YCut(0.5)
+  ,fkPhiSlope(10.6)
+  ,fkNMeanClusters(20.)
+  ,fkNSigmaClusters(2.)
+  ,fkNClusterNoise(0.)
+  ,fkNMeanTracklets(5.5)
   ,fkTrackLikelihood(-15.)
   ,fMinMaxCutSigma(4.)
   ,fMinLeftRightCutSigma(8.)
   ,fClusMaxThresh(4.5)
   ,fClusSigThresh(3.5)
-  ,fTCnexp(2)
+  ,fTCnexp(1)
   ,fNumberOfPresamples(0)
   ,fNumberOfPostsamples(0)
 {
@@ -68,8 +72,8 @@ AliTRDrecoParam::AliTRDrecoParam()
   SetVertexConstrained();
   SetImproveTracklets(kFALSE);
 
-  fSysCovMatrix[0] = 1.; // y direction (1 cm)
-  fSysCovMatrix[1] = 1.; // z direction (1 cm)
+  fSysCovMatrix[0] = 0.; // y direction (1 cm)
+  fSysCovMatrix[1] = 0.; // z direction (1 cm)
   fSysCovMatrix[2] = 0.; // snp
   fSysCovMatrix[3] = 0.; // tgl
   fSysCovMatrix[4] = 0.; // 1/pt
@@ -80,6 +84,7 @@ AliTRDrecoParam::AliTRDrecoParam()
 //______________________________________________________________
 AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
   :AliDetectorRecoParam(ref)
+  ,fkdNchdy(ref.fkdNchdy)
   ,fkMaxTheta(ref.fkMaxTheta)
   ,fkMaxPhi(ref.fkMaxPhi)
   ,fkRoad0y(ref.fkRoad0y)
@@ -92,11 +97,14 @@ AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
   ,fkFindable(ref.fkFindable)
   ,fkChi2Z(ref.fkChi2Z)
   ,fkChi2Y(ref.fkChi2Y)
-  ,fkChi2YCut(ref.fkChi2YCut)
-  ,fkChi2ZCut(ref.fkChi2ZCut)
-  ,fkPhiCut(ref.fkPhiCut)
-  ,fkMeanNclusters(ref.fkMeanNclusters)
-  ,fkSigmaNclusters(ref.fkSigmaNclusters)
+  ,fkChi2YSlope(ref.fkChi2YSlope)
+  ,fkChi2ZSlope(ref.fkChi2ZSlope)
+       ,fkChi2YCut(ref.fkChi2YCut)
+  ,fkPhiSlope(ref.fkPhiSlope)
+  ,fkNMeanClusters(ref.fkNMeanClusters)
+  ,fkNSigmaClusters(ref.fkNSigmaClusters)
+  ,fkNClusterNoise(ref.fkNClusterNoise)
+  ,fkNMeanTracklets(ref.fkNMeanTracklets)
   ,fkTrackLikelihood(ref.fkTrackLikelihood)
   ,fMinMaxCutSigma(ref.fMinMaxCutSigma)
   ,fMinLeftRightCutSigma(ref.fMinLeftRightCutSigma)
@@ -126,7 +134,9 @@ AliTRDrecoParam *AliTRDrecoParam::GetLowFluxParam()
   // Parameters for the low flux environment
   //
 
-  return new AliTRDrecoParam();
+  AliTRDrecoParam *rec = new AliTRDrecoParam();
+  rec->fkdNchdy = 12.; // pp in TRD
+  return rec;
 
 }
 
@@ -138,7 +148,9 @@ AliTRDrecoParam *AliTRDrecoParam::GetHighFluxParam()
   //
 
   AliTRDrecoParam *rec = new AliTRDrecoParam();
+  rec->fkdNchdy = 4000.; // PbPb in TRD
   rec->SetImproveTracklets(kTRUE);
+
   return rec;
 
 }
@@ -152,13 +164,28 @@ AliTRDrecoParam *AliTRDrecoParam::GetCosmicTestParam()
 
   AliTRDrawStreamBase::SetRawStreamVersion("TB");
   AliTRDrecoParam *par = new AliTRDrecoParam();
+  par->fSysCovMatrix[0] = 1.; // y direction (1 cm)
+  par->fSysCovMatrix[1] = 1.; // z direction (1 cm)
   par->SetVertexConstrained(kFALSE);
-  par->SetChi2YCut(1.136);
-  par->SetChi2ZCut(0.069);
-  par->SetMaxTheta(2.1445);
-  par->SetMaxPhi(2.7475);
-  par->SetMeanNclusters(48.1197);
-  par->SetSigmaNclusters(8.59347);
+  par->fkChi2YSlope       = 0.11853;
+  par->fkChi2ZSlope       = 0.04527;
+       par->fkChi2YCut                                 = 1.;
+       par->fkPhiSlope                                 = 3.17954;
+  par->fkMaxTheta       = 2.1445;
+  par->fkMaxPhi         = 2.7475;
+  par->fkNMeanClusters  = 12.89;
+  par->fkNSigmaClusters = 2.095;
   return par;
 
 }
+
+
+//______________________________________________________________
+Float_t AliTRDrecoParam::GetNClusters() const
+{
+  // Estimate the number of clusters in the TRD detector
+  
+  Float_t nclusters = (fkNMeanClusters + 2*fkNSigmaClusters)*fkNMeanTracklets*fkdNchdy;
+  nclusters *= 1.+fkNClusterNoise;
+  return nclusters;
+}