BmnSsdPhysics.h

Go to the documentation of this file.

 
#ifndef BMNSSDPHYSICS_H
#define BMNSSDPHYSICS_H 1
 
 
#include "Rtypes.h"
#include "TObject.h"
#include <iostream>
#include <map>
 
class BmnSsdPhysics : public TObject {
 
  public:
 
    virtual ~BmnSsdPhysics();
 
 
 
    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 Double_t GetSiCharge() { return fgkSiCharge; }
 
 
    static BmnSsdPhysics* 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:
 
    BmnSsdPhysics();
 
    BmnSsdPhysics(const BmnSsdPhysics&) = delete;
 
    BmnSsdPhysics operator=(const BmnSsdPhysics&) = delete;
 
    static BmnSsdPhysics* fgInstance;
 
    // --- Physical constants
    static const Double_t fgkSiCharge;     
    static const Double_t fgkSiDensity;    
    static const Double_t fgkProtonMass;   
 
    // --- Parameters for the Urban model
    Double_t fUrbanI;     
    Double_t fUrbanE1;    
    Double_t fUrbanE2;    
    Double_t fUrbanF1;    
    Double_t fUrbanF2;    
    Double_t fUrbanEmax;  
    Double_t fUrbanR;     
 
    // --- 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; 
 
    Double_t InterpolateDataTable(Double_t eKin, std::map<Double_t, Double_t>& table);
 
    void ReadDataTablesStoppingPower();
 
    void ReadDataTablesLandauWidth();
 
 
    void SetUrbanParameters(Double_t z);
 
 
    ClassDef(BmnSsdPhysics,1);
 
};
 
#endif /* BMNSSDPHYSICS_H */