Geant4 Cross Reference |
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 // Parameters comes from JQMD 28 // Niita et al., JAERI-Data/Code 99-042 29 // 30 // 230307 Skyrme-QMD parameters added by Y-H. Sato and A. Haga 31 32 #include "G4LightIonQMDParameters.hh" 33 #include "G4Pow.hh" 34 #include "G4PhysicalConstants.hh" 35 36 G4ThreadLocal G4LightIonQMDParameters* G4LightIonQMDParameters::parameters = NULL; 37 38 G4LightIonQMDParameters::G4LightIonQMDParameters() 39 { 40 G4Pow* pow=G4Pow::GetInstance(); 41 42 wl = 1.26; // width of wave packet [fm] 43 hbc = 0.19732857; // h-bar c in GeVfm 44 45 //Pauli 46 cpw = 1.0 / 2.0 / wl; 47 48 cph = 2.0 * wl / (hbc*hbc); 49 50 cpc = 4.0; 51 52 epsx = -20.0 ; 53 54 55 // JQMD 56 /* 57 rho0 = 0.168; // satulation density 58 G4double rpot = 1.0/3.0; 59 60 G4double ebinm = -16.0; // bounding energy [MeV] 61 G4double ebin = ebinm * 0.001; 62 63 G4double pfer = hbc * pow->A13 ( 3./2. *pi*pi * rho0 ); 64 65 G4double rmass = 0.938; 66 67 G4double efer = pfer*pfer / 2. / rmass; 68 69 G4double t3 = 8. / 3. / rpot / pow->powA( rho0 , ( 1.+rpot ) ) * ( efer / 5. - ebin ); 70 71 G4double t0 = -16./15. * efer / rho0 - ( 1.+rpot ) * t3 * pow->powA( rho0 , rpot ); 72 73 74 G4double aaa = 3./4. * t0 * rho0; 75 G4double bbb = 3./8. * t3 * ( 2.+rpot ) * pow->powA( rho0 , ( 1.+rpot ) ); 76 G4double esymm = 25 * 0.001; // symetric potential 25 [MeV] -> GeV 77 78 gamm = rpot + 1.0; 79 */ 80 81 82 // Skyrme-QMD 83 // Ref. Y. Zhang and Z. Li, Elliptic flow and system size dependence of transition energies at intermediate energies, Phys.Rev. C74 (2006) 014602. 84 85 // ImQMD-SLy4 86 /* 87 rho0 = 0.159546; 88 G4double aaa = -297.82 * 0.001; 89 G4double bbb = 219.21 * 0.001; 90 gamm = 7.0/6; 91 eta = 5.0/3; 92 kappas = 0.08; 93 g0 = 24.569/(2 * rho0 * pow->powA( 4 * pi * wl , 1.5 )) * 0.001; 94 g0iso = 4.557/(rho0 * pow->powA( 4 * pi * wl , 1.5 )) * 0.001; 95 gtau0 = 9.70/(pow->powA( rho0 , eta ) * pow->powA ( (4.0*pi*wl) , (1.5*eta) )) * 0.001; 96 G4double esymm = 32 * 0.001; 97 */ 98 99 // ImQMD-SkMstar 100 101 //rho0 = 0.165; 102 rho0 = 0.1603; // satulation density 103 G4double aaa = -318.0 * 0.001; 104 G4double bbb = 249.5 * 0.001; 105 gamm = 7.0/6; 106 eta = 5.0/3; 107 kappas = 0.08; 108 g0 = 21.86/(2 * rho0 * pow->powA( 4 * pi * wl , 1.5 )) * 0.001; 109 //g0iso = -5.485/(rho0 * pow->powA( 4 * pi * wl , 1.5 )) * 0.001; -> kappas 110 gtau0 = 5.9357/(pow->powA( rho0 , eta ) * pow->powA ( (4.0*pi*wl) , (1.5*eta) )) * 0.001; 111 G4double esymm = 32 * 0.001; 112 113 // ImQMD-SIII 114 /* 115 rho0 = 0.1452; // satulation density 116 G4double aaa = -122.921 * 0.001; 117 G4double bbb = 55.343 * 0.001; 118 gamm = 2; 119 eta = 5.0/3; 120 kappas = 0.08; 121 g0 = 18.286/(2 * rho0 * pow->powA( 4 * pi * wl , 1.5 )) * 0.001; 122 //g0iso = -5.485/(rho0 * pow->powA( 4 * pi * wl , 1.5 )) * 0.001; -> kappas 123 gtau0 = 6.439/(pow->powA( rho0 , eta ) * pow->powA ( (4.0*pi*wl) , (1.5*eta) )) * 0.001; 124 G4double esymm = 28.17 * 0.001; 125 */ 126 127 // Local Potenials 128 c0 = aaa / ( rho0 * pow->powA( 4 * pi * wl , 1.5 ) * 2.0 ); 129 130 c3 = bbb / ( pow->powA( rho0 , gamm ) * pow->powA ( (4.0*pi*wl) , (1.5*gamm) ) * ( gamm+1.0) ); 131 132 cs = esymm / ( rho0 * pow->powA( (4.0*pi*wl) , 1.5 ) * 2.0 ); 133 134 G4double ccoul = 0.001439767; 135 cl = ccoul/2.0 * 1; // Include Coulomb interaction 136 //cl = ccoul/2.0 * 0; // Not Include Coulomb interaction 137 138 139 140 // GroundStateNucleus 141 cdp = 1.0 / pow->powA ( ( 4.0 * pi * wl ) , 1.5 ); 142 c0p = c0 * 2.0; 143 c3p = c3 * ( gamm + 1.0 ); 144 csp = cs * 2.0; 145 clp = cl * 2.0; 146 147 g0p = g0 * 2.0; // Skyrme-QMD 148 g0isop = g0iso * 2.0; // Skyrme-QMD 149 gtau0p = gtau0 * ( eta + 1.0 ); // Skyrme-QMD 150 151 } 152 153 154 155 G4LightIonQMDParameters::~G4LightIonQMDParameters() 156 { 157 ; 158 } 159 160