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
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056 #ifndef G4MollerBhabhaModel_h
00057 #define G4MollerBhabhaModel_h 1
00058
00059 #include "G4VEmModel.hh"
00060
00061 class G4ParticleChangeForLoss;
00062
00063 class G4MollerBhabhaModel : public G4VEmModel
00064 {
00065
00066 public:
00067
00068 G4MollerBhabhaModel(const G4ParticleDefinition* p = 0,
00069 const G4String& nam = "MollerBhabha");
00070
00071 virtual ~G4MollerBhabhaModel();
00072
00073 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00074
00075 virtual G4double ComputeCrossSectionPerElectron(
00076 const G4ParticleDefinition*,
00077 G4double kineticEnergy,
00078 G4double cutEnergy,
00079 G4double maxEnergy);
00080
00081 virtual G4double ComputeCrossSectionPerAtom(
00082 const G4ParticleDefinition*,
00083 G4double kineticEnergy,
00084 G4double Z, G4double A,
00085 G4double cutEnergy,
00086 G4double maxEnergy);
00087
00088 virtual G4double CrossSectionPerVolume(const G4Material*,
00089 const G4ParticleDefinition*,
00090 G4double kineticEnergy,
00091 G4double cutEnergy,
00092 G4double maxEnergy);
00093
00094 virtual G4double ComputeDEDXPerVolume(const G4Material*,
00095 const G4ParticleDefinition*,
00096 G4double kineticEnergy,
00097 G4double cutEnergy);
00098
00099 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00100 const G4MaterialCutsCouple*,
00101 const G4DynamicParticle*,
00102 G4double tmin,
00103 G4double maxEnergy);
00104
00105 protected:
00106
00107 virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
00108 G4double kinEnergy);
00109
00110 inline void SetParticle(const G4ParticleDefinition* p);
00111
00112 const G4ParticleDefinition* particle;
00113 G4ParticleDefinition* theElectron;
00114 G4ParticleChangeForLoss* fParticleChange;
00115
00116 G4bool isElectron;
00117 G4double twoln10;
00118 G4double lowLimit;
00119
00120 private:
00121
00122
00123 G4MollerBhabhaModel & operator=(const G4MollerBhabhaModel &right);
00124 G4MollerBhabhaModel(const G4MollerBhabhaModel&);
00125
00126 G4bool isInitialised;
00127
00128 };
00129
00130
00131
00132 inline void G4MollerBhabhaModel::SetParticle(const G4ParticleDefinition* p)
00133 {
00134 particle = p;
00135 if(p != theElectron) { isElectron = false; }
00136 }
00137
00138
00139
00140 #endif