From 258306d442d871e8708652bcf7de8ac4185f1349 Mon Sep 17 00:00:00 2001 From: James Szalkie Date: Thu, 4 Jun 2026 12:16:54 -0400 Subject: [PATCH] sausage eliminated --- Armory/ClassSX3.h | 4 +- Armory/anasenMS.cpp | 56 ++- ELoss/E_vs_x_N-14.dat | 1000 ++++++++++++++++++------------------- ELoss/E_vs_x_N-17.dat | 1000 ++++++++++++++++++------------------- ELoss/E_vs_x_Ne-18.dat | 501 +++++++++++++++++++ ELoss/E_vs_x_alpha.dat | 1000 ++++++++++++++++++------------------- ELoss/E_vs_x_proton.dat | 1000 ++++++++++++++++++------------------- ELoss/PCEnergyAnalysis.py | 286 ++++------- 8 files changed, 2644 insertions(+), 2203 deletions(-) create mode 100644 ELoss/E_vs_x_Ne-18.dat diff --git a/Armory/ClassSX3.h b/Armory/ClassSX3.h index 3586292..9f43ba8 100644 --- a/Armory/ClassSX3.h +++ b/Armory/ClassSX3.h @@ -57,8 +57,8 @@ private: const int numDet = 12; const float radius = 88; const float width = 40; - const float length = 75; - const float gap = 46; + const float length = 174.3; // 75 + const float gap = 0; // 46 short id; // -1 when no hit short chUp; diff --git a/Armory/anasenMS.cpp b/Armory/anasenMS.cpp index faf0d54..93b5fe4 100644 --- a/Armory/anasenMS.cpp +++ b/Armory/anasenMS.cpp @@ -96,13 +96,15 @@ int main(int argc, char **argv){ if( argc >= 2 ) numEvent = atoi(argv[1]); TransferReaction transfer; - transfer.SetA(14, 7, 0); // 18Ne projectile - transfer.SetIncidentEnergyAngle((42.82/14), 0, 0); // KEA in MeV/u, theta and phi in rad + transfer.SetA(18, 10, 0); // 18Ne projectile + //transfer.SetIncidentEnergyAngle(0, 0, 0); // KEA in MeV/u, theta and phi in rad + TGraph* elossBeam = LoadELoss("../ELoss/E_vs_x_Ne-18.dat"); transfer.Seta(4, 2); // 4He target transfer.Setb(1, 1); // outgoing proton from the primary transfer - transfer.SetB(17, 8); // 21Na* heavy product + transfer.SetB(21, 11); // 21Na* heavy product + const double beamA = 18.0; // mass number of 18Ne beam - bool enableSequentialDecay = false; // turning to false to disable sequential decay for now, can be set to true to enable + bool enableSequentialDecay = true; // turning to false to disable sequential decay for now, can be set to true to enable const int decayDaughterA = 20; const int decayDaughterZ = 10; const int decayEjectA = 1; @@ -110,12 +112,14 @@ int main(int argc, char **argv){ // Excited state lists (projectile and heavy-product excitation states) std::vector ExAList = {0}; // 18Ne projectile excitations in MeV - std::vector ExList = {0}; // 21Na* excitation in MeV + std::vector ExList = {2.797, 2.829, 3.544, 3.862, 4.294}; // 21Na* excitation in MeV // define vertex position uniform distribution ranges (mm) - double vertexXRange[2] = { -5, 5}; // mm - double vertexYRange[2] = { -5, 5}; - double vertexZRange[2] = { -100, 100}; + double vertexXRange[2] = { -5, 5}; // mm - 5, 5 + double vertexYRange[2] = { -5, 5}; // -5, 5 + double vertexZRange[2] = { -100, 100}; // -174.3, 174.3 (full length of gas volume, centered at 0) + const double beamEntranceZ = -174.3; //vertexZRange[0]; // mm, assumed beam entrance into the gas + // detector resolution / uncertainty parameters double sigmaSX3_W = 0; // mm, if < 0 use mid-point (no spread in SX3 horizontal dimension) @@ -168,8 +172,17 @@ int main(int argc, char **argv){ // beam and CM variables saved in tree double KEA; double KEA2; - tree1->Branch("beamKEA", &KEA, "beamKEA/D"); - tree2->Branch("beamKEA", &KEA2, "beamKEA/D"); + double beamPath_cm; + double beamEnergy; + double beamEnergyLoss; + tree1->Branch("beamKEA", &KEA, "beamKEA/D"); + tree2->Branch("beamKEA", &KEA2, "beamKEA/D"); + tree1->Branch("beamPath_cm", &beamPath_cm, "beamPath_cm/D"); + tree2->Branch("beamPath_cm", &beamPath_cm, "beamPath_cm/D"); + tree1->Branch("beamEnergy", &beamEnergy, "beamEnergy/D"); + tree2->Branch("beamEnergy", &beamEnergy, "beamEnergy/D"); + tree1->Branch("beamEnergyLoss", &beamEnergyLoss, "beamEnergyLoss/D"); + tree2->Branch("beamEnergyLoss", &beamEnergyLoss, "beamEnergyLoss/D"); double thetaCM, phiCM; double thetaCM2, phiCM2; @@ -319,6 +332,22 @@ int main(int argc, char **argv){ // recalc kinematic constants for chosen states transfer.CalReactionConstant(); + // vertex position in target volume + vertexX = (vertexXRange[1]- vertexXRange[0])*gRandom->Rndm() + vertexXRange[0]; + vertexY = (vertexYRange[1]- vertexYRange[0])*gRandom->Rndm() + vertexYRange[0]; + vertexZ = (vertexZRange[1]- vertexZRange[0])*gRandom->Rndm() + vertexZRange[0]; + + TVector3 vertex(vertexX, vertexY, vertexZ); + + // compute beam energy at the event vertex from the gas path length + beamPath_cm = TVector3(vertexZ - beamEntranceZ, vertexX, vertexY).Mag() * 0.1; + if( beamPath_cm < 0 ) beamPath_cm = 0; + beamEnergy = elossBeam->Eval(beamPath_cm); // MeV + beamEnergyLoss = elossBeam->Eval(0.0) - beamEnergy; + KEA = beamEnergy / beamA; + transfer.SetIncidentEnergyAngle(KEA, 0, 0); + transfer.CalReactionConstant(); + // isotropic CM direction thetaCM = TMath::ACos(2 * gRandom->Rndm() - 1) ; // polar angle in CM frame phiCM = (gRandom->Rndm() - 0.5) * TMath::TwoPi(); @@ -364,13 +393,6 @@ int main(int argc, char **argv){ delete [] output; - // vertex position in target volume - vertexX = (vertexXRange[1]- vertexXRange[0])*gRandom->Rndm() + vertexXRange[0]; - vertexY = (vertexYRange[1]- vertexYRange[0])*gRandom->Rndm() + vertexYRange[0]; - vertexZ = (vertexZRange[1]- vertexZRange[0])*gRandom->Rndm() + vertexZRange[0]; - - TVector3 vertex(vertexX, vertexY, vertexZ); - // set direction vector from lab angle TVector3 dir(1, 0, 0); dir.SetTheta(thetab * TMath::DegToRad()); diff --git a/ELoss/E_vs_x_N-14.dat b/ELoss/E_vs_x_N-14.dat index d85584d..c86a859 100644 --- a/ELoss/E_vs_x_N-14.dat +++ b/ELoss/E_vs_x_N-14.dat @@ -1,501 +1,501 @@ Distance_cm Energy_MeV --0.0 30.0 -0.07180179111052108 29.9398997995992 -0.14352830066451855 29.879799599198396 -0.2151795230231735 29.819699398797596 -0.2867554524730064 29.759599198396792 -0.3582560832248319 29.699498997995992 -0.42968140941276456 29.63939879759519 -0.5010314250931679 29.579298597194388 -0.57230612424359 29.519198396793588 -0.6435055007617274 29.459098196392784 -0.7146295484643119 29.398997995991984 -0.7856782610860453 29.338897795591183 -0.8566516322784783 29.27879759519038 -0.9275496556088745 29.21869739478958 -0.9983723245590967 29.15859719438878 -1.0691196325244414 29.098496993987975 -1.1397915728124568 29.038396793587175 -1.210388138641787 28.97829659318637 -1.280909323140939 28.91819639278557 -1.3513551193470947 28.85809619238477 -1.4217255202048664 28.797995991983967 -1.4920205185650337 28.737895791583167 -1.5622401071833079 28.677795591182363 -1.6323842787190146 28.617695390781563 -1.7024530257338235 28.557595190380763 -1.7724463406904174 28.49749498997996 -1.8423642159511449 28.43739478957916 -1.9122066437766951 28.377294589178355 -1.9819736163246942 28.317194388777555 -2.0516651256483427 28.257094188376755 -2.1212811636949955 28.19699398797595 -2.19082172230472 28.13689378757515 -2.2602867932088757 28.07679358717435 -2.329676368028632 28.016693386773547 -2.3989904382734726 27.956593186372746 -2.4682289953397145 27.896492985971943 -2.5373920305089497 27.836392785571142 -2.6064795349465246 27.776292585170342 -2.6754914996999597 27.71619238476954 -2.7444279156973463 27.656092184368738 -2.813288773745763 27.595991983967934 -2.8820740645296095 27.535891783567134 -2.9507837786089777 27.475791583166334 -3.0194179064179623 27.41569138276553 -3.087976438262945 27.35559118236473 -3.1564593643208982 27.295490981963926 -3.224866674637605 27.235390781563126 -3.293198359125913 27.175290581162326 -3.361454407563926 27.115190380761522 -3.429634809593171 27.05509018036072 -3.497739554716773 26.99498997995992 -3.5657686322975746 26.934889779559118 -3.6337220315562218 26.874789579158318 -3.701599741569271 26.814689378757514 -3.7694017512672056 26.754589178356714 -3.8371280494324855 26.694488977955913 -3.9047786246975313 26.63438877755511 -3.9723534655426787 26.57428857715431 -4.039852560294146 26.514188376753506 -4.107275897121914 26.454088176352705 -4.174623464037633 26.393987975951905 -4.241895248892465 26.3338877755511 -4.309091239374895 26.2737875751503 -4.376211423008545 26.213687374749497 -4.443255787149914 26.153587174348697 -4.510224318986124 26.093486973947897 -4.577117005532616 26.033386773547093 -4.643933833630797 25.973286573146293 -4.7106747899457035 25.913186372745493 -4.777339860963585 25.85308617234469 -4.843929032989463 25.79298597194389 -4.910442292144685 25.732885771543085 -4.9768796243644 25.672785571142285 -5.043241015395039 25.612685370741485 -5.109526450791738 25.55258517034068 -5.175735915915712 25.49248496993988 -5.241869395931638 25.432384769539077 -5.307926875804942 25.372284569138277 -5.3739083402991 25.312184368737476 -5.4398137739728725 25.252084168336673 -5.505643161177487 25.191983967935872 -5.5713964860538345 25.13188376753507 -5.637073732529551 25.07178356713427 -5.70267488431613 25.011683366733468 -5.768199924905947 24.951583166332664 -5.833648837569244 24.891482965931864 -5.89902160535111 24.831382765531064 -5.964318211068373 24.77128256513026 -6.029538637306457 24.71118236472946 -6.094682866416231 24.651082164328656 -6.15975088051075 24.590981963927856 -6.224742661462017 24.530881763527056 -6.2896581908976446 24.470781563126252 -6.354497450197486 24.41068136272545 -6.419260420490242 24.350581162324648 -6.483947082649972 24.290480961923848 -6.548557417292601 24.230380761523048 -6.613091404772351 24.170280561122244 -6.6775490251781076 24.110180360721444 -6.741930258329784 24.05008016032064 -6.806235083774566 23.98997995991984 -6.8704634807831635 23.92987975951904 -6.934615428345972 23.869779559118236 -6.9986909051691715 23.809679358717435 -7.062689889670811 23.749579158316635 -7.126612359976794 23.68947895791583 -7.190458293916811 23.62937875751503 -7.254227669020249 23.569278557114227 -7.31792046251198 23.509178356713427 -7.381536651308153 23.449078156312627 -7.445076212011879 23.388977955911823 -7.508539120908857 23.328877755511023 -7.57192535396297 23.26877755511022 -7.635234886811764 23.20867735470942 -7.698467694761913 23.14857715430862 -7.761623752784579 23.088476953907815 -7.82470303551071 23.028376753507015 -7.887705517226298 22.96827655310621 -7.950631171867514 22.90817635270541 -8.013479973015828 22.84807615230461 -8.076251893893017 22.787975951903807 -8.138946907356104 22.727875751503007 -8.201564985892244 22.667775551102206 -8.264106101613512 22.607675350701403 -8.326570226251613 22.547575150300602 -8.388957331152536 22.4874749498998 -8.45126738727109 22.427374749499 -8.513500365165402 22.367274549098198 -8.575656234991307 22.307174348697394 -8.637734966496634 22.247074148296594 -8.69973652901547 22.18697394789579 -8.761660891462254 22.12687374749499 -8.823508022325866 22.06677354709419 -8.885277889663561 22.006673346693386 -8.946970461094837 21.946573146292586 -9.00858570379523 21.886472945891782 -9.070123584489966 21.826372745490982 -9.131584069447577 21.76627254509018 -9.19296712447337 21.706172344689378 -9.254272714902815 21.646072144288578 -9.315500805594846 21.585971943887778 -9.376651360925043 21.525871743486974 -9.437724344778701 21.465771543086174 -9.49871972054383 21.40567134268537 -9.559637451104003 21.34557114228457 -9.620477498831134 21.28547094188377 -9.681239825578121 21.225370741482966 -9.741924392671372 21.165270541082165 -9.80253116090325 21.10517034068136 -9.863060090524352 21.04507014028056 -9.92351114123572 20.98496993987976 -9.983884272180898 20.924869739478957 -10.04417944193787 20.864769539078157 -10.104396608510893 20.804669338677353 -10.164535729322173 20.744569138276553 -10.224596761203454 20.684468937875753 -10.284579660387434 20.62436873747495 -10.344484382499063 20.56426853707415 -10.404310882546723 20.50416833667335 -10.46405911491325 20.444068136272545 -10.523729033346811 20.383967935871745 -10.583320590951676 20.32386773547094 -10.642833740178787 20.26376753507014 -10.702268432816243 20.20366733466934 -10.76162461997959 20.143567134268537 -10.820902252101968 20.083466933867737 -10.880101278924126 20.023366733466933 -10.939221649484248 19.963266533066133 -10.998263312107644 19.903166332665332 -11.057226214396275 19.84306613226453 -11.11611030321808 19.78296593186373 -11.1749155246962 19.722865731462925 -11.233641824197953 19.662765531062124 -11.292289146323709 19.602665330661324 -11.350857434895538 19.54256513026052 -11.409346632945686 19.48246492985972 -11.467756682704895 19.42236472945892 -11.526087525590508 19.362264529058116 -11.58433910219439 19.302164328657316 -11.642511352270676 19.242064128256512 -11.70060421472329 19.181963927855712 -11.758617627593308 19.12186372745491 -11.816551528046077 19.061763527054108 -11.87440585235814 19.001663326653308 -11.932180535903985 18.941563126252504 -11.989875513142515 18.881462925851704 -12.047490717603374 18.821362725450903 -12.105026081873003 18.7612625250501 -12.162481537580478 18.7011623246493 -12.219857015383148 18.641062124248496 -12.277152444952002 18.580961923847696 -12.334367754956837 18.520861723446895 -12.391502873051166 18.46076152304609 -12.448557725856865 18.40066132264529 -12.505532238948621 18.34056112224449 -12.56242633683808 18.280460921843687 -12.61923994295775 18.220360721442887 -12.67597297964467 18.160260521042083 -12.732625368123767 18.100160320641283 -12.789197028490985 18.040060120240483 -12.84568787969612 17.97995991983968 -12.902097839525354 17.91985971943888 -12.958426824583556 17.859759519038075 -13.014674750276237 17.799659318637275 -13.07084153079126 17.739559118236475 -13.126927079080211 17.67945891783567 -13.182931306839478 17.61935871743487 -13.238854124491043 17.559258517034067 -13.294695441162911 17.499158316633267 -13.35045516466927 17.439058116232466 -13.406135057732337 17.378957915831663 -13.4617388016865 17.318857715430862 -13.51726839581522 17.258757515030062 -13.57272416217035 17.19865731462926 -13.628106295572564 17.13855711422846 -13.683414795676521 17.078456913827655 -13.73864966198905 17.018356713426854 -13.793810893866711 16.958256513026054 -13.848898490513296 16.89815631262525 -13.90391245097732 16.83805611222445 -13.958852774149522 16.777955911823646 -14.013719458760315 16.717855711422846 -14.068512503377256 16.657755511022046 -14.12323190640248 16.597655310621242 -14.177877666070119 16.537555110220442 -14.232449780443737 16.477454909819638 -14.286948247413727 16.417354709418838 -14.34137306469472 16.357254509018038 -14.395724229822974 16.297154308617234 -14.450001740153757 16.237054108216434 -14.504205592858739 16.176953907815633 -14.558335784923374 16.11685370741483 -14.612392313144266 16.05675350701403 -14.666375174126575 15.996653306613226 -14.720284364281374 15.936553106212424 -14.774119879823063 15.876452905811622 -14.82788171676675 15.81635270541082 -14.881569870925665 15.75625250501002 -14.935184337908577 15.696152304609218 -14.988725113117225 15.636052104208416 -15.042192191743764 15.575951903807614 -15.095585568768243 15.515851703406813 -15.148905238956097 15.455751503006011 -15.202151196855652 15.39565130260521 -15.255323436795692 15.335551102204407 -15.308421952883021 15.275450901803605 -15.361446739000097 15.215350701402805 -15.414397788802681 15.155250501002003 -15.467275095717538 15.095150300601201 -15.520078652940192 15.0350501002004 -15.572808453432728 14.974949899799599 -15.62546448992166 14.914849699398797 -15.678046754895849 14.854749498997995 -15.7305552406045 14.794649298597193 -15.782989939055238 14.734549098196391 -15.83535084201224 14.67444889779559 -15.887637940994487 14.614348697394789 -15.939851227274083 14.554248496993987 -15.991990691874676 14.494148296593185 -16.044056325569994 14.434048096192384 -16.0960481188825 14.373947895791582 -16.147966062082148 14.31384769539078 -16.199810145185282 14.253747494989979 -16.251580357953678 14.193647294589177 -16.303276689893725 14.133547094188376 -16.354899130255767 14.073446893787574 -16.406447668033614 14.013346693386772 -16.45792229196422 13.95324649298597 -16.509322990527565 13.89314629258517 -16.560649751946734 13.833046092184368 -16.611902564188192 13.772945891783566 -16.66308141496231 13.712845691382764 -16.714186291724115 13.652745490981962 -16.76521718167428 13.592645290581162 -16.81617407176042 13.53254509018036 -16.867056948678616 13.472444889779558 -16.917865798875297 13.412344689378756 -16.968600608549387 13.352244488977956 -17.01926136365482 13.292144288577154 -17.06984804990339 13.232044088176352 -17.120360652767964 13.17194388777555 -17.170799157486123 13.111843687374748 -17.221163549064173 13.051743486973947 -17.27145381228162 12.991643286573145 -17.321669931696093 12.931543086172343 -17.37181189164875 12.871442885771541 -17.421879676270205 12.811342685370741 -17.471873269486984 12.75124248496994 -17.521792655028538 12.691142284569137 -17.571637816434862 12.631042084168335 -17.621408737064744 12.570941883767533 -17.67110540010466 12.510841683366733 -17.72072778857838 12.450741482965931 -17.770275885357268 12.390641282565129 -17.81974967317139 12.330541082164327 -17.869149134621406 12.270440881763527 -17.91847425219131 12.210340681362725 -17.96772500826206 12.150240480961923 -18.016901385126154 12.090140280561121 -18.06600336500317 12.030040080160319 -18.115030930056353 11.969939879759519 -18.163984062410297 11.909839679358717 -18.21286274416973 11.849739478957915 -18.261666957439527 11.789639278557113 -18.31039668434597 11.729539078156312 -18.359051907059335 11.66943887775551 -18.407632607817842 11.609338677354708 -18.45613876895308 11.549238476953906 -18.504570372916948 11.489138276553104 -18.55292740231021 11.429038076152304 -18.601209839912713 11.368937875751502 -18.649417668715376 11.3088376753507 -18.697550871954036 11.248737474949898 -18.745609433145237 11.188637274549098 -18.793593336124054 11.128537074148296 -18.841502565084046 11.068436873747494 -18.88933710461947 11.008336673346692 -18.937096939769848 10.94823647294589 -18.984782056066983 10.88813627254509 -19.032392439584594 10.828036072144288 -19.079928076990655 10.767935871743486 -19.127388955602573 10.707835671342684 -19.174775063445388 10.647735470941884 -19.222086389313084 10.587635270541082 -19.26932292283319 10.52753507014028 -19.31648465453486 10.467434869739478 -19.363571575920545 10.407334669338676 -19.410583679541464 10.347234468937875 -19.45752095907709 10.287134268537073 -19.504383409418775 10.227034068136271 -19.551171026757775 10.16693386773547 -19.597883808677874 10.10683366733467 -19.644521754252835 10.046733466933867 -19.691084864148884 9.986633266533065 -19.73757314073253 9.926533066132263 -19.78398658818393 9.866432865731461 -19.830325212616096 9.806332665330661 -19.87658902220025 9.746232464929859 -19.922778027297554 9.686132264529057 -19.968892240597626 9.626032064128255 -20.014931677264105 9.565931863727455 -20.060896355087653 9.505831663326653 -20.106786294646717 9.44573146292585 -20.15260388337894 9.385631262525049 -20.19835324261332 9.325531062124247 -20.24403635438935 9.265430861723447 -20.289653709967453 9.205330661322645 -20.335205646669117 9.145230460921843 -20.380692334128472 9.08513026052104 -20.426113948586362 9.02503006012024 -20.471470673193984 8.964929859719438 -20.51676269833093 8.904829659318636 -20.56199022193833 8.844729458917834 -20.60715344986786 8.784629258517032 -20.65225259624735 8.724529058116232 -20.697287883863904 8.66442885771543 -20.742259544565325 8.604328657314628 -20.78716781968081 8.544228456913826 -20.832012960461878 8.484128256513026 -20.876795228544545 8.424028056112224 -20.9215148964338 8.363927855711422 -20.966172248011556 8.30382765531062 -21.010767579069267 8.243727454909818 -21.05530119786641 8.183627254509018 -21.099773425716265 8.123527054108216 -21.14418459760029 8.063426853707414 -21.18853506281266 8.003326653306614 -21.23282518563649 7.943226452905811 -21.27705534605338 7.88312625250501 -21.321225940488095 7.823026052104208 -21.365337382590127 7.7629258517034065 -21.409390104054193 7.7028256513026045 -21.453384555481648 7.6427254509018026 -21.497321207285044 7.5826252505010014 -21.541200550638127 7.5225250501001995 -21.585023098473723 7.462424849699398 -21.628789386532095 7.402324649298596 -21.672499974462514 7.342224448897795 -21.716155446980963 7.282124248496993 -21.75975641508704 7.222024048096192 -21.80330351734335 7.16192384769539 -21.84679742122083 7.101823647294588 -21.89023882451372 7.041723446893787 -21.93362845682805 6.981623246492985 -21.97696708114787 6.921523046092184 -22.020255495483557 6.861422845691382 -22.063494534606978 6.801322645290581 -22.106685071878456 6.741222444889779 -22.14982802117091 6.681122244488978 -22.192924338896823 6.621022044088176 -22.235975026144065 6.560921843687374 -22.278981130927104 6.500821643286573 -22.321943750560397 6.440721442885771 -22.36486403416138 6.3806212424849695 -22.407743185290894 6.3205210420841675 -22.450582464739462 6.260420841683366 -22.493383193468393 6.200320641282564 -22.536146755715347 6.140220440881763 -22.57886084602819 6.080120240480961 -22.621507998728188 6.020020040080159 -22.664083083613544 5.959919839679358 -22.706585162587967 5.899819639278556 -22.749014293993014 5.839719438877755 -22.7913711367681 5.779619238476953 -22.83365639649778 5.719519038076152 -22.875870828020997 5.65941883767535 -22.918015238182466 5.599318637274549 -22.9600904887345 5.539218436873747 -23.002097499398058 5.479118236472945 -23.04403725109238 5.419018036072144 -23.085910789343284 5.358917835671342 -23.12771922788074 5.298817635270541 -23.169463752437263 5.238717434869739 -23.211145624759308 5.178617234468938 -23.25276618684479 5.118517034068136 -23.294326865420835 5.0584168336673345 -23.335829176676782 4.9983166332665325 -23.37727473126868 4.9382164328657305 -23.418665239612682 4.878116232464929 -23.46000251748603 4.818016032064127 -23.50128849195584 4.757915831663326 -23.54252520765742 4.697815631262524 -23.58371483344552 4.637715430861723 -23.624859669443914 4.577615230460921 -23.665962154520678 4.51751503006012 -23.707024874218767 4.457414829659318 -23.748050569174072 4.397314629258516 -23.78904214405575 4.337214428857715 -23.83000267706666 4.277114228456913 -23.870935430045073 4.217014028056112 -23.911843859212375 4.15691382765531 -23.952731626615574 4.096813627254509 -23.99360261231784 4.036713426853707 -24.0344609273952 3.9766132264529053 -24.075310927803017 3.9165130260521037 -24.116157229181862 3.856412825651302 -24.157004722679194 3.7963126252505006 -24.197858591870666 3.736212424849699 -24.238724330873303 3.6761122244488975 -24.279607763752022 3.616012024048096 -24.32051506533147 3.555911823647294 -24.361452783536723 3.4958116232464924 -24.40242786339951 3.435711422845691 -24.443447672881266 3.3756112224448893 -24.484520030680844 3.3155110220440878 -24.525653236213284 3.255410821643286 -24.566856101967026 3.1953106212424847 -24.608137988470688 3.135210420841683 -24.64950884212734 3.0751102204408816 -24.690979236204775 3.0150100200400796 -24.732560415304846 2.954909819639278 -24.77426434367444 2.8948096192384765 -24.816103757765784 2.834709418837675 -24.85809222350526 2.7746092184368734 -24.900244198789036 2.714509018036072 -24.942575101791842 2.6544088176352703 -24.98510138575333 2.5943086172344687 -25.027840620997036 2.534208416833667 -25.070811585041426 2.474108216432865 -25.114034361784228 2.4140080160320636 -25.157530450882575 2.353907815631262 -25.201322888616833 2.2938076152304605 -25.245436381719497 2.233707414829659 -25.28989745587799 2.1736072144288574 -25.334734620888515 2.113507014028056 -25.379978554755372 2.0534068136272543 -25.42566230940736 1.9933066132264528 -25.47182154115234 1.9332064128256512 -25.518494769529443 1.8731062124248496 -25.565723668865896 1.813006012024048 -25.613553397627733 1.7529058116232463 -25.662032971603566 1.6928056112224448 -25.711215688119985 1.6327054108216432 -25.76115960991029 1.5726052104208417 -25.8119281190153 1.51250501002004 -25.863590553278588 1.4524048096192383 -25.916222940730588 1.3923046092184368 -25.969908850600266 1.3322044088176352 -26.024740384070153 1.2721042084168337 -26.080819333502628 1.212004008016032 -26.13825854613062 1.1519038076152304 -26.197183537712476 1.0918036072144288 -26.25773441423771 1.0317034068136273 -26.32006817665718 0.9716032064128256 -26.38436150660483 0.911503006012024 -26.450814162914718 0.8514028056112224 -26.51965316366565 0.7913026052104208 -26.59113799328843 0.7312024048096192 -26.665567170123857 0.6711022044088176 -26.743286655809747 0.611002004008016 -26.824700817961897 0.5509018036072144 -26.910287035101458 0.4908016032064128 -27.00061568214767 0.4307014028056112 -27.096378416212726 0.3706012024048096 -27.19891273784251 0.310501002004008 -27.31067456563479 0.2504008016032064 -27.43503429907699 0.1903006012024048 -27.576994110147854 0.1302004008016032 -27.74597338383011 0.0701002004008016 -27.950434782116325 0.01 +-0.0 42.82 +0.1253066290957655 42.73420841683367 +0.25046072905153527 42.64841683366734 +0.37546229782093876 42.56262525050101 +0.5003113333504398 42.476833667334674 +0.6250078335789916 42.39104208416834 +0.7495517964376327 42.30525050100201 +0.8739432198490851 42.219458917835674 +0.9981821017273447 42.13366733466934 +1.1222684399772636 42.04787575150301 +1.2462022324941247 41.96208416833667 +1.369983477163207 41.87629258517034 +1.493612171859333 41.790501002004014 +1.6170883144464583 41.70470941883768 +1.7404119027771503 41.61891783567135 +1.8635829346921604 41.533126252505014 +1.9866014080199357 41.44733466933868 +2.109467320576131 41.36154308617235 +2.2321806701631135 41.27575150300601 +2.3547414545694556 41.18995991983968 +2.47714967156942 41.10416833667335 +2.599405318922432 41.01837675350701 +2.721508394372536 40.93258517034069 +2.8434588956478883 40.84679358717435 +2.9652568204601404 40.76100200400802 +3.0869021665039122 40.67521042084169 +3.208394931456204 40.58941883767535 +3.3297351129758073 40.50362725450902 +3.450922708702706 40.41783567134269 +3.571957716257464 40.33204408817635 +3.692840133240601 40.24625250501002 +3.813569957231961 40.160460921843686 +3.934147185790053 40.07466933867736 +4.054571816451433 39.98887775551103 +4.174843846729972 39.90308617234469 +4.294963274116216 39.81729458917836 +4.414930096076671 39.731503006012026 +4.534744310053103 39.64571142284569 +4.654405913461808 39.55991983967936 +4.773914903692879 39.474128256513026 +4.893271278109454 39.38833667334669 +5.012475034046952 39.30254509018036 +5.131526168812287 39.21675350701403 +5.2504246796831096 39.1309619238477 +5.369170563906945 39.045170340681366 +5.487763818700406 38.95937875751503 +5.606204441248349 38.8735871743487 +5.724492428703016 38.787795591182366 +5.842627778183168 38.70200400801603 +5.9606104867731995 38.6162124248497 +6.078440551522235 38.530420841683366 +6.196117969443211 38.44462925851703 +6.313642737511934 38.358837675350706 +6.431014852666161 38.27304609218437 +6.548234311804575 38.18725450901804 +6.665301111785833 38.101462925851706 +6.782215249427555 38.01567134268537 +6.898976721505291 37.92987975951904 +7.015585524751485 37.844088176352706 +7.132041655854405 37.75829659318637 +7.248345111457067 37.67250501002004 +7.36449588815612 37.58671342685371 +7.480493982500761 37.50092184368738 +7.596339390991535 37.415130260521046 +7.712032110079209 37.32933867735471 +7.827572136163578 37.24354709418838 +7.942959465592258 37.157755511022046 +8.058194094659454 37.07196392785571 +8.17327601960471 36.98617234468938 +8.288205236611635 36.900380761523046 +8.402981741806602 36.81458917835671 +8.517605531257422 36.728797595190386 +8.63207660097203 36.64300601202405 +8.74639494689706 36.55721442885772 +8.860560564916492 36.471422845691386 +8.97457345085021 36.38563126252505 +9.088433600452566 36.29983967935872 +9.202141009410902 36.214048096192386 +9.315695673344049 36.12825651302605 +9.429097587800802 36.04246492985972 +9.54234674825836 35.956673346693385 +9.655443150120737 35.87088176352706 +9.768386788717184 35.785090180360726 +9.881177659300489 35.69929859719439 +9.993815757045347 35.61350701402806 +10.106301077046643 35.527715430861726 +10.21863361431772 35.44192384769539 +10.330813363788613 35.35613226452906 +10.44284032030425 35.270340681362725 +10.554714478622623 35.18454909819639 +10.666435833412924 35.09875751503006 +10.778004379253638 35.01296593186373 +10.889420110630642 34.9271743486974 +11.000683021935192 34.841382765531066 +11.111793107461944 34.75559118236473 +11.222750361406915 34.6697995991984 +11.333554777865395 34.584008016032065 +11.44420635082984 34.49821643286573 +11.55470507418771 34.4124248496994 +11.665050941719285 34.326633266533065 +11.775243947095417 34.24084168336673 +11.885284083875261 34.155050100200405 +11.995171345503984 34.06925851703407 +12.104905725310358 33.98346693386774 +12.21448721650439 33.897675350701405 +12.32391581217487 33.81188376753507 +12.433191505286882 33.72609218436874 +12.542314288679256 33.640300601202405 +12.651284155062005 33.55450901803607 +12.760101097013676 33.46871743486974 +12.868765106978683 33.382925851703405 +12.977276177264566 33.29713426853708 +13.085634300039247 33.211342685370745 +13.193839467328155 33.12555110220441 +13.301891671011374 33.03975951903808 +13.409790902820703 32.953967935871745 +13.517537154336667 32.86817635270541 +13.625130416985478 32.78238476953908 +13.73257068203593 32.696593186372745 +13.839857940596254 32.61080160320641 +13.946992183610881 32.525010020040085 +14.053973401857217 32.43921843687375 +14.160801585942249 32.35342685370742 +14.267476726299185 32.267635270541085 +14.373998813183993 32.18184368737475 +14.480367836671876 32.09605210420842 +14.586583786653682 32.010260521042085 +14.692646652832249 31.924468937875755 +14.798556424718699 31.838677354709425 +14.904313091628635 31.75288577154309 +15.009916642678277 31.66709418837676 +15.115367066780538 31.581302605210425 +15.220664352641016 31.49551102204409 +15.325808488753914 31.40971943887776 +15.430799463397895 31.32392785571143 +15.535637264631829 31.238136272545095 +15.640321880290495 31.15234468937876 +15.744853297980185 31.066553106212428 +15.849231505074231 30.980761523046098 +15.953456488708454 30.894969939879765 +16.057528235776505 30.80917835671343 +16.161446732925146 30.723386773547098 +16.26521196654943 30.637595190380765 +16.368823922787787 30.551803607214435 +16.472282587517043 30.4660120240481 +16.575587946347294 30.380220440881768 +16.67873998461674 30.294428857715435 +16.781738687386387 30.2086372745491 +16.884584039434653 30.12284569138277 +16.9872760252519 30.037054108216438 +17.089814629034823 29.951262525050105 +17.192199834680743 29.86547094188377 +17.29443162578182 29.77967935871744 +17.39650998561914 29.693887775551108 +17.498434897156663 29.608096192384775 +17.6002063430351 29.52230460921844 +17.701824305565648 29.436513026052108 +17.8032887667236 29.350721442885778 +17.90459970814188 29.264929859719444 +18.005757111104373 29.17913827655311 +18.10676095653921 29.093346693386778 +18.20761122501187 29.007555110220444 +18.308307896718183 28.921763527054114 +18.408850951477188 28.83597194388778 +18.509240368723834 28.750180360721448 +18.609476127501576 28.664388777555114 +18.709558206454805 28.57859719438878 +18.809486583821133 28.49280561122245 +18.909261237423564 28.407014028056118 +19.008882144662454 28.321222444889784 +19.108349282507362 28.23543086172345 +19.20766262748873 28.149639278557117 +19.3068221556894 28.063847695390788 +19.405827842735977 27.978056112224454 +19.504679663790007 27.89226452905812 +19.60337759353899 27.806472945891787 +19.701921606187234 27.720681362725454 +19.800311675446505 27.634889779559124 +19.89854777452653 27.54909819639279 +19.99662987612527 27.463306613226457 +20.09455795241903 27.377515030060124 +20.19233197505239 27.29172344689379 +20.28995191512789 27.20593186372746 +20.387417743195584 27.120140280561127 +20.48472942924231 27.034348697394794 +20.581886942680807 26.94855711422846 +20.67889025233859 26.862765531062127 +20.775739326446615 26.776973947895797 +20.87243413262773 26.691182364729464 +20.968974637884873 26.60539078156313 +21.06536080858905 26.519599198396797 +21.161592610467075 26.433807615230464 +21.257670008589077 26.348016032064134 +21.353592967355738 26.2622244488978 +21.449361450485277 26.176432865731467 +21.544975421000203 26.090641282565134 +21.640434841213764 26.0048496993988 +21.735739672716164 25.91905811623247 +21.830889876360477 25.833266533066137 +21.925885412248277 25.747474949899804 +22.020726239715 25.66168336673347 +22.11541231731499 25.575891783567137 +22.209943602806266 25.490100200400807 +22.304320053134962 25.404308617234474 +22.398541624419458 25.31851703406814 +22.492608271934188 25.232725450901807 +22.58651995009313 25.146933867735473 +22.680276612432948 25.061142284569144 +22.773878211595804 24.97535070140281 +22.8673246993118 24.889559118236477 +22.96061602638109 24.803767535070143 +23.05375214265558 24.717975951903814 +23.14673299702033 24.63218436873748 +23.23955853737449 24.546392785571147 +23.33222871061189 24.460601202404813 +23.42474346260124 24.37480961923848 +23.517102738165875 24.28901803607215 +23.609306481063157 24.203226452905817 +23.701354633963366 24.117434869739483 +23.793247138428203 24.03164328657315 +23.884983934888858 23.945851703406817 +23.976564962623577 23.860060120240487 +24.067990159734816 23.774268537074153 +24.159259463125863 23.68847695390782 +24.250372808477017 23.602685370741487 +24.341330130221245 23.516893787575153 +24.432131361519335 23.431102204408823 +24.522776434234544 23.34531062124249 +24.613265278906667 23.259519038076157 +24.703597824725612 23.173727454909823 +24.793773999504396 23.08793587174349 +24.883793729651554 23.00214428857716 +24.973656940143005 22.916352705410826 +25.063363554493275 22.830561122244493 +25.152913494726135 22.74476953907816 +25.24230668134462 22.658977955911826 +25.331543033300385 22.573186372745496 +25.42062246796244 22.487394789579163 +25.509544901085167 22.40160320641283 +25.598310246775718 22.315811623246496 +25.686918417460642 22.230020040080163 +25.775369323851855 22.144228456913833 +25.86366287491185 22.0584368737475 +25.951798977818136 21.972645290581166 +26.03977753792696 21.886853707414833 +26.12759845873619 21.8010621242485 +26.215261641847437 21.71527054108217 +26.30276698692732 21.629478957915836 +26.39011439166791 21.543687374749503 +26.477303751746305 21.45789579158317 +26.56433496078331 21.372104208416836 +26.65120791030123 21.286312625250506 +26.737922489680756 21.200521042084173 +26.824478586116832 21.11472945891784 +26.91087608457364 21.028937875751506 +26.997114867738535 20.943146292585173 +27.083194815974984 20.857354709418843 +27.169115807274473 20.77156312625251 +27.25487771720731 20.685771543086176 +27.3404804188724 20.599979959919843 +27.42592378284583 20.51418837675351 +27.51120767712838 20.42839679358718 +27.596331967091807 20.342605210420846 +27.681296515423963 20.256813627254513 +27.76610118207267 20.17102204408818 +27.850745824188344 20.085230460921846 +27.935230296065317 19.999438877755516 +28.019554449081866 19.913647294589182 +28.103718131638857 19.82785571142285 +28.187721189097022 19.742064128256516 +28.27156346371283 19.656272545090182 +28.355244794572872 19.570480961923852 +28.438765017526798 19.48468937875752 +28.522123965118695 19.398897795591186 +28.605321466516916 19.313106212424852 +28.68835734744232 19.22731462925852 +28.771231430094858 19.14152304609219 +28.853943533078493 19.055731462925856 +28.936493471324376 18.969939879759522 +29.01888105601229 18.88414829659319 +29.101106094490266 18.79835671342686 +29.183168390192378 18.712565130260526 +29.265067742554585 18.626773547094192 +29.3468039469287 18.54098196392786 +29.42837679449431 18.455190380761525 +29.509786072168705 18.369398797595196 +29.591031562514694 18.283607214428862 +29.67211304364626 18.19781563126253 +29.75303028913206 18.112024048096195 +29.833783067896654 18.026232464929862 +29.91437114411944 17.940440881763532 +29.994794277131216 17.8546492985972 +30.075052221308336 17.768857715430865 +30.15514472596438 17.683066132264532 +30.235071535239282 17.5972745490982 +30.31483238798586 17.51148296593187 +30.394427017653687 17.425691382765535 +30.473861297823554 17.339899799599202 +30.553142947210066 17.25410821643287 +30.63227411605424 17.168316633266535 +30.711255243600842 17.082525050100205 +30.790086328535562 16.996733466933872 +30.86876736892077 16.91094188376754 +30.947298362180778 16.825150300601205 +31.025679305086957 16.73935871743487 +31.103910193742728 16.653567134268542 +31.181991023568454 16.56777555110221 +31.259921789286196 16.481983967935875 +31.3377024849044 16.39619238476954 +31.415333103702515 16.31040080160321 +31.49281363821552 16.22460921843688 +31.570144080218462 16.138817635270545 +31.647324420710902 16.05302605210421 +31.72435464990142 15.967234468937876 +31.8012347571921 15.881442885771545 +31.877964731163075 15.795651302605211 +31.954544559557153 15.70985971943888 +32.03097422926451 15.624068136272546 +32.10725372630753 15.538276553106213 +32.183383035825855 15.452484969939881 +32.25936214206152 15.366693386773548 +32.33519102834444 15.280901803607216 +32.41086967707813 15.195110220440883 +32.48639806972576 15.10931863727455 +32.561776186796536 15.023527054108218 +32.6370040078326 14.937735470941885 +32.71208151139635 14.851943887775553 +32.787008675058296 14.76615230460922 +32.86178547538556 14.680360721442888 +32.93641188793105 14.594569138276555 +33.01088788722337 14.508777555110221 +33.08521344675756 14.42298597194389 +33.159388538986796 14.337194388777556 +33.233413135315075 14.251402805611225 +33.307287206091 14.165611222444891 +33.381010720602816 14.079819639278558 +33.45458364707478 13.994028056112226 +33.52800595266495 13.908236472945893 +33.601277603464624 13.822444889779561 +33.67439856449945 13.736653306613228 +33.747368799732435 13.650861723446894 +33.82018827206902 13.565070140280563 +33.892856943364336 13.47927855711423 +33.965374774432874 13.393486973947898 +34.037741725060734 13.307695390781564 +34.109957754020705 13.22190380761523 +34.18202281909034 13.1361122244489 +34.253936877073365 13.050320641282566 +34.32569988382453 12.964529058116234 +34.39731179427836 12.8787374749499 +34.46877256248196 12.792945891783567 +34.540082141632226 12.707154308617236 +34.61124048411786 12.621362725450902 +34.682247541566454 12.53557114228457 +34.75310326489713 12.449779559118237 +34.82380760437902 12.363987975951906 +34.894360509696206 12.278196392785572 +34.96476193001937 12.192404809619239 +35.03501181408488 12.106613226452907 +35.105110110281686 12.020821643286574 +35.17505676674668 11.935030060120242 +35.244851731469126 11.849238476953909 +35.314494952404786 11.763446893787576 +35.383986377600486 11.677655310621244 +35.453325955329866 11.59186372745491 +35.52251363424108 11.506072144288579 +35.591549363517345 11.420280561122246 +35.66043309305126 11.334488977955912 +35.72916477363385 11.24869739478958 +35.797744357159374 11.162905811623247 +35.866171796847105 11.077114228456916 +35.93444704748115 10.991322645290582 +36.00257006566971 10.905531062124249 +36.070540810125095 10.819739478957917 +36.138359241966 10.733947895791584 +36.20602532504354 10.648156312625252 +36.27353902629283 10.562364729458919 +36.34090031611183 10.476573146292585 +36.40810916876941 10.390781563126254 +36.47516556284471 10.30498997995992 +36.54206948169996 10.219198396793589 +36.60882091398918 10.133406813627255 +36.67541985420521 10.047615230460922 +36.74186630326787 9.96182364729459 +36.80816026915608 9.876032064128257 +36.87430176758709 9.790240480961925 +36.9402908227461 9.704448897795592 +37.006127468069906 9.618657314629258 +37.07181174708833 9.532865731462927 +37.1373437143276 9.447074148296593 +37.20272916905411 9.361282565130262 +37.26797675933494 9.275490981963928 +37.33309012847224 9.189699398797597 +37.398070285292086 9.103907815631263 +37.46291773240203 9.01811623246493 +37.527633000717685 8.932324649298598 +37.592216651340415 8.846533066132265 +37.65666927756368 8.760741482965933 +37.72099150701718 8.6749498997996 +37.78518400395867 8.589158316633267 +37.84924747172412 8.503366733466935 +37.91318265534775 8.417575150300602 +37.97699034436413 8.33178356713427 +38.04067137580594 8.245991983967937 +38.10422663741147 8.160200400801603 +38.16765707105759 8.074408817635272 +38.23096367643484 7.988617234468938 +38.29414751498274 7.902825651302606 +38.35720971410489 7.817034068136273 +38.420151471685 7.731242484969941 +38.482974060926736 7.645450901803608 +38.54567883554216 7.559659318637275 +38.6082672353157 7.473867735470942 +38.67074079207265 7.38807615230461 +38.733101136083974 7.302284569138277 +38.795350002941625 7.216492985971945 +38.85748924094182 7.130701402805612 +38.91952081901682 7.044909819639279 +38.9814468352595 6.959118236472946 +39.04326952608887 6.873326653306614 +39.10499127610925 6.787535070140281 +39.166614628720424 6.701743486973949 +39.22814229754168 6.615951903807615 +39.28957717871848 6.530160320641283 +39.35092236418705 6.44436873747495 +39.41218115597956 6.358577154308618 +39.47335708166071 6.272785571142285 +39.53445391099549 6.186993987975953 +39.5954653781818 6.101202404809619 +39.65635724114077 6.015410821643287 +39.71710183675489 5.929619238476954 +39.77769779796367 5.843827655310622 +39.83814637917637 5.758036072144289 +39.89844965476232 5.672244488977956 +39.9586099077926 5.5864529058116235 +40.018629646867886 5.500661322645291 +40.07851162427874 5.414869739478958 +40.138258855612065 5.329078156312626 +40.19787464092837 5.2432865731462925 +40.25736258764703 5.15749498997996 +40.3167266352907 5.0717034068136275 +40.375971082255596 4.985911823647295 +40.43510061479213 4.9001202404809625 +40.49412033840004 4.814328657314629 +40.55303581186471 4.728537074148297 +40.61185308418661 4.642745490981964 +40.670578734684575 4.556953907815632 +40.72921991658635 4.471162324649299 +40.78778440445705 4.385370741482967 +40.846280645858656 4.299579158316633 +40.90471781768247 4.213787575150301 +40.96310588765215 4.127995991983968 +41.02145568155941 4.042204408817636 +41.07977895686838 3.9564128256513027 +41.13808848341046 3.87062124248497 +41.19639813199078 3.7848296593186372 +41.25472297184284 3.6990380761523047 +41.31307937800252 3.6132464929859722 +41.371485149830036 3.5274549098196393 +41.42995964209279 3.4416633266533068 +41.488523910239174 3.3558717434869743 +41.54720087174925 3.2700801603206413 +41.60601548575131 3.184288577154309 +41.66499495345326 3.0984969939879763 +41.72416894236683 3.0127054108216433 +41.783569837816025 2.926913827655311 +41.843233025838174 2.841122244488978 +41.9031972123299 2.7553306613226454 +41.963504784192 2.669539078156313 +42.024202219324486 2.58374749498998 +42.08534055366527 2.4979559118236474 +42.14697591511581 2.4121643286573144 +42.209170136236764 2.326372745490982 +42.271991460131844 2.2405811623246494 +42.335515357109664 2.1547895791583165 +42.39982547370629 2.068997995991984 +42.46501474071347 1.9832064128256515 +42.53118667332464 1.8974148296593187 +42.598456904837846 1.8116232464929862 +42.6669550061728 1.7258316633266535 +42.7368266576443 1.6400400801603208 +42.80823625822887 1.5542484969939883 +42.88137008274485 1.4684569138276555 +42.95644013153293 1.3826653306613228 +43.03368886422277 1.29687374749499 +43.11339507482985 1.2110821643286573 +43.195881258766626 1.1252905811623248 +43.2815229577175 1.039498997995992 +43.37076076920711 0.9537074148296594 +43.46411601392595 0.8679158316633268 +43.56221153586788 0.7821242484969941 +43.66579989869254 0.6963326653306614 +43.7758025933236 0.6105410821643287 +43.89336632865798 0.5247494989979961 +44.019947285520644 0.4389579158316634 +44.15744455614658 0.3531663326653307 +44.31040818247419 0.26737474949899803 +44.48735181159537 0.18158316633266536 +44.70389791126603 0.09579158316633267 +44.98408614316132 0.01 diff --git a/ELoss/E_vs_x_N-17.dat b/ELoss/E_vs_x_N-17.dat index eb7faac..669043c 100644 --- a/ELoss/E_vs_x_N-17.dat +++ b/ELoss/E_vs_x_N-17.dat @@ -1,501 +1,501 @@ Distance_cm Energy_MeV --0.0 80.0 -0.310447080374236 79.83969939879759 -0.6204578294907753 79.6793987975952 -0.930032231257981 79.51909819639279 -1.2391702696195288 79.35879759519038 -1.5478719285548268 79.19849699398797 -1.8561371920792706 79.03819639278558 -2.1639660442446633 78.87789579158317 -2.4713584691393633 78.71759519038076 -2.7783144508887063 78.55729458917835 -3.084833973655291 78.39699398797595 -3.390917021639267 78.23669338677355 -3.6965635790787554 78.07639278557114 -4.001773630250005 77.91609218436874 -4.306547159467815 77.75579158316633 -4.610884151085797 77.59549098196393 -4.914784589496799 77.43519038076153 -5.218248459133066 77.27488977955912 -5.521275744466659 77.11458917835671 -5.823866430009726 76.95428857715432 -6.126020500314924 76.79398797595191 -6.427737939975579 76.6336873747495 -6.72901873362611 76.47338677354709 -7.0298628659423 76.3130861723447 -7.330270321641717 76.15278557114229 -7.630241085483879 75.99248496993988 -7.929775142270676 75.83218436873747 -8.228872476846638 75.67188376753508 -8.527533074099361 75.51158316633267 -8.825756918959671 75.35128256513026 -9.123543996402042 75.19098196392785 -9.420894291444894 75.03068136272545 -9.717807789150884 74.87038076152305 -10.01428447462733 74.71008016032064 -10.31032433302637 74.54977955911824 -10.60592734954539 74.38947895791583 -10.901093509427275 74.22917835671343 -11.195822797960844 74.06887775551102 -11.490115200481 73.90857715430862 -11.783970702369155 73.74827655310621 -12.077389289053482 73.58797595190381 -12.37037094600934 73.4276753507014 -12.662915658759424 73.267374749499 -12.955023412874183 73.10707414829659 -13.246694193972075 72.9467735470942 -13.537927987719975 72.78647294589179 -13.828724779833328 72.62617234468938 -14.119084556076556 72.46587174348697 -14.40900730226331 72.30557114228458 -14.698493004256871 72.14527054108217 -14.987541647970296 71.98496993987976 -15.276153219366819 71.82466933867735 -15.564327704460112 71.66436873747494 -15.85206508931456 71.50406813627255 -16.13936536004565 71.34376753507014 -16.426228502820102 71.18346693386773 -16.712654503856257 71.02316633266533 -16.99864334942431 70.86286573146293 -17.28419502584668 70.70256513026052 -17.569309519498155 70.54226452905812 -17.853986816806238 70.38196392785571 -18.138226904251383 70.22166332665331 -18.422029768367356 70.0613627254509 -18.705395395741345 69.9010621242485 -18.98832377301432 69.74076152304609 -19.27081488688123 69.5804609218437 -19.552868724091358 69.42016032064129 -19.834485271448415 69.25985971943888 -20.11566451581088 69.09955911823647 -20.396406444092182 68.93925851703408 -20.676711043261037 68.77895791583167 -20.95657830034151 68.61865731462926 -21.23600820241335 68.45835671342685 -21.515000736612173 68.29805611222444 -21.793555890129642 68.13775551102205 -22.071673650213775 67.97745490981964 -22.349354004168987 67.81715430861723 -22.626596939356382 67.65685370741483 -22.903402443193883 67.49655310621243 -23.17977050315651 67.33625250501002 -23.455701106776406 67.17595190380761 -23.731194241643077 67.0156513026052 -24.00624989540352 66.85535070140281 -24.280868055762454 66.6950501002004 -24.555048710482314 66.534749498998 -24.828791847383474 66.37444889779559 -25.102097454344328 66.2141482965932 -25.37496551930148 66.05384769539079 -25.64739603024972 65.89354709418838 -25.919388975242196 65.73324649298597 -26.190944342390463 65.57294589178358 -26.462062119864633 65.41264529058117 -26.732742295893306 65.25234468937876 -27.002984858763703 65.09204408817635 -27.272789796821694 64.93174348697394 -27.54215709847179 64.77144288577155 -27.811086752177253 64.61114228456914 -28.079578746459955 64.45084168336673 -28.347633069900464 64.29054108216432 -28.615249711137974 64.13024048096193 -28.882428658870364 63.969939879759515 -29.14916990185396 63.80963927855711 -29.41547342890362 63.649338677354706 -29.681339228892636 63.4890380761523 -29.94676729075255 63.3287374749499 -30.21175760347315 63.16843687374749 -30.476310156102226 63.00813627254509 -30.740424937745534 62.84783567134268 -31.004101937566528 62.68753507014028 -31.267341144786275 62.52723446893787 -31.530142548683152 62.36693386773547 -31.79250613859272 62.20663326653306 -32.0544319039074 62.04633266533066 -32.31591983407629 61.88603206412825 -32.5769699186048 61.72573146292585 -32.83758214705445 61.56543086172344 -33.097756509042426 61.405130260521034 -33.35749299424132 61.24482965931863 -33.616791592378746 61.084529058116225 -33.875652293236904 60.92422845691382 -34.13407508665221 60.763927855711415 -34.39205996251479 60.603627254509014 -34.6496069107681 60.443326653306606 -34.906715921408306 60.283026052104205 -35.16338698448388 60.1227254509018 -35.41962009009493 59.962424849699396 -35.67541522839274 59.80212424849699 -35.93077238957902 59.641823647294586 -36.18569156390539 59.48152304609218 -36.44017274167259 59.32122244488978 -36.694215913229876 59.16092184368737 -36.947821068974164 59.00062124248497 -37.200988199349375 58.84032064128256 -37.4537172948455 58.68002004008016 -37.70600834599785 58.51971943887775 -37.95786134338611 58.35941883767535 -38.209276277633435 58.19911823647294 -38.46025313940548 58.03881763527053 -38.71079191940938 57.87851703406813 -38.960892608392754 57.718216432865724 -39.21055519714253 57.55791583166332 -39.45977967648391 57.397615230460914 -39.708566037279084 57.23731462925851 -39.9569142704261 57.077014028056105 -40.20482436685751 56.916713426853704 -40.452296317539115 56.756412825651296 -40.699330113468534 56.596112224448895 -40.94592574567383 56.43581162324649 -41.19208320521197 56.275511022044086 -41.43780248316737 56.11521042084168 -41.6830835706502 55.954909819639276 -41.92792645879486 55.79460921843687 -42.17233113875816 55.63430861723447 -42.41629760171765 55.47400801603206 -42.65982583886971 55.31370741482966 -42.9029158414277 55.15340681362725 -43.14556760061997 54.99310621242485 -43.38778110768789 54.83280561122244 -43.629556353883686 54.67250501002003 -43.87089333046829 54.51220440881763 -44.11179202870915 54.35190380761522 -44.3522524398778 54.19160320641282 -44.59227455524763 54.03130260521041 -44.831858366091225 53.87100200400801 -45.07100386367798 53.710701402805604 -45.30971103927132 53.5504008016032 -45.5479798841261 53.390100200400795 -45.785810389485675 53.229799599198394 -46.0232025465791 53.069498997995986 -46.260156346618054 52.909198396793585 -46.496671780793804 52.748897795591176 -46.73274884027397 52.588597194388775 -46.96838751619929 52.42829659318637 -47.20358779968015 52.267995991983966 -47.43834968179316 52.10769539078156 -47.672673153577485 51.94739478957916 -47.90655820603118 51.78709418837675 -48.14000483010728 51.62679358717435 -48.37301301670994 51.46649298597194 -48.60558275669026 51.30619238476954 -48.83771404084219 51.14589178356713 -49.06940685989809 50.98559118236472 -49.300661204524346 50.82529058116232 -49.53147706531679 50.66498997995991 -49.7618544327959 50.50468937875751 -49.991793297402 50.3443887775551 -50.22129364949015 50.1840881763527 -50.4503554793251 50.023787575150294 -50.67897877707583 49.86348697394789 -50.90716353281019 49.703186372745485 -51.13490973648916 49.542885771543084 -51.362217377961144 49.382585170340676 -51.58908644695586 49.222284569138274 -51.81551693307831 49.061983967935866 -52.04150882580233 48.901683366733465 -52.26706211446418 48.74138276553106 -52.49217678825573 48.581082164328656 -52.71685283621764 48.42078156312625 -52.94109024723215 48.26048096192385 -53.16488901001588 48.10018036072144 -53.388249113112195 47.93987975951904 -53.61117054488357 47.77957915831663 -53.833653293503524 47.61927855711422 -54.055697346948456 47.45897795591182 -54.27730269298925 47.29867735470941 -54.49846931918255 47.13837675350701 -54.719197212861864 46.9780761523046 -54.93948636112836 46.8177755511022 -55.15933675084147 46.65747494989979 -55.37874836860911 46.49717434869739 -55.59772120077776 46.336873747494984 -55.81625523342214 46.17657314629258 -56.03435045233469 46.016272545090175 -56.25200684301464 45.85597194388777 -56.46922439065691 45.695671342685365 -56.686003080140544 45.535370741482964 -56.90234289601696 45.375070140280556 -57.11824382249773 45.214769539078155 -57.33370584344215 45.05446893787575 -57.548728942344304 44.894168336673346 -57.763313102319955 44.73386773547094 -57.97745830609284 44.57356713426854 -58.1911645359808 44.41326653306613 -58.40443177388134 44.25296593186372 -58.61726000125688 44.09266533066132 -58.82964919911958 43.93236472945891 -59.041599348015666 43.77206412825651 -59.25311042800946 43.6117635270541 -59.464182418666745 43.4514629258517 -59.67481529903792 43.29116232464929 -59.88500904764041 43.13086172344689 -60.094763642440824 42.97056112224448 -60.30407906083641 42.81026052104208 -60.51295527963614 42.649959919839674 -60.72139227504114 42.48965931863727 -60.9293900226247 42.329358717434864 -61.13694849731155 42.16905811623246 -61.34406767335673 42.008757515030055 -61.550747524323704 41.848456913827654 -61.75698802306201 41.688156312625246 -61.9627891416841 41.527855711422845 -62.16815085154173 41.36755511022044 -62.3730731232015 41.207254509018036 -62.57755592641982 41.04695390781563 -62.78159923011709 40.88665330661322 -62.985203002351206 40.72635270541082 -63.188367210290316 40.56605210420841 -63.39109182018473 40.40575150300601 -63.593376797338124 40.2454509018036 -63.795222106077865 40.0851503006012 -63.99662770972455 39.92484969939879 -64.19759357056058 39.76454909819639 -64.39811964979799 39.60424849699398 -64.59820590754516 39.44394789579158 -64.79785230277284 39.28364729458917 -64.9970587932789 39.12334669338677 -65.19582533565236 38.96304609218436 -65.39415188523618 38.80274549098196 -65.59203839608917 38.642444889779554 -65.78948482094654 38.48214428857715 -65.98649111117969 38.321843687374745 -66.18305721675445 38.161543086172344 -66.37918308618846 38.001242484969936 -66.57486866650709 37.840941883767535 -66.77011390319821 37.68064128256513 -66.9649187401656 37.52034068136272 -67.15928311968105 37.36004008016032 -67.35320698233508 37.19973947895791 -67.54669026698615 37.03943887775551 -67.73973291070853 36.8791382765531 -67.93233484873852 36.7188376753507 -68.12449601441925 36.55853707414829 -68.31621633914371 36.39823647294589 -68.50749575229626 36.23793587174348 -68.69833418119235 36.07763527054108 -68.8887315510165 35.91733466933867 -69.07868778475847 35.75703406813627 -69.26820280314756 35.59673346693386 -69.45727652458493 35.43643286573146 -69.64590886507402 35.27613226452905 -69.8340997381488 35.11583166332665 -70.0218490548 34.955531062124244 -70.20915672339913 34.79523046092184 -70.3960226496202 34.634929859719435 -70.58244673635917 34.474629258517034 -70.768428883651 34.314328657314626 -70.95396898858414 34.15402805611222 -71.13906694521262 33.993727454909816 -71.3237226444654 33.83342685370741 -71.507935974053 33.67312625250501 -71.69170681837136 33.5128256513026 -71.8750350584028 33.3525250501002 -72.05792057161402 33.19222444889779 -72.2403632318509 33.03192384769539 -72.42236290923033 32.87162324649298 -72.60391947002849 32.71132264529058 -72.78503277656596 32.55102204408817 -72.96570268708916 32.39072144288577 -73.14592905564838 32.23042084168336 -73.32571173197175 32.07012024048096 -73.5050505613357 31.909819639278556 -73.68394538443115 31.74951903807615 -73.8623960372258 31.589218436873747 -74.04040235082199 31.428917835671342 -74.21796415131024 31.268617234468937 -74.39508125961827 31.108316633266533 -74.57175349135525 30.948016032064128 -74.74798065665131 30.787715430861724 -74.92376255999187 30.62741482965932 -75.09909900004692 30.467114228456914 -75.2739897694949 30.30681362725451 -75.448434654841 30.146513026052105 -75.62243343622976 29.9862124248497 -75.79598588725176 29.825911823647296 -75.96909177474416 29.66561122244489 -76.14175085858483 29.505310621242486 -76.31396289148007 29.345010020040082 -76.48572761874544 29.184709418837677 -76.65704477807957 29.02440881763527 -76.82791409933073 28.864108216432864 -76.99833530425592 28.70380761523046 -77.16830810627215 28.543507014028055 -77.33783221019964 28.38320641282565 -77.50690731199671 28.222905811623246 -77.67553309848611 28.06260521042084 -77.84370924707233 27.902304609218437 -78.01143542544968 27.742004008016032 -78.17871129130087 27.581703406813627 -78.34553649198561 27.421402805611223 -78.51191066421895 27.261102204408818 -78.67783343373898 27.100801603206413 -78.84330441496364 26.94050100200401 -79.00832321063588 26.780200400801604 -79.17288941145723 26.6198997995992 -79.33700259570901 26.459599198396795 -79.50066232886087 26.29929859719439 -79.66386816316616 26.138997995991986 -79.82661963724372 25.97869739478958 -79.9889162756455 25.818396793587176 -80.15075758840949 25.65809619238477 -80.31214307059764 25.497795591182363 -80.47307220181781 25.33749498997996 -80.63354444572956 25.177194388777554 -80.79355924953296 25.01689378757515 -80.95311604343986 24.856593186372745 -81.11221424012695 24.69629258517034 -81.27085323416992 24.535991983967936 -81.4290324014581 24.37569138276553 -81.58675109858872 24.215390781563126 -81.74400866224012 24.05509018036072 -81.9008044085231 23.894789579158317 -82.05713763230955 23.734488977955913 -82.21300760653759 23.574188376753508 -82.3684135814922 23.413887775551103 -82.5233547840606 23.2535871743487 -82.6778304169612 23.093286573146294 -82.8318396579453 22.93298597194389 -82.98538165897044 22.772685370741485 -83.13845554534433 22.61238476953908 -83.29106041483816 22.452084168336675 -83.44319533676827 22.29178356713427 -83.59485935104485 22.131482965931863 -83.74605146718648 21.971182364729458 -83.89677066329916 21.810881763527053 -84.04701588501854 21.65058116232465 -84.19678604441374 21.490280561122244 -84.3460800188517 21.32997995991984 -84.49489664981998 21.169679358717435 -84.64325344981809 21.00937875751503 -84.7911735435064 20.849078156312626 -84.93866215834781 20.68877755511022 -85.08571949361519 20.528476953907816 -85.23234553831101 20.36817635270541 -85.37854027698381 20.207875751503007 -85.52430368957071 20.047575150300602 -85.66963575123822 19.887274549098198 -85.81453643222159 19.726973947895793 -85.95900569766322 19.56667334669339 -86.1030435074504 19.406372745490984 -86.24664981605297 19.24607214428858 -86.38982457236148 19.085771543086175 -86.53256771952638 18.92547094188377 -86.67487919479912 18.76517034068136 -86.81675892937572 18.604869739478957 -86.95820684824407 18.444569138276552 -87.09922287003563 18.284268537074148 -87.2398069068829 18.123967935871743 -87.3799588642839 17.96366733466934 -87.5196786409751 17.803366733466934 -87.6589661288143 17.64306613226453 -87.79782121267546 17.482765531062125 -87.93624377035732 17.32246492985972 -88.0742336725081 17.162164328657315 -88.21179078256873 17.00186372745491 -88.34891495673746 16.841563126252506 -88.48560604395878 16.6812625250501 -88.62186388594026 16.520961923847697 -88.7576883172008 16.360661322645292 -88.89307916515493 16.200360721442888 -89.02803625023736 16.040060120240483 -89.16255938607337 15.879759519038075 -89.29664837970046 15.71945891783567 -89.4303030318479 15.559158316633265 -89.56352313728112 15.39885771543086 -89.69630848521874 15.238557114228456 -89.82865885983098 15.078256513026052 -89.96057404082914 14.917955911823647 -90.09205380415666 14.757655310621242 -90.22309792279367 14.597354709418838 -90.3537061676881 14.437054108216431 -90.48387830882788 14.276753507014027 -90.61361411647022 14.116452905811622 -90.74291336254586 13.956152304609217 -90.87177582225802 13.795851703406813 -91.0002012758978 13.635551102204408 -91.12818951090043 13.475250501002003 -91.25574032416915 13.314949899799599 -91.3828535246965 13.154649298597194 -91.50952893651608 12.99434869739479 -91.63576640202146 12.834048096192385 -91.76156578569284 12.673747494989978 -91.88692697827655 12.513446893787574 -92.01184990146761 12.35314629258517 -92.13633451315076 12.192845691382765 -92.26038081326182 12.03254509018036 -92.38398885033806 11.872244488977955 -92.50715872883376 11.71194388777555 -92.6298906172858 11.551643286573146 -92.75219234963754 11.391342685370741 -92.87408919710457 11.231042084168337 -92.99560099136974 11.07074148296593 -93.11673115955776 10.910440881763526 -93.23748207018804 10.750140280561121 -93.35785630294139 10.589839679358716 -93.47785667066756 10.429539078156312 -93.59748624376778 10.269238476953907 -93.71674837722104 10.108937875751502 -93.83564674055611 9.948637274549098 -93.95418535110745 9.788336673346693 -94.07236861093571 9.628036072144289 -94.19020134784105 9.467735470941884 -94.3076888609519 9.307434869739478 -94.4248369714335 9.147134268537073 -94.54165207893212 8.986833667334668 -94.65814122445144 8.826533066132264 -94.77431216045173 8.666232464929859 -94.89017342906953 8.505931863727454 -95.00573444948039 8.34563126252505 -95.12100561557092 8.185330661322645 -95.2359984052543 8.02503006012024 -95.35072550295808 7.864729458917835 -95.46520093704179 7.70442885771543 -95.5794402341697 7.544128256513026 -95.69342306494791 7.383827655310621 -95.80705194267318 7.2235270541082155 -95.9202562979675 7.063226452905811 -96.03303907254428 6.902925851703406 -96.14541000628606 6.742625250501002 -96.25738036100483 6.582324649298597 -96.36896311300154 6.422024048096192 -96.48017316978807 6.261723446893787 -96.59102761428412 6.101422845691382 -96.70154598033344 5.9411222444889775 -96.81175056402161 5.780821643286573 -96.92166677604182 5.620521042084168 -97.03132354127928 5.460220440881763 -97.14075375290562 5.299919839679358 -97.24999478964256 5.1396192384769535 -97.35908910653238 4.979318637274549 -97.46808491162416 4.819018036072144 -97.5770369435547 4.658717434869739 -97.68600736821337 4.498416833667334 -97.79506681671315 4.338116232464929 -97.90429559199127 4.177815631262525 -98.01378507785854 4.01751503006012 -98.1236393926483 3.857214428857715 -98.23397734038403 3.6969138276553104 -98.34493472641476 3.5366132264529053 -98.45666712290719 3.3763126252505007 -98.56935319404431 3.216012024048096 -98.68319872355592 3.055711422845691 -98.7984415316025 2.8954108216432863 -98.91535752888466 2.7351102204408813 -99.03426824033879 2.5748096192384766 -99.15555024976442 2.414509018036072 -99.2796471869406 2.254208416833667 -99.40708512656974 2.0939078156312623 -99.53849263621014 1.9336072144288576 -99.6746272686095 1.7733066132264528 -99.81641116290147 1.6130060120240481 -99.96497981252479 1.4527054108216433 -100.12175037127554 1.2924048096192384 -100.28851987621947 1.1321042084168336 -100.46761108259498 0.9718036072144288 -100.66209790124877 0.8115030060120241 -100.8761730207967 0.6512024048096192 -101.11579468387575 0.4909018036072144 -101.39334197127955 0.3306012024048096 -101.74118762747285 0.1703006012024048 -102.22847195660619 0.01 +-0.0 42.82 +0.11805569609279197 42.73420841683367 +0.23597861826135275 42.64841683366734 +0.353768762631506 42.56262525050101 +0.4714261252746706 42.476833667334674 +0.58895070220695 42.39104208416834 +0.7063424893881578 42.30525050100201 +0.8236014827208391 42.219458917835674 +0.9407276780492773 42.13366733466934 +1.0577210711584852 42.04787575150301 +1.1745816577731822 41.96208416833667 +1.291309433556755 41.87629258517034 +1.4079043941101932 41.790501002004014 +1.524366534971058 41.70470941883768 +1.6406958516123378 41.61891783567135 +1.756892339441384 41.533126252505014 +1.8729559937987823 41.44733466933868 +1.988886809957216 41.36154308617235 +2.104684783120313 41.27575150300601 +2.220349908421475 41.18995991983968 +2.335882180922691 41.10416833667335 +2.4512815956133274 41.01837675350701 +2.5665481474088985 40.93258517034069 +2.681681831149859 40.84679358717435 +2.796682641600287 40.76100200400802 +2.9115505734466396 40.67521042084169 +3.0262856212964464 40.58941883767535 +3.1408877796769885 40.50362725450902 +3.255357043033959 40.41783567134269 +3.369693405730105 40.33204408817635 +3.4838968620438453 40.24625250501002 +3.597967406167872 40.160460921843686 +3.7119050322077167 40.07466933867736 +3.8257097341803474 39.98887775551103 +3.939381506012644 39.90308617234469 +4.0529203415399495 39.81729458917836 +4.166326234504551 39.731503006012026 +4.279599178554146 39.64571142284569 +4.392739167240287 39.55991983967936 +4.505746194016802 39.474128256513026 +4.618620252238188 39.38833667334669 +4.731361335157986 39.30254509018036 +4.843969435927121 39.21675350701403 +4.956444547592253 39.1309619238477 +5.068786663094028 39.045170340681366 +5.180995775265379 38.95937875751503 +5.293071876829765 38.8735871743487 +5.405014960399384 38.787795591182366 +5.516825018473369 38.70200400801603 +5.628502043435949 38.6162124248497 +5.740046027554583 38.530420841683366 +5.851456962978067 38.44462925851703 +5.962734841734602 38.358837675350706 +6.073879655729882 38.27304609218437 +6.184891396745052 38.18725450901804 +6.295770056434738 38.101462925851706 +6.406515626324992 38.01567134268537 +6.517128097811218 37.92987975951904 +6.627607462156071 37.844088176352706 +6.737953710487311 37.75829659318637 +6.848166833795639 37.67250501002004 +6.958246822932479 37.58671342685371 +7.068193668607782 37.50092184368738 +7.178007361387692 37.415130260521046 +7.287687891692281 37.32933867735471 +7.397235249793193 37.24354709418838 +7.506649425811266 37.157755511022046 +7.615930409714118 37.07196392785571 +7.725078191313693 36.98617234468938 +7.8340927602637676 36.900380761523046 +7.942974106057426 36.81458917835671 +8.05172221802448 36.728797595190386 +8.160337085328896 36.64300601202405 +8.268818696966093 36.55721442885772 +8.377167041760291 36.471422845691386 +8.485382108361767 36.38563126252505 +8.593463885244082 36.29983967935872 +8.701412360701267 36.214048096192386 +8.80922752284496 36.12825651302605 +8.916909359601497 36.04246492985972 +9.02445785870897 35.956673346693385 +9.131873007714212 35.87088176352706 +9.239154793969796 35.785090180360726 +9.346303204630885 35.69929859719439 +9.453318226652137 35.61350701402806 +9.560199846784496 35.527715430861726 +9.666948051571964 35.44192384769539 +9.773562827348307 35.35613226452906 +9.880044160233712 35.270340681362725 +9.9863920361314 35.18454909819639 +10.092606440724175 35.09875751503006 +10.198687359470917 35.01296593186373 +10.304634777603056 34.9271743486974 +10.410448680120913 34.841382765531066 +10.516129051790067 34.75559118236473 +10.621675877137612 34.6697995991984 +10.72708914044838 34.584008016032065 +10.832368825761087 34.49821643286573 +10.937514916864433 34.4124248496994 +11.04252739729313 34.326633266533065 +11.14740625032387 34.24084168336673 +11.252151458971225 34.155050100200405 +11.356763005983513 34.06925851703407 +11.461240873838525 33.98346693386774 +11.565585044739267 33.897675350701405 +11.66979550060958 33.81188376753507 +11.773872223089715 33.72609218436874 +11.877815193531829 33.640300601202405 +11.981624392995403 33.55450901803607 +12.0852998022426 33.46871743486974 +12.188841401733539 33.382925851703405 +12.292249171621487 33.29713426853708 +12.395523091748018 33.211342685370745 +12.498663141638003 33.12555110220441 +12.601669300494613 33.03975951903808 +12.704541547194198 32.953967935871745 +12.807279860281078 32.86817635270541 +12.90988421796227 32.78238476953908 +13.012354598102117 32.696593186372745 +13.114690978216833 32.61080160320641 +13.21689333546895 32.525010020040085 +13.318961646661716 32.43921843687375 +13.420895888233328 32.35342685370742 +13.522696036251139 32.267635270541085 +13.624362066405741 32.18184368737475 +13.725893954004952 32.09605210420842 +13.827291673967714 32.010260521042085 +13.928555200817874 31.924468937875755 +14.029684508677894 31.838677354709425 +14.130679571262426 31.75288577154309 +14.23154036187179 31.66709418837676 +14.332266853385356 31.581302605210425 +14.432859018254812 31.49551102204409 +14.533316828497309 31.40971943887776 +14.633640255688526 31.32392785571143 +14.733829270955571 31.238136272545095 +14.833883844969804 31.15234468937876 +14.93380394793953 31.066553106212428 +15.033589549602565 30.980761523046098 +15.133240619218693 30.894969939879765 +15.232757125561973 30.80917835671343 +15.33213903691294 30.723386773547098 +15.431386321050674 30.637595190380765 +15.53049894524472 30.551803607214435 +15.629476876246903 30.4660120240481 +15.72832008028297 30.380220440881768 +15.827028523044115 30.294428857715435 +15.925602169678363 30.2086372745491 +16.024040984781795 30.12284569138277 +16.122344932389645 30.037054108216438 +16.220513975967222 29.951262525050105 +16.318548078400703 29.86547094188377 +16.416447201987754 29.77967935871744 +16.514211308428013 29.693887775551108 +16.611840358813378 29.608096192384775 +16.70933431361816 29.52230460921844 +16.806693132689066 29.436513026052108 +16.90391677523499 29.350721442885778 +17.001005199816667 29.264929859719444 +17.097958364336108 29.17913827655311 +17.19477622602589 29.093346693386778 +17.291458741438245 29.007555110220444 +17.388005866433957 28.921763527054114 +17.48441755617111 28.83597194388778 +17.580693765093564 28.750180360721448 +17.676834446919322 28.664388777555114 +17.772839554628636 28.57859719438878 +17.86870904045192 28.49280561122245 +17.9644428558575 28.407014028056118 +18.060040951539058 28.321222444889784 +18.155503277402968 28.23543086172345 +18.25082978255532 28.149639278557117 +18.34602041528879 28.063847695390788 +18.44107512306924 27.978056112224454 +18.53599385252208 27.89226452905812 +18.630776549418428 27.806472945891787 +18.72542315866101 27.720681362725454 +18.81993362426979 27.634889779559124 +18.914307889367414 27.54909819639279 +19.008545896164307 27.463306613226457 +19.10264758594359 27.377515030060124 +19.19661289904569 27.29172344689379 +19.290441774852685 27.20593186372746 +19.3841341517724 27.120140280561127 +19.47768996722217 27.034348697394794 +19.571109157612355 26.94855711422846 +19.66439165832957 26.862765531062127 +19.757537403719585 26.776973947895797 +19.850546327069953 26.691182364729464 +19.9434183605923 26.60539078156313 +20.036153435404326 26.519599198396797 +20.12875148151147 26.433807615230464 +20.22121242778824 26.348016032064134 +20.31353620195927 26.2622244488978 +20.405722730579917 26.176432865731467 +20.49777193901664 26.090641282565134 +20.589683751426943 26.0048496993988 +20.681458090738985 25.91905811623247 +20.773094878630843 25.833266533066137 +20.864594035509352 25.747474949899804 +20.955955480488605 25.66168336673347 +21.047179131368047 25.575891783567137 +21.13826490461018 25.490100200400807 +21.22921271531786 25.404308617234474 +21.32002247721118 25.31851703406814 +21.410694102603916 25.232725450901807 +21.501227502379596 25.146933867735473 +21.59162258596707 25.061142284569144 +21.68187926131568 24.97535070140281 +21.77199743486994 24.889559118236477 +21.8619770115438 24.803767535070143 +21.95181789469437 24.717975951903814 +22.041519986095235 24.63218436873748 +22.131083185909223 24.546392785571147 +22.2205073926607 24.460601202404813 +22.309792503207348 24.37480961923848 +22.39893841271141 24.28901803607215 +22.487945014610435 24.203226452905817 +22.57681220058744 24.117434869739483 +22.66553986054054 24.03164328657315 +22.75412788255202 23.945851703406817 +22.84257615285681 23.860060120240487 +22.930884555810422 23.774268537074153 +23.019052973856205 23.68847695390782 +23.107081287492072 23.602685370741487 +23.194969375236553 23.516893787575153 +23.28271711359423 23.431102204408823 +23.370324377020548 23.34531062124249 +23.4577910378859 23.259519038076157 +23.54511696643911 23.173727454909823 +23.632302030770173 23.08793587174349 +23.719346096772327 23.00214428857716 +23.80624902810341 22.916352705410826 +23.89301068614644 22.830561122244493 +23.979630929969517 22.74476953907816 +24.06610961628491 22.658977955911826 +24.15244659940739 22.573186372745496 +24.238641731211803 22.487394789579163 +24.32469486108975 22.40160320641283 +24.410605835905535 22.315811623246496 +24.496374499951223 22.230020040080163 +24.58200069490085 22.144228456913833 +24.667484259763764 22.0584368737475 +24.752825030837062 21.972645290581166 +24.83802284165712 21.886853707414833 +24.923077522950216 21.8010621242485 +25.007988902582152 21.71527054108217 +25.09275680550698 21.629478957915836 +25.177381053714665 21.543687374749503 +25.261861466177805 21.45789579158317 +25.346197858797286 21.372104208416836 +25.430390044346886 21.286312625250506 +25.514437832416856 21.200521042084173 +25.598343372429525 21.11472945891784 +25.68211315407838 21.028937875751506 +25.76575174457193 20.943146292585173 +25.849259714096803 20.857354709418843 +25.932637309541242 20.77156312625251 +26.01588453001807 20.685771543086176 +26.09900137427483 20.599979959919843 +26.181987840686656 20.51418837675351 +26.26484392724908 20.42839679358718 +26.347569631570803 20.342605210420846 +26.4301649508664 20.256813627254513 +26.512629881948982 20.17102204408818 +26.594964421222837 20.085230460921846 +26.677168564676 19.999438877755516 +26.759242307872814 19.913647294589182 +26.841185645946457 19.82785571142285 +26.922998573591435 19.742064128256516 +27.004681085056077 19.656272545090182 +27.086233174135007 19.570480961923852 +27.167654834161624 19.48468937875752 +27.24894605800056 19.398897795591186 +27.33010683804017 19.313106212424852 +27.411137166185043 19.22731462925852 +27.492037033848522 19.14152304609219 +27.57280643194529 19.055731462925856 +27.653445350883956 18.969939879759522 +27.733953780559755 18.88414829659319 +27.814331710347265 18.79835671342686 +27.89457912909325 18.712565130260526 +27.97469602510955 18.626773547094192 +28.054682386166128 18.54098196392786 +28.134538199484183 18.455190380761525 +28.214263451729447 18.369398797595196 +28.293858129005628 18.283607214428862 +28.37332221684798 18.19781563126253 +28.45265570021712 18.112024048096195 +28.531858563493024 18.026232464929862 +28.61093079046924 17.940440881763532 +28.68987236434739 17.8546492985972 +28.76868326773191 17.768857715430865 +28.847363482625095 17.683066132264532 +28.92591299042251 17.5972745490982 +29.004331771908685 17.51148296593187 +29.08261980725326 17.425691382765535 +29.16077707600749 17.339899799599202 +29.238803557101242 17.25410821643287 +29.316699228840438 17.168316633266535 +29.394464068905037 17.082525050100205 +29.472098054347583 16.996733466933872 +29.54960116159232 16.91094188376754 +29.626973366434985 16.825150300601205 +29.70421464404327 16.73935871743487 +29.781324968958035 16.653567134268542 +29.85830431509532 16.56777555110221 +29.935152655749167 16.481983967935875 +30.011869963595387 16.39619238476954 +30.08845621069628 16.31040080160321 +30.16491136850636 16.22460921843688 +30.241235407879245 16.138817635270545 +30.317428299075633 16.05302605210421 +30.393490011772595 15.967234468937876 +30.469420515074162 15.881442885771545 +30.545219777523357 15.795651302605211 +30.62088776711569 15.70985971943888 +30.696424451314304 15.624068136272546 +30.77182979706679 15.538276553106213 +30.84710377082382 15.452484969939881 +30.922246338559727 15.366693386773548 +30.99725746579507 15.280901803607216 +31.07213711762144 15.195110220440883 +31.146885258728513 15.10931863727455 +31.22150185343356 15.023527054108218 +31.295986865713566 14.937735470941885 +31.37034025924007 14.851943887775553 +31.444561997416937 14.76615230460922 +31.518652043421206 14.680360721442888 +31.592610360247228 14.594569138276555 +31.66643691075425 14.508777555110221 +31.7401316577177 14.42298597194389 +31.813694563884358 14.337194388777556 +31.887125592031623 14.251402805611225 +31.960424705031212 14.165611222444891 +32.03359186591741 14.079819639278558 +32.10662703796027 13.994028056112226 +32.179530184744 13.908236472945893 +32.25230127025078 13.822444889779561 +32.324940258950456 13.736653306613228 +32.39744711589633 13.650861723446894 +32.469821806827504 13.565070140280563 +32.54206429827807 13.47927855711423 +32.614174557693595 13.393486973947898 +32.68615255355532 13.307695390781564 +32.75799825551248 13.22190380761523 +32.82971163452334 13.1361122244489 +32.901292663005215 13.050320641282566 +32.972741314994266 12.964529058116234 +33.04405756631545 12.8787374749499 +33.11524139476328 12.792945891783567 +33.186292780293996 12.707154308617236 +33.25721170522985 12.621362725450902 +33.327998154476134 12.53557114228457 +33.39865211575176 12.449779559118237 +33.46917357983412 12.363987975951906 +33.53956254081907 12.278196392785572 +33.60981899639692 12.192404809619239 +33.67994294814533 12.106613226452907 +33.74993440184014 12.020821643286574 +33.81979336778508 11.935030060120242 +33.88951986116157 11.849238476953909 +33.959113902399636 11.763446893787576 +34.02857551757133 11.677655310621244 +34.0979047388078 11.59186372745491 +34.167101604741454 11.506072144288579 +34.23616802338304 11.420280561122246 +34.30511127877856 11.334488977955912 +34.37393720363919 11.24869739478958 +34.44264671487873 11.162905811623247 +34.51124046849496 11.077114228456916 +34.57971882202794 10.991322645290582 +34.648082149486754 10.905531062124249 +34.71633084224327 10.819739478957917 +34.78446530997601 10.733947895791584 +34.852485981666845 10.648156312625252 +34.920393306653935 10.562364729458919 +34.98818775574391 10.476573146292585 +35.055869822386974 10.390781563126254 +35.123440023918626 10.30498997995992 +35.1908989028719 10.219198396793589 +35.25824702836442 10.133406813627255 +35.32548499756465 10.047615230460922 +35.39261343724221 9.96182364729459 +35.4596330054073 9.876032064128257 +35.52654439304457 9.790240480961925 +35.593348325947325 9.704448897795592 +35.660045566658106 9.618657314629258 +35.72663691652224 9.532865731462927 +35.79312321786134 9.447074148296593 +35.859505356274234 9.361282565130262 +35.92578426307318 9.275490981963928 +35.991960917864084 9.189699398797597 +36.05803635127954 9.103907815631263 +36.124011647874646 9.01811623246493 +36.18988794919584 8.932324649298598 +36.25566645703399 8.846533066132265 +36.32134843687352 8.760741482965933 +36.38693522155054 8.6749498997996 +36.45242821513342 8.589158316633267 +36.5178288970408 8.503366733466935 +36.58313882641247 8.417575150300602 +36.64835964675037 8.33178356713427 +36.7134930908477 8.245991983967937 +36.77854098602582 8.160200400801603 +36.84350525970006 8.074408817635272 +36.908387945297086 7.988617234468938 +36.97319118854833 7.902825651302606 +37.03791725418594 7.817034068136273 +37.102568533069814 7.731242484969941 +37.16714754977645 7.645450901803608 +37.23165697068315 7.559659318637275 +37.296099612583546 7.473867735470942 +37.3604588968638 7.38807615230461 +37.42470431162644 7.302284569138277 +37.4888218028035 7.216492985971945 +37.55280990353413 7.130701402805612 +37.61666939937871 7.044909819639279 +37.68040171662944 6.959118236472946 +37.74400839984512 6.873326653306614 +37.80749111965541 6.787535070140281 +37.87085168106867 6.701743486973949 +37.93409203231815 6.615951903807615 +37.99721427428394 6.530160320641283 +38.06022067053137 6.44436873747495 +38.123113658009466 6.358577154308618 +38.18589585845703 6.272785571142285 +38.24857009056772 6.186993987975953 +38.311139382969756 6.101202404809619 +38.37360698808099 6.015410821643287 +38.435976396905026 5.929619238476954 +38.49825135484027 5.843827655310622 +38.56043587857996 5.758036072144289 +38.62253427418869 5.672244488977956 +38.68455115644868 5.5864529058116235 +38.7464914695781 5.500661322645291 +38.80836050943348 5.414869739478958 +38.870163947319234 5.329078156312626 +38.93190785553981 5.2432865731462925 +38.99359873484339 5.15749498997996 +39.055243543921726 5.0717034068136275 +39.116849731147916 4.985911823647295 +39.17842526875315 4.9001202404809625 +39.239978689665676 4.814328657314629 +39.30151912725952 4.728537074148297 +39.3630563582887 4.642745490981964 +39.42460084931409 4.556953907815632 +39.4861638069658 4.471162324649299 +39.54775723242492 4.385370741482967 +39.60939398055459 4.299579158316633 +39.67108782416329 4.213787575150301 +39.732853523943845 4.127995991983968 +39.794706904700924 4.042204408817636 +39.85666493855927 3.9564128256513027 +39.918745835936875 3.87062124248497 +39.98096914517284 3.7848296593186372 +40.04335586182223 3.6990380761523047 +40.10592854877206 3.6132464929859722 +40.16871146849785 3.5274549098196393 +40.231730728972465 3.4416633266533068 +40.29501444496455 3.3558717434869743 +40.35859291672775 3.2700801603206413 +40.422498828393174 3.184288577154309 +40.486767468744674 3.0984969939879763 +40.551436977491726 3.0127054108216433 +40.616548620672134 2.926913827655311 +40.68214709943468 2.841122244488978 +40.74828089719228 2.7553306613226454 +40.81500267102752 2.669539078156313 +40.88236969431013 2.58374749498998 +40.95044435879488 2.4979559118236474 +41.01929474606661 2.4121643286573144 +41.08899528016014 2.326372745490982 +41.15962747560352 2.2405811623246494 +41.231280798138314 2.1547895791583165 +41.304053659125664 2.068997995991984 +41.378054569371855 1.9832064128256515 +41.453403484095276 1.8974148296593187 +41.530233378407864 1.8116232464929862 +41.608692102544246 1.7258316633266535 +41.68894457889956 1.6400400801603208 +41.77117541979916 1.5542484969939883 +41.855592067342734 1.4684569138276555 +41.94242858685405 1.3826653306613228 +42.031950286696514 1.29687374749499 +42.12445939442829 1.2110821643286573 +42.220302100117344 1.1252905811623248 +42.31987739427955 1.039498997995992 +42.423648300291056 0.9537074148296594 +42.53215636317217 0.8679158316633268 +42.646040668433116 0.7821242484969941 +42.76606333791741 0.6963326653306614 +42.89314460389153 0.6105410821643287 +43.02841265887782 0.5247494989979961 +43.17327756020756 0.4389579158316634 +43.330829725471716 0.3531663326653307 +43.5065660928558 0.26737474949899803 +43.70860855237439 0.18158316633266536 +43.95235326660341 0.09579158316633267 +44.25102111629877 0.01 diff --git a/ELoss/E_vs_x_Ne-18.dat b/ELoss/E_vs_x_Ne-18.dat new file mode 100644 index 0000000..cb9f4de --- /dev/null +++ b/ELoss/E_vs_x_Ne-18.dat @@ -0,0 +1,501 @@ +Distance_cm Energy_MeV +-0.0 80.0 +0.16713525463015264 79.83969939879759 +0.33406473805004305 79.6793987975952 +0.5007884444544184 79.51909819639279 +0.6673063680316278 79.35879759519038 +0.8336185029635974 79.19849699398797 +0.9997248434257149 79.03819639278558 +1.1656253835868005 78.87789579158317 +1.3313201176089287 78.71759519038076 +1.4968090396473945 78.55729458917835 +1.662092143850605 78.39699398797595 +1.8271694243599683 78.23669338677355 +1.9920408753098522 78.07639278557114 +2.1567064908273967 77.91609218436874 +2.3211662650324683 77.75579158316633 +2.485420192037525 77.59549098196393 +2.6494682659475672 77.43519038076153 +2.813310480859938 77.27488977955912 +2.97694683086427 77.11458917835671 +3.1403773100423384 76.95428857715432 +3.303601912468001 76.79398797595191 +3.466620632206991 76.6336873747495 +3.6294334633168486 76.47338677354709 +3.7920403998467678 76.3130861723447 +3.9544414358375213 76.15278557114229 +4.116636565321244 75.99248496993988 +4.278625782321351 75.83218436873747 +4.440409080852373 75.67188376753508 +4.601986454919871 75.51158316633267 +4.763357898520203 75.35128256513026 +4.924523405640434 75.19098196392785 +5.085482970258169 75.03068136272545 +5.246236586341382 74.87038076152305 +5.406784247848318 74.71008016032064 +5.567125948727237 74.54977955911824 +5.727261682916312 74.38947895791583 +5.8871914443434274 74.22917835671343 +6.0469152269260595 74.06887775551102 +6.206433024571017 73.90857715430862 +6.365744831174311 73.74827655310621 +6.524850640620946 73.58797595190381 +6.683750446784783 73.4276753507014 +6.842444243528258 73.267374749499 +7.0009320247022435 73.10707414829659 +7.1592137841458126 72.9467735470942 +7.31728951568609 72.78647294589179 +7.475159213137954 72.62617234468938 +7.632822870303875 72.46587174348697 +7.790280480973665 72.30557114228458 +7.947532038924305 72.14527054108217 +8.10457753791963 71.98496993987976 +8.261416971710144 71.82466933867735 +8.418050334032769 71.66436873747494 +8.574477618610576 71.50406813627255 +8.73069881915259 71.34376753507014 +8.886713929353444 71.18346693386773 +9.042522942893175 71.02316633266533 +9.198125853436917 70.86286573146293 +9.353522654634682 70.70256513026052 +9.508713340120996 70.54226452905812 +9.66369790351466 70.38196392785571 +9.818476338418428 70.22166332665331 +9.973048638418758 70.0613627254509 +10.127414797085422 69.9010621242485 +10.28157480797124 69.74076152304609 +10.43552866461173 69.5804609218437 +10.589276360524838 69.42016032064129 +10.742817889210507 69.25985971943888 +10.896153244150392 69.09955911823647 +11.049282418807483 68.93925851703408 +11.202205406625788 68.77895791583167 +11.354922201029893 68.61865731462926 +11.50743279542463 68.45835671342685 +11.659737183194682 68.29805611222444 +11.811835357704176 68.13775551102205 +11.96372731229633 67.97745490981964 +12.115413040292964 67.81715430861723 +12.266892534994133 67.65685370741483 +12.418165789677664 67.49655310621243 +12.569232797598776 67.33625250501002 +12.72009355198955 67.17595190380761 +12.870748046058516 67.0156513026052 +13.021196272990165 66.85535070140281 +13.17143822594452 66.6950501002004 +13.321473898056563 66.534749498998 +13.471303282435787 66.37444889779559 +13.62092637216566 66.2141482965932 +13.770343160303154 66.05384769539079 +13.919553639878124 65.89354709418838 +14.068557803892825 65.73324649298597 +14.217355645321316 65.57294589178358 +14.365947157108948 65.41264529058117 +14.514332332171698 65.25234468937876 +14.662511163395624 65.09204408817635 +14.810483643636248 64.93174348697394 +14.958249765717909 64.77144288577155 +15.105809522433189 64.61114228456914 +15.253162906542174 64.45084168336673 +15.400309910771849 64.29054108216432 +15.547250527815386 64.13024048096193 +15.69398475033151 63.969939879759515 +15.840512570943682 63.80963927855711 +15.986833982239466 63.649338677354706 +16.132948976769757 63.4890380761523 +16.278857547048005 63.3287374749499 +16.42455968554947 63.16843687374749 +16.570055384710404 63.00813627254509 +16.71534463692726 62.84783567134268 +16.86042743455584 62.68753507014028 +17.005303769910494 62.52723446893787 +17.149973635263198 62.36693386773547 +17.294437022842725 62.20663326653306 +17.438693924833697 62.04633266533066 +17.5827443333757 61.88603206412825 +17.726588240562307 61.72573146292585 +17.870225638440154 61.56543086172344 +18.013656519007906 61.405130260521034 +18.156880874215283 61.24482965931863 +18.299898695962035 61.084529058116225 +18.442709976096864 60.92422845691382 +18.58531470641638 60.763927855711415 +18.72771287866396 60.603627254509014 +18.869904484528682 60.443326653306606 +19.01188951564411 60.283026052104205 +19.15366796358719 60.1227254509018 +19.295239819876986 59.962424849699396 +19.436605075973514 59.80212424849699 +19.57776372327643 59.641823647294586 +19.718715753123806 59.48152304609218 +19.859461156790765 59.32122244488978 +19.999999925488197 59.16092184368737 +20.140332050361344 59.00062124248497 +20.28045752248844 58.84032064128256 +20.420376332879247 58.68002004008016 +20.560088472473627 58.51971943887775 +20.69959393214001 58.35941883767535 +20.83889270267392 58.19911823647294 +20.977984774796354 58.03881763527053 +21.116870139152223 57.87851703406813 +21.25554878630873 57.718216432865724 +21.394020706753658 57.55791583166332 +21.53228589089372 57.397615230460914 +21.670344329052767 57.23731462925851 +21.808196011470045 57.077014028056105 +21.945840928298335 56.916713426853704 +22.083279069602135 56.756412825651296 +22.220510425355712 56.596112224448895 +22.357534985441198 56.43581162324649 +22.49435273964655 56.275511022044086 +22.630963677663576 56.11521042084168 +22.767367789085785 55.954909819639276 +22.90356506340633 55.79460921843687 +23.03955549001576 55.63430861723447 +23.175339058199874 55.47400801603206 +23.310915757137373 55.31370741482966 +23.446285575897594 55.15340681362725 +23.58144850343808 54.99310621242485 +23.716404528602197 54.83280561122244 +23.851153640116614 54.67250501002003 +23.985695826588763 54.51220440881763 +24.12003107650427 54.35190380761522 +24.254159378224237 54.19160320641282 +24.388080719982597 54.03130260521041 +24.52179508988326 53.87100200400801 +24.655302475897322 53.710701402805604 +24.788602865860124 53.5504008016032 +24.921696247468297 53.390100200400795 +25.0545826082767 53.229799599198394 +25.187261935695332 53.069498997995986 +25.319734216986117 52.909198396793585 +25.45199943925968 52.748897795591176 +25.58405758947199 52.588597194388775 +25.715908654420975 52.42829659318637 +25.847552620743016 52.267995991983966 +25.978989474909405 52.10769539078156 +26.11021920322267 51.94739478957916 +26.24124179181288 51.78709418837675 +26.372057226633782 51.62679358717435 +26.502665493458952 51.46649298597194 +26.633066577877752 51.30619238476954 +26.763260465291278 51.14589178356713 +26.893247140908148 50.98559118236472 +27.023026589740244 50.82529058116232 +27.152598796598348 50.66498997995991 +27.281963746087623 50.50468937875751 +27.411121422603074 50.3443887775551 +27.540071810324818 50.1840881763527 +27.668814893213327 50.023787575150294 +27.797350655004475 49.86348697394789 +27.92567907920456 49.703186372745485 +28.05380014908511 49.542885771543084 +28.18171384767767 49.382585170340676 +28.309420157768376 49.222284569138274 +28.43691906189247 49.061983967935866 +28.564210542328635 48.901683366733465 +28.69129458109325 48.74138276553106 +28.818171159934437 48.581082164328656 +28.944840260326067 48.42078156312625 +29.071301863461503 48.26048096192385 +29.197555950247317 48.10018036072144 +29.323602501296758 47.93987975951904 +29.449441496923143 47.77957915831663 +29.57507291713302 47.61927855711422 +29.700496741619226 47.45897795591182 +29.82571294975377 47.29867735470941 +29.95072152058049 47.13837675350701 +30.075522432807634 46.9780761523046 +30.20011566480015 46.8177755511022 +30.32450119457191 46.65747494989979 +30.448678999777613 46.49717434869739 +30.57264905770465 46.336873747494984 +30.696411345264607 46.17657314629258 +30.819965838984718 46.016272545090175 +30.94331251499898 45.85597194388777 +31.066451349039156 45.695671342685365 +31.189382316425476 45.535370741482964 +31.3121053920572 45.375070140280556 +31.43462055040285 45.214769539078155 +31.556927765490304 45.05446893787575 +31.679027010896565 44.894168336673346 +31.800918259737347 44.73386773547094 +31.92260148465635 44.57356713426854 +32.04407665781432 44.41326653306613 +32.16534375087781 44.25296593186372 +32.28640273500765 44.09266533066132 +32.4072535808472 43.93236472945891 +32.52789625851022 43.77206412825651 +32.64833073756853 43.6117635270541 +32.76855698703928 43.4514629258517 +32.888574975372 43.29116232464929 +33.00838467043521 43.13086172344689 +33.12798603950282 42.97056112224448 +33.247379049240074 42.81026052104208 +33.36656366568924 42.649959919839674 +33.48553985425489 42.48965931863727 +33.60430757968881 42.329358717434864 +33.72286680607454 42.16905811623246 +33.841217496811566 42.008757515030055 +33.959359614598995 41.848456913827654 +34.07729312141897 41.688156312625246 +34.19501797851952 41.527855711422845 +34.31253414639713 41.36755511022044 +34.42984158477867 41.207254509018036 +34.54694025260308 41.04695390781563 +34.66383010800239 40.88665330661322 +34.78051110828239 40.72635270541082 +34.896983209902785 40.56605210420841 +35.013246368456734 40.40575150300601 +35.12930053865002 40.2454509018036 +35.245145674279534 40.0851503006012 +35.360781728211336 39.92484969939879 +35.47620865235801 39.76454909819639 +35.59142639765558 39.60424849699398 +35.706434914039676 39.44394789579158 +35.82123415042119 39.28364729458917 +35.93582405466124 39.12334669338677 +36.050204573545514 38.96304609218436 +36.16437565275787 38.80274549098196 +36.27833723685334 38.642444889779554 +36.39208926923033 38.48214428857715 +36.50563169210217 38.321843687374745 +36.61896444646781 38.161543086172344 +36.73208747208184 38.001242484969936 +36.84500070742364 37.840941883767535 +36.95770408966574 37.68064128256513 +37.070197554641325 37.52034068136272 +37.18248103681086 37.36004008016032 +37.294554469227855 37.19973947895791 +37.40641778350363 37.03943887775551 +37.51807090977123 36.8791382765531 +37.62951377664829 36.7188376753507 +37.740746311198905 36.55853707414829 +37.851768438894496 36.39823647294589 +37.96258008357358 36.23793587174348 +38.07318116740042 36.07763527054108 +38.18357161082264 35.91733466933867 +38.29375133252751 35.75703406813627 +38.40372024939723 35.59673346693386 +38.51347827646278 35.43643286573146 +38.623025326856656 35.27613226452905 +38.732361311764194 35.11583166332665 +38.841486140373604 34.955531062124244 +38.950399719824546 34.79523046092184 +39.05910195515535 34.634929859719435 +39.16759274924867 34.474629258517034 +39.275872002775714 34.314328657314626 +39.3839396141388 34.15402805611222 +39.491795479412346 33.993727454909816 +39.59943949228223 33.83342685370741 +39.70687154398332 33.67312625250501 +39.81409152323534 33.5128256513026 +39.92109931617683 33.3525250501002 +40.02789480629727 33.19222444889779 +40.13447787436722 33.03192384769539 +40.24084839836649 32.87162324649298 +40.34700625341019 32.71132264529058 +40.45295131167273 32.55102204408817 +40.55868344230947 32.39072144288577 +40.6642025113763 32.23042084168336 +40.76950838174669 32.07012024048096 +40.87460091302642 31.909819639278556 +40.97947996146579 31.74951903807615 +41.08414537986927 31.589218436873747 +41.18859701750243 31.428917835671342 +41.292834719996215 31.268617234468937 +41.396858329248325 31.108316633266533 +41.50066768332168 30.948016032064128 +41.604262616339874 30.787715430861724 +41.707642958379516 30.62741482965932 +41.81080853535932 30.467114228456914 +41.91375916892593 30.30681362725451 +42.01649467633623 30.146513026052105 +42.11901487033617 29.9862124248497 +42.22131955903592 29.825911823647296 +42.32340854578118 29.66561122244489 +42.42528162902067 29.505310621242486 +42.52693860216951 29.345010020040082 +42.628379253468445 29.184709418837677 +42.72960336583878 29.02440881763527 +42.830610716732785 28.864108216432864 +42.9314010779796 28.70380761523046 +43.03197421562627 28.543507014028055 +43.13232988977391 28.38320641282565 +43.23246785440877 28.222905811623246 +43.33238785722805 28.06260521042084 +43.43208963946024 27.902304609218437 +43.53157293567985 27.742004008016032 +43.63083747361634 27.581703406813627 +43.72988297395699 27.421402805611223 +43.828709150143595 27.261102204408818 +43.927315708162666 27.100801603206413 +44.02570234632903 26.94050100200401 +44.123868755062475 26.780200400801604 +44.221814616657355 26.6198997995992 +44.31953960504472 26.459599198396795 +44.4170433855469 26.29929859719439 +44.514325614624134 26.138997995991986 +44.611385939613086 25.97869739478958 +44.70822399845685 25.818396793587176 +44.804839419426244 25.65809619238477 +44.901231820832024 25.497795591182363 +44.99740081072773 25.33749498997996 +45.09334598660287 25.177194388777554 +45.189066935065966 25.01689378757515 +45.28456323151731 24.856593186372745 +45.37983443981088 24.69629258517034 +45.47488011190513 24.535991983967936 +45.56969978750229 24.37569138276553 +45.66429299367563 24.215390781563126 +45.75865924448445 24.05509018036072 +45.85279804057618 23.894789579158317 +45.94670886877532 23.734488977955913 +46.04039120165851 23.574188376753508 +46.13384449711553 23.413887775551103 +46.227068197895484 23.2535871743487 +46.320061731137756 23.093286573146294 +46.41282450788721 22.93298597194389 +46.50535592259301 22.772685370741485 +46.59765535259054 22.61238476953908 +46.68972828343022 22.452084168336675 +46.781591894934756 22.29178356713427 +46.873258589408394 22.131482965931863 +46.964729457462816 21.971182364729458 +47.05600463342356 21.810881763527053 +47.14708425265793 21.65058116232465 +47.2379684515823 21.490280561122244 +47.32865736767084 21.32997995991984 +47.41915113946602 21.169679358717435 +47.509449906590845 21.00937875751503 +47.59955380976344 20.849078156312626 +47.689462990813915 20.68877755511022 +47.779177592704045 20.528476953907816 +47.86869775955004 20.36817635270541 +47.95802363664879 20.207875751503007 +48.04715537050797 20.047575150300602 +48.13609310888054 19.887274549098198 +48.22483700080402 19.726973947895793 +48.31338719664531 19.56667334669339 +48.4017438481514 19.406372745490984 +48.48990710850689 19.24607214428858 +48.577877132398925 19.085771543086175 +48.665654076090405 18.92547094188377 +48.7532380975024 18.76517034068136 +48.84062935630667 18.604869739478957 +48.927828014029565 18.444569138276552 +49.01483423416827 18.284268537074148 +49.10164818232089 18.123967935871743 +49.18827002633179 17.96366733466934 +49.27469993645373 17.803366733466934 +49.36093808552857 17.64306613226453 +49.446984649188565 17.482765531062125 +49.53283980608019 17.32246492985972 +49.618503738113006 17.162164328657315 +49.70397663073599 17.00186372745491 +49.78925867324422 16.841563126252506 +49.87435005911901 16.6812625250501 +49.95925098640481 16.520961923847697 +50.043961658126726 16.360661322645292 +50.12848228275268 16.200360721442888 +50.212813074704805 16.040060120240483 +50.29695425492496 15.879759519038075 +50.380906051499935 15.71945891783567 +50.46466870035226 15.559158316633265 +50.54824244600334 15.39885771543086 +50.63162754241612 15.238557114228456 +50.714824253925336 15.078256513026052 +50.797832856264215 14.917955911823647 +50.88065363769728 14.757655310621242 +50.96328690027001 14.597354709418838 +51.04573296118719 14.437054108216431 +51.127992154332844 14.276753507014027 +51.21006483194625 14.116452905811622 +51.29195136646969 13.956152304609217 +51.373652152585464 13.795851703406813 +51.45516760946135 13.635551102204408 +51.53649818322577 13.475250501002003 +51.617644349696086 13.314949899799599 +51.69860661738592 13.154649298597194 +51.77938553082005 12.99434869739479 +51.859997377378576 12.834048096192385 +51.940470139903056 12.673747494989978 +52.020817710343564 12.513446893787574 +52.101043131727295 12.35314629258517 +52.18114871928441 12.192845691382765 +52.26113691037399 12.03254509018036 +52.341010274159125 11.872244488977955 +52.42077152218919 11.71194388777555 +52.50042351998291 11.551643286573146 +52.57996929971611 11.391342685370741 +52.659412074129975 11.231042084168337 +52.73875525178858 11.07074148296593 +52.8180024538295 10.910440881763526 +52.89715753236777 10.750140280561121 +52.9762245907325 10.589839679358716 +53.05520800573647 10.429539078156312 +53.1341124522032 10.269238476953907 +53.21294293000315 10.108937875751502 +53.291704793881614 9.948637274549098 +53.37040378639588 9.788336673346693 +53.44904607431941 9.628036072144289 +53.527638288916116 9.467735470941884 +53.60618757054037 9.307434869739478 +53.68470161807782 9.147134268537073 +53.76318874381103 8.986833667334668 +53.841657934373046 8.826533066132264 +53.920118918543174 8.666232464929859 +53.99858224274524 8.505931863727454 +54.07705935523119 8.34563126252505 +54.1555627000758 8.185330661322645 +54.23410582227539 8.02503006012024 +54.31270348543876 7.864729458917835 +54.39137180378868 7.70442885771543 +54.47012839046343 7.544128256513026 +54.54899252442924 7.383827655310621 +54.62798533869604 7.2235270541082155 +54.70713003298487 7.063226452905811 +54.78639077236708 6.902925851703406 +54.8656671752427 6.742625250501002 +54.94490987977844 6.582324649298597 +55.02413071285859 6.422024048096192 +55.10334829033552 6.261723446893787 +55.18258344018353 6.101422845691382 +55.26185947485493 5.9411222444889775 +55.34120250050544 5.780821643286573 +55.420641768914145 5.620521042084168 +55.50021007903053 5.460220440881763 +55.579944236439935 5.299919839679358 +55.65988558071507 5.1396192384769535 +55.74008059270073 4.979318637274549 +55.82058159637453 4.819018036072144 +55.90144757318593 4.658717434869739 +55.98274511089628 4.498416833667334 +56.06454951418689 4.338116232464929 +56.14694611102537 4.177815631262525 +56.23003179746609 4.01751503006012 +56.313916874873435 3.857214428857715 +56.39872724841941 3.6969138276553104 +56.48460707541674 3.5366132264529053 +56.57172197844761 3.3763126252505007 +56.66026297398809 3.216012024048096 +56.75045131617868 3.055711422845691 +56.842544523282776 2.8954108216432863 +56.936843949828585 2.7351102204408813 +57.03370440364967 2.5748096192384766 +57.133546504624924 2.414509018036072 +57.23687277369295 2.254208416833667 +57.344288880173586 2.0939078156312623 +57.4565321519782 1.9336072144288576 +57.57451052057477 1.7733066132264528 +57.69935680289275 1.6130060120240481 +57.832506115155404 1.4527054108216433 +57.975809220578434 1.2924048096192384 +58.13170362498847 1.1321042084168336 +58.30348118145025 0.9718036072144288 +58.49572437647055 0.8115030060120241 +58.715052349153865 0.6512024048096192 +58.97145636708799 0.4909018036072144 +59.276292273963826 0.3306012024048096 +59.64908939725906 0.1703006012024048 +60.08468062361392 0.01 diff --git a/ELoss/E_vs_x_alpha.dat b/ELoss/E_vs_x_alpha.dat index d357192..431dae8 100644 --- a/ELoss/E_vs_x_alpha.dat +++ b/ELoss/E_vs_x_alpha.dat @@ -1,501 +1,501 @@ Distance_cm Energy_MeV --0.0 30.0 -1.9446985513936768 29.9398997995992 -3.8860830606275085 29.879799599198396 -5.824152858570795 29.819699398797596 -7.75890727441193 29.759599198396792 -9.690345635651623 29.699498997995992 -11.618467268097803 29.63939879759519 -13.543271495859377 29.579298597194388 -15.464757641339961 29.519198396793588 -17.382925025232733 29.459098196392784 -19.29777296651357 29.398997995991984 -21.209300782435854 29.338897795591183 -23.11750778852414 29.27879759519038 -25.022393298567813 29.21869739478958 -26.923956624615716 29.15859719438878 -28.82219707696977 29.098496993987975 -30.71711396417857 29.038396793587175 -32.60870659303208 28.97829659318637 -34.496974268554645 28.91819639278557 -36.38191629399964 28.85809619238477 -38.26353197084298 28.797995991983967 -40.14182059877668 28.737895791583167 -42.016781475703404 28.677795591182363 -43.88841389772941 28.617695390781563 -45.75671715915911 28.557595190380763 -47.62169055248852 28.49749498997996 -49.483333368398675 28.43739478957916 -51.341644895750136 28.377294589178355 -53.19662442157583 28.317194388777555 -55.048271231075525 28.257094188376755 -56.896584607609185 28.19699398797595 -58.74156383269031 28.13689378757515 -60.58320818598024 28.07679358717435 -62.42151694528146 28.016693386773547 -64.25648938653089 27.956593186372746 -66.08812478379424 27.896492985971943 -67.91642240925876 27.836392785571142 -69.74138153322748 27.776292585170342 -71.56300142411253 27.71619238476954 -73.38128134842827 27.656092184368738 -75.19622057078563 27.595991983967934 -77.00781835388466 27.535891783567134 -78.81607395850885 27.475791583166334 -80.62098664351821 27.41569138276553 -82.42255566584241 27.35559118236473 -84.22078028047498 27.295490981963926 -86.01565974046585 27.235390781563126 -87.80719329691551 27.175290581162326 -89.59538019896806 27.115190380761522 -91.38021969380424 27.05509018036072 -93.16171102663546 26.99498997995992 -94.93985344069677 26.934889779559118 -96.71464617723989 26.874789579158318 -98.48608847552723 26.814689378757514 -100.25417957282433 26.754589178356714 -102.0189187043939 26.694488977955913 -103.78030510348871 26.63438877755511 -105.53833800134457 26.57428857715431 -107.2930166271742 26.514188376753506 -109.04434020815968 26.454088176352705 -110.79230796944634 26.393987975951905 -112.53691913413562 26.3338877755511 -114.2781729232779 26.2737875751503 -116.01606855586644 26.213687374749497 -117.75060524882959 26.153587174348697 -119.48178221702469 26.093486973947897 -121.20959867323084 26.033386773547093 -122.93405382814169 25.973286573146293 -124.65514689035908 25.913186372745493 -126.37287706638585 25.85308617234469 -128.08724356061853 25.79298597194389 -129.7982455753411 25.732885771543085 -131.5058823107172 25.672785571142285 -133.21015296478376 25.612685370741485 -134.9110567334438 25.55258517034068 -136.60859281045896 25.49248496993988 -138.3027603874433 25.432384769539077 -139.99355865385533 25.372284569138277 -141.68098679699168 25.312184368737476 -143.3650440019797 25.252084168336673 -145.04572945177006 25.191983967935872 -146.72304232713032 25.13188376753507 -148.39698180663706 25.07178356713427 -150.06754706666933 25.011683366733468 -151.7347372814013 24.951583166332664 -153.3985516227947 24.891482965931864 -155.05898926059237 24.831382765531064 -156.71604936231054 24.77128256513026 -158.36973109323154 24.71118236472946 -160.02003361639714 24.651082164328656 -161.66695609260057 24.590981963927856 -163.31049768038 24.530881763527056 -164.95065753601085 24.470781563126252 -166.58743481349836 24.41068136272545 -168.2208286645709 24.350581162324648 -169.8508382386719 24.290480961923848 -171.4774626829533 24.230380761523048 -173.10070114226775 24.170280561122244 -174.7205527591612 24.110180360721444 -176.33701667386597 24.05008016032064 -177.9500920242929 23.98997995991984 -179.5597779460244 23.92987975951904 -181.16607357230689 23.869779559118236 -182.76897803404304 23.809679358717435 -184.36849045978505 23.749579158316635 -185.96460997572675 23.68947895791583 -187.55733570569603 23.62937875751503 -189.146666771148 23.569278557114227 -190.7326022911567 23.509178356713427 -192.31514138240834 23.449078156312627 -193.89428315919355 23.388977955911823 -195.47002673339952 23.328877755511023 -197.04237121450322 23.26877755511022 -198.611315709563 23.20867735470942 -200.1768593232118 23.14857715430862 -201.73900115764928 23.088476953907815 -203.29774031263398 23.028376753507015 -204.8530758854763 22.96827655310621 -206.40500697103033 22.90817635270541 -207.95353266168675 22.84807615230461 -209.49865204736494 22.787975951903807 -211.04036421550512 22.727875751503007 -212.5786682510612 22.667775551102206 -214.11356323649287 22.607675350701403 -215.64504825175757 22.547575150300602 -217.17312237430357 22.4874749498998 -218.69778467906144 22.427374749499 -220.21903423843693 22.367274549098198 -221.73687012230306 22.307174348697394 -223.25129139799196 22.247074148296594 -224.76229713028792 22.18697394789579 -226.26988638141867 22.12687374749499 -227.77405821104836 22.06677354709419 -229.2748116762694 22.006673346693386 -230.7721458315944 21.946573146292586 -232.2660597289489 21.886472945891782 -233.75655241766282 21.826372745490982 -235.2436229444632 21.76627254509018 -236.72727035346594 21.706172344689378 -238.20749368616777 21.646072144288578 -239.6842919814387 21.585971943887778 -241.15766427551384 21.525871743486974 -242.62760960198517 21.465771543086174 -244.0941269917941 21.40567134268537 -245.5572154732228 21.34557114228457 -247.01687407188675 21.28547094188377 -248.47310181072638 21.225370741482966 -249.9258977099987 21.165270541082165 -251.3752607872699 21.10517034068136 -252.82119005740637 21.04507014028056 -254.26368453256725 20.98496993987976 -255.70274322219592 20.924869739478957 -257.1383651330115 20.864769539078157 -258.57054926900133 20.804669338677353 -259.9992946314118 20.744569138276553 -261.42460021874075 20.684468937875753 -262.846465026729 20.62436873747495 -264.2648880483515 20.56426853707415 -265.6798682738096 20.50416833667335 -267.09140469052215 20.444068136272545 -268.4994962831171 20.383967935871745 -269.9041420334233 20.32386773547094 -271.30534092046133 20.26376753507014 -272.70309192043567 20.20366733466934 -274.09739400672555 20.143567134268537 -275.4882461498762 20.083466933867737 -276.87564731759073 20.023366733466933 -278.25959647472064 19.963266533066133 -279.6400925832575 19.903166332665332 -281.0171346023241 19.84306613226453 -282.39072148816507 19.78296593186373 -283.76085219413847 19.722865731462925 -285.12752567070635 19.662765531062124 -286.490740865426 19.602665330661324 -287.85049672294076 19.54256513026052 -289.20679218497037 19.48246492985972 -290.5596261903024 19.42236472945892 -291.90899767478265 19.362264529058116 -293.2549055713053 19.302164328657316 -294.59734880980426 19.242064128256512 -295.9363263172429 19.181963927855712 -297.271837017605 19.12186372745491 -298.60387983188474 19.061763527054108 -299.9324536780768 19.001663326653308 -301.2575574711669 18.941563126252504 -302.57919012312135 18.881462925851704 -303.8973505428775 18.821362725450903 -305.2120376363334 18.7612625250501 -306.52325030633705 18.7011623246493 -307.8309874526771 18.641062124248496 -309.13524797207134 18.580961923847696 -310.43603075815685 18.520861723446895 -311.733334701479 18.46076152304609 -313.02715868948053 18.40066132264529 -314.3175016064909 18.34056112224449 -315.60436233371524 18.280460921843687 -316.8877397492227 18.220360721442887 -318.1676327279357 18.160260521042083 -319.44404014161796 18.100160320641283 -320.71696085886316 18.040060120240483 -321.98639374508315 17.97995991983968 -323.2523376624957 17.91985971943888 -324.5147914701128 17.859759519038075 -325.7737540237279 17.799659318637275 -327.0292241759039 17.739559118236475 -328.2812007759601 17.67945891783567 -329.5296826699594 17.61935871743487 -330.77466870069566 17.559258517034067 -332.0161577076796 17.499158316633267 -333.2541485271261 17.439058116232466 -334.4886399919402 17.378957915831663 -335.7196309317028 17.318857715430862 -336.9471201726569 17.258757515030062 -338.1711065376932 17.19865731462926 -339.39158884633486 17.13855711422846 -340.60856591472304 17.078456913827655 -341.8220365556013 17.018356713426854 -343.03199957830026 16.958256513026054 -344.2384537887217 16.89815631262525 -345.4413979893221 16.83805611222445 -346.64083097909673 16.777955911823646 -347.836751553562 16.717855711422846 -349.0291585047391 16.657755511022046 -350.218050621136 16.597655310621242 -351.4034266877294 16.537555110220442 -352.58528548594694 16.477454909819638 -353.76362579364786 16.417354709418838 -354.9384463851044 16.357254509018038 -356.10974603098197 16.297154308617234 -357.27752349831906 16.237054108216434 -358.441777550507 16.176953907815633 -359.6025069472691 16.11685370741483 -360.75971044463853 16.05675350701403 -361.91338679493725 15.996653306613226 -363.06353474675285 15.936553106212424 -364.21015304491584 15.876452905811622 -365.35324043047586 15.81635270541082 -366.49279564067746 15.75625250501002 -367.6288174089355 15.696152304609218 -368.76130446480937 15.636052104208416 -369.89025553397704 15.575951903807614 -371.01566933820817 15.515851703406813 -372.13754459533675 15.455751503006011 -373.2558800192327 15.39565130260521 -374.37067431977283 15.335551102204407 -375.4819262028113 15.275450901803605 -376.58963437014864 15.215350701402805 -377.6937975195009 15.155250501002003 -378.7944143444667 15.095150300601201 -379.8914835344945 15.0350501002004 -380.98500377484817 14.974949899799599 -382.0749737465722 14.914849699398797 -383.16139212645527 14.854749498997995 -384.24425758699334 14.794649298597193 -385.3235687963515 14.734549098196391 -386.3993244183244 14.67444889779559 -387.47152311229627 14.614348697394789 -388.54016353319884 14.554248496993987 -389.6052443314687 14.494148296593185 -390.66676415300327 14.434048096192384 -391.72472163911516 14.373947895791582 -392.7791154264856 14.31384769539078 -393.829944147116 14.253747494989979 -394.87720642827867 14.193647294589177 -395.9209008924653 14.133547094188376 -396.9610261573346 14.073446893787574 -397.9975808356577 14.013346693386772 -399.03056353526233 13.95324649298597 -400.0599728589752 13.89314629258517 -401.08580740456233 13.833046092184368 -402.1080657646679 13.772945891783566 -403.126746526751 13.712845691382764 -404.14184827302023 13.652745490981962 -405.1533695803667 13.592645290581162 -406.1613090202947 13.53254509018036 -407.16566515884995 13.472444889779558 -408.1664365565459 13.412344689378756 -409.1636217682876 13.352244488977956 -410.1572193432934 13.292144288577154 -411.14722782501326 13.232044088176352 -412.1336457510458 13.17194388777555 -413.11647165305146 13.111843687374748 -414.0957040566638 13.051743486973947 -415.0713414813974 12.991643286573145 -416.04338244055276 12.931543086172343 -417.01182544111856 12.871442885771541 -417.9766689836704 12.811342685370741 -418.9379115622664 12.75124248496994 -419.89555166433934 12.691142284569137 -420.84958777058534 12.631042084168335 -421.80001835484893 12.570941883767533 -422.7468418840042 12.510841683366733 -423.6900568178323 12.450741482965931 -424.62966160889493 12.390641282565129 -425.56565470240315 12.330541082164327 -426.49803453608274 12.270440881763527 -427.4267995400344 12.210340681362725 -428.35194813658967 12.150240480961923 -429.27347874016203 12.090140280561121 -430.1913897570932 12.030040080160319 -431.1056795854939 11.969939879759519 -432.01634661508 11.909839679358717 -432.9233892270028 11.849739478957915 -433.82680579367354 11.789639278557113 -434.72659467858244 11.729539078156312 -435.62275423611163 11.66943887775551 -436.5152828113416 11.609338677354708 -437.4041787398513 11.549238476953906 -438.2894403475119 11.489138276553104 -439.17106595027286 11.429038076152304 -440.04905385394176 11.368937875751502 -440.923402353956 11.3088376753507 -441.7941097351473 11.248737474949898 -442.6611742714979 11.188637274549098 -443.52459422588896 11.128537074148296 -444.3843678498403 11.068436873747494 -445.2404933832412 11.008336673346692 -446.0929690540722 10.94823647294589 -446.94179307811805 10.88813627254509 -447.7869636586699 10.828036072144288 -448.62847898621834 10.767935871743486 -449.4663372381356 10.707835671342684 -450.3005365783471 10.647735470941884 -451.1310751569918 10.587635270541082 -451.95795111007124 10.52753507014028 -452.78116255908617 10.467434869739478 -453.6007076106614 10.407334669338676 -454.4165843561575 10.347234468937875 -455.22879087126944 10.287134268537073 -456.03732521561153 10.227034068136271 -456.84218543228815 10.16693386773547 -457.6433695474499 10.10683366733467 -458.44087556983453 10.046733466933867 -459.23470149029237 9.986633266533065 -460.02484528129526 9.926533066132263 -460.8113048964288 9.866432865731461 -461.5940782698673 9.806332665330661 -462.37316331583065 9.746232464929859 -463.1485579280226 9.686132264529057 -463.9202599790499 9.626032064128255 -464.68826731982114 9.565931863727455 -465.4525777789255 9.505831663326653 -466.21318916198965 9.44573146292585 -466.9700992510129 9.385631262525049 -467.7233058036794 9.325531062124247 -468.47280655264655 9.265430861723447 -469.2185992048092 9.205330661322645 -469.96068144053834 9.145230460921843 -470.6990509128935 9.08513026052104 -471.43370524680824 9.02503006012024 -472.16464203824773 8.964929859719438 -472.89185885333677 8.904829659318636 -473.61535322745823 8.844729458917834 -474.33512266432024 8.784629258517032 -475.05116463499144 8.724529058116232 -475.7634765769027 8.66442885771543 -476.4720558928146 8.604328657314628 -477.1768999497494 8.544228456913826 -477.8780060778861 8.484128256513026 -478.57537156941754 8.424028056112224 -479.26899367736803 8.363927855711422 -479.9588696143704 8.30382765531062 -480.6449965514011 8.243727454909818 -481.3273716164715 8.183627254509018 -482.00599189327454 8.123527054108216 -482.68085441978474 8.063426853707414 -483.3519561868102 8.003326653306614 -484.01929413649515 7.943226452905811 -484.68286516077103 7.88312625250501 -485.34266609975475 7.823026052104208 -485.99869374009216 7.7629258517034065 -486.6509448132453 7.7028256513026045 -487.2994159937209 7.6427254509018026 -487.9441038972393 7.5826252505010014 -488.58500507884065 7.5225250501001995 -489.2221160309277 7.462424849699398 -489.8554331812419 7.402324649298596 -490.484952890772 7.342224448897795 -491.110671451592 7.282124248496993 -491.73258508462703 7.222024048096192 -492.3506899373447 7.16192384769539 -492.9649820813698 7.101823647294588 -493.5754575100198 7.041723446893787 -494.18211213575995 6.981623246492985 -494.78494178757404 6.921523046092184 -495.38394220825035 6.861422845691382 -495.9791090515788 6.801322645290581 -496.57043787945867 6.741222444889779 -497.15792415891303 6.681122244488978 -497.7415632590089 6.621022044088176 -498.32135044768006 6.560921843687374 -498.8972808884509 6.500821643286573 -499.46934963705877 6.440721442885771 -500.0375516379729 6.3806212424849695 -500.6018817208082 6.3205210420841675 -501.1623345966315 6.260420841683366 -501.7189048541595 6.200320641282564 -502.27158695584546 6.140220440881763 -502.8203752338549 6.080120240480961 -503.3652638859279 6.020020040080159 -503.90624697112776 5.959919839679358 -504.44331840547557 5.899819639278556 -504.97647195747 5.839719438877755 -505.5057012434929 5.779619238476953 -506.0309997231017 5.719519038076152 -506.55236069420863 5.65941883767535 -507.06977728815053 5.599318637274549 -507.58324246464963 5.539218436873747 -508.0927490066702 5.479118236472945 -508.5982895151747 5.419018036072144 -509.0998564037844 5.358917835671342 -509.5974418933509 5.298817635270541 -510.0910380064469 5.238717434869739 -510.580636561784 5.178617234468938 -511.0662291685689 5.118517034068136 -511.5478072208109 5.0584168336673345 -512.0253618915941 4.9983166332665325 -512.4988841273324 4.9382164328657305 -512.9683646420252 4.878116232464929 -513.4337939115383 4.818016032064127 -513.8951621679325 4.757915831663326 -514.3524593938729 4.697815631262524 -514.8056753171481 4.637715430861723 -515.2547994053415 4.577615230460921 -515.6998208606931 4.51751503006012 -516.1407286152053 4.457414829659318 -516.5775113260445 4.397314629258516 -517.0101573713032 4.337214428857715 -517.4386548461923 4.277114228456913 -517.8629915597443 4.217014028056112 -518.283155032117 4.15691382765531 -518.6991324926005 4.096813627254509 -519.1109108784422 4.036713426853707 -519.5184768346184 3.9766132264529053 -519.9218167147003 3.9165130260521037 -520.3209165829779 3.856412825651302 -520.7157622180277 3.7963126252505006 -521.1063391179338 3.736212424849699 -521.492632507397 3.6761122244488975 -521.8746273470002 3.616012024048096 -522.2523083449257 3.555911823647294 -522.625659971467 3.4958116232464924 -522.9946664767127 3.435711422845691 -523.3593119118344 3.3756112224448893 -523.7195801544641 3.3155110220440878 -524.0754549387094 3.255410821643286 -524.4269198904268 3.1953106212424847 -524.7739585684527 3.135210420841683 -525.116554512587 3.0751102204408816 -525.4546912992241 3.0150100200400796 -525.7883526056503 2.954909819639278 -526.117522284157 2.8948096192384765 -526.4421844472781 2.834709418837675 -526.7623235656347 2.7746092184368734 -527.0779245800711 2.714509018036072 -527.3889730299994 2.6544088176352703 -527.6954552001333 2.5943086172344687 -527.9973582880958 2.534208416833667 -528.2946705957371 2.474108216432865 -528.5873817473994 2.4140080160320636 -528.875482938834 2.353907815631262 -529.1589672210148 2.2938076152304605 -529.4378298237193 2.233707414829659 -529.7120685244831 2.1736072144288574 -529.9816840693895 2.113507014028056 -530.2466806531692 2.0534068136272543 -530.5070664672812 1.9933066132264528 -530.7628543260761 1.9332064128256512 -531.0140623828552 1.8731062124248496 -531.2607149497184 1.813006012024048 -531.5028434376286 1.7529058116232463 -531.7404874362488 1.6928056112224448 -531.9736959570038 1.6327054108216432 -532.202528867726 1.5726052104208417 -532.427058553486 1.51250501002004 -532.6473718462431 1.4524048096192383 -532.8635722764056 1.3923046092184368 -533.0757827131423 1.3322044088176352 -533.2841484785854 1.2721042084168337 -533.4888410456759 1.212004008016032 -533.6900624628815 1.1519038076152304 -533.8880506950984 1.0918036072144288 -534.0830861342573 1.0317034068136273 -534.2754996238341 0.9716032064128256 -534.4656824714835 0.911503006012024 -534.654099113704 0.8514028056112224 -534.8413033789114 0.7913026052104208 -535.0279597264604 0.7312024048096192 -535.2148715173979 0.6711022044088176 -535.4030194802259 0.611002004008016 -535.5936154291123 0.5509018036072144 -535.7881797238598 0.4908016032064128 -535.9886576384534 0.4307014028056112 -536.1976039873823 0.3706012024048096 -536.4184988965045 0.310501002004008 -536.656348262986 0.2504008016032064 -536.9190144190857 0.1903006012024048 -537.2209352294095 0.1302004008016032 -537.5927100947087 0.0701002004008016 -538.2452746759647 0.01 +-0.0 80.0 +12.749770620274157 79.83969939879759 +25.478020644395357 79.6793987975952 +38.18473677772476 79.51909819639279 +50.86990573754322 79.35879759519038 +63.53351425306336 79.19849699398797 +76.17554906543448 79.03819639278558 +88.79599692775395 78.87789579158317 +101.39484460506704 78.71759519038076 +113.97207887437784 78.55729458917835 +126.52768652465411 78.39699398797595 +139.06165435683195 78.23669338677355 +151.5739691838256 78.07639278557114 +164.064617830526 77.91609218436874 +176.53358713381016 77.75579158316633 +188.9808639425433 77.59549098196393 +201.40643511758762 77.43519038076153 +213.8102875317996 77.27488977955912 +226.19240807003786 77.11458917835671 +238.55278362916454 76.95428857715432 +250.89140111805244 76.79398797595191 +263.20824745758114 76.6336873747495 +275.5033095806437 76.47338677354709 +287.7765744321465 76.3130861723447 +300.0280289690152 76.15278557114229 +312.2576601601896 75.99248496993988 +324.4654549866288 75.83218436873747 +336.65140044131 75.67188376753508 +348.8154835292328 75.51158316633267 +360.95769126741277 75.35128256513026 +373.0780106848855 75.19098196392785 +385.1764288227046 75.03068136272545 +397.2529327339396 74.87038076152305 +409.30750948367853 74.71008016032064 +421.3401461490201 74.54977955911824 +433.3508298190753 74.38947895791583 +445.33954759496294 74.22917835671343 +457.30628658981095 74.06887775551102 +469.25103392874655 73.90857715430862 +481.173777644168 73.74827655310621 +493.0745057894051 73.58797595190381 +504.95320552294834 73.4276753507014 +516.8098639948078 73.267374749499 +528.6444683670481 73.10707414829659 +540.4570058137814 72.9467735470942 +552.2474635211648 72.78647294589179 +564.0158286873883 72.62617234468938 +575.7620885226721 72.46587174348697 +587.4862302492573 72.30557114228458 +599.1882411014028 72.14527054108217 +610.8681083253705 71.98496993987976 +622.5258191794214 71.82466933867735 +634.161360933806 71.66436873747494 +645.7747208707533 71.50406813627255 +657.3658862844655 71.34376753507014 +668.934844481102 71.18346693386773 +680.4815827787722 71.02316633266533 +692.0060885075228 70.86286573146293 +703.5083490093302 70.70256513026052 +714.9883516380825 70.54226452905812 +726.4460837595707 70.38196392785571 +737.8815327514743 70.22166332665331 +749.2946860033516 70.0613627254509 +760.6855309166194 69.9010621242485 +772.0540549045437 69.74076152304609 +783.4002453922214 69.5804609218437 +794.7240898165701 69.42016032064129 +806.0255756263053 69.25985971943888 +817.304690281929 69.09955911823647 +828.5614212557098 68.93925851703408 +839.7957560316703 68.77895791583167 +851.0076821055627 68.61865731462926 +862.1971869848546 68.45835671342685 +873.3642581887088 68.29805611222444 +884.5088833382628 68.13775551102205 +895.6310508066263 67.97745490981964 +906.730748875815 67.81715430861723 +917.8079650867644 67.65685370741483 +928.8626869920927 67.49655310621243 +939.8949021560828 67.33625250501002 +950.9045981546539 67.17595190380761 +961.891762575342 67.0156513026052 +972.8563830172745 66.85535070140281 +983.7984470911497 66.6950501002004 +994.7179424192058 66.534749498998 +1005.6148566351997 66.37444889779559 +1016.4891773843786 66.2141482965932 +1027.3408923234579 66.05384769539079 +1038.169989120587 65.89354709418838 +1048.976455455326 65.73324649298597 +1059.7602790186154 65.57294589178358 +1070.5214475127502 65.41264529058117 +1081.259948651345 65.25234468937876 +1091.9757701593064 65.09204408817635 +1102.6688997728024 64.93174348697394 +1113.3393252392289 64.77144288577155 +1123.9870343171808 64.61114228456914 +1134.6120147764145 64.45084168336673 +1145.214254397817 64.29054108216432 +1155.79374097337 64.13024048096193 +1166.3504623061187 63.969939879759515 +1176.88440621013 63.80963927855711 +1187.3955605104602 63.649338677354706 +1197.8839130431186 63.4890380761523 +1208.349451655026 63.3287374749499 +1218.7921642039805 63.16843687374749 +1229.2120385586143 63.00813627254509 +1239.6090626197222 62.84783567134268 +1249.9832250156726 62.68753507014028 +1260.3345143509216 62.52723446893787 +1270.6629185105687 62.36693386773547 +1280.9684253904566 62.20663326653306 +1291.2510228971246 62.04633266533066 +1301.510698947768 61.88603206412825 +1311.7474414701896 61.72573146292585 +1321.9612384027575 61.56543086172344 +1332.1520776943557 61.405130260521034 +1342.319947304338 61.24482965931863 +1352.464835202482 61.084529058116225 +1362.5867293689362 60.92422845691382 +1372.6856177941752 60.763927855711415 +1382.7614884789443 60.603627254509014 +1392.8143294342133 60.443326653306606 +1402.8441286811199 60.283026052104205 +1412.8508742509205 60.1227254509018 +1422.8345541849335 59.962424849699396 +1432.7951565344881 59.80212424849699 +1442.7326693608647 59.641823647294586 +1452.6470807352434 59.48152304609218 +1462.5383787386425 59.32122244488978 +1472.406551461864 59.16092184368737 +1482.2515870054317 59.00062124248497 +1492.073473479535 58.84032064128256 +1501.8721990039649 58.68002004008016 +1511.6477517080548 58.51971943887775 +1521.400119730616 58.35941883767535 +1531.1292912198778 58.19911823647294 +1540.8352543334192 58.03881763527053 +1550.5179975719552 57.87851703406813 +1560.1775097389607 57.718216432865724 +1569.8137792991565 57.55791583166332 +1579.4267944188239 57.397615230460914 +1589.016543273396 57.23731462925851 +1598.583014047393 57.077014028056105 +1608.1261949343466 56.916713426853704 +1617.6460741367334 56.756412825651296 +1627.1426398658982 56.596112224448895 +1636.6158803419844 56.43581162324649 +1646.0657837938559 56.275511022044086 +1655.4923384590252 56.11521042084168 +1664.8955325835736 55.954909819639276 +1674.2753544220775 55.79460921843687 +1683.6317922375263 55.63430861723447 +1692.9648343012466 55.47400801603206 +1702.2744688928183 55.31370741482966 +1711.5606842999964 55.15340681362725 +1720.8234688186258 54.99310621242485 +1730.0628107525604 54.83280561122244 +1739.2786984135762 54.67250501002003 +1748.4711201212867 54.51220440881763 +1757.640064203058 54.35190380761522 +1766.785518993918 54.19160320641282 +1775.9074728364706 54.03130260521041 +1785.0059140808034 53.87100200400801 +1794.0808310844002 53.710701402805604 +1803.132212212045 53.5504008016032 +1812.1600460405348 53.390100200400795 +1821.1643214623273 53.229799599198394 +1830.1450271572353 53.069498997995986 +1839.102151488665 52.909198396793585 +1848.035682827018 52.748897795591176 +1856.9456095495907 52.588597194388775 +1865.8319200404776 52.42829659318637 +1874.6946026904689 52.267995991983966 +1883.5336458969518 52.10769539078156 +1892.349038063806 51.94739478957916 +1901.1407676013023 51.78709418837675 +1909.9088229259964 51.62679358717435 +1918.6531924606254 51.46649298597194 +1927.3738646339987 51.30619238476954 +1936.0708278808934 51.14589178356713 +1944.7440706419427 50.98559118236472 +1953.3935813635278 50.82529058116232 +1962.019348497668 50.66498997995991 +1970.6213605019057 50.50468937875751 +1979.199605839197 50.3443887775551 +1987.7540729777947 50.1840881763527 +1996.2847503911364 50.023787575150294 +2004.7916265577248 49.86348697394789 +2013.2746899610142 49.703186372745485 +2021.7339290892896 49.542885771543084 +2030.16933243555 49.382585170340676 +2038.5808886743296 49.222284569138274 +2046.9685866915843 49.061983967935866 +2055.332415185043 48.901683366733465 +2063.6723626341836 48.74138276553106 +2071.988417522787 48.581082164328656 +2080.280568338814 48.42078156312625 +2088.5488035742787 48.26048096192385 +2096.7931117251237 48.10018036072144 +2105.01348129109 47.93987975951904 +2113.2099007755933 47.77957915831663 +2121.3823586855897 47.61927855711422 +2129.530843531449 47.45897795591182 +2137.655343826825 47.29867735470941 +2145.7558480885214 47.13837675350701 +2153.832344836361 46.9780761523046 +2161.884822593052 46.8177755511022 +2169.9132698840567 46.65747494989979 +2177.9176752374538 46.49717434869739 +2185.8980271838063 46.336873747494984 +2193.8543142560243 46.17657314629258 +2201.78652498923 46.016272545090175 +2209.69464792062 45.85597194388777 +2217.5786715893296 45.695671342685365 +2225.438584536293 45.535370741482964 +2233.274375435914 45.375070140280556 +2241.08603307195 45.214769539078155 +2248.8735460814087 45.05446893787575 +2256.6369029784987 44.894168336673346 +2264.376092278781 44.73386773547094 +2272.0911024990255 44.57356713426854 +2279.7819221570758 44.41326653306613 +2287.448539771708 44.25296593186372 +2295.0909438624926 44.09266533066132 +2302.709122949657 43.93236472945891 +2310.3030655539433 43.77206412825651 +2317.8727601964742 43.6117635270541 +2325.4181953986126 43.4514629258517 +2332.9393596818245 43.29116232464929 +2340.4362415675414 43.13086172344689 +2347.9088295770257 42.97056112224448 +2355.3571122312323 42.81026052104208 +2362.781078050676 42.649959919839674 +2370.1807155552956 42.48965931863727 +2377.556013264322 42.329358717434864 +2384.9069596961453 42.16905811623246 +2392.233543368186 42.008757515030055 +2399.535752840079 41.848456913827654 +2406.8135767087715 41.688156312625246 +2414.0670035095595 41.527855711422845 +2421.29602172089 41.36755511022044 +2428.500619819898 41.207254509018036 +2435.6807862822893 41.04695390781563 +2442.8365095822205 40.88665330661322 +2449.9677781921846 40.72635270541082 +2457.074580582899 40.56605210420841 +2464.156905223194 40.40575150300601 +2471.2147405799074 40.2454509018036 +2478.248075117779 40.0851503006012 +2485.2570013314544 39.92484969939879 +2492.2416735199654 39.76454909819639 +2499.2021431176217 39.60424849699398 +2506.1384008950845 39.44394789579158 +2513.0504375717605 39.28364729458917 +2519.93824381534 39.12334669338677 +2526.8018102413307 38.96304609218436 +2533.6411274125876 38.80274549098196 +2540.4561858388374 38.642444889779554 +2547.2469759761957 38.48214428857715 +2554.013488226684 38.321843687374745 +2560.755712937734 38.161543086172344 +2567.473640401697 38.001242484969936 +2574.1672608553354 37.840941883767535 +2580.836564479321 37.68064128256513 +2587.4815413977194 37.52034068136272 +2594.102181677472 37.36004008016032 +2600.6984753278743 37.19973947895791 +2607.2704123000453 37.03943887775551 +2613.8179824863937 36.8791382765531 +2620.3411757200756 36.7188376753507 +2626.839981774452 36.55853707414829 +2633.3143903625332 36.39823647294589 +2639.764391136425 36.23793587174348 +2646.189973686761 36.07763527054108 +2652.591127542138 35.91733466933867 +2658.967842168537 35.75703406813627 +2665.3201069687448 35.59673346693386 +2671.6479112817638 35.43643286573146 +2677.951244382222 35.27613226452905 +2684.23009547977 35.11583166332665 +2690.484453718481 34.955531062124244 +2696.714308176231 34.79523046092184 +2702.9196478640893 34.634929859719435 +2709.1004617256867 34.474629258517034 +2715.256738636591 34.314328657314626 +2721.388467403665 34.15402805611222 +2727.4956367644263 33.993727454909816 +2733.5782353863965 33.83342685370741 +2739.6362518664437 33.67312625250501 +2745.6696747301216 33.5128256513026 +2751.6784924309977 33.3525250501002 +2757.6626933499797 33.19222444889779 +2763.6222657946296 33.03192384769539 +2769.5571979984775 32.87162324649298 +2775.467478120323 32.71132264529058 +2781.3530942435345 32.55102204408817 +2787.214034375337 32.39072144288577 +2793.0502864460987 32.23042084168336 +2798.8618383086045 32.07012024048096 +2804.6486777373293 31.909819639278556 +2810.410792427698 31.74951903807615 +2816.1481699953447 31.589218436873747 +2821.860797975359 31.428917835671342 +2827.5486638215316 31.268617234468937 +2833.2117549055865 31.108316633266533 +2838.850058516413 30.948016032064128 +2844.4635618592843 30.787715430861724 +2850.0522520550744 30.62741482965932 +2855.616116139465 30.467114228456914 +2861.1551410621455 30.30681362725451 +2866.669313686009 30.146513026052105 +2872.1586207863375 29.9862124248497 +2877.623049049982 29.825911823647296 +2883.0625850745364 29.66561122244489 +2888.477215367503 29.505310621242486 +2893.8669263454526 29.345010020040082 +2899.231704333176 29.184709418837677 +2904.5715355628295 29.02440881763527 +2909.886406173075 28.864108216432864 +2915.1763022082105 28.70380761523046 +2920.441209617296 28.543507014028055 +2925.681114253271 28.38320641282565 +2930.8960018720686 28.222905811623246 +2936.085858131717 28.06260521042084 +2941.250668591442 27.902304609218437 +2946.390418710755 27.742004008016032 +2951.505093848542 27.581703406813627 +2956.5946792621403 27.421402805611223 +2961.6591601064115 27.261102204408818 +2966.698521432809 27.100801603206413 +2971.712748188436 26.94050100200401 +2976.701825215101 26.780200400801604 +2981.665737248365 26.6198997995992 +2986.604468916583 26.459599198396795 +2991.518004739939 26.29929859719439 +2996.406329129477 26.138997995991986 +3001.2694263861245 25.97869739478958 +3006.1072806997095 25.818396793587176 +3010.919876147975 25.65809619238477 +3015.707196695585 25.497795591182363 +3020.469226193125 25.33749498997996 +3025.2059483760986 25.177194388777554 +3029.917346863919 25.01689378757515 +3034.603405158894 24.856593186372745 +3039.264106645205 24.69629258517034 +3043.899434587883 24.535991983967936 +3048.5093721317794 24.37569138276553 +3053.0939023005285 24.215390781563126 +3057.65300799551 24.05509018036072 +3062.1866719948034 23.894789579158317 +3066.6948769521377 23.734488977955913 +3071.1776053958374 23.574188376753508 +3075.634839727763 23.413887775551103 +3080.0665622222464 23.2535871743487 +3084.4727550250227 23.093286573146294 +3088.853400152155 22.93298597194389 +3093.208479488956 22.772685370741485 +3097.5379747889046 22.61238476953908 +3101.841867672556 22.452084168336675 +3106.1201396264464 22.29178356713427 +3110.3727720019947 22.131482965931863 +3114.599746014395 21.971182364729458 +3118.8010427415056 21.810881763527053 +3122.9766431227313 21.65058116232465 +3127.1265279578975 21.490280561122244 +3131.250677906118 21.32997995991984 +3135.349073484657 21.169679358717435 +3139.421695067778 21.00937875751503 +3143.468522885589 20.849078156312626 +3147.4895370228764 20.68877755511022 +3151.4847174179254 20.528476953907816 +3155.4540438613335 20.36817635270541 +3159.3974959948073 20.207875751503007 +3163.315053309949 20.047575150300602 +3167.2066951470233 19.887274549098198 +3171.072400693714 19.726973947895793 +3174.9121489838544 19.56667334669339 +3178.725918896143 19.406372745490984 +3182.513689152835 19.24607214428858 +3186.2754383184074 19.085771543086175 +3190.0111447981994 18.92547094188377 +3193.72078683702 18.76517034068136 +3197.4043425177224 18.604869739478957 +3201.0617897597426 18.444569138276552 +3204.6931063175957 18.284268537074148 +3208.298269779327 18.123967935871743 +3211.8772575649127 17.96366733466934 +3215.4300469246073 17.803366733466934 +3218.9566149372263 17.64306613226453 +3222.4569385083646 17.482765531062125 +3225.93099436854 17.32246492985972 +3229.378759071252 17.162164328657315 +3232.800208990955 17.00186372745491 +3236.195320320929 16.841563126252506 +3239.5640690710393 16.6812625250501 +3242.90643106538 16.520961923847697 +3246.2223819397805 16.360661322645292 +3249.5118971391694 16.200360721442888 +3252.774951914776 16.040060120240483 +3256.0115213211566 15.879759519038075 +3259.221580213026 15.71945891783567 +3262.4051032418774 15.559158316633265 +3265.56206485237 15.39885771543086 +3268.692439278458 15.238557114228456 +3271.79620053924 15.078256513026052 +3274.8733224345 14.917955911823647 +3277.923778539912 14.757655310621242 +3280.94754220187 14.597354709418838 +3283.9445865319176 14.437054108216431 +3286.914884400729 14.276753507014027 +3289.8584084316026 14.116452905811622 +3292.7751309934233 13.956152304609217 +3295.6650241930374 13.795851703406813 +3298.528059866989 13.635551102204408 +3301.3642095725554 13.475250501002003 +3304.173444578012 13.314949899799599 +3306.9557358520633 13.154649298597194 +3309.711054052347 12.99434869739479 +3312.439369512937 12.834048096192385 +3315.140652230741 12.673747494989978 +3317.8148718506945 12.513446893787574 +3320.461997649633 12.35314629258517 +3323.0819985187204 12.192845691382765 +3325.674842944298 12.03254509018036 +3328.240498986997 11.872244488977955 +3330.7789342589645 11.71194388777555 +3333.2901158990085 11.551643286573146 +3335.7740105454773 11.391342685370741 +3338.2305843066506 11.231042084168337 +3340.6598027284113 11.07074148296593 +3343.061630758937 10.910440881763526 +3345.436032710129 10.750140280561121 +3347.7829722154656 10.589839679358716 +3350.1024121839464 10.429539078156312 +3352.394314749745 10.269238476953907 +3354.6586412171687 10.108937875751502 +3356.8953520004757 9.948637274549098 +3359.104406558051 9.788336673346693 +3361.2857633204126 9.628036072144289 +3363.4393796114523 9.467735470941884 +3365.5652115622615 9.307434869739478 +3367.6632140168394 9.147134268537073 +3369.7333404289034 8.986833667334668 +3371.775542748954 8.826533066132264 +3373.789771300669 8.666232464929859 +3375.775974645608 8.505931863727454 +3377.7340994351302 8.34563126252505 +3379.66409024831 8.185330661322645 +3381.565889414554 8.02503006012024 +3383.4394368194853 7.864729458917835 +3385.2846696925694 7.70442885771543 +3387.1015223748154 7.544128256513026 +3388.8899260647736 7.383827655310621 +3390.6498085409285 7.2235270541082155 +3392.3810938584616 7.063226452905811 +3394.0837020182585 6.902925851703406 +3395.757548605939 6.742625250501002 +3397.4025443986297 6.582324649298597 +3399.018594937181 6.422024048096192 +3400.605600061568 6.261723446893787 +3402.1634534073473 6.101422845691382 +3403.692041861288 5.9411222444889775 +3405.1912449747183 5.780821643286573 +3406.660934333761 5.620521042084168 +3408.100972886596 5.460220440881763 +3409.5112142292487 5.299919839679358 +3410.8915018533357 5.1396192384769535 +3412.241668361891 4.979318637274549 +3413.5615346630957 4.819018036072144 +3414.8509091568226 4.658717434869739 +3416.109586935824 4.498416833667334 +3417.3373490328404 4.338116232464929 +3418.5339617577183 4.177815631262525 +3419.699176186092 4.01751503006012 +3420.832727884929 3.857214428857715 +3421.9343369927 3.6969138276553104 +3423.0037088163135 3.5366132264529053 +3424.04053516801 3.3763126252505007 +3425.0444967495923 3.216012024048096 +3426.015267008151 3.055711422845691 +3426.9525180502146 2.8954108216432863 +3427.8559294294455 2.7351102204408813 +3428.725200945019 2.5748096192384766 +3429.560071045515 2.414509018036072 +3430.3603430897388 2.254208416833667 +3431.125922669185 2.0939078156312623 +3431.8568706050264 1.9336072144288576 +3432.553478365708 1.7733066132264528 +3433.2163760038984 1.6130060120240481 +3433.846688257485 1.4527054108216433 +3434.4462642565104 1.2924048096192384 +3435.01802496156 1.1321042084168336 +3435.5665111296275 0.9718036072144288 +3436.098801816051 0.8115030060120241 +3436.6261906188843 0.6512024048096192 +3437.167629836447 0.4909018036072144 +3437.7582548561163 0.3306012024048096 +3438.481013284422 0.1703006012024048 +3440.140107711027 0.01 diff --git a/ELoss/E_vs_x_proton.dat b/ELoss/E_vs_x_proton.dat index 1be5562..9a5762c 100644 --- a/ELoss/E_vs_x_proton.dat +++ b/ELoss/E_vs_x_proton.dat @@ -1,501 +1,501 @@ Distance_cm Energy_MeV --0.0 16.5 -8.406046275708023 16.466953907815633 -16.797597620189666 16.43390781563126 -25.17464772858279 16.400861723446894 -33.53719030136126 16.367815631262527 -41.88521904432812 16.334769539078156 -50.21872766860876 16.30172344689379 -58.53770989065286 16.268677354709418 -66.84215943222279 16.23563126252505 -75.13207002039528 16.202585170340683 -83.40743538755417 16.169539078156312 -91.66824927138298 16.136492985971945 -99.91450541486613 16.103446893787574 -108.146197566277 16.070400801603206 -116.36331947917888 16.03735470941884 -124.56586491241717 16.004308617234468 -132.75382763011174 15.971262525050099 -140.92720140165557 15.93821643286573 -149.08598000170676 15.90517034068136 -157.23015721018294 15.872124248496993 -165.3597268122572 15.839078156312624 -173.47468267334733 15.806032064128255 -181.57501871569755 15.772985971943886 -189.66072878991827 15.739939879759518 -197.7318066984019 15.706893787575149 -205.7882462487757 15.67384769539078 -213.83004125389695 15.64080160320641 -221.85718553184546 15.607755511022043 -229.86967290591846 15.574709418837674 -237.86749720462112 15.541663326653305 -245.85065226166117 15.508617234468936 -253.8191319159414 15.475571142284567 -261.7729300115518 15.4425250501002 -269.71204039776404 15.40947895791583 -277.63645692902105 15.376432865731461 -285.54617346493126 15.343386773547092 -293.44118387025986 15.310340681362725 -301.32148201492254 15.277294589178355 -309.1870617739749 15.244248496993986 -317.0379170276059 15.211202404809617 -324.8740416611286 15.17815631262525 -332.69542956497355 15.14511022044088 -340.5020746346774 15.112064128256511 -348.293970770876 15.079018036072142 -356.0711118792947 15.045971943887775 -363.8334918707409 15.012925851703406 -371.5811046610924 14.979879759519036 -379.31394417128996 14.946833667334667 -387.03200432732666 14.9137875751503 -394.7352790602404 14.88074148296593 -402.42376230610137 14.847695390781562 -410.09744800600424 14.814649298597192 -417.7563301060574 14.781603206412823 -425.4004025573727 14.748557114228456 -433.0296593160565 14.715511022044087 -440.64409434319725 14.682464929859718 -448.24370160485626 14.649418837675348 -455.82847507205634 14.616372745490981 -463.39840875669114 14.583326653306612 -470.9534967481723 14.550280561122243 -478.49373310294425 14.517234468937874 -486.019111807819 14.484188376753506 -493.529626854532 14.451142284569137 -501.02527223972817 14.418096192384768 -508.5060419649517 14.385050100200399 -515.9719300366336 14.352004008016031 -523.4229304660813 14.318957915831662 -530.8590372694642 14.285911823647293 -538.280244467803 14.252865731462924 -545.6865460869565 14.219819639278557 -553.0779361576107 14.186773547094187 -560.454408715264 14.153727454909818 -567.8159578002154 14.120681362725449 -575.1625774575518 14.08763527054108 -582.4942617371341 14.054589178356713 -589.8110046935857 14.021543086172343 -597.1128003862765 13.988496993987974 -604.3996428793113 13.955450901803605 -611.6715262415149 13.922404809619238 -618.9284445464206 13.889358717434868 -626.1703918722528 13.8563126252505 -633.3973623019155 13.82326653306613 -640.6093499229763 13.790220440881763 -647.8063488276546 13.757174348697394 -654.9883531128036 13.724128256513024 -662.155356879898 13.691082164328655 -669.307354235018 13.658036072144288 -676.444339288836 13.624989979959919 -683.5663061565989 13.59194388777555 -690.6732489581148 13.55889779559118 -697.7651618177364 13.525851703406813 -704.8420388643473 13.492805611222444 -711.9038742313438 13.459759519038075 -718.9506621426533 13.426713426853706 -725.9823968275381 13.393667334669336 -732.9990724316182 13.360621242484969 -740.0006831020465 13.3275751503006 -746.9872229904272 13.29452905811623 -753.9586862528006 13.261482965931862 -760.9150670496258 13.228436873747494 -767.8563595457656 13.195390781563125 -774.7825579104671 13.162344689378756 -781.693656317346 13.129298597194387 -788.5896489443693 13.09625250501002 -795.4705299738387 13.06320641282565 -802.3362935923715 13.030160320641281 -809.1869339908841 12.997114228456912 -816.022445364574 12.964068136272545 -822.8428219129031 12.931022044088175 -829.6480578395776 12.897975951903806 -836.4381473525316 12.864929859719437 -843.2130846639077 12.83188376753507 -849.9728639900408 12.7988376753507 -856.7174795514369 12.765791583166331 -863.4469255727564 12.732745490981962 -870.1611962827949 12.699699398797593 -876.860285914464 12.666653306613226 -883.5441887047743 12.633607214428856 -890.2128988948136 12.600561122244487 -896.8664107297298 12.567515030060118 -903.5047184587105 12.53446893787575 -910.1278163349657 12.501422845691382 -916.7356986157059 12.468376753507012 -923.3283595621243 12.435330661322643 -929.9057934583401 12.402284569138276 -936.4679946395635 12.369238476953907 -943.0149574237822 12.336192384769538 -949.546676083572 12.303146292585168 -956.0631448954147 12.270100200400801 -962.5643581396795 12.237054108216432 -969.0503101006015 12.204008016032063 -975.5209950662623 12.170961923847694 -981.9764073285701 12.137915831663326 -988.4165411832399 12.104869739478957 -994.8413909297722 12.071823647294588 -1001.2509508714334 12.038777555110219 -1007.6452153152351 12.005731462925851 -1014.0241785719155 11.972685370741482 -1020.3878349559164 11.939639278557113 -1026.7361787853647 11.906593186372744 -1033.069204382051 11.873547094188375 -1039.3869060714103 11.840501002004007 -1045.6892781825009 11.807454909819638 -1051.976315047984 11.774408817635269 -1058.2480110041033 11.7413627254509 -1064.5043603906647 11.708316633266532 -1070.7453575510172 11.675270541082163 -1076.97099683203 11.642224448897794 -1083.1812725840748 11.609178356713425 -1089.3761791610038 11.576132264529058 -1095.5557109201313 11.543086172344688 -1101.7198622222118 11.51004008016032 -1107.8686274314211 11.47699398797595 -1114.0020009179275 11.443947895791583 -1120.1199770948845 11.410901803607214 -1126.2225503739382 11.377855711422844 -1132.3097151282118 11.344809619238475 -1138.3814657341777 11.311763527054108 -1144.43779657164 11.278717434869739 -1150.4787020237138 11.24567134268537 -1156.5041764768068 11.212625250501 -1162.5142143206008 11.179579158316631 -1168.5088099480326 11.146533066132264 -1174.4879577552774 11.113486973947895 -1180.451652141728 11.080440881763526 -1186.3998875099778 11.047394789579156 -1192.3326582658044 11.014348697394789 -1198.2499588181515 10.98130260521042 -1204.151783579111 10.94825651302605 -1210.0381269639067 10.915210420841682 -1215.9089833908783 10.882164328657314 -1221.7643472814655 10.849118236472945 -1227.6042130601904 10.816072144288576 -1233.4285751546442 10.783026052104207 -1239.2374279954706 10.74997995991984 -1245.0307660163533 10.71693386773547 -1250.808583653999 10.683887775551101 -1256.570875348126 10.650841683366732 -1262.3176355414496 10.617795591182364 -1268.048858679671 10.584749498997995 -1273.7645392133095 10.551703406813626 -1279.4646716091677 10.518657314629257 -1285.1492503327906 10.485611222444888 -1290.8182698369958 10.45256513026052 -1296.4717245775637 10.419519038076151 -1302.109609013228 10.386472945891782 -1307.731917605668 10.353426853707413 -1313.3386448195004 10.320380761523046 -1318.9297851222732 10.287334669338676 -1324.5053329844588 10.254288577154307 -1330.0652828794498 10.221242484969938 -1335.6096292835537 10.18819639278557 -1341.1383666759907 10.155150300601202 -1346.65148953889 10.122104208416832 -1352.1489923572894 10.089058116232463 -1357.6309361544413 10.056012024048096 -1363.0974158721303 10.022965931863727 -1368.5484604553164 9.989919839679358 -1373.9840654610061 9.956873747494988 -1379.4042264274146 9.923827655310621 -1384.8089388738244 9.890781563126252 -1390.1981983004378 9.857735470941883 -1395.5720001882319 9.824689378757514 -1400.9303399988105 9.791643286573144 -1406.2732131742534 9.758597194388777 -1411.6006151369666 9.725551102204408 -1416.9125412895264 9.692505010020039 -1422.2089870145269 9.65945891783567 -1427.4899476744213 9.626412825651302 -1432.7554186113655 9.593366733466933 -1438.0053951470554 9.560320641282564 -1443.2398725825662 9.527274549098195 -1448.4588461981875 9.494228456913827 -1453.6623112532582 9.461182364729458 -1458.8502629859972 9.428136272545089 -1464.022696613335 9.39509018036072 -1469.17960733074 9.362044088176352 -1474.3209903120467 9.328997995991983 -1479.4468407092775 9.295951903807614 -1484.5571536524658 9.262905811623245 -1489.6519242494755 9.229859719438878 -1494.7311475858196 9.196813627254508 -1499.7948187244733 9.16376753507014 -1504.8429327056892 9.13072144288577 -1509.8754845468072 9.097675350701401 -1514.892469242063 9.064629258517034 -1519.893881762396 9.031583166332664 -1524.8797170552505 8.998537074148295 -1529.8499700443792 8.965490981963926 -1534.8046356296413 8.932444889779559 -1539.7437086868 8.89939879759519 -1544.6671840673143 8.86635270541082 -1549.5750565981332 8.833306613226451 -1554.4673210814817 8.800260521042084 -1559.3439722946496 8.767214428857715 -1564.2050049897723 8.734168336673346 -1569.050413893613 8.701122244488976 -1573.8801937073401 8.668076152304609 -1578.6943391063026 8.63503006012024 -1583.4928447398006 8.60198396793587 -1588.2757052308561 8.568937875751502 -1593.0429151759784 8.535891783567134 -1597.7944691449272 8.502845691382765 -1602.5303616804727 8.469799599198396 -1607.2505872981512 8.436753507014027 -1611.9551404860206 8.403707414829658 -1616.644015704409 8.37066132264529 -1621.3172073856629 8.337615230460921 -1625.9747099338895 8.304569138276552 -1630.6165177246974 8.271523046092183 -1635.242625104933 8.238476953907815 -1639.8530263924147 8.205430861723446 -1644.4477158756597 8.172384769539077 -1649.026687813612 8.139338677354708 -1653.5899364353622 8.10629258517034 -1658.1374559398685 8.073246492985971 -1662.669240495667 8.040200400801602 -1667.1852842405845 8.007154308617235 -1671.6855812814438 7.974108216432865 -1676.1701256937645 7.9410621242484964 -1680.638911521463 7.908016032064127 -1685.091932776543 7.874969939879759 -1689.529183438787 7.84192384769539 -1693.9506574554389 7.808877755511022 -1698.3563487408856 7.775831663326652 -1702.7462511763304 7.742785571142283 -1707.1203586094648 7.709739478957915 -1711.4786648541349 7.676693386773546 -1715.8211636900003 7.643647294589178 -1720.1478488621915 7.610601202404808 -1724.4587140809601 7.57755511022044 -1728.7537530213242 7.544509018036071 -1733.0329593227082 7.511462925851703 -1737.2963265885785 7.478416833667334 -1741.5438483860719 7.445370741482965 -1745.7755182456197 7.412324649298596 -1749.9913296605646 7.379278557114228 -1754.191276086774 7.346232464929859 -1758.375350942244 7.31318637274549 -1762.5435476067014 7.280140280561121 -1766.6958594211953 7.247094188376753 -1770.8322796876855 7.214048096192384 -1774.952801668622 7.181002004008016 -1779.05741858652 7.147955911823646 -1783.1461236235266 7.114909819639278 -1787.2189099209806 7.081863727454909 -1791.2757705789666 7.04881763527054 -1795.3166986558592 7.015771543086172 -1799.3416871678637 6.982725450901802 -1803.3507290885448 6.949679358717434 -1807.3438173483507 6.916633266533065 -1811.3209448341286 6.883587174348697 -1815.2821043886313 6.8505410821643276 -1819.2272888100158 6.817494989979959 -1823.1564908513349 6.78444889779559 -1827.0697032200173 6.751402805611222 -1830.966918577342 6.718356713426853 -1834.8481295379013 6.685310621242484 -1838.7133286690555 6.652264529058115 -1842.5625084903781 6.619218436873747 -1846.3956614730916 6.586172344689378 -1850.2127800394926 6.55312625250501 -1854.0138565623672 6.52008016032064 -1857.7988833643958 6.487034068136272 -1861.5678527175482 6.453987975951903 -1865.3207568424662 6.420941883767535 -1869.0575879078363 6.387895791583166 -1872.77833802975 6.354849699398796 -1876.4829992710536 6.321803607214428 -1880.1715636406848 6.288757515030059 -1883.8440230929975 6.255711422845691 -1887.5003695270746 6.2226653306613215 -1891.1405947860258 6.189619238476953 -1894.7646906562752 6.156573146292584 -1898.3726488668326 6.123527054108216 -1901.9644610885523 6.090480961923847 -1905.5401189333772 6.057434869739478 -1909.0996139535682 6.024388777555109 -1912.642937640918 5.991342685370741 -1916.1700814259511 5.958296593186372 -1919.6810366771058 5.925250501002004 -1923.1757946999014 5.892204408817634 -1926.6543467360877 5.859158316633266 -1930.1166839627788 5.826112224448897 -1933.5627974915667 5.793066132264529 -1936.9926783676196 5.7600200400801596 -1940.4063175687593 5.726973947895791 -1943.8037060045206 5.693927855711422 -1947.1848345151907 5.660881763527054 -1950.5496938708284 5.627835671342685 -1953.898274770262 5.5947895791583155 -1957.230567840067 5.561743486973947 -1960.5465636335193 5.528697394789578 -1963.846252629529 5.49565130260521 -1967.1296252315478 5.462605210420841 -1970.3966717664537 5.429559118236472 -1973.647382483411 5.396513026052103 -1976.8817475527053 5.363466933867735 -1980.0997570645495 5.330420841683366 -1983.3014010278666 5.297374749498998 -1986.4866693690412 5.264328657314628 -1989.655551930644 5.23128256513026 -1992.808038470126 5.198236472945891 -1995.9441186584822 5.165190380761523 -1999.0637820788843 5.1321442885771535 -2002.1670182252783 5.099098196392785 -2005.2538165009519 5.066052104208416 -2008.3241662170644 5.033006012024048 -2011.3780565911425 4.999959919839679 -2014.4154767455366 4.96691382765531 -2017.4364157058421 4.933867735470941 -2020.440862399278 4.900821643286572 -2023.428805653026 4.867775551102204 -2026.4002341925273 4.834729458917835 -2029.3551366397337 4.801683366733466 -2032.2935015113171 4.768637274549097 -2035.2153172168275 4.735591182364729 -2038.1205720568057 4.70254509018036 -2041.0092542208438 4.6694989979959916 -2043.8813517855945 4.636452905811622 -2046.7368527127257 4.603406813627254 -2049.5757448468194 4.570360721442885 -2052.3980159132116 4.537314629258517 -2055.203653515774 4.5042685370741475 -2057.992645134629 4.471222444889779 -2060.7649781238056 4.43817635270541 -2063.520639708821 4.405130260521042 -2066.2596169841977 4.372084168336673 -2068.9818969109056 4.339038076152304 -2071.6874663137273 4.305991983967935 -2074.3763118785464 4.272945891783567 -2077.0484201495524 4.239899799599198 -2079.703777526361 4.206853707414829 -2082.3423702610476 4.17380761523046 -2084.964184455083 4.140761523046091 -2087.5692060561787 4.107715430861723 -2090.1574208550287 4.074669338677354 -2092.7288144819468 4.0416232464929855 -2095.2833724033944 4.008577154308617 -2097.8210799183958 3.975531062124248 -2100.341922154833 3.9424849699398794 -2102.845884065617 3.9094388777555107 -2105.33295042473 3.8763927855711415 -2107.803105823127 3.843346693386773 -2110.2563346645047 3.810300601202404 -2112.6926211609107 3.7772545090180354 -2115.111949328206 3.7442084168336667 -2117.5143029813594 3.711162324649298 -2119.8996657295747 3.6781162324649292 -2122.2680209712344 3.6450701402805605 -2124.6193518886594 3.612024048096192 -2126.9536414426702 3.578977955911823 -2129.2708723669416 3.5459318637274544 -2131.5710271621383 3.5128857715430857 -2133.854088089825 3.479839679358717 -2136.120037166133 3.4467935871743483 -2138.368856155175 3.4137474949899795 -2140.6005265621943 3.380701402805611 -2142.815029626433 3.347655310621242 -2145.012346313704 3.3146092184368734 -2147.19245730865 3.2815631262525047 -2149.355343006677 3.248517034068136 -2151.500983505536 3.2154709418837673 -2153.62935859654 3.182424849699398 -2155.7404477553896 3.1493787575150294 -2157.834230132592 3.1163326653306607 -2159.9106845434417 3.083286573146292 -2161.969789457541 3.0502404809619232 -2164.011522987832 3.0171943887775545 -2166.03586287911 2.984148296593186 -2168.042786495988 2.951102204408817 -2170.032270810274 2.9180561122244484 -2172.004292387731 2.8850100200400797 -2173.9588273741706 2.851963927855711 -2175.895851480848 2.8189178356713422 -2177.8153399691037 2.7858717434869735 -2179.7172676342084 2.752825651302605 -2181.6016087883518 2.719779559118236 -2183.4683372427235 2.6867334669338674 -2185.317426288618 2.6536873747494987 -2187.148848677502 2.62064128256513 -2188.9625765999613 2.5875951903807612 -2190.7585816634605 2.5545490981963925 -2192.536834868817 2.521503006012024 -2194.2973065853053 2.488456913827655 -2196.039966524281 2.455410821643286 -2197.764783711226 2.4223647294589172 -2199.471726456078 2.3893186372745485 -2201.160762321733 2.35627254509018 -2202.8318580905557 2.323226452905811 -2204.484979728765 2.2901803607214424 -2206.1200923485058 2.2571342685370737 -2207.737160167428 2.224088176352705 -2209.3361464655745 2.1910420841683362 -2210.9170135393438 2.1579959919839675 -2212.479722652293 2.124949899799599 -2214.024233982508 2.09190380761523 -2215.550506566251 2.0588577154308614 -2217.0584982375563 2.0258116232464927 -2218.548165563428 1.9927655310621242 -2220.019463774245 1.9597194388777555 -2221.4723466889436 1.9266733466933865 -2222.9067666345086 1.8936272545090178 -2224.322674359253 1.860581162324649 -2225.7200189393116 1.8275350701402804 -2227.0987476777163 1.7944889779559117 -2228.4588059953644 1.761442885771543 -2229.800137313103 1.7283967935871742 -2231.1226829240954 1.6953507014028055 -2232.4263818555323 1.6623046092184368 -2233.711170718666 1.629258517034068 -2234.9769835460424 1.5962124248496992 -2236.2237516146974 1.5631663326653304 -2237.4514032539546 1.5301202404809617 -2238.659863636356 1.497074148296593 -2239.8490545501013 1.4640280561122243 -2241.0188941512515 1.4309819639278556 -2242.169296693807 1.3979358717434869 -2243.3001722356366 1.3648897795591182 -2244.411426318109 1.3318436873747495 -2245.502959617184 1.2987975951903807 -2246.5746675636447 1.265751503006012 -2247.6264399301567 1.232705410821643 -2248.658160382909 1.1996593186372744 -2249.6697059958124 1.1666132264529057 -2250.6609467256144 1.133567134268537 -2251.6317448469495 1.1005210420841682 -2252.581954347367 1.0674749498997995 -2253.5114202839122 1.0344288577154308 -2254.419978105077 1.001382765531062 -2255.307452945161 0.9683366733466933 -2256.1736589026286 0.9352905811623246 -2257.0183983204206 0.9022444889779558 -2257.84146109503 0.8691983967935871 -2258.6426240533374 0.8361523046092184 -2259.421650452904 0.8031062124248496 -2260.17828968417 0.7700601202404809 -2260.912277283807 0.7370140280561122 -2261.6233354099572 0.7039679358717434 -2262.311173985597 0.6709218436873747 -2262.9754927899785 0.637875751503006 -2263.615984875288 0.6048296593186372 -2264.232341812548 0.5717835671342685 -2264.824261434824 0.5387374749498998 -2265.3914589556684 0.505691382765531 -2265.9336826066656 0.4726452905811623 -2266.450735273437 0.43959919839679357 -2266.9425040355172 0.4065531062124248 -2267.4090000731435 0.3735070140280561 -2267.850412184595 0.34046092184368737 -2268.267178376049 0.3074148296593186 -2268.660082153422 0.2743687374749499 -2269.0303844912305 0.24132264529058114 -2269.3800120148485 0.2082765531062124 -2269.7118446943537 0.1752304609218437 -2270.030205654177 0.14218436873747495 -2270.3418319261446 0.10913827655310619 -2270.6582520480956 0.07609218436873746 -2271.0039656545337 0.04304609218436874 -2271.462498726208 0.01 +-0.0 80.0 +154.18305973173608 79.83969939879759 +308.12735926191283 79.6793987975952 +461.8327373234883 79.51909819639279 +615.2990316169373 79.35879759519038 +768.5260796125001 79.19849699398797 +921.5137185493666 79.03819639278558 +1074.2617854349387 78.87789579158317 +1226.770117043952 78.71759519038076 +1379.0385499177276 78.55729458917835 +1531.066920363351 78.39699398797595 +1682.855064452848 78.23669338677355 +1834.4028180224204 78.07639278557114 +1985.7100166715445 77.91609218436874 +2136.776495762197 77.75579158316633 +2287.6020904179973 77.59549098196393 +2438.1866355234224 77.43519038076153 +2588.5299657228848 77.27488977955912 +2738.631915419939 77.11458917835671 +2888.4923187764016 76.95428857715432 +3038.111009711546 76.79398797595191 +3187.4878219011594 76.6336873747495 +3336.622588776726 76.47338677354709 +3485.5151435245275 76.3130861723447 +3634.165319084815 76.15278557114229 +3782.5729481508456 75.99248496993988 +3930.737863168042 75.83218436873747 +4078.659896333071 75.67188376753508 +4226.338879592993 75.51158316633267 +4373.774644644274 75.35128256513026 +4520.967022931927 75.19098196392785 +4667.915845648577 75.03068136272545 +4814.620943733519 74.87038076152305 +4961.082147871845 74.71008016032064 +5107.299288493425 74.54977955911824 +5253.272195772012 74.38947895791583 +5399.000699624273 74.22917835671343 +5544.484629708876 74.06887775551102 +5689.72381542546 73.90857715430862 +5834.718085913709 73.74827655310621 +5979.467273154079 73.58797595190381 +6123.9712090511475 73.4276753507014 +6268.229722135261 73.267374749499 +6412.242640477294 73.10707414829659 +6556.00979188318 72.9467735470942 +6699.531003892951 72.78647294589179 +6842.806103779651 72.62617234468938 +6985.834918548363 72.46587174348697 +7128.617274935159 72.30557114228458 +7271.152999406111 72.14527054108217 +7413.441918156177 71.98496993987976 +7555.483857108203 71.82466933867735 +7697.278641911851 71.66436873747494 +7838.8260979425295 71.50406813627255 +7980.126050300369 71.34376753507014 +8121.178323809068 71.18346693386773 +8261.982743014867 71.02316633266533 +8402.539132185419 70.86286573146293 +8542.84731530875 70.70256513026052 +8682.907116092074 70.54226452905812 +8822.718357960724 70.38196392785571 +8962.280864057015 70.22166332665331 +9101.594457239149 70.0613627254509 +9240.658960080014 69.9010621242485 +9379.474194866087 69.74076152304609 +9518.039983596247 69.5804609218437 +9656.35614798067 69.42016032064129 +9794.422509439586 69.25985971943888 +9932.238889102144 69.09955911823647 +10069.805107805203 68.93925851703408 +10207.120986092188 68.77895791583167 +10344.186344211814 68.61865731462926 +10481.001002116918 68.45835671342685 +10617.56477946323 68.29805611222444 +10753.877495608132 68.13775551102205 +10889.938970125799 67.97745490981964 +11025.749025173951 67.81715430861723 +11161.307482081154 67.65685370741483 +11296.614158977436 67.49655310621243 +11431.66887368821 67.33625250501002 +11566.471443732924 67.17595190380761 +11701.021686323815 67.0156513026052 +11835.319418364583 66.85535070140281 +11969.364456449124 66.6950501002004 +12103.156616860146 66.534749498998 +12236.695715567874 66.37444889779559 +12369.981568228695 66.2141482965932 +12503.013990183852 66.05384769539079 +12635.792796458009 65.89354709418838 +12768.317801757934 65.73324649298597 +12900.588820471092 65.57294589178358 +13032.605666664298 65.41264529058117 +13164.36815408225 65.25234468937876 +13295.876096146163 65.09204408817635 +13427.129305952329 64.93174348697394 +13558.127596270679 64.77144288577155 +13688.870779543386 64.61114228456914 +13819.358667883333 64.45084168336673 +13949.591073072701 64.29054108216432 +14079.567806561465 64.13024048096193 +14209.288679465957 63.969939879759515 +14338.753502567275 63.80963927855711 +14467.962086309848 63.649338677354706 +14596.914240799899 63.4890380761523 +14725.609775803876 63.3287374749499 +14854.048500746952 63.16843687374749 +14982.230224711424 63.00813627254509 +15110.154756435179 62.84783567134268 +15237.821904828645 62.68753507014028 +15365.231481461387 62.52723446893787 +15492.383297022994 62.36693386773547 +15619.27715883334 62.20663326653306 +15745.912873859768 62.04633266533066 +15872.290248715499 61.88603206412825 +15998.409089657915 61.72573146292585 +16124.269202586942 61.56543086172344 +16249.870393043318 61.405130260521034 +16375.212466206913 61.24482965931863 +16500.295226895032 61.084529058116225 +16625.118479560646 60.92422845691382 +16749.682028290703 60.763927855711415 +16873.985676804325 60.603627254509014 +16998.029228451072 60.443326653306606 +17121.81248620914 60.283026052104205 +17245.335252683577 60.1227254509018 +17368.597330104447 59.962424849699396 +17491.598520325042 59.80212424849699 +17614.338624819986 59.641823647294586 +17736.81744468343 59.48152304609218 +17859.034780627135 59.32122244488978 +17980.990432978622 59.16092184368737 +18102.68420167922 59.00062124248497 +18224.11588628221 58.84032064128256 +18345.285285950806 58.68002004008016 +18466.192199456287 58.51971943887775 +18586.836425175956 58.35941883767535 +18707.21776109121 58.19911823647294 +18827.336004785502 58.03881763527053 +18947.190953442318 57.87851703406813 +19066.782406452363 57.718216432865724 +19186.110163856676 57.55791583166332 +19305.174022667197 57.397615230460914 +19423.97377841939 57.23731462925851 +19542.50922623837 57.077014028056105 +19660.78016083675 56.916713426853704 +19778.78637651254 56.756412825651296 +19896.527667146915 56.596112224448895 +20014.0038262021 56.43581162324649 +20131.214646719127 56.275511022044086 +20248.159921315633 56.11521042084168 +20364.839442183602 55.954909819639276 +20481.25300108713 55.79460921843687 +20597.40038936011 55.63430861723447 +20713.281397903975 55.47400801603206 +20828.895817185323 55.31370741482966 +20944.243437233632 55.15340681362725 +21059.324047638842 54.99310621242485 +21174.137437549023 54.83280561122244 +21288.683395667922 54.67250501002003 +21402.961710252555 54.51220440881763 +21516.972169110777 54.35190380761522 +21630.714559598768 54.19160320641282 +21744.18866861858 54.03130260521041 +21857.394282615576 53.87100200400801 +21970.331187575943 53.710701402805604 +22082.999169024064 53.5504008016032 +22195.39801201999 53.390100200400795 +22307.52750361747 53.229799599198394 +22419.3874317656 53.069498997995986 +22530.977581464762 52.909198396793585 +22642.29773585974 52.748897795591176 +22753.34767761759 52.588597194388775 +22864.12718892491 52.42829659318637 +22974.63605148507 52.267995991983966 +23084.874046515448 52.10769539078156 +23194.840954744563 51.94739478957916 +23304.53655640929 51.78709418837675 +23413.96063125194 51.62679358717435 +23523.112958517395 51.46649298597194 +23631.993316950164 51.30619238476954 +23740.60148479144 51.14589178356713 +23848.937239776104 50.98559118236472 +23957.00035912971 50.82529058116232 +24064.790619565472 50.66498997995991 +24172.30779728115 50.50468937875751 +24279.551667955988 50.3443887775551 +24386.522006747535 50.1840881763527 +24493.21858828853 50.023787575150294 +24599.641186683668 49.86348697394789 +24705.789575506402 49.703186372745485 +24811.66352779565 49.542885771543084 +24917.26281605253 49.382585170340676 +25022.587212237024 49.222284569138274 +25127.63649080643 49.061983967935866 +25232.410426636045 48.901683366733465 +25336.90879099011 48.74138276553106 +25441.13135359422 48.581082164328656 +25545.077883615642 48.42078156312625 +25648.748149659736 48.26048096192385 +25752.14191976643 48.10018036072144 +25855.258961406573 47.93987975951904 +25958.099041478315 47.77957915831663 +26060.661926303397 47.61927855711422 +26162.947381623442 47.45897795591182 +26264.955172596197 47.29867735470941 +26366.68506379171 47.13837675350701 +26468.136819188505 46.9780761523046 +26569.31020216968 46.8177755511022 +26670.20497551901 46.65747494989979 +26770.820901416937 46.49717434869739 +26871.157741436604 46.336873747494984 +26971.21525653976 46.17657314629258 +27070.993207072694 46.016272545090175 +27170.491352762063 45.85597194388777 +27269.70945271074 45.695671342685365 +27368.64726539354 45.535370741482964 +27467.30454865298 45.375070140280556 +27565.681063441407 45.214769539078155 +27663.776570538143 45.05446893787575 +27761.590826310017 44.894168336673346 +27859.123585985417 44.73386773547094 +27956.374604138182 44.57356713426854 +28053.343634683075 44.41326653306613 +28150.03043087112 44.25296593186372 +28246.434745284932 44.09266533066132 +28342.556329834002 43.93236472945891 +28438.39493574988 43.77206412825651 +28533.950313581357 43.6117635270541 +28629.22221318955 43.4514629258517 +28724.210383742968 43.29116232464929 +28818.91457371246 43.13086172344689 +28913.334530866203 42.97056112224448 +29007.47000226451 42.81026052104208 +29101.32073425469 42.649959919839674 +29194.886472465754 42.48965931863727 +29288.16696180315 42.329358717434864 +29381.161946443324 42.16905811623246 +29473.871169828344 42.008757515030055 +29566.294374660334 41.848456913827654 +29658.43130679031 41.688156312625246 +29750.28171187288 41.527855711422845 +29841.84533088793 41.36755511022044 +29933.12190347435 41.207254509018036 +30024.11116850251 41.04695390781563 +30114.812864068266 40.88665330661322 +30205.226727487014 40.72635270541082 +30295.352495287614 40.56605210420841 +30385.18990320623 40.40575150300601 +30474.73868618012 40.2454509018036 +30563.998578341336 40.0851503006012 +30652.969313010344 39.92484969939879 +30741.65062268956 39.76454909819639 +30830.042239056827 39.60424849699398 +30918.143892958757 39.44394789579158 +31005.95531440406 39.28364729458917 +31093.476232556706 39.12334669338677 +31180.706375729093 38.96304609218436 +31267.64547137501 38.80274549098196 +31354.293246082634 38.642444889779554 +31440.649428275603 38.48214428857715 +31526.713747520695 38.321843687374745 +31612.485929761842 38.161543086172344 +31697.96569799777 38.001242484969936 +31783.152774324688 37.840941883767535 +31868.046879928675 37.68064128256513 +31952.647735077913 37.52034068136272 +32036.95505911483 37.36004008016032 +32120.96857044816 37.19973947895791 +32204.687986544846 37.03943887775551 +32288.113023921895 36.8791382765531 +32371.24339813804 36.7188376753507 +32454.07882378539 36.55853707414829 +32536.61901448084 36.39823647294589 +32618.86368285749 36.23793587174348 +32700.81254055583 36.07763527054108 +32782.4652982149 35.91733466933867 +32863.82166546321 35.75703406813627 +32944.88135090969 35.59673346693386 +33025.64406641419 35.43643286573146 +33106.10952353573 35.27613226452905 +33186.27742847418 35.11583166332665 +33266.1474856031 34.955531062124244 +33345.7193982285 34.79523046092184 +33424.99286857898 34.634929859719435 +33503.9675977955 34.474629258517034 +33582.643285921215 34.314328657314626 +33661.019631891 34.15402805611222 +33739.09633352089 33.993727454909816 +33816.873087497355 33.83342685370741 +33894.34958936639 33.67312625250501 +33971.52553352248 33.5128256513026 +34048.40061319734 33.3525250501002 +34124.97452044852 33.19222444889779 +34201.246946147854 33.03192384769539 +34277.21757996967 32.87162324649298 +34352.88611352987 32.71132264529058 +34428.252239447604 32.55102204408817 +34503.31564592074 32.39072144288577 +34578.07601763951 32.23042084168336 +34652.53303803486 32.07012024048096 +34726.68638926552 31.909819639278556 +34800.53575220487 31.74951903807615 +34874.080806427635 31.589218436873747 +34947.32123019628 31.428917835671342 +35020.25670044725 31.268617234468937 +35092.88689277696 31.108316633266533 +35165.21148142754 30.948016032064128 +35237.23013927238 30.787715430861724 +35308.94253780138 30.62741482965932 +35380.34834710598 30.467114228456914 +35451.44723586397 30.30681362725451 +35522.245081721965 30.146513026052105 +35592.75268234916 29.9862124248497 +35662.974306361895 29.825911823647296 +35732.908983693706 29.66561122244489 +35802.55574571999 29.505310621242486 +35871.913625289664 29.345010020040082 +35940.981656756754 29.184709418837677 +36009.75887601199 29.02440881763527 +36078.24432051431 28.864108216432864 +36146.43702932237 28.70380761523046 +36214.33604312593 28.543507014028055 +36281.94040427729 28.38320641282565 +36349.24915682255 28.222905811623246 +36416.26134653289 28.06260521042084 +36482.97602093569 27.902304609218437 +36549.39223472073 27.742004008016032 +36615.509044379934 27.581703406813627 +36681.325502741965 27.421402805611223 +36746.840664394906 27.261102204408818 +36812.05358583889 27.100801603206413 +36876.9633255167 26.94050100200401 +36941.56894384422 26.780200400801604 +37005.86950324079 26.6198997995992 +37069.864068159404 26.459599198396795 +37133.551705116726 26.29929859719439 +37196.93148272303 26.138997995991986 +37260.002471711865 25.97869739478958 +37322.76374496967 25.818396793587176 +37385.21437986621 25.65809619238477 +37447.353458834885 25.497795591182363 +37509.18006409903 25.33749498997996 +37570.69327715053 25.177194388777554 +37631.89218178084 25.01689378757515 +37692.77586410935 24.856593186372745 +37753.34341261153 24.69629258517034 +37813.59391814675 24.535991983967936 +37873.52647398588 24.37569138276553 +37933.14017583867 24.215390781563126 +37992.43412188072 24.05509018036072 +38051.40741278028 23.894789579158317 +38110.05915172469 23.734488977955913 +38168.388449131824 23.574188376753508 +38226.39441807867 23.413887775551103 +38284.076169392734 23.2535871743487 +38341.43281627571 23.093286573146294 +38398.463474578086 22.93298597194389 +38455.16726282355 22.772685370741485 +38511.543302232916 22.61238476953908 +38567.59071674772 22.452084168336675 +38623.30863305334 22.29178356713427 +38678.69618060169 22.131482965931863 +38733.75249163343 21.971182364729458 +38788.4767032246 21.810881763527053 +38842.867957793664 21.65058116232465 +38896.9253984136 21.490280561122244 +38950.648168348 21.32997995991984 +39004.03541375611 21.169679358717435 +39057.08628371198 21.00937875751503 +39109.79993022304 20.849078156312626 +39162.17550824803 20.68877755511022 +39214.212175714354 20.528476953907816 +39265.909093534734 20.36817635270541 +39317.265425623234 20.207875751503007 +39368.280342404236 20.047575150300602 +39418.953017958884 19.887274549098198 +39469.282625735475 19.726973947895793 +39519.26834143274 19.56667334669339 +39568.909343824096 19.406372745490984 +39618.20481476919 19.24607214428858 +39667.15393922459 19.085771543086175 +39715.75590525363 18.92547094188377 +39764.00990403534 18.76517034068136 +39811.91512995735 18.604869739478957 +39859.47078435391 18.444569138276552 +39906.67607141289 18.284268537074148 +39953.530194452454 18.123967935871743 +40000.03235994172 17.96366733466934 +40046.18177750371 17.803366733466934 +40091.97765991722 17.64306613226453 +40137.41922311743 17.482765531062125 +40182.50568619546 17.32246492985972 +40227.236271396585 17.162164328657315 +40271.61020759174 17.00186372745491 +40315.62672662542 16.841563126252506 +40359.285059801325 16.6812625250501 +40402.58444152622 16.520961923847697 +40445.52410933893 16.360661322645292 +40488.10330390047 16.200360721442888 +40530.32126898268 16.040060120240483 +40572.177251455294 15.879759519038075 +40613.67050333278 15.71945891783567 +40654.80028056926 15.559158316633265 +40695.56583988652 15.39885771543086 +40735.96643994735 15.238557114228456 +40776.001342430856 15.078256513026052 +40815.669812009495 14.917955911823647 +40854.97111632417 14.757655310621242 +40893.90452595757 14.597354709418838 +40932.4693169715 14.437054108216431 +40970.66476813032 14.276753507014027 +41008.49015825209 14.116452905811622 +41045.94476892484 13.956152304609217 +41083.027884525174 13.795851703406813 +41119.73879217987 13.635551102204408 +41156.076781725504 13.475250501002003 +41192.04114756819 13.314949899799599 +41227.63118672216 13.154649298597194 +41262.84619663469 12.99434869739479 +41297.68547706074 12.834048096192385 +41332.14833024534 12.673747494989978 +41366.23406087208 12.513446893787574 +41399.941976606264 12.35314629258517 +41433.27138830625 12.192845691382765 +41466.22160825391 12.03254509018036 +41498.79194983882 11.872244488977955 +41530.981728625164 11.71194388777555 +41562.790262294024 11.551643286573146 +41594.216870937526 11.391342685370741 +41625.26087714943 11.231042084168337 +41655.92160485331 11.07074148296593 +41686.19837910929 10.910440881763526 +41716.09052683559 10.750140280561121 +41745.59737676145 10.589839679358716 +41774.718259621564 10.429539078156312 +41803.45250770008 10.269238476953907 +41831.79945448089 10.108937875751502 +41859.76045661166 9.948637274549098 +41887.337411157394 9.788336673346693 +41914.530250393924 9.628036072144289 +41941.33841792872 9.467735470941884 +41967.761345412386 9.307434869739478 +41993.79845206821 9.147134268537073 +42019.44914419414 8.986833667334668 +42044.71281463526 8.826533066132264 +42069.58884222413 8.666232464929859 +42094.07659118656 8.505931863727454 +42118.17541050981 8.34563126252505 +42141.88463327002 8.185330661322645 +42165.20357591537 8.02503006012024 +42188.13153750087 7.864729458917835 +42210.66779887049 7.70442885771543 +42232.81162178161 7.544128256513026 +42254.56224796631 7.383827655310621 +42275.91889812315 7.2235270541082155 +42296.880770832635 7.063226452905811 +42317.4470413882 6.902925851703406 +42337.61686053389 6.742625250501002 +42357.38935309841 6.582324649298597 +42376.76361651393 6.422024048096192 +42395.73871920611 6.261723446893787 +42414.31369884013 6.101422845691382 +42432.4875604047 5.9411222444889775 +42450.25927411347 5.780821643286573 +42467.62777309974 5.620521042084168 +42484.59195087608 5.460220440881763 +42501.15065852565 5.299919839679358 +42517.30270158569 5.1396192384769535 +42533.04683657607 4.979318637274549 +42548.38176711646 4.819018036072144 +42563.30613956394 4.658717434869739 +42577.81853808798 4.498416833667334 +42591.91747908106 4.338116232464929 +42605.60140477906 4.177815631262525 +42618.868675934544 4.01751503006012 +42631.71756334523 3.857214428857715 +42644.146237986286 3.6969138276553104 +42656.15275942319 3.5366132264529053 +42667.73506208455 3.3763126252505007 +42678.89093884083 3.216012024048096 +42689.6180211494 3.055711422845691 +42699.91375476411 2.8954108216432863 +42709.77536963178 2.7351102204408813 +42719.19984204987 2.5748096192384766 +42728.18384634584 2.414509018036072 +42736.723692108964 2.254208416833667 +42744.81524111143 2.0939078156312623 +42752.45379508432 1.9336072144288576 +42759.63394077371 1.7733066132264528 +42766.34933103669 1.6130060120240481 +42772.592368292775 1.4527054108216433 +42778.35373678859 1.2924048096192384 +42783.62170074818 1.1321042084168336 +42788.38105321053 0.9718036072144288 +42792.61162030574 0.8115030060120241 +42796.2865773871 0.6512024048096192 +42799.37275039113 0.4909018036072144 +42801.8433788484 0.3306012024048096 +42803.74630774719 0.1703006012024048 +42805.973184099894 0.01 diff --git a/ELoss/PCEnergyAnalysis.py b/ELoss/PCEnergyAnalysis.py index ce708ff..e0d8c9d 100644 --- a/ELoss/PCEnergyAnalysis.py +++ b/ELoss/PCEnergyAnalysis.py @@ -20,7 +20,6 @@ import os import periodictable as pt import re from matplotlib.colors import LinearSegmentedColormap -import matplotlib as mpl # ROOT-like styling plt.rcParams.update({ @@ -79,9 +78,8 @@ interp_cache = {} def make_E_vs_x(z, mass_u, emax_mev, label, npoints, P_TORR, TEMP_K): # GAS SETUP R = 8.3144 - P_TORR = 400 - TEMP_K = 293.15 - R = 8.3144 + P_TORR = P_TORR + TEMP_K = TEMP_K # Gas density p_pa = P_TORR * 133.322 @@ -221,8 +219,6 @@ def energy_loss(particle, Ei, dl): Ef = E_of_x(xf) return np.maximum(Ef, 0.0) - -#def energy_loss_angular(particle, Ei, theta): def energy_reconstruction(particle, Ef, dl): @@ -242,11 +238,10 @@ def energy_distance(particle, Ei, Ef): xi = x_of_E(Ei) - xf = x_of_E(Ef) + xf = x_of_E(Ef) * .686 return np.abs(xf - xi) - def resolve_particle(name): name = name.lower().strip().rstrip("s") if name in particles: @@ -273,8 +268,7 @@ def resolve_particle(name): return elem.number, elem.mass, 30.0, name except Exception: raise ValueError(f"Unknown particle/isotope: {name}") - - + def update_plot_data(name, values): for i, (existing_name, _) in enumerate(plot_data): if existing_name == name: @@ -296,19 +290,17 @@ def calculate_distance_tree2(vz, theta, z_max=34.86): cos_theta = np.where(np.abs(cos_theta) < 1e-10, 1e-10, cos_theta) return np.abs((z_max - vz) / cos_theta) - def load_tree_arrays(tree, treename, max_events=None): branches = ["Tb", "thetab", "sx3Z", "vX", "vY", "vZ", "sx3ID", "aID", "cID", "aDist", "cDist"] if treename == 'tree1': branches.extend(["sx3X", "sx3Y"]) return tree.arrays(branches, library="np", entry_stop=max_events) - def prepare_tree_data(tree, treename, particle, max_events=None, z_max=34.86): data = load_tree_arrays(tree, treename, max_events) - mask = data["thetab"] >= 0 - data = {key: value[mask] for key, value in data.items()} + #mask = data["thetab"] >= 0 + #data = {key: value[mask] for key, value in data.items()} Ei = data["Tb"] theta = np.radians(data["thetab"]) vX = data["vX"] @@ -340,7 +332,7 @@ def prepare_tree_data(tree, treename, particle, max_events=None, z_max=34.86): sin_theta = np.sin(theta) sin_theta = np.where(sin_theta != 0, sin_theta, 1e-10) - radii = np.array([3.7, 4.2]) + radii = np.array([3.7, 4.3]) dA = (radii[0] - np.sqrt((vX/10)**2 + (vY/10)**2))/ sin_theta dC = (radii[1] - np.sqrt((vX/10)**2 + (vY/10)**2))/ sin_theta @@ -392,10 +384,9 @@ def prepare_tree_data(tree, treename, particle, max_events=None, z_max=34.86): "Eprop": Eprop, "dA": dA, "dC": dC, - "thetab": data["thetab"] + "thetab": np.degrees(theta) } - def process_file(filename, treename, particle=None, max_events=None): tree = uproot.open(filename)[f"{treename}"] @@ -410,7 +401,6 @@ def process_file(filename, treename, particle=None, max_events=None): print(f"File {filename} particle {particle}, tree {treename}") return prepare_tree_data(tree, treename, particle, max_events=max_events) - def power_fit_and_plot(x, y, label, color=None): x = np.array(x) @@ -487,12 +477,10 @@ class MyInteractiveApp(cmd.Cmd): except: self.file = None print("\nATTENTION: Root file not found, continue without uproot functions or uproot file manually") - #intro = "Interactive Shell Started. Type 'help' to see commands." prompt = ">> " - def default(self, line): # Check if the command starts with our multi-word phrase if line.startswith("make table "): @@ -761,7 +749,6 @@ class MyInteractiveApp(cmd.Cmd): import readline import atexit import os - import awkward as ak histfile = os.path.expanduser("~/.uproot_shell_history") @@ -781,7 +768,6 @@ class MyInteractiveApp(cmd.Cmd): "file": self.file, "uproot": uproot, "np": np, - #"ak": ak, "plt": plt, "self": self } @@ -820,7 +806,7 @@ class MyInteractiveApp(cmd.Cmd): look up 'uproot' for more details""" self.run_command_line() - def do_energy_analysis(self, arg): + def do_energy_analysis(self, arg): #geometry needs update args = shlex.split(arg) @@ -840,7 +826,6 @@ class MyInteractiveApp(cmd.Cmd): print(f"Using TTree: {self.tree}") treename = self.tree.name if hasattr(self.tree, 'name') else "tree1" - branches = ["Tb", "thetab", "sx3Z", "vX", "vY", "vZ", "sx3X", "sx3Y"] Tb_key = "Tb" @@ -901,7 +886,7 @@ class MyInteractiveApp(cmd.Cmd): global EA EA = energy_loss(particle, Ei, dA) global EC - EC = energy_loss(particle, Ei, ClassSX3dC) + EC = energy_loss(particle, Ei, dC) global Esx3 Esx3 = energy_loss(particle, Ei, dsx3) global Eprop @@ -910,25 +895,15 @@ class MyInteractiveApp(cmd.Cmd): Elost = Ei - Esx3 print("Analysis complete") - print(f"Processed events: {len(Ei)}") - print(f"Anode average energy: {np.mean(EA):.3f} MeV") - print(f"Cathode average energy: {np.mean(EC):.3f} MeV") - print(f"sx3 average energy: {np.mean(Esx3):.3f} MeV") - print(f"Average total energy loClassSX3ss to sx3: {np.mean(Elost):.3f} MeV") - print(f"Maximum total energy loss to sx3: {np.max(Elost):.3f} MeV") - print(f"Minimum total energy loss to sx3: {np.min(Elost):.3f} MeV") - print(f"Proportion counter average energy difference: {np.mean(Eprop):.3f} MeV") - print(f"Maximum proportion counter energy difference: {np.max(Eprop):.3f} MeV") - print(f"Minimum proportion counter energy difference: {np.min(Eprop):.3f} MeV") @@ -937,46 +912,30 @@ class MyInteractiveApp(cmd.Cmd): print(f"Writing new tree to {output_filename}") # Load ALL original branches - all_data = self.tree.arrays( - library="np", - entry_stop=max_events - ) + all_data = self.tree.arrays(library="np",entry_stop=max_events) # Create full-length arrays initialized to NaN n_total = len(data["Tb"]) EA_full = np.full(n_total, np.nan) - EC_full = np.full(n_total, np.nan) - Esx3_full = np.full(n_total, np.nan) - Eprop_full = np.full(n_total, np.nan) - Elost_full = np.full(n_total, np.nan) # Put values back into valid entries EA_full[mask] = EA - EC_full[mask] = EC - Esx3_full[mask] = Esx3 - Eprop_full[mask] = Eprop - Elost_full[mask] = Elost # Add new branches all_data["EA"] = EA_full - all_data["EC"] = EC_full - all_data["Esx3"] = Esx3_full - all_data["Eprop"] = Eprop_full - all_data["Elost"] = Elost_full - # Write new ROOT file as a classic TTree with uproot.recreate(output_filename) as fout: @@ -1010,17 +969,15 @@ class MyInteractiveApp(cmd.Cmd): data = prepare_tree_data(self.tree, treename, particle, max_events=max_events) - mask = data["thetab"] >= 0 - - Ei = data["Ei"][mask] - sx3Z = data["sx3Z"][mask] - EA = data["EA"][mask] - EC = data["EC"][mask] - Esx3 = data["Esx3"][mask] - Elost = data["Elost"][mask] - Eprop = data["Eprop"][mask] - dA = data["dA"][mask] - dC = data["dC"][mask] + Ei = data["Ei"] + sx3Z = data["sx3Z"] + EA = data["EA"] + EC = data["EC"] + Esx3 = data["Esx3"] + Elost = data["Elost"] + Eprop = data["Eprop"] + dA = data["dA"] + dC = data["dC"] update_plot_data(f"{particle}_{treename}_Ei", Ei) update_plot_data(f"{particle}_{treename}_sx3Z", sx3Z) @@ -1170,7 +1127,7 @@ class MyInteractiveApp(cmd.Cmd): # Keep only events with thetab >= 45 #thetab_mask = all_branches["thetab"] >= 0 - thetab_mask = mask + for branch in branch_names: values = all_branches[branch] @@ -1186,7 +1143,7 @@ class MyInteractiveApp(cmd.Cmd): continue # Apply the thetab cut - values = values[thetab_mask] + values = values values = values[~np.isnan(values)] @@ -1212,10 +1169,6 @@ class MyInteractiveApp(cmd.Cmd): def do_hist_comp(self, arg): particle = arg - self.do_set_tree("tree1") - data1 = prepare_tree_data(self.tree, "tree1", particle, max_events=None) - self.do_set_tree("tree2") - data2 = prepare_tree_data(self.tree, "tree2", particle, max_events=None) base = f"{particle}_comp_plots" os.makedirs(base, exist_ok=True) @@ -1287,7 +1240,20 @@ class MyInteractiveApp(cmd.Cmd): #tree2_name = args[3] if len(args) > 3 else 'tree1' outdir = "dual_plots" - data1 = process_file(os.path.join("..", "Armory", file1), "tree1") + # load both trees for file1 and combine their arrays into a single dataset + data1_tree1 = process_file(os.path.join("..", "Armory", file1), "tree1") + data1_tree2 = process_file(os.path.join("..", "Armory", file1), "tree2") + # concatenate matching array entries + data1 = {"particle": f"{data1_tree1['particle']}_combined"} + for key in data1_tree1: + if key == "particle": + continue + try: + data1[key] = np.concatenate([data1_tree1[key], data1_tree2[key]]) + except Exception: + # fallback: prefer tree1 value if concat fails + data1[key] = data1_tree1[key] + data2 = process_file(os.path.join("..", "Armory", file2), "tree1") #print(f"File one {file1} ({tree1_name}) length {len(data1['Ei'])} \nFile two {file2} ({tree2_name}) length {len(data2['Ei'])}") @@ -1295,27 +1261,10 @@ class MyInteractiveApp(cmd.Cmd): print(f"Saving plots to: {outdir}") - #Overlay histogram: Elost plt.figure(figsize=(8,6)) - - plt.hist( - data1["Elost"], - bins=200, - histtype='step', - linewidth=2, - density=True, - label=data1["particle"] - ) - - plt.hist( - data2["Elost"], - bins=200, - histtype='step', - linewidth=2, - density=True, - label=data2["particle"] - ) + plt.hist(data1["Elost"],bins=200,histtype='step',linewidth=2,density=True,label=data1["particle"]) + plt.hist(data2["Elost"],bins=200,histtype='step',linewidth=2,density=True,label=data2["particle"]) plt.xlabel("Energy Loss (MeV)") plt.ylabel("Normalized Counts") @@ -1329,106 +1278,81 @@ class MyInteractiveApp(cmd.Cmd): #Overlay histogram: sx3Z plt.figure(figsize=(8,6)) - - plt.hist( - data1["sx3Z"], - bins=150, - histtype='step', - linewidth=2, - density=True, - label=data1["particle"] - ) - - plt.hist( - data2["sx3Z"], - bins=150, - histtype='step', - linewidth=2, - density=True, - label=data2["particle"] - ) - + plt.hist(data1["sx3Z"],bins=150,histtype='step',linewidth=2,density=True,label=data1["particle"]) + plt.hist(data2["sx3Z"],bins=150,histtype='step',linewidth=2,density=True,label=data2["particle"]) plt.xlabel("SX3 Z") plt.ylabel("Normalized Counts") plt.title("SX3 Position Comparison") plt.legend() plt.grid(True) plt.tight_layout() - plt.savefig(f"{outdir}/sx3Z_overlay.png", dpi=300) plt.show() - - #Side-by-side 2D plots - fig, axes = plt.subplots(1, 2, figsize=(14,6)) - - h1 = axes[0].hist2d( - data1["sx3Z"], - data1["Elost"], - bins=200 - ) - - axes[0].set_title(f'{data1["particle"]} Elost vs SX3') - axes[0].set_xlabel("SX3 Z") - axes[0].set_ylabel("Energy Loss (MeV)") - - h2 = axes[1].hist2d( - data2["sx3Z"], - data2["Elost"], - bins=200 - ) - - axes[1].set_title(f'{data2["particle"]} Elost vs SX3') - axes[1].set_xlabel("SX3 Z") - axes[1].set_ylabel("Energy Loss (MeV)") - - fig.colorbar(h1[3], ax=axes[0], label="Counts") - fig.colorbar(h2[3], ax=axes[1], label="Counts") - - plt.tight_layout() - - plt.savefig(f"{outdir}/Elost_vs_sx3_comparison.png", dpi=300) - plt.show() + try: + #Side-by-side 2D plots + fig, axes = plt.subplots(1, 2, figsize=(14,6)) + + h1 = axes[0].hist2d(data1["sx3Z"],data1["Elost"],bins=200) + + axes[0].set_title(f'{data1["particle"]} Elost vs SX3') + axes[0].set_xlabel("SX3 Z") + axes[0].set_ylabel("Energy Loss (MeV)") + + h2 = axes[1].hist2d(data2["sx3Z"],data2["Elost"],bins=200) + + axes[1].set_title(f'{data2["particle"]} Elost vs SX3') + axes[1].set_xlabel("SX3 Z") + axes[1].set_ylabel("Energy Loss (MeV)") + + fig.colorbar(h1[3], ax=axes[0], label="Counts") + fig.colorbar(h2[3], ax=axes[1], label="Counts") + + plt.tight_layout() + + plt.savefig(f"{outdir}/Elost_vs_sx3_comparison.png", dpi=300) + plt.show() + except: + print("Error with side-by-side plots") #EA vs Esx3 overlay scatter plt.figure(figsize=(8,6)) - - plt.scatter( - data1["EA"], - data1["Esx3"], - s=1, - alpha=0.3, - label=data1["particle"] - ) - - plt.scatter( - data2["EA"], - data2["Esx3"], - s=1, - alpha=0.3, - label=data2["particle"] - ) - + plt.scatter(data1["EA"],data1["Esx3"],s=1,alpha=0.3,label=data1["particle"]) + plt.scatter(data2["EA"],data2["Esx3"],s=1,alpha=0.3,label=data2["particle"]) plt.xlabel("EA (MeV)") plt.ylabel("Esx3 (MeV)") plt.title("Anode vs SX3 Energy") plt.legend() plt.grid(True) plt.tight_layout() - plt.savefig(f"{outdir}/EA_vs_Esx3_overlay.png", dpi=300) plt.show() #PCE vs SiE - mask1 = data1["Esx3"] > 0 - mask2 = data2["Esx3"] > 0 + mask1 = data1["Esx3"] > 1 + mask2 = data2["Esx3"] > 1 + + thetab1 = np.deg2rad(data1["thetab"][mask1]) + thetab2 = np.deg2rad(data2["thetab"][mask2]) + + #theta smear (spatial uncertainty) + if True: + sigma = 0.2 * np.sqrt(np.maximum(thetab1, 0)) + + thetab1 += np.random.normal(0,sigma,len(thetab1)) + thetab1 += np.random.normal(0,0.02*np.sqrt(thetab1),len(thetab1)) + + sigma = 0.2 * np.sqrt(np.maximum(thetab2, 0)) + thetab2 += np.random.normal(0,sigma,len(thetab2)) + thetab2 += np.random.normal(0,0.02*np.sqrt(thetab2),len(thetab2)) + combined_Esx3 = np.concatenate([ (data1["Esx3"][mask1]), data2["Esx3"][mask2] ]) combined_Eprop = np.concatenate([ - data1["Eprop"][mask1] * 3, - data2["Eprop"][mask2] * 2.7 + data1["Eprop"][mask1] * np.sin(thetab1) * 3, + data2["Eprop"][mask2] * np.sin(thetab2) * 3 ]) combined_Esx3 = ( @@ -1436,7 +1360,8 @@ class MyInteractiveApp(cmd.Cmd): + np.random.normal(0,0.08,len(combined_Esx3)) ) - if True: + #Artifical smear + if False: sigma = 0.02 * np.sqrt(np.maximum(combined_Eprop, 0)) combined_Eprop += np.random.normal( @@ -1451,13 +1376,9 @@ class MyInteractiveApp(cmd.Cmd): len(combined_Eprop) ) - mask = ( - np.isfinite(combined_Esx3) - & np.isfinite(combined_Eprop) - ) - + mask = (np.isfinite(combined_Esx3)& np.isfinite(combined_Eprop)) combined_Esx3 = combined_Esx3[mask] - combined_Eprop = combined_Eprop[mask]# * 3 + combined_Eprop = combined_Eprop[mask] #combined_Eprop = combined_Eprop * .686 @@ -1468,27 +1389,27 @@ class MyInteractiveApp(cmd.Cmd): bins=200 ) plt.xlabel("SX3 Energy (MeV)") - plt.ylabel("PCEnergy") - plt.xlim(0,30) - plt.ylim(0,0.45) + plt.ylabel("PCEnergy x Sin(theta)") + #plt.xlim(0,30) + #plt.ylim(0,0.45) plt.colorbar(label="Counts") plt.tight_layout() plt.savefig(f"{outdir}/Eprop_vs_Esx3.png", dpi=300) plt.show() - plt.figure(figsize=(11,6), facecolor='white') + plt.figure(figsize=(14,6), facecolor='white') plt.hist2d( combined_Esx3, combined_Eprop, bins=[500,500], #range=[[0,35],[0,0.6]], cmap=yellow_jet, - cmin=1 - ) + cmin=1) + plt.margins(0) plt.xlabel("SX3 Energy (MeV)") - plt.ylabel("PCEnergy") - plt.xlim(0,30) - plt.ylim(0,0.45) + plt.ylabel("PCEnergy x Sin(theta)") + plt.xlim(0,35) + plt.ylim(0,0.6) cbar = plt.colorbar() cbar.set_label("Counts") plt.minorticks_on() @@ -1513,14 +1434,14 @@ class MyInteractiveApp(cmd.Cmd): mask = data1["Esx3"] > 0 power_fit_and_plot( data1["Esx3"][mask], - data1["Eprop"][mask] * .686, + data1["Eprop"][mask], label=data1["particle"], color='red' ) mask2 = data2["Esx3"] > 0 power_fit_and_plot( data2["Esx3"][mask2], - data2["Eprop"][mask2] * .686, + data2["Eprop"][mask2], label=data2["particle"], color='cyan' ) @@ -1539,11 +1460,8 @@ class MyInteractiveApp(cmd.Cmd): plt.savefig(f"{outdir}/Combined_Eprop_vs_Esx3_fit.png",dpi=300) plt.show() - - print("Dual plotting complete.") - - + #exec(open("PCEnergyAnalysis.py").read()) if __name__ == "__main__":