M AliTPCcalibCalib.cxx - Disabled PointCorrection which is not yet...
authormarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Mar 2010 14:24:16 +0000 (14:24 +0000)
committermarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Mar 2010 14:24:16 +0000 (14:24 +0000)
M      AliTPCPointCorrection.cxx   - Access correction  parameters from the OCDB

Marian

TPC/AliTPCPointCorrection.cxx
TPC/AliTPCcalibCalib.cxx

index fb46aaa..25c0513 100644 (file)
@@ -28,7 +28,7 @@
     
 */
 
-
+#include "AliTPCcalibDB.h"
 #include "TLinearFitter.h"
 #include "Riostream.h"
 #include "TList.h"
@@ -229,7 +229,11 @@ Double_t AliTPCPointCorrection::RPhiCOGCorrection(Int_t sector, Int_t padrow, Fl
   //  qMax  - cluster max charge - cl-.GetMax()
   //  threshold - clusterer threshold
   //
-  AliTPCClusterParam * clparam = AliTPCClusterParam::Instance(); 
+  AliTPCClusterParam * clparam = AliTPCcalibDB::Instance()->GetClusterParam();
+  if (!clparam) {
+    AliFatal("TPC OCDB not initialized"); 
+    return 0;
+  }
   Int_t padtype=0;
   if (sector>=36) padtype = (padrow>64)?2:1;
   Double_t padwidth=(padtype==0)? 0.4:0.6;
index 088cef6..63edd33 100644 (file)
@@ -195,7 +195,7 @@ Bool_t  AliTPCcalibCalib::RefitTrack(AliESDtrack * track, AliTPCseed *seed, Floa
   //
   // First apply calibration
   //
-  AliTPCPointCorrection * corr =  AliTPCPointCorrection::Instance();
+  //  AliTPCPointCorrection * corr =  AliTPCPointCorrection::Instance();
   for (Int_t irow=0;irow<159;irow++) {
     AliTPCclusterMI *cluster=seed->GetClusterPointer(irow);
     if (!cluster) continue; 
@@ -229,22 +229,22 @@ Bool_t  AliTPCcalibCalib::RefitTrack(AliESDtrack * track, AliTPCseed *seed, Floa
       x[1]-=dyq;
       x[2]-=dzq;
     }
-    //
-    // Apply r-phi correction  - To be done on track level- knowing the track angle !!!
-    //
-    Double_t corrclY =  
-      corr->RPhiCOGCorrection(cluster->GetDetector(),cluster->GetRow(), cluster->GetPad(),
-                                 cluster->GetY(),cluster->GetY(), cluster->GetZ(), 0., cluster->GetMax(),2.5);
-    // R correction
-    Double_t corrR   = corr->CorrectionOutR0(kFALSE,kFALSE,cluster->GetX(),cluster->GetY(),cluster->GetZ(),cluster->GetDetector());
-
-    if (0&fApplyRPhiCorrection){
-      if (cluster->GetY()>0) x[1]+=corrclY;  // rphi correction
-      if (cluster->GetY()<0) x[1]-=corrclY;  // rphi correction
-    }
-    if (0&fApplyRCorrection){      
-      x[0]+=corrR;                           // radial correction
-    }
+//     //
+//     // Apply r-phi correction  - To be done on track level- knowing the track angle !!!
+//     //
+//     Double_t corrclY =  
+//       corr->RPhiCOGCorrection(cluster->GetDetector(),cluster->GetRow(), cluster->GetPad(),
+//                               cluster->GetY(),cluster->GetY(), cluster->GetZ(), 0., cluster->GetMax(),2.5);
+//     // R correction
+//     Double_t corrR   = corr->CorrectionOutR0(kFALSE,kFALSE,cluster->GetX(),cluster->GetY(),cluster->GetZ(),cluster->GetDetector());
+
+//     if (0&fApplyRPhiCorrection){
+//       if (cluster->GetY()>0) x[1]+=corrclY;  // rphi correction
+//       if (cluster->GetY()<0) x[1]-=corrclY;  // rphi correction
+//     }
+//     if (0&fApplyRCorrection){      
+//       x[0]+=corrR;                           // radial correction
+//     }
 
     //
     //
@@ -287,8 +287,8 @@ Bool_t  AliTPCcalibCalib::RefitTrack(AliESDtrack * track, AliTPCseed *seed, Floa
          "cl.="<<cluster<<
          "cy="<<dy<<
          "cz="<<dz<<
-         "cY="<<corrclY<<
-         "cR="<<corrR<<
+         //      "cY="<<corrclY<<
+         //      "cR="<<corrR<<
          "dxq="<<dxq<<
          "dyq="<<dyq<<
          "dzq="<<dzq<<