From: abercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 29 Sep 2010 10:40:48 +0000 (+0000)
Subject: modify systematic r-phi plotting
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=ea2e835e9f3c16e527ee28ea045e942e02c085bd;p=u%2Fmrichter%2FAliRoot.git

modify systematic r-phi plotting
---

diff --git a/PWG1/TRD/AliTRDclusterResolution.cxx b/PWG1/TRD/AliTRDclusterResolution.cxx
index 95d97d4a154..cb7ac0830ab 100644
--- a/PWG1/TRD/AliTRDclusterResolution.cxx
+++ b/PWG1/TRD/AliTRDclusterResolution.cxx
@@ -517,7 +517,7 @@ void AliTRDclusterResolution::UserExec(Option_t *)
   TH3S *h3 = NULL;
 
   // define limits around ExB for which x contribution is negligible
-  const Float_t kDtgPhi = 3.5e-2; //(+- 2 deg)
+  const Float_t kAroundZero = 3.5e-2; //(+- 2 deg)
 
   TObjArray *arr0 = (TObjArray*)fContainer->At(kCenter);
   TObjArray *arr1 = (TObjArray*)fContainer->At(kSigm);
@@ -540,10 +540,11 @@ void AliTRDclusterResolution::UserExec(Option_t *)
     AliDebug(4, Form("det[%d] tb[%2d] q[%4.0f Log[%6.4f]] dy[%7.2f][um] ypull[%5.2f]", det, t, q, TMath::Log(q), 1.e4*dy, dy/TMath::Sqrt(covcl[0])));
     
     cli->GetGlobalPosition(y, z, dydx, dzdx, &cov[0]);
+    Float_t tilt(cli->GetTilt());
 
     // resolution as a function of cluster charge
     // only for phi equal exB 
-    if(TMath::Abs(dydx-fExB) < kDtgPhi){
+    if(TMath::Abs(dydx-fExB-tilt*dzdx) < kAroundZero){
       h3 = (TH3S*)fContainer->At(kQRes);
       h3->Fill(TMath::Log(q), dy, dy/TMath::Sqrt(covcl[0]));
     }
@@ -556,7 +557,7 @@ void AliTRDclusterResolution::UserExec(Option_t *)
 
     // resolution as a function of y displacement from pad center
     // only for phi equal exB
-    if(TMath::Abs(dydx-fExB) < kDtgPhi){
+    if(TMath::Abs(dydx-fExB-tilt*dzdx) < kAroundZero){
       h3 = (TH3S*)arr0->At(0);
       h3->Fill(t, cli->GetYDisplacement(), dy);
       h3 = (TH3S*)arr0->At(1);
@@ -566,7 +567,6 @@ void AliTRDclusterResolution::UserExec(Option_t *)
     Int_t it(((TH3S*)arr0->At(0))->GetXaxis()->FindBin(t));
 
     // fill histo for resolution (sigma)
-    Float_t tilt(cli->GetTilt());
     ((TH3S*)arr1->At(it-1))->Fill(tilt*dzdx, dydx, dy);
 
     // fill histo for systematic (mean)