Changes for cosmics reconstruction
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 27 Aug 2009 19:02:52 +0000 (19:02 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 27 Aug 2009 19:02:52 +0000 (19:02 +0000)
TRD/AliTRDrecoParam.cxx
TRD/AliTRDrecoParam.h
TRD/AliTRDseedV1.cxx
TRD/AliTRDtransform.cxx

index 310f92a..c60fb53 100644 (file)
@@ -43,6 +43,7 @@ AliTRDrecoParam::AliTRDrecoParam()
   ,fkRoad2y(3.0)
   ,fkRoad2z(20.0)
   ,fkPlaneQualityThreshold(5.0)// 4.2? under Investigation
+  ,fkRoadzMultiplicator(1.5)
   ,fkFindable(.333)
   ,fkChi2Z(30./*14.*//*12.5*/)
   ,fkChi2Y(.25)
@@ -88,6 +89,7 @@ AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
   ,fkRoad2y(ref.fkRoad2y)
   ,fkRoad2z(ref.fkRoad2z)
   ,fkPlaneQualityThreshold(ref.fkPlaneQualityThreshold)
+  ,fkRoadzMultiplicator(ref.fkRoadzMultiplicator)
   ,fkFindable(ref.fkFindable)
   ,fkChi2Z(ref.fkChi2Z)
   ,fkChi2Y(ref.fkChi2Y)
@@ -161,6 +163,7 @@ AliTRDrecoParam *AliTRDrecoParam::GetCosmicTestParam()
   par->fkMaxPhi         = 2.7475;
   par->fkNMeanClusters  = 12.89;
   par->fkNSigmaClusters = 2.095;
+  par->fkRoadzMultiplicator = 3.;
   return par;
 
 }
index 1be90bb..0eba24f 100644 (file)
@@ -45,6 +45,7 @@ public:
   Double_t GetRoad1z() const                { return fkRoad1z;   }
   Double_t GetRoad2y() const                { return fkRoad2y;   }
   Double_t GetRoad2z() const                { return fkRoad2z;   }
+  Double_t GetRoadzMultiplicator() const    { return fkRoadzMultiplicator; }
   Double_t GetTrackLikelihood() const       { return fkTrackLikelihood;       }
   inline void GetSysCovMatrix(Double_t *sys) const;  
   Double_t GetMinMaxCutSigma() const        { return fMinMaxCutSigma;     };
@@ -97,6 +98,7 @@ private:
   Double_t  fkRoad2z;                // Road in z for extrapolated cluster
   
   Double_t  fkPlaneQualityThreshold; // Quality threshold
+  Double_t  fkRoadzMultiplicator;    // Multiplicator for the Roads in z 
   Double_t  fkFindable;              // minimum ratio of clusters per tracklet supposed to be attached.
   Double_t  fkChi2Z;                 // Max chi2 on the z direction for seeding clusters fit
   Double_t  fkChi2Y;                 // Max chi2 on the y direction for seeding clusters Rieman fit
index 91e50ca..123f868 100644 (file)
@@ -906,7 +906,7 @@ Bool_t      AliTRDseedV1::AttachClusters(AliTRDtrackingChamber *const chamber, Bool_t
            t2    = GetTilt()*GetTilt();
   //define roads
   Double_t kroady = 1., //fkReconstructor->GetRecoParam() ->GetRoad1y();
-           kroadz = GetPadLength() * 1.5 + 1.;
+           kroadz = GetPadLength() * fkReconstructor->GetRecoParam()->GetRoadzMultiplicator() + 1.;
   // define probing cluster (the perfect cluster) and default calibration
   Short_t sig[] = {0, 0, 10, 30, 10, 0,0};
   AliTRDcluster cp(fDet, 6, 75, 0, sig, 0);
index 77ca47c..f5d9b18 100644 (file)
@@ -221,6 +221,7 @@ Bool_t AliTRDtransform::Transform(AliTRDcluster *c)
   Double_t vd  = fCalVdriftDetValue * fCalVdriftROC->GetValue(col,row);
   // t0
   Double_t t0  = fCalT0DetValue     + fCalT0ROC->GetValue(col,row);
+  t0 /= fSamplingFrequency;
   // ExB correction
   Double_t exb = AliTRDCommonParam::Instance()->GetOmegaTau(vd);