satcorr

SOS satcorr parameterization 03/14/06

The parameterization of d(ssdelta), which is a function of ssxp_tar and spcentral, is shown below.

Replay conditions: default matrix, genable_sos_satcorr = 1, genable_hms_satcorr = 2000, Chungchen's angle corrections turned on.

This parameterization needs to be applied IN ADDITION TO the existing sos saturation correction (NB see important correction here)

s_satcorr.f
       ssdelta = ssdelta +
     > (0.0199-0.34515*ssxp_tar+12.8196*ssxp_tar**2+212.49192*ssxp_tar**3)
     > +(-0.0160165-0.277808*ssxp_tar-14.74545*ssxp_tar**2)*spcentral
     > +(-174.44034*exp(138.8557*(ssxp_tar-0.07550238))+0.1112637+6.1343792*ssxp_tar)
     > *(1.077034-1.56109*spcentral+0.387612*spcentral**3)*exp(0.475572*spcentral)

The plots of ssxptar vs. missing mass look much better now (new plots).

Hee'p setting Ebeam=4021 (electron in the sos, psos=-1.59 GeV)

H(e,e'pi) production data, electron in the sos, setting 1alh2, psos=-1.73 GeV