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 G4mplIonisationWithDeltaModel_h
00050 #define G4mplIonisationWithDeltaModel_h 1
00051
00052 #include "G4VEmModel.hh"
00053 #include "G4VEmFluctuationModel.hh"
00054
00055 class G4ParticleChangeForLoss;
00056
00057 class G4mplIonisationWithDeltaModel : public G4VEmModel, public G4VEmFluctuationModel
00058 {
00059
00060 public:
00061
00062 G4mplIonisationWithDeltaModel(G4double mCharge,
00063 const G4String& nam = "mplIonisationWithDelta");
00064
00065 virtual ~G4mplIonisationWithDeltaModel();
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 G4double ComputeCrossSectionPerElectron(
00075 const G4ParticleDefinition*,
00076 G4double kineticEnergy,
00077 G4double cutEnergy,
00078 G4double maxEnergy);
00079
00080 virtual G4double ComputeCrossSectionPerAtom(
00081 const G4ParticleDefinition*,
00082 G4double kineticEnergy,
00083 G4double Z, G4double A,
00084 G4double cutEnergy,
00085 G4double maxEnergy);
00086
00087 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00088 const G4MaterialCutsCouple*,
00089 const G4DynamicParticle*,
00090 G4double tmin,
00091 G4double maxEnergy);
00092
00093
00094 virtual G4double SampleFluctuations(const G4Material*,
00095 const G4DynamicParticle*,
00096 G4double& tmax,
00097 G4double& length,
00098 G4double& meanLoss);
00099
00100 virtual G4double Dispersion(const G4Material*,
00101 const G4DynamicParticle*,
00102 G4double& tmax,
00103 G4double& length);
00104
00105 void SetParticle(const G4ParticleDefinition* p);
00106
00107 protected:
00108
00109 virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
00110 G4double kinEnergy);
00111
00112 private:
00113
00114 G4double ComputeDEDXAhlen(const G4Material* material, G4double bg2, G4double cut);
00115
00116
00117 G4mplIonisationWithDeltaModel & operator=(const G4mplIonisationWithDeltaModel &right);
00118 G4mplIonisationWithDeltaModel(const G4mplIonisationWithDeltaModel&);
00119
00120 const G4ParticleDefinition* monopole;
00121 G4ParticleDefinition* theElectron;
00122 G4ParticleChangeForLoss* fParticleChange;
00123
00124 G4double mass;
00125 G4double magCharge;
00126 G4double twoln10;
00127 G4double betalow;
00128 G4double betalim;
00129 G4double beta2lim;
00130 G4double bg2lim;
00131 G4double chargeSquare;
00132 G4double dedxlim;
00133 G4int nmpl;
00134 G4double pi_hbarc2_over_mc2;
00135 };
00136
00137 #endif
00138
00139