KVLightEnergyCsI.h

#ifndef KV_LIGHT_ENERGY_CSI_H
#define KV_LIGHT_ENERGY_CSI_H
 
#include "KVCalibrator.h"
#include "KVValueRange.h"
 
class KVLightEnergyCsI: public KVCalibrator {
 
   Double_t CalculLumiere(Double_t*, Double_t*);
   Double_t InvCalc(Double_t*, Double_t*);
   std::unique_ptr<TF1> inv_fit_func;
   Double_t a3_fit = 6.;
   Double_t a4_fit = 0.385;
   std::optional<KVValueRange<double>> a3_range;
   std::optional<KVValueRange<double>> a4_range;
 
   Bool_t CheckParameters(const KVNameValueList& z_and_a) const;
 
protected:
   mutable Double_t Z;
   mutable Double_t A;
 
public:
   KVLightEnergyCsI(Bool_t make_func = kTRUE);
   Double_t Compute(Double_t chan, const KVNameValueList& z_and_a = "") const override;
   Double_t Invert(Double_t, const KVNameValueList& z_and_a = "") const override;
 
   void Fit(TH1*, int z, int a, std::optional<KVValueRange<Double_t>> fit_e_range = {});
   void invFit(TH1*, int z, int a, std::optional<KVValueRange<Double_t>> fit_lite_range = {});
 
   void SetZandA(double z, double a)
   {
      Z = z; A = a;
   }
   void fit_with_a3_fixed(double a3)
   {
      a3_fit = a3;
      a3_range.reset();//remove any previously set range
   }
   void fit_with_a3_range(double a3_min, double a3_max)
   {
      a3_range = {a3_min,a3_max};
   }
   void fit_with_a4_fixed(double a4)
   {
      a4_fit = a4;
      a4_range.reset();//remove any previously set range
   }
   void fit_with_a4_range(double a4_min, double a4_max)
   {
      a4_range = {a4_min,a4_max};
   }
 
   ClassDefOverride(KVLightEnergyCsI, 1)        //Light-energy calibration for CsI detectors
};
 
#endif