File indexing completed on 2025-01-18 09:58:38
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031 #ifndef G4LivermoreRayleighModel_h
0032 #define G4LivermoreRayleighModel_h 1
0033
0034 #include "G4ParticleChangeForGamma.hh"
0035 #include "G4PhysicsFreeVector.hh"
0036 #include "G4ProductionCutsTable.hh"
0037 #include "G4VEmModel.hh"
0038
0039 class G4LivermoreRayleighModel : public G4VEmModel
0040 {
0041 public:
0042 explicit G4LivermoreRayleighModel();
0043 ~G4LivermoreRayleighModel() override;
0044
0045 void Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
0046 void InitialiseLocal(const G4ParticleDefinition*, G4VEmModel* masterModel) override;
0047 void InitialiseForElement(const G4ParticleDefinition*, G4int Z) override;
0048
0049 G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*, G4double kinEnergy, G4double Z,
0050 G4double A = 0, G4double cut = 0,
0051 G4double emax = DBL_MAX) override;
0052
0053 void SampleSecondaries(std::vector<G4DynamicParticle*>*, const G4MaterialCutsCouple*,
0054 const G4DynamicParticle*, G4double tmin, G4double maxEnergy) override;
0055
0056 [[maybe_unused]] inline void SetLowEnergyThreshold(G4double);
0057
0058 G4LivermoreRayleighModel& operator=(const G4LivermoreRayleighModel& right) = delete;
0059 G4LivermoreRayleighModel(const G4LivermoreRayleighModel&) = delete;
0060
0061 private:
0062 void ReadData(const G4int ZZ);
0063 const G4String& FindDirectoryPath();
0064
0065 G4ParticleChangeForGamma* fParticleChange;
0066
0067 static G4PhysicsFreeVector* dataCS[101];
0068 static G4String gDataDirectory;
0069
0070 G4double lowEnergyLimit;
0071 G4int verboseLevel;
0072 G4int maxZ = 100;
0073 G4bool isInitialised = false;
0074 };
0075
0076 [[maybe_unused]] inline void G4LivermoreRayleighModel::SetLowEnergyThreshold(G4double val)
0077 {
0078 lowEnergyLimit = val;
0079 }
0080
0081
0082 #endif