The document describes the use of fr2in.for It was written in 2011-12 by Alex Brown with the guidance of Ian Thompson. It is used for (p,t), (t,p), (p,3He), (3He,p) and (3He,n) reactions on targets with J=0 (even=even). Any final J is allowed. If you want to make a *.eps figure you also need plotsig.for and cps.for. All of these programs can be compiled with gfortran for.zip contains the source codes for fr2in, plotsig and cps ca40.zip contains the input and output files for the example below ------------------------------------------------------- Sample input file j0.dai for 40Ca 0+ to 38Ca 0+ bbcm00_co00.tna re = pt el = 35. The first line must be the name of the *.tna file from oxbash, nushell or nushellx. The following lines must all be in the format xx = yy re is the type of reaction. The options are pt, tp, ph and hp el is the laboratory projectile energy in MeV With no other input the program uses the following (default) internal optical potentials (for details see below): proton: CH deuteron: DG triton: XL helium: BG To run fr2in type fr2in < j0.dai > j0.dao fr2in makes the fresco input files j0t1.in (direct only) and j0t2.in (direct plus sequential) for fresco. The windows command file fr2.bat given below runs fr2in and fresco and then makes a plot. Other input options are documented in the *.dao file made with fr2in. These include the following given by example: nf = 2 for the second 0+ state (the default nf = 1) ef = 2. excitation energy (MeV) of the final state (the default is 0.) sp = 0. to change effective single-particle separation energy, the default is (bef-bei)/2 rm = 30. to change rmatch in fresco (in fm) pp = CH, BG, ME or PE optical potentials for protons dp = DG, PP or LH optical potentials for deuterons tp = XL, BG or PA optical potentials for tritons hp = BG or PA optical potentials for helium proton potentials [BG] Bechetti-Greenlees (A>40 2040 E>10 MeV):Phys Rep 201, 57 (1991) [ME] Menet (3040 827 12 af = 38 <-> 40 | ! | zi <-> zf = 20 <-> 20 | ! | 2ji <-> 2jf = 0 <-> 0 | ! | pi <-> pf = + <-> + | ! ---------------------------------- ! ! one-body transition densities a(dj) = ! = /sqrt(2dj+1) ! for input to some dwba codes: z = a(dj)/sqrt(2ji+1) ! with edmonds (de-shalit talmi) reduced matrix element convention ! ! for n,l,2j = 1 2 3 label k = 5 ! for n,l,2j = 1 2 5 label k = 4 ! for n,l,2j = 2 0 1 label k = 6 ! ! ji, jf, ti, tf, ti, fz, ! dj, ni, nf, ef, ei, exi, exf, ! 0.0, 0.0, 1.0, 0.0, -1.0, 0.0, 0.0, 1., 1., -279.961, -250.864, 0.000, 0.000, 4, 4, 0.00000, 0.16403, ! k1,k2,a(pp),a(nn) 5, 5, 0.00000, 0.97237, 6, 6, 0.00000, 0.16610, 0, 0.0, 2., 1., -279.961, -244.836, 6.028, 0.000, 4, 4, 0.00000, 0.07794, 5, 5, 0.00000, -0.18065, 6, 6, 0.00000, 0.98051, 0, -------------------------------------- fr2.bat (given below) is a windows command file that runs fr2in and fresco and then makes a figure -------------------------------------- rem - run fr2in to make the fresco inputs rem - *t1.in for direct rem - *t2.in for direct plus sequential fr2in < %1.dai > %1.dao rem ----- run fresco for *t1.in copy babt1.in %1t1.in del babt1.in fresco < %1t1.in > %1t1.out rem ----- copy the fresco output into the *.elas, *.sig files if exist fort.200 goto aa copy fort.201 %1t1.elas copy fort.202 %1t1.sig del fort.201 del fort.202 goto bb :aa copy fort.200 %1t1.elas copy fort.201+fort.202 %1t1.sig del fort.200 del fort.201 del fort.202 :bb copy fort.13 %1t1.tot rem ----- run fresco for *t2.in copy babt2.in %1t2.in del babt2.in fresco < %1t2.in > %1t2.out rem ----- copy the fresco output into the *.elas, *.sig files if exist fort.200 goto aa copy fort.201 %1t2.elas copy fort.202 %1t2.sig del fort.201 del fort.202 goto bb :aa copy fort.200 %1t2.elas copy fort.201+fort.202 %1t2.sig del fort.200 del fort.201 del fort.202 :bb copy fort.13 %1t2.tot rem ----- convert the *.sig to *.plt plotsig %1 rem ----- make *.eps with log y axis - you must have cps.exe to use this cps %1 > temp.out rem ----- make *.eps with decade y axis - you must have cps.exe to use this cps %1x > temp.out del temp.out rem ----- delete the uneeded file del fort.13 del fort.16 del fort.3 del fort.35 del fort.38 del fort.39 del fort.4 del fort.40 del fort.45 del fort.46 del fort.48 del fort.56 del fort.7 del fort.75