The inputs to SFresco are
Card 1: name of search_file
Lines 2: repeatedly one of (in all upper or all lower case):
EX: exit
Q: query status of search variables
V: query status of unfixed search variables
CHA: query status of changed search variables
CHI: list status by dataset
SET var val: set variable number var to value val.
FIX var: fix variable number var (set step=0).
FXWID : fix all R-matrix widths variables.
FXRPE : fix all R-matrix pole energy variables.
FXPBG : fix all pole energies variables with 'BG' in their name.
FXBG : fix all variables with 'BG' in their name.
FXNOR : fix all dataset norm variables (kind=5)
FXSHF : fix all dataset energy shift variables (kind=6)
FREE : free all R-matrix energies and widths variables.
STEP var step: unfix variable number var with step step.
SCAN var val1 val2 step: scan variable number var
from value val1 to value val2 in steps of step.
SHOW: list all datasets with current predictions and values.
PLOT plotfile: write file (default: search.plot) for reading by xmgr, xvgr,
with data and theoretical curves.
DONE search_file: write single new search file (default: done.fresco) for rereading by SFRESCO as a combination of FRESCO input_file and current parameter values. Data specifications are however not included.
LINE plotfile: write file (default: search.plot) for reading by xmgr, xvgr,
with theoretical curves only.
READ file: read plot output file for further searches, if not:
READ snapfile: if the input string file contains the string 'snap',
read to the end of previous snap output snapfile, for further searches.
ELIM emin emax: restrict data scans to this energy range for plots and sums
POLES emin emax: vary only R-matrix poles in this energy range (others are fixed)
ESCAN emin emax estep: scan lab. energy in incident channel, to give
phase shifts in fort.71, fusion & reaction cross sections in fort.40,
and S-factors in fort.35 and 75.
This is to look for/at resonances, or at behaviour out of data's energy range,
not for fitting.
MIN: call minuit interactively, to read from stdin.
Return to SFresco by minuit command end.
A example of this search_file input is:
'6hep.in' '6hep.out' 7 1 &variable kind=1 name='V' kp=1 pline=2 col=1 / &variable kind=1 name='r0' kp=1 pline=2 col=2 valmin=.7 valmax=2.0/ &variable kind=1 name='a' kp=1 pline=2 col=3 / &variable kind=1 name='W' kp=1 pline=2 col=4 / &variable kind=1 name='rw' kp=1 pline=2 col=5 valmin=.7 valmax=2.0/ &variable kind=1 name='aw' kp=1 pline=2 col=6 / &variable kind=5 name='exptnorm' / &data type=0 iscale=2 idir=0 lab=F abserr=T/ 13.99 565.7 28.3 17.59 478.19 23.91 21.22 425.41 21.28 24.92 338.44 16.93 28.68 284.71 14.24 36.48 168.59 8.43 40.56 137.9 6.86 44.8 104.75 5.24 53. 30.0 9.14 53.94 38.45 2.71 59. 25.62 2.55 &A example of the full SFresco input, suitable for batch running, is:
search4.in min fix 4 migrad end q show plotwhere search4.in is the above search_file input file name.