From: cblume Date: Thu, 27 Aug 2009 19:02:52 +0000 (+0000) Subject: Changes for cosmics reconstruction X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=566bf88701f29a57ead95f17240d28f0d8c765d7;hp=d4cbb500156b7953228b25f39b9bac22d65ab679 Changes for cosmics reconstruction --- diff --git a/TRD/AliTRDrecoParam.cxx b/TRD/AliTRDrecoParam.cxx index 310f92a029a..c60fb5386a6 100644 --- a/TRD/AliTRDrecoParam.cxx +++ b/TRD/AliTRDrecoParam.cxx @@ -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; } diff --git a/TRD/AliTRDrecoParam.h b/TRD/AliTRDrecoParam.h index 1be90bbcaa8..0eba24f4560 100644 --- a/TRD/AliTRDrecoParam.h +++ b/TRD/AliTRDrecoParam.h @@ -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 diff --git a/TRD/AliTRDseedV1.cxx b/TRD/AliTRDseedV1.cxx index 91e50ca6917..123f8680da0 100644 --- a/TRD/AliTRDseedV1.cxx +++ b/TRD/AliTRDseedV1.cxx @@ -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); diff --git a/TRD/AliTRDtransform.cxx b/TRD/AliTRDtransform.cxx index 77ca47c1b5c..f5d9b18a6d5 100644 --- a/TRD/AliTRDtransform.cxx +++ b/TRD/AliTRDtransform.cxx @@ -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);