Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4NuclearLevelData.cc
Go to the documentation of this file.
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// -------------------------------------------------------------------
28//
29// GEANT4 source file
30//
31// File name: G4NuclearLevelData
32//
33// Author: V.Ivanchenko
34//
35// Creation date: 10 February 2015
36//
37// Modifications:
38//
39// -------------------------------------------------------------------
40
41#include "G4NuclearLevelData.hh"
42#include "G4LevelReader.hh"
43#include "G4LevelManager.hh"
44#include "G4Element.hh"
45#include "G4ElementTable.hh"
48#include "G4ShellCorrection.hh"
49#include "G4SystemOfUnits.hh"
50#include "G4AutoLock.hh"
51#include "G4Pow.hh"
52#include <iomanip>
53
54G4NuclearLevelData* G4NuclearLevelData::theInstance = nullptr;
55
56const G4int G4NuclearLevelData::AMIN[] = {1,
57 1,3,4,5,7,8,10,12,14,16, // Z= 1-10
58 20,19,21,22,25,27,29,31,35,35, // Z= 11-20
59 39,39,42,42,44,45,50,48,53,54, // Z= 21-30
60 60,60,63,64,69,69,73,73,76,78, // Z= 31-40
61 81,83,85,88,90,91,93,96,98,100, // Z= 41-50
62 104,105,108,109,112,114,117,121,121,125, // Z= 51-60
63 128,129,130,134,135,139,140,143,144,149, // Z= 61-70
64 150,153,155,158,159,161,164,166,170,171, // Z= 71-80
65 176,178,184,186,191,193,199,201,206,208, // Z= 81-90
66 212,217,225,228,230,233,234,237,241,241, // Z= 91-100
67 245,248,252,253,255,258,260,263,266,269, // Z=101-110
68 272,277,282,0,289,0,293}; // Z=111-117
69const G4int G4NuclearLevelData::AMAX[] = {1,
70 6,10,11,14,17,22,24,26,30,30, // Z= 1-10
71 37,40,42,44,46,48,51,53,56,58, // Z= 11-20
72 61,63,66,68,71,74,76,79,82,85, // Z= 21-30
73 87,90,92,95,98,101,103,107,109,112, // Z= 31-40
74 114,117,120,124,126,128,129,133,135,138, // Z= 41-50
75 140,143,145,148,151,152,153,155,155,161, // Z= 51-60
76 163,165,167,169,171,173,175,177,179,181, // Z= 61-70
77 184,189,192,194,198,202,204,205,210,216, // Z= 71-80,
78 217,220,224,227,229,229,233,234,235,238, // Z= 81-90,
79 239,242,244,247,248,252,253,256,257,259, // Z= 91-100
80 260,262,266,265,270,269,274,277,278,277, // Z=101-110
81 282,281,286,0,290,0,294}; // Z=111-117
82const G4int G4NuclearLevelData::LEVELIDX[] = {0,
83 0, 6, 14, 22, 32, 43, 58, 73, 88, 105, // Z= 1-10
84 120, 138, 160, 182, 205, 227, 249, 272, 295, 317, // Z= 11-20
85 341, 364, 389, 414, 441, 469, 499, 526, 558, 588, // Z= 21-30
86 620, 648, 679, 709, 741, 771, 804, 835, 870, 904, // Z= 31-40
87 939, 973, 1008, 1044, 1081, 1118, 1156, 1193, 1231, 1269, // Z= 41-50
88 1308, 1345, 1384, 1422, 1462, 1502, 1541, 1578, 1613, 1648, // Z= 51-60
89 1685, 1721, 1758, 1796, 1832, 1869, 1904, 1940, 1975, 2011, // Z= 61-70
90 2044, 2079, 2116, 2154, 2191, 2231, 2273, 2314, 2354, 2395, // Z= 71-80
91 2441, 2483, 2526, 2567, 2609, 2648, 2685, 2720, 2754, 2784, // Z= 81-90
92 2815, 2843, 2869, 2889, 2909, 2928, 2948, 2968, 2988, 3005, // Z= 91-100
93 3024, 3040, 3055, 3070, 3083, 3099, 3111, 3126, 3141, 3154, // Z= 101-110
94 3163, 3174, 3179, 3184, 3184, 3186, 3186}; // Z= 111-117
95
96// obtained from PhotonEvaporation5.2
97static const G4float LEVELMAX[3188] = {0.0f,
98 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1-10
99 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 5.366f, 4.63f, 2.255f, 2.691f, 0.0f, //11-20
100 0.0f, 0.0f, 0.0f,0.42908f, 27.494f, 17.495f, 7.371f,0.32004f, 0.0f, 0.0f, //21-30
101 0.0f, 0.0f, 2.32f, 17.076f, 8.895f, 9.2717f, 2.723f, 4.131f, 0.0f, 0.0f, //31-40
102 0.0f, 0.0f, 0.0f, 0.0f, 3.3537f, 10.679f, 16.57f, 15.108f, 11.306f, 0.74f, //41-50
103 4.142f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //51-60
104 15.3f, 11.05f, 13.39f,0.39727f, 5.772f, 2.6144f, 0.0f, 0.0f, 4.195f, 0.0f, //61-70
105 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 11.846f, 13.09f, 11.079f, 8.282f, 3.9449f, //71-80
106 6.555f, 3.073f, 6.936f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //81-90
107 11.193f, 7.584f, 10.411f, 6.6475f, 4.584f, 2.5722f, 6.873f, 1.8313f, 5.45f, 0.665f, //91-100
108 1.281f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.6164f, 4.635f, 18.43f, 11.989f, //101-110
109 11.533f, 4.4361f, 8.393f, 4.07f, 3.6913f, 0.885f, 3.904f, 0.931f, 2.235f, 0.0f, //111-120
110 4.4679f, 9.813f, 12.592f, 6.578f, 4.289f, 0.2336f, 9.1867f, 3.5125f, 4.166f, 2.1136f, //121-130
111 1.162f, 0.569f, 1.117f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //131-140
112 6.2542f, 9.596f, 15.15f, 13.332f, 12.479f, 5.9258f, 8.4394f, 3.9857f, 5.4129f, 3.8145f, //141-150
113 5.203f, 1.2424f, 2.12f, 0.67f, 0.66f, 0.0f, 2.016f, 0.0f, 0.0f, 0.0f, //151-160
114 0.0f, 0.55f, 0.4258f, 7.97f, 6.695f, 8.952f, 7.4088f, 5.9936f, 9.3731f, 5.1495f, //161-170
115 3.202f, 4.7304f, 0.657f, 1.02f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //171-180
116 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 5.677f, 8.3758f, 13.744f, 10.252f, 15.529f, //181-190
117 6.5874f, 6.705f, 4.931f, 6.0233f, 2.1682f, 3.692f, 1.437f, 1.084f, 0.0f, 1.831f, //191-200
118 0.0f, 0.77f, 0.0f, 0.0f, 0.0f, 0.0f, 1.12f, 2.104f, 8.379f, 8.0143f, //201-210
119 13.879f, 7.9356f, 7.9666f, 7.9201f, 6.2205f, 2.281f, 7.897f, 0.0f, 2.167f, 0.0f, //211-220
120 1.56f, 0.0f, 0.184f, 0.0f, 0.0f, 0.0f, 0.0f, 1.507f, 2.887f, 5.136f, //221-230
121 10.146f, 12.048f, 10.779f, 16.649f, 6.9861f, 7.2719f, 4.3036f, 6.0053f, 0.0f, 5.0094f, //231-240
122 0.904f, 0.89f, 0.94f, 3.257f, 0.0f, 0.952f, 0.0f, 0.0f, 0.0f, 0.0f, //241-250
123 0.0f, 5.0463f, 5.5484f, 8.3055f, 12.572f, 8.5798f, 17.009f, 4.836f, 4.3544f, 4.0871f, //251-260
124 2.451f, 3.0304f, 1.829f, 0.475f, 1.616f, 1.024f, 0.0f, 0.0f, 0.0f, 0.0f, //261-270
125 0.0f, 0.0f, 1.867f, 3.819f, 5.542f, 12.276f, 22.365f, 8.7912f, 21.662f, 5.5355f, //271-280
126 10.857f, 4.271f, 4.6339f, 4.2469f, 5.3515f, 4.3261f, 1.577f, 0.0f, 3.279f, 0.0f, //281-290
127 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.9179f, 5.0498f, 10.981f, 18.612f, 7.4724f, //291-300
128 9.7407f, 5.978f, 3.7142f, 2.5743f, 4.5691f, 1.9443f, 5.8424f, 2.177f, 2.1042f, 0.1714f, //301-310
129 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.045f, 3.842f, 4.3837f, //311-320
130 7.75f, 22.06f, 7.146f, 10.45f, 6.2237f, 12.188f, 4.9997f, 6.626f, 5.5034f, 9.5457f, //321-330
131 5.6802f, 6.8692f, 3.53f, 2.5631f, 2.2f, 3.699f, 0.0f, 0.0f, 0.0f, 0.0f, //331-340
132 0.0f, 0.0f, 6.148f, 6.2534f, 7.3549f, 9.1411f, 16.462f, 5.5415f, 10.31f, 6.6775f, //341-350
133 15.662f, 1.8478f, 6.184f, 4.2657f, 3.382f, 0.247f, 0.0f, 0.7749f, 0.0f, 0.0f, //351-360
134 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.744f, 3.0664f, 14.55f, 13.03f, //361-370
135 13.169f, 8.0051f, 9.977f, 8.8816f, 10.939f, 4.8821f, 9.0862f, 8.1317f, 6.4321f, 5.4631f, //371-380
136 4.2098f, 0.0f, 2.881f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //381-390
137 0.0f, 7.1595f, 13.631f, 15.258f, 13.282f, 9.1953f, 8.59f, 11.604f, 5.0389f, 4.0852f, //391-400
138 2.2979f, 2.1528f, 0.0f, 2.4756f, 0.114f, 0.0f, 0.2021f, 0.0f, 0.0f, 0.0f, //401-410
139 0.081f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.494f, 8.1625f, 10.022f, 17.378f, //411-420
140 10.7f, 17.958f, 9.2616f, 12.795f, 11.291f, 12.54f, 7.4333f, 5.6033f, 12.95f, 4.6797f, //421-430
141 0.503f, 4.6814f, 1.027f, 1.175f, 0.0f, 1.141f, 0.0f, 0.0f, 0.0f, 0.0f, //431-440
142 0.0f, 0.0f, 0.0f, 0.0f, 6.2931f, 10.724f, 8.2774f, 19.636f, 12.066f, 11.159f, //441-450
143 4.9981f, 10.559f, 7.4078f, 4.7103f, 9.8311f, 1.3507f, 9.235f, 2.3782f, 1.5f, 0.0f, //451-460
144 0.188f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //461-470
145 0.0f, 0.0f, 0.153f, 6.9944f, 12.65f, 7.3819f, 12.593f, 15.062f, 10.243f, 11.964f, //471-480
146 8.323f, 12.813f, 4.9852f, 17.956f, 0.861f, 5.474f, 2.2231f, 4.6284f, 2.2919f, 1.407f, //481-490
147 0.387f, 1.389f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //491-500
148 0.0f, 4.1151f, 10.507f, 14.882f, 18.6f, 11.292f, 8.6251f, 9.5531f, 8.1224f, 4.8042f, //501-510
149 2.3097f, 4.1671f, 1.561f, 3.2712f, 0.642f, 0.0f, 1.8165f, 0.0f, 0.0f, 0.0f, //511-520
150 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.247f, //521-530
151 1.4534f, 6.4574f, 4.4831f, 22.459f, 8.3457f, 42.007f, 22.24f, 22.997f, 7.82f, 10.596f, //531-540
152 4.461f, 8.24f, 4.0111f, 6.5798f, 2.155f, 5.7752f, 2.7f, 3.3612f, 1.2734f, 3.593f, //541-550
153 0.0f, 1.763f, 0.0f, 2.42f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //551-560
154 0.0f, 2.398f, 22.747f, 31.964f, 8.1327f, 9.4062f, 10.884f, 9.8643f, 4.7634f, 8.4843f, //561-570
155 5.0772f, 3.4636f, 0.9563f, 3.828f, 1.9801f, 5.3307f, 2.5968f, 2.386f, 0.0f, 0.0662f, //571-580
156 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //581-590
157 0.0f, 3.378f, 3.3862f, 33.899f, 32.617f, 42.521f, 9.774f, 24.869f, 19.887f, 7.6933f, //591-600
158 4.684f, 18.016f, 2.378f, 6.1158f, 2.377f, 4.7703f, 2.0089f, 5.628f, 5.022f, 2.9744f, //601-610
159 4.6053f, 3.106f, 3.304f, 1.492f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //611-620
160 2.903f, 6.8423f, 14.465f, 5.6291f, 15.981f, 5.1092f, 10.083f, 7.7252f, 10.257f, 2.8864f, //621-630
161 4.1652f,0.74133f, 2.9895f, 1.0857f, 4.1479f, 2.6024f, 3.9459f, 3.5539f, 3.3348f, 3.38f, //631-640
162 1.6216f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.285f, //641-650
163 0.0f, 9.2999f, 8.1524f, 18.08f, 4.848f, 22.959f, 9.1823f, 13.173f, 7.4159f, 8.8679f, //651-660
164 2.3614f, 8.3611f, 6.5056f, 6.0652f, 6.0714f, 5.0782f, 4.2777f, 6.1553f, 4.4705f, 4.402f, //661-670
165 1.9415f, 1.6704f, 2.3482f, 0.527f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //671-680
166 0.0f, 9.7929f, 7.7913f, 8.4997f, 18.181f, 8.9417f, 9.6844f, 6.088f, 7.8614f, 1.129f, //681-690
167 9.399f, 1.0645f, 4.4563f, 2.024f, 1.9648f, 1.8731f, 3.5627f, 1.092f, 4.8418f, 1.1983f, //691-700
168 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //701-710
169 0.929f, 0.0f, 21.201f, 13.363f, 20.245f, 6.0355f, 20.798f, 10.468f, 13.202f, 10.243f, //711-720
170 5.4325f, 7.4197f, 5.7838f, 3.6768f, 7.8185f, 3.2227f, 5.6886f, 4.9884f, 6.6046f, 5.1648f, //721-730
171 3.302f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //731-740
172 0.0f, 21.412f, 14.959f, 18.968f, 25.968f, 7.6143f, 13.682f, 9.415f, 11.344f, 6.0869f, //741-750
173 8.8115f, 4.4501f, 3.7987f, 2.2433f, 5.635f, 2.7416f, 5.3906f, 3.814f, 3.9871f, 3.1541f, //751-760
174 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //761-770
175 0.0f, 0.0f, 0.0f, 23.62f, 23.237f, 30.936f, 19.135f, 14.735f, 17.354f, 17.296f, //771-780
176 11.822f, 11.484f, 5.7594f, 6.0117f, 6.3734f, 7.6532f, 4.7906f, 10.116f, 5.8215f, 7.9695f, //781-790
177 4.7073f, 5.7301f, 4.4529f, 5.0788f, 3.1999f, 2.52f, 0.1955f, 0.5541f, 0.0f, 0.0f, //791-800
178 0.0f, 0.0f, 0.0f, 0.0f, 21.612f, 32.365f, 2.0191f, 18.376f, 9.7229f, 15.984f, //801-810
179 13.176f, 10.027f, 6.0137f, 9.342f, 7.3827f, 7.1071f, 8.651f, 3.3085f, 5.905f, 7.3913f, //811-820
180 3.8813f, 6.2389f, 4.7881f, 6.7256f, 3.9442f, 0.835f, 1.1346f, 0.0f, 0.0f, 0.0f, //821-830
181 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.043f, 0.0f, 14.23f, 13.943f, 17.764f, //831-840
182 8.1749f, 21.818f, 24.46f, 18.674f, 19.803f, 15.08f, 7.5549f, 10.874f, 7.0323f, 12.026f, //841-850
183 7.9844f, 5.8279f, 5.742f, 7.363f, 6.7074f, 6.0637f, 6.176f, 7.5265f, 6.3059f, 3.6224f, //851-860
184 2.3208f, 3.3456f, 1.3629f, 0.126f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //861-870
185 0.0f, 0.787f, 5.702f, 12.46f, 17.67f, 18.897f, 19.47f, 17.8f, 11.099f, 12.415f, //871-880
186 7.0166f, 9.6184f, 11.991f, 4.8351f, 6.8961f, 4.0479f, 4.314f, 2.9699f, 5.8102f, 1.9836f, //881-890
187 3.5226f, 1.897f, 3.3895f, 1.7f, 2.6959f, 1.6896f, 0.0f, 0.0f, 0.0f, 0.0f, //891-900
188 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.789f, 7.7086f, 16.12f, 23.223f, 36.877f, //901-910
189 11.141f, 20.532f, 10.092f, 11.2f, 9.6011f, 13.11f, 12.773f, 9.7222f, 7.2941f, 8.9806f, //911-920
190 6.5618f, 6.8213f, 7.2959f, 6.5436f, 6.567f, 7.6149f, 6.5062f, 7.4517f, 4.028f, 4.2239f, //921-930
191 0.0f, 0.607f, 0.0f, 0.5216f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.18f, //931-940
192 1.723f, 9.093f, 15.145f, 20.786f, 18.374f, 9.7372f, 7.2726f, 3.0718f, 10.137f, 3.7969f, //941-950
193 10.955f, 6.496f, 8.6952f, 1.6147f, 2.2475f, 0.0f, 2.3363f, 0.7036f, 3.3966f, 2.2864f, //951-960
194 3.6071f, 0.5143f, 1.0456f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //961-970
195 0.0f, 0.0f, 0.0f, 4.2764f, 5.5996f, 12.766f, 11.838f, 14.058f, 7.5901f, 14.486f, //971-980
196 10.749f, 11.215f, 9.67f, 11.588f, 10.509f, 9.8824f, 5.5025f, 4.6101f, 8.118f, 8.1137f, //981-990
197 5.031f, 6.2005f, 6.308f, 7.2832f, 6.0757f, 7.66f, 5.057f, 5.3475f, 1.4436f, 2.5319f, //991-1000
198 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.524f, 1.599f, //1001-1010
199 4.9129f, 20.291f, 11.246f, 20.862f, 7.8331f, 13.258f, 9.1864f, 14.461f, 10.047f, 8.3451f, //1011-1020
200 3.1294f, 4.1779f, 4.3581f, 4.231f,0.35986f, 4.069f, 1.0173f, 2.7664f, 2.3395f, 4.3206f, //1021-1030
201 2.359f, 4.2009f, 3.6771f, 3.9519f, 0.35f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1031-1040
202 0.0f, 0.0f, 0.0f, 3.4808f, 6.2392f, 6.3879f, 9.63f, 11.311f, 11.95f, 18.321f, //1041-1050
203 12.161f, 11.601f, 17.292f, 17.237f, 19.254f, 20.198f, 9.8262f, 12.222f, 5.127f, 5.357f, //1051-1060
204 3.6412f, 4.4457f, 3.7652f, 5.1544f, 2.7331f, 8.1598f, 7.5439f, 7.7493f, 3.4822f, 4.3195f, //1061-1070
205 0.0617f, 2.1664f, 0.0f, 1.4156f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1071-1080
206 0.0f, 7.0196f, 9.744f, 6.9245f, 10.426f, 13.876f, 5.56f, 9.2422f, 9.4964f, 10.282f, //1081-1090
207 14.519f, 17.234f, 15.978f, 7.9989f, 4.4759f, 4.7022f, 6.671f, 3.8021f, 2.2865f, 3.3456f, //1091-1100
208 1.7662f, 4.2482f, 2.7694f, 4.006f, 2.0203f, 2.2801f, 0.0f, 0.0f, 0.0f, 0.0f, //1101-1110
209 0.1572f, 0.0f, 0.271f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.5358f, 7.6629f, //1111-1120
210 7.7002f, 10.418f, 7.04f, 7.523f, 10.863f, 4.5178f, 16.108f, 6.4882f, 15.414f, 17.357f, //1121-1130
211 12.707f, 16.909f, 5.895f, 10.676f, 6.8279f, 4.958f, 4.4841f, 3.831f, 5.2219f, 4.5176f, //1131-1140
212 5.2557f, 4.0712f, 7.8188f, 3.886f, 3.3024f, 0.0f, 3.4705f, 0.135f, 1.164f, 0.0f, //1141-1150
213 1.3f, 0.0f, 2.1097f, 0.0f, 2.151f, 0.0f, 0.0f, 4.86f, 0.0f, 6.9485f, //1151-1160
214 4.4757f, 7.7702f, 8.699f, 10.706f, 6.1075f, 8.2579f, 7.159f, 12.619f, 6.7608f, 8.7181f, //1161-1170
215 4.3106f, 5.4145f, 2.6661f, 2.3527f, 1.8189f,0.78316f, 1.8858f, 3.239f, 2.672f, 0.3239f, //1171-1180
2160.72042f, 0.326f, 0.8058f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1181-1190
217 0.0f, 0.0f, 0.0f, 0.0f, 2.4275f, 7.0998f, 9.3882f, 7.1789f, 10.514f, 11.878f, //1191-1200
218 18.243f, 10.851f, 15.863f, 11.852f, 20.979f, 11.133f, 14.206f, 9.4072f, 9.3942f, 4.2015f, //1201-1210
219 4.6052f, 3.8328f, 4.9689f, 4.5389f, 4.367f, 3.9881f, 5.5227f, 3.3902f, 3.1702f, 2.9102f, //1211-1220
220 3.6636f, 3.6065f, 3.151f, 3.2399f, 1.429f, 1.94f, 2.13f, 0.0f, 0.0f, 0.0f, //1221-1230
221 0.0f, 0.0f, 2.721f, 6.0742f, 9.0586f, 6.3034f, 6.1527f, 15.283f, 5.4832f, 6.989f, //1231-1240
222 8.7931f, 6.3857f, 12.745f, 13.503f, 8.328f, 11.96f, 5.485f, 4.5831f, 1.4511f, 4.1875f, //1241-1250
223 0.2f, 4.0259f, 0.0f, 3.89f, 0.0f, 2.7235f,0.24268f, 2.8633f,0.68823f, 3.5895f, //1251-1260
224 1.1729f, 4.5789f, 5.3908f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.1717f, //1261-1270
225 1.969f, 1.7846f, 9.998f, 10.288f, 19.932f, 6.0447f, 25.56f, 24.742f, 11.516f, 17.813f, //1271-1280
226 12.965f, 18.22f, 17.871f, 9.0696f, 9.3219f, 4.0436f, 7.01f, 3.1033f, 8.993f, 3.228f, //1281-1290
227 4.2838f, 3.306f, 8.4335f, 2.6235f, 4.9902f, 3.8995f, 4.898f, 3.9925f, 5.2628f, 6.6536f, //1291-1300
228 7.2441f, 2.0046f, 2.5089f, 0.0f, 1.295f, 0.0f, 1.344f, 0.0f, 3.9738f, 6.7841f, //1301-1310
229 14.217f, 7.2145f, 25.458f, 27.548f, 23.013f, 18.439f, 21.104f, 16.268f, 10.306f, 4.2456f, //1311-1320
230 8.1846f, 2.9184f, 6.1139f, 2.8847f, 3.9728f,0.92057f, 7.3097f, 1.0596f, 5.3651f, 0.1279f, //1321-1330
231 5.3543f, 0.833f, 3.2807f, 2.9831f, 3.5687f, 5.1092f, 6.9509f, 5.324f, 3.263f, 0.0f, //1331-1340
232 0.0f, 0.0f, 0.0f, 0.0f, 0.15f, 3.127f, 0.721f, 10.546f, 13.752f, 20.637f, //1341-1350
233 9.9578f, 28.646f, 8.7643f, 32.926f, 14.783f, 23.721f, 10.426f, 11.851f, 9.6987f, 6.0391f, //1351-1360
234 8.5832f, 6.9141f, 8.0304f, 5.7514f, 3.5642f, 7.9153f, 6.2876f, 7.724f, 6.0824f, 7.6365f, //1361-1370
235 5.1722f, 4.8903f, 6.1635f, 7.7224f, 6.6692f, 5.16f, 3.6273f, 2.7603f, 2.2105f, 0.0f, //1371-1380
236 0.0f, 0.0f, 0.0f, 0.0f, 4.6208f, 7.075f, 13.017f, 8.7122f, 31.621f, 8.576f, //1381-1390
237 19.449f, 11.403f, 21.128f, 2.8059f, 11.823f, 6.3618f, 11.143f, 8.3232f, 8.9855f, 5.5555f, //1391-1400
238 3.8679f, 1.1022f, 3.3503f, 4.1498f, 3.408f, 1.079f, 2.3522f,0.27786f, 4.0473f, 1.1065f, //1401-1410
239 5.8492f, 4.0745f, 3.8528f, 0.0f, 3.3323f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1411-1420
240 0.0f, 0.0f, 1.8897f, 0.0f, 4.4691f, 13.218f, 34.591f, 11.183f, 23.483f, 24.255f, //1421-1430
241 15.399f, 20.542f, 12.15f, 13.651f, 13.473f, 9.4266f, 16.53f, 37.618f, 8.7448f, 9.523f, //1431-1440
242 5.5734f, 3.8291f, 6.9713f, 3.1859f, 4.1884f, 2.1237f, 3.4923f, 3.1699f, 6.6241f, 5.1702f, //1441-1450
243 5.814f, 5.0965f, 5.5036f, 3.3657f, 4.5112f, 1.9627f, 2.3027f, 0.0f, 0.0f, 0.0f, //1451-1460
244 0.0f, 2.1434f, 8.29f, 6.1544f, 0.0f, 10.381f, 14.396f, 0.0f, 4.891f, 7.126f, //1461-1470
245 12.184f, 7.873f, 11.233f, 6.127f, 11.31f, 8.687f, 5.8654f, 4.9467f, 8.0995f, 4.0405f, //1471-1480
246 5.7351f, 5.6965f, 2.8342f, 2.2537f, 1.6329f, 0.0f, 5.4942f, 2.5084f, 4.2277f, 2.9303f, //1481-1490
247 4.3457f, 4.0013f, 3.1181f, 1.426f, 1.5108f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1491-1500
248 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 6.446f, 19.082f, 20.601f, 9.0649f, 10.75f, //1501-1510
249 13.551f, 25.371f, 11.179f, 16.895f, 13.477f, 18.217f, 10.388f, 8.0228f, 6.4409f, 8.3104f, //1511-1520
250 8.052f, 5.2309f, 5.8502f, 4.1374f, 4.7644f, 6.4343f, 3.9508f, 5.7653f, 4.9317f, 5.2804f, //1521-1530
251 3.9443f, 5.0276f, 4.623f, 4.0718f, 2.3652f, 2.867f, 0.0f, 0.0f, 0.0f, 0.0f, //1531-1540
252 0.0f, 0.0f, 0.0f, 0.0f, 7.718f, 0.0f, 17.663f, 14.185f, 14.391f, 6.4758f, //1541-1550
253 21.268f, 9.255f, 16.478f, 9.7134f, 17.252f, 5.2166f, 6.2833f, 7.179f, 7.1596f, 5.479f, //1551-1560
254 5.9853f, 1.7885f, 8.5821f, 3.0101f, 2.7443f, 1.5393f, 4.35f, 1.2402f, 4.1525f, 2.1657f, //1561-1570
255 3.677f, 2.226f, 0.8939f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.6427f, 5.8331f, //1571-1580
256 11.735f, 15.066f, 16.603f, 9.8191f, 17.839f, 18.905f, 25.065f, 28.556f, 38.188f, 22.002f, //1581-1590
257 17.039f, 8.0346f, 8.6258f, 7.6612f, 8.8731f, 3.7014f, 5.6601f, 2.5231f, 6.8799f, 2.8251f, //1591-1600
258 3.9736f, 4.5905f, 4.69f, 3.47f, 5.3112f, 3.634f, 4.3673f, 2.5818f, 2.9249f, 0.0f, //1601-1610
259 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 12.279f, 17.151f, 8.429f, 9.552f, //1611-1620
260 7.858f, 12.411f, 12.757f, 15.734f, 16.356f, 13.656f, 10.745f, 7.8242f, 7.434f, 4.617f, //1621-1630
261 4.8626f, 4.718f, 8.8805f, 1.1757f, 1.3974f,0.13352f, 1.6088f,0.50295f, 2.2496f,0.76546f, //1631-1640
262 3.1859f, 0.897f, 2.5919f, 3.0022f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1641-1650
263 2.6538f, 15.765f, 13.081f, 16.107f, 21.847f, 19.906f, 21.007f, 21.889f, 27.517f, 23.972f, //1651-1660
264 21.439f, 6.508f, 26.69f, 9.3653f, 12.159f, 11.788f, 8.946f, 6.081f, 7.364f, 5.2072f, //1661-1670
265 3.8609f, 1.1665f, 3.2796f, 3.2208f, 4.001f, 2.9363f, 3.3993f, 0.0f, 2.737f, 0.0f, //1671-1680
266 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 6.537f, 4.9648f, 9.8777f, //1681-1690
267 7.4389f, 6.8243f, 12.955f, 6.2838f, 6.8641f, 10.592f, 5.0213f, 7.3537f, 7.0303f, 10.535f, //1691-1700
268 8.5958f, 7.2168f, 6.8764f, 7.9774f, 0.8006f, 2.1223f, 0.3693f, 2.434f, 0.6599f, 2.1725f, //1701-1710
269 1.69f, 2.7706f, 0.6187f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1711-1720
270 0.0f, 0.122f, 0.0f, 4.245f, 7.0015f, 2.594f, 12.599f, 19.618f, 9.54f, 13.31f, //1721-1730
271 11.115f, 8.1008f, 11.238f, 23.37f, 12.248f, 12.739f, 15.283f, 6.1768f, 3.4678f, 11.525f, //1731-1740
272 5.802f, 9.7369f, 4.574f, 4.7496f, 3.4675f, 6.4652f, 1.3621f, 3.335f, 1.6945f, 2.3349f, //1741-1750
273 3.142f, 2.2584f, 0.0f, 0.0f, 0.0f, 1.4855f, 0.0f, 0.0f, 0.0f, 0.0f, //1751-1760
274 0.0f, 0.0f, 0.0f, 0.0f, 4.645f, 6.6324f, 7.2656f, 7.3347f, 9.0354f, 16.899f, //1761-1770
275 26.393f, 18.711f, 11.234f, 3.4709f, 21.43f, 18.293f, 7.6842f, 6.1919f, 5.7769f, 2.5922f, //1771-1780
276 4.9286f, 0.952f, 2.1987f,0.51329f, 1.4631f, 0.0f, 1.0518f, 0.0f, 0.0f, 0.0f, //1781-1790
277 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 7.443f, //1791-1800
278 10.498f, 6.4229f, 13.32f, 14.58f, 14.532f, 22.262f, 19.34f, 17.885f, 22.1f, 22.449f, //1801-1810
279 35.996f, 23.397f, 24.919f, 18.722f, 7.0352f, 7.3531f, 6.2406f, 6.5826f, 3.842f, 3.9233f, //1811-1820
280 1.5209f, 3.55f, 1.49f, 2.8571f, 0.0f, 2.3769f, 0.0f, 1.6015f, 0.0f, 0.0f, //1821-1830
281 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.0429f, 6.712f, 6.7863f, 10.126f, //1831-1840
282 5.3799f, 14.119f, 2.1568f, 14.344f, 14.091f, 19.595f, 21.805f, 35.544f, 18.003f, 7.447f, //1841-1850
283 7.45f, 17.07f, 15.907f, 6.5819f, 0.3884f, 3.794f,0.76797f, 1.857f,0.44211f, 0.373f, //1851-1860
284 0.0f, 0.0f, 1.0155f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.1661f, //1861-1870
285 0.0f, 3.2412f, 8.4975f, 4.6021f, 2.7442f, 6.095f, 3.6506f, 14.235f, 14.233f, 16.291f, //1871-1880
286 21.825f, 34.66f, 21.565f, 22.435f, 18.45f, 30.241f, 12.72f, 7.454f, 13.897f, 7.2306f, //1881-1890
287 5.799f, 7.276f, 5.0034f, 4.2123f, 3.9791f, 3.119f, 0.0978f, 1.315f, 0.0f, 0.235f, //1891-1900
288 0.0f, 0.0f, 0.0f, 0.0f, 4.0846f, 0.0f, 0.318f, 4.6911f, 2.6816f, 0.0f, //1901-1910
289 8.6613f, 3.7894f, 5.0987f, 10.753f, 9.9587f, 12.939f, 18.466f, 5.435f, 8.5161f, 10.248f, //1911-1920
290 15.876f, 13.422f, 5.331f, 5.77f, 5.949f, 5.882f, 4.3429f, 0.581f, 4.1095f, 1.1901f, //1921-1930
291 1.919f,0.63041f, 1.5782f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //1931-1940
292 0.33f, 0.0f, 0.0f, 5.008f, 9.5907f, 3.2427f, 9.5086f, 10.287f, 14.945f, 8.4096f, //1941-1950
293 16.03f, 17.086f, 16.583f, 19.505f, 24.115f, 17.888f, 22.788f, 26.143f, 27.581f, 13.865f, //1951-1960
294 11.549f, 1.6218f, 3.498f, 3.7561f, 3.921f, 4.549f, 7.5314f, 2.1041f, 1.3967f, 0.0f, //1961-1970
295 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.3115f, 0.0f, 3.5235f, 0.0f, 0.0f, //1971-1980
296 0.6713f, 6.9085f, 2.555f, 9.1688f, 6.1412f, 0.6142f, 4.773f, 13.242f, 10.357f, 8.1322f, //1981-1990
297 8.635f, 11.582f, 15.865f, 14.882f, 5.8891f, 10.051f, 9.3388f, 2.799f, 4.951f, 5.5983f, //1991-2000
298 3.8641f, 1.4006f,0.79742f, 1.2129f, 0.7729f, 1.3677f, 0.0f, 0.0f, 0.0f, 0.0f, //2001-2010
299 0.0f, 0.0f, 2.448f, 2.7445f, 2.5782f, 7.6099f, 2.2813f, 10.232f, 6.9884f, 13.96f, //2011-2020
300 11.971f, 15.403f, 15.709f, 12.392f, 6.081f, 12.931f, 11.848f, 12.716f, 16.767f, 16.846f, //2021-2030
301 10.961f, 5.0848f, 4.8797f, 3.863f, 2.018f, 4.61f, 4.425f, 4.729f, 2.9788f, 2.77f, //2031-2040
302 0.0f, 0.0f, 0.0f, 0.0f, 3.9984f, 0.0f, 2.6329f, 2.4313f, 4.528f, 5.2867f, //2041-2050
303 8.5435f, 0.0f, 6.4758f, 5.1397f, 15.943f, 9.3446f, 18.436f, 19.112f, 16.462f, 4.6733f, //2051-2060
304 16.821f, 6.293f, 6.9628f, 8.3751f, 10.907f, 0.602f, 1.8369f, 1.8828f, 3.524f, 3.0213f, //2061-2070
305 2.4978f, 0.3908f, 1.5862f, 0.9816f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.671f, //2071-2080
306 0.0f, 5.0193f, 6.4994f, 9.1066f, 7.0098f, 7.7475f, 11.599f, 14.178f, 11.972f, 11.219f, //2081-2090
307 9.9276f, 16.006f, 11.941f, 19.176f, 15.295f, 19.1f, 14.368f, 12.644f, 15.62f, 16.46f, //2091-2100
308 26.178f, 4.8635f, 5.0647f, 5.3886f, 4.2047f, 5.0575f, 3.0969f, 4.7339f, 1.6048f, 1.2722f, //2101-2110
309 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.9116f, //2111-2120
310 0.0f, 6.764f, 0.0f, 7.2134f, 3.3507f, 6.1682f, 3.9721f, 14.483f, 5.177f, 11.993f, //2121-2130
311 4.9183f, 9.817f, 4.7857f, 4.7272f, 7.5245f, 5.6101f, 6.562f, 6.4373f, 5.8217f, 5.9537f, //2131-2140
312 4.1868f, 6.4177f, 6.0629f, 1.7843f, 0.6172f, 1.2585f, 0.0f, 0.0f, 0.0f, 0.0f, //2141-2150
313 0.175f, 0.0f, 0.0f, 0.0f, 0.0f, 4.8611f, 0.0f, 0.45f, 7.737f, 9.2967f, //2151-2160
314 7.4704f, 10.076f, 9.662f, 11.128f, 7.224f, 11.37f, 9.7766f, 9.0844f, 3.8879f, 4.6062f, //2161-2170
315 4.0387f, 6.9217f, 7.2047f, 13.394f, 6.7925f, 7.632f, 8.655f, 5.6669f, 3.6661f, 6.7601f, //2171-2180
316 3.2769f, 6.417f, 3.3438f, 1.5382f, 0.0f, 2.381f, 0.0f, 0.219f, 0.0f, 0.0f, //2181-2190
317 0.0f, 0.0f, 4.4008f, 0.0f, 0.0f, 0.0f, 0.0f, 0.175f, 0.092f, 0.3492f, //2191-2200
318 4.7279f, 0.3701f, 7.3785f, 0.9063f, 7.566f, 6.042f, 12.733f, 7.8227f, 10.824f, 5.9987f, //2201-2210
319 5.408f, 4.8181f, 7.5907f, 5.858f, 2.3926f, 2.9915f, 4.799f,0.99784f, 2.1996f,0.75766f, //2211-2220
320 0.4824f, 0.3197f, 0.0f, 0.267f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2221-2230
321 0.0f, 2.1895f, 0.0f, 2.8393f, 2.6094f, 3.9108f, 3.9845f, 5.1584f, 3.625f, 6.7068f, //2231-2240
322 6.2602f, 8.6901f, 7.541f, 10.545f, 7.0075f, 6.683f, 5.612f, 13.518f, 11.02f, 11.147f, //2241-2250
323 12.741f, 13.482f, 4.1802f, 10.671f, 7.3585f, 7.7774f, 1.2104f, 4.5075f, 1.4723f, 5.1304f, //2251-2260
324 2.64f, 4.1138f, 1.4595f, 2.541f, 1.4647f, 0.0f, 0.0f, 3.1985f, 0.0f, 0.0f, //2261-2270
325 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.072f, 3.4415f, 2.0874f, 3.6317f, //2271-2280
326 0.1684f, 6.254f, 0.0f, 6.076f, 7.14f, 9.1252f, 9.2007f, 6.8299f, 5.0893f, 7.0521f, //2281-2290
327 3.8694f, 6.987f, 3.1382f, 5.0544f, 4.7856f, 5.497f, 0.9235f, 2.3332f, 1.181f, 2.3121f, //2291-2300
328 6.1981f, 2.404f,0.54259f,0.94375f,0.52237f, 0.0f, 0.1164f, 0.0f, 0.0f, 0.0f, //2301-2310
329 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.0332f, 0.545f, 2.6366f, 3.5166f, 4.218f, //2311-2320
330 4.085f, 3.449f, 4.231f, 6.8786f, 4.6838f, 6.6012f, 8.1748f, 7.434f, 7.0384f, 7.3969f, //2321-2330
331 5.2578f, 7.7332f, 4.9136f, 7.4078f, 5.5103f, 6.5531f, 2.8648f, 4.2143f, 6.1489f, 4.9507f, //2331-2340
332 3.1292f, 4.5298f, 2.5926f, 3.4243f,0.74782f, 3.017f, 0.424f, 3.2266f, 1.4555f, 1.7885f, //2341-2350
333 0.0f, 3.153f, 0.0f, 0.0f, 0.0f, 1.8274f, 2.1985f, 0.0f, 5.4301f, 1.7992f, //2351-2360
334 6.1583f, 0.0f, 7.1722f, 0.0f, 5.7177f, 4.5791f, 7.8795f, 6.3223f, 7.038f, 3.954f, //2361-2370
335 7.22f, 5.3252f, 9.5804f, 5.3789f, 10.935f, 4.6355f, 4.0632f, 2.1848f, 2.792f, 1.673f, //2371-2380
336 1.242f, 1.5604f, 1.7705f,0.46871f, 1.76f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2381-2390
337 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.2786f, 0.0f, 0.0f, 2.523f, 2.8744f, //2391-2400
338 1.7955f, 5.5345f, 5.2669f, 5.8034f, 2.9527f, 5.1081f, 4.7019f, 7.3176f, 4.7027f, 8.873f, //2401-2410
339 4.8702f, 7.9412f, 6.8197f, 10.031f, 8.7932f, 11.427f, 11.406f, 18.664f, 11.574f, 8.6523f, //2411-2420
340 4.0612f, 5.2842f, 3.0684f, 6.1621f, 1.1878f, 4.922f, 1.4883f, 3.6893f, 2.5795f, 6.2766f, //2421-2430
341 0.0f, 1.2969f, 0.0f, 1.366f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2431-2440
342 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.657f, 0.52f, 3.915f, 0.3198f, 4.6724f, //2441-2450
343 0.9716f, 5.608f, 0.8768f, 6.3331f, 1.8545f, 6.9284f, 7.4519f, 13.772f, 6.815f, 8.8743f, //2451-2460
344 3.6298f, 3.9657f, 4.0664f, 2.6432f, 2.5481f, 2.7474f, 2.0448f, 6.4188f, 6.6548f, 8.4086f, //2461-2470
345 2.6431f, 4.418f, 1.728f, 0.8683f, 0.2536f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2471-2480
346 0.0f, 0.0f, 0.0f, 0.0f, 2.138f, 0.0f, 2.8122f, 0.0f, 1.6634f, 0.278f, //2481-2490
347 3.3139f, 1.723f, 4.1665f, 3.1024f, 5.3303f, 5.142f, 9.6693f, 7.932f, 12.37f, 7.538f, //2491-2500
348 13.516f, 10.406f, 12.699f, 10.66f, 8.6305f, 8.6548f, 9.2054f, 5.5718f, 8.1261f, 6.8658f, //2501-2510
349 8.7056f, 7.306f, 13.675f, 4.142f, 4.133f, 5.5574f, 0.8049f, 0.0f, 1.365f, 0.0f, //2511-2520
350 1.459f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.063f, 0.0f, //2521-2530
351 2.973f, 0.273f, 1.7366f, 0.0f, 3.401f, 0.396f, 2.9226f, 5.4847f, 4.0248f, 3.9636f, //2531-2540
352 4.2969f, 3.5777f, 5.2835f, 7.2753f, 5.8247f, 3.809f, 6.7188f, 4.3053f, 4.5599f, 5.6266f, //2541-2550
353 7.632f, 5.996f, 6.191f, 0.4945f, 0.7589f,0.88804f, 1.9598f, 0.0f, 1.436f, 0.0f, //2551-2560
354 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.194f, //2561-2570
355 3.0402f, 1.671f, 2.2946f, 1.1297f, 3.3265f, 2.0207f, 3.6467f, 1.6862f, 5.1133f, 3.931f, //2571-2580
356 5.2269f, 3.7108f, 5.5289f, 4.3583f, 5.9111f, 4.6288f, 9.7247f, 4.3463f, 9.6437f, 4.3542f, //2581-2590
357 10.084f, 4.911f, 2.922f, 1.3282f, 3.2624f, 2.1595f, 2.2338f, 0.2647f, 2.2868f, 0.0f, //2591-2600
358 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2601-2610
359 0.0f, 0.0f, 1.636f, 0.1579f, 0.324f, 0.0f, 0.573f, 0.344f, 2.0774f, 0.3917f, //2611-2620
360 2.7011f, 2.3462f, 4.546f, 2.9023f, 2.1495f, 4.3669f, 4.695f, 9.0673f, 10.016f, 4.7716f, //2621-2630
361 2.926f, 0.7284f, 0.517f, 0.4793f, 0.8919f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2631-2640
362 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2641-2650
363 0.0f, 0.0f, 2.3838f, 0.0f, 2.7762f, 0.0f, 2.7755f, 4.7208f, 4.5832f, 4.0594f, //2651-2660
364 4.1298f, 3.013f, 5.9306f, 7.3077f, 12.026f, 11.232f, 10.612f, 8.8319f, 5.051f, 2.2871f, //2661-2670
365 4.2997f, 0.6189f, 5.1678f, 0.7328f, 3.9615f, 0.0f, 3.6961f, 0.0f, 0.0f, 0.0f, //2671-2680
366 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2681-2690
367 0.0f, 0.531f, 1.8852f, 0.827f, 0.0f, 5.292f, 8.23f, 8.5334f, 8.0948f, 8.658f, //2691-2700
368 3.4623f, 1.973f, 3.0023f, 3.0451f, 0.778f, 0.5111f, 0.8242f, 0.0f, 1.6954f, 0.0f, //2701-2710
369 1.7498f, 0.0f,0.95506f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2711-2720
370 0.0f, 0.0f, 0.0f, 0.0f, 2.0101f, 0.0f, 2.1474f, 0.0f, 4.257f, 1.1981f, //2721-2730
371 3.9294f, 1.77f, 7.87f, 6.2832f, 6.2661f, 5.0f, 6.6788f, 4.3458f, 6.2555f, 1.3758f, //2731-2740
372 3.2877f, 1.0289f, 2.3687f, 0.7241f, 2.382f, 1.4748f, 3.4189f, 1.7205f, 2.4183f, 1.7739f, //2741-2750
373 1.05f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.1336f, //2751-2760
374 1.884f, 0.0918f, 2.438f, 0.322f, 2.013f, 2.5077f, 3.7208f, 3.4056f, 2.5144f, 0.0f, //2761-2770
375 0.1675f, 0.452f,0.25704f, 0.0f, 0.8747f,0.03307f, 0.3355f,0.63577f, 2.9997f, 0.4785f, //2771-2780
376 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 2.181f, //2781-2790
377 1.4213f, 3.6814f, 0.6738f, 2.104f, 2.322f, 4.892f, 2.4215f, 5.0979f, 1.7025f, 2.864f, //2791-2800
378 2.4944f, 3.0971f, 2.1663f, 3.2834f,0.77929f, 3.819f, 1.173f, 5.164f, 4.512f, 3.7751f, //2801-2810
379 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2811-2820
380 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2821-2830
381 0.0f, 0.8408f, 0.0f, 2.14f, 0.0f, 4.0515f,0.18673f,0.75575f,0.67811f, 0.0f, //2831-2840
382 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2841-2850
383 0.0f, 1.847f, 0.0f, 0.0f, 0.0f, 3.243f, 1.2681f, 2.6598f, 3.174f, 4.807f, //2851-2860
384 4.0407f, 5.077f, 4.835f, 5.513f, 3.107f, 1.2761f, 0.0f, 0.0f, 0.0f, 0.0f, //2861-2870
385 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.2433f, 0.0f, 3.401f, 0.0f, //2871-2880
386 4.6146f,0.64676f, 1.1971f, 0.2998f, 0.0f, 0.6406f, 0.0f, 0.0f, 0.0f, 0.0f, //2881-2890
387 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.1188f, 1.786f, 1.0006f, 4.2652f, 4.413f, //2891-2900
388 6.0963f, 1.3577f, 4.172f, 5.0342f, 6.123f, 0.766f, 1.1408f, 0.0f, 0.0f, 0.0f, //2901-2910
389 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.557f, 0.0f, 5.9809f, //2911-2920
390 2.2f, 0.5324f, 0.8789f, 1.1856f,0.29937f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2921-2930
391 0.0f, 0.0f, 0.0f, 0.0f, 0.038f, 0.4202f, 0.288f, 0.0874f, 1.1059f, 1.1023f, //2931-2940
392 4.0332f, 0.5817f, 5.6807f, 4.7135f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2941-2950
393 0.0f, 0.0f, 0.0f, 0.24f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2951-2960
394 0.4894f, 0.59f, 1.2275f, 0.471f, 0.9784f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //2961-2970
395 0.0f, 0.0f, 0.0f, 0.315f, 0.0f, 0.055f, 0.0f, 0.738f, 0.7373f, 3.315f, //2971-2980
396 1.6952f, 1.25f, 0.9698f, 0.417f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.2096f, //2981-2990
397 0.0f, 0.2532f, 0.0f, 0.0f, 0.0f, 0.3532f, 0.0f, 1.357f, 0.6f, 0.435f, //2991-3000
398 0.4692f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3001-3010
399 0.0f, 2.372f, 0.669f, 0.0f,0.55843f, 0.0f, 0.398f,0.73365f, 0.2311f, 1.5602f, //3011-3020
400 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 3.1278f, //3021-3030
401 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3031-3040
402 0.0f, 0.0f, 1.6997f, 2.8792f, 3.465f, 2.839f, 0.0f, 0.0f, 0.0f, 0.0f, //3041-3050
403 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 1.8999f, 0.0f, 0.0f, //3051-3060
404 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3061-3070
405 0.0f, 0.656f, 2.868f, 1.157f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3071-3080
406 0.0f, 0.0f, 0.0f, 0.0f, 0.6f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3081-3090
407 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3091-3100
408 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3101-3110
409 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3111-3120
410 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3121-3130
411 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3131-3140
412 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3141-3150
413 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3151-3160
414 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3161-3170
415 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, //3171-3180
416 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f};
417
418namespace
419{
420 G4Mutex nuclearLevelDataMutex = G4MUTEX_INITIALIZER;
421}
422
424{
425 if (nullptr == theInstance) {
426 G4AutoLock l(&nuclearLevelDataMutex);
427 if (nullptr == theInstance) {
428 static G4NuclearLevelData theData;
429 theInstance = &theData;
430 }
431 l.unlock();
432 }
433 return theInstance;
434}
435
436G4NuclearLevelData::G4NuclearLevelData()
437{
438 fDeexPrecoParameters = new G4DeexPrecoParameters();
439 fLevelReader = new G4LevelReader(this);
440 for(G4int Z=0; Z<ZMAX; ++Z) {
441 (fLevelManagers[Z]).resize(AMAX[Z]-AMIN[Z]+1,nullptr);
442 (fLevelManagerFlags[Z]).resize(AMAX[Z]-AMIN[Z]+1,false);
443 }
444 fShellCorrection = new G4ShellCorrection();
445 fPairingCorrection = new G4PairingCorrection();
446 fG4calc = G4Pow::GetInstance();
447}
448
450{
451 delete fLevelReader;
452 delete fDeexPrecoParameters;
453 delete fShellCorrection;
454 delete fPairingCorrection;
455 for(G4int Z=1; Z<ZMAX; ++Z) {
456 size_t nn = (fLevelManagers[Z]).size();
457 for(size_t j=0; j<nn; ++j) {
458 delete (fLevelManagers[Z])[j];
459 }
460 }
461}
462
463const G4LevelManager*
465{
466 if(Z < 1 || Z >= ZMAX || A < AMIN[Z] || A > AMAX[Z]) { return nullptr; }
467 const G4int idx = A - AMIN[Z];
468 if( !(fLevelManagerFlags[Z])[idx] ) {
469 G4AutoLock l(&nuclearLevelDataMutex);
470 if( !(fLevelManagerFlags[Z])[idx] ) {
471 (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
472 (fLevelManagerFlags[Z])[idx] = true;
473 }
474 l.unlock();
475 }
476 return (fLevelManagers[Z])[idx];
477}
478
479G4bool
481{
482 G4bool res = false;
483 if(Z > 0 && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) {
484 G4AutoLock l(&nuclearLevelDataMutex);
485 const G4LevelManager* newman = fLevelReader->MakeLevelManager(Z, A, fname);
486 // if file is corrupted G4LevelReader should make the exception
487 if(newman != nullptr) {
488 res = true;
489 if(0 < fDeexPrecoParameters->GetVerbose()) {
490 G4cout << "G4NuclearLevelData::AddPrivateData for Z= " << Z
491 << " A= " << A << " from <" << fname
492 << "> is done" << G4endl;
493 }
494 const G4int idx = A - AMIN[Z];
495 delete (fLevelManagers[Z])[idx];
496 (fLevelManagers[Z])[idx] = newman;
497 (fLevelManagerFlags[Z])[idx] = true;
498 }
499 l.unlock();
500 } else {
502 ed << "private nuclear level data for Z= " << Z << " A= " << A
503 << " outside allowed limits ";
504 G4Exception("G4NuclearLevelData::AddPrivateData","had0433",FatalException,
505 ed,"Stop execution");
506 }
507 return res;
508}
509
511{
512 return (Z >= 0 && Z < ZMAX) ? AMIN[Z] : 0;
513}
514
516{
517 return (Z >= 0 && Z < ZMAX) ? AMAX[Z] : 0;
518}
519
521{
522 if(fInitialized) return;
523 G4AutoLock l(&nuclearLevelDataMutex);
524 if(!fInitialized) {
525 fInitialized = true;
526 G4int mZ = Zlim + 1;
527 if(mZ > ZMAX) { mZ = ZMAX; }
528 for(G4int Z=1; Z<mZ; ++Z) {
529 for(G4int A=AMIN[Z]; A<=AMAX[Z]; ++A) {
530 G4int idx = A - AMIN[Z];
531 if( !(fLevelManagerFlags[Z])[idx] ) {
532 (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
533 (fLevelManagerFlags[Z])[idx] = true;
534 }
535 }
536 }
537 }
538 l.unlock();
539}
540
542{
543 return (0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) ?
544 (G4double)(LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0.0;
545}
546
548{
549 return (0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) ?
550 (LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0.0f;
551}
552
554{
555 G4double e = energy;
556 if(e <= GetMaxLevelEnergy(Z, A)) {
557 const G4LevelManager* man = GetLevelManager(Z, A);
558 if(man != nullptr) {
559 e = (G4double)man->NearestLevelEnergy(e, man->NumberOfTransitions());
560 }
561 }
562 return e;
563}
564
567{
569 if(energy < e) {
570 const G4LevelManager* man = GetLevelManager(Z, A);
571 if(man != nullptr) {
572 e = (G4double)man->NearestLowEdgeLevelEnergy(energy);
573 }
574 }
575 return e;
576}
577
580 G4double Mass, G4double partMass, G4double ekin)
581{
582 G4double T = ekin;
583 G4double E2 = (Mass - partMass)*(Mass - partMass);
584 G4double Eex = std::sqrt(E2 - 2.0*Mass*T) - resMass;
585 if(Eex <= GetMaxLevelEnergy(Z, A)) {
586 if(Eex <= 0.0) {
587 Eex = 0.0;
588 T = (E2 - resMass*resMass)*0.5/Mass;
589 } else {
590 const G4LevelManager* man = GetLevelManager(Z, A);
591 if(man != nullptr) {
592 size_t idx = man->NearestLevelIndex(Eex, man->NumberOfTransitions());
593 for(;;) {
594 Eex = (G4double)man->LevelEnergy(idx);
595 T = (E2 - (resMass+Eex)*(resMass+Eex))*0.5/Mass;
596 if(T >= 0.0 || 0 == idx) { break; }
597 --idx;
598 }
599 }
600 }
601 T = std::max(T, 0.0);
602 }
603 return T;
604}
605
607{
608 return fDeexPrecoParameters;
609}
610
612{
613 return fPairingCorrection;
614}
615
617{
618 return fShellCorrection;
619}
620
622{
623 if(fDeexPrecoParameters->GetLevelDensityFlag()) {
624 return A*fDeexPrecoParameters->GetLevelDensity();
625 }
626 const G4LevelManager* man = GetLevelManager(Z, A);
627 return (man) ? man->LevelDensity(U)
628 : 0.058025*A*(1.0 + 5.9059/fG4calc->Z13(A));
629}
630
632{
633 if(fDeexPrecoParameters->GetLevelDensityFlag()) {
634 return fPairingCorrection->GetPairingCorrection(A, Z);
635 }
636 G4int N = A - Z;
637 const G4double par = 12.*CLHEP::MeV;
638 G4double x = (A <= 36) ? 6.0 : std::sqrt((G4double)A);
639 return (2 - Z + (Z/2)*2 - N + (N/2)*2)*par/x;
640}
641
642void G4NuclearLevelData::StreamLevels(std::ostream& os,
643 G4int Z, G4int A)
644{
645 const G4LevelManager* man = GetLevelManager(Z, A);
646 if(man) {
647 os << "Level data for Z= " << Z << " A= " << A << " "
648 << man->NumberOfTransitions() + 1 << " levels \n";
649 man->StreamInfo(os);
650 }
651}
@ FatalException
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
std::ostringstream G4ExceptionDescription
#define G4MUTEX_INITIALIZER
std::mutex G4Mutex
float G4float
Definition G4Types.hh:84
double G4double
Definition G4Types.hh:83
bool G4bool
Definition G4Types.hh:86
int G4int
Definition G4Types.hh:85
const G4double A[17]
#define G4endl
Definition G4ios.hh:67
G4GLOB_DLL std::ostream G4cout
void StreamInfo(std::ostream &os) const
G4double NearestLowEdgeLevelEnergy(const G4double energy) const
std::size_t NearestLevelIndex(const G4double energy, const std::size_t index=0) const
G4double LevelEnergy(const std::size_t i) const
G4double LevelDensity(const G4double U) const
std::size_t NumberOfTransitions() const
G4double NearestLevelEnergy(const G4double energy, const std::size_t index=0) const
const G4LevelManager * MakeLevelManager(G4int Z, G4int A, const G4String &filename)
const G4LevelManager * CreateLevelManager(G4int Z, G4int A)
G4float MaxLevelEnergy(G4int Z, G4int A) const
G4double GetLevelDensity(G4int Z, G4int A, G4double U)
G4double GetLevelEnergy(G4int Z, G4int A, G4double energy)
G4double FindLevel(G4int Z, G4int A, G4double resMass, G4double Mass, G4double partMass, G4double T)
G4DeexPrecoParameters * GetParameters()
G4bool AddPrivateData(G4int Z, G4int A, const G4String &filename)
G4int GetMinA(G4int Z) const
void UploadNuclearLevelData(G4int Zlim)
G4double GetLowEdgeLevelEnergy(G4int Z, G4int A, G4double energy)
G4ShellCorrection * GetShellCorrection()
const G4LevelManager * GetLevelManager(G4int Z, G4int A)
G4PairingCorrection * GetPairingCorrection()
static G4NuclearLevelData * GetInstance()
G4double GetMaxLevelEnergy(G4int Z, G4int A) const
void StreamLevels(std::ostream &os, G4int Z, G4int A)
G4int GetMaxA(G4int Z) const
G4double GetPairingCorrection(G4int A, G4int Z) const
static G4Pow * GetInstance()
Definition G4Pow.cc:41
G4double Z13(G4int Z) const
Definition G4Pow.hh:123
#define N
Definition crc32.c:57