Skip to content

Commit

Permalink
add legacy folder for og implementation
Browse files Browse the repository at this point in the history
  • Loading branch information
KedoKudo committed Jan 1, 2025
1 parent cfb252a commit 1dd4264
Show file tree
Hide file tree
Showing 34 changed files with 6,700 additions and 5 deletions.
3 changes: 0 additions & 3 deletions .gitignore
Original file line number Diff line number Diff line change
Expand Up @@ -201,6 +201,3 @@ cython_debug/

# Ignore vscode stuff.
.vscode/

# Ignore legacy stuff.
legacy/
4 changes: 2 additions & 2 deletions .pre-commit-config.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -17,9 +17,9 @@ repos:
hooks:
- id: ruff
args: [--fix, --exit-non-zero-on-fix]
exclude: "tests/data/.*"
exclude: "(tests/data/.*|legacy/.*)"
- id: ruff-format
exclude: "tests/data/.*"
exclude: "(tests/data/.*|legacy/.*)"
- repo: https://github.com/codespell-project/codespell
rev: v2.3.0
hooks:
Expand Down
56 changes: 56 additions & 0 deletions legacy/Notes/Pleiades.notes
Original file line number Diff line number Diff line change
@@ -0,0 +1,56 @@
General structure of git repo

[PLEIADES]
docs/
- Makefile
- _build/
- api.rst
- conf.py
- examples.rst
- generated
- index.rst
- make.bat
- requirements.in
- requirements.txt
- source/
- modules.rst
- pleiades.rst
- usage.rst
notes/
- Pleiades.notes [this file]
- SammyInputFile.notes
- SammyOutputFile.notes
- SammyParameterFile.notes
examples/
- configFiles/
- nat_europium.ini
- uranium.ini
- data/
- natEu.twenty
- u235-u238.twenty
- Notebooks/
- natEu_fit_tutorial.ipynb
- u_fit_tutorial.ipynb
- scripts/
- makeCompoundFit.py
nucDataLibs/
- isotopeInfo/
- isotopes.info
- mass.mas20
- neutrons.list
- resonanceTables/
- res_endf8.endf
- xSections/
- [cross-section files for select isotopes]
pleiades/
- __pycache__
- __init__.py
- sammyInput.py
- sammyPlotter.py
- nucData.py
- simData.py
tests/
- test_nucData.py
- test_sammyInput.py
- test_sammyRunner.py
- test_sammyUtils.py
123 changes: 123 additions & 0 deletions legacy/Notes/SammyInputFile.notes
Original file line number Diff line number Diff line change
@@ -0,0 +1,123 @@
Max Number of Columns in the fortran input file = 80

Data input types
“A” can be numbers or letters or symbols
“F” requires a number with a decimal point
“I” implies a number without a decimal point like an integer and requires that the number be in the right-most columns

Classes:
Card 1 -
Line 1:
1-80 columns, A "TITLE"
Card 2 -
Line 1:
01-10 columns, A "ELMNT"; # Sample element’s name
11-20 columns, F "AW"; # Atomic weight (amu)
21-30 columns, F "EMIN"; # Minimum energy for this data set (eV)
31-40 columns, F "EMAX"; # Maximum energy for this data set (eV)
41-45 columns, I "NEPNTS"; # Number of points to be used in generating artificial energy grid (default = 10001)
46-50 columns, I "ITMAX"; # Number of iterations (default = 2)
51-52 columns, I "ICORR"; # Correlations smaller than this value (divided by 100) are not to be printed
53-55 columns, I "NXTRA"; # Number of extra points to be added between each pair of data points for auxiliary energy grid
56-57 columns, I "IPTDOP"; # Number of points to be added to auxiliary energy grid across small resonances
59-60 columns, I "IPTWID"; # Determines the number of points to be added to auxiliary grid in tails of small resonances
61-70 columns, I "IXXCHN"; # Number of energy channels in ODF-type data file to be ignored
71-72 columns, I "NDIGIT"; # Number of digits forcompact format for covariance matrix
73-74 columns, I "IDROPP"; # The input resonance parameter covariance matrix will be modified before being used in the fitting procedure
75-80 columns, I "MATNUM" # ENDF Material Number
Card 3 -
Note: Alphanumeric information concerning program options.
Each command is its own line. This is followed by a blank line.
Lines: various command statements.

Card 4 -
Note: This card set is present only when the phrase “KEY WORD PARTICLE PAir
definitions are given” occurs in card set 3 above
Line 1:
[FILL IN LATER]

Card 5 -
Line 1:
01-10 columns, F "TEMP"; # Effective temperature of the sample (K);
11-20 columns, F "DIST"; # Flight-path length (m)
21-30 columns, F "DELTAL"; # That portion of the Gaussian resolution width attributed to the flight-path length (m).
31-40 columns, F "DELTAE"; # E-folding width of exponential resolution function (μs);
41-50 columns, F "DELTAG"; # The absolute value of DELTAG is the FWHM (μs) of a Gaussian resolution function representing the burst width
51-60 columns, F "DELTTT"; # Approximate step size for the Effective Temperature for the Leal-Hwang method of Doppler broadening (K).
61-70 columns, F "ELOWBR"; # The energy below which no broadening is wanted, for highenergy Gaussian approximation to Doppler broadening (eV)
71-80 columns, F "DDDEEE"; # Range for straight-line energy average resolution function (eV)

