release 1.4 code

History

@@ -1,5 +1,22 @@
 
 
+================================================================================
+                               October 2024
+--------------------------------------------------------------------------------
+E. Mendoza
+  - SaG4n.01.04
+     + We have included the possibility of using the alphas of a given
+     radioactive decay in the source, without having to write all the alpha
+     energies and intensities in the input. For this, a database has been added
+     in this release, in the new directory “data/AlphaSpectra”. The database has
+     been obtained mainly from the ENDF/B-VIII.0 decay data library.
+     + A new keyword is introduced in the input file in order to provide the path
+     to the mentioned database.
+     + A new example with makes use of this database is provided in examples/fuel.
+     + A new output histogram is provided by the code, “Source_00”, with the
+     energy spectrum of the alpha particle source used in the calculation.
+================================================================================
+
 ================================================================================
                                March 2022
 --------------------------------------------------------------------------------

README.md

@@ -33,15 +33,13 @@ The SaG4n code, together with different (alpha,xn) data libraries can be downloa
 
 ## Directories
 
-- `CMakeLists.txt`: needed to compie SaG4n
-- `SaG4n.cc`: SaG4n source code
 - `include`: SaG4n source code
 - `src`: SaG4n source code
 - `PhysicsListFiles`: SaG4n source code (modified G4ParticleHP)
 - `inputs/examples`: Some SaG4n input examples
 - `outputs/examples_reference`: reference outputs of the input examples
+- `data`: SaG4n database
 - `docs`: SaG4n manual
-- `History`: list of changes along the SaG4n releases
 
 
 ## Installation
@@ -112,7 +110,10 @@ Where ENDFLIBING4FORMAT must be one of the (alpha,xn) cross section data librari
 
 ## Examples
 
-Some input examples are distributed together with the code. They are located inside the inputs/examples directory. Outputs generated with SaG4n using these inputs can be found in `outputs/examples_reference/`. They have been generated using Geant4.11.0.0 and the JENDLTENDL01 data library.
+Some input examples are distributed together with the code. They are located inside the `inputs/examples` directory. Outputs generated with SaG4n using these inputs can be found in `outputs/examples_reference/`. They have been generated using Geant4.11.2.1 and the JENDLTENDL01 data library.
+
+We have found some differences in the results of the examples compared to those distributed in previous releases of SaG4n. Specifically, in the results of `yields/yield01.inp` and `yields/yield03.inp`, which corresponds to the neutron production in Li and B.  The differences come from changes in the Geant4, and not in SaG4n. If the inputs are run using Geant4.11.00.0, instead of Geant4.11.2.1, then the same results are obtained as those distributed in previous SaG4n releases.
+
 
 The `outputs/examples` directory does not content any output file. The idea is to place there the output files generated by SaG4n, and compare then with the reference output files.
 
@@ -200,7 +201,7 @@ Again, to run them:
 ```
 
 
-- `source01.inp`, `source02.inp`, `source03.inp` and `source04.inp`: they are the same input but with different sources. The results from `source01.inp` have to be the same as the sum of the results from `source02.inp`, `source03.inp` and `source04.inp`.
+ - `source01.inp`, `source02.inp`, `source03.inp` and `source04.inp`: they are the same input but with different sources. The results from `source01.inp` have to be the same as the sum of the results from `source02.inp`, `source03.inp` and `source04.inp`.
 
 ```sh
 ./build/SaG4n inputs/examples/basic01/source01.inp
@@ -218,6 +219,23 @@ Again, to run them:
       + momentum direction: they are are emittedin the (0,0,1) direction.
       + position: they are emitted from a 2 cm radius flat disk perpendicular to the Z axis and centered in (0,0,-1) cm. The position of the alphas are sampled according to a 2D-gaussian distribution centered in the center of the disk and with 0.7 cm standard deviations in X and Y, and always inside the disk (i.e. 2 cm is the maximum distance from the Z axis).
 
