added verbosity to QA histograms (Yves)

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

diff --git a/TRD/AliTRDrecoParam.cxx b/TRD/AliTRDrecoParam.cxx
index 4ce768dcb683b516809af748292b28c5e85cff7f..63b61420f18c325f3aa360627c7a292003a4180a 100644 (file)
--- a/TRD/AliTRDrecoParam.cxx
+++ b/TRD/AliTRDrecoParam.cxx
@@ -26,16 +26,16 @@
 ///////////////////////////////////////////////////////////////////////////////
 
 #include "AliTRDrecoParam.h"
+#include "AliTRDrawStreamBase.h"
 
 ClassImp(AliTRDrecoParam)
 
 //______________________________________________________________
 AliTRDrecoParam::AliTRDrecoParam()
   :AliDetectorRecoParam()
-  ,fkClusterSharing(0)
-  ,fkPIDMethod(1) // LQ PID
-  ,fkMaxTheta(1.0)
-  ,fkMaxPhi(2.0)
+  ,fkdNchdy(12.)
+  ,fkMaxTheta(1.0)     
+  ,fkMaxPhi(2.0) 
   ,fkRoad0y(6.0)
   ,fkRoad0z(8.5) 
   ,fkRoad1y(2.0)
@@ -46,12 +46,73 @@ AliTRDrecoParam::AliTRDrecoParam()
   ,fkFindable(.333)
   ,fkChi2Z(30./*14.*//*12.5*/)
   ,fkChi2Y(.25)
+  ,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(1)
+  ,fNumberOfPresamples(0)
+  ,fNumberOfPostsamples(0)
 {
   //
   // Default constructor
   //
+  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
 
+  memset(fPIDThreshold, 0, AliTRDCalPID::kNMom*sizeof(Double_t));
+}
+
+//______________________________________________________________
+AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
+  :AliDetectorRecoParam(ref)
+  ,fkdNchdy(ref.fkdNchdy)
+  ,fkMaxTheta(ref.fkMaxTheta)
+  ,fkMaxPhi(ref.fkMaxPhi)
+  ,fkRoad0y(ref.fkRoad0y)
+  ,fkRoad0z(ref.fkRoad0z) 
+  ,fkRoad1y(ref.fkRoad1y)
+  ,fkRoad1z(ref.fkRoad1z)      
+  ,fkRoad2y(ref.fkRoad2y)
+  ,fkRoad2z(ref.fkRoad2z)
+  ,fkPlaneQualityThreshold(ref.fkPlaneQualityThreshold)
+  ,fkFindable(ref.fkFindable)
+  ,fkChi2Z(ref.fkChi2Z)
+  ,fkChi2Y(ref.fkChi2Y)
+  ,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)
+  ,fClusMaxThresh(ref.fClusMaxThresh)
+  ,fClusSigThresh(ref.fClusSigThresh)
+  ,fTCnexp(ref.fTCnexp)
+  ,fNumberOfPresamples(ref.fNumberOfPresamples)
+  ,fNumberOfPostsamples(ref.fNumberOfPostsamples)
+{
+  //
+  // Copy constructor
+  //
+  memcpy(fSysCovMatrix, ref.fSysCovMatrix, 5*sizeof(Double_t));
+  memcpy(fPIDThreshold, ref.fPIDThreshold, AliTRDCalPID::kNMom*sizeof(Double_t));
 }
 
 //______________________________________________________________
@@ -61,7 +122,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;
 
 }
 
@@ -72,6 +135,43 @@ AliTRDrecoParam *AliTRDrecoParam::GetHighFluxParam()
   // Parameters for the high flux environment
   //
 
-  return new AliTRDrecoParam();
+  AliTRDrecoParam *rec = new AliTRDrecoParam();
+  rec->fkdNchdy = 4000.; // PbPb in TRD
+
+  return rec;
+
+}
+
+//______________________________________________________________
+AliTRDrecoParam *AliTRDrecoParam::GetCosmicTestParam()
+{
+  //
+  // Parameters for the cosmics data
+  //
+
+  AliTRDrawStreamBase::SetRawStreamVersion("TB");
+  AliTRDrecoParam *par = new AliTRDrecoParam();
+  par->fSysCovMatrix[0] = 1.; // y direction (1 cm)
+  par->fSysCovMatrix[1] = 1.; // z direction (1 cm)
+  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;
+}