This directory concerns the narrow resonances, the resonance ranges,
the BW shapes physical and bare (undressed fro vacuum polarization).
TH parameters are in the main directory resonances_dat.f procedures:

resonances0n.f writes BW resonance shapes in plot files {\tt
rbwxxx.dap} where xxx=name of resonances omega, phi, psi1...psi6,
ypsilon1...ypsilon6; {\tt rbwxxx.dau} are plot files of the
unrenormalized BW plot data.
For check: x-value rasters: preset XRdat-nonres.f written to fort.102,
with resonance grids merged written to fort.101 as actual raster
including resonance regions should agree with xRdat-extended.

resonances1n.f writes BW resonance shapes in FORTRAN data tables
{\tt rbwxxx.dat} where xxx=name of resonances omega, phi, psi1...psi6,
ypsilon1...ypsilon6

resonances_count.f writes physical and bare BW data renormalized via
the alphaQEDc26 routines of complex alpha(s) on preset raster
xRdat-extended.f; includes function BWRENOfac(s,fracerror) complex alpha
renormalization via alphaQEDc26; output:

plotfiles:
fort.210 for omega, fort.220 for phi, and  
fort.211...216 for psi1...psi6; fort.221...226 for ypsilon1...ypsilon6

FORTRAN data arrays:
fort.102 omega, fort.103 phi, fort.104 psi1, ....

all merge into reso_reno.dat by running via resonances_count.sh

gnuplot shows: reson0_gnu.sh for resnances0n.sh and reson1_gnu.sh for resonances_count.sh

Input and routines utlized: on main directory

resonances_dat.f -- input parameters

functions included in Rdat_fun.f:
BW(s,M,G,P) Breit-Wigner function providing physical (dressed) or bare versions
   depending on IREN switch: IREN=0,1,2 for physical,undressed by spacelike,
   or by timelike HVP subtraction.
BWRENO(s,fracerror) renormalization factor (vacuum polarization subtrsction)
   using tabels from reson_reno.dat
BWdressed(s,M,G,P) physical BW shape based on PDG parameters

