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 #ifndef G4eeToHadronsModel_h
00052 #define G4eeToHadronsModel_h 1
00053
00054 #include "G4VEmModel.hh"
00055
00056 class G4PhysicsVector;
00057 class G4Vee2hadrons;
00058
00059 class G4eeToHadronsModel : public G4VEmModel
00060 {
00061
00062 public:
00063
00064 G4eeToHadronsModel(G4Vee2hadrons*, G4int ver=0,
00065 const G4String& nam = "eeToHadrons");
00066
00067 virtual ~G4eeToHadronsModel();
00068
00069 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
00070
00071 virtual G4double CrossSectionPerVolume(const G4Material*,
00072 const G4ParticleDefinition*,
00073 G4double kineticEnergy,
00074 G4double cutEnergy,
00075 G4double maxEnergy);
00076
00077 virtual G4double ComputeCrossSectionPerAtom(
00078 const G4ParticleDefinition*,
00079 G4double kineticEnergy,
00080 G4double Z, G4double A,
00081 G4double cutEnergy = 0.0,
00082 G4double maxEnergy = DBL_MAX);
00083
00084 virtual G4double ComputeCrossSectionPerElectron(
00085 const G4ParticleDefinition*,
00086 G4double kineticEnergy,
00087 G4double cutEnergy = 0.0,
00088 G4double maxEnergy = DBL_MAX);
00089
00090 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
00091 const G4MaterialCutsCouple*,
00092 const G4DynamicParticle*,
00093 G4double tmin = 0.0,
00094 G4double maxEnergy = DBL_MAX);
00095
00096 G4DynamicParticle* GenerateCMPhoton(G4double);
00097
00098 inline G4double PeakEnergy() const;
00099
00100 private:
00101
00102 void ComputeCMCrossSectionPerElectron();
00103
00104
00105 G4eeToHadronsModel & operator=(const G4eeToHadronsModel &right);
00106 G4eeToHadronsModel(const G4eeToHadronsModel&);
00107
00108 G4Vee2hadrons* model;
00109 G4ParticleDefinition* theGamma;
00110 G4PhysicsVector* crossPerElectron;
00111 G4PhysicsVector* crossBornPerElectron;
00112 G4bool isInitialised;
00113 G4int nbins;
00114 G4int verbose;
00115
00116 G4double lowKinEnergy;
00117 G4double peakKinEnergy;
00118 G4double highKinEnergy;
00119
00120 G4double emin;
00121 G4double epeak;
00122 G4double emax;
00123 };
00124
00125
00126
00127 inline G4double G4eeToHadronsModel::PeakEnergy() const
00128 {
00129 return peakKinEnergy;
00130 }
00131
00132
00133
00134 #endif