-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathSigmaT.hxx
2120 lines (1707 loc) · 162 KB
/
SigmaT.hxx
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
#include <iostream>
#include <fstream>
#include <stdio.h>
#include <cstdlib>
#include <stdlib.h>
#include <iomanip>
#include <stdio.h>
#include <math.h>
#include <cmath>
using std::setw;
using std::setprecision;
using std::cout;
using std::cin;
using std::endl;
using namespace std;
/* double fSigmaT; */
double MySigmaT( double qsqgev, double tgev, double wgev ) {
double s = 0.0;
/* double wgev = par[0]; */
/* double qsqgev = par[1]; */
/* double tgev = x[0]; */
//---------------------------------------------------------------------------------------------//
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
//---------------------------------------------------------------------------------------------//
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( wgev <= 2.1 ) {
s = ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 2.546 + pow( qsqgev , 1 ) * -1.074 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -45.04 + pow( qsqgev , 1 ) * 60.66 + pow( qsqgev , 2 ) * -22.69 ) ) +
( pow( qsqgev , 0 ) * 6.504 + pow( qsqgev , 1 ) * -7.199 + pow( qsqgev , 2 ) * 2.295 ) +
( tgev * ( pow( qsqgev , 0 ) * -4.284 + pow( qsqgev , 1 ) * 5.635 + pow( qsqgev , 2 ) * -1.911 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.1
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) ) {
s = ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 2.332 + pow( qsqgev , 1 ) * -1.054 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -40.39 + pow( qsqgev , 1 ) * 51.49 + pow( qsqgev , 2 ) * -18.97 ) ) +
( pow( qsqgev , 0 ) * 4.724 + pow( qsqgev , 1 ) * -5.13 + pow( qsqgev , 2 ) * 1.628 ) +
( tgev * ( pow( qsqgev , 0 ) * -3.525 + pow( qsqgev , 1 ) * 4.489 + pow( qsqgev , 2 ) * -1.495 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.2
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) ) {
s = ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 2.12 + pow( qsqgev , 1 ) * -1.023 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -34.61 + pow( qsqgev , 1 ) * 40.38 + pow( qsqgev , 2 ) * -14.55 ) ) +
( pow( qsqgev , 0 ) * 3.461 + pow( qsqgev , 1 ) * -3.685 + pow( qsqgev , 2 ) * 1.166 ) +
( tgev * ( pow( qsqgev , 0 ) * -2.799 + pow( qsqgev , 1 ) * 3.458 + pow( qsqgev , 2 ) * -1.135 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.3
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) ) {
s = ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 1.91 + pow( qsqgev , 1 ) * -0.9857 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -27.26 + pow( qsqgev , 1 ) * 26.45 + pow( qsqgev , 2 ) * -9.059 ) ) +
( pow( qsqgev , 0 ) * 2.538 + pow( qsqgev , 1 ) * -2.63 + pow( qsqgev , 2 ) * 0.8264 ) +
( tgev * ( pow( qsqgev , 0 ) * -2.161 + pow( qsqgev , 1 ) * 2.574 + pow( qsqgev , 2 ) * -0.831 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.4
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) ) {
s = ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 1.694 + pow( qsqgev , 1 ) * -0.9393 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -18.64 + pow( qsqgev , 1 ) * 10.18 + pow( qsqgev , 2 ) * -2.648 ) ) +
( pow( qsqgev , 0 ) * 1.771 + pow( qsqgev , 1 ) * -1.691 + pow( qsqgev , 2 ) * 0.5114 ) +
( tgev * ( pow( qsqgev , 0 ) * -1.533 + pow( qsqgev , 1 ) * 1.671 + pow( qsqgev , 2 ) * -0.5125 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.5
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) ) {
s = ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 1.521 + pow( qsqgev , 1 ) * -0.9198 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -17.9 + pow( qsqgev , 1 ) * 8.338 + pow( qsqgev , 2 ) * -1.914 ) ) +
( pow( qsqgev , 0 ) * 1.326 + pow( qsqgev , 1 ) * -1.224 + pow( qsqgev , 2 ) * 0.3666 ) +
( tgev * ( pow( qsqgev , 0 ) * -1.184 + pow( qsqgev , 1 ) * 1.236 + pow( qsqgev , 2 ) * -0.3718 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.6
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) ) {
s = ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 1.358 + pow( qsqgev , 1 ) * -0.9031 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -17.73 + pow( qsqgev , 1 ) * 7.539 + pow( qsqgev , 2 ) * -1.597 ) ) +
( pow( qsqgev , 0 ) * 1.016 + pow( qsqgev , 1 ) * -0.9168 + pow( qsqgev , 2 ) * 0.2741 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.9309 + pow( qsqgev , 1 ) * 0.9422 + pow( qsqgev , 2 ) * -0.28 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.7
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) ) {
s = ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 1.202 + pow( qsqgev , 1 ) * -0.8882 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -17.69 + pow( qsqgev , 1 ) * 6.986 + pow( qsqgev , 2 ) * -1.374 ) ) +
( pow( qsqgev , 0 ) * 0.7906 + pow( qsqgev , 1 ) * -0.7023 + pow( qsqgev , 2 ) * 0.2104 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.7403 + pow( qsqgev , 1 ) * 0.7301 + pow( qsqgev , 2 ) * -0.2154 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.8
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) ) {
s = ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 1.052 + pow( qsqgev , 1 ) * -0.8745 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -17.73 + pow( qsqgev , 1 ) * 6.568 + pow( qsqgev , 2 ) * -1.2 ) ) +
( pow( qsqgev , 0 ) * 0.6239 + pow( qsqgev , 1 ) * -0.5481 + pow( qsqgev , 2 ) * 0.165 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.5943 + pow( qsqgev , 1 ) * 0.5737 + pow( qsqgev , 2 ) * -0.1686 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 2.9
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) ) {
s = ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 0.9084 + pow( qsqgev , 1 ) * -0.8628 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -17.94 + pow( qsqgev , 1 ) * 6.43 + pow( qsqgev , 2 ) * -1.13 ) ) +
( pow( qsqgev , 0 ) * 0.4999 + pow( qsqgev , 1 ) * -0.4371 + pow( qsqgev , 2 ) * 0.1326 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.483 + pow( qsqgev , 1 ) * 0.4591 + pow( qsqgev , 2 ) * -0.1349 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 3
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) ) {
s = ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 0.7705 + pow( qsqgev , 1 ) * -0.852 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -18.28 + pow( qsqgev , 1 ) * 6.522 + pow( qsqgev , 2 ) * -1.153 ) ) +
( pow( qsqgev , 0 ) * 0.4062 + pow( qsqgev , 1 ) * -0.3561 + pow( qsqgev , 2 ) * 0.1092 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.3972 + pow( qsqgev , 1 ) * 0.3742 + pow( qsqgev , 2 ) * -0.1106 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 3.1
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) ) {
s = ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 0.6384 + pow( qsqgev , 1 ) * -0.8428 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -18.69 + pow( qsqgev , 1 ) * 6.698 + pow( qsqgev , 2 ) * -1.192 ) ) +
( pow( qsqgev , 0 ) * 0.3342 + pow( qsqgev , 1 ) * -0.2948 + pow( qsqgev , 2 ) * 0.09148 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.33 + pow( qsqgev , 1 ) * 0.3095 + pow( qsqgev , 2 ) * -0.09206 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 3.2
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) ) {
s = ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 0.5101 + pow( qsqgev , 1 ) * -0.8344 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -19.11 + pow( qsqgev , 1 ) * 6.918 + pow( qsqgev , 2 ) * -1.25 ) ) +
( pow( qsqgev , 0 ) * 0.2778 + pow( qsqgev , 1 ) * -0.2475 + pow( qsqgev , 2 ) * 0.07769 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.2765 + pow( qsqgev , 1 ) * 0.259 + pow( qsqgev , 2 ) * -0.07767 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 3.3
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) ) {
s = ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 0.3852 + pow( qsqgev , 1 ) * -0.8252 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -19.58 + pow( qsqgev , 1 ) * 7.256 + pow( qsqgev , 2 ) * -1.365 ) ) +
( pow( qsqgev , 0 ) * 0.2334 + pow( qsqgev , 1 ) * -0.2109 + pow( qsqgev , 2 ) * 0.06704 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.2338 + pow( qsqgev , 1 ) * 0.2199 + pow( qsqgev , 2 ) * -0.06668 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 0.1 to 2 at w = 3.4
//---------------------------------------------------------------------------------------------
if ( qsqgev <= 2 ) {
if ( wgev >= 3.4 ) {
s = ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) *
( exp( ( pow( qsqgev , 0 ) * 0.2654 + pow( qsqgev , 1 ) * -0.8181 + pow( qsqgev , 2 ) * 0 ) +
tgev * ( pow( qsqgev , 0 ) * -20.09 + pow( qsqgev , 1 ) * 7.614 + pow( qsqgev , 2 ) * -1.474 ) ) +
( pow( qsqgev , 0 ) * 0.1979 + pow( qsqgev , 1 ) * -0.1816 + pow( qsqgev , 2 ) * 0.05842 ) +
( tgev * ( pow( qsqgev , 0 ) * -0.1993 + pow( qsqgev , 1 ) * 0.1887 + pow( qsqgev , 2 ) * -0.05786 ) ) );
}
}
//---------------------------------------------------------------------------------------------
//---------------------------------------------------------------------------------------------//
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
//---------------------------------------------------------------------------------------------//
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( wgev <= 2.1 ) {
s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.1
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) ) {
s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) ) {
s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.3
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) ) {
s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.4
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) ) {
s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.5
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) ) {
s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.6
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) ) {
s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.7
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) ) {
s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.8
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) ) {
s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 2.9
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) ) {
s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 3
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) ) {
s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 3.1
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) ) {
s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 3.2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) ) {
s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 3.3
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) ) {
s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 2 to 3 at w = 3.4
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 2 ) && ( qsqgev <= 3 ) ) {
if ( wgev >= 3.4 ) {
s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//---------------------------------------------------------------------------------------------//
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
//---------------------------------------------------------------------------------------------//
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( wgev <= 2.1 ) {
s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.1
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) ) {
s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) ) {
s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.3
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.3 ) && ( wgev <= 2.4 ) ) {
s = ( ( ( pow( ( pow( 2.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.4
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.4 ) && ( wgev <= 2.5 ) ) {
s = ( ( ( pow( ( pow( 2.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.5
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.5 ) && ( wgev <= 2.6 ) ) {
s = ( ( ( pow( ( pow( 2.5 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.6
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.6 ) && ( wgev <= 2.7 ) ) {
s = ( ( ( pow( ( pow( 2.6 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.7
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.7 ) && ( wgev <= 2.8 ) ) {
s = ( ( ( pow( ( pow( 2.7 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.8
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.8 ) && ( wgev <= 2.9 ) ) {
s = ( ( ( pow( ( pow( 2.8 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 2.9
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 2.9 ) && ( wgev <= 3 ) ) {
s = ( ( ( pow( ( pow( 2.9 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 3
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 3 ) && ( wgev <= 3.1 ) ) {
s = ( ( ( pow( ( pow( 3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 3.1
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 3.1 ) && ( wgev <= 3.2 ) ) {
s = ( ( ( pow( ( pow( 3.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 3.2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 3.2 ) && ( wgev <= 3.3 ) ) {
s = ( ( ( pow( ( pow( 3.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 3.3
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( ( wgev >= 3.3 ) && ( wgev <= 3.4 ) ) {
s = ( ( ( pow( ( pow( 3.3 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 3 to 4 at w = 3.4
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 3 ) && ( qsqgev <= 4 ) ) {
if ( wgev >= 3.4 ) {
s = ( ( ( pow( ( pow( 3.4 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//---------------------------------------------------------------------------------------------//
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
// //
//---------------------------------------------------------------------------------------------//
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 4 to 5 at w = 2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) ) {
if ( wgev <= 2.1 ) {
s = ( ( ( pow( ( pow( 2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 4 to 5 at w = 2.1
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) ) {
if ( ( wgev >= 2.1 ) && ( wgev <= 2.2 ) ) {
s = ( ( ( pow( ( pow( 2.1 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------
//Paramaterization fo Sigma_T for Qsq: 4 to 5 at w = 2.2
//---------------------------------------------------------------------------------------------
if ( ( qsqgev >= 4 ) && ( qsqgev <= 5 ) ) {
if ( ( wgev >= 2.2 ) && ( wgev <= 2.3 ) ) {
s = ( ( ( pow( ( pow( 2.2 , 2 ) - pow( 0.938 , 2 ) ) , 2 ) /
pow( ( pow( wgev , 2 ) - pow( 0.938 , 2 ) ) , 2 ) ) ) *
( (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 ) ) +
( 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 ) ) +
( 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 ) ) +
( 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 ) ) +
( pow( tgev , 4 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 5 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) +
( pow( tgev , 6 ) * ( pow( qsqgev , 0 ) * 0 + pow( qsqgev , 1 ) * 0 + pow( qsqgev , 2 ) * 0 + pow( qsqgev , 3 ) * 0 + pow( qsqgev , 4 ) * 0 ) ) ) );
}
}
//---------------------------------------------------------------------------------------------