| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // G4RToEConvForGamma class implementation <<
27 // 26 //
28 // Author: H.Kurashige, 05 October 2002 - Firs << 27 // $Id: G4RToEConvForGamma.cc 70745 2013-06-05 10:54:00Z gcosmo $
29 // ------------------------------------------- << 28 //
>> 29 //
>> 30 // --------------------------------------------------------------
>> 31 // GEANT 4 class implementation file/ History:
>> 32 // 5 Oct. 2002, H.Kuirashige : Structure created based on object model
>> 33 // --------------------------------------------------------------
30 34
31 #include "G4RToEConvForGamma.hh" 35 #include "G4RToEConvForGamma.hh"
32 #include "G4ParticleDefinition.hh" 36 #include "G4ParticleDefinition.hh"
33 #include "G4ParticleTable.hh" 37 #include "G4ParticleTable.hh"
>> 38 #include "G4Material.hh"
>> 39 #include "G4PhysicsLogVector.hh"
>> 40
>> 41 #include "G4ios.hh"
34 #include "G4SystemOfUnits.hh" 42 #include "G4SystemOfUnits.hh"
35 #include "G4Log.hh" <<
36 #include "G4Exp.hh" <<
37 #include "G4Pow.hh" <<
38 43
39 // ------------------------------------------- <<
40 G4RToEConvForGamma::G4RToEConvForGamma() 44 G4RToEConvForGamma::G4RToEConvForGamma()
41 : G4VRangeToEnergyConverter() << 45 : G4VRangeToEnergyConverter(),
>> 46 Z(-1),
>> 47 s200keV(0.), s1keV(0.),
>> 48 tmin(0.), tlow(0.),
>> 49 smin(0.), slow(0.),
>> 50 cmin(0.), clow(0.), chigh(0.)
42 { 51 {
43 theParticle = G4ParticleTable::GetParticleTa << 52 theParticle = G4ParticleTable::GetParticleTable()->FindParticle("gamma");
44 if (theParticle == nullptr) << 53 if (theParticle ==0) {
45 { <<
46 #ifdef G4VERBOSE 54 #ifdef G4VERBOSE
47 if (GetVerboseLevel()>0) << 55 if (GetVerboseLevel()>0) {
48 { << 56 G4cout << " G4RToEConvForGamma::G4RToEConvForGamma() ";
49 G4cout << " G4RToEConvForGamma::G4RToECo << 57 G4cout << " Gamma is not defined !!" << G4endl;
50 G4cout << "Gamma is not defined !!" << G <<
51 } 58 }
52 #endif 59 #endif
53 } 60 }
54 else << 61 }
55 { << 62
56 fPDG = theParticle->GetPDGEncoding(); << 63 G4RToEConvForGamma::~G4RToEConvForGamma()
>> 64 {
>> 65 }
>> 66
>> 67
>> 68 // ***********************************************************************
>> 69 // ******************* BuildAbsorptionLengthVector ***********************
>> 70 // ***********************************************************************
>> 71 void G4RToEConvForGamma::BuildAbsorptionLengthVector(
>> 72 const G4Material* aMaterial,
>> 73 G4RangeVector* absorptionLengthVector )
>> 74 {
>> 75 // fill the absorption length vector for this material
>> 76 // absorption length is defined here as
>> 77 //
>> 78 // absorption length = 5./ macroscopic absorption cross section
>> 79 //
>> 80 const G4CrossSectionTable* aCrossSectionTable = (G4CrossSectionTable*)(theLossTable);
>> 81 const G4ElementVector* elementVector = aMaterial->GetElementVector();
>> 82 const G4double* atomicNumDensityVector = aMaterial->GetAtomicNumDensityVector();
>> 83
>> 84 // fill absorption length vector
>> 85 G4int NumEl = aMaterial->GetNumberOfElements();
>> 86 G4double absorptionLengthMax = 0.0;
>> 87 for (size_t ibin=0; ibin<size_t(TotBin); ibin++) {
>> 88 G4double SIGMA = 0. ;
>> 89 for (size_t iel=0; iel<size_t(NumEl); iel++) {
>> 90 G4int IndEl = (*elementVector)[iel]->GetIndex();
>> 91 SIGMA += atomicNumDensityVector[iel]*
>> 92 (*((*aCrossSectionTable)[IndEl]))[ibin];
>> 93 }
>> 94 // absorption length=5./SIGMA
>> 95 absorptionLengthVector->PutValue(ibin, 5./SIGMA);
>> 96 if (absorptionLengthMax < 5./SIGMA ) absorptionLengthMax = 5./SIGMA;
57 } 97 }
58 } 98 }
59 99
60 // ------------------------------------------- << 100
61 G4RToEConvForGamma::~G4RToEConvForGamma() << 101
62 {} << 102 // ***********************************************************************
63 << 103 // ********************** ComputeCrossSection ****************************
64 // ------------------------------------------- << 104 // ***********************************************************************
65 G4double G4RToEConvForGamma::ComputeValue(cons << 105 G4double G4RToEConvForGamma::ComputeCrossSection(G4double AtomicNumber,
66 cons << 106 G4double KineticEnergy)
67 { 107 {
68 // Compute the "absorption" cross-section of << 108 // Compute the "absorption" cross section of the photon "absorption"
69 // Cross-section means here the sum of the c << 109 // cross section means here the sum of the cross sections of the
70 // pair production, Compton scattering and p << 110 // pair production, Compton scattering and photoelectric processes
71 << 111 const G4double t1keV = 1.*keV;
72 const G4double t1keV = 1.*CLHEP::keV; << 112 const G4double t200keV = 200.*keV;
73 const G4double t200keV = 200.*CLHEP::keV; << 113 const G4double t100MeV = 100.*MeV;
74 const G4double t100MeV = 100.*CLHEP::MeV; << 114
75 << 115 // compute Z dependent quantities in the case of a new AtomicNumber
76 G4double Zsquare = Z*Z; << 116 if(std::abs(AtomicNumber-Z)>0.1) {
77 G4double Zlog = G4Pow::GetInstance()->logZ(Z << 117 Z = AtomicNumber;
78 G4double Zlogsquare = Zlog*Zlog; << 118 G4double Zsquare = Z*Z;
79 << 119 G4double Zlog = std::log(Z);
80 G4double tmin = (0.552+218.5/Z+557.17/Zsquar << 120 G4double Zlogsquare = Zlog*Zlog;
81 G4double tlow = 0.2*G4Exp(-7.355/std::sqrt(Z << 121
82 << 122 s200keV = (0.2651-0.1501*Zlog+0.02283*Zlogsquare)*Zsquare;
83 G4double smin = (0.01239+0.005585*Zlog-0.000 << 123 tmin = (0.552+218.5/Z+557.17/Zsquare)*MeV;
84 G4double s200keV = (0.2651-0.1501*Zlog+0.022 << 124 smin = (0.01239+0.005585*Zlog-0.000923*Zlogsquare)*std::exp(1.5*Zlog);
85 << 125 cmin=std::log(s200keV/smin)/(std::log(tmin/t200keV)*std::log(tmin/t200keV));
86 G4double cminlog = G4Log(tmin/t200keV); << 126 tlow = 0.2*std::exp(-7.355/std::sqrt(Z))*MeV;
87 G4double cmin = G4Log(s200keV/smin)/(cminlog << 127 slow = s200keV*std::exp(0.042*Z*std::log(t200keV/tlow)*std::log(t200keV/tlow));
88 << 128 s1keV = 300.*Zsquare;
89 G4double slowlog = G4Log(t200keV/tlow); << 129 clow =std::log(s1keV/slow)/std::log(tlow/t1keV);
90 G4double slow = s200keV * G4Exp(0.042*Z*slow <<
91 G4double logtlow = G4Log(tlow/t1keV); <<
92 G4double clow = G4Log(300.*Zsquare/slow)/log <<
93 G4double chigh = (7.55e-5 - 0.0542e-5*Z)*Zsq <<
94 130
95 // Calculate the cross-section (using an app << 131 chigh=(7.55e-5-0.0542e-5*Z)*Zsquare*Z/std::log(t100MeV/tmin);
96 G4double xs; <<
97 if ( energy < tlow ) <<
98 { <<
99 xs = (energy < t1keV) ? slow*G4Exp(clow*lo <<
100 slow*G4Exp(clow*G4Log(tlow/energy)); <<
101 } <<
102 else if ( energy < t200keV ) <<
103 { <<
104 G4double x = G4Log(t200keV/energy); <<
105 xs = s200keV * G4Exp(0.042*Z*x*x); <<
106 } <<
107 else if( energy<tmin ) <<
108 { <<
109 const G4double x = G4Log(tmin/energy); <<
110 xs = smin * G4Exp(cmin*x*x); <<
111 } 132 }
112 else << 133
113 { << 134 // calculate the cross section (using an approximate empirical formula)
114 xs = smin + chigh*G4Log(energy/tmin); << 135 G4double xs;
>> 136 if ( KineticEnergy<tlow ) {
>> 137 if(KineticEnergy<t1keV) xs = slow*std::exp(clow*std::log(tlow/t1keV));
>> 138 else xs = slow*std::exp(clow*std::log(tlow/KineticEnergy));
>> 139
>> 140 } else if ( KineticEnergy<t200keV ) {
>> 141 xs = s200keV
>> 142 * std::exp(0.042*Z*std::log(t200keV/KineticEnergy)*std::log(t200keV/KineticEnergy));
>> 143
>> 144 } else if( KineticEnergy<tmin ){
>> 145 xs = smin
>> 146 * std::exp(cmin*std::log(tmin/KineticEnergy)*std::log(tmin/KineticEnergy));
>> 147
>> 148 } else {
>> 149 xs = smin + chigh*std::log(KineticEnergy/tmin);
>> 150
115 } 151 }
116 return xs * CLHEP::barn; << 152 return xs * barn;
117 } 153 }
118 154
119 // ------------------------------------------- <<
120 155