| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 // -------------------------------------------
28 //
29 // GEANT4 Class file
30 //
31 //
32 // File name: G4hBremsstrahlungModel
33 //
34 // Author: Vladimir Ivanchenko on base
35 //
36 // Creation date: 28.02.2008
37 //
38 // Modifications:
39 //
40
41 //
42 // Class Description:
43 //
44 //
45 // -------------------------------------------
46 //
47 //....oooOO0OOooo........oooOO0OOooo........oo
48 //....oooOO0OOooo........oooOO0OOooo........oo
49
50 #include "G4hBremsstrahlungModel.hh"
51 #include "G4PhysicalConstants.hh"
52 #include "G4SystemOfUnits.hh"
53 #include "G4Log.hh"
54 #include "G4NistManager.hh"
55
56 using namespace std;
57
58 G4hBremsstrahlungModel::G4hBremsstrahlungModel
59 const G4String& nam)
60 : G4MuBremsstrahlungModel(p, nam)
61 {}
62
63 //....oooOO0OOooo........oooOO0OOooo........oo
64
65 G4hBremsstrahlungModel::~G4hBremsstrahlungMode
66 {}
67
68 //....oooOO0OOooo........oooOO0OOooo........oo
69
70 G4double G4hBremsstrahlungModel::ComputeDMicro
71 G4d
72 G4d
73 G4d
74 // differential cross section
75 {
76 G4double dxsection = 0.;
77
78 if(gammaEnergy > tkin) return dxsection;
79 // G4cout << "G4hBremsstrahlungModel m= " <
80 // << " " << particle->GetParticleName()
81 G4double E = tkin + mass ;
82 G4double v = gammaEnergy/E ;
83 G4double delta = 0.5*mass*mass*v/(E-gammaEne
84 G4double rab0=delta*sqrte ;
85
86 G4int iz = std::max(G4lrint(Z), 1);
87
88 G4double z13 = 1.0/nist->GetZ13(iz);
89 G4double dn = mass*nist->GetA27(iz)/(70.*Me
90
91 G4double b = (1 == iz) ? bh : btf;
92
93 // nucleus contribution logarithm
94 G4double rab1=b*z13;
95 G4double fn=G4Log(rab1/(dn*(electron_mass_c2
96 (mass+delta*(dn*sqrte-2.))) ;
97 fn = std::max(fn, 0.0);
98
99 G4double x = 1.0 - v;
100 if(particle->GetPDGSpin() != 0) { x += 0.75*
101
102 dxsection = coeff*x*Z*Z*fn/gammaEnergy;
103
104 return dxsection;
105 }
106
107 //....oooOO0OOooo........oooOO0OOooo........oo
108