| 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: G4Generator2BS
33 //
34 // Author: Andreia Trindade (andreia@lip.p
35 // Pedro Rodrigues (psilva@lip.pt
36 // Luis Peralta (luis@lip.pt)
37 //
38 // Creation date: 2 June 2003
39 //
40 // Modifications:
41 // 02 Jun 2003 First implementation acording
42 // 12 Oct 2010 V.Ivanchenko moved RejectionFu
43 //
44 //
45 // Class Description:
46 //
47 // Concrete class for Bremsstrahlung Angular D
48 // 2BS Distribution
49 //
50
51 // -------------------------------------------
52 //
53
54 #ifndef G4Generator2BS_h
55 #define G4Generator2BS_h 1
56
57 #include "G4ios.hh"
58 #include "globals.hh"
59 #include "G4VEmAngularDistribution.hh"
60 #include "G4Log.hh"
61
62 class G4Pow;
63
64 class G4Generator2BS : public G4VEmAngularDist
65 {
66
67 public:
68
69 explicit G4Generator2BS(const G4String& name
70
71 virtual ~G4Generator2BS();
72
73 G4ThreeVector& SampleDirection(const G4Dynam
74 G4double out_energy,
75 G4int Z,
76 const G4Material* mat = nullptr) over
77
78 void PrintGeneratorInformation() const overr
79
80 G4Generator2BS & operator=(const G4Generato
81 G4Generator2BS(const G4Generator2BS&) = del
82
83 protected:
84 inline G4double RejectionFunction(G4double v
85
86 private:
87 G4Pow* g4pow;
88
89 G4double fz;
90 G4double ratio;
91 G4double ratio1;
92 G4double ratio2;
93 G4double delta;
94
95 G4int nwarn;
96
97 };
98
99 inline G4double G4Generator2BS::RejectionFunct
100 {
101 G4double y2 = (1 + y)*(1 + y);
102 G4double x = 4*y*ratio/y2;
103 return 4*x - ratio1 - (ratio2 - x)*G4Log(del
104 }
105
106 #endif
107
108