source: trunk/source/processes/electromagnetic/xrays/src/G4TransitionRadiation.cc@ 966

Last change on this file since 966 was 961, checked in by garnier, 17 years ago

update processes
File size: 8.4 KB
Line 
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
56const G4int G4TransitionRadiation::fSympsonNumber = 100 ;
57const G4int G4TransitionRadiation::fGammaNumber = 15 ;
58const G4int G4TransitionRadiation::fPointNumber = 100 ;
59
60using namespace std;
61
62///////////////////////////////////////////////////////////////////////
63//
64// Constructor for selected couple of materials
65//
66
67G4TransitionRadiation::
68G4TransitionRadiation( 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
80G4TransitionRadiation::~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
91G4double
92G4TransitionRadiation::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
117G4double
118G4TransitionRadiation::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
143G4double G4TransitionRadiation::
144AngleIntegralDistribution( 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
193G4double G4TransitionRadiation::
194EnergyIntegralDistribution( 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 --------------------------
Note: See TracBrowser for help on using the repository browser.