+```sh
+./build/SaG4n inputs/examples/beam/beam01.inp
+```
+
+### fuel
+
+`inputs/examples/fuel/`: here we calculate the neutron emission in PuO2. The geometry is a shpere made of PuO2, with a certain Pu isotopic composition. The alpha particles come from the decay of the different Pu isotopoes. We use the database introduced in SaG4n.01.04 to obtain the alpha spectra. The results are "per Bq", i.e. normalized to the number of alpha particles. There are two different inputs, which should lead to the same result:
+   - fuel01.inp: the intensities of the alpha spectra are provided for each Pu isotope.
+   - fuel02.inp: the intensities of the alpha spectra are not provided directly, but the half-lifes and isotopic abundances. 
+
+
+```sh
+./build/SaG4n inputs/examples/fuel/fuel01.inp
+./build/SaG4n inputs/examples/fuel/fuel02.inp
+```
+
+
 ### decay series
 
 Additionally, sources corresponding to the Th-232, U-235 and U-238 decay series in secular equilibrium are located in `inputs/examples/sources/`

data/AlphaSpectra/Ac207

@@ -0,0 +1,2 @@
+1
+89207  7.6935             1

data/AlphaSpectra/Ac208

@@ -0,0 +1,2 @@
+1
+89208  7.5725          0.99

data/AlphaSpectra/Ac209

@@ -0,0 +1,2 @@
+1
+89209  7.5775             1

data/AlphaSpectra/Ac211

@@ -0,0 +1,2 @@
+1
+89211  7.4775             1

data/AlphaSpectra/Ac213

@@ -0,0 +1,2 @@
+1
+89213  7.3645             1

data/AlphaSpectra/Ac215

@@ -0,0 +1,5 @@
+4
+89215  6.9625    0.00139874
+89215  7.0285     0.0019982
+89215  7.2165    0.00459586
+89215  7.6045      0.991107

data/AlphaSpectra/Ac220

@@ -0,0 +1,16 @@
+15
+89220  7.6225          0.04
+89220  7.6355          0.04
+89220  7.6525          0.09
+89220  7.6645          0.04
+89220  7.6705          0.08
+89220  7.7095          0.11
+89220  7.7925           0.1
+89220  7.8505          0.05
+89220  7.8555          0.26
+89220  7.9445          0.02
+89220  7.9715          0.04
+89220  8.0065          0.03
+89220  8.0555          0.04
+89220  8.0635          0.02
+89220  8.1945          0.04

data/AlphaSpectra/Ac221

@@ -0,0 +1,5 @@
+4
+89221  7.1705         0.024
+89221  7.3775          0.09
+89221  7.4405          0.21
+89221  7.6455          0.68

data/AlphaSpectra/Ac222

@@ -0,0 +1,3 @@
+2
+89222  6.9635        0.0594
+89222  7.0085        0.9306

data/AlphaSpectra/Ac223

@@ -0,0 +1,30 @@
+29
+89223  5.8925       9.9e-05
+89223  5.9665      0.000297
+89223  6.0235       9.9e-05
+89223  6.0825      0.000297
+89223  6.1345      0.001188
+89223  6.1405      0.000297
+89223  6.1635      0.000495
+89223  6.1775      0.009405
+89223  6.2055      0.000297
+89223  6.2235      5.94e-05
+89223  6.2355      0.000891
+89223  6.2815      0.000495
+89223  6.2925      0.004653
+89223  6.2935      0.000594
+89223  6.3255       0.00297
+89223  6.3325      0.001386
+89223  6.3415      0.000495
+89223  6.3605      0.002178
+89223  6.3965      0.001287
+89223  6.4485       0.00198
+89223  6.4545      0.000594
+89223  6.4735       0.03069
+89223  6.5235       0.00594
+89223  6.5285       0.03069
+89223  6.5635       0.13563
+89223  6.5815       0.00297
+89223  6.6065       0.00495
+89223  6.6465       0.44154
+89223  6.6615       0.31482

data/AlphaSpectra/Ac224

