CbmStsPhysics.h

Go to the documentation of this file.

/* Copyright (C) 2014-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
#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
#include <utility>
#include <vector>
 
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);
 
  std::pair<std::vector<double>, double> GetLandauWidthTable() const;
 
  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_ */