G4NuMuXQ2Integration.cc & G4NuMuXYIntegration.cc
================================================
Makes integration of d(sig)/dxdQ2 (DIS) for neutrino/antineutrino over x and then over Q2,
having as an intermediate result the d(sig)/dQ2 distribution and as a final result sig_tot.
Separation of QES (W<M_N+m_pi) is planned at this time.
  Output for 20 points equidistant from ln(E_min) to ln(E_max=300 (GeV)) (A=12):
  ------------------------------------------------------------------------------
  >>>>Emin=0.111651,lEi=-2.19238,lEa=5.70378,dlE=0.394808 (sig_t, 10^-38 cm^2 GeV^-1)
-----------------+----------+
I(Idx)dQ2,exp=.67| I(Idx)dy | ! NOT UP TO DATE !
-----------------+----------+
 E_nu   sig_t/E  | sig_t/E  |
-----------------+----------+
.136017 0.158940 | 0.188611 |
.201863 0.347557 | 0.354376 |
.299584 0.482431 | 0.484008 |
.444613 0.587544 | 0.587829 |
.659850 0.663442 | 0.663354 |
.979282 0.710895 | 0.710521 |
1.45335 0.734348 | 0.733853 |
2.15692 0.741216 | 0.740403 |
3.20108 0.737759 | 0.736918 |
4.75073 0.729512 | 0.728337 |
7.05055 0.718947 | 0.717549 |
10.4637 0.707357 | 0.705886 |
15.5292 0.695349 | 0.693841 |
23.0469 0.683188 | 0.681529 |
34.2038 0.671010 | 0.669100 |
50.7619 0.658915 | 0.656679 |
75.3357 0.646951 | 0.644425 |
111.806 0.635081 | 0.632419 |
165.931 0.622685 | 0.620626 |
246.258 0.610883 | 0.608874 |
-----------------+----------+
I(Idx)dQ2,exp=.34| I(Idx)dy |
-----------------+----------+
 E_anu  sig_t/E  | sig_t/E  |
-----------------+----------+
.136017 0.134604 | 0.162247 |
.201863 0.267721 | 0.274488 |
.299584 0.335296 | 0.337052 |
.444613 0.368757 | 0.369132 |
.659850 0.379554 | 0.379583 |
.979282 0.377011 | 0.376908 |
1.45335 0.368407 | 0.368181 |
2.15692 0.358362 | 0.357970 |
3.20108 0.349276 | 0.348610 |
4.75073 0.341937 | 0.341011 |
7.05055 0.336524 | 0.335304 |
10.4637 0.332704 | 0.331202 |
15.5292 0.330011 | 0.328298 |
23.0469 0.328097 | 0.326188 |
34.2038 0.326474 | 0.324666 |
50.7619 0.325437 | 0.323606 |
75.3357 0.324840 | 0.322993 |
111.806 0.324660 | 0.322805 |
165.931 0.324814 | 0.322974 |
246.258 0.325131 | 0.323362 |
-----------------+----------+
Conclusion (05.10.2005):
 a) at high energies sig(nu) and sig(anu) look to be close to the experimental values
 b) fL OFF increase less then by 10% -> the only hope is to modify C(Q2), which is bigger
 c) both waysof calculation give the same result (within the accuracy of calculation)(?)

















































































