12 include
"mcRadCor.inc"
19 real phrad(4),q2true,nutrue,radweight
20 real nu,q2,
phi,yys,xxs
28 call
radgen(mcset_enebeam,q2,
nu,yys,xxs,
phi,phrad,q2true,nutrue,
40 mcradcor_nutrue=nutrue
41 mcradcor_q2true=q2true
42 mcradcor_ytrue=nutrue/mcset_enebeam
43 mcradcor_xtrue=q2true/mcradcor_ytrue/(4.*ebeam*pbeam)
44 mcradcor_w2true=amp2 + (q2true*(1./mcradcor_xtrue-1.))
49 mcradcor_ebrems=phrad(4)
50 mcradcor_thetabrems=0.
52 + mcradcor_thetabrems = acos(phrad(3)/phrad(4))
54 if (.not.(phrad(1).eq.0..and.phrad(2).eq.0.))
then
55 mcradcor_phibrems = atan2(phrad(2),phrad(1))
56 if (mcradcor_phibrems.lt.0.)
57 + mcradcor_phibrems = mcradcor_phibrems +
twopi
66 mcradcor_sigrad=sigrad
67 mcradcor_sigcor=sigcor
75 mcradcor_redfac=redfac
81 else if (ita.eq.3)
then