TFCSEnergyInterpolationSpline.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_FASTCALOSIMEVENT_TFCSEnergyInterpolationSpline_h
#define ISF_FASTCALOSIMEVENT_TFCSEnergyInterpolationSpline_h
 
#include "ISF_FastCaloSimEvent/TFCSParametrization.h"
#include "TSpline.h"
 
class TGraph;
 
class TFCSEnergyInterpolationSpline : public TFCSParametrization {
public:
  TFCSEnergyInterpolationSpline(const char *name = nullptr,
                                const char *title = nullptr);

  enum FCSEnergyInitializationStatusBits {
    kOnlyScaleEnergy =
        BIT(15) 
  };
 
  bool OnlyScaleEnergy() const { return TestBit(kOnlyScaleEnergy); };
  void set_OnlyScaleEnergy() { SetBit(kOnlyScaleEnergy); };
  void reset_OnlyScaleEnergy() { ResetBit(kOnlyScaleEnergy); };
 
  virtual bool is_match_Ekin_bin(int /*Ekin_bin*/) const override {
    return true;
  };
  virtual bool is_match_calosample(int /*calosample*/) const override {
    return true;
  };

  void InitFromSpline(const TSpline3 &spline) { m_spline = spline; };

  void InitFromArrayInLogEkin(Int_t np, Double_t logEkin[], Double_t response[],
                              const char *opt = nullptr, Double_t valbeg = 0,
                              Double_t valend = 0);

  void InitFromArrayInEkin(Int_t np, Double_t Ekin[], Double_t response[],
                           const char *opt = nullptr, Double_t valbeg = 0,
                           Double_t valend = 0);
 
  const TSpline3 &spline() const { return m_spline; };

  virtual FCSReturnCode
  simulate(TFCSSimulationState &simulstate, const TFCSTruthState *truth,
           const TFCSExtrapolationState *extrapol) const override;
 
  void Print(Option_t *option = "") const override;
 
  static void unit_test(TFCSSimulationState *simulstate = nullptr,
                        TFCSTruthState *truth = nullptr,
                        const TFCSExtrapolationState *extrapol = nullptr,
                        TGraph *grspline = nullptr);
 
private:
  TSpline3 m_spline;
 
  ClassDefOverride(TFCSEnergyInterpolationSpline,
                   1) // TFCSEnergyInterpolationSpline
};
 
#endif