Adopt high level TPCFastTransform methods

Author: shahor02

GPU/TPCFastTransformation/TPCFastTransform.h

@@ -250,6 +250,9 @@ class TPCFastTransform : public FlatObject
   /// Return TOF correction (vdrift / C)
   GPUd() float getTOFCorr() const { return mLdriftCorr; }
 
+  /// Return nominal PV Z position correction (vdrift / C)
+  GPUd() float getPrimVtxZ() const { return mPrimVtxZ; }
+
   /// Return map lumi
   GPUd() float getLumi() const { return mLumi; }
 

GPU/TPCFastTransformation/TPCFastTransformPOD.cxx

@@ -15,6 +15,7 @@
 /// \author  ruben.shahoayn@cern.ch
 
 #include "TPCFastTransformPOD.h"
+#include "GPUDebugStreamer.h"
 #if !defined(GPUCA_GPUCODE)
 #include <TRandom.h>
 #endif
@@ -23,6 +24,7 @@ namespace GPUCA_NAMESPACE
 {
 namespace gpu
 {
+
 #if !defined(GPUCA_GPUCODE)
 
 size_t TPCFastTransformPOD::estimateSize(const TPCFastTransform& src)
@@ -63,6 +65,7 @@ const TPCFastTransformPOD& TPCFastTransformPOD::create(char* buff, size_t buffSi
   assert(buffSize > sizeof(TPCFastTransformPOD));
   const auto& origCorr = src.getCorrection();
   auto& podMap = getNonConst(buff);
+  podMap.mApplyCorrections = true; // by default always apply corrections
 
   // copy fixed size data --- start
   podMap.mNumberOfScenarios = origCorr.mNumberOfScenarios;
@@ -78,6 +81,14 @@ const TPCFastTransformPOD& TPCFastTransformPOD::create(char* buff, size_t buffSi
   }
   podMap.mInterpolationSafetyMargin = origCorr.fInterpolationSafetyMargin;
   podMap.mTimeStamp = origCorr.mTimeStamp;
+
+  podMap.mVdrift = src.getVDrift();
+  podMap.mT0 = src.getT0();
+  podMap.mVdriftCorrY = src.getVdriftCorrY();
+  podMap.mLdriftCorr = src.getLdriftCorr();
+  podMap.mTOFcorr = src.getTOFCorr();
+  podMap.mPrimVtxZ = src.getPrimVtxZ();
+
   // copy fixed size data --- end
 
   size_t nextDynOffs = alignOffset(sizeof(TPCFastTransformPOD));
@@ -207,10 +218,10 @@ bool TPCFastTransformPOD::test(const TPCFastTransform& src, int npoints) const
     }
   }
   //
-  LOGP(info, " (ns per call)              original       this      Nmissmatch");
-  LOGP(info, "getCorrection               {:.3e}    {:.3e}   {}", double(origEnd[0] - origStart[0]) / npoints, double(thisEnd[0] - thisStart[0]) / npoints, ndiff[0]);
-  LOGP(info, "getCorrectionInvCorrectedX  {:.3e}    {:.3e}   {}", double(origEnd[1] - origStart[1]) / npoints, double(thisEnd[1] - thisStart[1]) / npoints, ndiff[1]);
-  LOGP(info, "getCorrectionInvUV          {:.3e}    {:.3e}   {}", double(origEnd[2] - origStart[2]) / npoints, double(thisEnd[2] - thisStart[2]) / npoints, ndiff[2]);
+  LOGP(info, " (ns per call)              original        this     Nmissmatch");
+  LOGP(info, "getCorrection               {:.3e}    {:.3e}   {}", double(origEnd[0] - origStart[0]) / npoints * 1000., double(thisEnd[0] - thisStart[0]) / npoints * 1000., ndiff[0]);
+  LOGP(info, "getCorrectionInvCorrectedX  {:.3e}    {:.3e}   {}", double(origEnd[1] - origStart[1]) / npoints * 1000., double(thisEnd[1] - thisStart[1]) / npoints * 1000., ndiff[1]);
+  LOGP(info, "getCorrectionInvUV          {:.3e}    {:.3e}   {}", double(origEnd[2] - origStart[2]) / npoints * 1000., double(thisEnd[2] - thisStart[2]) / npoints * 1000., ndiff[2]);
   return ndiff[0] == 0 && ndiff[1] == 0 && ndiff[2] == 0;
 }