o fix integration

diff --git a/TPC/Base/AliTPCCorrection.cxx b/TPC/Base/AliTPCCorrection.cxx
index 511800e7b051abdea2d91021b2b56f33e937b0c2..2ba1b2a3ca48714d82c30775a204e5b3c7a76db1 100644 (file)
--- a/TPC/Base/AliTPCCorrection.cxx
+++ b/TPC/Base/AliTPCCorrection.cxx
@@ -455,7 +455,8 @@ void AliTPCCorrection::GetCorrectionIntegralDz(const Float_t x[],const Short_t r
     //  the slopes will be positive.
     // but since we chose deltaZ opposite sign the singn of the corretion should be fine
     
-    GetCorrectionDz(xyz,roc,dxyz,delta);
+    Float_t xyz2[3]={xyz[0],xyz[1],xyz[2]+deltaZ/2.};
+    GetCorrectionDz(xyz2,roc,dxyz,delta/2.);
     xyz[0]+=deltaZ*dxyz[0];
     xyz[1]+=deltaZ*dxyz[1];
     xyz[2]+=deltaZ;           //
@@ -497,7 +498,8 @@ void AliTPCCorrection::GetDistortionIntegralDz(const Float_t x[],const Short_t r
     // and since we are moving towards the read-out plane the deltaZ for
     //   weighting the dK/dz should have the opposite sign
     deltaZ*=sign;
-    GetDistortionDz(xyz,roc,dxyz,delta);
+    Float_t xyz2[3]={xyz[0],xyz[1],xyz[2]+deltaZ/2.};
+    GetDistortionDz(xyz2,roc,dxyz,delta/2.);
     xyz[0]+=-deltaZ*dxyz[0];
     xyz[1]+=-deltaZ*dxyz[1];
     xyz[2]+=deltaZ;           //TODO: Should this also be corrected for the dxyz[2]

diff --git a/TPC/Upgrade/AliToyMCReconstruction.cxx b/TPC/Upgrade/AliToyMCReconstruction.cxx
index 505e3954174649156165691f0e77c6667eef263b..b1f01dbcf5451775cceafe99d99a0105397cd2ff 100644 (file)
--- a/TPC/Upgrade/AliToyMCReconstruction.cxx
+++ b/TPC/Upgrade/AliToyMCReconstruction.cxx
@@ -312,7 +312,7 @@ void AliToyMCReconstruction::RunReco(const char* file, Int_t nmaxEv)
         if (arrClustRes) {
           const Int_t nCl=arrClustRes->GetEntriesFast();
           // fracktion of outliers from track extrapolation
-          // for 1, 1.5, 2, 2.5 and 3 sigma of the cluster resolution (~1mm)
+          // for 3, 3.5, 4, 4.5 and 5 sigma of the cluster resolution (~1mm)
           Float_t fracY[5]={0.};
           Float_t fracZ[5]={0.};
           
@@ -321,8 +321,8 @@ void AliToyMCReconstruction::RunReco(const char* file, Int_t nmaxEv)
             const Float_t sigmaY=TMath::Sqrt(cl->GetSigmaY2());
             const Float_t sigmaZ=TMath::Sqrt(cl->GetSigmaZ2());
             for (Int_t inSig=0; inSig<5; ++inSig) {
-              fracY[inSig] += cl->GetY()>(1+inSig*.5)*sigmaY;
-              fracZ[inSig] += cl->GetZ()>(1+inSig*.5)*sigmaZ;
+              fracY[inSig] += cl->GetY()>(3+inSig*.5)*sigmaY;
+              fracZ[inSig] += cl->GetZ()>(3+inSig*.5)*sigmaZ;
             }
           }
           
@@ -336,8 +336,8 @@ void AliToyMCReconstruction::RunReco(const char* file, Int_t nmaxEv)
           (*fStreamer) << "Tracks" <<
           "clustRes.=" << arrClustRes;
           for (Int_t inSig=0; inSig<5; ++inSig) {
-            const char* fracYname=Form("clFracY%02d=", 10+inSig*5);
-            const char* fracZname=Form("clFracZ%02d=", 10+inSig*5);
+            const char* fracYname=Form("clFracY%02d=", 30+inSig*5);
+            const char* fracZname=Form("clFracZ%02d=", 30+inSig*5);
             (*fStreamer) << "Tracks" <<
             fracYname << fracY[inSig] <<
             fracZname << fracZ[inSig];