| Geant4 Cross Reference |
>> 1 // This code implementation is the intellectual property of
>> 2 // the GEANT4 collaboration.
1 // 3 //
2 // ******************************************* << 4 // By copying, distributing or modifying the Program (or any work
3 // * License and Disclaimer << 5 // based on the Program) you indicate your acceptance of this statement,
4 // * << 6 // and all its terms.
5 // * The Geant4 software is copyright of th << 7 //
6 // * the Geant4 Collaboration. It is provided << 8 // $Id: G4TransitionRadiation.hh,v 1.2.8.1 1999/12/07 20:51:24 gunter Exp $
7 // * conditions of the Geant4 Software License << 9 // GEANT4 tag $Name: geant4-01-01 $
8 // * LICENSE and available at http://cern.ch/ << 10 //
9 // * include a list of copyright holders. << 11 // G4TransitionRadiation -- header file
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 // 12 //
26 // Class for description of transition radiat 13 // Class for description of transition radiation generated
27 // by charged particle crossed interface betw 14 // by charged particle crossed interface between material 1
28 // and material 2 (1 -> 2). Transition radiati 15 // and material 2 (1 -> 2). Transition radiation could be of kind:
29 // - optical back 16 // - optical back
30 // - optical forward 17 // - optical forward
31 // - X-ray forward (for relativistic case Tk 18 // - X-ray forward (for relativistic case Tkin/mass >= 10^2)
32 // 19 //
33 // GEANT 4 class header file --- Copyright CER 20 // GEANT 4 class header file --- Copyright CERN 1995
>> 21 // CERB Geneva Switzerland
34 // 22 //
>> 23 // for information related to this code, please, contact
>> 24 // CERN, CN Division, ASD Group
35 // History: 25 // History:
36 // 18.12.97, V. Grichine (Vladimir.Grichine@ce 26 // 18.12.97, V. Grichine (Vladimir.Grichine@cern.ch)
37 // 02.02.00, V.Grichine, new data fEnergy and << 27
38 // numerical integration <<
39 // 03.06.03, V.Ivanchenko fix compilation warn <<
40 // 28.07.05, P.Gumplinger add G4ProcessType to <<
41 28
42 #ifndef G4TransitionRadiation_h 29 #ifndef G4TransitionRadiation_h
43 #define G4TransitionRadiation_h 30 #define G4TransitionRadiation_h
44 31
45 #include "globals.hh" << 32
46 #include "G4ParticleDefinition.hh" <<
47 #include "G4Step.hh" <<
48 #include "G4Track.hh" <<
49 #include "G4VDiscreteProcess.hh" 33 #include "G4VDiscreteProcess.hh"
50 #include "G4VParticleChange.hh" << 34 #include "G4Material.hh"
>> 35 // #include "G4OpBoundaryProcess.hh"
51 36
52 class G4TransitionRadiation : public G4VDiscre << 37 class G4TransitionRadiation : public G4VDiscreteProcess
53 { 38 {
54 public: << 39 public:
55 explicit G4TransitionRadiation(const G4Strin << 40
56 G4ProcessType << 41 // Constructors
>> 42
>> 43
>> 44 G4TransitionRadiation( const G4String& processName = "TR") ;
>> 45
>> 46
>> 47 // G4TransitionRadiation(const G4TransitionRadiation& right) ;
>> 48
>> 49 // Destructor
>> 50
>> 51 ~G4TransitionRadiation() ;
>> 52
>> 53 // Operators
>> 54 // G4TransitionRadiation& operator=(const G4TransitionRadiation& right) ;
>> 55 // G4int operator==(const G4TransitionRadiation& right)const ;
>> 56 // G4int operator!=(const G4TransitionRadiation& right)const ;
>> 57
>> 58 // Methods
>> 59
>> 60 G4bool IsApplicable(const G4ParticleDefinition& aParticleType)
>> 61 {
>> 62 return ( aParticleType.GetPDGCharge() != 0.0 );
>> 63 }
>> 64
>> 65 G4double GetMeanFreePath(const G4Track& aTrack,
>> 66 G4double previousStepSize,
>> 67 G4ForceCondition* condition)
>> 68 {
>> 69 *condition = Forced;
>> 70 return DBL_MAX; // so TR doesn't limit mean free path
>> 71 }
>> 72
>> 73 G4VParticleChange* PostStepDoIt(const G4Track& aTrack,
>> 74 const G4Step& aStep)
>> 75 {
>> 76 ClearNumberOfInteractionLengthLeft();
>> 77 return &aParticleChange;
>> 78 }
>> 79
>> 80
>> 81
>> 82
>> 83 virtual
>> 84 G4double SpectralAngleTRdensity( G4double energy,
>> 85 G4double varAngle ) const = 0 ;
57 86
58 virtual ~G4TransitionRadiation(); << 87 G4double IntegralOverEnergy( G4double energy1,
>> 88 G4double energy2,
>> 89 G4double varAngle ) const ;
59 90
60 G4TransitionRadiation(const G4TransitionRadi << 91 G4double IntegralOverAngle( G4double energy,
61 G4TransitionRadiation& operator=(const G4Tra << 92 G4double varAngle1,
>> 93 G4double varAngle2 ) const ;
62 94
63 // Methods << 95 G4double AngleIntegralDistribution( G4double varAngle1,
>> 96 G4double varAngle2 ) const ;
64 97
65 G4bool IsApplicable(const G4ParticleDefiniti << 98 G4double EnergyIntegralDistribution( G4double energy1,
>> 99 G4double energy2 ) const ;
66 100
67 virtual G4double GetMeanFreePath(const G4Tra <<
68 G4ForceCond <<
69 101
70 virtual G4VParticleChange* PostStepDoIt(cons <<
71 cons <<
72 102
73 virtual void ProcessDescription(std::ostream << 103 // Access functions
74 virtual void DumpInfo() const override { Pro <<
75 104
76 virtual G4double SpectralAngleTRdensity(G4do <<
77 G4do <<
78 105
79 G4double IntegralOverEnergy(G4double energy1 << 106 protected :
80 G4double varAngl <<
81 107
82 G4double IntegralOverAngle(G4double energy, << 108 G4int fMatIndex1 ; // index of the 1st material
83 G4double varAngle << 109 G4int fMatIndex2 ; // index of the 2nd material
84 110
85 G4double AngleIntegralDistribution(G4double << 111 private :
86 G4double <<
87 112
88 G4double EnergyIntegralDistribution(G4double << 113 G4double fGamma ;
89 114
90 protected: << 115 // Local constants
91 // Local constants << 116 static const G4int fSympsonNumber ; // Accuracy of Sympson integration 10
92 // Accuracy of Sympson integration << 117 static const G4int fGammaNumber ; // = 15
93 static constexpr G4int fSympsonNumber = 100; << 118 static const G4int fPointNumber ; // = 100
94 static constexpr G4int fGammaNumber = 15; <<
95 static constexpr G4int fPointNumber = 100; <<
96 119
97 G4double fGamma; << 120 G4double fMinEnergy ; // min TR energy
98 G4double fEnergy; << 121 G4double fMaxEnergy ; // max TR energy
99 G4double fVarAngle; << 122 G4double fMaxTheta ; // max theta of TR quanta
100 123
101 G4double fMinEnergy; // min TR energy << 124 G4double fSigma1 ; // plasma energy Sq of matter1
102 G4double fMaxEnergy; // max TR energy << 125 G4double fSigma2 ; // plasma energy Sq of matter2
103 G4double fMaxTheta; // max theta of TR qu <<
104 126
105 G4double fSigma1; // plasma energy Sq of ma <<
106 G4double fSigma2; // plasma energy Sq of ma <<
107 127
108 G4int fMatIndex1; // index of the 1st mater << 128 } ;
109 G4int fMatIndex2; // index of the 2nd mater <<
110 }; <<
111 129
112 #endif // G4TransitionRadiation_h << 130 #endif // G4TransitionRadiation_h
113 131