Bugfix to use cylindrical coordinates consitently (in Marian creation of the residual...

diff --git a/TPC/AliTPCCorrection.cxx b/TPC/AliTPCCorrection.cxx
index 24e5ae3b72e3ceb9f5b7b34f1de95294e1396f0e..07b8b0fc67abc40c64ad6fc13363946425de9b3e 100644 (file)
--- a/TPC/AliTPCCorrection.cxx
+++ b/TPC/AliTPCCorrection.cxx
@@ -1927,7 +1927,7 @@ void AliTPCCorrection::AddVisualCorrection(AliTPCCorrection* corr, Int_t positio
 
 
 
-Double_t AliTPCCorrection::GetCorrSector(Double_t sector, Double_t localX, Double_t kZ, Int_t axisType, Int_t corrType){
+Double_t AliTPCCorrection::GetCorrSector(Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType){
   //
   // calculate the correction at given position - check the geffCorr
   //
@@ -1935,9 +1935,9 @@ Double_t AliTPCCorrection::GetCorrSector(Double_t sector, Double_t localX, Doubl
   AliTPCCorrection *corr = (AliTPCCorrection*)fgVisualCorrection->At(corrType);
   if (!corr) return 0;
   Double_t phi=sector*TMath::Pi()/9.;
-  Double_t gx = localX*TMath::Cos(phi);
-  Double_t gy = localX*TMath::Sin(phi);
-  Double_t gz = localX*kZ;
+  Double_t gx = r*TMath::Cos(phi);
+  Double_t gy = r*TMath::Sin(phi);
+  Double_t gz = r*kZ;
   Int_t nsector=(gz>0) ? 0:18; 
   //
   //

diff --git a/TPC/AliTPCCorrection.h b/TPC/AliTPCCorrection.h
index b1fd7f3f82ba8fa6057611cdee413ad539ef0dbc..571ac5843a5d78d7b8e71dc863661bfcd98b3e88 100644 (file)
--- a/TPC/AliTPCCorrection.h
+++ b/TPC/AliTPCCorrection.h
@@ -90,7 +90,7 @@ public:
   void FastSimDistortedVertex(Double_t orgVertex[3], Int_t nTracks, AliESDVertex &aV, AliESDVertex &avOrg, AliESDVertex &cV, AliESDVertex &cvOrg, TTreeSRedirector * const pcstream, Double_t etaCuts);
 
   static void AddVisualCorrection(AliTPCCorrection* corr, Int_t position);
-  static Double_t GetCorrSector(Double_t sector, Double_t localX, Double_t kZ, Int_t axisType, Int_t corrType=0);
+  static Double_t GetCorrSector(Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType=0);
   static Double_t GetCorrXYZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0);
 
 protected:

diff --git a/TPC/AliTPCcalibAlign.cxx b/TPC/AliTPCcalibAlign.cxx
index a2c30064f832899fafabeaef662ae0323453f1d8..27cd6fdaaec739502ca193dff580ad913901a965 100644 (file)
--- a/TPC/AliTPCcalibAlign.cxx
+++ b/TPC/AliTPCcalibAlign.cxx
@@ -1676,7 +1676,7 @@ void AliTPCcalibAlign::MakeResidualHistos(){
   axisName[1]="sector";   axisTitle[1]="Sector Number"; 
   binsTrack[1]=180;       xminTrack[1]=0;        xmaxTrack[1]=18; 
   //
-  axisName[2]="localX";   axisTitle[2]="x (cm)"; 
+  axisName[2]="R";   axisTitle[2]="r (cm)"; 
   binsTrack[2]=53;       xminTrack[2]=85.;        xmaxTrack[2]=245.; 
   //
   //
@@ -2722,7 +2722,7 @@ void AliTPCcalibAlign::UpdateClusterDeltaField(const AliTPCseed * seed){
     trackOut.GetXYZ(xyz);
     resVector[1]= 9.*TMath::ATan2(xyz[1],xyz[0])/TMath::Pi();
     if (resVector[1]<0) resVector[1]+=18;
-    resVector[2]= cl->GetX();
+    resVector[2]= TMath::Sqrt(cl->GetX()*cl->GetX()+cl->GetY()*cl->GetY());
     resVector[3]= cl->GetZ()/cl->GetX();
     //
     resVector[0]= cl->GetY()-trackOut.GetY();
@@ -2772,7 +2772,7 @@ void AliTPCcalibAlign::UpdateClusterDeltaField(const AliTPCseed * seed){
     trackIn.GetXYZ(xyz);
     resVector[1]= 9.*TMath::ATan2(xyz[1],xyz[0])/TMath::Pi();
     if (resVector[1]<0) resVector[1]+=18;
-    resVector[2]= cl->GetX();
+    resVector[2]= TMath::Sqrt(cl->GetX()*cl->GetX()+cl->GetY()*cl->GetY());
     resVector[3]= cl->GetZ()/cl->GetX();
     //
     resVector[0]= cl->GetY()-trackIn.GetY();
@@ -2963,7 +2963,7 @@ void  AliTPCcalibAlign::UpdateAlignSector(const AliTPCseed * track,Int_t isec){
     if (fCurrentTrack->IsOn(0x4)&&TMath::Abs(vImpact[0])<1&&TMath::Abs(vImpact[1])<1){    
       Double_t resVector[5];
       resVector[1]= 9.*gphi/TMath::Pi();
-      resVector[2]= c->GetX();
+      resVector[2]= TMath::Sqrt(c->GetX()*c->GetX()+c->GetY()*c->GetY());
       resVector[3]= c->GetZ()/c->GetX();
       //
       //