First, the "integral" folder should be with the epion-tch-v2.f90 file in the same folder.  

Select channel

1 (is selected for demonstration) 


Choose the type of nucleon form factor between  GPD parameterization (1) Conventional dipole form factor (2)

1 (is selected)

If (1) is chosen, input parameter alpha' is required; 
alpha'=1.105 yields a similar result from the nucleon dipole FF of Lambda=0.843 GeV,
We adopt alpha'=0.3 which corresponds to the result from the dipole FF of Lambda=1.55 GeV

0.3 (is selected)

- The next option asks on/off of nucleon magnetic FF F_2; 
  
1 (is selected. In the GPD FF we set on the F_2.)

- The next option needs the cutoff Lambda of pion charge FF of the monopole type;
  For DESY data Lambda=710 MeV, 
  For JLab Horn data, Lambda=680 MeV
  For Cornell, CEA, & JLab Lambda=755 MeV 

775 (is selected)
 
- The next asks the sign of the axial meson a_1(1260) which affects a change meaningful to the shape of differential and beam SSA.
  In Alicia's experiment we used the +sign of a_1. But, -sign of a_1 better improves the differential and SSA.  

2 (is selected)

Up to these options we choose the parameters for nucleon and pion FF with the sign of a_1.
In the following we choose the kinematic inputs (1) for components of differential cross sections, the ratio of ds_L/ds_T, AU_L (3) for total cross section  

1 (is selected)

In case of choosing (1);

The range of pion production angle theta is selected as small (below 60 deg), or full (180 deg)

1 (is selected)


 CHOOSE X-AXIS & FIX VARs. : W,QQ&t(1), QQ,W&t(2), QQ,xb&t(3), QQ,xb&theta(4), QQ,W&theta(5)
    W,QQ&cos(theta)(6), QQ,W&cos(theta)(7),  -t,xb&QQ(8) ,  -t,W&QQ(9),  OR  cos(theta),W&QQ(10)

1 (is selected)

**Then, choose kinematic variables; if W, QQ&t(1) is selected, for instance, you get the cross section as y-axis vs. total energy W as the x-axis with Q**2=? and -t=? as the inputs.**

INPUT: FINAL W(GEV) 

4 (is selected) -->   Select the last value on the x-axis with energy W=4 GeV 

INPUT: QQ(GEV^2)

2.9 (is selected) --> Q**2=2.9 GeV^2

INPUT: -T(GEV^2)   (These numbers are automatically estimated as the W & QQ are input.)

0.75 (is selected) --> e=0.75

INPUT: PHOTON POL EL(0<EL<1)?
Separated dcs needs not PHOTON POL EL(0<EL<1) and PUT ANLY VALUE, e.g. 0.5, IF NOT RELEVANT

0.45 (is selected) 
**photon polarization epsilon is required only when ds_U and total cross section are calculated. 
So the case is irrelevant to these cross sections, say, the components of diff. crs sections, you may put any value, e.g., 0.5.**    

OUTPUT UNIT: MICROBARN(1), NANOBARN(2)

1 (is selected)
 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

In the option for nucleon FF above1

2. Choose the type of nucleon form factor between  GPD parameterization (1) Conventional dipole form factor (2)