AAARRRRRRRGGGHHHHH!!!!!! More details to come.
using_rad = 1 ; (ONE = TRUE)
use_expon = 0 ; (LEAVE AT 0)
one_tail = -3 ; 0=all, 1=e, 2=e', 3=p, -3=all but p
intcor_mode = 1 ; (LEAVE AT 1)
Egamma_res_limit = 1. ; Egamma resolution limit (MeV)
spect_mode = 0 ; 0=e+p arms, -1=p arm, -2=e arm only, 1=none
cuts.Em.min = 0. ; (Em.min=Em.max=0.0 gives wide open cuts)
cuts.Em.max = 200. ; Must be wider than cuts in analysis(elastic or e,e'p)
using_Eloss = 1 ; (ONE = TRUE)
correct_Eloss = 1 ; ONE = correct reconstructed events for eloss.
correct_raster = 1 ; ONE = Reconstruct events using 'raster' matrix elements.
mc_smear = 1 ; ONE = target & hut mult scatt AND DC smearing.
deForest_flag = 0 ; 0=sigcc1, 1=sigcc2, -1=sigcc1 ONSHELL
rad_flag = 0 ; (radiative option #1...see init.f)
extrad_flag = 2 ; (rad. option #2...see init.f)
lambda(1) = 0.0 ; if rad_flag.eq.4 then lambda(1) = {TF}
lambda(2) = 0.0 ; if rad_flag.eq.4 then lambda(2) = {TF}
lambda(3) = 0.0 ; if rad_flag.eq.4 then lambda(3) = {TF}
Nntu = 1 ; ONE = generate ntuples
using_Coulomb = 1 ; (ONE = TRUE)
dE_edge_test = 0. ; (move around energy edges)
use_offshell_rad = 1 ; (ONE = TRUE)
Egamma_gen_max = 0. ; Set $>$0 to hardwire the Egamma limits.
---------