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 #ifndef G4mplIonisationModel_h
00050 #define G4mplIonisationModel_h 1
00051
00052 #include "G4VEmModel.hh"
00053 #include "G4VEmFluctuationModel.hh"
00054
00055 class G4ParticleChangeForLoss;
00056
00057 class G4mplIonisationModel : public G4VEmModel, public G4VEmFluctuationModel
00058 {
00059
00060 public:
00061
00062 G4mplIonisationModel(G4double mCharge,
00063 const G4String& nam = "mplIonisation");
00064
00065 virtual ~G4mplIonisationModel();
00066
00067 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00068
00069 virtual G4double ComputeDEDXPerVolume(const G4Material*,
00070 const G4ParticleDefinition*,
00071 G4double kineticEnergy,
00072 G4double cutEnergy);
00073
00074 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00075 const G4MaterialCutsCouple*,
00076 const G4DynamicParticle*,
00077 G4double tmin,
00078 G4double maxEnergy);
00079
00080
00081 virtual G4double SampleFluctuations(const G4Material*,
00082 const G4DynamicParticle*,
00083 G4double& tmax,
00084 G4double& length,
00085 G4double& meanLoss);
00086
00087 virtual G4double Dispersion(const G4Material*,
00088 const G4DynamicParticle*,
00089 G4double& tmax,
00090 G4double& length);
00091
00092 void SetParticle(const G4ParticleDefinition* p);
00093
00094 private:
00095
00096 G4double ComputeDEDXAhlen(const G4Material* material, G4double bg2);
00097
00098
00099 G4mplIonisationModel & operator=(const G4mplIonisationModel &right);
00100 G4mplIonisationModel(const G4mplIonisationModel&);
00101
00102 const G4ParticleDefinition* monopole;
00103 G4ParticleChangeForLoss* fParticleChange;
00104
00105 G4double mass;
00106 G4double magCharge;
00107 G4double twoln10;
00108 G4double betalow;
00109 G4double betalim;
00110 G4double beta2lim;
00111 G4double bg2lim;
00112 G4double chargeSquare;
00113 G4double dedxlim;
00114 G4int nmpl;
00115 G4double pi_hbarc2_over_mc2;
00116 };
00117
00118 #endif
00119
00120