@@ -0,0 +1,28 @@
+27
+89224  5.6445      5.46e-05
+89224  5.7125     0.0001092
+89224  5.7215     0.0001092
+89224  5.7425      4.55e-05
+89224  5.7715     0.0002184
+89224  5.7795      5.46e-05
+89224  5.8065      1.82e-05
+89224  5.8405     0.0002366
+89224  5.8445     0.0005005
+89224  5.8545      6.37e-05
+89224  5.8575     0.0002275
+89224  5.8645     0.0006825
+89224  5.8715       9.1e-05
+89224  5.8805      0.001547
+89224  5.9055     0.0001274
+89224  5.9125     0.0001365
+89224  5.9195     0.0008736
+89224  5.9455      0.004004
+89224  5.9615      2.73e-05
+89224  5.9715     0.0001638
+89224  6.0035      0.006097
+89224  6.0175      0.001274
+89224  6.0595      0.019929
+89224  6.1415      0.023296
+89224  6.1585     0.0009373
+89224  6.2075      0.010829
+89224  6.2135      0.018564

data/AlphaSpectra/Ac225

@@ -0,0 +1,30 @@
+29
+89225  5.0215         5e-06
+89225  5.0675         3e-05
+89225  5.0915         6e-05
+89225  5.1315         4e-05
+89225  5.2025         2e-05
+89225  5.2115        0.0003
+89225  5.2395         3e-05
+89225  5.2705       0.00014
+89225  5.2865        0.0023
+89225  5.3215        0.0007
+89225  5.4365        0.0007
+89225  5.4435        0.0014
+89225  5.5405       0.00015
+89225  5.5455        0.0003
+89225  5.5545         0.001
+89225  5.5635       0.00034
+89225  5.5805         0.012
+89225  5.5985        0.0004
+89225  5.6095         0.011
+89225  5.6375         0.044
+89225  5.6825         0.013
+89225  5.7245         0.031
+89225  5.7315        0.0087
+89225  5.7325        0.0932
+89225  5.7905         0.086
+89225  5.7925         0.181
+89225  5.8055         0.003
+89225  5.8295       0.00023
+89225  5.8305         0.507

data/AlphaSpectra/Ac226

@@ -0,0 +1,2 @@
+1
+89226  5.3995         6e-05

data/AlphaSpectra/Ac227

@@ -0,0 +1,20 @@
+19
+89227  4.3635      4.14e-07
+89227  4.4235      8.28e-07
+89227  4.4455       6.9e-06
+89227  4.4585       6.9e-07
+89227  4.5115      4.14e-07
+89227  4.5815      4.14e-07
+89227  4.5895      1.38e-06
+89227  4.5945      2.76e-06
+89227  4.7155      5.52e-05
+89227  4.7385     1.242e-05
+89227  4.7685     0.0002484
+89227  4.7855     1.104e-05
+89227  4.7965      0.000138
+89227  4.8225      9.66e-06
+89227  4.8555      0.000828
+89227  4.8725     0.0008694
+89227  4.8985     1.518e-05
+89227  4.9405     0.0054648
+89227  4.9535     0.0065826

data/AlphaSpectra/Am237

@@ -0,0 +1,2 @@
+1
+95237  6.0425       0.00025

data/AlphaSpectra/Am238

@@ -0,0 +1,2 @@
+1
+95238  5.9405         1e-06

data/AlphaSpectra/Am239

@@ -0,0 +1,5 @@
+4
+95239  5.6805      1.98e-06
+95239  5.7345     1.375e-05
+95239  5.7745      8.37e-05
+95239  5.8255       3.3e-07

data/AlphaSpectra/Am240

@@ -0,0 +1,4 @@
+3
+95240  5.2865     2.337e-08
+95240  5.3375      2.28e-07
+95240  5.3785    1.6492e-06

data/AlphaSpectra/Am241

@@ -0,0 +1,26 @@
+25
+95241  4.7585         5e-08
+95241  4.8005       8.6e-07
+95241  4.8335         7e-06
+95241  5.0035         1e-06
+95241  5.0675       1.4e-06
+95241  5.0895         4e-06
+95241  5.0965         4e-06
+95241  5.1145         4e-06
+95241  5.1365       3.2e-06
+95241  5.1555         7e-06
+95241  5.1785         3e-06
+95241  5.1825         9e-06
+95241  5.1925         6e-06
+95241  5.2165         1e-07
+95241  5.2235       1.3e-05
+95241  5.2445       2.4e-05
+95241  5.2795         5e-06
+95241  5.3225       0.00015
+95241  5.3885        0.0166
+95241  5.4165        0.0001
+95241  5.4425         0.131
+95241  5.4695        0.0002
+95241  5.4855         0.848
+95241  5.5115       0.00225
+95241  5.5445        0.0037

