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
00032
00033 #ifndef G4LivermorePolarizedRayleighModel_h
00034 #define G4LivermorePolarizedRayleighModel_h 1
00035
00036 #include "G4VEmModel.hh"
00037 #include "G4ParticleChangeForGamma.hh"
00038 #include "G4CrossSectionHandler.hh"
00039 #include "G4LogLogInterpolation.hh"
00040 #include "G4CompositeEMDataSet.hh"
00041 #include "G4Gamma.hh"
00042
00043 class G4LivermorePolarizedRayleighModel : public G4VEmModel
00044 {
00045
00046 public:
00047
00048 G4LivermorePolarizedRayleighModel(const G4ParticleDefinition* p = 0,
00049 const G4String& nam = "LivermorePolarizedRayleigh");
00050
00051 virtual ~G4LivermorePolarizedRayleighModel();
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 private:
00074
00075 G4double lowEnergyLimit;
00076 G4double highEnergyLimit;
00077
00078 G4int verboseLevel;
00079 G4bool isInitialised;
00080
00081 G4VCrossSectionHandler* crossSectionHandler;
00082 G4VEMDataSet* formFactorData;
00083
00084
00085
00086
00087
00088 G4double GenerateCosTheta(G4double incomingPhotonEnergy, G4int zAtom) const;
00089
00090
00091
00092
00093 G4double GeneratePhi(G4double cosTheta) const;
00094
00095
00096
00097 G4double GeneratePolarizationAngle(void) const;
00098
00099 G4ThreeVector GetPhotonPolarization(const G4DynamicParticle& photon);
00100
00101 G4LivermorePolarizedRayleighModel & operator=(const G4LivermorePolarizedRayleighModel &right);
00102 G4LivermorePolarizedRayleighModel(const G4LivermorePolarizedRayleighModel&);
00103
00104 };
00105
00106
00107
00108 #endif