CbmStsPhysics.h

Go to the documentation of this file.

 
#ifndef CBMSTSPHYSICS_H
#define CBMSTSPHYSICS_H 1
 
#include <Rtypes.h>      // for THashConsistencyHolder, ClassDef
#include <RtypesCore.h>  // for Double_t, Int_t, Bool_t
 
#include <map>  // for map
 
class CbmStsPhysics {
 
public:
  virtual ~CbmStsPhysics();
 
 
  static Double_t DiffusionWidth(Double_t z,
                                 Double_t d,
                                 Double_t vBias,
                                 Double_t vFd,
                                 Double_t temperature,
                                 Int_t chargeType);
 
 
  static Double_t
  ElectricField(Double_t vBias, Double_t vFd, Double_t dZ, Double_t z);
 
 
  Double_t
  EnergyLoss(Double_t dz, Double_t mass, Double_t eKin, Double_t dedx) const;
 
 
  static CbmStsPhysics* Instance();
 
 
  Double_t LandauWidth(Double_t mostProbableCharge);
 
 
  static Double_t PairCreationEnergy() { return 3.57142e-9; }
 
 
  static Double_t ParticleCharge(Int_t pid);
 
 
  static Double_t ParticleMass(Int_t pid);
 
 
  Double_t StoppingPower(Double_t eKin, Int_t pid);
 
 
  Double_t StoppingPower(Double_t energy,
                         Double_t mass,
                         Double_t charge,
                         Bool_t isElectron);
 
 
private:
  CbmStsPhysics();
 
 
  CbmStsPhysics(const CbmStsPhysics&) = delete;
  ;
 
 
  CbmStsPhysics operator=(const CbmStsPhysics&) = delete;
 
 
  Double_t InterpolateDataTable(Double_t eKin,
                                std::map<Double_t, Double_t>& table);
 
 
  void ReadDataTablesStoppingPower();
 
 
  void ReadDataTablesLandauWidth();
 
 
  void SetUrbanParameters(Double_t z);
 
 
private:
  static CbmStsPhysics* fgInstance;  
 
  // --- Parameters for the Urban model
  Double_t fUrbanI = 0.;  
  Double_t fUrbanE1 = 0.;  
  Double_t fUrbanE2 = 0.;  
  Double_t fUrbanF1 = 0.;  
  Double_t fUrbanF2 = 0.;  
  Double_t fUrbanEmax =
    0.;  
  Double_t fUrbanR =
    0.;  
 
  // --- Data tables for stopping power
  std::map<Double_t, Double_t>
    fStoppingElectron {};  
  std::map<Double_t, Double_t>
    fStoppingProton {};  
 
  // --- Data tables for width of Landau distribution
  std::map<Double_t, Double_t> fLandauWidth;  
 
 
  ClassDef(CbmStsPhysics, 2);
};
 
#endif /* CBMSTSPHYSICS_H_ */