data/AlphaSpectra/Am242m

@@ -0,0 +1,16 @@
+15
+95242  4.9745      9.18e-08
+95242  5.0275      9.18e-07
+95242  5.0675    1.1475e-05
+95242  5.0825     1.377e-06
+95242  5.0915     9.639e-06
+95242  5.1435    0.00025704
+95242  5.1535      9.18e-07
+95242  5.2075     0.0040851
+95242  5.2145     1.377e-06
+95242  5.2725     5.049e-05
+95242  5.3145     2.754e-05
+95242  5.3315     6.885e-06
+95242  5.3675     5.049e-05
+95242  5.4105      4.59e-05
+95242  5.4585     6.426e-06

data/AlphaSpectra/Am243

@@ -0,0 +1,16 @@
+15
+95243  4.6955       1.7e-05
+95243  4.9185       8.5e-07
+95243  4.9305       1.8e-06
+95243  4.9465       3.4e-06
+95243  4.9975         8e-06
+95243  5.0085         8e-06
+95243  5.0295       1.5e-05
+95243  5.0355       1.5e-05
+95243  5.0885         5e-05
+95243  5.1135         8e-05
+95243  5.1815        0.0136
+95243  5.2335         0.112
+95243  5.2755         0.871
+95243  5.3215        0.0016
+95243  5.3495        0.0016

data/AlphaSpectra/At204

@@ -0,0 +1,2 @@
+1
+85204  5.9505        0.0391

data/AlphaSpectra/At205

@@ -0,0 +1,2 @@
+1
+85205  5.9025           0.1

data/AlphaSpectra/At206

@@ -0,0 +1,5 @@
+4
+85206  5.7025      0.008613
+85206  5.7345       9.9e-05
+85206  5.7675      0.000207
+85206  5.7745       8.1e-05

data/AlphaSpectra/At207

@@ -0,0 +1,2 @@
+1
+85207  5.7585         0.086

data/AlphaSpectra/At208

@@ -0,0 +1,4 @@
+3
+85208  5.5865     4.785e-05
+85208  5.6265     0.0001177
+85208  5.6405     0.0053295

data/AlphaSpectra/At209

@@ -0,0 +1,2 @@
+1
+85209  5.6475      0.040959

data/AlphaSpectra/At210

@@ -0,0 +1,9 @@
+8
+85210  5.1755     3.675e-06
+85210  5.2425     1.575e-05
+85210  5.3615     0.0004865
+85210  5.3865      8.05e-05
+85210  5.4425      0.000497
+85210  5.4565         7e-06
+85210  5.4655      0.000126
+85210  5.5245    0.00053375

data/AlphaSpectra/At211

@@ -0,0 +1,5 @@
+4
+85211  4.9975      4.18e-06
+85211  5.1415     1.045e-05
+85211  5.2105    3.8456e-05
+85211  5.8695         0.418

data/AlphaSpectra/At212

@@ -0,0 +1,9 @@
+8
+85212  6.6215       0.00135
+85212  6.7585        0.0007
+85212  6.7985       0.00048
+85212  7.0495         0.004
+85212  7.0795        0.0059
+85212  7.1685        0.0015
+85212  7.6075         0.154
+85212  7.6695         0.832

data/AlphaSpectra/At212m

@@ -0,0 +1,10 @@
+9
+85212  6.8125        0.0036
+85212  6.9475       0.00052
+85212  7.0225        0.0013
+85212  7.2515        0.0038
+85212  7.2725        0.0036
+85212  7.2985        0.0007
+85212  7.3885        0.0035
+85212  7.8265         0.676
+85212  7.8885         0.307

data/AlphaSpectra/At217

@@ -0,0 +1,5 @@
+4
+85217  6.3225   4.99965e-05
+85217  6.4845   0.000209985
+85217  6.8135   0.000359975
+85217  7.0665       0.99893

data/AlphaSpectra/At218

@@ -0,0 +1,4 @@
+3
+85218  6.6535      0.063936
+85218  6.6935        0.8991
+85218  6.7565      0.035964

data/AlphaSpectra/At219

@@ -0,0 +1,2 @@
+1
+85219  6.2085          0.97