KVedaLoss.h

 
#ifndef __KVEDALOSS_H
#define __KVEDALOSS_H
 
#include "KVIonRangeTable.h"
#include "KVHashList.h"
 
class KVedaLossMaterial;
class TGeoMaterial;
 
class KVedaLoss : public KVIonRangeTable {
    static KVHashList* fMaterials;
   TString fLocalMaterialsDirectory;
 
   Bool_t init_materials() const;
   Bool_t CheckMaterialsList() const
   {
      if (!fMaterials) return init_materials();
      return kTRUE;
   };
   KVIonRangeTableMaterial* GetMaterialWithNameOrType(const Char_t* material) const override;
    static Bool_t fgNewRangeInversion;
 
   void AddMaterial(KVIonRangeTableMaterial*) const;
 
public:
   KVedaLoss();
   virtual ~KVedaLoss();
 
 
   void Print(Option_t* = "") const override;
   TObjArray* GetListOfMaterials() override;
 
   static void SetIgnoreEnergyLimits(Bool_t yes = kTRUE);
 
   static void SetUseNewRangeInversion(Bool_t yes = kTRUE)
   {
      fgNewRangeInversion = yes;
   }
   static Bool_t IsUseNewRangeInversion()
   {
      return fgNewRangeInversion;
   }
 
   Bool_t CheckIon(Int_t Z, Int_t A) const override;
   Bool_t ReadMaterials(const Char_t* path) const override;
 
   KVIonRangeTableMaterial* AddElementalMaterial(Int_t Z, Int_t A = 0) const override;
   Bool_t AddRANGEMaterial(const Char_t* name) const;
   KVIonRangeTableMaterial* AddCompoundMaterial(
      const Char_t* /*name*/, const Char_t* /* symbol */,
      const std::vector<CompoundFormulaElement>& /* elements */, Double_t /* density */ = -1.0) const override;
   KVIonRangeTableMaterial* AddMixedMaterial(
      const Char_t* /* name */, const Char_t* /* symbol */,
      const std::vector<MixtureFormulaElement>& /* elements */,
      Double_t /* density */ = -1.0) const override;
   ClassDefOverride(KVedaLoss, 1) //C++ implementation of VEDALOSS stopping power calculation
};
 
#endif