KVIonRangeTable.h

 
#ifndef KVIONRANGETABLE_H
#define KVIONRANGETABLE_H
 
#include "KVBase.h"
#include "KVUnits.h"
 
class TGeoMaterial;
class KVIonRangeTableMaterial;
class TGeoManager;
class TVector3;
 
class KVIonRangeTable : public KVBase {
 
protected:
   virtual KVIonRangeTableMaterial* GetMaterialWithPointer(TGeoMaterial*) const;
   virtual KVIonRangeTableMaterial* GetMaterialWithNameOrType(const Char_t* material) const = 0;
 
public:
   enum SolType {
      kEmax,
      kEmin
   };
    std::optional<double> get_element_density(int z) const;
    bool is_gas(int z) const;
    std::string get_element_name(int z) const;
    std::string get_element_symbol(int z) const;
 
 
   KVIonRangeTable(const Char_t* name = "", const Char_t* title = "");
   virtual ~KVIonRangeTable();
 
   static KVIonRangeTable* GetRangeTable(const Char_t* name);
 
   virtual KVIonRangeTableMaterial* AddElementalMaterial(Int_t z, Int_t a = 0) const
   {
 
      IGNORE_UNUSED(z);
      IGNORE_UNUSED(a);
      return nullptr;
   }
   struct CompoundFormulaElement {
      Int_t Z;
      Int_t A;
      Int_t natoms;
   };
 
   virtual KVIonRangeTableMaterial* AddCompoundMaterial(
      const Char_t* name, const Char_t* symbol,
      const std::vector<CompoundFormulaElement>& elements, Double_t density = -1.0) const
   {
 
      IGNORE_UNUSED(name);
      IGNORE_UNUSED(symbol);
      IGNORE_UNUSED(elements);
      IGNORE_UNUSED(density);
      return nullptr;
   }
   struct MixtureFormulaElement : CompoundFormulaElement {
      Double_t proportion;
   };
   virtual KVIonRangeTableMaterial* AddMixedMaterial(
      const Char_t*  name, const Char_t*  symbol,
      const std::vector<MixtureFormulaElement>& elements,
      Double_t  density  = -1.0) const
   {
      IGNORE_UNUSED(name);
      IGNORE_UNUSED(symbol);
      IGNORE_UNUSED(elements);
      IGNORE_UNUSED(density);
      return nullptr;
   }
 
   virtual TObjArray* GetListOfMaterials() = 0;
 
   virtual TGeoMaterial* GetTGeoMaterial(const Char_t* material);
 
   virtual Double_t GetEmaxValid(const Char_t* material, Int_t Z, Int_t A);
 
   virtual Double_t GetDensity(const Char_t*);
 
   virtual void SetDensity(const Char_t*, Double_t);
 
   virtual void SetTemperatureAndPressure(const Char_t*, Double_t temperature, Double_t pressure);
 
   virtual Double_t GetZ(const Char_t*);
 
   virtual Double_t GetAtomicMass(const Char_t*);
 
   KVIonRangeTableMaterial* GetMaterial(const Char_t* material) const;
   KVIonRangeTableMaterial* GetMaterial(TGeoMaterial*) const;
 
   virtual Bool_t IsMaterialKnown(const Char_t*);
 
   virtual Bool_t IsMaterialKnown(TGeoMaterial*);
 
   virtual Bool_t IsMaterialGas(const Char_t*);
 
   const Char_t* GetMaterialName(const Char_t*);
 
   const Char_t* GetMaterialType(const Char_t*);
 
   virtual Double_t GetRangeOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t E,
                                  Double_t Amat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetLinearRangeOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t E,
                                        Double_t Amat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t E, Double_t r,
                                   Double_t Amat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetLinearDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t E, Double_t d,
                                         Double_t Amat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetEResOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t E, Double_t r,
                                 Double_t Amat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearEResOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t E, Double_t d,
                                       Double_t Amat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetEIncFromEResOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t Eres, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearEIncFromEResOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t Eres, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetEIncFromDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t DeltaE, Double_t e, enum SolType type = kEmax, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearEIncFromDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t DeltaE, Double_t e, enum SolType type = kEmax, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetDeltaEFromEResOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t ERes, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearDeltaEFromEResOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t ERes, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetPunchThroughEnergy(const Char_t* mat, Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearPunchThroughEnergy(const Char_t* mat, Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetMaxDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetEIncOfMaxDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearMaxDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
   virtual Double_t GetLinearEIncOfMaxDeltaEOfIon(const Char_t* mat, Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   void Print(Option_t* = "") const override;
 
   virtual Bool_t CheckIon(Int_t, Int_t) const
   {
      AbstractMethod("CheckIon");
      return kTRUE;
   }
   virtual Bool_t ReadMaterials(const Char_t*) const = 0;
 
   ClassDefOverride(KVIonRangeTable, 1) //Abstract base class for calculation of range & energy loss of charged particles in matter
};
 
#endif