Card 6 -
Note: This card set is included only if DELTAG < 0 in card set 5.
Line 1:
01-10 columns, F "DELTAB" # Minimum channel width for interval from EMIN to EMAX (μs)
11-15 columns, I "NCF" # Number of crunch boundaries (energy intervals)
Line 2:
01-10 columns, F "BCF(I)" # maximum energy for this crunch factor (read in order of increasing energy) (eV)
11-20 columns, F "CF(I)" # crunch factors (integer units of DELTAB)

Card 7 -
Line 1:
01-10 columns, F "CRFN" # Channel radius (fermi)
11-20 columns, F "THICK" # Sample thickness (atoms/barn)
21-30 columns, F "DCOVA" # Constant term in data covariance (default = 0.0)
31-40 columns, F "DCOVB" # Linear term in data covariance (default = 0.0)
41-50 columns, F "VMIN" # Minimum absolute uncertainty (standard deviation) on experimental data

Card 8 -
Note: Line 1 will contain a "command" string that represents the type
of data, which must begin in first column on the line. Only those
characters in CAPITAL letters are required. Lines 2 and greater
are not needed unless "DIFFErential elastic cross sections" or
"DIFFERENTIAL REaction" is selected.
Line 1:
1-80 columns, A "CROSS"
Line 2-N:

Card 9 -
Note: No longer used.

Card 10.1 -
Note: It is no longer necessary to have “USE NEW SPIN GROUP Format” in
card set 3 in order to use this alternative; card set 10.1 is now
the default.
Line 2 is repeated for a total of (NENT + NEXT) lines!
Line 1 followed by multiple line 2s is repeated for all spin groups
Line 1:
1-3 column, I "JJ" # Spin group number: 1,2,3, ...
5 column, A "EXCL" # Flag to exclude this spin group from calculation of cross section
6-10 column, I "NENT(JJ)" # Number of entrance channels (must be > 0)
11-15 column, I "NEXT(JJ)" # Number of exit channels, excluding those that are also entrance channels (can be zero)
16-20 column, F "SPINJ(JJ)" # Integer or half-integer spin for resonances in group JJ
21-30 column, F "ABNDNC(JJ)" # Isotopic abundance for this spin group (These values may be overwritten in the PARameter file)
31-35 column, F "SPINI" # Ground-state spin for the nuclide (i.e., for the target particle)
35-80 column, A "Comments" # Use this to state which nuclei is covered by this spin group!

Line 2:
4-5, I "N" # Channel number (One line is needed for each channel!)
6-8, I "KZ1(N,JJ)" # Charge number Z1 for the first particle in this channel
10, I "LPENT(N,JJ)" # Is SAMMY to calculate the penetrability for channel N, spin group JJ?
11-13, I "KZ2(N,JJ)" # Charge number Z2 for the second particle in this channel
15, I "ISHIFT(N,JJ)" # 1 if calculate shift; 0 Otherwise
18, I "IFEXCL" # 1 if exclude this channel from final-state calculation for reaction cross section; 0 if include
19-20, I "LSPIN" # Orbital angular momentum for channel
21-30, F "CHSPIN(N,JJ)" # Channel spin (integer or halfinteger, sign implies parity)
31-40, F "ENBND(N,JJ)" # Boundary condition (eV). Used only if ISHIFT = 1
41-50, F "ECHAN(N,JJ)" # Excitation energy of residual nucleus in an inelastic channel (eV)
51-55, F "RDEFF(N,JJ)" # Effective radius (fermi)
56-60, F "RDTRU(N,JJ)" # True radius (fermi) (Value may be overwritten in the PARameter file.)
61-70, F "EMM1(N,JJ)" # Mass of first particle in this channel (that is, of the target nuclide, in the case of an incident channel) (amu)
71-80, F "EMM2(N,JJ)" # Mass of second particle in the channel (that is, of the neutron or other incident particle, for entrance channels) (amu)



Run time notes:
Need to use the lowest mass in a compound nuclei fit otherwise you get this:
***************************************************************************
* Warning -- Mass used to determine limits for auxiliary grid (for Dopplr *
* and multiple-scattering corrections) may be too big. SAMMY is using the *
* mass from nuclide number 2 (238.050770), because that is closest to the *
* mass from Line 2 of the INPut file (238.050800). For more accurate *
* results, use the mass of the smallest nuclide (235.043940). (Expect *
* longer runtime if you change to the smaller mass.) *
***************************************************************************
14 changes: 14 additions & 0 deletions legacy/Notes/SammyOutputLstFile.notes
Original file line number Diff line number Diff line change
@@ -0,0 +1,14 @@
LST FILE for plotting
Column 1: Energy (Units are as specified by input command. If not specified, then units are eV if appropriate, otherwise keV.)
Column 2: Experimental cross section (barns)
Column 3: Absolute uncertainty in experimental cross section (barns)
Column 4: Zeroth-order theoretical cross section as evaluated by SAMMY (barns)
Column 5: Final theoretical cross section as evaluated by SAMMY (barns)
Column 6: Experimental transmission (dimensionless)
Column 7: Absolute uncertainty in experimental transmission
Column 8: Zeroth-order theoretical transmission as evaluated by SAMMY
Column 9: Final theoretical transmission as evaluated by SAMMY
Column 10: Theoretical uncertainty on section 4 if data were total cross section, or on section 8 if data were transmission
Column 11: Theoretical uncertainty on section 5 or section 9
Column 12: Adjusted energy initially, when t0 and L0 are to be varied. See Section III.E.8 of this manual.
Column 13: Adjusted energy
Loading

0 comments on commit 1dd4264

Please sign in to comment.