| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 // ------------------------------------------- 27 // -------------------------------------------------------------------
28 // 28 //
29 // GEANT4 source file 29 // GEANT4 source file
30 // 30 //
31 // File name: G4NuclearLevelData 31 // File name: G4NuclearLevelData
32 // 32 //
33 // Author: V.Ivanchenko 33 // Author: V.Ivanchenko
34 // 34 //
35 // Creation date: 10 February 2015 35 // Creation date: 10 February 2015
36 // 36 //
37 // Modifications: 37 // Modifications:
38 // 38 //
39 // ------------------------------------------- 39 // -------------------------------------------------------------------
40 40
41 #include "G4NuclearLevelData.hh" 41 #include "G4NuclearLevelData.hh"
42 #include "G4LevelReader.hh" 42 #include "G4LevelReader.hh"
43 #include "G4LevelManager.hh" 43 #include "G4LevelManager.hh"
44 #include "G4Element.hh" 44 #include "G4Element.hh"
45 #include "G4ElementTable.hh" 45 #include "G4ElementTable.hh"
46 #include "G4DeexPrecoParameters.hh" 46 #include "G4DeexPrecoParameters.hh"
47 #include "G4PairingCorrection.hh" 47 #include "G4PairingCorrection.hh"
48 #include "G4ShellCorrection.hh" 48 #include "G4ShellCorrection.hh"
49 #include "G4SystemOfUnits.hh" 49 #include "G4SystemOfUnits.hh"
50 #include "G4AutoLock.hh" <<
51 #include "G4Pow.hh" 50 #include "G4Pow.hh"
52 #include <iomanip> 51 #include <iomanip>
53 52
54 G4NuclearLevelData* G4NuclearLevelData::theIns 53 G4NuclearLevelData* G4NuclearLevelData::theInstance = nullptr;
55 54
56 const G4int G4NuclearLevelData::AMIN[] = {1, 55 const G4int G4NuclearLevelData::AMIN[] = {1,
57 1,3,4,5,7,8,10,12,14,16, // 56 1,3,4,5,7,8,10,12,14,16, // Z= 1-10
58 20,19,21,22,25,27,29,31,35,35, // 57 20,19,21,22,25,27,29,31,35,35, // Z= 11-20
59 39,39,42,42,44,45,50,48,53,54, // 58 39,39,42,42,44,45,50,48,53,54, // Z= 21-30
60 60,60,63,64,69,69,73,73,76,78, // 59 60,60,63,64,69,69,73,73,76,78, // Z= 31-40
61 81,83,85,88,90,91,93,96,98,100, // 60 81,83,85,88,90,91,93,96,98,100, // Z= 41-50
62 104,105,108,109,112,114,117,121,121,125, // 61 104,105,108,109,112,114,117,121,121,125, // Z= 51-60
63 128,129,130,134,135,139,140,143,144,149, // 62 128,129,130,134,135,139,140,143,144,149, // Z= 61-70
64 150,153,155,158,159,161,164,166,170,171, // 63 150,153,155,158,159,161,164,166,170,171, // Z= 71-80
65 176,178,184,186,191,193,199,201,206,208, // 64 176,178,184,186,191,193,199,201,206,208, // Z= 81-90
66 212,217,225,228,230,233,234,237,241,241, // 65 212,217,225,228,230,233,234,237,241,241, // Z= 91-100
67 245,248,252,253,255,258,260,263,266,269, // 66 245,248,252,253,255,258,260,263,266,269, // Z=101-110
68 272,277,282,0,289,0,293}; // 67 272,277,282,0,289,0,293}; // Z=111-117
69 const G4int G4NuclearLevelData::AMAX[] = {1, 68 const G4int G4NuclearLevelData::AMAX[] = {1,
70 6,10,11,14,17,22,24,26,30,30, // 69 6,10,11,14,17,22,24,26,30,30, // Z= 1-10
71 37,40,42,44,46,48,51,53,56,58, // 70 37,40,42,44,46,48,51,53,56,58, // Z= 11-20
72 61,63,66,68,71,74,76,79,82,85, // 71 61,63,66,68,71,74,76,79,82,85, // Z= 21-30
73 87,90,92,95,98,101,103,107,109,112, // 72 87,90,92,95,98,101,103,107,109,112, // Z= 31-40
74 114,117,120,124,126,128,129,133,135,138, // 73 114,117,120,124,126,128,129,133,135,138, // Z= 41-50
75 140,143,145,148,151,152,153,155,155,161, // 74 140,143,145,148,151,152,153,155,155,161, // Z= 51-60
76 163,165,167,169,171,173,175,177,179,181, // 75 163,165,167,169,171,173,175,177,179,181, // Z= 61-70
77 184,189,192,194,198,202,204,205,210,216, // 76 184,189,192,194,198,202,204,205,210,216, // Z= 71-80,
78 217,220,224,227,229,229,233,234,235,238, // 77 217,220,224,227,229,229,233,234,235,238, // Z= 81-90,
79 239,242,244,247,248,252,253,256,257,259, // 78 239,242,244,247,248,252,253,256,257,259, // Z= 91-100
80 260,262,266,265,270,269,274,277,278,277, // 79 260,262,266,265,270,269,274,277,278,277, // Z=101-110
81 282,281,286,0,290,0,294}; // 80 282,281,286,0,290,0,294}; // Z=111-117
82 const G4int G4NuclearLevelData::LEVELIDX[] = { 81 const G4int G4NuclearLevelData::LEVELIDX[] = {0,
83 0, 6, 14, 22, 32, 43 82 0, 6, 14, 22, 32, 43, 58, 73, 88, 105, // Z= 1-10
84 120, 138, 160, 182, 205, 227 83 120, 138, 160, 182, 205, 227, 249, 272, 295, 317, // Z= 11-20
85 341, 364, 389, 414, 441, 469 84 341, 364, 389, 414, 441, 469, 499, 526, 558, 588, // Z= 21-30
86 620, 648, 679, 709, 741, 771 85 620, 648, 679, 709, 741, 771, 804, 835, 870, 904, // Z= 31-40
87 939, 973, 1008, 1044, 1081, 1118 86 939, 973, 1008, 1044, 1081, 1118, 1156, 1193, 1231, 1269, // Z= 41-50
88 1308, 1345, 1384, 1422, 1462, 1502 87 1308, 1345, 1384, 1422, 1462, 1502, 1541, 1578, 1613, 1648, // Z= 51-60
89 1685, 1721, 1758, 1796, 1832, 1869 88 1685, 1721, 1758, 1796, 1832, 1869, 1904, 1940, 1975, 2011, // Z= 61-70
90 2044, 2079, 2116, 2154, 2191, 2231 89 2044, 2079, 2116, 2154, 2191, 2231, 2273, 2314, 2354, 2395, // Z= 71-80
91 2441, 2483, 2526, 2567, 2609, 2648 90 2441, 2483, 2526, 2567, 2609, 2648, 2685, 2720, 2754, 2784, // Z= 81-90
92 2815, 2843, 2869, 2889, 2909, 2928 91 2815, 2843, 2869, 2889, 2909, 2928, 2948, 2968, 2988, 3005, // Z= 91-100
93 3024, 3040, 3055, 3070, 3083, 3099 92 3024, 3040, 3055, 3070, 3083, 3099, 3111, 3126, 3141, 3154, // Z= 101-110
94 3163, 3174, 3179, 3184, 3184, 3186 93 3163, 3174, 3179, 3184, 3184, 3186, 3186}; // Z= 111-117
95 94
96 // obtained from PhotonEvaporation5.2 95 // obtained from PhotonEvaporation5.2
97 static const G4float LEVELMAX[3188] = {0.0f, 96 static const G4float LEVELMAX[3188] = {0.0f,
98 0.0f, 0.0f, 0.0f, 0.0f, 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
99 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 98 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, 99 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, 100 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, 101 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, 102 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, 103 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, 104 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, 105 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, 106 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, 107 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, 108 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, 109 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, 110 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, 111 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, 112 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, 113 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, 114 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, 115 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, 116 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, 117 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, 118 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, 119 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, 120 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, 121 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, 122 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, 123 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, 124 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, 125 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, 126 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, 127 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, 128 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, 129 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, 130 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, 131 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, 132 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, 133 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, 134 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, 135 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, 136 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, 137 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, 138 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, 139 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, 140 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, 141 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, 142 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, 143 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, 144 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, 145 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, 146 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, 147 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, 148 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, 149 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, 150 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, 151 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, 152 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, 153 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, 154 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, 155 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, 156 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, 157 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, 158 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, 159 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, 160 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, 161 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, 162 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, 163 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, 164 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, 165 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, 166 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, 167 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, 168 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, 169 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, 170 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, 171 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, 172 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, 173 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, 174 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, 175 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, 176 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, 177 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, 178 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, 179 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, 180 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, 181 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, 182 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, 183 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, 184 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, 185 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, 186 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, 187 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, 188 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, 189 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, 190 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, 191 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, 192 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, 193 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, 194 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, 195 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, 196 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, 197 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, 198 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, 199 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, 200 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, 201 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, 202 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, 203 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, 204 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, 205 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, 206 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, 207 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, 208 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, 209 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, 210 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, 211 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, 212 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, 213 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 214 4.3106f, 5.4145f, 2.6661f, 2.3527f, 1.8189f,0.78316f, 1.8858f, 3.239f, 2.672f, 0.3239f, //1171-1180
216 0.72042f, 0.326f, 0.8058f, 0.0f, 0.0f, 215 0.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, 216 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, 217 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, 218 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, 219 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, 220 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, 221 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, 222 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, 223 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, 224 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, 225 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, 226 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, 227 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, 228 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 229 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, 230 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, 231 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, 232 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, 233 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, 234 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, 235 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, 236 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, 237 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, 238 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, 239 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, 240 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, 241 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, 242 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, 243 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, 244 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, 245 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, 246 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, 247 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, 248 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, 249 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, 250 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, 251 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, 252 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, 253 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, 254 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, 255 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, 256 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, 257 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, 258 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, 259 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 260 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, 261 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, 262 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, 263 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, 264 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, 265 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, 266 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, 267 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, 268 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, 269 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, 270 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, 271 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, 272 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, 273 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, 274 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, 275 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, 276 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, 277 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, 278 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, 279 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, 280 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, 281 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, 282 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, 283 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, 284 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, 285 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, 286 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, 287 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, 288 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, 289 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, 290 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, 291 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, 292 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, 293 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, 294 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, 295 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, 296 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, 297 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, 298 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, 299 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, 300 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, 301 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, 302 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, 303 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, 304 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, 305 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, 306 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, 307 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, 308 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, 309 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, 310 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, 311 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, 312 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, 313 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, 314 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, 315 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, 316 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, 317 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, 318 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, 319 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, 320 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, 321 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, 322 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, 323 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, 324 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, 325 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, 326 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, 327 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, 328 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, 329 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, 330 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, 331 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, 332 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, 333 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, 334 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, 335 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, 336 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, 337 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, 338 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, 339 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, 340 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, 341 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, 342 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, 343 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, 344 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, 345 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, 346 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, 347 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, 348 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, 349 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, 350 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, 351 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 352 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, 353 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, 354 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, 355 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, 356 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, 357 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, 358 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, 359 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, 360 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, 361 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, 362 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, 363 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, 364 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, 365 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, 366 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, 367 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, 368 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, 369 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, 370 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, 371 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, 372 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, 373 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 374 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, 375 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, 376 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, 377 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, 378 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, 379 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, 380 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, 381 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, 382 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, 383 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, 384 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, 385 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, 386 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, 387 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, 388 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, 389 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, 390 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, 391 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, 392 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, 393 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, 394 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, 395 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, 396 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, 397 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, 398 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, 399 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, 400 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, 401 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, 402 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, 403 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, 404 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, 405 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, 406 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, 407 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, 408 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, 409 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, 410 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, 411 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, 412 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, 413 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, 414 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, 415 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f};
417 416
418 namespace << 417 #ifdef G4MULTITHREADED
419 { << 418 G4Mutex G4NuclearLevelData::nuclearLevelDataMutex = G4MUTEX_INITIALIZER;
420 G4Mutex nuclearLevelDataMutex = G4MUTEX_INIT << 419 #endif
421 } <<
422 420
423 G4NuclearLevelData* G4NuclearLevelData::GetIns 421 G4NuclearLevelData* G4NuclearLevelData::GetInstance()
424 { 422 {
425 if (nullptr == theInstance) { << 423 if (!theInstance) {
426 G4AutoLock l(&nuclearLevelDataMutex); << 424 #ifdef G4MULTITHREADED
427 if (nullptr == theInstance) { << 425 G4MUTEXLOCK(&nuclearLevelDataMutex);
>> 426 if (!theInstance) {
>> 427 #endif
428 static G4NuclearLevelData theData; 428 static G4NuclearLevelData theData;
429 theInstance = &theData; 429 theInstance = &theData;
>> 430 #ifdef G4MULTITHREADED
430 } 431 }
431 l.unlock(); << 432 G4MUTEXUNLOCK(&nuclearLevelDataMutex);
>> 433 #endif
432 } 434 }
433 return theInstance; 435 return theInstance;
434 } 436 }
435 437
436 G4NuclearLevelData::G4NuclearLevelData() 438 G4NuclearLevelData::G4NuclearLevelData()
437 { 439 {
438 fDeexPrecoParameters = new G4DeexPrecoParame 440 fDeexPrecoParameters = new G4DeexPrecoParameters();
439 fLevelReader = new G4LevelReader(this); 441 fLevelReader = new G4LevelReader(this);
440 for(G4int Z=0; Z<ZMAX; ++Z) { 442 for(G4int Z=0; Z<ZMAX; ++Z) {
441 (fLevelManagers[Z]).resize(AMAX[Z]-AMIN[Z] 443 (fLevelManagers[Z]).resize(AMAX[Z]-AMIN[Z]+1,nullptr);
442 (fLevelManagerFlags[Z]).resize(AMAX[Z]-AMI 444 (fLevelManagerFlags[Z]).resize(AMAX[Z]-AMIN[Z]+1,false);
443 } 445 }
444 fShellCorrection = new G4ShellCorrection(); 446 fShellCorrection = new G4ShellCorrection();
445 fPairingCorrection = new G4PairingCorrection 447 fPairingCorrection = new G4PairingCorrection();
446 fG4calc = G4Pow::GetInstance(); 448 fG4calc = G4Pow::GetInstance();
>> 449 fInitialized = false;
447 } 450 }
448 451
449 G4NuclearLevelData::~G4NuclearLevelData() 452 G4NuclearLevelData::~G4NuclearLevelData()
450 { 453 {
451 delete fLevelReader; 454 delete fLevelReader;
452 delete fDeexPrecoParameters; 455 delete fDeexPrecoParameters;
453 delete fShellCorrection; 456 delete fShellCorrection;
454 delete fPairingCorrection; 457 delete fPairingCorrection;
455 for(G4int Z=1; Z<ZMAX; ++Z) { 458 for(G4int Z=1; Z<ZMAX; ++Z) {
456 size_t nn = (fLevelManagers[Z]).size(); 459 size_t nn = (fLevelManagers[Z]).size();
457 for(size_t j=0; j<nn; ++j) { 460 for(size_t j=0; j<nn; ++j) {
>> 461 //G4cout << " G4NuclearLevelData delete Z= " << Z
>> 462 // << " A= " << AMIN[Z]+j << G4endl;
458 delete (fLevelManagers[Z])[j]; 463 delete (fLevelManagers[Z])[j];
459 } 464 }
460 } 465 }
461 } 466 }
462 467
463 const G4LevelManager* 468 const G4LevelManager*
464 G4NuclearLevelData::GetLevelManager(G4int Z, G << 469 G4NuclearLevelData::GetLevelManager(G4int Z, G4int A, G4bool isLocked)
465 { 470 {
466 if(Z < 1 || Z >= ZMAX || A < AMIN[Z] || A > << 471 const G4LevelManager* man = nullptr;
467 const G4int idx = A - AMIN[Z]; << 472 if(0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) {
468 if( !(fLevelManagerFlags[Z])[idx] ) { << 473 const G4int idx = A - AMIN[Z];
469 G4AutoLock l(&nuclearLevelDataMutex); << 474 if(!(fLevelManagerFlags[Z])[idx]) {
470 if( !(fLevelManagerFlags[Z])[idx] ) { << 475 if(isLocked) {
471 (fLevelManagers[Z])[idx] = fLevelReader- << 476 (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
472 (fLevelManagerFlags[Z])[idx] = true; << 477 (fLevelManagerFlags[Z])[idx] = true;
>> 478 } else {
>> 479 InitialiseForIsotope(Z, A);
>> 480 }
473 } 481 }
474 l.unlock(); << 482 man = (fLevelManagers[Z])[idx];
475 } 483 }
476 return (fLevelManagers[Z])[idx]; << 484 return man;
477 } 485 }
478 486
479 G4bool << 487 G4bool
480 G4NuclearLevelData::AddPrivateData(G4int Z, G4 << 488 G4NuclearLevelData::AddPrivateData(G4int Z, G4int A, const G4String& filename)
481 { 489 {
482 G4bool res = false; 490 G4bool res = false;
>> 491 #ifdef G4MULTITHREADED
>> 492 G4MUTEXLOCK(&nuclearLevelDataMutex);
>> 493 #endif
483 if(Z > 0 && Z < ZMAX && A >= AMIN[Z] && A <= 494 if(Z > 0 && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) {
484 G4AutoLock l(&nuclearLevelDataMutex); << 495 const G4LevelManager* newman =
485 const G4LevelManager* newman = fLevelReade << 496 fLevelReader->MakeLevelManager(Z, A, filename);
486 // if file is corrupted G4LevelReader shou << 497 if(newman) {
487 if(newman != nullptr) { <<
488 res = true; 498 res = true;
489 if(0 < fDeexPrecoParameters->GetVerbose( << 499 G4cout << "G4NuclearLevelData::AddPrivateData for Z= " << Z
490 G4cout << "G4NuclearLevelData::AddPrivateDat << 500 << " A= " << A << " from <" << filename
491 << " A= " << A << " from <" << fname << 501 << "> is done" << G4endl;
492 << "> is done" << G4endl; <<
493 } <<
494 const G4int idx = A - AMIN[Z]; 502 const G4int idx = A - AMIN[Z];
495 delete (fLevelManagers[Z])[idx]; 503 delete (fLevelManagers[Z])[idx];
496 (fLevelManagers[Z])[idx] = newman; 504 (fLevelManagers[Z])[idx] = newman;
497 (fLevelManagerFlags[Z])[idx] = true; 505 (fLevelManagerFlags[Z])[idx] = true;
498 } 506 }
499 l.unlock(); <<
500 } else { 507 } else {
501 G4ExceptionDescription ed; 508 G4ExceptionDescription ed;
502 ed << "private nuclear level data for Z= " 509 ed << "private nuclear level data for Z= " << Z << " A= " << A
503 << " outside allowed limits "; 510 << " outside allowed limits ";
504 G4Exception("G4NuclearLevelData::AddPrivat 511 G4Exception("G4NuclearLevelData::AddPrivateData","had0433",FatalException,
505 ed,"Stop execution"); 512 ed,"Stop execution");
506 } 513 }
>> 514 #ifdef G4MULTITHREADED
>> 515 G4MUTEXUNLOCK(&nuclearLevelDataMutex);
>> 516 #endif
507 return res; 517 return res;
508 } 518 }
509 519
510 G4int G4NuclearLevelData::GetMinA(G4int Z) con 520 G4int G4NuclearLevelData::GetMinA(G4int Z) const
511 { 521 {
512 return (Z >= 0 && Z < ZMAX) ? AMIN[Z] : 0; 522 return (Z >= 0 && Z < ZMAX) ? AMIN[Z] : 0;
513 } 523 }
514 524
515 G4int G4NuclearLevelData::GetMaxA(G4int Z) con 525 G4int G4NuclearLevelData::GetMaxA(G4int Z) const
516 { 526 {
517 return (Z >= 0 && Z < ZMAX) ? AMAX[Z] : 0; 527 return (Z >= 0 && Z < ZMAX) ? AMAX[Z] : 0;
518 } 528 }
519 529
520 void G4NuclearLevelData::UploadNuclearLevelDat << 530 void G4NuclearLevelData::InitialiseForIsotope(G4int Z, G4int A)
521 { 531 {
522 if(fInitialized) return; << 532 if(Z < 1 || Z >= ZMAX || A < AMIN[Z] || A > AMAX[Z]) { return; }
523 G4AutoLock l(&nuclearLevelDataMutex); << 533 const G4int idx = A - AMIN[Z];
524 if(!fInitialized) { << 534 #ifdef G4MULTITHREADED
525 fInitialized = true; << 535 if(!(fLevelManagerFlags[Z])[idx]) {
526 G4int mZ = Zlim + 1; << 536 G4MUTEXLOCK(&nuclearLevelDataMutex);
527 if(mZ > ZMAX) { mZ = ZMAX; } << 537 #endif
528 for(G4int Z=1; Z<mZ; ++Z) { << 538 // initialise only once
529 for(G4int A=AMIN[Z]; A<=AMAX[Z]; ++A) { << 539 // before 1st event fragments Z < zmax are initialized
530 G4int idx = A - AMIN[Z]; << 540 if(!fInitialized) {
531 if( !(fLevelManagerFlags[Z])[idx] ) { << 541 fInitialized = true;
532 (fLevelManagers[Z])[idx] = fLevelReader->C << 542 InitialiseUp(fDeexPrecoParameters->GetUploadZ());
533 (fLevelManagerFlags[Z])[idx] = true; << 543 }
534 } << 544 if(!(fLevelManagerFlags[Z])[idx]) {
>> 545 (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
>> 546 (fLevelManagerFlags[Z])[idx] = true;
>> 547 }
>> 548 #ifdef G4MULTITHREADED
>> 549 G4MUTEXUNLOCK(&nuclearLevelDataMutex);
>> 550 }
>> 551 #endif
>> 552 }
>> 553
>> 554 void G4NuclearLevelData::UploadNuclearLevelData(G4int Z)
>> 555 {
>> 556 #ifdef G4MULTITHREADED
>> 557 G4MUTEXLOCK(&nuclearLevelDataMutex);
>> 558 #endif
>> 559 fDeexPrecoParameters->SetUploadZ(Z);
>> 560 InitialiseUp(Z);
>> 561 #ifdef G4MULTITHREADED
>> 562 G4MUTEXUNLOCK(&nuclearLevelDataMutex);
>> 563 #endif
>> 564 }
>> 565
>> 566 void G4NuclearLevelData::InitialiseUp(G4int ZZ)
>> 567 {
>> 568 G4int mZ = ZZ;
>> 569 if(mZ >= ZMAX) { mZ = ZMAX; }
>> 570 for(G4int Z=1; Z<mZ; ++Z) {
>> 571 for(G4int A=AMIN[Z]; A<=AMAX[Z]; ++A) {
>> 572 G4int idx = A - AMIN[Z];
>> 573 if(!(fLevelManagerFlags[Z])[idx]) {
>> 574 (fLevelManagers[Z])[idx] = fLevelReader->CreateLevelManager(Z, A);
>> 575 (fLevelManagerFlags[Z])[idx] = true;
535 } 576 }
536 } 577 }
537 } 578 }
538 l.unlock(); <<
539 } 579 }
540 580
541 G4double G4NuclearLevelData::GetMaxLevelEnergy 581 G4double G4NuclearLevelData::GetMaxLevelEnergy(G4int Z, G4int A) const
542 { 582 {
543 return (0 < Z && Z < ZMAX && A >= AMIN[Z] && 583 return (0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) ?
544 (G4double)(LEVELMAX[LEVELIDX[Z] + A - AMIN 584 (G4double)(LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0.0;
545 } 585 }
546 586
547 G4float G4NuclearLevelData::MaxLevelEnergy(G4i 587 G4float G4NuclearLevelData::MaxLevelEnergy(G4int Z, G4int A) const
548 { 588 {
549 return (0 < Z && Z < ZMAX && A >= AMIN[Z] && 589 return (0 < Z && Z < ZMAX && A >= AMIN[Z] && A <= AMAX[Z]) ?
550 (LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0. 590 (LEVELMAX[LEVELIDX[Z] + A - AMIN[Z]]) : 0.0f;
551 } 591 }
552 592
553 G4double G4NuclearLevelData::GetLevelEnergy(G4 593 G4double G4NuclearLevelData::GetLevelEnergy(G4int Z, G4int A, G4double energy)
554 { 594 {
555 G4double e = energy; 595 G4double e = energy;
556 if(e <= GetMaxLevelEnergy(Z, A)) { 596 if(e <= GetMaxLevelEnergy(Z, A)) {
557 const G4LevelManager* man = GetLevelManage 597 const G4LevelManager* man = GetLevelManager(Z, A);
558 if(man != nullptr) { 598 if(man != nullptr) {
559 e = (G4double)man->NearestLevelEnergy(e, 599 e = (G4double)man->NearestLevelEnergy(e, man->NumberOfTransitions());
560 } 600 }
561 } 601 }
562 return e; 602 return e;
563 } 603 }
564 604
565 G4double 605 G4double
566 G4NuclearLevelData::GetLowEdgeLevelEnergy(G4in 606 G4NuclearLevelData::GetLowEdgeLevelEnergy(G4int Z, G4int A, G4double energy)
567 { 607 {
568 G4double e = GetMaxLevelEnergy(Z, A); 608 G4double e = GetMaxLevelEnergy(Z, A);
569 if(energy < e) { 609 if(energy < e) {
570 const G4LevelManager* man = GetLevelManage 610 const G4LevelManager* man = GetLevelManager(Z, A);
571 if(man != nullptr) { 611 if(man != nullptr) {
572 e = (G4double)man->NearestLowEdgeLevelEn 612 e = (G4double)man->NearestLowEdgeLevelEnergy(energy);
573 } 613 }
574 } 614 }
575 return e; 615 return e;
576 } 616 }
577 617
578 G4double 618 G4double
579 G4NuclearLevelData::FindLevel(G4int Z, G4int A 619 G4NuclearLevelData::FindLevel(G4int Z, G4int A, G4double resMass,
580 G4double Mass, G 620 G4double Mass, G4double partMass, G4double ekin)
581 { 621 {
582 G4double T = ekin; 622 G4double T = ekin;
583 G4double E2 = (Mass - partMass)*(Mass - part 623 G4double E2 = (Mass - partMass)*(Mass - partMass);
584 G4double Eex = std::sqrt(E2 - 2.0*Mass*T) - 624 G4double Eex = std::sqrt(E2 - 2.0*Mass*T) - resMass;
585 if(Eex <= GetMaxLevelEnergy(Z, A)) { 625 if(Eex <= GetMaxLevelEnergy(Z, A)) {
586 if(Eex <= 0.0) { 626 if(Eex <= 0.0) {
587 Eex = 0.0; 627 Eex = 0.0;
588 T = (E2 - resMass*resMass)*0.5/Mass; 628 T = (E2 - resMass*resMass)*0.5/Mass;
589 } else { 629 } else {
590 const G4LevelManager* man = GetLevelMana 630 const G4LevelManager* man = GetLevelManager(Z, A);
591 if(man != nullptr) { 631 if(man != nullptr) {
592 size_t idx = man->NearestLevelIndex(Ee 632 size_t idx = man->NearestLevelIndex(Eex, man->NumberOfTransitions());
593 for(;;) { 633 for(;;) {
594 Eex = (G4double)man->LevelEnergy(idx 634 Eex = (G4double)man->LevelEnergy(idx);
595 T = (E2 - (resMass+Eex)*(resMass+Eex 635 T = (E2 - (resMass+Eex)*(resMass+Eex))*0.5/Mass;
596 if(T >= 0.0 || 0 == idx) { break; } 636 if(T >= 0.0 || 0 == idx) { break; }
597 --idx; 637 --idx;
598 } 638 }
599 } 639 }
600 } 640 }
601 T = std::max(T, 0.0); 641 T = std::max(T, 0.0);
602 } 642 }
603 return T; 643 return T;
604 } 644 }
605 645
606 G4DeexPrecoParameters* G4NuclearLevelData::Get 646 G4DeexPrecoParameters* G4NuclearLevelData::GetParameters()
607 { 647 {
608 return fDeexPrecoParameters; 648 return fDeexPrecoParameters;
609 } 649 }
610 650
611 G4PairingCorrection* G4NuclearLevelData::GetPa 651 G4PairingCorrection* G4NuclearLevelData::GetPairingCorrection()
612 { 652 {
613 return fPairingCorrection; 653 return fPairingCorrection;
614 } 654 }
615 655
616 G4ShellCorrection* G4NuclearLevelData::GetShel 656 G4ShellCorrection* G4NuclearLevelData::GetShellCorrection()
617 { 657 {
618 return fShellCorrection; 658 return fShellCorrection;
619 } 659 }
620 660
621 G4double G4NuclearLevelData::GetLevelDensity(G 661 G4double G4NuclearLevelData::GetLevelDensity(G4int Z, G4int A, G4double U)
622 { 662 {
623 if (fDeexPrecoParameters->GetLevelDensityFla << 663 if(fDeexPrecoParameters->GetLevelDensityFlag()) {
624 return A*fDeexPrecoParameters->GetLevelDen 664 return A*fDeexPrecoParameters->GetLevelDensity();
625 } 665 }
626 const G4LevelManager* man = GetLevelManager( 666 const G4LevelManager* man = GetLevelManager(Z, A);
627 return (nullptr != man) ? man->LevelDensity( << 667 return (man) ? man->LevelDensity(U)
628 : 0.058025*A*(1.0 + 5.9059/fG4calc->Z13(A) 668 : 0.058025*A*(1.0 + 5.9059/fG4calc->Z13(A));
629 } 669 }
630 670
631 G4double G4NuclearLevelData::GetPairingCorrect 671 G4double G4NuclearLevelData::GetPairingCorrection(G4int Z, G4int A)
632 { 672 {
633 return fPairingCorrection->GetPairingCorrect << 673 if(fDeexPrecoParameters->GetLevelDensityFlag()) {
>> 674 return fPairingCorrection->GetPairingCorrection(A, Z);
>> 675 }
>> 676 G4int N = A - Z;
>> 677 const G4double par = 12.*CLHEP::MeV;
>> 678 G4double x = (A <= 36) ? 6.0 : std::sqrt((G4double)A);
>> 679 return (2 - Z + (Z/2)*2 - N + (N/2)*2)*par/x;
634 } 680 }
635 681
636 void G4NuclearLevelData::StreamLevels(std::ost 682 void G4NuclearLevelData::StreamLevels(std::ostream& os,
637 G4int Z, 683 G4int Z, G4int A)
638 { 684 {
639 const G4LevelManager* man = GetLevelManager( 685 const G4LevelManager* man = GetLevelManager(Z, A);
640 if (man) { << 686 if(man) {
641 os << "Level data for Z= " << Z << " A= " 687 os << "Level data for Z= " << Z << " A= " << A << " "
642 << man->NumberOfTransitions() + 1 << " 688 << man->NumberOfTransitions() + 1 << " levels \n";
643 man->StreamInfo(os); 689 man->StreamInfo(os);
644 } 690 }
645 } 691 }
646 692