}
+//_____________________________________________________________________________
+Double_t SigmaY2(Double_t r, Double_t tgl, Double_t pt)
+{
+ //
+ // Parametrised error of the cluster reconstruction (pad direction)
+ //
+ // Sigma rphi
+ const Float_t kArphi=0.41818e-2;
+ const Float_t kBrphi=0.17460e-4;
+ const Float_t kCrphi=0.30993e-2;
+ const Float_t kDrphi=0.41061e-3;
+
+ pt=TMath::Abs(pt)*1000.;
+ Double_t x=r/pt;
+ tgl=TMath::Abs(tgl);
+ Double_t s=kArphi - kBrphi*r*tgl + kCrphi*x*x + kDrphi*x;
+ if (s<0.4e-3) s=0.4e-3;
+ s*=1.3; //Iouri Belikov
+
+ return s;
+}
+
+//_____________________________________________________________________________
+Double_t SigmaZ2(Double_t r, Double_t tgl)
+{
+ //
+ // Parametrised error of the cluster reconstruction (drift direction)
+ //
+ // Sigma z
+ const Float_t kAz=0.39614e-2;
+ const Float_t kBz=0.22443e-4;
+ const Float_t kCz=0.51504e-1;
+
+
+ tgl=TMath::Abs(tgl);
+ Double_t s=kAz - kBz*r*tgl + kCz*tgl*tgl;
+ if (s<0.4e-3) s=0.4e-3;
+ s*=1.3; //Iouri Belikov
+
+ return s;
+}
/*
$Log$
+Revision 1.38 2003/10/17 12:01:16 kowal2
+Removed compiler warning.
+
Revision 1.37 2003/07/22 15:56:14 hristov
Implementing ESD functionality in the NewIO (Yu.Belikov)
#include "AliTPCParam.h"
#include "AliClusters.h"
+extern Double_t SigmaY2(Double_t, Double_t, Double_t);
+extern Double_t SigmaZ2(Double_t, Double_t);
+
ClassImp(AliTPCtracker)
//_____________________________________________________________________________
fBarrelTrack->SetRefPlane(refPlane, isIn);
}
-//_____________________________________________________________________________
-Double_t SigmaY2(Double_t r, Double_t tgl, Double_t pt)
-{
- //
- // Parametrised error of the cluster reconstruction (pad direction)
- //
- // Sigma rphi
- const Float_t kArphi=0.41818e-2;
- const Float_t kBrphi=0.17460e-4;
- const Float_t kCrphi=0.30993e-2;
- const Float_t kDrphi=0.41061e-3;
-
- pt=TMath::Abs(pt)*1000.;
- Double_t x=r/pt;
- tgl=TMath::Abs(tgl);
- Double_t s=kArphi - kBrphi*r*tgl + kCrphi*x*x + kDrphi*x;
- if (s<0.4e-3) s=0.4e-3;
- s*=1.3; //Iouri Belikov
-
- return s;
-}
-
-//_____________________________________________________________________________
-Double_t SigmaZ2(Double_t r, Double_t tgl)
-{
- //
- // Parametrised error of the cluster reconstruction (drift direction)
- //
- // Sigma z
- const Float_t kAz=0.39614e-2;
- const Float_t kBz=0.22443e-4;
- const Float_t kCz=0.51504e-1;
-
-
- tgl=TMath::Abs(tgl);
- Double_t s=kAz - kBz*r*tgl + kCz*tgl*tgl;
- if (s<0.4e-3) s=0.4e-3;
- s*=1.3; //Iouri Belikov
-
- return s;
-}
//_____________________________________________________________________________
Double_t f1(Double_t x1,Double_t y1,