00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #ifndef G4LivermoreRayleighModel_h
00032 #define G4LivermoreRayleighModel_h 1
00033
00034 #include "G4VEmModel.hh"
00035 #include "G4ParticleChangeForGamma.hh"
00036 #include "G4LPhysicsFreeVector.hh"
00037 #include "G4ProductionCutsTable.hh"
00038
00039 class G4LivermoreRayleighModel : public G4VEmModel
00040 {
00041
00042 public:
00043
00044 G4LivermoreRayleighModel();
00045
00046 virtual ~G4LivermoreRayleighModel();
00047
00048 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00049
00050 virtual G4double ComputeCrossSectionPerAtom(
00051 const G4ParticleDefinition*,
00052 G4double kinEnergy,
00053 G4double Z,
00054 G4double A=0,
00055 G4double cut=0,
00056 G4double emax=DBL_MAX);
00057
00058 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00059 const G4MaterialCutsCouple*,
00060 const G4DynamicParticle*,
00061 G4double tmin,
00062 G4double maxEnergy);
00063
00064 inline void SetLowEnergyThreshold(G4double);
00065
00066 private:
00067
00068 void ReadData(size_t Z, const char* path = 0);
00069
00070 G4LivermoreRayleighModel & operator=(const G4LivermoreRayleighModel &right);
00071 G4LivermoreRayleighModel(const G4LivermoreRayleighModel&);
00072
00073 G4bool isInitialised;
00074 G4int maxZ;
00075 G4int verboseLevel;
00076
00077 G4double lowEnergyLimit;
00078
00079 std::vector<G4LPhysicsFreeVector*> dataCS;
00080
00081 G4ParticleChangeForGamma* fParticleChange;
00082
00083 };
00084
00085 inline void G4LivermoreRayleighModel::SetLowEnergyThreshold(G4double val)
00086 {
00087 lowEnergyLimit = val;
00088 }
00089
00090
00091
00092 #endif