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 G4LivermoreNuclearGammaConversionModel_h
00032 #define G4LivermoreNuclearGammaConversionModel_h 1
00033
00034 #include "G4VEmModel.hh"
00035 #include "G4Electron.hh"
00036 #include "G4Positron.hh"
00037 #include "G4ParticleChangeForGamma.hh"
00038 #include "G4CrossSectionHandler.hh"
00039 #include "G4ForceCondition.hh"
00040 #include "G4CompositeEMDataSet.hh"
00041 #include "G4Gamma.hh"
00042
00043 class G4LivermoreNuclearGammaConversionModel : public G4VEmModel
00044 {
00045
00046 public:
00047
00048 G4LivermoreNuclearGammaConversionModel(const G4ParticleDefinition* p = 0,
00049 const G4String& nam = "LivermoreNuclearGammaConversion");
00050
00051 virtual ~G4LivermoreNuclearGammaConversionModel();
00052
00053 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00054
00055 virtual G4double ComputeCrossSectionPerAtom(
00056 const G4ParticleDefinition*,
00057 G4double kinEnergy,
00058 G4double Z,
00059 G4double A=0,
00060 G4double cut=0,
00061 G4double emax=DBL_MAX);
00062
00063 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00064 const G4MaterialCutsCouple*,
00065 const G4DynamicParticle*,
00066 G4double tmin,
00067 G4double maxEnergy);
00068
00069 protected:
00070
00071 G4ParticleChangeForGamma* fParticleChange;
00072
00073 G4double GetMeanFreePath(const G4Track& aTrack,
00074 G4double previousStepSize,
00075 G4ForceCondition* condition);
00076 private:
00077
00078 const G4double smallEnergy;
00079 G4double lowEnergyLimit;
00080 G4double highEnergyLimit;
00081
00082 G4int verboseLevel;
00083 G4bool isInitialised;
00084
00085 G4double ScreenFunction1(G4double screenVariable);
00086 G4double ScreenFunction2(G4double screenVariable);
00087 G4VCrossSectionHandler* crossSectionHandler;
00088 G4VEMDataSet* meanFreePathTable;
00089
00090 G4LivermoreNuclearGammaConversionModel & operator=(const G4LivermoreNuclearGammaConversionModel &right);
00091 G4LivermoreNuclearGammaConversionModel(const G4LivermoreNuclearGammaConversionModel&);
00092
00093
00094 };
00095
00096
00097
00098 #endif