KVRangeYanezMaterial.h

 
#ifndef __KVRANGEYANEZMATERIAL_H
#define __KVRANGEYANEZMATERIAL_H
#include "KVIonRangeTableMaterial.h"
 
class range;  // PIMPL C++ implementation of range library
 
class KVRangeYanezMaterial : public KVIonRangeTableMaterial {
    Int_t fTableType;
 
    Int_t fNelem;
    Int_t iabso;
   struct Elem {
      Double_t w;
      Int_t z;
      Int_t a;
    } fAbsorb[10];
    Double_t thickness;
    Int_t Zp, Ap;
    Double_t error;
 
    std::unique_ptr<range> range_lib;
protected:
   void PrepareRangeLibVariables(Int_t Z, Int_t A);
   Double_t DeltaEFunc(Double_t*, Double_t*);
   Double_t RangeFunc(Double_t*, Double_t*);
   Double_t EResFunc(Double_t*, Double_t*);
 
   void MakeFunctionObjects();
 
public:
   ~KVRangeYanezMaterial();
   KVRangeYanezMaterial();
   KVRangeYanezMaterial(const KVRangeYanezMaterial&) ;
   KVRangeYanezMaterial(const KVIonRangeTable* t, const Char_t* name, const Char_t* symbol, const Char_t* state,
                        Double_t density = -1, Int_t Z = -1, Int_t A = -1);
 
   void Copy(TObject&) const override;
   void Initialize() override;
 
   TF1* GetDeltaEFunction(Double_t e, Int_t Z, Int_t A, Double_t isoAmat = 0) override;
   TF1* GetRangeFunction(Int_t Z, Int_t A, Double_t isoAmat = 0) override;
   TF1* GetEResFunction(Double_t e, Int_t Z, Int_t A, Double_t isoAmat = 0) override;
   TF1* GetStoppingFunction(Int_t, Int_t, Double_t = 0) override
   {
      KVError::Warning(this, "GetStoppingFunction", "Not available for Yanez Range tables");
      return 0;
   };
   void SetTableType(int type)
   {
      fTableType = type;
   };
 
   Double_t GetEIncFromEResOfIon(Int_t Z, Int_t A, Double_t Eres, Double_t e, Double_t isoAmat = 0.) override;
 
   Int_t GetNElem() const
   {
      return fNelem;
   }
   void SaveMaterial(std::ofstream& matfile);
 
   ClassDefOverride(KVRangeYanezMaterial, 1) //Description of absorber for the Range dE/dx and range library (Ricardo Yanez)
};
 
#endif