| 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 #include "G4GammaXTRadiator.hh"
28
29 #include "G4Gamma.hh"
30
31 //////////////////////////////////////////////
32 // Constructor, destructor
33 G4GammaXTRadiator::G4GammaXTRadiator(G4Logical
34 G4double
35 G4Materia
36 G4double
37 const G4S
38 : G4VXTRenergyLoss(anEnvelope, foilMat, gasM
39 {
40 G4cout << "Gamma distributed X-ray TR radiat
41
42 // Build energy and angular integral spectra
43 // a radiator
44
45 fAlphaPlate = alphaPlate;
46 fAlphaGas = alphaGas;
47 G4cout << "fAlphaPlate = " << fAlphaPlate <<
48 << G4endl;
49 }
50
51 //////////////////////////////////////////////
52 G4GammaXTRadiator::~G4GammaXTRadiator() = defa
53
54 void G4GammaXTRadiator::ProcessDescription(std
55 {
56 out
57 << "Rough approximation describing a radia
58 "radiation.\n"
59 "Thicknesses of plates and gas gaps are
60 "description.\n";
61 }
62
63 //////////////////////////////////////////////
64 // Rough approximation for radiator interferen
65 // fully GamDistr radiator. The plate and gas
66 // according to exponent. The mean values of t
67 // are supposed to be about XTR formation zone
68 // mean absorption length of XTR photons in co
69 G4double G4GammaXTRadiator::GetStackFactor(G4d
70 G4d
71 {
72 G4double result, Za, Zb, Ma, Mb;
73
74 Za = GetPlateFormationZone(energy, gamma, va
75 Zb = GetGasFormationZone(energy, gamma, varA
76
77 Ma = GetPlateLinearPhotoAbs(energy);
78 Mb = GetGasLinearPhotoAbs(energy);
79
80 G4complex Ca(1.0 + 0.5 * fPlateThick * Ma /
81 fPlateThick / Za / fAlphaPlate)
82 G4complex Cb(1.0 + 0.5 * fGasThick * Mb / fA
83 fGasThick / Zb / fAlphaGas);
84
85 G4complex Ha = std::pow(Ca, -fAlphaPlate);
86 G4complex Hb = std::pow(Cb, -fAlphaGas);
87 G4complex H = Ha * Hb;
88
89 G4complex F1 = (1.0 - Ha) * (1.0 - Hb) / (1.
90
91 G4complex F2 = (1.0 - Ha) * (1.0 - Ha) * Hb
92 (1.0 - std::pow(H, fPlateNumb
93
94 G4complex R = (F1 + F2) * OneInterfaceXTRdEd
95
96 result = 2.0 * std::real(R);
97
98 return result;
99 }
100