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