double CMToPLab(double sNN = 4.5) { double M_NUCL = 0.938; double sN = sNN/2.; double pN = TMath::Sqrt(sN*sN - M_NUCL*M_NUCL); TLorentzVector tMom4(0, 0, -pN, sN); TVector3 vBoost{tMom4.BoostVector()}; TLorentzVector pMom4(0, 0, pN, sN); pMom4.Boost(-vBoost); return (double)pMom4.Pz(); } double SigmaInelNN_Plab(double plab, double A1, double Z1, double A2, double Z2) { double sigma_inel_NN = -1.; double N1 = A1-Z1; double N2 = A2-Z2; // Data points from PDG std::vector 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}; std::vector 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}; std::vector 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}; std::vector v_plab_el_np = {1.25000 ,1.73000 ,2.23000 ,3.00000 ,3.00000 ,3.83000 ,5.10000 ,5.50000 ,8.00000 ,17.00000}; std::vector 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}; std::vector 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}; std::vector 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}; std::vector 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 }; std::vector 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}; std::vector 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}; std::vector 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}; std::vector v_plabMin_el_np = {0.95000 ,1.68000 ,2.16000 ,2.50000 ,2.50000 ,3.71000 ,4.93000 ,1.00000 ,8.00000 ,10.00000}; std::vector 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}; std::vector 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}; std::vector 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}; std::vector v_plabMax_el_np = {1.55000 ,1.78000 ,2.30000 ,3.50000 ,3.50000 ,3.95000 ,5.27000 ,10.00000,8.00000 ,24.00000}; std::vector 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}; std::vector 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}; std::vector 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}; std::vector 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 }; std::vector 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}; std::vector 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}; std::vector 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}; std::vector 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 }; std::vector v_plabErrUp_tot_pp, v_plabErrUp_el_pp, v_plabErrUp_tot_np, v_plabErrUp_el_np; std::vector v_plabErrLo_tot_pp, v_plabErrLo_el_pp, v_plabErrLo_tot_np, v_plabErrLo_el_np; for (int i=0; i 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; 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; 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; 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; 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; 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; 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; double sigma_tot_pp = gr_tot_pp->Eval(plab); double sigma_el_pp = gr_el_pp->Eval(plab); double sigma_tot_np = gr_tot_np->Eval(plab); double sigma_el_np = gr_el_np->Eval(plab); double sigma_inel_pp = sigma_tot_pp - sigma_el_pp; double sigma_inel_np = sigma_tot_np - sigma_el_np; // most probably no reliable data extrapolation can be made if (sigma_inel_np<0){ sigma_inel_np = sigma_inel_pp; 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; } 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; 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; sigma_inel_NN = ( N1*N2*sigma_inel_pp + Z1*Z2*sigma_inel_pp + (Z1*N2+Z2*N1)*sigma_inel_np )/( A1*A2 ); return sigma_inel_NN; } double SigmaInelNN_CM(double sNN, double A1, double Z1, double A2, double Z2) { double plab = CMToPLab(sNN); return SigmaInelNN_Plab(plab, A1, Z1, A2, Z2); } double SigmaInelNN_Ekin(double Ekin, double A1, double Z1, double A2, double Z2) { double mp=0.938; double El=Ekin+mp; double plab=sqrt(El*El-mp*mp); return SigmaInelNN_Plab(plab, A1, Z1, A2, Z2); }