Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2025-01-18 09:15:03

0001 #include <iostream>
0002 #include <fstream>
0003 #include <stdio.h>
0004 #include <cstdlib>
0005 #include <stdlib.h>
0006 #include <iomanip>
0007 #include <stdio.h>
0008 #include <math.h>
0009 #include <cmath>
0010 
0011 using std::setw;
0012 using std::setprecision;
0013 using std::cout;
0014 using std::cin;
0015 using std::endl;
0016 using namespace std;
0017 
0018 
0019 double MySigmaL( double qsqgev, double tgev, double wgev ) { 
0020 
0021    double s =  0.0;
0022    /* double wgev   = par[0]; */
0023    /* double qsqgev = par[1]; */
0024    /* double tgev   = x[0]; */
0025 
0026    //---------------------------------------------------------------------------------------------//
0027    //                                                                                             //
0028    //                                                                                             //
0029    //                                                                                             //
0030    //                                                                                             //
0031    //                                                                                             //
0032    //                                                                                             //
0033    //                                                                                             //
0034    //                                                                                             //
0035    //                                                                                             //
0036    //                                                                                             //
0037    //                                                                                             //
0038    //                                                                                             //
0039    //---------------------------------------------------------------------------------------------//
0040    //---------------------------------------------------------------------------------------------
0041    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2
0042    //---------------------------------------------------------------------------------------------
0043    if ( ( qsqgev <= 2 ) )  { 
0044      if ( ( wgev <= 2.1 ) )  { 
0045        s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0046      ( exp( 3.664 + qsqgev * -0.6073 + tgev * ( -16.22 + qsqgev * 3.759 ) ) + 
0047        exp( -4.753 + pow(qsqgev,1) * 11.27 + pow(qsqgev,2) * -4.634 + tgev * ( 4.881 + pow(qsqgev,1) * -16.22 + pow(qsqgev,2) * 6.862 ) ) );
0048      }
0049    }
0050    //---------------------------------------------------------------------------------------------
0051 
0052 
0053    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.1
0054    //---------------------------------------------------------------------------------------------
0055    if ( ( qsqgev <= 2 ) )  { 
0056      if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0057 
0058        s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0059          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0060      ( exp( 3.311 + qsqgev * -0.489 + 
0061         tgev * ( -16.94 + qsqgev * 2.873 ) ) + exp( -6.04 + pow(qsqgev,1) * 13.43 + pow(qsqgev,2) * -5.412 + 
0062                                 tgev * ( 5.984 + pow(qsqgev,1) * -19.85 + pow(qsqgev,2) * 8.324 ) ) );
0063 
0064      }
0065    }
0066    //---------------------------------------------------------------------------------------------
0067 
0068 
0069    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.2
0070    //---------------------------------------------------------------------------------------------
0071    if ( ( qsqgev <= 2 ) )  { 
0072      if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0073 
0074        s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0075          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0076      ( exp( 2.988 + qsqgev * -0.3947 + 
0077         tgev * ( -16.72 + qsqgev * 1.332 ) ) + exp( -7.208 + pow(qsqgev,1) * 14.99 + pow(qsqgev,2) * -5.879 + 
0078                                 tgev * ( 7.025 + pow(qsqgev,1) * -22.66 + pow(qsqgev,2) * 9.306 ) ) );
0079 
0080      }
0081    }
0082    //---------------------------------------------------------------------------------------------
0083 
0084 
0085    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.3
0086    //---------------------------------------------------------------------------------------------
0087    if ( ( qsqgev <= 2 ) )  { 
0088      if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0089 
0090        s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0091          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0092      ( exp( 2.702 + qsqgev * -0.327 + 
0093         tgev * ( -15.71 + qsqgev * -0.673 ) ) + exp( -8.252 + pow(qsqgev,1) * 16.09 + pow(qsqgev,2) * -6.124 + 
0094                                  tgev * ( 7.913 + pow(qsqgev,1) * -24.58 + pow(qsqgev,2) * 9.823 ) ) );
0095 
0096      }
0097    }
0098    //---------------------------------------------------------------------------------------------
0099 
0100 
0101    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.4
0102    //---------------------------------------------------------------------------------------------
0103    if ( ( qsqgev <= 2 ) )  { 
0104      if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0105 
0106        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0107          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0108      ( exp( 2.467 + qsqgev * -0.2932 + 
0109         tgev * ( -14 + qsqgev * -2.955 ) ) + exp( -9.39 + pow(qsqgev,1) * 16.6 + pow(qsqgev,2) * -6.004 + 
0110                               tgev * ( 8.998 + pow(qsqgev,1) * -24.97 + pow(qsqgev,2) * 9.404 ) ) );
0111 
0112      }
0113    }
0114    //---------------------------------------------------------------------------------------------
0115 
0116 
0117    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.5
0118    //---------------------------------------------------------------------------------------------
0119    if ( ( qsqgev <= 2 ) )  { 
0120      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0121 
0122        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0123          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0124      ( exp( 2.234 + qsqgev * -0.2634 + 
0125         tgev * ( -12.53 + qsqgev * -4.85 ) ) + exp( -10.29 + pow(qsqgev,1) * 17.08 + pow(qsqgev,2) * -5.98 + 
0126                                 tgev * ( 9.618 + pow(qsqgev,1) * -25.49 + pow(qsqgev,2) * 9.26 ) ) );
0127 
0128      }
0129    }
0130    //---------------------------------------------------------------------------------------------
0131 
0132 
0133    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.6
0134    //---------------------------------------------------------------------------------------------
0135    if ( ( qsqgev <= 2 ) )  { 
0136      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0137 
0138        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0139          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0140      ( exp( 2.038 + qsqgev * -0.2547 + 
0141         tgev * ( -10.93 + qsqgev * -6.522 ) ) + exp( -11.14 + pow(qsqgev,1) * 16.89 + pow(qsqgev,2) * -5.599 + 
0142                                  tgev * ( 10.12 + pow(qsqgev,1) * -24.34 + pow(qsqgev,2) * 8.237 ) ) );
0143 
0144      }
0145    }
0146    //---------------------------------------------------------------------------------------------
0147 
0148 
0149    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.7
0150    //---------------------------------------------------------------------------------------------
0151    if ( ( qsqgev <= 2 ) )  { 
0152      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0153 
0154        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0155          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0156      ( exp( 1.849 + qsqgev * -0.2489 + 
0157         tgev * ( -9.653 + qsqgev * -7.727 ) ) + exp( -11.85 + pow(qsqgev,1) * 16.67 + pow(qsqgev,2) * -5.259 + 
0158                                  tgev * ( 10.44 + pow(qsqgev,1) * -23.56 + pow(qsqgev,2) * 7.519 ) ) );
0159 
0160      }
0161    }
0162    //---------------------------------------------------------------------------------------------
0163 
0164 
0165    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.8
0166    //---------------------------------------------------------------------------------------------
0167    if ( ( qsqgev <= 2 ) )  { 
0168      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0169 
0170        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0171          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0172      ( exp( 1.68 + qsqgev * -0.2521 + 
0173         tgev * ( -8.554 + qsqgev * -8.54 ) ) + exp( -12.41 + pow(qsqgev,1) * 15.78 + pow(qsqgev,2) * -4.597 + 
0174                                 tgev * ( 10.49 + pow(qsqgev,1) * -21.37 + pow(qsqgev,2) * 6.119 ) ) );
0175 
0176      }
0177    }
0178    //---------------------------------------------------------------------------------------------
0179 
0180 
0181    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 2.9
0182    //---------------------------------------------------------------------------------------------
0183    if ( ( qsqgev <= 2 ) )  { 
0184      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0185 
0186        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0187          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0188      ( exp( 1.515 + qsqgev * -0.2524 + 
0189         tgev * ( -7.748 + qsqgev * -8.954 ) ) + exp( -12.86 + pow(qsqgev,1) * 14.71 + pow(qsqgev,2) * -3.879 + 
0190                                  tgev * ( 10.45 + pow(qsqgev,1) * -19.41 + pow(qsqgev,2) * 4.911 ) ) );
0191 
0192      }
0193    }
0194    //---------------------------------------------------------------------------------------------
0195 
0196 
0197    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 3
0198    //---------------------------------------------------------------------------------------------
0199    if ( ( qsqgev <= 2 ) )  { 
0200      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0201 
0202        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0203          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0204      ( exp( 1.358 + qsqgev * -0.2512 + 
0205         tgev * ( -7.204 + qsqgev * -8.984 ) ) + exp( -13.13 + pow(qsqgev,1) * 13.06 + pow(qsqgev,2) * -2.906 + 
0206                                  tgev * ( 10.17 + pow(qsqgev,1) * -16.58 + pow(qsqgev,2) * 3.327 ) ) );
0207 
0208      }
0209    }
0210    //---------------------------------------------------------------------------------------------
0211 
0212 
0213    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 3.1
0214    //---------------------------------------------------------------------------------------------
0215    if ( ( qsqgev <= 2 ) )  { 
0216      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0217 
0218        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0219          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0220      ( exp( 1.203 + qsqgev * -0.2438 + 
0221         tgev * ( -6.909 + qsqgev * -8.714 ) ) + exp( -13.29 + pow(qsqgev,1) * 11.19 + pow(qsqgev,2) * -1.85 + 
0222                                  tgev * ( 9.888 + pow(qsqgev,1) * -14.19 + pow(qsqgev,2) * 2.024 ) ) );
0223 
0224      }
0225    }
0226    //---------------------------------------------------------------------------------------------
0227 
0228 
0229    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 3.2
0230    //---------------------------------------------------------------------------------------------
0231    if ( ( qsqgev <= 2 ) )  { 
0232      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0233 
0234        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0235          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0236      ( exp( 1.047 + qsqgev * -0.228 + 
0237         tgev * ( -6.859 + qsqgev * -8.162 ) ) + exp( -14.74 + pow(qsqgev,1) * 11.49 + pow(qsqgev,2) * -1.638 + 
0238                                  tgev * ( 6.205 + pow(qsqgev,1) * -5.893 + pow(qsqgev,2) * -1.5 ) ) );
0239 
0240      }
0241    }
0242    //---------------------------------------------------------------------------------------------
0243 
0244 
0245    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 3.3
0246    //---------------------------------------------------------------------------------------------
0247    if ( ( qsqgev <= 2 ) )  { 
0248      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
0249 
0250        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0251          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0252      ( exp( 0.891 + qsqgev * -0.2052 + 
0253         tgev * ( -6.992 + qsqgev * -7.436 ) ) + exp( -16.86 + pow(qsqgev,1) * 12.36 + pow(qsqgev,2) * -1.541 + 
0254                                  tgev * ( 1.811 + pow(qsqgev,1) * 2.739 + pow(qsqgev,2) * -4.907 ) ) );
0255 
0256      }
0257    }
0258    //---------------------------------------------------------------------------------------------
0259 
0260 
0261    //Paramaterization fo Sigma_L for Qsq: 0.1  to  2  at w = 3.4
0262    //---------------------------------------------------------------------------------------------
0263    if ( ( qsqgev <= 2 ) )  { 
0264      if ( ( wgev >= 3.4 ) )  { 
0265 
0266        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0267          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0268      ( exp( 0.7328 + qsqgev * -0.1743 + 
0269         tgev * ( -7.255 + qsqgev * -6.599 ) ) + exp( -15.92 + pow(qsqgev,1) * 8.21 + pow(qsqgev,2) * 0.3937 + 
0270                                  tgev * ( -1.547 + pow(qsqgev,1) * 9.437 + pow(qsqgev,2) * -7.492 ) ) );
0271 
0272      }
0273    }
0274    //---------------------------------------------------------------------------------------------
0275 
0276 
0277    //---------------------------------------------------------------------------------------------//
0278    //                                                                                             //
0279    //                                                                                             //
0280    //                                                                                             //
0281    //                                                                                             //
0282    //                                                                                             //
0283    //                                                                                             //
0284    //                                                                                             //
0285    //                                                                                             //
0286    //                                                                                             //
0287    //                                                                                             //
0288    //                                                                                             //
0289    //                                                                                             //
0290    //---------------------------------------------------------------------------------------------//
0291    //---------------------------------------------------------------------------------------------
0292    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2
0293    //---------------------------------------------------------------------------------------------
0294    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0295      if ( wgev <= 2.1 )  { 
0296 
0297        s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0298          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0299      ( exp( 3.797 + qsqgev * -0.7708 + 
0300         tgev * ( -11.4 + qsqgev * 2.043 ) ) + exp( 2.818 + pow(qsqgev,1) * -1.488 + pow(qsqgev,2) * 0.1288 + 
0301                                tgev * ( -5.239 + pow(qsqgev,1) * 2.237 + pow(qsqgev,2) * -0.2923 ) ) );
0302 
0303      }
0304    }
0305    //---------------------------------------------------------------------------------------------
0306 
0307 
0308    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.1
0309    //---------------------------------------------------------------------------------------------
0310    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0311      if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0312 
0313        s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0314          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0315      ( exp( 3.638 + qsqgev * -0.7541 + 
0316         tgev * ( -13.35 + qsqgev * 2.48 ) ) + exp( 3.451 + pow(qsqgev,1) * -1.955 + pow(qsqgev,2) * 0.2047 + 
0317                                tgev * ( -6.971 + pow(qsqgev,1) * 3.092 + pow(qsqgev,2) * -0.4229 ) ) );
0318 
0319      }
0320    }
0321    //---------------------------------------------------------------------------------------------
0322 
0323 
0324    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.2
0325    //---------------------------------------------------------------------------------------------
0326    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0327      if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0328 
0329        s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0330          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0331      ( exp( 3.482 + qsqgev * -0.738 + 
0332         tgev * ( -16.03 + qsqgev * 3.16 ) ) + exp( 3.698 + pow(qsqgev,1) * -2.017 + pow(qsqgev,2) * 0.189 + 
0333                                tgev * ( -8.654 + pow(qsqgev,1) * 3.772 + pow(qsqgev,2) * -0.5019 ) ) );
0334 
0335      }
0336    }
0337    //---------------------------------------------------------------------------------------------
0338 
0339 
0340    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.3
0341    //---------------------------------------------------------------------------------------------
0342    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0343      if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0344 
0345        s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0346          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0347      ( exp( 3.285 + qsqgev * -0.7044 + 
0348         tgev * ( -19.27 + qsqgev * 3.977 ) ) + exp( 4.278 + pow(qsqgev,1) * -2.353 + pow(qsqgev,2) * 0.2382 + 
0349                                 tgev * ( -10.97 + pow(qsqgev,1) * 4.897 + pow(qsqgev,2) * -0.6705 ) ) );
0350 
0351      }
0352    }
0353    //---------------------------------------------------------------------------------------------
0354 
0355 
0356    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.4
0357    //---------------------------------------------------------------------------------------------
0358    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0359      if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0360 
0361        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0362          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0363      ( exp( 3.072 + qsqgev * -0.6655 + 
0364         tgev * ( -23.49 + qsqgev * 5.085 ) ) + exp( 4.402 + pow(qsqgev,1) * -2.327 + pow(qsqgev,2) * 0.2203 + 
0365                                 tgev * ( -12.97 + pow(qsqgev,1) * 5.712 + pow(qsqgev,2) * -0.7753 ) ) );
0366 
0367      }
0368    }
0369    //---------------------------------------------------------------------------------------------
0370 
0371 
0372    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.5
0373    //---------------------------------------------------------------------------------------------
0374    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0375      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0376 
0377        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0378          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0379      ( exp( 2.84 + qsqgev * -0.6208 + 
0380         tgev * ( -27.9 + qsqgev * 6.183 ) ) + exp( 3.524 + pow(qsqgev,1) * -1.573 + pow(qsqgev,2) * 0.07332 + 
0381                                tgev * ( -13.29 + pow(qsqgev,1) * 5.239 + pow(qsqgev,2) * -0.641 ) ) );
0382 
0383      }
0384    }
0385    //---------------------------------------------------------------------------------------------
0386 
0387 
0388    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.6
0389    //---------------------------------------------------------------------------------------------
0390    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0391      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0392 
0393        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0394          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0395      ( exp( 2.59 + qsqgev * -0.5701 + 
0396         tgev * ( -31.51 + qsqgev * 6.912 ) ) + exp( 2.71 + pow(qsqgev,1) * -0.9737 + pow(qsqgev,2) * -0.02709 + 
0397                                 tgev * ( -12.99 + pow(qsqgev,1) * 4.427 + pow(qsqgev,2) * -0.4684 ) ) );
0398 
0399      }
0400    }
0401    //---------------------------------------------------------------------------------------------
0402 
0403 
0404    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.7
0405    //---------------------------------------------------------------------------------------------
0406    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0407      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0408 
0409        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0410          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0411      ( exp( 2.326 + qsqgev * -0.5167 + 
0412         tgev * ( -33.95 + qsqgev * 7.176 ) ) + exp( 2.217 + pow(qsqgev,1) * -0.6905 + pow(qsqgev,2) * -0.05912 + 
0413                                 tgev * ( -12.77 + pow(qsqgev,1) * 3.818 + pow(qsqgev,2) * -0.3557 ) ) );
0414 
0415      }
0416    }
0417    //---------------------------------------------------------------------------------------------
0418 
0419 
0420    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.8
0421    //---------------------------------------------------------------------------------------------
0422    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0423      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0424 
0425        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0426          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0427      ( exp( 2.039 + qsqgev * -0.4559 + 
0428         tgev * ( -35.11 + qsqgev * 6.981 ) ) + exp( 1.796 + pow(qsqgev,1) * -0.4831 + pow(qsqgev,2) * -0.07818 + 
0429                                 tgev * ( -12.3 + pow(qsqgev,1) * 3.097 + pow(qsqgev,2) * -0.2287 ) ) );
0430 
0431      }
0432    }
0433    //---------------------------------------------------------------------------------------------
0434 
0435 
0436    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 2.9
0437    //---------------------------------------------------------------------------------------------
0438    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0439      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0440 
0441        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0442          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0443      ( exp( 1.761 + qsqgev * -0.4007 + 
0444         tgev * ( -35.32 + qsqgev * 6.47 ) ) + exp( 1.536 + pow(qsqgev,1) * -0.4128 + pow(qsqgev,2) * -0.07247 + 
0445                                tgev * ( -12.23 + pow(qsqgev,1) * 2.787 + pow(qsqgev,2) * -0.1901 ) ) );
0446 
0447      }
0448    }
0449    //---------------------------------------------------------------------------------------------
0450 
0451 
0452    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 3
0453    //---------------------------------------------------------------------------------------------
0454    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0455      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0456 
0457        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0458          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0459      ( exp( 1.479 + qsqgev * -0.3454 + 
0460         tgev * ( -34.8 + qsqgev * 5.795 ) ) + exp( 1.397 + pow(qsqgev,1) * -0.4326 + pow(qsqgev,2) * -0.05434 + 
0461                                tgev * ( -12.26 + pow(qsqgev,1) * 2.575 + pow(qsqgev,2) * -0.1656 ) ) );
0462 
0463      }
0464    }
0465    //---------------------------------------------------------------------------------------------
0466 
0467 
0468    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 3.1
0469    //---------------------------------------------------------------------------------------------
0470    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0471      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0472 
0473        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0474          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0475      ( exp( 1.23 + qsqgev * -0.3052 + 
0476         tgev * ( -33.13 + qsqgev * 4.742 ) ) + exp( 1.172 + pow(qsqgev,1) * -0.4012 + pow(qsqgev,2) * -0.04431 + 
0477                                 tgev * ( -11.94 + pow(qsqgev,1) * 2.188 + pow(qsqgev,2) * -0.12 ) ) );
0478 
0479      }
0480    }
0481    //---------------------------------------------------------------------------------------------
0482 
0483 
0484    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 3.2
0485    //---------------------------------------------------------------------------------------------
0486    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0487      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0488 
0489        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0490          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0491      ( exp( 0.9782 + qsqgev * -0.2629 + 
0492         tgev * ( -30.96 + qsqgev * 3.591 ) ) + exp( 0.7577 + pow(qsqgev,1) * -0.2057 + pow(qsqgev,2) * -0.07127 + 
0493                                 tgev * ( -11.23 + pow(qsqgev,1) * 1.496 + pow(qsqgev,2) * -0.006339 ) ) );
0494 
0495      }
0496    }
0497    //---------------------------------------------------------------------------------------------
0498 
0499 
0500    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 3.3
0501    //---------------------------------------------------------------------------------------------
0502    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0503      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
0504 
0505        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0506          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0507      ( exp( 0.7846 + qsqgev * -0.244 + 
0508         tgev * ( -27.97 + qsqgev * 2.235 ) ) + exp( 0.231 + pow(qsqgev,1) * 0.06126 + pow(qsqgev,2) * -0.1103 + 
0509                                 tgev * ( -10.13 + pow(qsqgev,1) * 0.5808 + pow(qsqgev,2) * 0.1404 ) ) );
0510 
0511      }
0512    }
0513    //---------------------------------------------------------------------------------------------
0514 
0515 
0516    //Paramaterization fo Sigma_L for Qsq: 2  to  3  at w = 3.4
0517    //---------------------------------------------------------------------------------------------
0518    if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0519      if ( wgev >= 3.4 )  { 
0520 
0521        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0522          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0523      ( exp( 0.5923 + qsqgev * -0.2236 + 
0524         tgev * ( -25.23 + qsqgev * 1.105 ) ) + exp( -0.3767 + pow(qsqgev,1) * 0.407 + pow(qsqgev,2) * -0.1697 + 
0525                                 tgev * ( -8.837 + pow(qsqgev,1) * -0.5062 + pow(qsqgev,2) * 0.3347 ) ) );
0526 
0527      }
0528    }
0529    //---------------------------------------------------------------------------------------------
0530 
0531 
0532 
0533    //---------------------------------------------------------------------------------------------//
0534    //                                                                                             //
0535    //                                                                                             //
0536    //                                                                                             //
0537    //                                                                                             //
0538    //                                                                                             //
0539    //                                                                                             //
0540    //                                                                                             //
0541    //                                                                                             //
0542    //                                                                                             //
0543    //                                                                                             //
0544    //                                                                                             //
0545    //                                                                                             //
0546    //---------------------------------------------------------------------------------------------//
0547    //---------------------------------------------------------------------------------------------
0548    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2
0549    //---------------------------------------------------------------------------------------------
0550    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0551      if ( wgev <= 2.1 )  { 
0552 
0553        s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0554          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0555      ( exp( 3.329 + qsqgev * -0.6126 + 
0556         tgev * ( -8.417 + qsqgev * 1.02 ) ) + exp( 2.116 + pow(qsqgev,1) * -1.009 + pow(qsqgev,2) * 0.04834 + 
0557                                tgev * ( -3.688 + pow(qsqgev,1) * 1.171 + pow(qsqgev,2) * -0.1098 ) ) );
0558 
0559      }
0560    }
0561    //---------------------------------------------------------------------------------------------
0562 
0563 
0564    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.1
0565    //---------------------------------------------------------------------------------------------
0566    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0567      if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0568 
0569        s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0570          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0571      ( exp( 3.163 + qsqgev * -0.5951 + 
0572         tgev * ( -9.325 + qsqgev * 1.108 ) ) + exp( 3.366 + pow(qsqgev,1) * -1.81 + pow(qsqgev,2) * 0.1664 + 
0573                                 tgev * ( -5.35 + pow(qsqgev,1) * 1.91 + pow(qsqgev,2) * -0.2088 ) ) );
0574 
0575      }
0576    }
0577    //---------------------------------------------------------------------------------------------
0578 
0579 
0580    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.2
0581    //---------------------------------------------------------------------------------------------
0582    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0583      if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0584 
0585        s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0586          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0587      ( exp( 3.09 + qsqgev * -0.6028 + 
0588         tgev * ( -11.15 + qsqgev * 1.441 ) ) + exp( 3.798 + pow(qsqgev,1) * -1.989 + pow(qsqgev,2) * 0.1772 + 
0589                                 tgev * ( -6.615 + pow(qsqgev,1) * 2.283 + pow(qsqgev,2) * -0.2394 ) ) );
0590 
0591      }
0592    }
0593    //---------------------------------------------------------------------------------------------
0594 
0595 
0596    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.3
0597    //---------------------------------------------------------------------------------------------
0598    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0599      if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0600 
0601        s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0602          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0603      ( exp( 2.577 + qsqgev * -0.4781 + 
0604         tgev * ( -9.33 + qsqgev * 0.6788 ) ) + exp( 4.136 + pow(qsqgev,1) * -2.551 + pow(qsqgev,2) * 0.3115 + 
0605                                 tgev * ( -7.668 + pow(qsqgev,1) * 2.924 + pow(qsqgev,2) * -0.3701 ) ) );
0606 
0607      }
0608    }
0609    //---------------------------------------------------------------------------------------------
0610 
0611 
0612    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.4
0613    //---------------------------------------------------------------------------------------------
0614    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0615      if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0616 
0617        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0618          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0619      ( exp( 2.722 + qsqgev * -0.5537 + 
0620         tgev * ( -13.4 + qsqgev * 1.683 ) ) + exp( 11.57 + pow(qsqgev,1) * -6.541 + pow(qsqgev,2) * 0.8329 + 
0621                                tgev * ( -16.17 + pow(qsqgev,1) * 7.235 + pow(qsqgev,2) * -0.9303 ) ) );
0622 
0623      }
0624    }
0625    //---------------------------------------------------------------------------------------------
0626 
0627 
0628    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.5
0629    //---------------------------------------------------------------------------------------------
0630    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0631      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0632 
0633        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0634          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0635      ( exp( 2.724 + qsqgev * -0.5834 + 
0636         tgev * ( -17.67 + qsqgev * 2.68 ) ) + exp( 7.375 + pow(qsqgev,1) * -3.822 + pow(qsqgev,2) * 0.4003 + 
0637                                tgev * ( -14.44 + pow(qsqgev,1) * 5.604 + pow(qsqgev,2) * -0.64 ) ) );
0638 
0639      }
0640    }
0641    //---------------------------------------------------------------------------------------------
0642 
0643 
0644    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.6
0645    //---------------------------------------------------------------------------------------------
0646    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0647      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0648 
0649        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0650          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0651      ( exp( 2.607 + qsqgev * -0.5766 + 
0652         tgev * ( -21.76 + qsqgev * 3.562 ) ) + exp( 0.2743 + pow(qsqgev,1) * 0.4692 + pow(qsqgev,2) * -0.2373 + 
0653                                 tgev * ( -8.353 + pow(qsqgev,1) * 1.509 + pow(qsqgev,2) * -0.01137 ) ) );
0654 
0655      }
0656    }
0657    //---------------------------------------------------------------------------------------------
0658 
0659 
0660    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.7
0661    //---------------------------------------------------------------------------------------------
0662    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0663      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0664 
0665        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0666          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0667      ( exp( 2.381 + qsqgev * -0.5375 + 
0668         tgev * ( -24.64 + qsqgev * 4.025 ) ) + exp( -2.131 + pow(qsqgev,1) * 1.836 + pow(qsqgev,2) * -0.422 + 
0669                                 tgev * ( -6.084 + pow(qsqgev,1) * -0.1815 + pow(qsqgev,2) * 0.2391 ) ) );
0670 
0671      }
0672    }
0673    //---------------------------------------------------------------------------------------------
0674 
0675 
0676    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.8
0677    //---------------------------------------------------------------------------------------------
0678    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0679      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0680 
0681        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0682          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0683      ( exp( 2.171 + qsqgev * -0.5044 + 
0684         tgev * ( -27.24 + qsqgev * 4.37 ) ) + exp( 2.613 + pow(qsqgev,1) * -1.006 + pow(qsqgev,2) * 0.006567 + 
0685                                tgev * ( -12.8 + pow(qsqgev,1) * 3.424 + pow(qsqgev,2) * -0.2841 ) ) );
0686 
0687      }
0688    }
0689    //---------------------------------------------------------------------------------------------
0690 
0691 
0692    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 2.9
0693    //---------------------------------------------------------------------------------------------
0694    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0695      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0696 
0697        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0698          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0699      ( exp( 1.999 + qsqgev * -0.4829 + 
0700         tgev * ( -30.2 + qsqgev * 4.82 ) ) + exp( 1.658 + pow(qsqgev,1) * -0.5098 + pow(qsqgev,2) * -0.05409 + 
0701                               tgev * ( -11.69 + pow(qsqgev,1) * 2.486 + pow(qsqgev,2) * -0.149 ) ) );
0702 
0703      }
0704    }
0705    //---------------------------------------------------------------------------------------------
0706 
0707 
0708    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 3
0709    //---------------------------------------------------------------------------------------------
0710    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0711      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0712 
0713        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0714          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0715      ( exp( 1.793 + qsqgev * -0.4533 + 
0716         tgev * ( -32.39 + qsqgev * 5.049 ) ) + exp( 1.01 + pow(qsqgev,1) * -0.2116 + pow(qsqgev,2) * -0.08537 + 
0717                                 tgev * ( -11.19 + pow(qsqgev,1) * 1.954 + pow(qsqgev,2) * -0.07677 ) ) );
0718 
0719      }
0720    }
0721    //---------------------------------------------------------------------------------------------
0722 
0723 
0724    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 3.1
0725    //---------------------------------------------------------------------------------------------
0726    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0727      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0728 
0729        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0730          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0731      ( exp( 1.571 + qsqgev * -0.42 + 
0732         tgev * ( -34.91 + qsqgev * 5.405 ) ) + exp( 1.182 + pow(qsqgev,1) * -0.3622 + pow(qsqgev,2) * -0.05847 + 
0733                                 tgev * ( -12.8 + pow(qsqgev,1) * 2.62 + pow(qsqgev,2) * -0.1682 ) ) );
0734 
0735      }
0736    }
0737    //---------------------------------------------------------------------------------------------
0738 
0739 
0740    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 3.2
0741    //---------------------------------------------------------------------------------------------
0742    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0743      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0744 
0745        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0746          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0747      ( exp( 1.317 + qsqgev * -0.3785 + 
0748         tgev * ( -34.7 + qsqgev * 5.003 ) ) + exp( 1.798 + pow(qsqgev,1) * -0.83 + pow(qsqgev,2) * 0.02141 + 
0749                                tgev * ( -14.01 + pow(qsqgev,1) * 3.196 + pow(qsqgev,2) * -0.265 ) ) );
0750 
0751      }
0752    }
0753    //---------------------------------------------------------------------------------------------
0754 
0755 
0756    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 3.3
0757    //---------------------------------------------------------------------------------------------
0758    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0759      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
0760 
0761        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0762          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0763      ( exp( 1.089 + qsqgev * -0.3455 + 
0764         tgev * ( -35.43 + qsqgev * 4.889 ) ) + exp( 0.6923 + pow(qsqgev,1) * -0.2664 + pow(qsqgev,2) * -0.05271 + 
0765                                 tgev * ( -11.5 + pow(qsqgev,1) * 1.583 + pow(qsqgev,2) * -0.04004 ) ) );
0766 
0767      }
0768    }
0769    //---------------------------------------------------------------------------------------------
0770 
0771 
0772    //Paramaterization fo Sigma_L for Qsq: 3  to  4  at w = 3.4
0773    //---------------------------------------------------------------------------------------------
0774    if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0775      if ( wgev >= 3.4 )  { 
0776 
0777        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0778          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0779      ( exp( 0.8895 + qsqgev * -0.3216 + 
0780         tgev * ( -36.38 + qsqgev * 4.903 ) ) + exp( 0.02066 + pow(qsqgev,1) * 0.05668 + pow(qsqgev,2) * -0.09543 + 
0781                                 tgev * ( -10.45 + pow(qsqgev,1) * 0.8147 + pow(qsqgev,2) * 0.06867 ) ) );
0782 
0783      }
0784    }
0785    //---------------------------------------------------------------------------------------------
0786 
0787 
0788    //---------------------------------------------------------------------------------------------//
0789    //                                                                                             //
0790    //                                                                                             //
0791    //                                                                                             //
0792    //                                                                                             //
0793    //                                                                                             //
0794    //                                                                                             //
0795    //                                                                                             //
0796    //                                                                                             //
0797    //                                                                                             //
0798    //                                                                                             //
0799    //                                                                                             //
0800    //                                                                                             //
0801    //---------------------------------------------------------------------------------------------//
0802    //---------------------------------------------------------------------------------------------
0803    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2
0804    //---------------------------------------------------------------------------------------------
0805    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0806      if ( wgev <= 2.1 )  { 
0807 
0808        s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0809          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0810      ( exp( 6.316 + qsqgev * -1.314 + 
0811         tgev * ( -11.25 + qsqgev * 1.708 ) ) + exp( 40.8 + pow(qsqgev,1) * -18.61 + pow(qsqgev,2) * 2.034 + 
0812                                 tgev * ( -46.21 + pow(qsqgev,1) * 20.54 + pow(qsqgev,2) * -2.302 ) ) );
0813 
0814      }
0815    }
0816    //---------------------------------------------------------------------------------------------
0817 
0818 
0819    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.1
0820    //---------------------------------------------------------------------------------------------
0821    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0822      if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0823 
0824        s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0825          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0826      ( exp( 3.49 + qsqgev * -0.6669 + 
0827         tgev * ( -10.47 + qsqgev * 1.354 ) ) + exp( 16.4 + pow(qsqgev,1) * -7.129 + pow(qsqgev,2) * 0.6894 + 
0828                                 tgev * ( -14.02 + pow(qsqgev,1) * 5.372 + pow(qsqgev,2) * -0.5381 ) ) );
0829 
0830      }
0831    }
0832    //---------------------------------------------------------------------------------------------
0833 
0834 
0835    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.2
0836    //---------------------------------------------------------------------------------------------
0837    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0838      if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0839 
0840        s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0841          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0842      ( exp( 3.175 + qsqgev * -0.6146 + 
0843         tgev * ( -10.9 + qsqgev * 1.323 ) ) + exp( -0.7067 + pow(qsqgev,1) * 0.4637 + pow(qsqgev,2) * -0.151 + 
0844                                tgev * ( -4.909 + pow(qsqgev,1) * 1.148 + pow(qsqgev,2) * -0.06541 ) ) );
0845 
0846      }
0847    }
0848    //---------------------------------------------------------------------------------------------
0849 
0850 
0851    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.3
0852    //---------------------------------------------------------------------------------------------
0853    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0854      if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0855 
0856        s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0857          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0858      ( exp( 3.24 + qsqgev * -0.6235 + 
0859         tgev * ( -14.17 + qsqgev * 1.798 ) ) + exp( -27.9 + pow(qsqgev,1) * 12.85 + pow(qsqgev,2) * -1.545 + 
0860                                 tgev * ( 17.48 + pow(qsqgev,1) * -9.275 + pow(qsqgev,2) * 1.116 ) ) );
0861 
0862      }
0863    }
0864    //---------------------------------------------------------------------------------------------
0865 
0866 
0867    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.4
0868    //---------------------------------------------------------------------------------------------
0869    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0870      if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0871 
0872        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0873          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0874      ( exp( 2.264 + qsqgev * -0.4265 + 
0875         tgev * ( -10.69 + qsqgev * 0.8795 ) ) + exp( -24.03 + pow(qsqgev,1) * 10.69 + pow(qsqgev,2) * -1.258 + 
0876                                  tgev * ( 16.76 + pow(qsqgev,1) * -8.826 + pow(qsqgev,2) * 1.034 ) ) );
0877 
0878      }
0879    }
0880    //---------------------------------------------------------------------------------------------
0881 
0882 
0883    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.5
0884    //---------------------------------------------------------------------------------------------
0885    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0886      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0887 
0888        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0889          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0890      ( exp( 1.827 + qsqgev * -0.3574 + 
0891         tgev * ( -8.378 + qsqgev * 0.2808 ) ) + exp( -0.9108 + pow(qsqgev,1) * -0.03561 + pow(qsqgev,2) * -0.02563 + 
0892                                  tgev * ( -3.982 + pow(qsqgev,1) * 0.5961 + pow(qsqgev,2) * -0.04425 ) ) );
0893 
0894      }
0895    }
0896    //---------------------------------------------------------------------------------------------
0897 
0898 
0899    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.6
0900    //---------------------------------------------------------------------------------------------
0901    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0902      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0903 
0904        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0905          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0906      ( exp( 1.878 + qsqgev * -0.3962 + 
0907         tgev * ( -9.64 + qsqgev * 0.4949 ) ) + exp( 1.547 + pow(qsqgev,1) * -1.173 + pow(qsqgev,2) * 0.09803 + 
0908                                 tgev * ( -5.869 + pow(qsqgev,1) * 1.292 + pow(qsqgev,2) * -0.1167 ) ) );
0909 
0910      }
0911    }
0912    //---------------------------------------------------------------------------------------------
0913 
0914 
0915    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.7
0916    //---------------------------------------------------------------------------------------------
0917    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0918      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0919 
0920        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0921          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0922      ( exp( 1.9 + qsqgev * -0.4238 + 
0923         tgev * ( -11.57 + qsqgev * 0.834 ) ) + exp( 3.609 + pow(qsqgev,1) * -2.043 + pow(qsqgev,2) * 0.1837 + 
0924                                 tgev * ( -8.755 + pow(qsqgev,1) * 2.348 + pow(qsqgev,2) * -0.2206 ) ) );
0925 
0926      }
0927    }
0928    //---------------------------------------------------------------------------------------------
0929 
0930 
0931    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.8
0932    //---------------------------------------------------------------------------------------------
0933    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0934      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0935 
0936        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0937          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0938      ( exp( 1.979 + qsqgev * -0.461 + 
0939         tgev * ( -16.62 + qsqgev * 1.817 ) ) + exp( 20.15 + pow(qsqgev,1) * -9.121 + pow(qsqgev,2) * 0.9355 + 
0940                                 tgev * ( -30.24 + pow(qsqgev,1) * 11.4 + pow(qsqgev,2) * -1.184 ) ) );
0941 
0942      }
0943    }
0944    //---------------------------------------------------------------------------------------------
0945 
0946 
0947    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 2.9
0948    //---------------------------------------------------------------------------------------------
0949    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0950      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0951 
0952        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0953          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0954      ( exp( 1.928 + qsqgev * -0.4668 + 
0955         tgev * ( -23.08 + qsqgev * 3.081 ) ) + exp( 25.13 + pow(qsqgev,1) * -10.96 + pow(qsqgev,2) * 1.097 + 
0956                                 tgev * ( -38.23 + pow(qsqgev,1) * 14.29 + pow(qsqgev,2) * -1.448 ) ) );
0957 
0958      }
0959    }
0960    //---------------------------------------------------------------------------------------------
0961 
0962 
0963    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 3
0964    //---------------------------------------------------------------------------------------------
0965    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0966      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0967 
0968        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0969          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0970      ( exp( 1.874 + qsqgev * -0.4722 + 
0971         tgev * ( -29.14 + qsqgev * 4.189 ) ) + exp( -2.064 + pow(qsqgev,1) * 1.461 + pow(qsqgev,2) * -0.309 + 
0972                                 tgev * ( -9.52 + pow(qsqgev,1) * 0.8933 + pow(qsqgev,2) * 0.08102 ) ) );
0973 
0974      }
0975    }
0976    //---------------------------------------------------------------------------------------------
0977 
0978 
0979    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 3.1
0980    //---------------------------------------------------------------------------------------------
0981    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0982      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0983 
0984        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0985          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0986      ( exp( 1.701 + qsqgev * -0.4519 + 
0987         tgev * ( -32.98 + qsqgev * 4.807 ) ) + exp( -22.96 + pow(qsqgev,1) * 10.86 + pow(qsqgev,2) * -1.358 + 
0988                                 tgev * ( 19.76 + pow(qsqgev,1) * -12.55 + pow(qsqgev,2) * 1.593 ) ) );
0989 
0990      }
0991    }
0992    //---------------------------------------------------------------------------------------------
0993 
0994 
0995    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 3.2
0996    //---------------------------------------------------------------------------------------------
0997    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0998      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0999 
1000        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1001          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1002      ( exp( 1.489 + qsqgev * -0.4227 + 
1003         tgev * ( -36.23 + qsqgev * 5.287 ) ) + exp( -23.54 + pow(qsqgev,1) * 11.05 + pow(qsqgev,2) * -1.37 + 
1004                                 tgev * ( 19.16 + pow(qsqgev,1) * -12.49 + pow(qsqgev,2) * 1.588 ) ) );
1005 
1006      }
1007    }
1008    //---------------------------------------------------------------------------------------------
1009 
1010 
1011    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 3.3
1012    //---------------------------------------------------------------------------------------------
1013    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1014      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1015 
1016        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1017          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1018      ( exp( 1.231 + qsqgev * -0.3836 + 
1019         tgev * ( -32.34 + qsqgev * 4.145 ) ) + exp( -1.517 + pow(qsqgev,1) * 0.757 + pow(qsqgev,2) * -0.1712 + 
1020                                 tgev * ( -7.345 + pow(qsqgev,1) * -0.2729 + pow(qsqgev,2) * 0.1655 ) ) );
1021 
1022      }
1023    }
1024    //---------------------------------------------------------------------------------------------
1025 
1026 
1027    //Paramaterization fo Sigma_L for Qsq: 4  to  5  at w = 3.4
1028    //---------------------------------------------------------------------------------------------
1029    if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1030      if ( wgev >= 3.4 )  { 
1031 
1032        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1033          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1034      ( exp( 1.033 + qsqgev * -0.3587 + 
1035         tgev * ( -33.35 + qsqgev * 4.115 ) ) + exp( -1.837 + pow(qsqgev,1) * 0.8157 + pow(qsqgev,2) * -0.1682 + 
1036                                 tgev * ( -8.106 + pow(qsqgev,1) * -0.02679 + pow(qsqgev,2) * 0.1305 ) ) );
1037 
1038      }
1039    }
1040    //---------------------------------------------------------------------------------------------
1041 
1042 
1043 
1044    //---------------------------------------------------------------------------------------------//
1045    //                                                                                             //
1046    //                                                                                             //
1047    //                                                                                             //
1048    //                                                                                             //
1049    //                                                                                             //
1050    //                                                                                             //
1051    //                                                                                             //
1052    //                                                                                             //
1053    //                                                                                             //
1054    //                                                                                             //
1055    //                                                                                             //
1056    //                                                                                             //
1057    //---------------------------------------------------------------------------------------------//
1058    //---------------------------------------------------------------------------------------------
1059    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2
1060    //---------------------------------------------------------------------------------------------
1061    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1062      if ( wgev <= 2.1 )  { 
1063 
1064        s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1065          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1066      ( exp( 1.427 + qsqgev * -0.5186 + 
1067         tgev * ( -2.977 + qsqgev * 0.266 ) ) + exp( 2.695 + pow(qsqgev,1) * -0.9802 + pow(qsqgev,2) * 0.04141 + 
1068                                 tgev * ( -4.763 + pow(qsqgev,1) * 1.058 + pow(qsqgev,2) * -0.07383 ) ) );
1069 
1070      }
1071    }
1072    //---------------------------------------------------------------------------------------------
1073 
1074 
1075    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.1
1076    //---------------------------------------------------------------------------------------------
1077    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1078      if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
1079 
1080        s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1081          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1082      ( exp( 12.1 + qsqgev * -2.383 + 
1083         tgev * ( -17.38 + qsqgev * 2.714 ) ) + exp( -48.11 + pow(qsqgev,1) * 16.37 + pow(qsqgev,2) * -1.44 + 
1084                                 tgev * ( 55.68 + pow(qsqgev,1) * -19.67 + pow(qsqgev,2) * 1.695 ) ) );
1085 
1086      }
1087    }
1088    //---------------------------------------------------------------------------------------------
1089 
1090 
1091    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.2
1092    //---------------------------------------------------------------------------------------------
1093    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1094      if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
1095 
1096        s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1097          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1098      ( exp( 12.02 + qsqgev * -2.295 + 
1099         tgev * ( -20.26 + qsqgev * 3.073 ) ) + exp( 213.6 + pow(qsqgev,1) * -78.36 + pow(qsqgev,2) * 7.078 + 
1100                                 tgev * ( -187.2 + pow(qsqgev,1) * 68.37 + pow(qsqgev,2) * -6.246 ) ) );
1101 
1102      }
1103    }
1104    //---------------------------------------------------------------------------------------------
1105 
1106 
1107    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.3
1108    //---------------------------------------------------------------------------------------------
1109    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1110      if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
1111 
1112        s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1113          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1114      ( exp( 2.214 + qsqgev * -0.4352 + 
1115         tgev * ( -7.533 + qsqgev * 0.5409 ) ) + exp( -3.134 + pow(qsqgev,1) * 0.7513 + pow(qsqgev,2) * -0.1125 + 
1116                                  tgev * ( -0.5375 + pow(qsqgev,1) * -0.386 + pow(qsqgev,2) * 0.05568 ) ) );
1117 
1118      }
1119    }
1120    //---------------------------------------------------------------------------------------------
1121 
1122 
1123    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.4
1124    //---------------------------------------------------------------------------------------------
1125    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1126      if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
1127 
1128        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1129          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1130      ( exp( 1.96 + qsqgev * -0.4031 + 
1131         tgev * ( -7.466 + qsqgev * 0.4516 ) ) + exp( -1.466 + pow(qsqgev,1) * 0.006673 + pow(qsqgev,2) * -0.03515 + 
1132                                  tgev * ( -1.846 + pow(qsqgev,1) * 0.05675 + pow(qsqgev,2) * 0.01095 ) ) );
1133 
1134      }
1135    }
1136    //---------------------------------------------------------------------------------------------
1137 
1138 
1139    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.5
1140    //---------------------------------------------------------------------------------------------
1141    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1142      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1143 
1144        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1145          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1146      ( exp( 2.47 + qsqgev * -0.5071 + 
1147         tgev * ( -12.36 + qsqgev * 1.243 ) ) + exp( 41.1 + pow(qsqgev,1) * -15.16 + pow(qsqgev,2) * 1.311 + 
1148                                 tgev * ( -39.67 + pow(qsqgev,1) * 13.42 + pow(qsqgev,2) * -1.173 ) ) );
1149 
1150      }
1151    }
1152    //---------------------------------------------------------------------------------------------
1153 
1154 
1155    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.6
1156    //---------------------------------------------------------------------------------------------
1157    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1158      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1159 
1160        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1161          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1162      ( exp( 2.613 + qsqgev * -0.5439 + 
1163         tgev * ( -16.32 + qsqgev * 1.855 ) ) + exp( 41.34 + pow(qsqgev,1) * -14.92 + pow(qsqgev,2) * 1.259 + 
1164                                 tgev * ( -38.92 + pow(qsqgev,1) * 12.7 + pow(qsqgev,2) * -1.078 ) ) );
1165 
1166      }
1167    }
1168    //---------------------------------------------------------------------------------------------
1169 
1170 
1171    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.7
1172    //---------------------------------------------------------------------------------------------
1173    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1174      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1175 
1176        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1177          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1178      ( exp( 2.53 + qsqgev * -0.5403 + 
1179         tgev * ( -18.41 + qsqgev * 2.122 ) ) + exp( -11.24 + pow(qsqgev,1) * 4.446 + pow(qsqgev,2) * -0.5165 + 
1180                                 tgev * ( -0.8509 + pow(qsqgev,1) * -1.504 + pow(qsqgev,2) * 0.23 ) ) );
1181 
1182      }
1183    }
1184    //---------------------------------------------------------------------------------------------
1185 
1186 
1187    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.8
1188    //---------------------------------------------------------------------------------------------
1189    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1190      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1191 
1192        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1193          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1194      ( exp( 3.903 + qsqgev * -0.8267 + 
1195         tgev * ( -28.39 + qsqgev * 3.974 ) ) + exp( -126.7 + pow(qsqgev,1) * 49.54 + pow(qsqgev,2) * -4.923 + 
1196                                 tgev * ( 92.64 + pow(qsqgev,1) * -38.13 + pow(qsqgev,2) * 3.806 ) ) );
1197 
1198      }
1199    }
1200    //---------------------------------------------------------------------------------------------
1201 
1202 
1203    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 2.9
1204    //---------------------------------------------------------------------------------------------
1205    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1206      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1207 
1208        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1209          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1210      ( exp( -0.3224 + qsqgev * -0.02778 + 
1211         tgev * ( 6.073 + qsqgev * -2.743 ) ) + exp( -13.23 + pow(qsqgev,1) * 3.455 + pow(qsqgev,2) * -0.2587 + 
1212                                 tgev * ( 3.159 + pow(qsqgev,1) * -1.551 + pow(qsqgev,2) * 0.07253 ) ) );
1213 
1214      }
1215    }
1216    //---------------------------------------------------------------------------------------------
1217 
1218 
1219    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 3
1220    //---------------------------------------------------------------------------------------------
1221    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1222      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
1223 
1224        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1225          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1226      ( exp( 1.86 + qsqgev * -0.4661 + 
1227         tgev * ( -20.37 + qsqgev * 2.315 ) ) + exp( 3.289 + pow(qsqgev,1) * -0.835 + pow(qsqgev,2) * -0.05204 + 
1228                                 tgev * ( -18.28 + pow(qsqgev,1) * 4.379 + pow(qsqgev,2) * -0.279 ) ) );
1229 
1230      }
1231    }
1232    //---------------------------------------------------------------------------------------------
1233 
1234 
1235    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 3.1
1236    //---------------------------------------------------------------------------------------------
1237    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1238      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
1239 
1240        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1241          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1242      ( exp( 1.638 + qsqgev * -0.4375 + 
1243         tgev * ( -20.3 + qsqgev * 2.215 ) ) + exp( -11.4 + pow(qsqgev,1) * 4.389 + pow(qsqgev,2) * -0.5162 + 
1244                                tgev * ( -3.89 + pow(qsqgev,1) * -0.8682 + pow(qsqgev,2) * 0.1893 ) ) );
1245 
1246      }
1247    }
1248    //---------------------------------------------------------------------------------------------
1249 
1250 
1251    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 3.2
1252    //---------------------------------------------------------------------------------------------
1253    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1254      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
1255 
1256        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1257          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1258      ( exp( 1.488 + qsqgev * -0.4216 + 
1259         tgev * ( -20.79 + qsqgev * 2.22 ) ) + exp( -26.52 + pow(qsqgev,1) * 9.834 + pow(qsqgev,2) * -1.006 + 
1260                                tgev * ( 16.87 + pow(qsqgev,1) * -8.514 + pow(qsqgev,2) * 0.8836 ) ) );
1261 
1262      }
1263    }
1264    //---------------------------------------------------------------------------------------------
1265 
1266 
1267    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 3.3
1268    //---------------------------------------------------------------------------------------------
1269    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1270      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1271 
1272        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1273          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1274      ( exp( 1.164 + qsqgev * -0.3716 + 
1275         tgev * ( -14.44 + qsqgev * 0.8443 ) ) + exp( -54.31 + pow(qsqgev,1) * 19.33 + pow(qsqgev,2) * -1.803 + 
1276                                  tgev * ( 57.38 + pow(qsqgev,1) * -22.72 + pow(qsqgev,2) * 2.103 ) ) );
1277 
1278      }
1279    }
1280    //---------------------------------------------------------------------------------------------
1281 
1282 
1283    //Paramaterization fo Sigma_L for Qsq: 5  to  6  at w = 3.4
1284    //---------------------------------------------------------------------------------------------
1285    if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1286      if ( wgev >= 3.4 )  { 
1287 
1288        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1289          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1290      ( exp( 1.156 + qsqgev * -0.3844 + 
1291         tgev * ( -17.55 + qsqgev * 1.314 ) ) + exp( 126.3 + pow(qsqgev,1) * -46.33 + pow(qsqgev,2) * 4.142 + 
1292                                 tgev * ( -175.3 + pow(qsqgev,1) * 61.86 + pow(qsqgev,2) * -5.561 ) ) );
1293 
1294      }
1295    }
1296    //---------------------------------------------------------------------------------------------
1297 
1298 
1299    //---------------------------------------------------------------------------------------------//
1300    //                                                                                             //
1301    //                                                                                             //
1302    //                                                                                             //
1303    //                                                                                             //
1304    //                                                                                             //
1305    //                                                                                             //
1306    //                                                                                             //
1307    //                                                                                             //
1308    //                                                                                             //
1309    //                                                                                             //
1310    //                                                                                             //
1311    //                                                                                             //
1312    //---------------------------------------------------------------------------------------------//
1313    //---------------------------------------------------------------------------------------------
1314    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2
1315    //---------------------------------------------------------------------------------------------
1316    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1317      if ( wgev <= 2.1 )  { 
1318        
1319      s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1320            pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1321        ( exp( -0.3944 + qsqgev * -0.2263 + 
1322           tgev * ( -2.257 + qsqgev * 0.1679 ) ) + exp( -118.6 + pow(qsqgev,1) * 37.12 + pow(qsqgev,2) * -2.947 + 
1323                                    tgev * ( -66.87 + pow(qsqgev,1) * 20.47 + pow(qsqgev,2) * -1.59 ) ) );
1324 
1325      }
1326    }
1327    //---------------------------------------------------------------------------------------------
1328    
1329 
1330    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.1
1331    //---------------------------------------------------------------------------------------------
1332    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1333      if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
1334 
1335        s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1336          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1337      ( exp( 0.7064 + qsqgev * -0.4166 + 
1338         tgev * ( -3.098 + qsqgev * 0.2583 ) ) + exp( -0.5318 + pow(qsqgev,1) * 0.04343 + pow(qsqgev,2) * -0.04229 + 
1339                                  tgev * ( -5.817 + pow(qsqgev,1) * 1.236 + pow(qsqgev,2) * -0.08258 ) ) );
1340 
1341      }
1342    }
1343    //---------------------------------------------------------------------------------------------
1344 
1345 
1346    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.2
1347    //---------------------------------------------------------------------------------------------
1348    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1349      if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
1350 
1351        s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1352          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1353      ( exp( 0.9339 + qsqgev * -0.4636 + 
1354         tgev * ( -3.212 + qsqgev * 0.233 ) ) + exp( 2.21 + pow(qsqgev,1) * -0.8574 + pow(qsqgev,2) * 0.02999 + 
1355                                 tgev * ( -5.793 + pow(qsqgev,1) * 1.07 + pow(qsqgev,2) * -0.06357 ) ) );
1356 
1357      }
1358    }
1359    //---------------------------------------------------------------------------------------------
1360 
1361 
1362    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.3
1363    //---------------------------------------------------------------------------------------------
1364    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1365      if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
1366 
1367        s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1368          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1369      ( exp( 1.558 + qsqgev * -0.5619 + 
1370         tgev * ( -5.179 + qsqgev * 0.4807 ) ) + exp( -0.08504 + pow(qsqgev,1) * -0.1783 + pow(qsqgev,2) * -0.02124 + 
1371                                  tgev * ( 1.176 + pow(qsqgev,1) * -1.095 + pow(qsqgev,2) * 0.09832 ) ) );
1372 
1373      }
1374    }
1375    //---------------------------------------------------------------------------------------------
1376 
1377 
1378    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.4
1379    //---------------------------------------------------------------------------------------------
1380    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1381      if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
1382 
1383        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1384          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1385      ( exp( 14.18 + qsqgev * -2.406 + 
1386         tgev * ( -23.41 + qsqgev * 3.089 ) ) + exp( -23.34 + pow(qsqgev,1) * 5.621 + pow(qsqgev,2) * -0.3728 + 
1387                                 tgev * ( 13.88 + pow(qsqgev,1) * -3.361 + pow(qsqgev,2) * 0.1534 ) ) );
1388 
1389      }
1390    }
1391    //---------------------------------------------------------------------------------------------
1392 
1393 
1394    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.5
1395    //---------------------------------------------------------------------------------------------
1396    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1397      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1398 
1399        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1400          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1401      ( exp( 10.58 + qsqgev * -1.794 + 
1402         tgev * ( -20.34 + qsqgev * 2.462 ) ) + exp( 80.53 + pow(qsqgev,1) * -26.13 + pow(qsqgev,2) * 2.042 + 
1403                                 tgev * ( -110 + pow(qsqgev,1) * 34.41 + pow(qsqgev,2) * -2.717 ) ) );
1404 
1405      }
1406    }
1407    //---------------------------------------------------------------------------------------------
1408 
1409 
1410    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.6
1411    //---------------------------------------------------------------------------------------------
1412    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1413      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1414 
1415        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1416          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1417      ( exp( 1.665 + qsqgev * -0.3821 + 
1418         tgev * ( -8.2 + qsqgev * 0.4645 ) ) + exp( 5.524 + pow(qsqgev,1) * -2.303 + pow(qsqgev,2) * 0.1493 + 
1419                                tgev * ( -6.979 + pow(qsqgev,1) * 1.548 + pow(qsqgev,2) * -0.1059 ) ) );
1420 
1421      }
1422    }
1423    //---------------------------------------------------------------------------------------------
1424 
1425 
1426    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.7
1427    //---------------------------------------------------------------------------------------------
1428    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1429      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1430 
1431        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1432          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1433      ( exp( 9.849 + qsqgev * -1.716 + 
1434         tgev * ( -26.44 + qsqgev * 3.375 ) ) + exp( 443.1 + pow(qsqgev,1) * -137.8 + pow(qsqgev,2) * 10.61 + 
1435                                 tgev * ( -347.7 + pow(qsqgev,1) * 108 + pow(qsqgev,2) * -8.408 ) ) );
1436 
1437      }
1438    }
1439    //---------------------------------------------------------------------------------------------
1440 
1441 
1442    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.8
1443    //---------------------------------------------------------------------------------------------
1444    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1445      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1446 
1447        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1448          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1449      ( exp( 2.619 + qsqgev * -0.5775 + 
1450         tgev * ( -14.33 + qsqgev * 1.442 ) ) + exp( -338.5 + pow(qsqgev,1) * 105 + pow(qsqgev,2) * -8.248 + 
1451                                 tgev * ( 211 + pow(qsqgev,1) * -66.64 + pow(qsqgev,2) * 5.243 ) ) );
1452 
1453      }
1454    }
1455    //---------------------------------------------------------------------------------------------
1456 
1457 
1458    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 2.9
1459    //---------------------------------------------------------------------------------------------
1460    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1461      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1462 
1463        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1464          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1465      ( exp( 5.105 + qsqgev * -0.9407 + 
1466         tgev * ( -41.86 + qsqgev * 5.409 ) ) + exp( 160 + pow(qsqgev,1) * -47.19 + pow(qsqgev,2) * 3.371 + 
1467                                 tgev * ( -157.9 + pow(qsqgev,1) * 45.34 + pow(qsqgev,2) * -3.271 ) ) );
1468 
1469      }
1470    }
1471    //---------------------------------------------------------------------------------------------
1472 
1473 
1474    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 3
1475    //---------------------------------------------------------------------------------------------
1476    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1477      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
1478 
1479        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1480          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1481      ( exp( 1.319 + qsqgev * -0.3737 + 
1482         tgev * ( -10.55 + qsqgev * 0.6363 ) ) + exp( 1.184 + pow(qsqgev,1) * -1.116 + pow(qsqgev,2) * 0.05812 + 
1483                                  tgev * ( -6.962 + pow(qsqgev,1) * 1.295 + pow(qsqgev,2) * -0.08359 ) ) );
1484 
1485      }
1486    }
1487    //---------------------------------------------------------------------------------------------
1488 
1489 
1490    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 3.1
1491    //---------------------------------------------------------------------------------------------
1492    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1493      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
1494 
1495        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1496          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1497      ( exp( 1.307 + qsqgev * -0.3807 + 
1498         tgev * ( -13.12 + qsqgev * 0.9742 ) ) + exp( -11.66 + pow(qsqgev,1) * 3.171 + pow(qsqgev,2) * -0.2985 + 
1499                                  tgev * ( 4.021 + pow(qsqgev,1) * -2.5 + pow(qsqgev,2) * 0.2348 ) ) );
1500 
1501      }
1502    }
1503    //---------------------------------------------------------------------------------------------
1504 
1505 
1506    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 3.2
1507    //---------------------------------------------------------------------------------------------
1508    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1509      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
1510 
1511        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1512          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1513      ( exp( 1.126 + qsqgev * -0.3601 + 
1514         tgev * ( -12.68 + qsqgev * 0.8358 ) ) + exp( 1.442 + pow(qsqgev,1) * -1.093 + pow(qsqgev,2) * 0.04699 + 
1515                                  tgev * ( -7.147 + pow(qsqgev,1) * 1.054 + pow(qsqgev,2) * -0.05304 ) ) );
1516 
1517      }
1518    }
1519    //---------------------------------------------------------------------------------------------
1520 
1521 
1522    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 3.3
1523    //---------------------------------------------------------------------------------------------
1524    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1525      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1526 
1527        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1528          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1529      ( exp( 1.071 + qsqgev * -0.3599 + 
1530         tgev * ( -13.37 + qsqgev * 0.8772 ) ) + exp( -16.13 + pow(qsqgev,1) * 4.267 + pow(qsqgev,2) * -0.3616 + 
1531                                  tgev * ( 10.33 + pow(qsqgev,1) * -4.322 + pow(qsqgev,2) * 0.3544 ) ) );
1532 
1533      }
1534    }
1535    //---------------------------------------------------------------------------------------------
1536 
1537 
1538    //Paramaterization fo Sigma_L for Qsq: 6  to  7  at w = 3.4
1539    //---------------------------------------------------------------------------------------------
1540    if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1541      if ( wgev >= 3.4 )  { 
1542 
1543        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1544          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1545      ( exp( 0.514 + qsqgev * -0.2839 + 
1546         tgev * ( -7.008 + qsqgev * -0.1449 ) ) + exp( -146.2 + pow(qsqgev,1) * 43.26 + pow(qsqgev,2) * -3.279 + 
1547                                   tgev * ( 132.8 + pow(qsqgev,1) * -41.07 + pow(qsqgev,2) * 3.104 ) ) );
1548 
1549      }
1550    }
1551    //---------------------------------------------------------------------------------------------
1552 
1553 
1554 
1555    //---------------------------------------------------------------------------------------------//
1556    //                                                                                             //
1557    //                                                                                             //
1558    //                                                                                             //
1559    //                                                                                             //
1560    //                                                                                             //
1561    //                                                                                             //
1562    //                                                                                             //
1563    //                                                                                             //
1564    //                                                                                             //
1565    //                                                                                             //
1566    //                                                                                             //
1567    //                                                                                             //
1568    //---------------------------------------------------------------------------------------------//
1569    //---------------------------------------------------------------------------------------------
1570    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 2.4
1571    //---------------------------------------------------------------------------------------------
1572    if ( qsqgev >= 7 )  {
1573      if ( wgev < 2.5 )  { 
1574 
1575        s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1576          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1577      ( exp( 0.4744 + qsqgev * -0.4071 + 
1578         tgev * ( -4.841 + qsqgev * 0.3756 ) ) + exp( -2.159 + pow(qsqgev,1) * 0.3347 + pow(qsqgev,2) * -0.05295 + 
1579                                  tgev * ( -1.705 + pow(qsqgev,1) * -0.165 + pow(qsqgev,2) * 0.01938 ) ) );
1580 
1581      }
1582    }
1583    //---------------------------------------------------------------------------------------------
1584 
1585 
1586    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 2.5
1587    //---------------------------------------------------------------------------------------------
1588    if ( qsqgev >= 7  )  { 
1589      if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1590 
1591        s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1592          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1593      ( exp( 1.016 + qsqgev * -0.4725 + 
1594         tgev * ( -7.169 + qsqgev * 0.6343 ) ) + exp( 2.248 + pow(qsqgev,1) * -0.8677 + pow(qsqgev,2) * 0.0279 + 
1595                                  tgev * ( 1.721 + pow(qsqgev,1) * -1.04 + pow(qsqgev,2) * 0.0706 ) ) );
1596 
1597      }
1598    }
1599    //---------------------------------------------------------------------------------------------
1600 
1601 
1602    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 2.6
1603    //---------------------------------------------------------------------------------------------
1604    if ( qsqgev >= 7  )  { 
1605      if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1606 
1607        s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1608          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1609      ( exp( 3.906 + qsqgev * -0.8198 + 
1610         tgev * ( -10.5 + qsqgev * 0.9657 ) ) + exp( -86.64 + pow(qsqgev,1) * 22.59 + pow(qsqgev,2) * -1.52 + 
1611                                 tgev * ( 105.4 + pow(qsqgev,1) * -28.66 + pow(qsqgev,2) * 1.908 ) ) );
1612 
1613      }
1614    }
1615    //---------------------------------------------------------------------------------------------
1616 
1617 
1618    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 2.7
1619    //---------------------------------------------------------------------------------------------
1620    if ( qsqgev >= 7  )  { 
1621      if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1622 
1623        s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1624          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1625      ( exp( -3.804 + qsqgev * 0.2026 + 
1626         tgev * ( 5.106 + qsqgev * -1.111 ) ) + exp( -123.3 + pow(qsqgev,1) * 32.75 + pow(qsqgev,2) * -2.222 + 
1627                                 tgev * ( 243.7 + pow(qsqgev,1) * -66.94 + pow(qsqgev,2) * 4.537 ) ) );
1628 
1629      }
1630    }
1631    //---------------------------------------------------------------------------------------------
1632 
1633 
1634    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 2.8
1635    //---------------------------------------------------------------------------------------------
1636    if ( qsqgev >= 7  )  { 
1637      if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1638 
1639        s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1640          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1641      ( exp( 5.326 + qsqgev * -0.9887 + 
1642         tgev * ( -8.259 + qsqgev * 0.6319 ) ) + exp( -624.5 + pow(qsqgev,1) * 161.8 + pow(qsqgev,2) * -10.54 + 
1643                                  tgev * ( 662.6 + pow(qsqgev,1) * -173.7 + pow(qsqgev,2) * 11.33 ) ) );
1644 
1645      }
1646    }
1647    //---------------------------------------------------------------------------------------------
1648 
1649 
1650    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 2.9
1651    //---------------------------------------------------------------------------------------------
1652    if ( qsqgev >= 7  )  { 
1653      if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1654 
1655        s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1656          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1657      ( exp( 4.206 + qsqgev * -0.7693 + 
1658         tgev * ( -22.87 + qsqgev * 2.401 ) ) + exp( 24.26 + pow(qsqgev,1) * -5.183 + pow(qsqgev,2) * 0.1576 + 
1659                                 tgev * ( 19.94 + pow(qsqgev,1) * -7.34 + pow(qsqgev,2) * 0.6158 ) ) );
1660 
1661      }
1662    }
1663    //---------------------------------------------------------------------------------------------
1664 
1665 
1666    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 3
1667    //---------------------------------------------------------------------------------------------
1668    if ( qsqgev >= 7  )  { 
1669      if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
1670 
1671        s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1672          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1673      ( exp( 0.3563 + qsqgev * -0.2464 + 
1674         tgev * ( -5.658 + qsqgev * 0.002434 ) ) + exp( 23.49 + pow(qsqgev,1) * -8.036 + pow(qsqgev,2) * 0.5752 + 
1675                                    tgev * ( -42.58 + pow(qsqgev,1) * 11.5 + pow(qsqgev,2) * -0.8028 ) ) );
1676 
1677      }
1678    }
1679    //---------------------------------------------------------------------------------------------
1680 
1681 
1682    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 3.1
1683    //---------------------------------------------------------------------------------------------
1684    if ( qsqgev >= 7  )  { 
1685      if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
1686 
1687        s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1688          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1689      ( exp( 1.076 + qsqgev * -0.3468 + 
1690         tgev * ( -10.44 + qsqgev * 0.574 ) ) + exp( -41.96 + pow(qsqgev,1) * 10.53 + pow(qsqgev,2) * -0.7292 + 
1691                                 tgev * ( 33.06 + pow(qsqgev,1) * -9.616 + pow(qsqgev,2) * 0.6573 ) ) );
1692 
1693      }
1694    }
1695    //---------------------------------------------------------------------------------------------
1696 
1697 
1698    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 3.2
1699    //---------------------------------------------------------------------------------------------
1700    if ( qsqgev >= 7  )  { 
1701      if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
1702 
1703        s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1704          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1705      ( exp( 0.9067 + qsqgev * -0.3293 + 
1706         tgev * ( -10.8 + qsqgev * 0.5637 ) ) + exp( -14.35 + pow(qsqgev,1) * 3.071 + pow(qsqgev,2) * -0.2259 + 
1707                                 tgev * ( 12.25 + pow(qsqgev,1) * -4.053 + pow(qsqgev,2) * 0.2812 ) ) );
1708 
1709      }
1710    }
1711    //---------------------------------------------------------------------------------------------
1712 
1713 
1714    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 3.3
1715    //---------------------------------------------------------------------------------------------
1716    if ( qsqgev >= 7  )  { 
1717      if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1718 
1719        s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1720          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1721      ( exp( 0.7918 + qsqgev * -0.3213 + 
1722         tgev * ( -11.24 + qsqgev * 0.5752 ) ) + exp( 2.169 + pow(qsqgev,1) * -1.341 + pow(qsqgev,2) * 0.06773 + 
1723                                  tgev * ( -9.616 + pow(qsqgev,1) * 1.704 + pow(qsqgev,2) * -0.101 ) ) );
1724 
1725      }
1726    }
1727    //---------------------------------------------------------------------------------------------
1728 
1729 
1730    //Paramaterization fo Sigma_L for Qsq: 7  to  8  at w = 3.4
1731    //---------------------------------------------------------------------------------------------
1732    if ( qsqgev >= 7  )  { 
1733      if ( wgev >= 3.4 )  { 
1734 
1735        s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1736          pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
1737      ( exp( 3.385 + qsqgev * -0.6811 + 
1738         tgev * ( -18.29 + qsqgev * 1.481 ) ) + exp( 89.47 + pow(qsqgev,1) * -25.23 + pow(qsqgev,2) * 1.699 + 
1739                                 tgev * ( -217.8 + pow(qsqgev,1) * 58.41 + pow(qsqgev,2) * -3.961 ) ) );
1740 
1741      }
1742    }
1743    //---------------------------------------------------------------------------------------------
1744 
1745 
1746  return s;
1747     
1748 
1749 }