store ntuple in same file as histos (Renu)

[u/mrichter/AliRoot.git] / TPC / AliTPCcluster.cxx

diff --git a/TPC/AliTPCcluster.cxx b/TPC/AliTPCcluster.cxx
index 107417824f5d4978a0f6cff592f4fdd7bf468b36..22232c020a68354f42d7129d927b04c86ef5cf0b 100644 (file)
--- a/TPC/AliTPCcluster.cxx
+++ b/TPC/AliTPCcluster.cxx
@@ -24,4 +24,57 @@
 #include "AliTPCcluster.h"
 
 ClassImp(AliTPCcluster)
+//
+AliTPCcluster::AliTPCcluster(Int_t *lab, Float_t *hit):
+  AliCluster(0,hit,0.,0.,lab),
+  fQ(0.)
+{
+  //
+  // constructor
+  //
+  fQ=hit[4];
+}
+ 
+//____________________________________________________________________________
+Double_t AliTPCcluster::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 AliTPCcluster::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;
+}