KVIonRangeTableMaterial.h

 
#ifndef __KVIONRANGETABLEMATERIAL_H
#define __KVIONRANGETABLEMATERIAL_H
 
#include "KVBase.h"
#include <TString.h>
#include <KVList.h>
#include "KVIonRangeTable.h"
 
class TGeoMaterial;
class TF1;
 
#define RTT  62.36367e+03  // cm^3.Torr.K^-1.mol^-1
#define ZERO_KELVIN  273.15
 
class KVIonRangeTableMaterial : public KVBase {
protected:
    const KVIonRangeTable* fTable;
    TString fState;
    KVList* fComposition;
   Bool_t fCompound;
   Bool_t fMixture;
    Double_t fDens;
    Double_t fZmat;
    Double_t fAmat;
    Double_t fMoleWt;
 
    Double_t fRangeOfLastDE;
 
    TF1* fDeltaE;
    TF1* fEres;
    TF1* fRange;
    TF1* fStopping;
 
public:
   KVIonRangeTableMaterial();
   KVIonRangeTableMaterial(const KVIonRangeTable*, const Char_t* name, const Char_t* symbol, const Char_t* state,
                           Double_t density = -1, Double_t Z = -1, Double_t A = -1);
   KVIonRangeTableMaterial(const KVIonRangeTableMaterial&) ;
   virtual ~KVIonRangeTableMaterial();
 
   virtual void Initialize();
   const KVIonRangeTable* GetTable() const
   {
      return fTable;
   }
 
   void ls(Option_t* = "") const;
   void Print(Option_t* = "") const;
   virtual Float_t GetEmaxValid(Int_t /* Z */, Int_t /* A */) const
   {
      return 1.e+4;
   };
   virtual Float_t GetEminValid(Int_t /* Z */, Int_t /* A */) const
   {
      return 0.0;
   };
 
   void Copy(TObject&) const;
   const Char_t* GetSymbol() const
   {
      return GetTitle();
   };
   Double_t GetMass() const
   {
      return fAmat;
   };
   Double_t GetZ() const
   {
      return fZmat;
   };
   Double_t GetMoleWt() const
   {
      return fMoleWt;
   };
   Double_t GetDensity() const
   {
      return fDens;
   };
   void SetDensity(Double_t d)
   {
      fDens = d;
   };
   void SetTemperatureAndPressure(Double_t T, Double_t P)
   {
      if (IsGas() && P == 0) fDens = 0;
      if (IsGas() && T > 0 && P > 0) fDens = (fMoleWt * P) / ((T + ZERO_KELVIN) * RTT);
   };
   void SetState(const Char_t* st)
   {
      fState = st;
   };
   Bool_t IsGas() const
   {
      return (fState == "gas");
   };
 
   Bool_t IsCompound() const
   {
      return (fCompound);
   };
   Bool_t IsMixture() const
   {
      return (fMixture);
   };
 
   KVList* GetComposition() const
   {
      return fComposition;
   };
 
   void AddCompoundElement(Int_t Z, Int_t A, Int_t Natoms);
   void AddMixtureElement(Int_t Z, Int_t A, Int_t Natoms, Double_t Proportion);
 
   TGeoMaterial* GetTGeoMaterial() const;
 
   virtual TF1* GetRangeFunction(Int_t Z, Int_t A, Double_t isoAmat = 0) = 0;
   virtual TF1* GetDeltaEFunction(Double_t e, Int_t Z, Int_t A, Double_t isoAmat = 0) = 0;
   virtual TF1* GetEResFunction(Double_t e, Int_t Z, Int_t A, Double_t isoAmat = 0) = 0;
   virtual TF1* GetStoppingFunction(Int_t Z, Int_t A, Double_t isoAmat = 0) = 0;
 
   void PrintRangeTable(Int_t Z, Int_t A, Double_t isoAmat = 0, Double_t units = KVUnits::cm, Double_t T = -1, Double_t P = -1);
   void PrintComposition(std::ostream&) const;
 
   virtual Double_t GetRangeOfIon(Int_t Z, Int_t A, Double_t E, Double_t isoAmat = 0.);
   virtual Double_t GetLinearRangeOfIon(Int_t Z, Int_t A, Double_t E, Double_t isoAmat = 0, Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetDeltaEOfIon(Int_t Z, Int_t A, Double_t E, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetLinearDeltaEOfIon(Int_t Z, Int_t A, Double_t E, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetEResOfIon(Int_t Z, Int_t A, Double_t E, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetLinearEResOfIon(Int_t Z, Int_t A, Double_t E, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   virtual Double_t GetEIncFromEResOfIon(Int_t Z, Int_t A, Double_t Eres, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetLinearEIncFromEResOfIon(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(Int_t Z, Int_t A, Double_t DeltaE, Double_t e, enum KVIonRangeTable::SolType type = KVIonRangeTable::kEmax, Double_t isoAmat = 0.);
   virtual Double_t GetLinearEIncFromDeltaEOfIon(Int_t Z, Int_t A, Double_t DeltaE, Double_t e, enum KVIonRangeTable::SolType type = KVIonRangeTable::kEmax, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   Double_t GetDeltaEFromEResOfIon(Int_t Z, Int_t A, Double_t ERes, Double_t e, Double_t isoAmat = 0.)
   {
      return GetEIncFromEResOfIon(Z, A, ERes, e, isoAmat) - ERes;
   }
 
   Double_t GetLinearDeltaEFromEResOfIon(Int_t Z, Int_t A, Double_t ERes, Double_t e,
                                         Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.)
   {
      return GetLinearEIncFromEResOfIon(Z, A, ERes, e, isoAmat, T, P) - ERes;
   }
 
   virtual Double_t GetPunchThroughEnergy(Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0.) ;
   virtual Double_t GetLinearPunchThroughEnergy(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(Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetEIncOfMaxDeltaEOfIon(Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetLinearMaxDeltaEOfIon(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(Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0., Double_t T = -1., Double_t P = -1.);
 
   Double_t GetRangeOfLastDE() const
   {
      return fRangeOfLastDE;
   }
 
   ClassDef(KVIonRangeTableMaterial, 1) //Material for use in energy loss & range calculations
};
 
#endif