1 | // |
---|
2 | // ******************************************************************** |
---|
3 | // * License and Disclaimer * |
---|
4 | // * * |
---|
5 | // * The Geant4 software is copyright of the Copyright Holders of * |
---|
6 | // * the Geant4 Collaboration. It is provided under the terms and * |
---|
7 | // * conditions of the Geant4 Software License, included in the file * |
---|
8 | // * LICENSE and available at http://cern.ch/geant4/license . These * |
---|
9 | // * include a list of copyright holders. * |
---|
10 | // * * |
---|
11 | // * Neither the authors of this software system, nor their employing * |
---|
12 | // * institutes,nor the agencies providing financial support for this * |
---|
13 | // * work make any representation or warranty, express or implied, * |
---|
14 | // * regarding this software system or assume any liability for its * |
---|
15 | // * use. Please see the license in the file LICENSE and URL above * |
---|
16 | // * for the full disclaimer and the limitation of liability. * |
---|
17 | // * * |
---|
18 | // * This code implementation is the result of the scientific and * |
---|
19 | // * technical work of the GEANT4 collaboration. * |
---|
20 | // * By using, copying, modifying or distributing the software (or * |
---|
21 | // * any work based on the software) you agree to acknowledge its * |
---|
22 | // * use in resulting scientific publications, and indicate your * |
---|
23 | // * acceptance of all terms of the Geant4 Software license. * |
---|
24 | // ******************************************************************** |
---|
25 | // |
---|
26 | // |
---|
27 | // $Id: G4TransitionRadiation.cc,v 1.7 2006/06/29 19:56:19 gunter Exp $ |
---|
28 | // GEANT4 tag $Name: geant4-09-02-ref-02 $ |
---|
29 | // |
---|
30 | // G4TransitionRadiation class -- implementation file |
---|
31 | |
---|
32 | // GEANT 4 class implementation file --- Copyright CERN 1995 |
---|
33 | // CERN Geneva Switzerland |
---|
34 | |
---|
35 | // For information related to this code, please, contact |
---|
36 | // CERN, CN Division, ASD Group |
---|
37 | // History: |
---|
38 | // 1st version 11.09.97 V. Grichine (Vladimir.Grichine@cern.ch ) |
---|
39 | // 2nd version 16.12.97 V. Grichine |
---|
40 | // 3rd version 28.07.05, P.Gumplinger add G4ProcessType to constructor |
---|
41 | |
---|
42 | |
---|
43 | #include <cmath> |
---|
44 | // #include "G4ios.hh" |
---|
45 | // #include <fstream.h> |
---|
46 | // #include <stdlib.h> |
---|
47 | |
---|
48 | #include "G4TransitionRadiation.hh" |
---|
49 | #include "G4Material.hh" |
---|
50 | |
---|
51 | // Init gamma array |
---|
52 | |
---|
53 | |
---|
54 | // Local constants |
---|
55 | |
---|
56 | const G4int G4TransitionRadiation::fSympsonNumber = 100 ; |
---|
57 | const G4int G4TransitionRadiation::fGammaNumber = 15 ; |
---|
58 | const G4int G4TransitionRadiation::fPointNumber = 100 ; |
---|
59 | |
---|
60 | using namespace std; |
---|
61 | |
---|
62 | /////////////////////////////////////////////////////////////////////// |
---|
63 | // |
---|
64 | // Constructor for selected couple of materials |
---|
65 | // |
---|
66 | |
---|
67 | G4TransitionRadiation:: |
---|
68 | G4TransitionRadiation( const G4String& processName, G4ProcessType type ) |
---|
69 | : G4VDiscreteProcess(processName, type) |
---|
70 | { |
---|
71 | // fMatIndex1 = pMat1->GetIndex() ; |
---|
72 | // fMatIndex2 = pMat2->GetIndex() ; |
---|
73 | } |
---|
74 | |
---|
75 | ////////////////////////////////////////////////////////////////////// |
---|
76 | // |
---|
77 | // Destructor |
---|
78 | // |
---|
79 | |
---|
80 | G4TransitionRadiation::~G4TransitionRadiation() |
---|
81 | { |
---|
82 | ; |
---|
83 | } |
---|
84 | |
---|
85 | |
---|
86 | /////////////////////////////////////////////////////////////////// |
---|
87 | // |
---|
88 | // Sympson integral of TR spectral-angle density over energy between |
---|
89 | // the limits energy 1 and energy2 at fixed varAngle = 1 - cos(Theta) |
---|
90 | |
---|
91 | G4double |
---|
92 | G4TransitionRadiation::IntegralOverEnergy( G4double energy1, |
---|
93 | G4double energy2, |
---|
94 | G4double varAngle ) const |
---|
95 | { |
---|
96 | G4int i ; |
---|
97 | G4double h , sumEven = 0.0 , sumOdd = 0.0 ; |
---|
98 | h = 0.5*(energy2 - energy1)/fSympsonNumber ; |
---|
99 | for(i=1;i<fSympsonNumber;i++) |
---|
100 | { |
---|
101 | sumEven += SpectralAngleTRdensity(energy1 + 2*i*h,varAngle) ; |
---|
102 | sumOdd += SpectralAngleTRdensity(energy1 + (2*i - 1)*h,varAngle) ; |
---|
103 | } |
---|
104 | sumOdd += SpectralAngleTRdensity(energy1 + (2*fSympsonNumber - 1)*h,varAngle) ; |
---|
105 | return h*( SpectralAngleTRdensity(energy1,varAngle) |
---|
106 | + SpectralAngleTRdensity(energy2,varAngle) |
---|
107 | + 4.0*sumOdd + 2.0*sumEven )/3.0 ; |
---|
108 | } |
---|
109 | |
---|
110 | |
---|
111 | |
---|
112 | /////////////////////////////////////////////////////////////////// |
---|
113 | // |
---|
114 | // Sympson integral of TR spectral-angle density over energy between |
---|
115 | // the limits varAngle1 and varAngle2 at fixed energy |
---|
116 | |
---|
117 | G4double |
---|
118 | G4TransitionRadiation::IntegralOverAngle( G4double energy, |
---|
119 | G4double varAngle1, |
---|
120 | G4double varAngle2 ) const |
---|
121 | { |
---|
122 | G4int i ; |
---|
123 | G4double h , sumEven = 0.0 , sumOdd = 0.0 ; |
---|
124 | h = 0.5*(varAngle2 - varAngle1)/fSympsonNumber ; |
---|
125 | for(i=1;i<fSympsonNumber;i++) |
---|
126 | { |
---|
127 | sumEven += SpectralAngleTRdensity(energy,varAngle1 + 2*i*h) ; |
---|
128 | sumOdd += SpectralAngleTRdensity(energy,varAngle1 + (2*i - 1)*h) ; |
---|
129 | } |
---|
130 | sumOdd += SpectralAngleTRdensity(energy,varAngle1 + (2*fSympsonNumber - 1)*h) ; |
---|
131 | |
---|
132 | return h*( SpectralAngleTRdensity(energy,varAngle1) |
---|
133 | + SpectralAngleTRdensity(energy,varAngle2) |
---|
134 | + 4.0*sumOdd + 2.0*sumEven )/3.0 ; |
---|
135 | } |
---|
136 | |
---|
137 | /////////////////////////////////////////////////////////////////// |
---|
138 | // |
---|
139 | // The number of transition radiation photons generated in the |
---|
140 | // angle interval between varAngle1 and varAngle2 |
---|
141 | // |
---|
142 | |
---|
143 | G4double G4TransitionRadiation:: |
---|
144 | AngleIntegralDistribution( G4double varAngle1, |
---|
145 | G4double varAngle2 ) const |
---|
146 | { |
---|
147 | G4int i ; |
---|
148 | G4double h , sumEven = 0.0 , sumOdd = 0.0 ; |
---|
149 | h = 0.5*(varAngle2 - varAngle1)/fSympsonNumber ; |
---|
150 | for(i=1;i<fSympsonNumber;i++) |
---|
151 | { |
---|
152 | sumEven += IntegralOverEnergy(fMinEnergy, |
---|
153 | fMinEnergy +0.3*(fMaxEnergy-fMinEnergy), |
---|
154 | varAngle1 + 2*i*h) |
---|
155 | + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
156 | fMaxEnergy, |
---|
157 | varAngle1 + 2*i*h); |
---|
158 | sumOdd += IntegralOverEnergy(fMinEnergy, |
---|
159 | fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
160 | varAngle1 + (2*i - 1)*h) |
---|
161 | + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
162 | fMaxEnergy, |
---|
163 | varAngle1 + (2*i - 1)*h) ; |
---|
164 | } |
---|
165 | sumOdd += IntegralOverEnergy(fMinEnergy, |
---|
166 | fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
167 | varAngle1 + (2*fSympsonNumber - 1)*h) |
---|
168 | + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
169 | fMaxEnergy, |
---|
170 | varAngle1 + (2*fSympsonNumber - 1)*h) ; |
---|
171 | |
---|
172 | return h*(IntegralOverEnergy(fMinEnergy, |
---|
173 | fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
174 | varAngle1) |
---|
175 | + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
176 | fMaxEnergy, |
---|
177 | varAngle1) |
---|
178 | + IntegralOverEnergy(fMinEnergy, |
---|
179 | fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
180 | varAngle2) |
---|
181 | + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy), |
---|
182 | fMaxEnergy, |
---|
183 | varAngle2) |
---|
184 | + 4.0*sumOdd + 2.0*sumEven )/3.0 ; |
---|
185 | } |
---|
186 | |
---|
187 | /////////////////////////////////////////////////////////////////// |
---|
188 | // |
---|
189 | // The number of transition radiation photons, generated in the |
---|
190 | // energy interval between energy1 and energy2 |
---|
191 | // |
---|
192 | |
---|
193 | G4double G4TransitionRadiation:: |
---|
194 | EnergyIntegralDistribution( G4double energy1, |
---|
195 | G4double energy2 ) const |
---|
196 | { |
---|
197 | G4int i ; |
---|
198 | G4double h , sumEven = 0.0 , sumOdd = 0.0 ; |
---|
199 | h = 0.5*(energy2 - energy1)/fSympsonNumber ; |
---|
200 | for(i=1;i<fSympsonNumber;i++) |
---|
201 | { |
---|
202 | sumEven += IntegralOverAngle(energy1 + 2*i*h,0.0,0.01*fMaxTheta ) |
---|
203 | + IntegralOverAngle(energy1 + 2*i*h,0.01*fMaxTheta,fMaxTheta); |
---|
204 | sumOdd += IntegralOverAngle(energy1 + (2*i - 1)*h,0.0,0.01*fMaxTheta) |
---|
205 | + IntegralOverAngle(energy1 + (2*i - 1)*h,0.01*fMaxTheta,fMaxTheta) ; |
---|
206 | } |
---|
207 | sumOdd += IntegralOverAngle(energy1 + (2*fSympsonNumber - 1)*h, |
---|
208 | 0.0,0.01*fMaxTheta) |
---|
209 | + IntegralOverAngle(energy1 + (2*fSympsonNumber - 1)*h, |
---|
210 | 0.01*fMaxTheta,fMaxTheta) ; |
---|
211 | |
---|
212 | return h*(IntegralOverAngle(energy1,0.0,0.01*fMaxTheta) |
---|
213 | + IntegralOverAngle(energy1,0.01*fMaxTheta,fMaxTheta) |
---|
214 | + IntegralOverAngle(energy2,0.0,0.01*fMaxTheta) |
---|
215 | + IntegralOverAngle(energy2,0.01*fMaxTheta,fMaxTheta) |
---|
216 | + 4.0*sumOdd + 2.0*sumEven )/3.0 ; |
---|
217 | } |
---|
218 | |
---|
219 | |
---|
220 | |
---|
221 | |
---|
222 | // end of G4TransitionRadiation implementation file -------------------------- |
---|