Home Explore contacts
Sign In
PEParfenov
/
usefulmacros
Watch 1
Star 0
Fork 0
Files Issues 0 Pull Requests 0 Wiki
Branch: master
Branches Tags
usefulmacros
/
SigmaInelNN.C

SigmaInelNN.C 62 KB
Permalink History Raw

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137 
  1. double CMToPLab(double sNN = 4.5) {
    2. 

  2. 

  3.   double M_NUCL = 0.938;
    4. 

  4.   double sN = sNN/2.;
    5. 

  5.   double pN = TMath::Sqrt(sN*sN - M_NUCL*M_NUCL);
    6. 

  6. 

  7.   TLorentzVector tMom4(0, 0, -pN, sN);
    8. 

  8.   TVector3 vBoost{tMom4.BoostVector()};
    9. 

  9.   TLorentzVector pMom4(0, 0, pN, sN);
    10. 

  10.   pMom4.Boost(-vBoost);
    11. 

  11. 

  12.   return (double)pMom4.Pz();
    13. 

  13. }
    14. 

  14. 

  15. double SigmaInelNN_Plab(double plab, double A1, double Z1, double A2, double Z2)
    16. 

  16. {
    17. 

  17. 

  18.   double sigma_inel_NN = -1.;
    19. 

  19. 

  20.   double N1 = A1-Z1;
    21. 

  21.   double N2 = A2-Z2;
    22. 

  22. 

  23.   // Data points from PDG
    24. 

  24.   std::vector<double> v_plab_tot_pp       = {0.14000    ,0.19000    ,0.24000    ,0.28000    ,0.31000    ,0.35000    ,0.37000    ,0.39000    ,0.43000    ,0.44000    ,0.49000    ,0.54000    ,0.57000    ,0.59000    ,0.60658    ,0.60700    ,0.65847    ,0.69000    ,0.72000    ,0.75000    ,0.75700    ,0.75732    ,0.83000    ,0.83086    ,0.84600    ,0.85000    ,0.87179    ,0.87200    ,0.88000    ,0.93000    ,0.93700    ,0.93736    ,0.96300    ,0.96336    ,0.96545    ,0.96824    ,0.96824    ,0.99452    ,1.00000    ,1.00549    ,1.00891    ,1.00900    ,1.00959    ,1.03673    ,1.03673    ,1.09008    ,1.09300    ,1.09337    ,1.10783    ,1.10800    ,1.11100    ,1.11100    ,1.12700    ,1.13586    ,1.14234    ,1.16168    ,1.16200    ,1.19365    ,1.20634    ,1.21500    ,1.21898    ,1.23786    ,1.24413    ,1.26909    ,1.28140    ,1.28900    ,1.29017    ,1.29388    ,1.32300    ,1.39149    ,1.40800    ,1.42000    ,1.42400    ,1.46329    ,1.49881    ,1.52235    ,1.53100    ,1.54580    ,1.59240    ,1.60000    ,1.60700    ,1.62800    ,1.66000    ,1.68700    ,1.69604    ,1.69604    ,1.73000    ,1.78000    ,1.78127    ,1.78500    ,1.85800    ,1.88600    ,1.89000    ,1.94000    ,1.95200    ,1.99000    ,2.00455    ,2.02661    ,2.05000    ,2.07900    ,2.21200    ,2.23968    ,2.28000    ,2.41900    ,2.45000    ,2.47000    ,2.59200    ,2.68000    ,2.70400    ,2.78444    ,2.81900    ,2.85700    ,2.88976    ,2.95800    ,2.97000    ,2.99400    ,3.05400    ,3.11000    ,3.13100    ,3.14200    ,3.27000    ,3.27700    ,3.30300    ,3.41160    ,3.44400    ,3.54600    ,3.58000    ,3.67024    ,3.73100    ,3.90800    ,4.00000    ,4.00000    ,4.03700    ,4.15000    ,4.20000    ,4.26500    ,4.50990    ,4.55200    ,4.78300    ,4.96600    ,5.22100    ,5.52000    ,5.52600    ,5.82400    ,5.83020    ,5.96493    ,6.00000    ,7.00000    ,7.75020    ,7.83500    ,7.85000    ,8.00000    ,8.10000    ,9.90020    ,10.00000   ,10.00000   ,10.00000   ,10.01000   ,10.11000   ,12.00000   ,12.00000   ,12.40000   ,14.00000   ,15.00000   ,15.80000   ,16.00000   ,17.70000   ,18.00000   ,19.33000   ,19.40000   ,19.80000   ,20.00000   ,20.00000   ,20.00000   ,21.40000   ,22.00000   ,23.00000   ,23.50000   ,24.00000   ,24.20000   ,24.50000   ,25.00000   ,26.42000   ,28.40000   ,30.00000   ,30.00000   ,32.10000   ,35.00000   ,35.00000   ,40.00000   ,40.00000   ,42.50000   ,45.00000   ,50.00000   ,50.00000   ,50.00000   ,50.00000   ,50.00000   ,52.20000   ,55.00000   ,60.00000   ,65.00000   ,69.00000   ,70.00000   ,72.99398   ,97.99551   ,100.00000  ,100.00000  ,100.00000  ,102.00000  ,120.00000  ,121.99639  ,146.99700  ,147.00000  ,150.00000  ,150.00000  ,170.00000  ,171.99743  ,175.00000  ,195.99775  ,200.00000  ,200.00000  ,205.00000  ,240.00000  ,280.00000  ,293.35186  ,293.35186  ,293.35186  ,293.35186  ,293.35186  ,295.86181  ,300.00000  ,300.00000  ,303.00000  ,310.00000  ,340.00000  ,370.00000  ,405.00000  ,498.04302  ,498.04302  ,498.04302  ,498.04302  ,498.04302  ,498.04998  ,501.30972  ,504.58697  ,659.34014  ,1063.80004 ,1063.83483 ,1073.38439 ,1087.78857 ,1484.68920 ,1484.68920 ,1484.68920 ,1484.68920 ,1484.68920 ,1490.30004 ,1507.28405 ,2067.38403 ,2074.10009 ,2080.68530 ,2094.02905 ,2094.02905 ,2100.71679 ,2114.12451 ,4.0594E+06 ,3982.00000 ,1.0001E+04 ,1.9954E+04 ,3.9812E+04 ,7.9433E+04 ,2.0248E+07 ,2.6112E+07 ,2.6112E+07 ,2.6112E+07 ,2.6112E+07 ,3.4104E+07 ,3.4104E+07 ,3.4104E+07 ,3.5186E+07 ,6.1147E+07 ,9.0046E+07 ,9.0046E+07 ,1.0626E+08 ,1.8466E+08 ,3.2090E+08 ,4.7961E+08 ,1.7314E+09 ,4.8092E+09};
    25. 

  25.   std::vector<double> v_plab_el_pp        = {0.84600   ,0.97935   ,1.00000   ,1.00959   ,1.03673   ,1.11100   ,1.12700   ,1.16810   ,1.20634   ,1.21500   ,1.23200   ,1.28140   ,1.29017   ,1.29400   ,1.32300   ,1.42400   ,1.46329   ,1.47516   ,1.53100   ,1.62800   ,1.66167   ,1.66167   ,1.68700   ,1.69000   ,1.69604   ,1.78500   ,1.88600   ,1.99000   ,2.20000   ,2.22883   ,2.80000   ,2.80553   ,3.00000   ,3.03662   ,3.65000   ,3.67024   ,4.00000   ,4.15000   ,4.80000   ,5.00000   ,5.25517   ,5.52000   ,5.96493   ,6.00000   ,6.60000   ,6.79990   ,6.80000   ,6.92000   ,7.02590   ,7.07634   ,7.07634   ,8.10000   ,8.50000   ,8.80000   ,8.80000   ,8.90000   ,9.00000   ,9.11008   ,9.39140   ,9.89388   ,9.90000   ,10.01000  ,10.80000  ,10.90000  ,12.00000  ,12.00000  ,12.10000  ,12.80000  ,13.20000  ,14.80000  ,14.90877  ,15.00000  ,15.10000  ,15.50000  ,15.50000  ,16.20000  ,16.70000  ,18.60000  ,18.90000  ,19.00000  ,19.20000  ,19.60000  ,20.00000  ,21.40000  ,21.70000  ,23.50000  ,24.00000  ,24.50000  ,24.60000  ,26.20000  ,27.50000  ,29.70000  ,30.00000  ,30.50000  ,33.30000  ,36.20000  ,38.00000  ,40.60000  ,45.20000  ,50.00000  ,50.00000  ,50.00000  ,50.60000  ,52.10000  ,54.40000  ,57.00000  ,58.00000  ,60.00000  ,60.20000  ,63.50000  ,66.10000  ,69.20000  ,69.80001  ,70.00000  ,70.00000  ,100.00000 ,100.00000 ,102.00000 ,140.00000 ,147.00000 ,150.00000 ,175.00000 ,200.00000 ,200.00000 ,205.00000 ,205.00000 ,250.00000 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,300.00000 ,300.00000 ,303.00000 ,491.54168 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,1073.38439,1073.38439,1484.68920,1484.68920,1484.68920,1484.68920,2067.38403,2080.68530,4.0594E+06,2.6112E+07,2.6112E+07,3.4104E+07,3.4104E+07,9.0046E+07,9.0046E+07};
    26. 

  26.   std::vector<double> v_plab_tot_np       = {0.00385  ,0.03043  ,0.03873  ,0.04347  ,0.04502  ,0.04973  ,0.05000  ,0.05020  ,0.05100  ,0.05200  ,0.05300  ,0.05311  ,0.05400  ,0.05448  ,0.05500  ,0.05600  ,0.05700  ,0.05800  ,0.05900  ,0.05970  ,0.06000  ,0.06100  ,0.06200  ,0.06300  ,0.06400  ,0.06500  ,0.06571  ,0.06600  ,0.06700  ,0.06800  ,0.06900  ,0.06907  ,0.06913  ,0.07000  ,0.07100  ,0.07121  ,0.07200  ,0.07300  ,0.07400  ,0.07500  ,0.07600  ,0.07633  ,0.07700  ,0.07767  ,0.07800  ,0.07900  ,0.08000  ,0.08100  ,0.08112  ,0.08200  ,0.08300  ,0.08400  ,0.08500  ,0.08570  ,0.08600  ,0.08700  ,0.08787  ,0.08800  ,0.08900  ,0.08999  ,0.09000  ,0.09100  ,0.09200  ,0.09300  ,0.09400  ,0.09410  ,0.09459  ,0.09500  ,0.09600  ,0.09700  ,0.09798  ,0.09800  ,0.09900  ,0.10000  ,0.10100  ,0.10176  ,0.10200  ,0.10300  ,0.10400  ,0.10500  ,0.10541  ,0.10600  ,0.10700  ,0.10800  ,0.10900  ,0.11000  ,0.11100  ,0.11200  ,0.11245  ,0.11300  ,0.11400  ,0.11500  ,0.11564  ,0.11600  ,0.11630  ,0.11700  ,0.11800  ,0.11900  ,0.12000  ,0.12100  ,0.12200  ,0.12287  ,0.12300  ,0.12400  ,0.12500  ,0.12600  ,0.12700  ,0.12800  ,0.12867  ,0.12900  ,0.13000  ,0.13047  ,0.13100  ,0.13200  ,0.13400  ,0.13500  ,0.13600  ,0.13700  ,0.13745  ,0.13800  ,0.13900  ,0.14000  ,0.14032  ,0.14100  ,0.14200  ,0.14300  ,0.14400  ,0.14419  ,0.14500  ,0.14505  ,0.14600  ,0.14700  ,0.14800  ,0.14900  ,0.15000  ,0.15064  ,0.15100  ,0.15200  ,0.15300  ,0.15377  ,0.15400  ,0.15500  ,0.15600  ,0.15681  ,0.15700  ,0.15762  ,0.15800  ,0.15900  ,0.15985  ,0.16000  ,0.16100  ,0.16200  ,0.16280  ,0.16280  ,0.16292  ,0.16300  ,0.16339  ,0.16339  ,0.16400  ,0.16500  ,0.16600  ,0.16700  ,0.16800  ,0.16853  ,0.16856  ,0.16900  ,0.17000  ,0.17100  ,0.17137  ,0.17200  ,0.17300  ,0.17300  ,0.17400  ,0.17413  ,0.17500  ,0.17500  ,0.17600  ,0.17686  ,0.17700  ,0.17700  ,0.17800  ,0.17800  ,0.17900  ,0.17954  ,0.18000  ,0.18000  ,0.18100  ,0.18200  ,0.18200  ,0.18218  ,0.18300  ,0.18375  ,0.18400  ,0.18400  ,0.18479  ,0.18500  ,0.18600  ,0.18600  ,0.18700  ,0.18800  ,0.18800  ,0.18900  ,0.19000  ,0.19000  ,0.19100  ,0.19200  ,0.19274  ,0.19281  ,0.19300  ,0.19300  ,0.19324  ,0.19324  ,0.19400  ,0.19455  ,0.19500  ,0.19500  ,0.19519  ,0.19600  ,0.19700  ,0.19700  ,0.19782  ,0.19800  ,0.19884  ,0.19900  ,0.19900  ,0.20000  ,0.20100  ,0.20200  ,0.20200  ,0.20243  ,0.20300  ,0.20400  ,0.20400  ,0.20418  ,0.20500  ,0.20600  ,0.20700  ,0.20700  ,0.20800  ,0.20900  ,0.20900  ,0.20923  ,0.21000  ,0.21100  ,0.21100  ,0.21200  ,0.21300  ,0.21350  ,0.21400  ,0.21400  ,0.21500  ,0.21600  ,0.21609  ,0.21654  ,0.21700  ,0.21700  ,0.21800  ,0.21900  ,0.21950  ,0.22000  ,0.22000  ,0.22100  ,0.22200  ,0.22200  ,0.22263  ,0.22300  ,0.22395  ,0.22400  ,0.22500  ,0.22500  ,0.22600  ,0.22700  ,0.22800  ,0.22800  ,0.22900  ,0.22905  ,0.23000  ,0.23087  ,0.23100  ,0.23100  ,0.23200  ,0.23234  ,0.23300  ,0.23400  ,0.23400  ,0.23500  ,0.23524  ,0.23600  ,0.23626  ,0.23700  ,0.23800  ,0.23800  ,0.23900  ,0.24000  ,0.24100  ,0.24100  ,0.24118  ,0.24200  ,0.24300  ,0.24400  ,0.24400  ,0.24500  ,0.24600  ,0.24700  ,0.24800  ,0.24800  ,0.24900  ,0.25000  ,0.25100  ,0.25100  ,0.25200  ,0.25500  ,0.25900  ,0.26300  ,0.26383  ,0.26700  ,0.26991  ,0.27200  ,0.27351  ,0.27551  ,0.27600  ,0.27707  ,0.28100  ,0.28406  ,0.28406  ,0.28600  ,0.29100  ,0.29426  ,0.29600  ,0.29818  ,0.30088  ,0.30100  ,0.30414  ,0.30700  ,0.30757  ,0.31200  ,0.31375  ,0.31800  ,0.32309  ,0.32423  ,0.32500  ,0.33100  ,0.33220  ,0.33800  ,0.33919  ,0.34110  ,0.34400  ,0.34979  ,0.34979  ,0.35122  ,0.35830  ,0.35858  ,0.36600  ,0.36664  ,0.37481  ,0.37500  ,0.38283  ,0.38400  ,0.39071  ,0.39300  ,0.39846  ,0.40300  ,0.40608  ,0.41300  ,0.41359  ,0.41606  ,0.41700  ,0.42098  ,0.42098  ,0.42098  ,0.42400  ,0.42826  ,0.42923  ,0.43306  ,0.43306  ,0.43500  ,0.43545  ,0.43782  ,0.43829  ,0.44000  ,0.44253  ,0.44700  ,0.44744  ,0.44953  ,0.45644  ,0.45900  ,0.46054  ,0.46326  ,0.46800  ,0.47001  ,0.47100  ,0.47667  ,0.48327  ,0.48327  ,0.48900  ,0.48979  ,0.49625  ,0.50264  ,0.50264  ,0.50897  ,0.51000  ,0.51524  ,0.52145  ,0.52760  ,0.53168  ,0.53300  ,0.53370  ,0.53975  ,0.54575  ,0.55170  ,0.55300  ,0.55760  ,0.55760  ,0.56346  ,0.56346  ,0.56927  ,0.57119  ,0.57503  ,0.58076  ,0.58645  ,0.58833  ,0.59209  ,0.59770  ,0.60327  ,0.60881  ,0.60881  ,0.61431  ,0.61977  ,0.62520  ,0.63060  ,0.63597  ,0.64131  ,0.64485  ,0.64501  ,0.66236  ,0.66582  ,0.67957  ,0.67957  ,0.67957  ,0.69649  ,0.71316  ,0.72500  ,0.72958  ,0.74577  ,0.75857  ,0.76175  ,0.76175  ,0.76175  ,0.77500  ,0.77753  ,0.79313  ,0.80854  ,0.82379  ,0.82500  ,0.83738  ,0.83888  ,0.85382  ,0.86862  ,0.87500  ,0.88329  ,0.89783  ,0.91080  ,0.91224  ,0.92500  ,0.92654  ,0.94073  ,0.95481  ,0.96879  ,0.96879  ,0.97500  ,0.97851  ,0.98267  ,0.99646  ,1.01016  ,1.02377  ,1.02400  ,1.03731  ,1.05076  ,1.06414  ,1.07500  ,1.07744  ,1.08407  ,1.09067  ,1.10384  ,1.11694  ,1.12500  ,1.12998  ,1.14295  ,1.15587  ,1.16873  ,1.17500  ,1.18153  ,1.19428  ,1.20698  ,1.21962  ,1.22500  ,1.23222  ,1.24478  ,1.25728  ,1.25728  ,1.26974  ,1.27500  ,1.28216  ,1.29454  ,1.30688  ,1.31917  ,1.32500  ,1.33143  ,1.34365  ,1.35584  ,1.36798  ,1.37500  ,1.38010  ,1.39218  ,1.40423  ,1.41624  ,1.42500  ,1.42823  ,1.50000  ,1.60000  ,1.70000  ,1.80000  ,1.90000  ,2.02200  ,2.14261  ,2.19300  ,2.39100  ,2.62700  ,2.94400  ,3.00000  ,3.27700  ,3.41260  ,4.00000  ,4.30000  ,4.74760  ,5.70000  ,5.86478  ,6.37060  ,6.37066  ,6.50000  ,7.78310  ,8.00000  ,8.24680  ,9.19190  ,10.00000 ,11.00000 ,12.00000 ,13.00000 ,14.00000 ,14.60000 ,16.00000 ,17.00000 ,17.80000 ,18.00000 ,21.00000 ,21.60000 ,26.28321 ,26.50000 ,27.00000 ,34.00000 ,34.00000 ,36.58794 ,52.69162 ,80.00000 ,80.00000 ,131.00000,131.00000,133.00000,175.00000,180.00000,180.00000,215.00000,215.00000,240.00000,240.00000,273.00000,273.00000,310.00000};
    27. 

  27.   std::vector<double> v_plab_el_np        = {1.25000 ,1.73000 ,2.23000 ,3.00000 ,3.00000 ,3.83000 ,5.10000 ,5.50000 ,8.00000 ,17.00000};
    28. 

  28.   std::vector<double> v_sigm_tot_pp       = {314.00,155.00,92.000,70.000,52.800,42.500,37.400,33.900,28.500,27.700,24.800,25.200,26.100,23.200,24.200,24.550,25.800,22.400,22.400,22.600,23.850,23.850,24.300,24.300,24.400,23.400,24.550,24.450,23.200,24.400,25.650,25.700,26.300,26.500,24.000,26.500,26.900,24.200,28.000,23.800,27.950,27.750,27.000,28.000,27.600,29.900,30.950,31.250,31.950,31.650,32.900,34.029,30.000,29.800,32.100,34.800,34.500,35.600,36.000,39.000,36.600,37.700,38.600,39.800,41.800,43.234,41.000,41.400,43.000,44.400,46.487,46.200,46.000,47.000,47.800,47.600,47.400,47.940,46.100,47.500,47.476,49.000,47.553,48.000,47.400,48.000,46.200,47.490,48.300,47.000,47.455,47.400,46.800,47.357,47.409,47.000,47.500,49.400,45.300,47.224,46.985,47.200,46.669,46.130,45.827,45.100,45.533,45.331,45.174,41.400,45.008,44.928,45.540,44.651,44.500,44.466,44.401,44.188,44.156,44.114,47.100,43.610,43.669,41.600,43.138,42.978,43.200,42.100,42.680,42.316,43.000,41.600,42.136,39.700,42.200,41.765,42.100,41.457,41.377,41.165,41.171,41.600,40.878,40.848,41.600,38.300,40.600,40.600,41.600,40.075,40.000,40.000,40.100,39.400,39.900,39.600,40.200,41.100,40.000,39.400,39.600,39.000,39.100,39.290,38.700,38.700,39.700,38.700,38.900,39.700,38.600,38.400,39.060,38.850,39.400,38.300,39.390,39.700,38.900,38.700,39.300,38.800,38.800,39.900,38.590,38.550,38.400,38.490,38.460,38.500,38.850,38.420,38.450,38.460,38.200,38.140,38.420,37.740,37.870,38.430,38.440,38.650,36.680,38.280,37.600,37.600,38.700,38.460,38.390,38.900,38.580,37.000,37.600,38.470,38.690,38.620,38.830,37.400,39.600,38.500,38.980,38.900,39.000,39.240,39.420,39.650,39.400,39.130,38.800,38.900,38.700,39.210,40.680,39.000,39.590,39.690,39.770,40.600,40.220,40.110,39.910,40.070,40.200,40.100,40.100,40.000,40.420,41.700,41.900,41.890,42.500,43.010,42.380,42.100,42.850,42.710,42.400,42.900,43.550,43.100,44.000,43.820,42.960,44.100,42.200,84.700,43.000,48.000,54.000,56.000,54.000,93.000,98.300,98.580,98.000,95.350,101.70, 96.07,102.90,101.00,117.00,110.6 ,109.5 ,104.00,100.00,124.00,120.00,133.00,170.00};
    29. 

  29.   std::vector<double> v_sigm_el_pp        = {23.200,24.000,25.800,24.000,22.000,25.800,24.000,25.200,25.000,25.300,24.200,25.100,23.000,24.700,24.600,24.300,21.000,24.000,24.100,26.000,26.800,24.800,24.700,28.200,19.000,22.700,22.500,22.100,22.200,19.860,16.300,19.210,17.200,17.000,15.200,15.320,13.500,11.600,14.400,12.700,10.000,11.990,10.000,11.500,11.470,11.790,10.600,11.400,8.0000,9.7000,9.8000,10.800,8.7400,11.710,9.8000,10.100,10.840,10.800,8.6000,10.000,10.200,10.200,11.040,9.9000,9.8700,9.8500,10.400,10.890,8.8700,10.480,11.000,8.1300,9.7000,9.2000,8.7500,9.3600,9.7400,10.200,8.5900,8.7000,9.4000,9.6400,9.0000,8.0000,8.1500,8.3000,8.3000,8.8000,8.0200,9.8000,7.9600,8.2000,7.7000,7.8700,7.6600,7.7000,7.6000,7.5200,7.4000,7.6100,7.0000,7.5600,7.4800,7.3300,7.2300,7.2100,7.4900,6.6000,7.2500,6.8900,7.0700,6.8600,6.8600,7.4100,7.1000,7.0700,7.0800,6.9000,7.0000,6.8500,6.9700,7.0600,6.8700,6.9500,6.9200,6.9200,6.7800,6.7000,6.8200,6.8000,6.8100,7.2900,7.8900,7.2000,6.8000,6.9000,7.3900,7.0000,6.7500,7.4500,7.5000,7.7900,7.5600,7.6000,7.1700,7.5100,7.7700,21.800,25.430,24.000,27.100,24.330,31.0  ,30.7};
    30. 

  30.   std::vector<double> v_sigm_tot_np       = {20360.,6202.0,4700.0,4228.0,4060.0,3675.0,3609.5,3630.0,3542.9,3472.1,3398.6,3447.0,3326.5,3330.0,3266.6,3206.6,3138.6,3079.6,3015.5,3004.0,2963.6,2907.2,2851.3,2792.6,2754.0,2705.5,2677.0,2651.5,2607.4,2558.2,2507.0,2536.0,2525.0,2464.0,2433.5,2417.0,2380.4,2341.5,2295.9,2258.5,2225.0,2274.0,2185.9,2206.0,2151.0,2123.9,2091.0,2049.9,2067.0,2017.0,1988.9,1945.4,1920.0,1948.0,1893.5,1864.0,1730.0,1845.0,1812.5,1829.0,1782.0,1744.4,1727.0,1701.8,1675.8,1679.0,1690.0,1653.5,1626.2,1601.1,1620.0,1577.3,1577.2,1533.1,1516.1,1517.0,1489.4,1468.5,1447.1,1427.2,1454.0,1402.0,1388.6,1367.7,1346.4,1329.8,1308.6,1290.3,1327.0,1274.1,1255.6,1238.1,1263.0,1218.0,1240.0,1203.0,1188.0,1170.8,1155.3,1138.8,1127.5,1134.0,1109.7,1099.6,1083.1,1066.9,1052.4,1038.0,1040.0,1024.3,1013.8,1026.0,1002.6,990.21,965.19,948.24,938.60,927.03,945.00,913.80,902.60,892.87,940.00,883.72,874.51,863.56,853.04,871.00,844.37,880.00,832.26,820.73,811.12,800.64,793.00,809.00,782.71,774.13,766.97,690.00,755.83,745.52,736.17,749.00,727.95,760.00,719.31,712.13,694.00,703.26,695.43,686.08,690.00,686.00,720.00,678.44,690.00,689.00,670.72,662.77,656.77,650.32,642.11,640.00,662.00,635.37,628.45,618.76,648.00,615.16,537.00,608.56,601.03,582.00,594.86,634.00,588.04,609.00,581.75,580.00,576.73,597.00,569.24,612.00,563.15,561.00,558.70,544.00,552.03,587.00,547.79,560.00,540.36,533.00,539.00,534.22,517.00,526.15,522.88,523.00,518.45,513.93,498.00,507.95,502.90,497.93,493.30,498.00,492.44,522.00,494.20,495.00,487.73,504.00,479.00,484.43,479.00,478.55,479.00,476.00,470.00,470.06,467.00,449.00,464.20,459.50,455.90,447.00,450.55,448.00,446.96,420.00,439.34,432.00,436.40,430.35,408.00,426.59,425.55,433.00,421.40,411.00,418.45,397.70,414.63,408.09,403.45,397.20,401.20,393.00,395.10,393.40,388.00,385.50,393.50,389.59,384.68,382.49,388.00,362.90,379.59,373.55,377.60,369.60,367.00,369.64,329.00,365.89,345.70,363.53,358.20,354.29,335.40,351.40,348.00,347.00,345.91,337.60,321.50,341.61,333.50,340.00,334.30,312.20,330.70,328.90,331.00,329.40,315.80,323.45,309.10,320.45,316.10,311.24,281.70,313.00,287.00,306.52,311.60,286.40,304.10,297.20,302.50,295.94,288.70,294.50,290.19,289.30,276.00,287.00,285.50,260.10,245.80,220.00,247.50,226.80,223.00,223.50,227.23,223.70,219.10,170.00,204.60,203.00,203.25,196.00,189.10,190.31,187.30,185.40,196.00,170.30,173.38,166.80,167.80,152.00,163.10,152.00,150.19,151.50,144.80,141.50,142.64,135.80,133.30,134.25,124.30,125.21,126.00,116.80,118.86,114.90,109.60,112.51,107.77,100.10,102.88,100.30,97.320,94.700,93.544,85.300,90.375,82.100,86.193,84.500,85.700,83.000,76.000,82.783,77.700,80.316,77.000,73.900,73.000,72.100,77.676,74.000,76.000,76.600,75.247,71.500,80.000,73.636,70.841,69.400,59.000,68.958,68.300,67.234,68.300,65.355,63.200,61.500,62.800,61.820,60.388,56.900,59.424,58.410,58.000,57.021,56.165,54.876,48.500,54.500,53.991,53.367,51.768,51.017,51.700,50.768,46.400,50.425,50.500,49.544,51.200,48.831,47.936,47.638,49.200,46.318,45.961,45.272,44.000,45.050,44.501,43.976,43.446,43.254,42.870,42.022,41.918,42.700,41.040,42.460,41.300,41.000,39.966,38.963,38.085,42.800,37.144,36.819,38.130,38.000,33.000,36.355,38.630,35.582,35.106,34.563,34.560,37.320,36.140,34.100,33.938,33.487,33.970,33.749,33.163,35.050,32.680,34.090,32.685,32.771,32.417,33.700,32.408,33.900,34.530,32.548,32.914,32.410,32.634,33.570,32.714,32.952,32.813,34.160,33.154,35.150,33.308,33.139,33.231,34.940,33.055,33.623,34.021,33.732,35.470,34.662,34.354,34.609,35.318,36.040,34.665,35.347,35.875,35.200,35.745,36.850,36.117,36.486,36.953,37.415,37.600,37.190,37.357,37.848,38.197,37.850,38.493,38.263,38.241,38.656,37.980,38.318,37.820,38.120,38.210,38.850,39.350,39.990,42.400,40.500,40.670,40.750,40.760,43.000,40.830,38.100,43.100,40.400,43.400,42.500,33.600,41.200,41.200,38.700,39.300,39.700,41.200,40.800,39.500,39.400,39.300,39.000,38.700,37.100,39.100,38.500,37.500,38.700,38.500,37.700,38.940,39.300,38.900,38.200,38.200,38.850,38.200,38.980,38.980,39.170,39.140,39.210,39.210,39.520,39.620,39.790,39.790,39.660,39.660,40.320,40.320,39.600};
    31. 

  31.   std::vector<double> v_sigm_el_np        = {28.580    ,24.300    ,21.140    ,14.200    ,1.3000    ,15.680    ,12.920    ,12.500    ,6.0000E-02,10.100    };
    32. 

  32.   std::vector<double> v_plabMin_tot_pp    = {0.14000   ,0.19000   ,0.24000   ,0.28000   ,0.31000   ,0.35000   ,0.37000   ,0.39000   ,0.43000   ,0.44000   ,0.49000   ,0.54000   ,0.57000   ,0.59000   ,0.60289   ,0.58700   ,0.65148   ,0.69000   ,0.72000   ,0.75000   ,0.73700   ,0.75413   ,0.83000   ,0.81875   ,0.83500   ,0.85000   ,0.86885   ,0.85200   ,0.88000   ,0.93000   ,0.91700   ,0.93453   ,0.94300   ,0.96056   ,0.95146   ,0.95678   ,0.96824   ,0.99452   ,0.98900   ,1.00549   ,1.00618   ,0.98900   ,0.99452   ,1.03552   ,1.03673   ,1.09008   ,1.07300   ,1.09074   ,1.10521   ,1.08800   ,1.10000   ,1.10878   ,1.12700   ,1.12194   ,1.14234   ,1.15911   ,1.14200   ,1.19365   ,1.18760   ,1.20500   ,1.21898   ,1.22239   ,1.24413   ,1.26909   ,1.28140   ,1.28642   ,1.28770   ,1.29388   ,1.31400   ,1.37360   ,1.40518   ,1.38450   ,1.41500   ,1.43948   ,1.47919   ,1.50235   ,1.52000   ,1.50494   ,1.59240   ,1.56000   ,1.60379   ,1.62800   ,1.65668   ,1.67500   ,1.69604   ,1.66741   ,1.68675   ,1.77644   ,1.75693   ,1.77300   ,1.85428   ,1.87400   ,1.84275   ,1.93612   ,1.94810   ,1.97900   ,1.97636   ,1.99803   ,1.99875   ,2.07484   ,2.20758   ,2.20733   ,2.27544   ,2.41416   ,2.44510   ,2.40825   ,2.58682   ,2.67464   ,2.69859   ,2.74219   ,2.81336   ,2.85129   ,2.80131   ,2.95208   ,2.89575   ,2.98801   ,3.04789   ,3.10378   ,3.12474   ,3.13572   ,3.27000   ,3.27045   ,3.29639   ,3.35764   ,3.43711   ,3.53891   ,3.49050   ,3.67024   ,3.72354   ,3.90018   ,3.96000   ,3.90000   ,4.02893   ,4.10000   ,4.15800   ,4.25647   ,4.50990   ,4.54290   ,4.77343   ,4.95607   ,5.21056   ,5.51000   ,5.51495   ,5.81235   ,5.83020   ,5.96493   ,5.89800   ,6.93000   ,7.75020   ,7.81933   ,7.83200   ,7.86400   ,8.00000   ,9.39990   ,9.83000   ,9.90000   ,10.00000  ,10.00000  ,10.11000  ,11.79600  ,12.00000  ,11.80000  ,13.76200  ,15.00000  ,15.05000  ,15.72800  ,16.80000  ,17.69400  ,19.33000  ,18.40000  ,19.80000  ,19.66000  ,20.00000  ,20.00000  ,20.30000  ,21.62600  ,23.00000  ,23.50000  ,24.00000  ,22.80000  ,24.50000  ,25.00000  ,26.42000  ,27.00000  ,30.00000  ,30.00000  ,31.90000  ,35.00000  ,35.00000  ,40.00000  ,40.00000  ,42.50000  ,45.00000  ,50.00000  ,50.00000  ,50.00000  ,50.00000  ,50.00000  ,52.20000  ,55.00000  ,60.00000  ,65.00000  ,69.00000  ,70.00000  ,72.99398  ,97.99551  ,99.00000  ,100.00000 ,100.00000 ,102.00000 ,120.00000 ,121.99639 ,146.99700 ,147.00000 ,150.00000 ,150.00000 ,170.00000 ,171.99743 ,173.25000 ,195.99775 ,200.00000 ,200.00000 ,205.00000 ,240.00000 ,280.00000 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,295.86181 ,300.00000 ,298.50000 ,303.00000 ,310.00000 ,340.00000 ,370.00000 ,405.00000 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,498.04998 ,501.30972 ,504.58697 ,659.34014 ,1063.80004,1063.83483,1073.38439,1087.78857,1484.68920,1484.68920,1484.68920,1484.68920,1484.68920,1490.30004,1507.28405,2067.38403,2074.10009,2080.68530,2094.02905,2094.02905,2100.71679,2114.12451,4.0594E+06,1996.19995,6310.50000,1.2590E+04,2.5120E+04,3.9812E+04,1.4791E+07,2.6112E+07,2.6112E+07,2.6112E+07,2.6112E+07,3.4104E+07,3.4104E+07,3.4104E+07,2.5704E+07,4.4668E+07,9.0046E+07,9.0046E+07,7.7625E+07,1.3490E+08,2.3442E+08,4.7961E+08,1.3861E+09,4.0333e+09};
    33. 

  33.   std::vector<double> v_plabMin_el_pp     = {0.83500   ,0.95846   ,0.98900   ,0.99452   ,1.03552   ,1.10000   ,1.12700   ,1.16810   ,1.18760   ,1.20500   ,1.23200   ,1.28140   ,1.28770   ,1.29400   ,1.31400   ,1.41500   ,1.43948   ,1.35515   ,1.52000   ,1.62800   ,1.66167   ,1.64442   ,1.67500   ,1.69000   ,1.66741   ,1.77300   ,1.87400   ,1.97900   ,2.20000   ,2.22883   ,2.80000   ,2.80553   ,3.00000   ,3.01317   ,3.65000   ,3.64081   ,3.96000   ,4.10000   ,4.80000   ,5.00000   ,5.21046   ,5.51000   ,5.96493   ,6.00000   ,6.59010   ,6.79990   ,6.80000   ,6.84500   ,6.96385   ,7.07634   ,7.07634   ,8.00000   ,8.50000   ,8.80000   ,8.80000   ,8.90000   ,9.00000   ,9.10003   ,9.39140   ,9.89388   ,9.80100   ,10.00000  ,10.80000  ,10.90000  ,12.00000  ,12.00000  ,12.10000  ,12.80000  ,13.20000  ,14.80000  ,14.90877  ,15.00000  ,14.94900  ,15.50000  ,15.50000  ,16.20000  ,16.70000  ,18.60000  ,18.90000  ,19.00000  ,19.20000  ,19.60000  ,19.80000  ,21.40000  ,21.70000  ,23.50000  ,24.00000  ,24.50000  ,24.60000  ,26.20000  ,27.50000  ,29.40300  ,30.00000  ,30.50000  ,33.30000  ,36.20000  ,38.00000  ,40.60000  ,45.20000  ,49.70000  ,50.00000  ,50.00000  ,50.60000  ,52.10000  ,54.40000  ,57.00000  ,58.00000  ,60.00000  ,60.20000  ,63.50000  ,66.10000  ,69.20000  ,69.80001  ,69.58000  ,70.00000  ,99.40000  ,100.00000 ,102.00000 ,139.16000 ,147.00000 ,150.00000 ,173.94999 ,200.00000 ,200.00000 ,205.00000 ,205.00000 ,250.00000 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,300.00000 ,298.50000 ,303.00000 ,491.54168 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,1073.38439,1073.38439,1484.68920,1484.68920,1484.68920,1484.68920,2067.38403,2080.68530,4.0594E+06,2.6112E+07,2.6112E+07,3.4104E+07,3.4104E+07,9.0046E+07,9.0046E+07};
    34. 

  34.   std::vector<double> v_plabMin_tot_np    = {0.00123  ,0.03043  ,0.03873  ,0.04347  ,0.04502  ,0.04973  ,0.05000  ,0.05020  ,0.05100  ,0.05200  ,0.05300  ,0.05311  ,0.05400  ,0.05448  ,0.05500  ,0.05600  ,0.05700  ,0.05800  ,0.05900  ,0.05970  ,0.06000  ,0.06100  ,0.06200  ,0.06300  ,0.06400  ,0.06500  ,0.06571  ,0.06600  ,0.06700  ,0.06800  ,0.06900  ,0.06907  ,0.06913  ,0.07000  ,0.07100  ,0.07121  ,0.07200  ,0.07300  ,0.07400  ,0.07500  ,0.07600  ,0.07633  ,0.07700  ,0.07767  ,0.07800  ,0.07900  ,0.08000  ,0.08100  ,0.08112  ,0.08200  ,0.08300  ,0.08400  ,0.08500  ,0.08570  ,0.08600  ,0.08700  ,0.08787  ,0.08800  ,0.08900  ,0.08999  ,0.09000  ,0.09100  ,0.09200  ,0.09300  ,0.09400  ,0.09410  ,0.09450  ,0.09500  ,0.09600  ,0.09700  ,0.09798  ,0.09800  ,0.09900  ,0.10000  ,0.10100  ,0.10176  ,0.10200  ,0.10300  ,0.10400  ,0.10500  ,0.10541  ,0.10600  ,0.10700  ,0.10800  ,0.10900  ,0.11000  ,0.11100  ,0.11200  ,0.11245  ,0.11300  ,0.11400  ,0.11500  ,0.11564  ,0.11600  ,0.11630  ,0.11700  ,0.11800  ,0.11900  ,0.12000  ,0.12100  ,0.12200  ,0.12287  ,0.12300  ,0.12400  ,0.12500  ,0.12600  ,0.12700  ,0.12800  ,0.12867  ,0.12900  ,0.13000  ,0.13047  ,0.13100  ,0.13200  ,0.13400  ,0.13500  ,0.13600  ,0.13700  ,0.13745  ,0.13800  ,0.13900  ,0.14000  ,0.14032  ,0.14100  ,0.14200  ,0.14300  ,0.14400  ,0.14419  ,0.14500  ,0.14505  ,0.14600  ,0.14700  ,0.14800  ,0.14900  ,0.15000  ,0.15064  ,0.15100  ,0.15200  ,0.15300  ,0.15377  ,0.15400  ,0.15500  ,0.15600  ,0.15681  ,0.15700  ,0.15762  ,0.15800  ,0.15900  ,0.15985  ,0.16000  ,0.16100  ,0.16200  ,0.16280  ,0.16280  ,0.16292  ,0.16300  ,0.16339  ,0.16044  ,0.16400  ,0.16500  ,0.16600  ,0.16700  ,0.16800  ,0.16853  ,0.16856  ,0.16900  ,0.17000  ,0.17100  ,0.17137  ,0.17200  ,0.17300  ,0.17300  ,0.17400  ,0.17413  ,0.17500  ,0.17500  ,0.17600  ,0.17686  ,0.17700  ,0.17700  ,0.17800  ,0.17800  ,0.17900  ,0.17954  ,0.18000  ,0.18000  ,0.18100  ,0.18200  ,0.18200  ,0.18218  ,0.18300  ,0.18375  ,0.18400  ,0.18400  ,0.18479  ,0.18500  ,0.18600  ,0.18600  ,0.18700  ,0.18800  ,0.18800  ,0.18900  ,0.19000  ,0.19000  ,0.19100  ,0.19200  ,0.19261  ,0.19281  ,0.19300  ,0.19300  ,0.19311  ,0.19324  ,0.19400  ,0.19390  ,0.19500  ,0.19500  ,0.19519  ,0.19600  ,0.19700  ,0.19700  ,0.19782  ,0.19800  ,0.19884  ,0.19900  ,0.19900  ,0.20000  ,0.20100  ,0.20200  ,0.20200  ,0.20243  ,0.20300  ,0.20400  ,0.20400  ,0.20418  ,0.20500  ,0.20600  ,0.20700  ,0.20700  ,0.20800  ,0.20900  ,0.20900  ,0.20923  ,0.21000  ,0.21100  ,0.21100  ,0.21200  ,0.21300  ,0.21342  ,0.21400  ,0.21400  ,0.21500  ,0.21600  ,0.21609  ,0.21654  ,0.21700  ,0.21700  ,0.21800  ,0.21900  ,0.21950  ,0.22000  ,0.22000  ,0.22100  ,0.22200  ,0.22200  ,0.22263  ,0.22300  ,0.22395  ,0.22400  ,0.22500  ,0.22500  ,0.22600  ,0.22700  ,0.22800  ,0.22800  ,0.22900  ,0.22905  ,0.23000  ,0.23081  ,0.23100  ,0.23100  ,0.23200  ,0.23234  ,0.23300  ,0.23400  ,0.23400  ,0.23500  ,0.23524  ,0.23600  ,0.23626  ,0.23700  ,0.23800  ,0.23800  ,0.23900  ,0.24000  ,0.24100  ,0.24100  ,0.24118  ,0.24200  ,0.24300  ,0.24400  ,0.24400  ,0.24500  ,0.24600  ,0.24700  ,0.24800  ,0.24800  ,0.24900  ,0.25000  ,0.25100  ,0.25100  ,0.25200  ,0.25500  ,0.25900  ,0.26300  ,0.26383  ,0.26700  ,0.26627  ,0.27200  ,0.26257  ,0.27551  ,0.27600  ,0.27707  ,0.28100  ,0.28406  ,0.27351  ,0.28600  ,0.29100  ,0.28406  ,0.29600  ,0.29818  ,0.28406  ,0.30100  ,0.29426  ,0.30700  ,0.30757  ,0.31200  ,0.30414  ,0.31800  ,0.31375  ,0.32423  ,0.32500  ,0.33100  ,0.32309  ,0.33800  ,0.33919  ,0.33220  ,0.34400  ,0.34110  ,0.34518  ,0.35122  ,0.34979  ,0.35858  ,0.36600  ,0.35830  ,0.36664  ,0.37500  ,0.37481  ,0.38400  ,0.38283  ,0.39300  ,0.39071  ,0.40300  ,0.39846  ,0.41300  ,0.40608  ,0.41110  ,0.41600  ,0.42098  ,0.42098  ,0.41359  ,0.42400  ,0.42098  ,0.42802  ,0.42947  ,0.43306  ,0.43500  ,0.42826  ,0.42585  ,0.43592  ,0.43900  ,0.43545  ,0.44700  ,0.44628  ,0.44253  ,0.44953  ,0.45900  ,0.45598  ,0.45644  ,0.46700  ,0.46326  ,0.47100  ,0.47001  ,0.47667  ,0.47224  ,0.48800  ,0.48327  ,0.48979  ,0.49839  ,0.49625  ,0.50264  ,0.50900  ,0.50897  ,0.51524  ,0.52145  ,0.52145  ,0.53200  ,0.52760  ,0.53370  ,0.53975  ,0.54575  ,0.55200  ,0.55170  ,0.55170  ,0.55760  ,0.55367  ,0.56346  ,0.57119  ,0.56927  ,0.57503  ,0.58076  ,0.58076  ,0.58645  ,0.59209  ,0.59770  ,0.59584  ,0.60327  ,0.60881  ,0.61431  ,0.61977  ,0.62520  ,0.63060  ,0.63597  ,0.63597  ,0.59040  ,0.65364  ,0.66582  ,0.66236  ,0.67957  ,0.67100  ,0.68806  ,0.70486  ,0.70000  ,0.72140  ,0.73770  ,0.75857  ,0.71316  ,0.73770  ,0.75379  ,0.75000  ,0.76967  ,0.78535  ,0.80086  ,0.81619  ,0.80000  ,0.83738  ,0.83136  ,0.84637  ,0.86124  ,0.85000  ,0.87597  ,0.89057  ,0.91080  ,0.90505  ,0.90000  ,0.91940  ,0.93365  ,0.94778  ,0.94073  ,0.96181  ,0.95000  ,0.97851  ,0.97574  ,0.98958  ,1.00332  ,1.01698  ,1.00000  ,1.03055  ,1.04404  ,1.05746  ,1.05000  ,1.07080  ,1.08407  ,1.08407  ,1.09727  ,1.11040  ,1.10000  ,1.12346  ,1.13647  ,1.14942  ,1.16230  ,1.15000  ,1.17513  ,1.18791  ,1.20063  ,1.21331  ,1.20000  ,1.22593  ,1.23851  ,1.25103  ,1.25728  ,1.26352  ,1.25000  ,1.27596  ,1.28836  ,1.30071  ,1.31303  ,1.30000  ,1.32531  ,1.33755  ,1.34975  ,1.36191  ,1.35000  ,1.37405  ,1.38614  ,1.39821  ,1.41024  ,1.40000  ,1.42224  ,1.45000  ,1.55000  ,1.65000  ,1.75000  ,1.85000  ,1.95000  ,1.92223  ,2.10000  ,2.30000  ,2.50000  ,2.80000  ,2.50000  ,3.20000  ,2.99560  ,3.40000  ,4.30000  ,4.13410  ,5.10000  ,5.05296  ,5.66210  ,5.66214  ,6.50000  ,7.02680  ,5.00000  ,6.93650  ,8.33680  ,10.00000 ,7.00000  ,12.00000 ,13.00000 ,9.00000  ,14.60000 ,13.00000 ,10.00000 ,17.80000 ,18.00000 ,17.00000 ,21.60000 ,26.28321 ,26.50000 ,27.00000 ,34.00000 ,14.00000 ,36.58794 ,52.69162 ,80.00000 ,52.00000 ,131.00000,104.00000,104.00000,154.00000,180.00000,154.00000,215.00000,206.00000,240.00000,231.00000,273.00000,252.00000,200.00000};
    35. 

  35.   std::vector<double> v_plabMin_el_np     = {0.95000 ,1.68000 ,2.16000 ,2.50000 ,2.50000 ,3.71000 ,4.93000 ,1.00000 ,8.00000 ,10.00000};
    36. 

  36.   std::vector<double> v_plabMax_tot_pp    = {0.14000   ,0.19000   ,0.24000   ,0.28000   ,0.31000   ,0.35000   ,0.37000   ,0.39000   ,0.43000   ,0.44000   ,0.49000   ,0.54000   ,0.57000   ,0.59000   ,0.61026   ,0.62700   ,0.66541   ,0.69000   ,0.72000   ,0.75000   ,0.77700   ,0.76050   ,0.83000   ,0.84288   ,0.85700   ,0.85000   ,0.87473   ,0.89200   ,0.88000   ,0.93000   ,0.95700   ,0.94019   ,0.98300   ,0.96615   ,0.97935   ,0.97962   ,0.96824   ,0.99452   ,1.01100   ,1.00549   ,1.01164   ,1.02900   ,1.02456   ,1.03794   ,1.03673   ,1.09008   ,1.11300   ,1.09601   ,1.11045   ,1.12800   ,1.12200   ,1.11322   ,1.12700   ,1.14970   ,1.14234   ,1.16425   ,1.18200   ,1.19365   ,1.22497   ,1.22500   ,1.21898   ,1.25326   ,1.24413   ,1.26909   ,1.28140   ,1.29158   ,1.29264   ,1.29388   ,1.33200   ,1.40930   ,1.41082   ,1.45550   ,1.43300   ,1.48700   ,1.51836   ,1.54229   ,1.54200   ,1.58637   ,1.59240   ,1.64000   ,1.61021   ,1.62800   ,1.66332   ,1.69900   ,1.69604   ,1.72455   ,1.77325   ,1.78356   ,1.80553   ,1.79700   ,1.86172   ,1.89800   ,1.93725   ,1.94388   ,1.95590   ,2.00100   ,2.03267   ,2.05512   ,2.10125   ,2.08316   ,2.21642   ,2.27195   ,2.28456   ,2.42384   ,2.45490   ,2.53175   ,2.59718   ,2.68536   ,2.70941   ,2.82661   ,2.82464   ,2.86271   ,2.97796   ,2.96392   ,3.04425   ,2.99999   ,3.06011   ,3.11622   ,3.13726   ,3.14828   ,3.27000   ,3.28355   ,3.30961   ,3.46550   ,3.45089   ,3.55309   ,3.66950   ,3.67024   ,3.73846   ,3.91582   ,4.04000   ,4.10000   ,4.04507   ,4.20000   ,4.24200   ,4.27353   ,4.50990   ,4.56110   ,4.79257   ,4.97593   ,5.23144   ,5.53000   ,5.53705   ,5.83565   ,5.83020   ,5.96493   ,6.10200   ,7.07000   ,7.75020   ,7.85067   ,7.86800   ,8.13600   ,8.20000   ,10.40000  ,10.17000  ,10.10000  ,10.00000  ,10.02000  ,10.11000  ,12.20400  ,12.00000  ,13.00000  ,14.23800  ,15.00000  ,16.55000  ,16.27200  ,18.60000  ,18.30600  ,19.33000  ,20.40000  ,19.80000  ,20.34000  ,20.00000  ,20.00000  ,22.50000  ,22.37400  ,23.00000  ,23.50000  ,24.00000  ,25.60000  ,24.50000  ,25.00000  ,26.42000  ,29.80000  ,30.00000  ,30.00000  ,32.30000  ,35.00000  ,35.00000  ,40.00000  ,40.00000  ,42.50000  ,45.00000  ,50.00000  ,50.00000  ,50.00000  ,50.00000  ,50.00000  ,52.20000  ,55.00000  ,60.00000  ,65.00000  ,69.00000  ,70.00000  ,72.99398  ,97.99551  ,101.00000 ,100.00000 ,100.00000 ,102.00000 ,120.00000 ,121.99639 ,146.99700 ,147.00000 ,150.00000 ,150.00000 ,170.00000 ,171.99743 ,176.75000 ,195.99775 ,200.00000 ,200.00000 ,205.00000 ,240.00000 ,280.00000 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,295.86181 ,300.00000 ,301.50000 ,303.00000 ,310.00000 ,340.00000 ,370.00000 ,405.00000 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,498.04998 ,501.30972 ,504.58697 ,659.34014 ,1063.80004,1063.83483,1073.38439,1087.78857,1484.68920,1484.68920,1484.68920,1484.68920,1484.68920,1490.30004,1507.28405,2067.38403,2074.10009,2080.68530,2094.02905,2094.02905,2100.71679,2114.12451,4.0594E+06,7944.20019,1.5850E+04,3.1624E+04,6.3097E+04,1.5849E+05,2.5704E+07,2.6112E+07,2.6112E+07,2.6112E+07,2.6112E+07,3.4104E+07,3.4104E+07,3.4104E+07,4.4668E+07,7.7625E+07,9.0046E+07,9.0046E+07,1.3490E+08,2.3442E+08,4.0738E+08,4.7961E+08,2.1151E+09,5.3288E+09};
    37. 

  37.   std::vector<double> v_plabMax_el_pp     = {0.85700   ,1.00002   ,1.01100   ,1.02456   ,1.03794   ,1.12200   ,1.12700   ,1.16810   ,1.22497   ,1.22500   ,1.23200   ,1.28140   ,1.29264   ,1.29400   ,1.33200   ,1.43300   ,1.48700   ,1.59240   ,1.54200   ,1.62800   ,1.66167   ,1.67887   ,1.69900   ,1.69000   ,1.72455   ,1.79700   ,1.89800   ,2.00100   ,2.20000   ,2.22883   ,2.80000   ,2.80553   ,3.00000   ,3.06006   ,3.65000   ,3.69965   ,4.04000   ,4.20000   ,4.80000   ,5.00000   ,5.29986   ,5.53000   ,5.96493   ,6.00000   ,6.60990   ,6.79990   ,6.80000   ,6.99500   ,7.08794   ,7.07634   ,7.07634   ,8.20000   ,8.50000   ,8.80000   ,8.80000   ,8.90000   ,9.00000   ,9.12013   ,9.39140   ,9.89388   ,9.99900   ,10.02000  ,10.80000  ,10.90000  ,12.00000  ,12.00000  ,12.10000  ,12.80000  ,13.20000  ,14.80000  ,14.90877  ,15.00000  ,15.25100  ,15.50000  ,15.50000  ,16.20000  ,16.70000  ,18.60000  ,18.90000  ,19.00000  ,19.20000  ,19.60000  ,20.20000  ,21.40000  ,21.70000  ,23.50000  ,24.00000  ,24.50000  ,24.60000  ,26.20000  ,27.50000  ,29.99700  ,30.00000  ,30.50000  ,33.30000  ,36.20000  ,38.00000  ,40.60000  ,45.20000  ,50.30000  ,50.00000  ,50.00000  ,50.60000  ,52.10000  ,54.40000  ,57.00000  ,58.00000  ,60.00000  ,60.20000  ,63.50000  ,66.10000  ,69.20000  ,69.80001  ,70.42002  ,70.00000  ,100.60000 ,100.00000 ,102.00000 ,140.83999 ,147.00000 ,150.00000 ,176.05000 ,200.00000 ,200.00000 ,205.00000 ,205.00000 ,250.00000 ,293.35186 ,293.35186 ,293.35186 ,293.35186 ,300.00000 ,301.50000 ,303.00000 ,491.54168 ,498.04302 ,498.04302 ,498.04302 ,498.04302 ,1073.38439,1073.38439,1484.68920,1484.68920,1484.68920,1484.68920,2067.38403,2080.68530,4.0594E+06,2.6112E+07,2.6112E+07,3.4104E+07,3.4104E+07,9.0046E+07,9.0046E+07};
    38. 

  38.   std::vector<double> v_plabMax_tot_np    = {0.00531  ,0.03043  ,0.03873  ,0.04347  ,0.04502  ,0.04973  ,0.05000  ,0.05020  ,0.05100  ,0.05200  ,0.05300  ,0.05311  ,0.05400  ,0.05448  ,0.05500  ,0.05600  ,0.05700  ,0.05800  ,0.05900  ,0.05970  ,0.06000  ,0.06100  ,0.06200  ,0.06300  ,0.06400  ,0.06500  ,0.06571  ,0.06600  ,0.06700  ,0.06800  ,0.06900  ,0.06907  ,0.06913  ,0.07000  ,0.07100  ,0.07121  ,0.07200  ,0.07300  ,0.07400  ,0.07500  ,0.07600  ,0.07633  ,0.07700  ,0.07767  ,0.07800  ,0.07900  ,0.08000  ,0.08100  ,0.08112  ,0.08200  ,0.08300  ,0.08400  ,0.08500  ,0.08570  ,0.08600  ,0.08700  ,0.08787  ,0.08800  ,0.08900  ,0.08999  ,0.09000  ,0.09100  ,0.09200  ,0.09300  ,0.09400  ,0.09410  ,0.09468  ,0.09500  ,0.09600  ,0.09700  ,0.09798  ,0.09800  ,0.09900  ,0.10000  ,0.10100  ,0.10176  ,0.10200  ,0.10300  ,0.10400  ,0.10500  ,0.10541  ,0.10600  ,0.10700  ,0.10800  ,0.10900  ,0.11000  ,0.11100  ,0.11200  ,0.11245  ,0.11300  ,0.11400  ,0.11500  ,0.11564  ,0.11600  ,0.11630  ,0.11700  ,0.11800  ,0.11900  ,0.12000  ,0.12100  ,0.12200  ,0.12287  ,0.12300  ,0.12400  ,0.12500  ,0.12600  ,0.12700  ,0.12800  ,0.12867  ,0.12900  ,0.13000  ,0.13047  ,0.13100  ,0.13200  ,0.13400  ,0.13500  ,0.13600  ,0.13700  ,0.13745  ,0.13800  ,0.13900  ,0.14000  ,0.14032  ,0.14100  ,0.14200  ,0.14300  ,0.14400  ,0.14419  ,0.14500  ,0.14505  ,0.14600  ,0.14700  ,0.14800  ,0.14900  ,0.15000  ,0.15064  ,0.15100  ,0.15200  ,0.15300  ,0.15377  ,0.15400  ,0.15500  ,0.15600  ,0.15681  ,0.15700  ,0.15762  ,0.15800  ,0.15900  ,0.15985  ,0.16000  ,0.16100  ,0.16200  ,0.16280  ,0.16280  ,0.16292  ,0.16300  ,0.16339  ,0.16628  ,0.16400  ,0.16500  ,0.16600  ,0.16700  ,0.16800  ,0.16853  ,0.16856  ,0.16900  ,0.17000  ,0.17100  ,0.17137  ,0.17200  ,0.17300  ,0.17300  ,0.17400  ,0.17413  ,0.17500  ,0.17500  ,0.17600  ,0.17686  ,0.17700  ,0.17700  ,0.17800  ,0.17800  ,0.17900  ,0.17954  ,0.18000  ,0.18000  ,0.18100  ,0.18200  ,0.18200  ,0.18218  ,0.18300  ,0.18375  ,0.18400  ,0.18400  ,0.18479  ,0.18500  ,0.18600  ,0.18600  ,0.18700  ,0.18800  ,0.18800  ,0.18900  ,0.19000  ,0.19000  ,0.19100  ,0.19200  ,0.19287  ,0.19281  ,0.19300  ,0.19300  ,0.19336  ,0.19324  ,0.19400  ,0.19519  ,0.19500  ,0.19500  ,0.19519  ,0.19600  ,0.19700  ,0.19700  ,0.19782  ,0.19800  ,0.19884  ,0.19900  ,0.19900  ,0.20000  ,0.20100  ,0.20200  ,0.20200  ,0.20243  ,0.20300  ,0.20400  ,0.20400  ,0.20418  ,0.20500  ,0.20600  ,0.20700  ,0.20700  ,0.20800  ,0.20900  ,0.20900  ,0.20923  ,0.21000  ,0.21100  ,0.21100  ,0.21200  ,0.21300  ,0.21357  ,0.21400  ,0.21400  ,0.21500  ,0.21600  ,0.21609  ,0.21654  ,0.21700  ,0.21700  ,0.21800  ,0.21900  ,0.21950  ,0.22000  ,0.22000  ,0.22100  ,0.22200  ,0.22200  ,0.22263  ,0.22300  ,0.22395  ,0.22400  ,0.22500  ,0.22500  ,0.22600  ,0.22700  ,0.22800  ,0.22800  ,0.22900  ,0.22905  ,0.23000  ,0.23094  ,0.23100  ,0.23100  ,0.23200  ,0.23234  ,0.23300  ,0.23400  ,0.23400  ,0.23500  ,0.23524  ,0.23600  ,0.23626  ,0.23700  ,0.23800  ,0.23800  ,0.23900  ,0.24000  ,0.24100  ,0.24100  ,0.24118  ,0.24200  ,0.24300  ,0.24400  ,0.24400  ,0.24500  ,0.24600  ,0.24700  ,0.24800  ,0.24800  ,0.24900  ,0.25000  ,0.25100  ,0.25100  ,0.25200  ,0.25500  ,0.25900  ,0.26300  ,0.26383  ,0.26700  ,0.27351  ,0.27200  ,0.28406  ,0.27551  ,0.27600  ,0.27707  ,0.28100  ,0.28406  ,0.29426  ,0.28600  ,0.29100  ,0.30414  ,0.29600  ,0.29818  ,0.31689  ,0.30100  ,0.31375  ,0.30700  ,0.30757  ,0.31200  ,0.32309  ,0.31800  ,0.33220  ,0.32423  ,0.32500  ,0.33100  ,0.34110  ,0.33800  ,0.33919  ,0.34979  ,0.34400  ,0.35830  ,0.35435  ,0.35122  ,0.36664  ,0.35858  ,0.36600  ,0.37481  ,0.38283  ,0.37500  ,0.39071  ,0.38400  ,0.39846  ,0.39300  ,0.40608  ,0.40300  ,0.41359  ,0.41300  ,0.42098  ,0.42098  ,0.41800  ,0.42098  ,0.42098  ,0.42826  ,0.42400  ,0.43545  ,0.43043  ,0.43663  ,0.43306  ,0.43500  ,0.44253  ,0.44953  ,0.44065  ,0.44100  ,0.44953  ,0.44700  ,0.44860  ,0.45644  ,0.46326  ,0.45900  ,0.46507  ,0.47001  ,0.46900  ,0.47667  ,0.47100  ,0.48327  ,0.48979  ,0.49410  ,0.49000  ,0.49625  ,0.50264  ,0.50687  ,0.50897  ,0.51524  ,0.51100  ,0.52145  ,0.52760  ,0.53370  ,0.54176  ,0.53400  ,0.53975  ,0.54575  ,0.55170  ,0.55760  ,0.55400  ,0.56346  ,0.56346  ,0.56927  ,0.57312  ,0.57503  ,0.57119  ,0.58076  ,0.58645  ,0.59209  ,0.59584  ,0.59770  ,0.60327  ,0.60881  ,0.62159  ,0.61431  ,0.61977  ,0.62520  ,0.63060  ,0.63597  ,0.64131  ,0.64661  ,0.65364  ,0.69665  ,0.67100  ,0.66582  ,0.69649  ,0.67957  ,0.68806  ,0.70486  ,0.72140  ,0.75000  ,0.73770  ,0.75379  ,0.75857  ,0.80854  ,0.78535  ,0.76967  ,0.80000  ,0.78535  ,0.80086  ,0.81619  ,0.83136  ,0.85000  ,0.83738  ,0.84637  ,0.86124  ,0.87597  ,0.90000  ,0.89057  ,0.90505  ,0.91080  ,0.91940  ,0.95000  ,0.93365  ,0.94778  ,0.96181  ,0.99646  ,0.97574  ,1.00000  ,0.97851  ,0.98958  ,1.00332  ,1.01698  ,1.03055  ,1.05000  ,1.04404  ,1.05746  ,1.07080  ,1.10000  ,1.08407  ,1.08407  ,1.09727  ,1.11040  ,1.12346  ,1.15000  ,1.13647  ,1.14942  ,1.16230  ,1.17513  ,1.20000  ,1.18791  ,1.20063  ,1.21331  ,1.22593  ,1.25000  ,1.23851  ,1.25103  ,1.26352  ,1.25728  ,1.27596  ,1.30000  ,1.28836  ,1.30071  ,1.31303  ,1.32531  ,1.35000  ,1.33755  ,1.34975  ,1.36191  ,1.37405  ,1.40000  ,1.38614  ,1.39821  ,1.41024  ,1.42224  ,1.45000  ,1.43421  ,1.55000  ,1.65000  ,1.75000  ,1.85000  ,1.95000  ,2.10000  ,2.35937  ,2.30000  ,2.50000  ,2.80000  ,3.20000  ,3.50000  ,3.60000  ,3.92870  ,4.60000  ,4.30000  ,5.35780  ,6.30000  ,6.77473  ,7.07740  ,7.07748  ,6.50000  ,8.53800  ,11.00000 ,9.45340  ,10.04600 ,10.00000 ,15.00000 ,12.00000 ,13.00000 ,19.00000 ,14.60000 ,19.00000 ,24.00000 ,17.80000 ,18.00000 ,25.00000 ,21.60000 ,26.28321 ,26.50000 ,27.00000 ,34.00000 ,52.00000 ,36.58794 ,52.69162 ,80.00000 ,104.00000,131.00000,154.00000,154.00000,206.00000,180.00000,206.00000,215.00000,231.00000,240.00000,252.00000,273.00000,300.00000,400.00000};
    39. 

  39.   std::vector<double> v_plabMax_el_np     = {1.55000 ,1.78000 ,2.30000 ,3.50000 ,3.50000 ,3.95000 ,5.27000 ,10.00000,8.00000 ,24.00000};
    40. 

  40.   std::vector<double> v_sigmErrUp_tot_pp  = {13.000  ,2.0000  ,1.0000  ,1.0000  ,6.0000  ,0.40000 ,2.3000  ,2.0000  ,1.3000  ,1.3000  ,0.80000 ,1.2000  ,1.0000  ,1.9000  ,0.24000 ,0.12000 ,2.0000  ,0.90000 ,1.8000  ,1.3000  ,0.10000 ,0.23000 ,1.0000  ,1.0000  ,0.40000 ,0.90000 ,0.23000 ,0.10000 ,0.30000 ,0.40000 ,0.10000 ,0.23000 ,0.10000 ,0.23000 ,1.0000  ,1.4000  ,0.70000 ,0.      ,0.40000 ,1.2000  ,0.29000 ,0.14000 ,2.0000  ,2.0000  ,0.40000 ,0.40000 ,0.15000 ,0.34000 ,0.36000 ,0.19000 ,0.50000 ,0.17000 ,7.0000  ,1.3000  ,0.50000 ,0.39000 ,0.24000 ,0.50000 ,3.0000  ,0.50000 ,0.50000 ,1.4000  ,0.50000 ,0.60000 ,1.1000  ,0.11300 ,2.0000  ,0.60000 ,0.60000 ,2.8000  ,5.2E-02 ,0.50000 ,0.60000 ,2.0000  ,1.6000  ,1.6000  ,0.60000 ,0.10000 ,0.50000 ,1.0200  ,5.8E-02 ,5.0000  ,5.8E-02 ,0.60000 ,3.0000  ,2.0000  ,0.82000 ,4.6E-02 ,1.6000  ,0.60000 ,4.1E-02 ,0.60000 ,1.5100  ,4.6E-02 ,4.1E-02 ,0.60000 ,1.6000  ,1.6000  ,1.1200  ,4.1E-02 ,4.6E-02 ,2.6000  ,4.1E-02 ,4.1E-02 ,4.1E-02 ,0.83000 ,4.1E-02 ,4.1E-02 ,4.1E-02 ,3.2000  ,4.1E-02 ,4.1E-02 ,0.10000 ,4.1E-02 ,0.46000 ,4.1E-02 ,4.1E-02 ,4.1E-02 ,4.1E-02 ,4.1E-02 ,0.90000 ,4.1E-02 ,4.1E-02 ,4.0000  ,4.1E-02 ,3.7E-02 ,0.43000 ,1.2000  ,4.1E-02 ,4.1E-02 ,1.5000  ,0.62000 ,4.1E-02 ,4.0000  ,0.30000 ,4.1E-02 ,0.70000 ,4.1E-02 ,3.7E-02 ,4.1E-02 ,3.2E-02 ,1.4000  ,4.1E-02 ,4.1E-02 ,0.60000 ,4.2000  ,0.60000 ,0.30000 ,1.1000  ,5.2E-02 ,0.60000 ,0.60000 ,0.20000 ,1.5000  ,0.60000 ,0.30000 ,0.30000 ,1.7000  ,0.30000 ,0.60000 ,0.10000 ,1.5000  ,0.60000 ,0.12000 ,1.5000  ,0.60000 ,1.5000  ,0.60000 ,0.30000 ,1.5000  ,0.60000 ,0.60000 ,0.12000 ,0.15000 ,1.5000  ,0.60000 ,0.42000 ,1.5000  ,0.10000 ,1.5000  ,0.80000 ,0.12000 ,0.30000 ,1.5000  ,0.12000 ,0.15000 ,0.60000 ,0.12000 ,5.0E-02 ,0.12000 ,0.15000 ,9.0E-02 ,0.12000 ,0.12000 ,5.0E-02 ,7.0E-02 ,0.15000 ,1.1800  ,0.12000 ,0.12000 ,0.12000 ,0.18000 ,0.53000 ,4.0E-02 ,1.1000  ,1.1000  ,1.3000  ,4.0E-02 ,6.0E-02 ,0.80000 ,4.0E-02 ,1.1000  ,1.1000  ,0.58000 ,4.0E-02 ,6.0E-02 ,4.0E-02 ,1.1000  ,1.1000  ,1.2000  ,4.0E-02 ,6.0E-02 ,1.0000  ,4.0E-02 ,4.0E-02 ,0.22000 ,0.30000 ,0.32000 ,0.25000 ,0.70000 ,0.70000 ,0.92000 ,0.55000 ,1.0000  ,6.0E-02 ,7.0E-02 ,6.0E-02 ,1.2000  ,0.14000 ,0.17000 ,0.33000 ,0.24000 ,0.80000 ,0.40000 ,0.20000 ,0.60000 ,0.47000 ,0.40000 ,0.24000 ,0.32000 ,0.80000 ,0.20000 ,0.15000 ,0.20000 ,0.33000 ,0.35000 ,0.40000 ,0.70000 ,0.31000 ,0.40000 ,0.36000 ,0.23000 ,0.38000 ,0.90000 ,3.5000  ,3.3000  ,3.0000  ,3.0000  ,4.0000  ,4.0000  ,4.0000  ,14.000  ,0.2000  ,0.0     ,0.0     ,0.38    ,0.0     ,0.18    ,0.0     ,16.000  ,18.000  ,0.0     ,0.0     ,26.000  ,27.000  ,34.000  ,15.000  ,13.000  ,48.000};
    41. 

  41.   std::vector<double> v_sigmErrUp_el_pp   = {0.50000,1.0000 ,0.50000,2.0000 ,2.0000 ,0.50000,5.0000 ,0.80000,2.0000 ,0.50000,1.6000 ,0.80000,2.0000 ,1.0000 ,0.50000,0.60000,2.0000 ,3.0000 ,0.50000,3.0000 ,2.3000 ,0.90000,0.50000,2.1000 ,3.0000 ,0.50000,0.50000,0.50000,3.4000 ,0.73000,1.0000 ,0.48000,0.70000,3.0000 ,0.60000,0.76000,0.30000,2.6000 ,1.2000 ,0.50000,2.0000 ,0.25000,2.1000 ,0.50000,0.33000,0.22000,0.60000,0.50000,2.0000 ,1.0000 ,0.90000,0.40000,0.40000,0.22000,0.30000,0.50000,0.32000,0.80000,0.80000,3.0000 ,0.50000,0.60000,0.22000,0.50000,0.23000,0.20000,1.7000 ,0.30000,0.29000,0.32000,4.0000 ,0.30000,0.50000,1.4000 ,0.29000,0.49000,0.37000,1.8000 ,0.17000,0.50000,1.3000 ,0.44000,0.50000,1.6000 ,0.16000,1.2000 ,0.20000,0.30000,0.16000,2.2000 ,0.15000,0.60000,0.20000,0.14000,0.14000,0.11000,0.10000,0.11000,0.11000,0.29000,0.20000,0.12000,0.12000,0.12000,0.11000,0.10000,8.0E-02,0.70000,0.10000,9.0E-02,9.0E-02,9.0E-02,0.10000,0.31000,0.20000,0.35000,9.0E-02,1.0000 ,0.28000,0.24000,0.11000,0.28000,0.13000,8.0E-02,0.44000,0.44000,0.23000,0.30000,8.0E-02,0.20000,0.19000,0.16000,0.52000,0.40000,0.60000,0.40000,8.0E-02,0.20000,0.17000,8.0E-02,0.30000,0.13000,8.0E-02,0.30000,9.0E-02,0.19000,0.10000,1.40000,1.10000,0.19000,1.40000,0.18000,0.0    ,0.0};
    42. 

  42.   std::vector<double> v_sigmErrUp_tot_np  = {100.00 ,11.200 ,40.000 ,18.000 ,30.000 ,16.000 ,4.5100 ,40.000 ,3.4400 ,3.8000 ,3.7400 ,22.000 ,3.6900 ,20.000 ,4.2100 ,3.6000 ,4.1000 ,3.5100 ,3.9900 ,24.000 ,3.9500 ,3.9000 ,3.8600 ,3.8000 ,4.6200 ,3.7400 ,22.000 ,3.6900 ,4.4800 ,3.6200 ,4.3800 ,4.0000 ,9.0000 ,4.3300 ,4.3100 ,27.000 ,3.4700 ,4.2200 ,4.1700 ,4.1300 ,4.1000 ,19.000 ,4.0700 ,6.8000 ,5.7000 ,4.0000 ,3.9800 ,3.9400 ,20.000 ,5.5200 ,3.8700 ,3.8300 ,5.3900 ,11.000 ,3.7900 ,5.3200 ,60.000 ,5.2900 ,3.7100 ,13.000 ,5.2000 ,2.1000 ,1.6200 ,1.5400 ,1.6000 ,12.000 ,6.6000 ,1.5900 ,1.5800 ,1.6500 ,11.000 ,1.6400 ,1.6300 ,1.6200 ,1.6200 ,15.000 ,1.6000 ,1.6900 ,1.6800 ,1.6700 ,13.000 ,1.6600 ,1.6200 ,1.6500 ,1.7500 ,1.6300 ,1.7300 ,1.7200 ,10.000 ,1.7100 ,1.7000 ,1.6900 ,10.000 ,1.8100 ,20.000 ,1.6700 ,1.8000 ,1.6500 ,1.7800 ,1.7700 ,1.7700 ,10.000 ,1.7500 ,1.9200 ,1.7400 ,1.7300 ,1.8900 ,1.7100 ,10.000 ,1.8700 ,1.8600 ,9.0000 ,1.8600 ,1.6900 ,1.8400 ,2.0300 ,1.8100 ,1.8100 ,9.0000 ,1.8000 ,2.0000 ,1.7900 ,20.000 ,1.9900 ,1.7800 ,1.9800 ,1.9100 ,11.000 ,1.9700 ,20.000 ,1.7500 ,1.9500 ,1.9400 ,1.9300 ,1.9300 ,11.000 ,1.9200 ,2.2100 ,1.9100 ,11.000 ,1.9000 ,1.8900 ,2.1700 ,12.000 ,1.8700 ,20.000 ,2.1600 ,1.8600 ,19.000 ,1.8600 ,2.1400 ,2.1300 ,60.000 ,7.0000 ,10.000 ,1.8400 ,12.000 ,5.0000 ,2.1100 ,2.1000 ,2.1000 ,2.1000 ,1.8100 ,10.000 ,10.000 ,2.0800 ,2.0800 ,2.0600 ,12.000 ,2.0700 ,30.000 ,2.0600 ,2.0500 ,30.000 ,2.0500 ,34.000 ,2.4900 ,15.000 ,2.0300 ,33.000 ,2.0300 ,32.000 ,2.0200 ,30.000 ,2.4600 ,29.000 ,2.0100 ,28.000 ,2.0100 ,20.000 ,2.4500 ,9.0000 ,1.9900 ,27.000 ,15.000 ,1.9800 ,27.000 ,2.4100 ,1.9800 ,28.000 ,2.4200 ,1.9700 ,40.000 ,2.4000 ,3.3900 ,1.9500 ,2.4000 ,12.000 ,2.3800 ,26.000 ,1.0900 ,3.0000 ,1.9400 ,10.000 ,24.000 ,2.3800 ,7.0000 ,2.3700 ,22.000 ,2.3700 ,10.000 ,1.9300 ,12.000 ,21.000 ,2.3500 ,2.3400 ,2.3400 ,20.000 ,2.3300 ,5.0000 ,1.9000 ,20.000 ,2.3100 ,11.000 ,2.3100 ,2.2900 ,20.000 ,2.2900 ,2.3000 ,19.000 ,3.2500 ,4.0000 ,2.2900 ,16.000 ,2.2900 ,2.2700 ,2.2600 ,1.7000 ,2.2600 ,15.000 ,3.1700 ,2.2400 ,6.0000 ,2.3000 ,14.000 ,2.2400 ,2.2300 ,2.2300 ,9.0000 ,13.500 ,2.2300 ,2.2100 ,12.700 ,3.1200 ,5.0000 ,2.2200 ,12.000 ,2.2100 ,12.300 ,2.2200 ,3.1100 ,2.1800 ,11.500 ,2.1800 ,3.0800 ,5.0000 ,2.1800 ,2.1000 ,11.000 ,2.1700 ,3.0100 ,10.000 ,2.1500 ,10.500 ,3.0300 ,2.1400 ,8.0000 ,3.0500 ,2.6000 ,2.1300 ,9.5000 ,2.1300 ,3.0000 ,2.1000 ,8.8000 ,3.0100 ,11.000 ,2.1000 ,3.0300 ,8.3000 ,2.1100 ,2.9300 ,3.0000 ,2.1000 ,7.4000 ,2.9500 ,2.0700 ,2.9500 ,6.8000 ,2.9400 ,2.9300 ,7.2000 ,6.7000 ,9.9000 ,1.8000 ,6.3000 ,7.6000 ,6.0000 ,4.1500 ,1.7000 ,5.7000 ,17.000 ,5.4000 ,7.0000 ,2.7200 ,5.2000 ,4.9000 ,2.0700 ,4.6000 ,1.6000 ,10.000 ,4.0000 ,1.6300 ,4.1000 ,1.5000 ,3.8000 ,1.3600 ,3.6000 ,1.1300 ,1.2000 ,3.3000 ,3.1000 ,0.97000,2.9000 ,2.2000 ,0.85000,2.8000 ,0.75000,3.0000 ,2.6000 ,0.69000,2.5000 ,2.3000 ,0.63000,0.58000,2.3000 ,0.54000,3.6000 ,0.50600,2.1000 ,0.47600,1.9000 ,0.45200,1.9000 ,0.43100,1.4000 ,1.7000 ,4.0000 ,7.6000 ,0.41300,1.7000 ,0.39700,5.0000 ,3.0000 ,1.5000 ,1.6000 ,0.38400,10.000 ,3.0000 ,1.5000 ,0.37300,1.7000 ,7.0000 ,0.36300,0.35300,1.6000 ,16.000 ,0.34600,1.4000 ,0.33900,2.1000 ,0.33300,0.32800,8.6000 ,1.3000 ,0.32400,0.32000,1.8000 ,0.31700,0.31400,1.2000 ,0.31100,0.30900,0.30800,5.6000 ,1.1000 ,0.30600,0.30500,0.30400,0.30300,1.0000 ,0.30200,1.2000 ,0.30200,8.3000 ,0.30200,2.6000 ,0.30200,0.30100,0.30100,1.6000 ,0.30100,0.30100,0.30100,12.000 ,0.30200,0.30300,0.30200,0.30300,0.30300,0.30400,0.30400,0.16700,0.90000,0.16700,1.2000 ,3.5000 ,4.0000 ,0.16800,0.16900,0.17000,1.0200 ,0.17100,0.17100,0.48000,1.5000 ,3.0000 ,0.17200,0.72000,0.17300,0.17400,0.17400,0.17400,0.60000,0.27000,0.17400,0.17300,0.17200,0.57000,0.17100,0.16900,0.65000,0.16700,0.26000,0.16500,0.16200,0.16000,1.3000 ,0.15800,0.24000,0.73000,0.15600,0.15500,0.15300,0.15200,0.24000,0.15200,0.15200,0.15300,0.24000,0.15400,0.59000,0.15600,0.16000,0.16400,0.24000,0.16900,0.17600,0.18300,0.19100,0.23000,0.19900,0.20600,0.21300,0.21700,0.23000,0.21900,0.21800,0.21500,0.90000,0.20800,0.22000,0.19900,0.18700,0.17200,0.15400,0.22000,0.13400,0.11900,0.11200,0.11200,0.22000,0.11500,0.12000,0.13100,0.15500,0.21000,0.21600,0.16000,0.16000,0.16000,0.16000,0.17000,0.15000,1.8000 ,0.14000,0.16000,0.17000,0.25000,1.0000 ,0.46000,2.6000 ,0.60000,1.9000 ,1.6000 ,0.60000,1.6000 ,1.7000 ,1.7000 ,1.5000 ,1.7000 ,0.50000,2.6000 ,1.9000 ,0.50000,0.50000,0.60000,0.60000,0.50000,1.2000 ,0.60000,0.     ,1.2000 ,0.50000,0.70000,0.80000,0.20000,1.0000 ,0.60000,0.90000,0.90000,0.20000,0.15000,0.33000,0.33000,0.19000,0.24000,0.30000,0.36000,0.18000,0.20000,0.24000,0.24000,0.24000,0.24000,0.23000,0.23000,0.60000};
    43. 

  43.   std::vector<double> v_sigmErrUp_el_np   = {1.3600 ,1.3100 ,0.97000,2.3000 ,0.60000,0.83000,0.75000,0.80000,3.0E-02,1.3000 };
    44. 

  44.   std::vector<double> v_sigmErrLo_tot_pp  = {13.000 ,2.0000 ,1.0000 ,1.0000 ,6.0000 ,0.40000,2.3000 ,2.0000 ,1.3000 ,1.3000 ,0.80000,1.2000 ,1.0000 ,1.9000 ,0.24000,0.12000,2.0000 ,0.90000,1.8000 ,1.3000 ,0.10000,0.23000,1.0000 ,1.0000 ,0.40000,0.90000,0.23000,0.10000,0.30000,0.40000,0.10000,0.23000,0.10000,0.23000,1.0000 ,1.3000 ,0.70000,0.     ,0.40000,1.2000 ,0.29000,0.14000,2.0000 ,2.0000 ,0.40000,0.40000,0.15000,0.34000,0.36000,0.19000,0.50000,0.17000,7.0000 ,1.1000 ,0.50000,0.39000,0.24000,0.50000,3.0000 ,0.50000,0.50000,1.0000 ,0.50000,0.60000,1.1000 ,0.11300,2.0000 ,0.60000,0.60000,2.6000 ,5.2E-02,0.45000,0.60000,2.0000 ,1.2000 ,1.2000 ,0.60000,0.10000,0.50000,0.61000,5.8E-02,5.0000 ,5.8E-02,0.60000,3.0000 ,2.0000 ,0.46000,4.6E-02,1.1000 ,0.60000,4.1E-02,0.60000,0.68000,4.6E-02,4.1E-02,0.60000,1.2000 ,1.1000 ,0.47000,4.1E-02,4.6E-02,1.2000 ,4.1E-02,4.1E-02,4.1E-02,0.45000,4.1E-02,4.1E-02,4.1E-02,1.4000 ,4.1E-02,4.1E-02,0.10000,4.1E-02,0.42000,4.1E-02,4.1E-02,4.1E-02,4.1E-02,4.1E-02,0.90000,4.1E-02,4.1E-02,1.6000 ,4.1E-02,3.7E-02,0.43000,1.2000 ,4.1E-02,4.1E-02,1.5000 ,0.62000,4.1E-02,4.0000 ,0.30000,4.1E-02,0.70000,4.1E-02,3.7E-02,4.1E-02,3.2E-02,1.4000 ,4.1E-02,4.1E-02,0.60000,4.2000 ,0.60000,0.30000,1.1000 ,5.2E-02,0.60000,0.60000,0.20000,1.5000 ,0.60000,0.30000,0.30000,1.7000 ,0.30000,0.60000,0.10000,1.5000 ,0.60000,0.12000,1.5000 ,0.60000,1.5000 ,0.60000,0.30000,1.5000 ,0.60000,0.60000,0.12000,0.15000,1.5000 ,0.60000,0.42000,1.5000 ,0.10000,1.5000 ,0.80000,0.12000,0.30000,1.5000 ,0.12000,0.15000,0.60000,0.12000,5.0E-02,0.12000,0.15000,9.0E-02,0.12000,0.12000,5.0E-02,7.0E-02,0.15000,1.1800 ,0.12000,0.12000,0.12000,0.18000,0.53000,4.0E-02,1.1000 ,1.1000 ,1.3000 ,4.0E-02,6.0E-02,0.80000,4.0E-02,1.1000 ,1.1000 ,0.58000,4.0E-02,6.0E-02,4.0E-02,1.1000 ,1.1000 ,1.2000 ,4.0E-02,6.0E-02,1.0000 ,4.0E-02,4.0E-02,0.22000,0.30000,0.32000,0.25000,0.70000,0.70000,0.92000,0.55000,1.0000 ,6.0E-02,7.0E-02,6.0E-02,1.2000 ,0.14000,0.17000,0.33000,0.24000,0.80000,0.40000,0.20000,0.60000,0.47000,0.40000,0.24000,0.32000,0.80000,0.20000,0.15000,0.20000,0.33000,0.35000,0.40000,0.70000,0.31000,0.40000,0.36000,0.23000,0.38000,0.90000,3.5000 ,3.3000 ,3.0000 ,3.0000 ,4.0000 ,4.0000 ,4.0000 ,14.000 ,0.2000 ,0.0    ,0.0    ,0.38   ,0.0    ,0.18   ,0.0    ,16.000 ,18.000 ,0.0    ,0.0    ,26.000 ,27.000 ,34.000 ,15.000 ,13.000 ,44.000};
    45. 

  45.   std::vector<double> v_sigmErrLo_el_pp   = {0.50000,1.0000 ,0.50000,2.0000 ,2.0000 ,0.50000,5.0000 ,0.80000,2.0000 ,0.50000,1.6000 ,0.80000,2.0000 ,1.0000 ,0.50000,0.60000,2.0000 ,3.0000 ,0.50000,3.0000 ,2.3000 ,0.90000,0.50000,2.1000 ,3.0000 ,0.50000,0.50000,0.50000,3.4000 ,0.64000,1.0000 ,0.48000,0.70000,3.0000 ,0.60000,0.76000,0.30000,2.6000 ,1.2000 ,0.50000,2.0000 ,0.25000,2.1000 ,0.50000,0.33000,0.22000,0.60000,0.50000,2.0000 ,1.0000 ,0.90000,0.40000,0.40000,0.22000,0.30000,0.50000,0.32000,0.80000,0.80000,3.0000 ,0.50000,0.60000,0.22000,0.50000,0.23000,0.20000,1.7000 ,0.30000,0.29000,0.32000,4.0000 ,0.30000,0.50000,1.4000 ,0.29000,0.49000,0.37000,1.8000 ,0.17000,0.50000,1.3000 ,0.44000,0.50000,1.6000 ,0.16000,0.70000,0.20000,0.30000,0.16000,2.2000 ,0.15000,0.60000,0.20000,0.14000,0.14000,0.11000,0.10000,0.11000,0.11000,0.29000,0.20000,0.12000,0.12000,0.12000,0.11000,0.10000,8.0E-02,0.70000,0.10000,9.0E-02,9.0E-02,9.0E-02,0.10000,0.31000,0.20000,0.35000,9.0E-02,1.0000 ,0.28000,0.24000,0.11000,0.28000,0.13000,8.0E-02,0.44000,0.44000,0.23000,0.30000,8.0E-02,0.20000,0.19000,0.16000,0.52000,0.40000,0.60000,0.40000,8.0E-02,0.20000,0.17000,8.0E-02,0.30000,0.13000,8.0E-02,0.30000,9.0E-02,0.19000,0.10000,1.40000,1.10000,0.19000,1.40000,0.18000,0.0    ,0.0};
    46. 

  46.   std::vector<double> v_sigmErrLo_tot_np  = {100.00 ,11.200 ,40.000 ,18.000 ,30.000 ,16.000 ,4.5100 ,40.000 ,3.4400 ,3.8000 ,3.7400 ,22.000 ,3.6900 ,20.000 ,4.2100 ,3.6000 ,4.1000 ,3.5100 ,3.9900 ,24.000 ,3.9500 ,3.9000 ,3.8600 ,3.8000 ,4.6200 ,3.7400 ,22.000 ,3.6900 ,4.4800 ,3.6200 ,4.3800 ,4.0000 ,9.0000 ,4.3300 ,4.3100 ,27.000 ,3.4700 ,4.2200 ,4.1700 ,4.1300 ,4.1000 ,19.000 ,4.0700 ,6.8000 ,5.7000 ,4.0000 ,3.9800 ,3.9400 ,20.000 ,5.5200 ,3.8700 ,3.8300 ,5.3900 ,11.000 ,3.7900 ,5.3200 ,60.000 ,5.2900 ,3.7100 ,13.000 ,5.2000 ,2.1000 ,1.6200 ,1.5400 ,1.6000 ,12.000 ,6.6000 ,1.5900 ,1.5800 ,1.6500 ,11.000 ,1.6400 ,1.6300 ,1.6200 ,1.6200 ,15.000 ,1.6000 ,1.6900 ,1.6800 ,1.6700 ,13.000 ,1.6600 ,1.6200 ,1.6500 ,1.7500 ,1.6300 ,1.7300 ,1.7200 ,10.000 ,1.7100 ,1.7000 ,1.6900 ,10.000 ,1.8100 ,20.000 ,1.6700 ,1.8000 ,1.6500 ,1.7800 ,1.7700 ,1.7700 ,10.000 ,1.7500 ,1.9200 ,1.7400 ,1.7300 ,1.8900 ,1.7100 ,10.000 ,1.8700 ,1.8600 ,9.0000 ,1.8600 ,1.6900 ,1.8400 ,2.0300 ,1.8100 ,1.8100 ,9.0000 ,1.8000 ,2.0000 ,1.7900 ,20.000 ,1.9900 ,1.7800 ,1.9800 ,1.9100 ,11.000 ,1.9700 ,20.000 ,1.7500 ,1.9500 ,1.9400 ,1.9300 ,1.9300 ,11.000 ,1.9200 ,2.2100 ,1.9100 ,11.000 ,1.9000 ,1.8900 ,2.1700 ,12.000 ,1.8700 ,20.000 ,2.1600 ,1.8600 ,19.000 ,1.8600 ,2.1400 ,2.1300 ,60.000 ,7.0000 ,10.000 ,1.8400 ,12.000 ,5.0000 ,2.1100 ,2.1000 ,2.1000 ,2.1000 ,1.8100 ,10.000 ,10.000 ,2.0800 ,2.0800 ,2.0600 ,12.000 ,2.0700 ,30.000 ,2.0600 ,2.0500 ,30.000 ,2.0500 ,34.000 ,2.4900 ,15.000 ,2.0300 ,33.000 ,2.0300 ,32.000 ,2.0200 ,30.000 ,2.4600 ,29.000 ,2.0100 ,28.000 ,2.0100 ,20.000 ,2.4500 ,9.0000 ,1.9900 ,27.000 ,15.000 ,1.9800 ,27.000 ,2.4100 ,1.9800 ,28.000 ,2.4200 ,1.9700 ,40.000 ,2.4000 ,3.3900 ,1.9500 ,2.4000 ,12.000 ,2.3800 ,26.000 ,1.0900 ,3.0000 ,1.9400 ,10.000 ,24.000 ,2.3800 ,7.0000 ,2.3700 ,22.000 ,2.3700 ,10.000 ,1.9300 ,12.000 ,21.000 ,2.3500 ,2.3400 ,2.3400 ,20.000 ,2.3300 ,5.0000 ,1.9000 ,20.000 ,2.3100 ,11.000 ,2.3100 ,2.2900 ,20.000 ,2.2900 ,2.3000 ,19.000 ,3.2500 ,4.0000 ,2.2900 ,16.000 ,2.2900 ,2.2700 ,2.2600 ,1.7000 ,2.2600 ,15.000 ,3.1700 ,2.2400 ,6.0000 ,2.3000 ,14.000 ,2.2400 ,2.2300 ,2.2300 ,9.0000 ,13.500 ,2.2300 ,2.2100 ,12.700 ,3.1200 ,5.0000 ,2.2200 ,12.000 ,2.2100 ,12.300 ,2.2200 ,3.1100 ,2.1800 ,11.500 ,2.1800 ,3.0800 ,5.0000 ,2.1800 ,2.1000 ,11.000 ,2.1700 ,3.0100 ,10.000 ,2.1500 ,10.500 ,3.0300 ,2.1400 ,8.0000 ,3.0500 ,2.6000 ,2.1300 ,9.5000 ,2.1300 ,3.0000 ,2.1000 ,8.8000 ,3.0100 ,11.000 ,2.1000 ,3.0300 ,8.3000 ,2.1100 ,2.9300 ,3.0000 ,2.1000 ,7.4000 ,2.9500 ,2.0700 ,2.9500 ,6.8000 ,2.9400 ,2.9300 ,7.2000 ,6.7000 ,9.9000 ,1.8000 ,6.3000 ,7.6000 ,6.0000 ,4.1500 ,1.7000 ,5.7000 ,17.000 ,5.4000 ,7.0000 ,2.7200 ,5.2000 ,4.9000 ,2.0700 ,4.6000 ,1.6000 ,10.000 ,4.0000 ,1.6300 ,4.1000 ,1.5000 ,3.8000 ,1.3600 ,3.6000 ,1.1300 ,1.2000 ,3.3000 ,3.1000 ,0.97000,2.9000 ,2.2000 ,0.85000,2.8000 ,0.75000,3.0000 ,2.6000 ,0.69000,2.5000 ,2.3000 ,0.63000,0.58000,2.3000 ,0.54000,3.6000 ,0.50600,2.1000 ,0.47600,1.9000 ,0.45200,1.9000 ,0.43100,1.4000 ,1.7000 ,4.0000 ,7.6000 ,0.41300,1.7000 ,0.39700,5.0000 ,3.0000 ,1.5000 ,1.6000 ,0.38400,10.000 ,3.0000 ,1.5000 ,0.37300,1.7000 ,7.0000 ,0.36300,0.35300,1.6000 ,16.000 ,0.34600,1.4000 ,0.33900,2.1000 ,0.33300,0.32800,8.6000 ,1.3000 ,0.32400,0.32000,1.8000 ,0.31700,0.31400,1.2000 ,0.31100,0.30900,0.30800,5.6000 ,1.1000 ,0.30600,0.30500,0.30400,0.30300,1.0000 ,0.30200,1.2000 ,0.30200,8.3000 ,0.30200,2.6000 ,0.30200,0.30100,0.30100,1.6000 ,0.30100,0.30100,0.30100,12.000 ,0.30200,0.30300,0.30200,0.30300,0.30300,0.30400,0.30400,0.16700,0.90000,0.16700,1.2000 ,3.5000 ,4.0000 ,0.16800,0.16900,0.17000,1.0200 ,0.17100,0.17100,0.48000,1.5000 ,3.0000 ,0.17200,0.72000,0.17300,0.17400,0.17400,0.17400,0.60000,0.27000,0.17400,0.17300,0.17200,0.57000,0.17100,0.16900,0.65000,0.16700,0.26000,0.16500,0.16200,0.16000,1.3000 ,0.15800,0.24000,0.73000,0.15600,0.15500,0.15300,0.15200,0.24000,0.15200,0.15200,0.15300,0.24000,0.15400,0.59000,0.15600,0.16000,0.16400,0.24000,0.16900,0.17600,0.18300,0.19100,0.23000,0.19900,0.20600,0.21300,0.21700,0.23000,0.21900,0.21800,0.21500,0.90000,0.20800,0.22000,0.19900,0.18700,0.17200,0.15400,0.22000,0.13400,0.11900,0.11200,0.11200,0.22000,0.11500,0.12000,0.13100,0.15500,0.21000,0.21600,0.16000,0.16000,0.16000,0.16000,0.17000,0.15000,1.8000 ,0.14000,0.16000,0.17000,0.25000,1.0000 ,0.46000,2.6000 ,0.60000,1.9000 ,1.6000 ,0.60000,1.6000 ,1.7000 ,1.7000 ,1.5000 ,1.7000 ,0.50000,2.6000 ,1.9000 ,0.50000,0.50000,0.60000,0.60000,0.50000,1.2000 ,0.60000,0.     ,1.2000 ,0.50000,0.70000,0.80000,0.20000,1.0000 ,0.60000,0.90000,0.90000,0.20000,0.15000,0.33000,0.33000,0.19000,0.24000,0.30000,0.36000,0.18000,0.20000,0.24000,0.24000,0.24000,0.24000,0.23000,0.23000,0.60000};
    47. 

  47.   std::vector<double> v_sigmErrLo_el_np   = {1.3600 ,1.3100 ,0.97000,2.3000 ,0.60000,0.83000,0.75000,0.80000,3.0E-02,1.3000 };
    48. 

  48. 

  49.   std::vector<double> v_plabErrUp_tot_pp, v_plabErrUp_el_pp, v_plabErrUp_tot_np, v_plabErrUp_el_np;
    50. 

  50.   std::vector<double> v_plabErrLo_tot_pp, v_plabErrLo_el_pp, v_plabErrLo_tot_np, v_plabErrLo_el_np;
    51. 

  51. 

  52.   for (int i=0; i<v_plab_tot_pp.size(); i++){
    53. 

  53.     v_plabErrUp_tot_pp.push_back(v_plabMax_tot_pp.at(i) - v_plab_tot_pp.at(i));
    54. 

  54.     v_plabErrLo_tot_pp.push_back(v_plab_tot_pp.at(i) - v_plabMin_tot_pp.at(i));
    55. 

  55.   }
    56. 

  56.   for (int i=0; i<v_plab_el_pp.size(); ++i){
    57. 

  57.     v_plabErrUp_el_pp.push_back(v_plabMax_el_pp.at(i) - v_plab_el_pp.at(i));
    58. 

  58.     v_plabErrLo_el_pp.push_back(v_plab_el_pp.at(i) - v_plabMin_el_pp.at(i));
    59. 

  59.   }
    60. 

  60.   for (int i=0; i<v_plab_tot_np.size(); ++i){
    61. 

  61.     v_plabErrUp_tot_np.push_back(v_plabMax_tot_np.at(i) - v_plab_tot_np.at(i));
    62. 

  62.     v_plabErrLo_tot_np.push_back(v_plab_tot_np.at(i) - v_plabMin_tot_np.at(i));
    63. 

  63.   }
    64. 

  64.   for (int i=0; i<v_plab_el_np.size(); ++i){
    65. 

  65.     v_plabErrUp_el_np.push_back(v_plabMax_el_np.at(i) - v_plab_el_np.at(i));
    66. 

  66.     v_plabErrLo_el_np.push_back(v_plab_el_np.at(i) - v_plabMin_el_np.at(i));
    67. 

  67.   }
    68. 

  68. 

  69.   TGraphAsymmErrors *gr_tot_pp, *gr_el_pp, *gr_tot_np, *gr_el_np;
    70. 

  70. 

  71.   gr_tot_pp = new TGraphAsymmErrors(v_plab_tot_pp.size(), v_plab_tot_pp.data(), v_sigm_tot_pp.data(), v_plabErrLo_tot_pp.data(), v_plabErrUp_tot_pp.data(), v_sigmErrLo_tot_pp.data(), v_sigmErrUp_tot_pp.data());
    72. 

  72.   gr_el_pp = new TGraphAsymmErrors(v_plab_el_pp.size(), v_plab_el_pp.data(), v_sigm_el_pp.data(), v_plabErrLo_el_pp.data(), v_plabErrUp_el_pp.data(), v_sigmErrLo_el_pp.data(), v_sigmErrUp_el_pp.data());
    73. 

  73.   gr_tot_np = new TGraphAsymmErrors(v_plab_tot_np.size(), v_plab_tot_np.data(), v_sigm_tot_np.data(), v_plabErrLo_tot_np.data(), v_plabErrUp_tot_np.data(), v_sigmErrLo_tot_np.data(), v_sigmErrUp_tot_np.data());
    74. 

  74.   gr_el_np = new TGraphAsymmErrors(v_plab_el_np.size(), v_plab_el_np.data(), v_sigm_el_np.data(), v_plabErrLo_el_np.data(), v_plabErrUp_el_np.data(), v_sigmErrLo_el_np.data(), v_sigmErrUp_el_np.data());
    75. 

  75. 

  76.   // Check beam energy ranges
    77. 

  77.   auto it_plab_min_tot_pp = std::min_element(v_plab_tot_pp.begin(), v_plab_tot_pp.end());
    78. 

  78.   double plab_min_tot_pp = *it_plab_min_tot_pp;
    79. 

  79.   auto it_plab_max_tot_pp = std::max_element(v_plab_tot_pp.begin(), v_plab_tot_pp.end());
    80. 

  80.   double plab_max_tot_pp = *it_plab_max_tot_pp;
    81. 

  81.   auto it_plab_min_el_pp = std::min_element(v_plab_el_pp.begin(), v_plab_el_pp.end());
    82. 

  82.   double plab_min_el_pp = *it_plab_min_el_pp;
    83. 

  83.   auto it_plab_max_el_pp = std::max_element(v_plab_el_pp.begin(), v_plab_el_pp.end());
    84. 

  84.   double plab_max_el_pp = *it_plab_max_el_pp;
    85. 

  85.   auto it_plab_min_tot_np = std::min_element(v_plab_tot_np.begin(), v_plab_tot_np.end());
    86. 

  86.   double plab_min_tot_np = *it_plab_min_tot_np;
    87. 

  87.   auto it_plab_max_tot_np = std::max_element(v_plab_tot_np.begin(), v_plab_tot_np.end());
    88. 

  88.   double plab_max_tot_np = *it_plab_max_tot_np;
    89. 

  89.   auto it_plab_min_el_np = std::min_element(v_plab_el_np.begin(), v_plab_el_np.end());
    90. 

  90.   double plab_min_el_np = *it_plab_min_el_np;
    91. 

  91.   auto it_plab_max_el_np = std::max_element(v_plab_el_np.begin(), v_plab_el_np.end());
    92. 

  92.   double plab_max_el_np = *it_plab_max_el_np;
    93. 

  93. 

  94.   if (plab < plab_min_tot_pp) std::cerr << "WARNING: beam energy is SMALLER than available data set for sigma_tot_pp!!! Use at own risk." << std::endl;
    95. 

  95.   if (plab > plab_max_tot_pp) std::cerr << "WARNING: beam energy is LARGER than available data set for sigma_tot_pp!!! Use at own risk." << std::endl;
    96. 

  96.   if (plab < plab_min_el_pp) std::cerr << "WARNING: beam energy is SMALLER than available data set for sigma_el_pp!!! Use at own risk." << std::endl;
    97. 

  97.   if (plab > plab_max_el_pp) std::cerr << "WARNING: beam energy is LARGER than available data set for sigma_el_pp!!! Use at own risk." << std::endl;
    98. 

  98.   if (plab < plab_min_tot_np) std::cerr << "WARNING: beam energy is SMALLER than available data set for sigma_tot_np!!! Use at own risk." << std::endl;
    99. 

  99.   if (plab > plab_max_tot_np) std::cerr << "WARNING: beam energy is LARGER than available data set for sigma_tot_np!!! Use at own risk." << std::endl;
    100. 

  100.   if (plab < plab_min_el_np) std::cerr << "WARNING: beam energy is SMALLER than available data set for sigma_el_np!!! Use at own risk." << std::endl;
    101. 

  101.   if (plab > plab_max_el_np) std::cerr << "WARNING: beam energy is LARGER than available data set for sigma_el_np!!! Use at own risk." << std::endl;
    102. 

  102. 

  103.   double sigma_tot_pp = gr_tot_pp->Eval(plab);
    104. 

  104.   double sigma_el_pp = gr_el_pp->Eval(plab);
    105. 

  105.   double sigma_tot_np = gr_tot_np->Eval(plab);
    106. 

  106.   double sigma_el_np = gr_el_np->Eval(plab);
    107. 

  107. 

  108.   double sigma_inel_pp = sigma_tot_pp - sigma_el_pp;
    109. 

  109.   double sigma_inel_np = sigma_tot_np - sigma_el_np;
    110. 

  110.   // most probably no reliable data extrapolation can be made
    111. 

  111.   if (sigma_inel_np<0){ 
    112. 

  112.     sigma_inel_np = sigma_inel_pp;
    113. 

  113.     std::cerr << "WARNING: sigma_tot_np and sigma_el_np produce unreliable extrapolation at this beam energy!!! Setting sigma_inel_np = sigma_inel_pp." << std::endl;
    114. 

  114.   }
    115. 

  115. 

  116.   std::cout << "sigma_tot_pp = " << sigma_tot_pp << " mb, sigma_el_pp = " << sigma_el_pp << " mb, sigma_inel_pp = " << sigma_inel_pp << " mb" << std::endl;
    117. 

  117.   std::cout << "sigma_tot_np = " << sigma_tot_np << " mb, sigma_el_np = " << sigma_el_np << " mb, sigma_inel_np = " << sigma_inel_np << " mb" << std::endl;
    118. 

  118. 

  119.   sigma_inel_NN = ( N1*N2*sigma_inel_pp + Z1*Z2*sigma_inel_pp + (Z1*N2+Z2*N1)*sigma_inel_np )/( A1*A2 );
    120. 

  120. 

  121.   return sigma_inel_NN;
    122. 

  122. }
    123. 

  123. 

  124. double SigmaInelNN_CM(double sNN, double A1, double Z1, double A2, double Z2)
    125. 

  125. {
    126. 

  126.   double plab = CMToPLab(sNN);
    127. 

  127.   return SigmaInelNN_Plab(plab, A1, Z1, A2, Z2);
    128. 

  128. }
    129. 

  129. 

  130. double SigmaInelNN_Ekin(double Ekin, double A1, double Z1, double A2, double Z2)
    131. 

  131. {
    132. 

  132. 

  133.   double mp=0.938;
    134. 

  134.   double El=Ekin+mp;
    135. 

  135.   double plab=sqrt(El*El-mp*mp);
    136. 

  136.   return SigmaInelNN_Plab(plab, A1, Z1, A2, Z2);
    137. 

  137. }
    138.