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>>>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)(?)