| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
>> 26 // $Id: G4eeToHadronsModel.hh,v 1.6 2007/05/22 17:37:30 vnivanch Exp $
>> 27 // GEANT4 tag $Name: geant4-09-01-patch-03 $
26 // 28 //
27 // ------------------------------------------- 29 // -------------------------------------------------------------------
28 // 30 //
29 // GEANT4 Class header file 31 // GEANT4 Class header file
30 // 32 //
31 // 33 //
32 // File name: G4eeToHadronsModel 34 // File name: G4eeToHadronsModel
33 // 35 //
34 // Author: Vladimir Ivanchenko 36 // Author: Vladimir Ivanchenko
35 // 37 //
36 // Creation date: 25.10.2003 38 // Creation date: 25.10.2003
37 // 39 //
38 // Modifications: 40 // Modifications:
39 // 08-04-05 Major optimisation of internal int 41 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
40 // 18-05-05 Use optimized interfaces (V.Ivantc 42 // 18-05-05 Use optimized interfaces (V.Ivantchenko)
41 // 43 //
42 44
43 // 45 //
44 // Class Description: 46 // Class Description:
45 // 47 //
46 48
47 // ------------------------------------------- 49 // -------------------------------------------------------------------
48 // 50 //
49 51
50 #ifndef G4eeToHadronsModel_h 52 #ifndef G4eeToHadronsModel_h
51 #define G4eeToHadronsModel_h 1 53 #define G4eeToHadronsModel_h 1
52 54
53 #include "G4VEmModel.hh" 55 #include "G4VEmModel.hh"
54 56
55 class G4PhysicsVector; 57 class G4PhysicsVector;
56 class G4Vee2hadrons; 58 class G4Vee2hadrons;
57 59
58 class G4eeToHadronsModel : public G4VEmModel 60 class G4eeToHadronsModel : public G4VEmModel
59 { 61 {
60 62
61 public: 63 public:
62 64
63 explicit G4eeToHadronsModel(G4Vee2hadrons*, << 65 G4eeToHadronsModel(const G4Vee2hadrons*, G4int ver=0,
64 const G4String& << 66 const G4String& nam = "eeToHadrons");
65 67
66 ~G4eeToHadronsModel() override; << 68 virtual ~G4eeToHadronsModel();
67 69
68 void Initialise(const G4ParticleDefinition*, << 70 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
69 71
70 G4double CrossSectionPerVolume(const G4Mater << 72 G4double PeakEnergy() const;
71 const G4ParticleDefinition*, <<
72 G4double kineticEnergy, <<
73 G4double cutEnergy, <<
74 G4double maxEnergy) override; <<
75 <<
76 G4double ComputeCrossSectionPerAtom(const G4 <<
77 G4double <<
78 G4double <<
79 G4double <<
80 G4double <<
81 <<
82 G4double ComputeCrossSectionPerElectron(cons <<
83 G4do <<
84 G4do <<
85 G4do <<
86 <<
87 void SampleSecondaries(std::vector<G4Dynamic <<
88 const G4MaterialCutsCouple*, <<
89 const G4DynamicParticle*, <<
90 G4double tmin = 0.0, <<
91 G4double maxEnergy = DBL_MAX) override; <<
92 73
93 G4DynamicParticle* GenerateCMPhoton(G4double << 74 virtual G4double CrossSectionPerVolume(const G4Material*,
94 << 75 const G4ParticleDefinition*,
95 inline G4double PeakEnergy() const; << 76 G4double kineticEnergy,
>> 77 G4double cutEnergy,
>> 78 G4double maxEnergy);
>> 79
>> 80 virtual G4double ComputeCrossSectionPerAtom(
>> 81 const G4ParticleDefinition*,
>> 82 G4double kineticEnergy,
>> 83 G4double Z, G4double A,
>> 84 G4double cutEnergy = 0.0,
>> 85 G4double maxEnergy = DBL_MAX);
>> 86
>> 87 virtual G4double ComputeCrossSectionPerElectron(
>> 88 const G4ParticleDefinition*,
>> 89 G4double kineticEnergy,
>> 90 G4double cutEnergy = 0.0,
>> 91 G4double maxEnergy = DBL_MAX);
>> 92
>> 93 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
>> 94 const G4MaterialCutsCouple*,
>> 95 const G4DynamicParticle*,
>> 96 G4double tmin = 0.0,
>> 97 G4double maxEnergy = DBL_MAX);
96 98
97 // hide assignment operator << 99 G4DynamicParticle* GenerateCMPhoton(G4double);
98 G4eeToHadronsModel & operator=(const G4eeTo <<
99 G4eeToHadronsModel(const G4eeToHadronsModel <<
100 100
101 private: 101 private:
102 102
103 void ComputeCMCrossSectionPerElectron(); 103 void ComputeCMCrossSectionPerElectron();
104 104
105 const G4ParticleDefinition* theGamma; << 105 // hide assignment operator
>> 106 G4eeToHadronsModel & operator=(const G4eeToHadronsModel &right);
>> 107 G4eeToHadronsModel(const G4eeToHadronsModel&);
106 108
107 G4Vee2hadrons* model; << 109 const G4Vee2hadrons* model;
108 G4PhysicsVector* crossPerElectron = nul << 110 G4ParticleDefinition* theGamma;
109 G4PhysicsVector* crossBornPerElectron = << 111 G4PhysicsVector* crossPerElectron;
>> 112 G4PhysicsVector* crossBornPerElectron;
>> 113 G4bool isInitialised;
>> 114 G4int nbins;
>> 115 G4int verbose;
110 116
111 G4double lowKinEnergy; 117 G4double lowKinEnergy;
112 G4double peakKinEnergy; 118 G4double peakKinEnergy;
113 G4double highKinEnergy; 119 G4double highKinEnergy;
114 120
115 G4double emin; 121 G4double emin;
116 G4double epeak; 122 G4double epeak;
117 G4double emax; 123 G4double emax;
118 <<
119 G4int nbins = 100; <<
120 G4int verbose; <<
121 G4bool isInitialised = false; <<
122 }; 124 };
123 125
124 //....oooOO0OOooo........oooOO0OOooo........oo 126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
125 127
126 inline G4double G4eeToHadronsModel::PeakEnergy 128 inline G4double G4eeToHadronsModel::PeakEnergy() const
127 { 129 {
128 return peakKinEnergy; 130 return peakKinEnergy;
>> 131 }
>> 132
>> 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 134
>> 135 inline G4double G4eeToHadronsModel::CrossSectionPerVolume(
>> 136 const G4Material* mat,
>> 137 const G4ParticleDefinition* p,
>> 138 G4double kineticEnergy,
>> 139 G4double, G4double)
>> 140 {
>> 141 return mat->GetElectronDensity()*
>> 142 ComputeCrossSectionPerElectron(p, kineticEnergy);
>> 143 }
>> 144
>> 145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 146
>> 147 inline G4double G4eeToHadronsModel::ComputeCrossSectionPerAtom(
>> 148 const G4ParticleDefinition* p,
>> 149 G4double kineticEnergy,
>> 150 G4double Z, G4double,
>> 151 G4double, G4double)
>> 152 {
>> 153 return Z*ComputeCrossSectionPerElectron(p, kineticEnergy);
129 } 154 }
130 155
131 //....oooOO0OOooo........oooOO0OOooo........oo 156 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
132 157
133 #endif 158 #endif
134 159