< Notations and definitions for output.dat > WRITE(6,*)'WHICH ONE PRINT TO OUT: DST, DSL, DSTT, DSLT, DSLTP, dSU, L/T_ratio, A_LU (1) OR DDS (3)' READ(*,*) OUT OUT=1 ---> You choose to calculate partial cross sections and SSA, e.g., ds_T, ds_L, ... ,ds_U, and A_LU, which are all printing out at the same time CASE(1) IF(XAXIS.EQ.8.OR.XAXIS.EQ.9) THEN !for choosing -t for the x-axis with inputs W & QQ or with inputs x_B & QQ. WRITE(40,*)REAL(-T)/10.**6,DST*PI/(K*P)*10.**6 !ds/dt WRITE(41,*)REAL(-T)/10.**6,DSL*PI/(K*P)*10.**6 WRITE(42,*)REAL(-T)/10.**6,DSTT*PI/(K*P)*10.**6 WRITE(43,*)REAL(-T)/10.**6,DSLT*PI/(K*P)*10.**6 WRITE(44,*)REAL(-T)/10.**6,DSLTP*PI/(K*P)*10.**6 WRITE(46,*)-REAL(T)/10.**6,DSU*PI/(K*P)*10.**6 WRITE(32,*)-REAL(T)/10.**6,DST !ds/d(Omega) WRITE(33,*)-REAL(T)/10.**6,DSL WRITE(34,*)REAL(-T)/10.**6,DSTT WRITE(35,*)REAL(-T)/10.**6,DSLT WRITE(36,*)REAL(-T)/10.**6,DSLTP WRITE(37,*)-REAL(T)/10.**6,DSU WRITE(50,*)REAL(-T+T0)/10.**6,DST*PI/(K*P)*10.**6 !ds/dt' ( t'=-(t-t_0) ) WRITE(51,*)REAL(-T+T0)/10.**6,DSL*PI/(K*P)*10.**6 WRITE(52,*)REAL(-T+T0)/10.**6,DSTT*PI/(K*P)*10.**6 WRITE(53,*)REAL(-T+T0)/10.**6,DSLT*PI/(K*P)*10.**6 WRITE(54,*)REAL(-T+T0)/10.**6,DSLTP*PI/(K*P)*10.**6 WRITE(56,*)-REAL(T-T0)/10.**6,DSU*PI/(K*P)*10.**6 WRITE(27,*)-REAL(T)/10.**6,DSL/DST WRITE(48,*)-REAL(T)/10.**6,ALU ELSEIF(XAXIS.EQ.10) THEN ! for choosing cos(theta) as the x-axis with inputs W & QQ or x_B & QQ. WRITE(40,*)COS(THETA),DST*PI/(K*P)*10.**6 WRITE(41,*)COS(THETA),DSL*PI/(K*P)*10.**6 WRITE(42,*)COS(THETA),DSTT*PI/(K*P)*10.**6 WRITE(43,*)COS(THETA),DSLT*PI/(K*P)*10.**6 WRITE(44,*)COS(THETA),DSLTP*PI/(K*P)*10.**6 WRITE(46,*)COS(THETA),DSU*PI/(K*P)*10.**6 WRITE(32,*)COS(THETA),DST WRITE(33,*)COS(THETA),DSL WRITE(34,*)COS(THETA),DSTT WRITE(35,*)COS(THETA),DSLT WRITE(36,*)COS(THETA),DSLTP WRITE(37,*)COS(THETA),DSU WRITE(27,*)COS(THETA),DSL/DST WRITE(48,*)COS(THETA),ALU ELSEIF(XAXIS.EQ.1.OR.XAXIS.EQ.6) THEN !for choosing the total energy W as the x-axis with inputs QQ & THETA. WRITE(40,*)E/10.**3,DST*PI/(K*P)*10.**6 !E=W=total energy, invariant energy WRITE(41,*)E/10.**3,DSL*PI/(K*P)*10.**6 WRITE(42,*)E/10.**3,DSTT*PI/(K*P)*10.**6 WRITE(43,*)E/10.**3,DSLT*PI/(K*P)*10.**6 WRITE(44,*)E/10.**3,DSLTP*PI/(K*P)*10.**6 WRITE(46,*)E/10.**3,DSU*PI/(K*P)*10.**6 WRITE(32,*)E/10.**3,DST WRITE(33,*)E/10.**3,DSL WRITE(34,*)E/10.**3,DSTT WRITE(35,*)E/10.**3,DSLT WRITE(36,*)E/10.**3,DSLTP WRITE(37,*)E/10.**3,DSU WRITE(27,*)E/10.**3,DSL/DST WRITE(48,*)E/10.**3,ALU ELSE !for choosing Q**2 as the x-axis with inputs W & FIXED THETA OR -T WRITE(40,*)QQ/10.**6,DST*PI/(K*P)*10.**6 WRITE(41,*)QQ/10.**6,DSL*PI/(K*P)*10.**6 WRITE(42,*)QQ/10.**6,DSTT*PI/(K*P)*10.**6 WRITE(43,*)QQ/10.**6,DSLT*PI/(K*P)*10.**6 WRITE(44,*)QQ/10.**6,DSLTP*PI/(K*P)*10.**6 WRITE(46,*)QQ/10.**6,DSU*PI/(K*P)*10.**6 WRITE(32,*)QQ/10.**6,DST WRITE(33,*)QQ/10.**6,DSL WRITE(34,*)QQ/10.**6,DSTT WRITE(35,*)QQ/10.**6,DSLT WRITE(36,*)QQ/10.**6,DSLTP WRITE(37,*)QQ/10.**6,DSU WRITE(27,*)QQ/10.**6,DSL/DST WRITE(48,*)QQ/10.**6,ALU ENDIF WRITE(6,*)'WHICH ONE PRINT TO OUT: DST, DSL, DSTT, DSLT, DSLTP, dSU, L/T_ratio, A_LU (1) OR DDS (3)' READ(*,*) OUT OUT=3 ---> You choose to calculate the total cross section 2pi*dds/dtd(phi) CASE(3) !Then you choose the type of the DDS below for the case of the PHI angle given; IF(DDS_TYPE.EQ.1) THEN !DESY IF(GTYPE.EQ.1.OR.GTYPE.EQ.3.OR.GTYPE.EQ.4) THEN !Fixed Phi=225deg FOR PROTON TARGET DS=DSU+EL*DSTT*COS(2.*PHI*PI/180)+SQRT(2.*EL*(1.+EL))*DSLT*COS(PHI*PI/180) ELSEIF(GTYPE.EQ.2) THEN !phi=120-240deg INTERATED FOR NEUTRON TARGET DS=(DSU*2.*PI/3.+EL*DSTT*(0.866)+SQRT(2.*EL*(1.+EL))*DSLT*(-1.7321))*3./2./PI !best fit ENDIF ELSEIF(DDS_TYPE.EQ.2) THEN !Cornell fixed DS=DSU+EL*DSTT*COS(2.*PHI*PI/180)+SQRT(0.5*EL*(1.+EL))*DSLT*COS(PHI*PI/180) ELSEIF(DDS_TYPE.EQ.3) THEN !Cornell average phi=-45, 45deg INTEGRATED. DS=(DSU*PI/2.+EL*DSTT*(1.)+SQRT(0.5*EL*(1.+EL))*DSLT*SQRT(2.))*2./PI ELSEIF(DDS_TYPE.EQ.4) THEN !CEA DS=DSU+EL*DSTT*COS_2PHI_INPUT+SQRT(0.5*EL*(1.+EL))*DSLT*COS_PHI_INPUT !CEA Brown !JLab ELSEIF(DDS_TYPE.EQ.5) THEN !JLab DS=DSU+EL*DSTT*COS(2.*PHI*PI/180)+SQRT(2.*EL*(1.+EL))*DSLT*COS(PHI*PI/180) ENDIF IF(XAXIS3.EQ.5) THEN !Choose -t as the x-axis for Phi angle given WRITE(45,*)-REAL(T)/10.D0**6,DS*PI/(K*P)*10.D0**6 !ddS_T.dat WRITE(49,*)-REAL(T)/10.D0**6,DS !ddS_W(omega).dat ELSEIF(XAXIS3.EQ.6) THEN !Choose W as the x-axis with inputs COS_2*PHI, COS_PHI for CEA by Brown !DS=DST+EL*DSL+EL*DSTT*COS_2PHI_INPUT+SQRT(0.5*EL*(1.+EL))*DSLT*COS_PHI_INPUT WRITE(49,*)E/10.D0**3,DS !ddS_W.dat ELSEIF(XAXIS3.EQ.7) THEN !Choose Q**2 as the x-axis with inputs COS_2*PHI, COS_PHI for CEA by Brown !DS=DST+EL*DSL+EL*DSTT*COS_2PHI_INPUT+SQRT(0.5*EL*(1.+EL))*DSLT*COS_PHI_INPUT WRITE(49,*)QQ/10.D0**6,DS !ddS_QQ.dat !WRITE(59,*)QQ/10.D0**6,DS !ddS_QQ.dat !For the case of choosing the PHI-angle as the x-axis, ELSE DO I10=1,360,1 PHI=I10*PI/180. IF(DDS_TYPE.EQ.1.OR.DDS_TYPE.EQ.5) THEN DS=DSU+EL*DSTT*COS(2.*PHI)+SQRT(2*EL*(1.+EL))*DSLT*COS(PHI) ELSE DS=DSU+EL*DSTT*COS(2.*PHI)+SQRT(0.5*EL*(1.+EL))*DSLT*COS(PHI) ENDIF WRITE(47,*)I10,DS*PI/(K*P)*10.**6 !ddS_PHI.dat WRITE(49,*)I10,DS !ddS_PHI.dat END DO END IF