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 // * 21 // * Parts of this code which have been devel 22 // * under contract to the European Space Agen 23 // * intellectual property of ESA. Rights to u 24 // * redistribute this software for general pu 25 // * in compliance with any licensing, distrib 26 // * policy adopted by the Geant4 Collaboratio 27 // * written by QinetiQ Ltd for the European S 28 // * contract 17191/03/NL/LvH (Aurora Programm 29 // * 30 // * By using, copying, modifying or distri 31 // * any work based on the software) you ag 32 // * use in resulting scientific publicati 33 // * acceptance of all terms of the Geant4 Sof 34 // ******************************************* 35 // 36 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 37 // 38 // MODULE: G4EMDissociationSpectr 39 // 40 // Version: B.1 41 // Date: 15/04/04 42 // Author: P R Truscott 43 // Organisation: QinetiQ Ltd, UK 44 // Customer: ESA/ESTEC, NOORDWIJK 45 // Contract: 17191/03/NL/LvH 46 // 47 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 48 // 49 // CHANGE HISTORY 50 // -------------- 51 // 52 // 6 October 2003, P R Truscott, QinetiQ Ltd, 53 // Created. 54 // 55 // 15 March 2004, P R Truscott, QinetiQ Ltd, U 56 // Beta release 57 // 58 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 59 ////////////////////////////////////////////// 60 // 61 #include "G4EMDissociationSpectrum.hh" 62 #include "globals.hh" 63 #include "G4PhysicalConstants.hh" 64 #include "G4SystemOfUnits.hh" 65 #include "G4Pow.hh" 66 67 ////////////////////////////////////////////// 68 // 69 G4EMDissociationSpectrum::G4EMDissociationSpec 70 { 71 bessel = new G4Bessel(); 72 } 73 ////////////////////////////////////////////// 74 // 75 G4EMDissociationSpectrum::~G4EMDissociationSpe 76 { 77 delete bessel; 78 } 79 ////////////////////////////////////////////// 80 // 81 G4double G4EMDissociationSpectrum::GetGeneralE 82 (G4double Eg, G4double b, G4double bmin) 83 { 84 G4double b2 = b*b; 85 G4double gg = 1.0/std::sqrt(1.0-b2); 86 G4double xi = Eg * bmin / gg / b / hbarc; 87 G4double K0 = bessel->K0(xi); 88 G4double K1 = bessel->K1(xi); 89 G4double n = 2.0 * fine_structure_const / p 90 (xi*K0*K1 - xi*xi*b2/2.0*(K1*K1-K0*K0)); 91 return n; 92 } 93 ////////////////////////////////////////////// 94 // 95 G4double G4EMDissociationSpectrum::GetGeneralE 96 (G4double Eg, G4double b, G4double bmin) 97 { 98 G4double b2 = b * b; 99 G4double b4 = b2 * b2; 100 G4double gg = 1.0/std::sqrt(1.0-b2); 101 G4double xi = Eg * bmin / gg / b / hbarc; 102 G4double K0 = bessel->K0(xi); 103 G4double K1 = bessel->K1(xi); 104 G4double n = 2.0 * fine_structure_const / p 105 (2.0*(1.0-b2)*K1*K1 + xi*G4Pow::GetInstanc 106 xi*xi*b4/2.0*(K1*K1-K0*K0)); 107 return n; 108 } 109 ////////////////////////////////////////////// 110 // 111 G4double G4EMDissociationSpectrum::GetClosestA 112 (const G4double AP, const G4double ZP, G4dou 113 { 114 G4double bsq = b * b; 115 G4double gg = 1.0/std::sqrt(1-bsq); 116 G4double AProot3 = G4Pow::GetInstance()->pow 117 G4double ATroot3 = G4Pow::GetInstance()->pow 118 G4double bc = 1.34 * fermi * (AProot3+A 119 // G4double a0 = ZP * ZT * classic_elect 120 G4double a0 = ZP * ZT * elm_coupling / 121 G4double bmin = 1.25 * bc + halfpi*a0/gg; 122 return bmin; 123 } 124 ////////////////////////////////////////////// 125 // 126