*/
-
+#include "AliTPCcalibDB.h"
#include "TLinearFitter.h"
#include "Riostream.h"
#include "TList.h"
// 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;
//
// 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;
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
+// }
//
//
"cl.="<<cluster<<
"cy="<<dy<<
"cz="<<dz<<
- "cY="<<corrclY<<
- "cR="<<corrR<<
+ // "cY="<<corrclY<<
+ // "cR="<<corrR<<
"dxq="<<dxq<<
"dyq="<<dyq<<
"dzq="<<dzq<<