Back to home page

EIC code displayed by LXR
 
 

    

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

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 /* double fSigmaT; */
0019 
0020 double MySigmaT( double qsqgev, double tgev, double wgev ) { 
0021 
0022 
0023   double s =  0.0;
0024   /* double wgev   = par[0]; */
0025   /* double qsqgev = par[1]; */
0026   /* double tgev   = x[0]; */
0027 
0028    //---------------------------------------------------------------------------------------------//
0029    //                                                                                             //
0030    //                                                                                             //
0031    //                                                                                             //
0032    //                                                                                             //
0033    //                                                                                             //
0034    //                                                                                             //
0035    //                                                                                             //
0036    //                                                                                             //
0037    //                                                                                             //
0038    //                                                                                             //
0039    //                                                                                             //
0040    //                                                                                             //
0041    //---------------------------------------------------------------------------------------------//
0042   //---------------------------------------------------------------------------------------------
0043   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2
0044   //---------------------------------------------------------------------------------------------
0045   if ( qsqgev <= 2  )  { 
0046     if ( wgev <= 2.1 )  { 
0047 
0048       s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0049         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0050     ( exp( (  pow( qsqgev , 0 ) * 2.546 + pow( qsqgev , 1 ) * -1.074 + pow( qsqgev , 2 ) * 0 ) + 
0051            tgev * ( pow( qsqgev , 0 ) * -45.04  + pow( qsqgev , 1 ) * 60.66 + pow( qsqgev , 2 ) * -22.69 ) ) + 
0052       (  pow( qsqgev , 0 ) * 6.504 + pow( qsqgev , 1 ) * -7.199 + pow( qsqgev , 2 ) * 2.295 ) + 
0053       ( tgev * ( pow( qsqgev , 0 ) * -4.284  + pow( qsqgev , 1 ) * 5.635 + pow( qsqgev , 2 ) * -1.911 ) ) );
0054 
0055     }
0056   }
0057   //---------------------------------------------------------------------------------------------
0058 
0059 
0060   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.1
0061   //---------------------------------------------------------------------------------------------
0062   if ( qsqgev <= 2  )  { 
0063     if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0064 
0065       s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0066         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0067     ( exp( (  pow( qsqgev , 0 ) * 2.332 + pow( qsqgev , 1 ) * -1.054 + pow( qsqgev , 2 ) * 0 ) + 
0068            tgev * ( pow( qsqgev , 0 ) * -40.39  + pow( qsqgev , 1 ) * 51.49 + pow( qsqgev , 2 ) * -18.97 ) ) + 
0069       (  pow( qsqgev , 0 ) * 4.724 + pow( qsqgev , 1 ) * -5.13 + pow( qsqgev , 2 ) * 1.628 ) + 
0070       ( tgev * ( pow( qsqgev , 0 ) * -3.525  + pow( qsqgev , 1 ) * 4.489 + pow( qsqgev , 2 ) * -1.495 ) ) );
0071 
0072     }
0073   }
0074   //---------------------------------------------------------------------------------------------
0075 
0076 
0077   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.2
0078   //---------------------------------------------------------------------------------------------
0079   if ( qsqgev <= 2  )  { 
0080     if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0081 
0082       s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0083         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0084     ( exp( (  pow( qsqgev , 0 ) * 2.12 + pow( qsqgev , 1 ) * -1.023 + pow( qsqgev , 2 ) * 0 ) + 
0085            tgev * ( pow( qsqgev , 0 ) * -34.61  + pow( qsqgev , 1 ) * 40.38 + pow( qsqgev , 2 ) * -14.55 ) ) + 
0086       (  pow( qsqgev , 0 ) * 3.461 + pow( qsqgev , 1 ) * -3.685 + pow( qsqgev , 2 ) * 1.166 ) + 
0087       ( tgev * ( pow( qsqgev , 0 ) * -2.799  + pow( qsqgev , 1 ) * 3.458 + pow( qsqgev , 2 ) * -1.135 ) ) );
0088 
0089     }
0090   }
0091   //---------------------------------------------------------------------------------------------
0092 
0093 
0094   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.3
0095   //---------------------------------------------------------------------------------------------
0096   if ( qsqgev <= 2  )  { 
0097     if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0098 
0099       s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0100         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0101     ( exp( (  pow( qsqgev , 0 ) * 1.91 + pow( qsqgev , 1 ) * -0.9857 + pow( qsqgev , 2 ) * 0 ) + 
0102            tgev * ( pow( qsqgev , 0 ) * -27.26  + pow( qsqgev , 1 ) * 26.45 + pow( qsqgev , 2 ) * -9.059 ) ) + 
0103       (  pow( qsqgev , 0 ) * 2.538 + pow( qsqgev , 1 ) * -2.63 + pow( qsqgev , 2 ) * 0.8264 ) + 
0104       ( tgev * ( pow( qsqgev , 0 ) * -2.161  + pow( qsqgev , 1 ) * 2.574 + pow( qsqgev , 2 ) * -0.831 ) ) );
0105 
0106     }
0107   }
0108   //---------------------------------------------------------------------------------------------
0109 
0110 
0111   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.4
0112   //---------------------------------------------------------------------------------------------
0113   if ( qsqgev <= 2  )  { 
0114     if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0115 
0116       s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0117         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0118     ( exp( (  pow( qsqgev , 0 ) * 1.694 + pow( qsqgev , 1 ) * -0.9393 + pow( qsqgev , 2 ) * 0 ) + 
0119            tgev * ( pow( qsqgev , 0 ) * -18.64  + pow( qsqgev , 1 ) * 10.18 + pow( qsqgev , 2 ) * -2.648 ) ) + 
0120       (  pow( qsqgev , 0 ) * 1.771 + pow( qsqgev , 1 ) * -1.691 + pow( qsqgev , 2 ) * 0.5114 ) + 
0121       ( tgev * ( pow( qsqgev , 0 ) * -1.533  + pow( qsqgev , 1 ) * 1.671 + pow( qsqgev , 2 ) * -0.5125 ) ) );
0122 
0123     }
0124   }
0125   //---------------------------------------------------------------------------------------------
0126 
0127 
0128   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.5
0129   //---------------------------------------------------------------------------------------------
0130   if ( qsqgev <= 2  )  { 
0131     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0132 
0133       s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0134         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0135     ( exp( (  pow( qsqgev , 0 ) * 1.521 + pow( qsqgev , 1 ) * -0.9198 + pow( qsqgev , 2 ) * 0 ) + 
0136            tgev * ( pow( qsqgev , 0 ) * -17.9  + pow( qsqgev , 1 ) * 8.338 + pow( qsqgev , 2 ) * -1.914 ) ) + 
0137       (  pow( qsqgev , 0 ) * 1.326 + pow( qsqgev , 1 ) * -1.224 + pow( qsqgev , 2 ) * 0.3666 ) + 
0138       ( tgev * ( pow( qsqgev , 0 ) * -1.184  + pow( qsqgev , 1 ) * 1.236 + pow( qsqgev , 2 ) * -0.3718 ) ) );
0139 
0140     }
0141   }
0142   //---------------------------------------------------------------------------------------------
0143 
0144 
0145   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.6
0146   //---------------------------------------------------------------------------------------------
0147   if ( qsqgev <= 2  )  { 
0148     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0149 
0150       s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0151         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0152     ( exp( (  pow( qsqgev , 0 ) * 1.358 + pow( qsqgev , 1 ) * -0.9031 + pow( qsqgev , 2 ) * 0 ) + 
0153            tgev * ( pow( qsqgev , 0 ) * -17.73  + pow( qsqgev , 1 ) * 7.539 + pow( qsqgev , 2 ) * -1.597 ) ) + 
0154       (  pow( qsqgev , 0 ) * 1.016 + pow( qsqgev , 1 ) * -0.9168 + pow( qsqgev , 2 ) * 0.2741 ) + 
0155       ( tgev * ( pow( qsqgev , 0 ) * -0.9309  + pow( qsqgev , 1 ) * 0.9422 + pow( qsqgev , 2 ) * -0.28 ) ) );
0156 
0157     }
0158   }
0159   //---------------------------------------------------------------------------------------------
0160 
0161 
0162   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.7
0163   //---------------------------------------------------------------------------------------------
0164   if ( qsqgev <= 2  )  { 
0165     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0166 
0167       s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0168         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0169     ( exp( (  pow( qsqgev , 0 ) * 1.202 + pow( qsqgev , 1 ) * -0.8882 + pow( qsqgev , 2 ) * 0 ) + 
0170            tgev * ( pow( qsqgev , 0 ) * -17.69  + pow( qsqgev , 1 ) * 6.986 + pow( qsqgev , 2 ) * -1.374 ) ) + 
0171       (  pow( qsqgev , 0 ) * 0.7906 + pow( qsqgev , 1 ) * -0.7023 + pow( qsqgev , 2 ) * 0.2104 ) + 
0172       ( tgev * ( pow( qsqgev , 0 ) * -0.7403  + pow( qsqgev , 1 ) * 0.7301 + pow( qsqgev , 2 ) * -0.2154 ) ) );
0173 
0174     }
0175   }
0176   //---------------------------------------------------------------------------------------------
0177 
0178 
0179   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.8
0180   //---------------------------------------------------------------------------------------------
0181   if ( qsqgev <= 2  )  { 
0182     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0183 
0184       s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0185         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0186     ( exp( (  pow( qsqgev , 0 ) * 1.052 + pow( qsqgev , 1 ) * -0.8745 + pow( qsqgev , 2 ) * 0 ) + 
0187            tgev * ( pow( qsqgev , 0 ) * -17.73  + pow( qsqgev , 1 ) * 6.568 + pow( qsqgev , 2 ) * -1.2 ) ) + 
0188       (  pow( qsqgev , 0 ) * 0.6239 + pow( qsqgev , 1 ) * -0.5481 + pow( qsqgev , 2 ) * 0.165 ) + 
0189       ( tgev * ( pow( qsqgev , 0 ) * -0.5943  + pow( qsqgev , 1 ) * 0.5737 + pow( qsqgev , 2 ) * -0.1686 ) ) );
0190 
0191     }
0192   }
0193   //---------------------------------------------------------------------------------------------
0194 
0195 
0196   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 2.9
0197   //---------------------------------------------------------------------------------------------
0198   if ( qsqgev <= 2  )  { 
0199     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0200 
0201       s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0202         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0203     ( exp( (  pow( qsqgev , 0 ) * 0.9084 + pow( qsqgev , 1 ) * -0.8628 + pow( qsqgev , 2 ) * 0 ) + 
0204            tgev * ( pow( qsqgev , 0 ) * -17.94  + pow( qsqgev , 1 ) * 6.43 + pow( qsqgev , 2 ) * -1.13 ) ) + 
0205       (  pow( qsqgev , 0 ) * 0.4999 + pow( qsqgev , 1 ) * -0.4371 + pow( qsqgev , 2 ) * 0.1326 ) + 
0206       ( tgev * ( pow( qsqgev , 0 ) * -0.483  + pow( qsqgev , 1 ) * 0.4591 + pow( qsqgev , 2 ) * -0.1349 ) ) );
0207 
0208     }
0209   }
0210   //---------------------------------------------------------------------------------------------
0211 
0212 
0213   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 3
0214   //---------------------------------------------------------------------------------------------
0215   if ( qsqgev <= 2  )  { 
0216     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0217 
0218       s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0219         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0220     ( exp( (  pow( qsqgev , 0 ) * 0.7705 + pow( qsqgev , 1 ) * -0.852 + pow( qsqgev , 2 ) * 0 ) + 
0221            tgev * ( pow( qsqgev , 0 ) * -18.28  + pow( qsqgev , 1 ) * 6.522 + pow( qsqgev , 2 ) * -1.153 ) ) + 
0222       (  pow( qsqgev , 0 ) * 0.4062 + pow( qsqgev , 1 ) * -0.3561 + pow( qsqgev , 2 ) * 0.1092 ) + 
0223       ( tgev * ( pow( qsqgev , 0 ) * -0.3972  + pow( qsqgev , 1 ) * 0.3742 + pow( qsqgev , 2 ) * -0.1106 ) ) );
0224 
0225     }
0226   }
0227   //---------------------------------------------------------------------------------------------
0228 
0229 
0230   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 3.1
0231   //---------------------------------------------------------------------------------------------
0232   if ( qsqgev <= 2  )  { 
0233     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0234 
0235       s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0236         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0237     ( exp( (  pow( qsqgev , 0 ) * 0.6384 + pow( qsqgev , 1 ) * -0.8428 + pow( qsqgev , 2 ) * 0 ) + 
0238            tgev * ( pow( qsqgev , 0 ) * -18.69  + pow( qsqgev , 1 ) * 6.698 + pow( qsqgev , 2 ) * -1.192 ) ) + 
0239       (  pow( qsqgev , 0 ) * 0.3342 + pow( qsqgev , 1 ) * -0.2948 + pow( qsqgev , 2 ) * 0.09148 ) + 
0240       ( tgev * ( pow( qsqgev , 0 ) * -0.33  + pow( qsqgev , 1 ) * 0.3095 + pow( qsqgev , 2 ) * -0.09206 ) ) );
0241 
0242     }
0243   }
0244   //---------------------------------------------------------------------------------------------
0245 
0246 
0247   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 3.2
0248   //---------------------------------------------------------------------------------------------
0249   if ( qsqgev <= 2  )  { 
0250     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0251 
0252       s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0253         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0254     ( exp( (  pow( qsqgev , 0 ) * 0.5101 + pow( qsqgev , 1 ) * -0.8344 + pow( qsqgev , 2 ) * 0 ) + 
0255            tgev * ( pow( qsqgev , 0 ) * -19.11  + pow( qsqgev , 1 ) * 6.918 + pow( qsqgev , 2 ) * -1.25 ) ) + 
0256       (  pow( qsqgev , 0 ) * 0.2778 + pow( qsqgev , 1 ) * -0.2475 + pow( qsqgev , 2 ) * 0.07769 ) + 
0257       ( tgev * ( pow( qsqgev , 0 ) * -0.2765  + pow( qsqgev , 1 ) * 0.259 + pow( qsqgev , 2 ) * -0.07767 ) ) );
0258 
0259     }
0260   }
0261   //---------------------------------------------------------------------------------------------
0262 
0263 
0264   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 3.3
0265   //---------------------------------------------------------------------------------------------
0266   if ( qsqgev <= 2  )  { 
0267     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
0268 
0269       s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0270         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0271     ( exp( (  pow( qsqgev , 0 ) * 0.3852 + pow( qsqgev , 1 ) * -0.8252 + pow( qsqgev , 2 ) * 0 ) + 
0272            tgev * ( pow( qsqgev , 0 ) * -19.58  + pow( qsqgev , 1 ) * 7.256 + pow( qsqgev , 2 ) * -1.365 ) ) + 
0273       (  pow( qsqgev , 0 ) * 0.2334 + pow( qsqgev , 1 ) * -0.2109 + pow( qsqgev , 2 ) * 0.06704 ) + 
0274       ( tgev * ( pow( qsqgev , 0 ) * -0.2338  + pow( qsqgev , 1 ) * 0.2199 + pow( qsqgev , 2 ) * -0.06668 ) ) );
0275 
0276     }
0277   }
0278   //---------------------------------------------------------------------------------------------
0279 
0280 
0281   //Paramaterization fo Sigma_T for Qsq: 0.1  to  2  at w = 3.4
0282   //---------------------------------------------------------------------------------------------
0283   if ( qsqgev <= 2  )  { 
0284     if ( wgev >= 3.4 )  { 
0285 
0286       s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0287         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) * 
0288     ( exp( (  pow( qsqgev , 0 ) * 0.2654 + pow( qsqgev , 1 ) * -0.8181 + pow( qsqgev , 2 ) * 0 ) + 
0289            tgev * ( pow( qsqgev , 0 ) * -20.09  + pow( qsqgev , 1 ) * 7.614 + pow( qsqgev , 2 ) * -1.474 ) ) + 
0290       (  pow( qsqgev , 0 ) * 0.1979 + pow( qsqgev , 1 ) * -0.1816 + pow( qsqgev , 2 ) * 0.05842 ) + 
0291       ( tgev * ( pow( qsqgev , 0 ) * -0.1993  + pow( qsqgev , 1 ) * 0.1887 + pow( qsqgev , 2 ) * -0.05786 ) ) );
0292 
0293     }
0294   }
0295   //---------------------------------------------------------------------------------------------
0296 
0297 
0298    //---------------------------------------------------------------------------------------------//
0299    //                                                                                             //
0300    //                                                                                             //
0301    //                                                                                             //
0302    //                                                                                             //
0303    //                                                                                             //
0304    //                                                                                             //
0305    //                                                                                             //
0306    //                                                                                             //
0307    //                                                                                             //
0308    //                                                                                             //
0309    //                                                                                             //
0310    //                                                                                             //
0311    //---------------------------------------------------------------------------------------------//
0312   //---------------------------------------------------------------------------------------------
0313   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2
0314   //---------------------------------------------------------------------------------------------
0315   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0316     if ( wgev <= 2.1 )  { 
0317 
0318       s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0319         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0320         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 3.089 + pow( qsqgev , 1 ) * -0.8796 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0321           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -6.194 + pow( qsqgev , 1 ) * 1.932 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0322           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 6.944 + pow( qsqgev , 1 ) * -2.2 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0323           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -2.724 + pow( qsqgev , 1 ) * 0.866 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0324           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0325           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0326           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0327 
0328     }
0329   }
0330   //---------------------------------------------------------------------------------------------
0331 
0332 
0333   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.1
0334   //---------------------------------------------------------------------------------------------
0335   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0336     if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0337 
0338       s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0339         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0340         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 2.612 + pow( qsqgev , 1 ) * -0.7368 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0341           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -5.789 + pow( qsqgev , 1 ) * 1.774 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0342           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 6.629 + pow( qsqgev , 1 ) * -2.064 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0343           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -2.635 + pow( qsqgev , 1 ) * 0.8255 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0344           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0345           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0346           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0347 
0348     }
0349   }
0350   //---------------------------------------------------------------------------------------------
0351 
0352 
0353   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.2
0354   //---------------------------------------------------------------------------------------------
0355   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0356     if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0357 
0358       s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0359         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0360         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 2.219 + pow( qsqgev , 1 ) * -0.6202 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0361           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -5.349 + pow( qsqgev , 1 ) * 1.616 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0362           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 6.25 + pow( qsqgev , 1 ) * -1.918 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0363           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -2.516 + pow( qsqgev , 1 ) * 0.778 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0364           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0365           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0366           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0367 
0368     }
0369   }
0370   //---------------------------------------------------------------------------------------------
0371 
0372 
0373   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.3
0374   //---------------------------------------------------------------------------------------------
0375   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0376     if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0377 
0378       s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0379         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0380         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.893 + pow( qsqgev , 1 ) * -0.5244 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0381           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -4.903 + pow( qsqgev , 1 ) * 1.463 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0382           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 5.836 + pow( qsqgev , 1 ) * -1.768 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0383           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -2.377 + pow( qsqgev , 1 ) * 0.7259 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0384           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0385           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0386           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0387 
0388     }
0389   }
0390   //---------------------------------------------------------------------------------------------
0391 
0392 
0393   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.4
0394   //---------------------------------------------------------------------------------------------
0395   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0396     if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0397 
0398       s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0399         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0400         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.621 + pow( qsqgev , 1 ) * -0.4451 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0401           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -4.468 + pow( qsqgev , 1 ) * 1.317 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0402           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 5.409 + pow( qsqgev , 1 ) * -1.62 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0403           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -2.224 + pow( qsqgev , 1 ) * 0.6717 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0404           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0405           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0406           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0407 
0408     }
0409   }
0410   //---------------------------------------------------------------------------------------------
0411 
0412 
0413   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.5
0414   //---------------------------------------------------------------------------------------------
0415   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0416     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0417 
0418       s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0419         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0420         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.394 + pow( qsqgev , 1 ) * -0.3794 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0421           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -4.055 + pow( qsqgev , 1 ) * 1.183 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0422           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 4.984 + pow( qsqgev , 1 ) * -1.476 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0423           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -2.067 + pow( qsqgev , 1 ) * 0.6174 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0424           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0425           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0426           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0427 
0428     }
0429   }
0430   //---------------------------------------------------------------------------------------------
0431 
0432 
0433   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.6
0434   //---------------------------------------------------------------------------------------------
0435   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0436     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0437 
0438       s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0439         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0440         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.202 + pow( qsqgev , 1 ) * -0.3246 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0441           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -3.67 + pow( qsqgev , 1 ) * 1.059 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0442           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 4.573 + pow( qsqgev , 1 ) * -1.341 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0443           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.91 + pow( qsqgev , 1 ) * 0.5647 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0444           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0445           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0446           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0447 
0448     }
0449   }
0450   //---------------------------------------------------------------------------------------------
0451 
0452 
0453   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.7
0454   //---------------------------------------------------------------------------------------------
0455   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0456     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0457 
0458       s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0459         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0460         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.041 + pow( qsqgev , 1 ) * -0.2787 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0461           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -3.315 + pow( qsqgev , 1 ) * 0.9472 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0462           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 4.182 + pow( qsqgev , 1 ) * -1.214 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0463           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.758 + pow( qsqgev , 1 ) * 0.5144 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0464           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0465           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0466           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0467 
0468     }
0469   }
0470   //---------------------------------------------------------------------------------------------
0471 
0472 
0473   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.8
0474   //---------------------------------------------------------------------------------------------
0475   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0476     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0477 
0478       s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0479         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0480         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.9044 + pow( qsqgev , 1 ) * -0.2402 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0481           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.991 + pow( qsqgev , 1 ) * 0.8465 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0482           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 3.817 + pow( qsqgev , 1 ) * -1.097 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0483           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.613 + pow( qsqgev , 1 ) * 0.4674 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0484           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0485           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0486           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0487 
0488     }
0489   }
0490   //---------------------------------------------------------------------------------------------
0491 
0492 
0493   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 2.9
0494   //---------------------------------------------------------------------------------------------
0495   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0496     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0497 
0498       s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0499         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0500         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.7885 + pow( qsqgev , 1 ) * -0.2078 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0501           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.697 + pow( qsqgev , 1 ) * 0.7564 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0502           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 3.478 + pow( qsqgev , 1 ) * -0.9906 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0503           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.477 + pow( qsqgev , 1 ) * 0.4239 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0504           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0505           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0506           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0507 
0508     }
0509   }
0510   //---------------------------------------------------------------------------------------------
0511 
0512 
0513   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 3
0514   //---------------------------------------------------------------------------------------------
0515   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0516     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0517 
0518       s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0519         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0520         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.6898 + pow( qsqgev , 1 ) * -0.1805 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0521           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.433 + pow( qsqgev , 1 ) * 0.6762 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0522           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 3.167 + pow( qsqgev , 1 ) * -0.8939 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0523           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.351 + pow( qsqgev , 1 ) * 0.3842 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0524           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0525           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0526           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0527 
0528     }
0529   }
0530   //---------------------------------------------------------------------------------------------
0531 
0532 
0533   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 3.1
0534   //---------------------------------------------------------------------------------------------
0535   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0536     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0537 
0538       s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0539         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0540         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.6054 + pow( qsqgev , 1 ) * -0.1573 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0541           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.195 + pow( qsqgev , 1 ) * 0.6049 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0542           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.883 + pow( qsqgev , 1 ) * -0.8068 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0543           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.235 + pow( qsqgev , 1 ) * 0.348 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0544           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0545           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0546           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0547 
0548     }
0549   }
0550   //---------------------------------------------------------------------------------------------
0551 
0552 
0553   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 3.2
0554   //---------------------------------------------------------------------------------------------
0555   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0556     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0557 
0558       s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0559         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0560         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.533 + pow( qsqgev , 1 ) * -0.1375 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0561           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.982 + pow( qsqgev , 1 ) * 0.5417 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0562           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.625 + pow( qsqgev , 1 ) * -0.7284 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0563           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.128 + pow( qsqgev , 1 ) * 0.3153 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0564           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0565           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0566           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0567 
0568     }
0569   }
0570   //---------------------------------------------------------------------------------------------
0571 
0572 
0573   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 3.3
0574   //---------------------------------------------------------------------------------------------
0575   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0576     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
0577 
0578       s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0579         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0580         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.4708 + pow( qsqgev , 1 ) * -0.1207 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0581           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.792 + pow( qsqgev , 1 ) * 0.4857 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0582           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.391 + pow( qsqgev , 1 ) * -0.6582 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0583           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -1.031 + pow( qsqgev , 1 ) * 0.2858 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0584           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0585           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0586           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0587 
0588     }
0589   }
0590   //---------------------------------------------------------------------------------------------
0591 
0592 
0593   //Paramaterization fo Sigma_T for Qsq: 2  to  3  at w = 3.4
0594   //---------------------------------------------------------------------------------------------
0595   if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) )  { 
0596     if ( wgev >= 3.4 )  { 
0597 
0598       s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0599         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0600         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.4171 + pow( qsqgev , 1 ) * -0.1063 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0601           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.621 + pow( qsqgev , 1 ) * 0.4362 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0602           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.179 + pow( qsqgev , 1 ) * -0.5953 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0603           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.9425 + pow( qsqgev , 1 ) * 0.2592 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0604           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0605           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0606           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0607 
0608     }
0609   }
0610   //---------------------------------------------------------------------------------------------
0611 
0612    //---------------------------------------------------------------------------------------------//
0613    //                                                                                             //
0614    //                                                                                             //
0615    //                                                                                             //
0616    //                                                                                             //
0617    //                                                                                             //
0618    //                                                                                             //
0619    //                                                                                             //
0620    //                                                                                             //
0621    //                                                                                             //
0622    //                                                                                             //
0623    //                                                                                             //
0624    //                                                                                             //
0625    //---------------------------------------------------------------------------------------------//
0626   //---------------------------------------------------------------------------------------------
0627   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2
0628   //---------------------------------------------------------------------------------------------
0629   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0630     if ( wgev <= 2.1 )  { 
0631 
0632       s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0633         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0634         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.45 + pow( qsqgev , 1 ) * -0.3202 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0635           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.038 + pow( qsqgev , 1 ) * 0.5048 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0636           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.034 + pow( qsqgev , 1 ) * -0.5094 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0637           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.7555 + pow( qsqgev , 1 ) * 0.1876 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0638           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0639           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0640           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0641 
0642     }
0643   }
0644   //---------------------------------------------------------------------------------------------
0645 
0646 
0647   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.1
0648   //---------------------------------------------------------------------------------------------
0649   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0650     if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0651 
0652       s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0653         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0654         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.259 + pow( qsqgev , 1 ) * -0.2752 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0655           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.05 + pow( qsqgev , 1 ) * 0.4916 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0656           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.118 + pow( qsqgev , 1 ) * -0.5127 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0657           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.7953 + pow( qsqgev , 1 ) * 0.1921 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0658           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0659           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0660           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0661 
0662     }
0663   }
0664   //---------------------------------------------------------------------------------------------
0665 
0666 
0667   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.2
0668   //---------------------------------------------------------------------------------------------
0669   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0670     if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0671 
0672       s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0673         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0674         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 1.097 + pow( qsqgev , 1 ) * -0.2378 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0675           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -2.01 + pow( qsqgev , 1 ) * 0.472 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0676           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.142 + pow( qsqgev , 1 ) * -0.5072 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0677           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.8164 + pow( qsqgev , 1 ) * 0.1936 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0678           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0679           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0680           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0681 
0682     }
0683   }
0684   //---------------------------------------------------------------------------------------------
0685 
0686 
0687   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.3
0688   //---------------------------------------------------------------------------------------------
0689   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0690     if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
0691 
0692       s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0693         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0694         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.9593 + pow( qsqgev , 1 ) * -0.2063 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0695           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.939 + pow( qsqgev , 1 ) * 0.4481 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0696           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.122 + pow( qsqgev , 1 ) * -0.4946 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0697           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.8215 + pow( qsqgev , 1 ) * 0.1921 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0698           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0699           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0700           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0701 
0702     }
0703   }
0704   //---------------------------------------------------------------------------------------------
0705 
0706 
0707   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.4
0708   //---------------------------------------------------------------------------------------------
0709   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0710     if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
0711 
0712       s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0713         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0714         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.8411 + pow( qsqgev , 1 ) * -0.1796 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0715           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.849 + pow( qsqgev , 1 ) * 0.4218 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0716           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.072 + pow( qsqgev , 1 ) * -0.4768 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0717           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.8139 + pow( qsqgev , 1 ) * 0.1881 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0718           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0719           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0720           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0721 
0722     }
0723   }
0724   //---------------------------------------------------------------------------------------------
0725 
0726 
0727   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.5
0728   //---------------------------------------------------------------------------------------------
0729   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0730     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
0731 
0732       s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0733         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0734         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.7392 + pow( qsqgev , 1 ) * -0.1568 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0735           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.748 + pow( qsqgev , 1 ) * 0.3944 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0736           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 2.001 + pow( qsqgev , 1 ) * -0.4555 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0737           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.7962 + pow( qsqgev , 1 ) * 0.1822 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0738           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0739           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0740           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0741 
0742     }
0743   }
0744   //---------------------------------------------------------------------------------------------
0745 
0746 
0747   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.6
0748   //---------------------------------------------------------------------------------------------
0749   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0750     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
0751 
0752       s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0753         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0754         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.6512 + pow( qsqgev , 1 ) * -0.1372 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0755           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.641 + pow( qsqgev , 1 ) * 0.3667 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0756           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.915 + pow( qsqgev , 1 ) * -0.4318 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0757           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.7707 + pow( qsqgev , 1 ) * 0.1747 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0758           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0759           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0760           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0761 
0762     }
0763   }
0764   //---------------------------------------------------------------------------------------------
0765 
0766 
0767   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.7
0768   //---------------------------------------------------------------------------------------------
0769   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0770     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
0771 
0772       s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0773         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0774         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.5749 + pow( qsqgev , 1 ) * -0.1203 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0775           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.533 + pow( qsqgev , 1 ) * 0.3395 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0776           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.82 + pow( qsqgev , 1 ) * -0.4068 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0777           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.7398 + pow( qsqgev , 1 ) * 0.1663 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0778           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0779           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0780           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0781 
0782     }
0783   }
0784   //---------------------------------------------------------------------------------------------
0785 
0786 
0787   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.8
0788   //---------------------------------------------------------------------------------------------
0789   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0790     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
0791 
0792       s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0793         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0794         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.5085 + pow( qsqgev , 1 ) * -0.1058 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0795           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.426 + pow( qsqgev , 1 ) * 0.3132 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0796           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.719 + pow( qsqgev , 1 ) * -0.3813 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0797           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.7052 + pow( qsqgev , 1 ) * 0.1572 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0798           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0799           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0800           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0801 
0802     }
0803   }
0804   //---------------------------------------------------------------------------------------------
0805 
0806 
0807   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 2.9
0808   //---------------------------------------------------------------------------------------------
0809   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0810     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
0811 
0812       s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0813         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0814         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.4508 + pow( qsqgev , 1 ) * -0.09318 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0815           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.323 + pow( qsqgev , 1 ) * 0.2881 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0816           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.617 + pow( qsqgev , 1 ) * -0.3558 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0817           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.6684 + pow( qsqgev , 1 ) * 0.1478 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0818           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0819           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0820           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0821 
0822     }
0823   }
0824   //---------------------------------------------------------------------------------------------
0825 
0826 
0827   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 3
0828   //---------------------------------------------------------------------------------------------
0829   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0830     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
0831 
0832       s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0833         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0834         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.4004 + pow( qsqgev , 1 ) * -0.08226 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0835           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.224 + pow( qsqgev , 1 ) * 0.2645 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0836           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.516 + pow( qsqgev , 1 ) * -0.3309 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0837           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.6306 + pow( qsqgev , 1 ) * 0.1384 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0838           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0839           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0840           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0841 
0842     }
0843   }
0844   //---------------------------------------------------------------------------------------------
0845 
0846 
0847   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 3.1
0848   //---------------------------------------------------------------------------------------------
0849   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0850     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
0851 
0852       s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0853         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0854         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.3564 + pow( qsqgev , 1 ) * -0.07279 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0855           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.13 + pow( qsqgev , 1 ) * 0.2424 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0856           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.416 + pow( qsqgev , 1 ) * -0.3069 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0857           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.5927 + pow( qsqgev , 1 ) * 0.1291 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0858           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0859           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0860           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0861 
0862     }
0863   }
0864   //---------------------------------------------------------------------------------------------
0865 
0866 
0867   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 3.2
0868   //---------------------------------------------------------------------------------------------
0869   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0870     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
0871 
0872       s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0873         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0874         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.3178 + pow( qsqgev , 1 ) * -0.06454 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0875           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -1.043 + pow( qsqgev , 1 ) * 0.222 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0876           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.32 + pow( qsqgev , 1 ) * -0.2841 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0877           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.5554 + pow( qsqgev , 1 ) * 0.1201 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0878           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0879           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0880           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0881 
0882     }
0883   }
0884   //---------------------------------------------------------------------------------------------
0885 
0886 
0887   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 3.3
0888   //---------------------------------------------------------------------------------------------
0889   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0890     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
0891 
0892       s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0893         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0894         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.284 + pow( qsqgev , 1 ) * -0.05735 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0895           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.9605 + pow( qsqgev , 1 ) * 0.2031 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0896           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.229 + pow( qsqgev , 1 ) * -0.2625 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0897           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.5193 + pow( qsqgev , 1 ) * 0.1115 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0898           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0899           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0900           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0901 
0902     }
0903   }
0904   //---------------------------------------------------------------------------------------------
0905 
0906 
0907   //Paramaterization fo Sigma_T for Qsq: 3  to  4  at w = 3.4
0908   //---------------------------------------------------------------------------------------------
0909   if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) )  { 
0910     if ( wgev >= 3.4 )  { 
0911 
0912       s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0913         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0914         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2543 + pow( qsqgev , 1 ) * -0.05108 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0915           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.8844 + pow( qsqgev , 1 ) * 0.1857 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0916           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 1.142 + pow( qsqgev , 1 ) * -0.2424 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0917           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.4847 + pow( qsqgev , 1 ) * 0.1034 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0918           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0919           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0920           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0921 
0922     }
0923   }
0924   //---------------------------------------------------------------------------------------------
0925 
0926 
0927    //---------------------------------------------------------------------------------------------//
0928    //                                                                                             //
0929    //                                                                                             //
0930    //                                                                                             //
0931    //                                                                                             //
0932    //                                                                                             //
0933    //                                                                                             //
0934    //                                                                                             //
0935    //                                                                                             //
0936    //                                                                                             //
0937    //                                                                                             //
0938    //                                                                                             //
0939    //                                                                                             //
0940    //---------------------------------------------------------------------------------------------//
0941   //---------------------------------------------------------------------------------------------
0942   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2
0943   //---------------------------------------------------------------------------------------------
0944   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0945     if ( wgev <= 2.1 )  { 
0946 
0947       s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0948         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0949         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.6873 + pow( qsqgev , 1 ) * -0.1265 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0950           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.5864 + pow( qsqgev , 1 ) * 0.1344 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0951           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.504 + pow( qsqgev , 1 ) * -0.1185 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0952           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.185 + pow( qsqgev , 1 ) * 0.04164 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0953           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0954           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0955           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0956 
0957     }
0958   }
0959   //---------------------------------------------------------------------------------------------
0960 
0961 
0962   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.1
0963   //---------------------------------------------------------------------------------------------
0964   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0965     if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
0966 
0967       s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0968         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0969         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.6116 + pow( qsqgev , 1 ) * -0.1108 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0970           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.6799 + pow( qsqgev , 1 ) * 0.1423 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0971           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.632 + pow( qsqgev , 1 ) * -0.1333 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0972           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2289 + pow( qsqgev , 1 ) * 0.04735 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0973           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0974           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0975           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0976 
0977     }
0978   }
0979   //---------------------------------------------------------------------------------------------
0980 
0981 
0982   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.2
0983   //---------------------------------------------------------------------------------------------
0984   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
0985     if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
0986 
0987       s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
0988         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
0989         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.545 + pow( qsqgev , 1 ) * -0.09761 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0990           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.7289 + pow( qsqgev , 1 ) * 0.1454 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0991           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7123 + pow( qsqgev , 1 ) * -0.1425 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0992           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2602 + pow( qsqgev , 1 ) * 0.05152 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0993           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0994           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
0995           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
0996 
0997     }
0998   }
0999   //---------------------------------------------------------------------------------------------
1000 
1001 
1002   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.3
1003   //---------------------------------------------------------------------------------------------
1004   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1005     if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
1006 
1007       s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1008         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1009         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.4866 + pow( qsqgev , 1 ) * -0.08635 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1010           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.75 + pow( qsqgev , 1 ) * 0.1453 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1011           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7611 + pow( qsqgev , 1 ) * -0.1475 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1012           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2822 + pow( qsqgev , 1 ) * 0.05443 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1013           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1014           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1015           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1016 
1017     }
1018   }
1019   //---------------------------------------------------------------------------------------------
1020 
1021 
1022   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.4
1023   //---------------------------------------------------------------------------------------------
1024   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1025     if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
1026 
1027       s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1028         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1029         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.4352 + pow( qsqgev , 1 ) * -0.07663 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1030           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.7528 + pow( qsqgev , 1 ) * 0.1428 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1031           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7879 + pow( qsqgev , 1 ) * -0.1496 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1032           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.297 + pow( qsqgev , 1 ) * 0.05628 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1033           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1034           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1035           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1036 
1037     }
1038   }
1039   //---------------------------------------------------------------------------------------------
1040 
1041 
1042   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.5
1043   //---------------------------------------------------------------------------------------------
1044   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1045     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1046 
1047       s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1048         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1049         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.3898 + pow( qsqgev , 1 ) * -0.06816 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1050           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.743 + pow( qsqgev , 1 ) * 0.1388 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1051           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7985 + pow( qsqgev , 1 ) * -0.1494 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1052           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.3058 + pow( qsqgev , 1 ) * 0.05719 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1053           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1054           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1055           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1056 
1057     }
1058   }
1059   //---------------------------------------------------------------------------------------------
1060 
1061 
1062   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.6
1063   //---------------------------------------------------------------------------------------------
1064   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1065     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1066 
1067       s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1068         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1069         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.3497 + pow( qsqgev , 1 ) * -0.06075 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1070           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.7247 + pow( qsqgev , 1 ) * 0.1337 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1071           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7972 + pow( qsqgev , 1 ) * -0.1474 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1072           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.3098 + pow( qsqgev , 1 ) * 0.05733 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1073           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1074           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1075           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1076 
1077     }
1078   }
1079   //---------------------------------------------------------------------------------------------
1080 
1081 
1082   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.7
1083   //---------------------------------------------------------------------------------------------
1084   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1085     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1086 
1087       s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1088         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1089         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.3141 + pow( qsqgev , 1 ) * -0.05423 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1090           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.7004 + pow( qsqgev , 1 ) * 0.1279 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1091           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7867 + pow( qsqgev , 1 ) * -0.1441 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1092           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.3097 + pow( qsqgev , 1 ) * 0.0568 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1093           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1094           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1095           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1096 
1097     }
1098   }
1099   //---------------------------------------------------------------------------------------------
1100 
1101 
1102   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.8
1103   //---------------------------------------------------------------------------------------------
1104   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1105     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1106 
1107       s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1108         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1109         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2824 + pow( qsqgev , 1 ) * -0.04848 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1110           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.6721 + pow( qsqgev , 1 ) * 0.1216 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1111           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7692 + pow( qsqgev , 1 ) * -0.1397 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1112           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.3063 + pow( qsqgev , 1 ) * 0.05573 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1113           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1114           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1115           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1116 
1117     }
1118   }
1119   //---------------------------------------------------------------------------------------------
1120 
1121 
1122   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 2.9
1123   //---------------------------------------------------------------------------------------------
1124   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1125     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1126 
1127       s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1128         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1129         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2542 + pow( qsqgev , 1 ) * -0.04339 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1130           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.6411 + pow( qsqgev , 1 ) * 0.115 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1131           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7464 + pow( qsqgev , 1 ) * -0.1345 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1132           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.3003 + pow( qsqgev , 1 ) * 0.05424 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1133           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1134           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1135           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1136 
1137     }
1138   }
1139   //---------------------------------------------------------------------------------------------
1140 
1141 
1142   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 3
1143   //---------------------------------------------------------------------------------------------
1144   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1145     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
1146 
1147       s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1148         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1149         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2291 + pow( qsqgev , 1 ) * -0.03889 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1150           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.6087 + pow( qsqgev , 1 ) * 0.1083 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1151           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.7199 + pow( qsqgev , 1 ) * -0.1288 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1152           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2922 + pow( qsqgev , 1 ) * 0.05241 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1153           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1154           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1155           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1156 
1157     }
1158   }
1159   //---------------------------------------------------------------------------------------------
1160 
1161 
1162   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 3.1
1163   //---------------------------------------------------------------------------------------------
1164   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1165     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
1166 
1167       s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1168         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1169         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2067 + pow( qsqgev , 1 ) * -0.0349 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1170           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.5756 + pow( qsqgev , 1 ) * 0.1017 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1171           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.6907 + pow( qsqgev , 1 ) * -0.1227 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1172           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2826 + pow( qsqgev , 1 ) * 0.05035 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1173           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1174           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1175           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1176 
1177     }
1178   }
1179   //---------------------------------------------------------------------------------------------
1180 
1181 
1182   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 3.2
1183   //---------------------------------------------------------------------------------------------
1184   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1185     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
1186 
1187       s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1188         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1189         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1867 + pow( qsqgev , 1 ) * -0.03136 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1190           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.5426 + pow( qsqgev , 1 ) * 0.09518 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1191           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.6597 + pow( qsqgev , 1 ) * -0.1164 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1192           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2719 + pow( qsqgev , 1 ) * 0.04813 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1193           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1194           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1195           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1196 
1197     }
1198   }
1199   //---------------------------------------------------------------------------------------------
1200 
1201 
1202   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 3.3
1203   //---------------------------------------------------------------------------------------------
1204   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1205     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1206 
1207       s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1208         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1209         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1688 + pow( qsqgev , 1 ) * -0.02821 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1210           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.5101 + pow( qsqgev , 1 ) * 0.08887 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1211           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.6278 + pow( qsqgev , 1 ) * -0.1101 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1212           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2604 + pow( qsqgev , 1 ) * 0.0458 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1213           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1214           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1215           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1216 
1217     }
1218   }
1219   //---------------------------------------------------------------------------------------------
1220 
1221 
1222   //Paramaterization fo Sigma_T for Qsq: 4  to  5  at w = 3.4
1223   //---------------------------------------------------------------------------------------------
1224   if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) )  { 
1225     if ( wgev >= 3.4 )  { 
1226 
1227       s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1228         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1229         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1528 + pow( qsqgev , 1 ) * -0.02541 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1230           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.4784 + pow( qsqgev , 1 ) * 0.08282 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1231           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.5956 + pow( qsqgev , 1 ) * -0.1038 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1232           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.2484 + pow( qsqgev , 1 ) * 0.04344 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1233           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1234           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1235           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1236 
1237     }
1238   }
1239   //---------------------------------------------------------------------------------------------
1240 
1241 
1242    //---------------------------------------------------------------------------------------------//
1243    //                                                                                             //
1244    //                                                                                             //
1245    //                                                                                             //
1246    //                                                                                             //
1247    //                                                                                             //
1248    //                                                                                             //
1249    //                                                                                             //
1250    //                                                                                             //
1251    //                                                                                             //
1252    //                                                                                             //
1253    //                                                                                             //
1254    //                                                                                             //
1255    //---------------------------------------------------------------------------------------------//
1256   //---------------------------------------------------------------------------------------------
1257   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2
1258   //---------------------------------------------------------------------------------------------
1259   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1260     if ( wgev <= 2.1 )  { 
1261 
1262       s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1263         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1264         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.3192 + pow( qsqgev , 1 ) * -0.05197 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1265           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.05751 + pow( qsqgev , 1 ) * 0.02695 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1266           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.004748 + pow( qsqgev , 1 ) * -0.017 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1267           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.01052 + pow( qsqgev , 1 ) * 0.006129 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1268           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1269           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1270           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1271 
1272     }
1273   }
1274   //---------------------------------------------------------------------------------------------
1275 
1276 
1277   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.1
1278   //---------------------------------------------------------------------------------------------
1279   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1280     if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
1281 
1282       s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1283         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1284         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2946 + pow( qsqgev , 1 ) * -0.04664 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1285           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1614 + pow( qsqgev , 1 ) * 0.03711 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1286           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1228 + pow( qsqgev , 1 ) * -0.02987 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1287           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.04678 + pow( qsqgev , 1 ) * 0.01031 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1288           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1289           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1290           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1291 
1292     }
1293   }
1294   //---------------------------------------------------------------------------------------------
1295 
1296 
1297   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.2
1298   //---------------------------------------------------------------------------------------------
1299   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1300     if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
1301 
1302       s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1303         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1304         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2697 + pow( qsqgev , 1 ) * -0.04193 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1305           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2249 + pow( qsqgev , 1 ) * 0.04321 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1306           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1995 + pow( qsqgev , 1 ) * -0.03831 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1307           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.0723 + pow( qsqgev , 1 ) * 0.01333 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1308           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1309           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1310           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1311 
1312     }
1313   }
1314   //---------------------------------------------------------------------------------------------
1315 
1316 
1317   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.3
1318   //---------------------------------------------------------------------------------------------
1319   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1320     if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
1321 
1322       s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1323         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1324         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2464 + pow( qsqgev , 1 ) * -0.03778 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1325           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2647 + pow( qsqgev , 1 ) * 0.04686 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1326           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.251 + pow( qsqgev , 1 ) * -0.044 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1327           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.09101 + pow( qsqgev , 1 ) * 0.01559 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1328           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1329           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1330           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1331 
1332     }
1333   }
1334   //---------------------------------------------------------------------------------------------
1335 
1336 
1337   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.4
1338   //---------------------------------------------------------------------------------------------
1339   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1340     if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
1341 
1342       s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1343         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1344         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2248 + pow( qsqgev , 1 ) * -0.0341 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1345           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2892 + pow( qsqgev , 1 ) * 0.04888 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1346           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.2861 + pow( qsqgev , 1 ) * -0.04782 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1347           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.105 + pow( qsqgev , 1 ) * 0.01729 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1348           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1349           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1350           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1351 
1352     }
1353   }
1354   //---------------------------------------------------------------------------------------------
1355 
1356 
1357   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.5
1358   //---------------------------------------------------------------------------------------------
1359   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1360     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1361 
1362       s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1363         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1364         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.2052 + pow( qsqgev , 1 ) * -0.03084 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1365           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.3037 + pow( qsqgev , 1 ) * 0.0498 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1366           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3102 + pow( qsqgev , 1 ) * -0.05037 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1367           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1155 + pow( qsqgev , 1 ) * 0.01857 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1368           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1369           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1370           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1371 
1372     }
1373   }
1374   //---------------------------------------------------------------------------------------------
1375 
1376 
1377   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.6
1378   //---------------------------------------------------------------------------------------------
1379   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1380     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1381 
1382       s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1383         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1384         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1873 + pow( qsqgev , 1 ) * -0.02794 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1385           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.3109 + pow( qsqgev , 1 ) * 0.04992 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1386           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3261 + pow( qsqgev , 1 ) * -0.05193 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1387           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1232 + pow( qsqgev , 1 ) * 0.01949 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1388           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1389           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1390           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1391 
1392     }
1393   }
1394   //---------------------------------------------------------------------------------------------
1395 
1396 
1397   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.7
1398   //---------------------------------------------------------------------------------------------
1399   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1400     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1401 
1402       s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1403         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1404         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1711 + pow( qsqgev , 1 ) * -0.02534 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1405           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.313 + pow( qsqgev , 1 ) * 0.04945 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1406           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3359 + pow( qsqgev , 1 ) * -0.05273 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1407           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1287 + pow( qsqgev , 1 ) * 0.02012 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1408           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1409           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1410           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1411 
1412     }
1413   }
1414   //---------------------------------------------------------------------------------------------
1415 
1416 
1417   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.8
1418   //---------------------------------------------------------------------------------------------
1419   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1420     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1421 
1422       s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1423         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1424         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1563 + pow( qsqgev , 1 ) * -0.023 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1425           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.3111 + pow( qsqgev , 1 ) * 0.04851 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1426           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3408 + pow( qsqgev , 1 ) * -0.05291 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1427           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1323 + pow( qsqgev , 1 ) * 0.02048 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1428           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1429           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1430           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1431 
1432     }
1433   }
1434   //---------------------------------------------------------------------------------------------
1435 
1436 
1437   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 2.9
1438   //---------------------------------------------------------------------------------------------
1439   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1440     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1441 
1442       s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1443         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1444         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1428 + pow( qsqgev , 1 ) * -0.02089 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1445           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.3062 + pow( qsqgev , 1 ) * 0.04723 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1446           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3418 + pow( qsqgev , 1 ) * -0.05257 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1447           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1343 + pow( qsqgev , 1 ) * 0.02061 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1448           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1449           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1450           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1451 
1452     }
1453   }
1454   //---------------------------------------------------------------------------------------------
1455 
1456 
1457   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 3
1458   //---------------------------------------------------------------------------------------------
1459   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1460     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
1461 
1462       s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1463         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1464         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1306 + pow( qsqgev , 1 ) * -0.01899 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1465           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.299 + pow( qsqgev , 1 ) * 0.04568 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1466           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3397 + pow( qsqgev , 1 ) * -0.05182 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1467           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1349 + pow( qsqgev , 1 ) * 0.02055 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1468           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1469           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1470           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1471 
1472     }
1473   }
1474   //---------------------------------------------------------------------------------------------
1475 
1476 
1477   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 3.1
1478   //---------------------------------------------------------------------------------------------
1479   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1480     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
1481 
1482       s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1483         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1484         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1194 + pow( qsqgev , 1 ) * -0.01727 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1485           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2901 + pow( qsqgev , 1 ) * 0.04394 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1486           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3349 + pow( qsqgev , 1 ) * -0.05072 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1487           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1342 + pow( qsqgev , 1 ) * 0.02032 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1488           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1489           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1490           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1491 
1492     }
1493   }
1494   //---------------------------------------------------------------------------------------------
1495 
1496 
1497   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 3.2
1498   //---------------------------------------------------------------------------------------------
1499   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1500     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
1501 
1502       s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1503         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1504         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1092 + pow( qsqgev , 1 ) * -0.01572 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1505           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2799 + pow( qsqgev , 1 ) * 0.04207 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1506           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.328 + pow( qsqgev , 1 ) * -0.04935 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1507           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1326 + pow( qsqgev , 1 ) * 0.01995 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1508           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1509           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1510           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1511 
1512     }
1513   }
1514   //---------------------------------------------------------------------------------------------
1515 
1516 
1517   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 3.3
1518   //---------------------------------------------------------------------------------------------
1519   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1520     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1521 
1522       s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1523         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1524         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.09998 + pow( qsqgev , 1 ) * -0.01431 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1525           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2689 + pow( qsqgev , 1 ) * 0.04012 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1526           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3196 + pow( qsqgev , 1 ) * -0.04777 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1527           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1302 + pow( qsqgev , 1 ) * 0.01947 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1528           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1529           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1530           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1531 
1532     }
1533   }
1534   //---------------------------------------------------------------------------------------------
1535 
1536 
1537   //Paramaterization fo Sigma_T for Qsq: 5  to  6  at w = 3.4
1538   //---------------------------------------------------------------------------------------------
1539   if ( ( qsqgev >= 5 ) && ( qsqgev <= 6 ) )  { 
1540     if ( wgev >= 3.4 )  { 
1541 
1542       s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1543         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1544         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.09153 + pow( qsqgev , 1 ) * -0.01304 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1545           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.2573 + pow( qsqgev , 1 ) * 0.03813 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1546           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.3098 + pow( qsqgev , 1 ) * -0.04603 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1547           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.1271 + pow( qsqgev , 1 ) * 0.01889 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1548           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1549           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1550           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1551 
1552     }
1553   }
1554   //---------------------------------------------------------------------------------------------
1555 
1556 
1557    //---------------------------------------------------------------------------------------------//
1558    //                                                                                             //
1559    //                                                                                             //
1560    //                                                                                             //
1561    //                                                                                             //
1562    //                                                                                             //
1563    //                                                                                             //
1564    //                                                                                             //
1565    //                                                                                             //
1566    //                                                                                             //
1567    //                                                                                             //
1568    //                                                                                             //
1569    //                                                                                             //
1570    //---------------------------------------------------------------------------------------------//
1571   //---------------------------------------------------------------------------------------------
1572   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2
1573   //---------------------------------------------------------------------------------------------
1574   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1575     if ( wgev <= 2.1 )  { 
1576 
1577       s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1578         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1579         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1287 + pow( qsqgev , 1 ) * -0.01998 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1580           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * 0.1505 + pow( qsqgev , 1 ) * -0.008036 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1581           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * -0.1683 + pow( qsqgev , 1 ) * 0.01216 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1582           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * 0.0459 + pow( qsqgev , 1 ) * -0.003386 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1583           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1584           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1585           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1586 
1587     }
1588   }
1589   //---------------------------------------------------------------------------------------------
1590 
1591 
1592   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.1
1593   //---------------------------------------------------------------------------------------------
1594   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1595     if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) )  { 
1596 
1597       s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1598         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1599         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1341 + pow( qsqgev , 1 ) * -0.01966 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1600           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * 0.0373 + pow( qsqgev , 1 ) * 0.003599 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1601           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * -0.05297 + pow( qsqgev , 1 ) * -0.0001854 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1602           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * 0.01242 + pow( qsqgev , 1 ) * 0.0003032 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1603           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1604           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1605           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1606 
1607     }
1608   }
1609   //---------------------------------------------------------------------------------------------
1610 
1611 
1612   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.2
1613   //---------------------------------------------------------------------------------------------
1614   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1615     if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) )  { 
1616 
1617       s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1618         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1619         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.129 + pow( qsqgev , 1 ) * -0.01827 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1620           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.02419 + pow( qsqgev , 1 ) * 0.009363 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1621           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.01269 + pow( qsqgev , 1 ) * -0.006786 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1622           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.007635 + pow( qsqgev , 1 ) * 0.002411 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1623           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1624           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1625           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1626 
1627     }
1628   }
1629   //---------------------------------------------------------------------------------------------
1630 
1631 
1632   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.3
1633   //---------------------------------------------------------------------------------------------
1634   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1635     if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) )  { 
1636 
1637       s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1638         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1639         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1219 + pow( qsqgev , 1 ) * -0.01685 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1640           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.06392 + pow( qsqgev , 1 ) * 0.01302 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1641           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.05636 + pow( qsqgev , 1 ) * -0.01116 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1642           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.02181 + pow( qsqgev , 1 ) * 0.003912 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1643           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1644           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1645           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1646 
1647     }
1648   }
1649   //---------------------------------------------------------------------------------------------
1650 
1651 
1652   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.4
1653   //---------------------------------------------------------------------------------------------
1654   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1655     if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) )  { 
1656 
1657       s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1658         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1659         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1143 + pow( qsqgev , 1 ) * -0.01551 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1660           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.09088 + pow( qsqgev , 1 ) * 0.01546 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1661           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.08695 + pow( qsqgev , 1 ) * -0.01424 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1662           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.0324 + pow( qsqgev , 1 ) * 0.005049 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1663           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1664           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1665           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1666 
1667     }
1668   }
1669   //---------------------------------------------------------------------------------------------
1670 
1671 
1672   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.5
1673   //---------------------------------------------------------------------------------------------
1674   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1675     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1676 
1677       s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1678         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1679         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.1066 + pow( qsqgev , 1 ) * -0.01428 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1680           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1095 + pow( qsqgev , 1 ) * 0.0171 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1681           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.109 + pow( qsqgev , 1 ) * -0.01646 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1682           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.04052 + pow( qsqgev , 1 ) * 0.005931 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1683           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1684           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1685           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1686 
1687     }
1688   }
1689   //---------------------------------------------------------------------------------------------
1690 
1691 
1692   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.6
1693   //---------------------------------------------------------------------------------------------
1694   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1695     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1696 
1697       s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1698         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1699         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.09933 + pow( qsqgev , 1 ) * -0.01315 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1700           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1226 + pow( qsqgev , 1 ) * 0.01821 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1701           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1253 + pow( qsqgev , 1 ) * -0.0181 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1702           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.04697 + pow( qsqgev , 1 ) * 0.006637 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1703           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1704           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1705           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1706 
1707     }
1708   }
1709   //---------------------------------------------------------------------------------------------
1710 
1711 
1712   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.7
1713   //---------------------------------------------------------------------------------------------
1714   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1715     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1716 
1717       s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1718         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1719         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.09237 + pow( qsqgev , 1 ) * -0.01211 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1720           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1316 + pow( qsqgev , 1 ) * 0.01891 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1721           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1374 + pow( qsqgev , 1 ) * -0.0193 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1722           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.05202 + pow( qsqgev , 1 ) * 0.00719 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1723           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1724           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1725           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1726 
1727     }
1728   }
1729   //---------------------------------------------------------------------------------------------
1730 
1731 
1732   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.8
1733   //---------------------------------------------------------------------------------------------
1734   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1735     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1736 
1737       s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1738         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1739         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.08583 + pow( qsqgev , 1 ) * -0.01116 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1740           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1375 + pow( qsqgev , 1 ) * 0.01931 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1741           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1464 + pow( qsqgev , 1 ) * -0.02018 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1742           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.05605 + pow( qsqgev , 1 ) * 0.007628 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1743           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1744           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1745           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1746 
1747     }
1748   }
1749   //---------------------------------------------------------------------------------------------
1750 
1751 
1752   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 2.9
1753   //---------------------------------------------------------------------------------------------
1754   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1755     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1756 
1757       s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1758         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1759         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.07969 + pow( qsqgev , 1 ) * -0.01028 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1760           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1411 + pow( qsqgev , 1 ) * 0.01945 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1761           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1528 + pow( qsqgev , 1 ) * -0.02076 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1762           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.05914 + pow( qsqgev , 1 ) * 0.007956 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1763           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1764           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1765           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1766 
1767     }
1768   }
1769   //---------------------------------------------------------------------------------------------
1770 
1771 
1772   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 3
1773   //---------------------------------------------------------------------------------------------
1774   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1775     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
1776 
1777       s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1778         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1779         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.07395 + pow( qsqgev , 1 ) * -0.009479 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1780           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1427 + pow( qsqgev , 1 ) * 0.01939 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1781           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1572 + pow( qsqgev , 1 ) * -0.02112 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1782           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.06146 + pow( qsqgev , 1 ) * 0.008194 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1783           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1784           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1785           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1786 
1787     }
1788   }
1789   //---------------------------------------------------------------------------------------------
1790 
1791 
1792   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 3.1
1793   //---------------------------------------------------------------------------------------------
1794   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1795     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
1796 
1797       s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1798         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1799         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.06859 + pow( qsqgev , 1 ) * -0.008737 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1800           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1428 + pow( qsqgev , 1 ) * 0.01917 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1801           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1598 + pow( qsqgev , 1 ) * -0.02127 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1802           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.06307 + pow( qsqgev , 1 ) * 0.008344 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1803           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1804           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1805           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1806 
1807     }
1808   }
1809   //---------------------------------------------------------------------------------------------
1810 
1811 
1812   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 3.2
1813   //---------------------------------------------------------------------------------------------
1814   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1815     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
1816 
1817       s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1818         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1819         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.06359 + pow( qsqgev , 1 ) * -0.008053 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1820           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1417 + pow( qsqgev , 1 ) * 0.01883 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1821           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.161 + pow( qsqgev , 1 ) * -0.02125 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1822           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.06408 + pow( qsqgev , 1 ) * 0.008421 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1823           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1824           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1825           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1826 
1827     }
1828   }
1829   //---------------------------------------------------------------------------------------------
1830 
1831 
1832   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 3.3
1833   //---------------------------------------------------------------------------------------------
1834   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1835     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
1836 
1837       s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1838         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1839         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.05894 + pow( qsqgev , 1 ) * -0.007422 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1840           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1395 + pow( qsqgev , 1 ) * 0.01837 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1841           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1608 + pow( qsqgev , 1 ) * -0.02108 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1842           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.06455 + pow( qsqgev , 1 ) * 0.00843 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1843           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1844           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1845           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1846 
1847     }
1848   }
1849   //---------------------------------------------------------------------------------------------
1850 
1851 
1852   //Paramaterization fo Sigma_T for Qsq: 6  to  7  at w = 3.4
1853   //---------------------------------------------------------------------------------------------
1854   if ( ( qsqgev >= 6 ) && ( qsqgev <= 7 ) )  { 
1855     if ( wgev >= 3.4 )  { 
1856 
1857       s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1858         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1859         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.05461 + pow( qsqgev , 1 ) * -0.006841 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1860           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.1365 + pow( qsqgev , 1 ) * 0.01783 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1861           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.1596 + pow( qsqgev , 1 ) * -0.02078 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1862           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.06454 + pow( qsqgev , 1 ) * 0.008382 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1863           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1864           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1865           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1866 
1867     }
1868   }
1869   //---------------------------------------------------------------------------------------------
1870 
1871 
1872    //---------------------------------------------------------------------------------------------//
1873    //                                                                                             //
1874    //                                                                                             //
1875    //                                                                                             //
1876    //                                                                                             //
1877    //                                                                                             //
1878    //                                                                                             //
1879    //                                                                                             //
1880    //                                                                                             //
1881    //                                                                                             //
1882    //                                                                                             //
1883    //                                                                                             //
1884    //                                                                                             //
1885    //---------------------------------------------------------------------------------------------//
1886   //---------------------------------------------------------------------------------------------
1887   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 2.4
1888   //---------------------------------------------------------------------------------------------
1889   if ( qsqgev >= 7 )  { 
1890     if ( wgev < 2.5 )  { 
1891 
1892       s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1893         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1894         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.05511 + pow( qsqgev , 1 ) * -0.007004 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1895           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.006249 + pow( qsqgev , 1 ) * 0.00326 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1896           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.007846 + pow( qsqgev , 1 ) * -0.002827 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1897           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.004861 + pow( qsqgev , 1 ) * 0.001075 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1898           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1899           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1900           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1901 
1902     }
1903   }
1904   //---------------------------------------------------------------------------------------------
1905 
1906 
1907   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 2.5
1908   //---------------------------------------------------------------------------------------------
1909   if ( qsqgev >= 7 )  { 
1910     if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) )  { 
1911 
1912       s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1913         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1914         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.05332 + pow( qsqgev , 1 ) * -0.006608 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1915           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.0239 + pow( qsqgev , 1 ) * 0.004764 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1916           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.02592 + pow( qsqgev , 1 ) * -0.004478 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1917           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.01087 + pow( qsqgev , 1 ) * 0.001653 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1918           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1919           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1920           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1921 
1922     }
1923   }
1924   //---------------------------------------------------------------------------------------------
1925 
1926 
1927   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 2.6
1928   //---------------------------------------------------------------------------------------------
1929   if ( qsqgev >= 7 )  { 
1930     if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) )  { 
1931 
1932       s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1933         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1934         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.05109 + pow( qsqgev , 1 ) * -0.006211 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1935           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.03685 + pow( qsqgev , 1 ) * 0.005857 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1936           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.03938 + pow( qsqgev , 1 ) * -0.005712 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1937           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.01557 + pow( qsqgev , 1 ) * 0.002109 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1938           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1939           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1940           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1941 
1942     }
1943   }
1944   //---------------------------------------------------------------------------------------------
1945 
1946 
1947   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 2.7
1948   //---------------------------------------------------------------------------------------------
1949   if ( qsqgev >= 7 )  { 
1950     if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) )  { 
1951 
1952       s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1953         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1954         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.04867 + pow( qsqgev , 1 ) * -0.005826 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1955           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.04658 + pow( qsqgev , 1 ) * 0.006668 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1956           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.04974 + pow( qsqgev , 1 ) * -0.006666 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1957           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.01935 + pow( qsqgev , 1 ) * 0.00248 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1958           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1959           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1960           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1961 
1962     }
1963   }
1964   //---------------------------------------------------------------------------------------------
1965 
1966 
1967   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 2.8
1968   //---------------------------------------------------------------------------------------------
1969   if ( qsqgev >= 7 )  { 
1970     if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) )  { 
1971 
1972       s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1973         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1974         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.04618 + pow( qsqgev , 1 ) * -0.005458 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1975           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.05396 + pow( qsqgev , 1 ) * 0.00727 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1976           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.05787 + pow( qsqgev , 1 ) * -0.007417 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1977           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.02245 + pow( qsqgev , 1 ) * 0.002787 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1978           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1979           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1980           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
1981 
1982     }
1983   }
1984   //---------------------------------------------------------------------------------------------
1985 
1986 
1987   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 2.9
1988   //---------------------------------------------------------------------------------------------
1989   if ( qsqgev >= 7 )  { 
1990     if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) )  { 
1991 
1992       s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
1993         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
1994         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.0437 + pow( qsqgev , 1 ) * -0.005109 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1995           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.05954 + pow( qsqgev , 1 ) * 0.007709 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1996           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.06433 + pow( qsqgev , 1 ) * -0.008012 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1997           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.02503 + pow( qsqgev , 1 ) * 0.003042 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1998           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
1999           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2000           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
2001 
2002     }
2003   }
2004   //---------------------------------------------------------------------------------------------
2005 
2006 
2007   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 3
2008   //---------------------------------------------------------------------------------------------
2009   if ( qsqgev >= 7 )  { 
2010     if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) )  { 
2011 
2012       s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
2013         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
2014         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.04126 + pow( qsqgev , 1 ) * -0.004778 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2015           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.06369 + pow( qsqgev , 1 ) * 0.008018 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2016           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.06947 + pow( qsqgev , 1 ) * -0.008481 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2017           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.02716 + pow( qsqgev , 1 ) * 0.003253 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2018           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2019           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2020           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
2021 
2022     }
2023   }
2024   //---------------------------------------------------------------------------------------------
2025 
2026 
2027   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 3.1
2028   //---------------------------------------------------------------------------------------------
2029   if ( qsqgev >= 7 )  { 
2030     if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) )  { 
2031 
2032       s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
2033         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
2034         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.03888 + pow( qsqgev , 1 ) * -0.004466 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2035           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.06666 + pow( qsqgev , 1 ) * 0.008217 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2036           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.07349 + pow( qsqgev , 1 ) * -0.008841 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2037           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.02891 + pow( qsqgev , 1 ) * 0.003424 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2038           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2039           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2040           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
2041 
2042     }
2043   }
2044   //---------------------------------------------------------------------------------------------
2045 
2046 
2047   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 3.2
2048   //---------------------------------------------------------------------------------------------
2049   if ( qsqgev >= 7 )  { 
2050     if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) )  { 
2051 
2052       s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
2053         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
2054         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.03659 + pow( qsqgev , 1 ) * -0.004172 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2055           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.06867 + pow( qsqgev , 1 ) * 0.008324 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2056           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.07659 + pow( qsqgev , 1 ) * -0.009106 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2057           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.03032 + pow( qsqgev , 1 ) * 0.00356 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2058           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2059           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2060           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
2061 
2062     }
2063   }
2064   //---------------------------------------------------------------------------------------------
2065 
2066 
2067   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 3.3
2068   //---------------------------------------------------------------------------------------------
2069   if ( qsqgev >= 7 )  { 
2070     if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) )  { 
2071 
2072       s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
2073         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
2074         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.03439 + pow( qsqgev , 1 ) * -0.003895 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2075           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.06986 + pow( qsqgev , 1 ) * 0.008353 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2076           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.07887 + pow( qsqgev , 1 ) * -0.009287 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2077           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.03142 + pow( qsqgev , 1 ) * 0.003663 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2078           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2079           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2080           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
2081 
2082     }
2083   }
2084   //---------------------------------------------------------------------------------------------
2085 
2086 
2087   //Paramaterization fo Sigma_T for Qsq: 7  to  8  at w = 3.4
2088   //---------------------------------------------------------------------------------------------
2089   if ( qsqgev >= 7 )  { 
2090     if ( wgev >= 3.4 )  { 
2091 
2092       s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) / 
2093         pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) * 
2094         ( (pow( tgev  , 0 ) * (  pow( qsqgev , 0 ) * 0.03229 + pow( qsqgev , 1 ) * -0.003634 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2095           (  pow( tgev  , 1 ) * (  pow( qsqgev , 0 ) * -0.07035 + pow( qsqgev , 1 ) * 0.008315 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2096           (  pow( tgev  , 2 ) * (  pow( qsqgev , 0 ) * 0.08041 + pow( qsqgev , 1 ) * -0.00939 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2097           (  pow( tgev  , 3 ) * (  pow( qsqgev , 0 ) * -0.03225 + pow( qsqgev , 1 ) * 0.003735 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2098           (  pow( tgev  , 4 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2099           (  pow( tgev  , 5 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) + 
2100           (  pow( tgev  , 6 ) * (  pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
2101 
2102     }
2103   }
2104   //---------------------------------------------------------------------------------------------
2105 
2106 
2107   /*****************************************************************************************************************************************/
2108   /*****************************************************************************************************************************************/
2109   /*****************************************************************************************************************************************/
2110   /*****************************************************************************************************************************************/
2111   /*****************************************************************************************************************************************/
2112   /*****************************************************************************************************************************************/
2113   /*****************************************************************************************************************************************/
2114   /*****************************************************************************************************************************************/
2115 
2116 
2117   return s; 
2118 }
2119 
2120
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!