isonav/TODO at master · ffavela/isonav · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
Solve this BUG!!:
$ isonav 197Au --Elab=64 --material=Si --thick=30
-1.000

But

$ isonav 197Au --Elab=0.001 --material=Si --thick=300
5884582.250

This makes no sense!!! BUG solved

Improve the documentation. Make a latex doc for the theory.

Make the interactive part of the program, less prone to errors.

Parse que extra data in the db (parity etc).

Add the relativistic scattering case.

Decays with excited states.

Improve the database, add half lifes, decay modes. Beta decay.
Add a get reference option.
Add a wikipedia parser for HF?
Add subnuclear particles?

Spectroscopic functions. Plotting stuff.

Testing cases for debugging.

# Use the database for masses and other properties from (suggested by
#    Silpion):

# https://www-nds.iaea.org/amdc/ame2012/mass.mas12
DONE!! (only for the masses)


Make a man page: example;  $ man isonav

Add coulomb energy for binary exit channels. DONE

Include an Elab option for looking up to which channels get opened.

Add additionally a sorting with respect of coulomb for the opened channels.

Add capability of printing higher order reactions (ternary, quaternary ...)
BUGS:

    isotopes with names 1950s and 1940s appear sometimes, for example:

        $ python isonav.py 235U --decay -v | head
         #Given an isotope, it gives all the possible decay schemes
         #res    daught          eRes    eDaugh  Q
          40S     195Os           77.659  15.923  93.582
          41S     194Os           76.083  16.075  92.157
          37P     198Ir           72.564  13.554  86.118
          34Si    201Pt           73.201  12.376  85.576
          38P     197Ir           70.896  13.672  84.568

BUG

$ isonav a -l
Error

$ isonav 4He -l
1	0.0
2	20.21
3	21.01
4	21.84
5	23.33
6	23.64
7	24.25
8	25.28
9	25.95
10	27.42
11	28.31
12	28.37
13	28.39
14	28.64
15	28.67
16	29.89

No error... this also happens with p,t,d (not that there are but
still). Note the listed levels that are not 0.0 are theoretical, put
some flag for this.