KVedaLossInverseRangeFunction.h

 
#ifndef __KVEDALOSSINVERSERANGEFUNCTION_H
#define __KVEDALOSSINVERSERANGEFUNCTION_H
 
#include "TObject.h"
#include "KVConfig.h"
#include "TSpline.h"
class TF1;
 
class KVedaLossInverseRangeFunction : public TObject {
 
#ifdef USING_ROOT5
    TSpline3* fInterpol;
#else
    std::unique_ptr<TSpline3> fInterpol;
#endif
 
public:
   KVedaLossInverseRangeFunction()
#ifdef USING_ROOT5
      : fInterpol(nullptr)
#endif
   {}
   KVedaLossInverseRangeFunction(TF1* range_func,
                                 Int_t A,
                                 Double_t riso,
                                 Int_t ninter = 50);
   virtual ~KVedaLossInverseRangeFunction()
   {
#ifdef USING_ROOT5
      SafeDelete(fInterpol);
#endif
   }
 
   Double_t GetEnergyPerNucleon(Double_t range, Double_t riso);
 
   ClassDef(KVedaLossInverseRangeFunction, 0) //Dedicated optimised inversion of range-energy function
};
 
#endif