KVedaLossMaterial.h

 
 
#ifndef __KVEDALOSSMATERIAL_H
#define __KVEDALOSSMATERIAL_H
 
#include "KVIonRangeTableMaterial.h"
#include <Riostream.h>
#include "TObjArray.h"
#include <vector>
 
class TGeoMaterial;
class KVedaLoss;
 
#define ZMAX_VEDALOSS 100
#define VERY_BIG_ENERGY  1.e+6
 
class KVedaLossMaterial : public KVIonRangeTableMaterial {
 
private:
   static KVedaLoss* fgTable;
   Double_t RF_Z;
   Double_t RF_A;
   std::vector<Double_t>* par;
   Double_t ran, adm, dleps, adn, riso, eps, DLEP, drande;
    Double_t thickness;
    TObjArray fInvRange;
 
protected:
    std::vector<Double_t> fEmin;
    std::vector<Double_t> fEmax;
    std::vector< std::vector<Double_t> > fCoeff;
 
    static Bool_t fNoLimits;
 
   Double_t DeltaEFunc(Double_t*, Double_t*);
   Double_t EResFunc(Double_t*, Double_t*);
   Double_t RangeFunc(Double_t*, Double_t*);
   Double_t StoppingFunc(Double_t*, Double_t*);
 
public:
   KVedaLossMaterial();
   KVedaLossMaterial(const KVIonRangeTable*, const Char_t* name, const Char_t* type, const Char_t* state,
                     Double_t density, Double_t Z, Double_t A, Double_t MoleWt = 0.0);
   virtual ~KVedaLossMaterial();
 
   Bool_t ReadRangeTable(FILE* fp);
   Float_t GetEmaxValid(Int_t Z, Int_t A) const
   {
      return (CheckIon(Z) ? A * fEmax[Z - 1] : 0.0);
   };
   Float_t GetEminValid(Int_t Z, Int_t A) const
   {
      return (CheckIon(Z) ? A * fEmin[Z - 1] : 0.0);
   };
 
   virtual TF1* GetRangeFunction(Int_t Z, Int_t A, Double_t isoAmat = 0);
   virtual TF1* GetDeltaEFunction(Double_t e, Int_t Z, Int_t A, Double_t isoAmat = 0);
   virtual TF1* GetEResFunction(Double_t e, Int_t Z, Int_t A, Double_t isoAmat = 0);
   virtual TF1* GetStoppingFunction(Int_t Z, Int_t A, Double_t isoAmat = 0);
 
   virtual Double_t GetRangeOfIon(Int_t Z, Int_t A, Double_t E, Double_t isoAmat = 0.);
   virtual Double_t GetDeltaEOfIon(Int_t Z, Int_t A, Double_t E, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetEResOfIon(Int_t Z, Int_t A, Double_t E, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetPunchThroughEnergy(Int_t Z, Int_t A, Double_t e, Double_t isoAmat = 0.);
   virtual Double_t GetEIncFromEResOfIon(Int_t Z, Int_t A, Double_t Eres, Double_t e, Double_t isoAmat = 0.);
 
   static void SetNoLimits(Bool_t on = kTRUE)
   {
      fNoLimits = on;
   };
 
   void GetParameters(Int_t Zion, Int_t& Aion, std::vector<Double_t>& rangepar);
   static Bool_t CheckIon(Int_t Z)
   {
      return (Z > 0 && Z <= ZMAX_VEDALOSS);
   }
 
   ClassDef(KVedaLossMaterial, 4) //Description of material properties used by KVedaLoss range calculation
};
 
#endif