]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/AliTPCCorrection.h
Test for Coverity
[u/mrichter/AliRoot.git] / TPC / AliTPCCorrection.h
index 679da2a8007a043a7e568302b80e7e00ea9d6eb4..d3518e60359db15f4319ed498634872c5b4c17d1 100644 (file)
@@ -12,6 +12,7 @@
 
 #include <TNamed.h>
 #include "TMatrixD.h"
+#include "TMatrixF.h"
 class TH2F;
 class TTimeStamp;
 class TCollection;
@@ -60,7 +61,7 @@ public:
 
 
   TTree* CreateDistortionTree(Double_t step=5);
-  static void  MakeDistortionMap(THnSparse * his0, TTreeSRedirector *pcstream, const char* hname, Int_t run,  Float_t refX, Int_t type);
+  static void  MakeDistortionMap(THnSparse * his0, TTreeSRedirector *pcstream, const char* hname, Int_t run,  Float_t refX, Int_t type, Int_t integ=1);
   static void  MakeDistortionMapCosmic(THnSparse * his0, TTreeSRedirector *pcstream, const char* hname, Int_t run,  Float_t refX, Int_t type);
   static void  MakeDistortionMapSector(THnSparse * his0, TTreeSRedirector *pcstream, const char* hname, Int_t run, Int_t type);
   // normally called directly in the correction classes which inherit from this class
@@ -105,12 +106,14 @@ protected:
 
   // Simple Interpolation functions: e.g. with tricubic interpolation (not yet in TH3)
   Int_t fILow, fJLow, fKLow;          // variable to help in the interpolation 
+  // Double_t versions
   void Interpolate2DEdistortion( const Int_t order, const Double_t r, const Double_t z, 
                                 const Double_t er[kNZ][kNR], Double_t &erValue );
   void Interpolate3DEdistortion( const Int_t order, const Double_t r, const Float_t phi, const Double_t z, 
                                 const Double_t er[kNZ][kNPhi][kNR], const Double_t ephi[kNZ][kNPhi][kNR], 
                                 const Double_t ez[kNZ][kNPhi][kNR],
                                 Double_t &erValue, Double_t &ephiValue, Double_t &ezValue);
+  // TMatrixD versions (for e.g. Poisson relaxation)
   Double_t Interpolate2DTable( const Int_t order, const Double_t x, const Double_t y, 
                              const Int_t nx,  const Int_t ny, const Double_t xv[], const Double_t yv[], 
                              const TMatrixD &array );
@@ -121,6 +124,18 @@ protected:
   Double_t Interpolate( const Double_t xArray[], const Double_t yArray[], 
                        const Int_t order, const Double_t x );
   void Search( const Int_t n, const Double_t xArray[], const Double_t x, Int_t &low );
+  // TMatrixF versions (smaller size, e.g. for final look up table)
+  Float_t Interpolate2DTable( const Int_t order, const Double_t x, const Double_t y, 
+                             const Int_t nx,  const Int_t ny, const Double_t xv[], const Double_t yv[], 
+                             const TMatrixF &array );
+  Float_t Interpolate3DTable( const Int_t order, const Double_t x,   const Double_t y,   const Double_t z,
+                             const Int_t  nx,    const Int_t  ny,    const Int_t  nz,
+                             const Double_t xv[], const Double_t yv[], const Double_t zv[],
+                              TMatrixF **arrayofArrays ); 
+  Float_t Interpolate( const Double_t xArray[], const Float_t yArray[], 
+                       const Int_t order, const Double_t x );
+
   virtual Int_t IsPowerOfTwo ( Int_t i ) const  ;