diff --git a/lstchain/data/LST1_backg_irf_gheffi_0.40_theffi_0.40.csv b/lstchain/data/LST1_backg_irf_gheffi_0.40_theffi_0.40.csv new file mode 100644 index 0000000000..b1f1986950 --- /dev/null +++ b/lstchain/data/LST1_backg_irf_gheffi_0.40_theffi_0.40.csv @@ -0,0 +1,191 @@ +ZD_deg,Ereco_min_TeV,Ereco_max_TeV,BckgRate_per_second,Theta_cut_deg,Gammaness_cut +6.0,0.012589254117941675,0.0199526231496888,0.08517203990519402,0.32,0.14980508270919676 +6.0,0.0199526231496888,0.03162277660168379,3.4842739444524904,0.25416790232400216,0.2232647189688829 +6.0,0.03162277660168379,0.05011872336272725,4.478140055010664,0.22836642194149148,0.24282859738626 +6.0,0.05011872336272725,0.07943282347242814,1.4772598658268499,0.17930820879898315,0.3082568921420486 +6.0,0.07943282347242814,0.12589254117941673,0.1278754571286665,0.11869128089021427,0.5205997890139312 +6.0,0.12589254117941673,0.19952623149688797,0.015107381355225944,0.1,0.6802999880849291 +6.0,0.19952623149688797,0.31622776601683794,0.0049059443594680376,0.1,0.7519160350778854 +6.0,0.31622776601683794,0.5011872336272722,0.0016774052026938128,0.1,0.8187521024602308 +6.0,0.5011872336272722,0.7943282347242814,0.0007139553959658439,0.1,0.8703749276810416 +6.0,0.7943282347242814,1.2589254117941668,0.00030388143938634826,0.1,0.9016936131710661 +6.0,1.2589254117941668,1.9952623149688788,0.0001293413142127672,0.1,0.9316904703903903 +6.0,1.9952623149688788,3.1622776601683795,5.505165302648401e-05,0.1,0.9396606170144983 +6.0,3.1622776601683795,5.01187233627272,2.343168166640779e-05,0.1,0.9500643114259029 +6.0,5.01187233627272,7.943282347242813,9.973246497280277e-06,0.1,0.9510363823031566 +6.0,7.943282347242813,12.589254117941662,4.244921346729869e-06,0.1,0.9401611478282427 +6.0,12.589254117941662,19.952623149688787,1.8067694651723323e-06,0.1,0.9567406906007412 +6.0,19.952623149688787,31.622776601683793,7.947470144435928e-07,0.1016591644784053,0.9477693994941487 +6.0,31.622776601683793,50.118723362727195,3.273172783976941e-07,0.1,0.9310066460490418 +6.0,50.118723362727195,79.43282347242814,1.42659587511951e-06,0.32,0.1 +9.579,0.012589254117941675,0.0199526231496888,0.07496221511832982,0.32,0.14954631782891498 +9.579,0.0199526231496888,0.03162277660168379,3.1143848906358587,0.2489388571568171,0.21880920839480064 +9.579,0.03162277660168379,0.05011872336272725,4.613326231585207,0.23467156439645145,0.24001297834835228 +9.579,0.05011872336272725,0.07943282347242814,1.4503731886193367,0.1782758956117673,0.30538889532484204 +9.579,0.07943282347242814,0.12589254117941673,0.12516326928576404,0.11670020873040216,0.5199252262542218 +9.579,0.12589254117941673,0.19952623149688797,0.015784489258984336,0.1,0.6741285296388073 +9.579,0.19952623149688797,0.31622776601683794,0.005171553265496122,0.1,0.7502508793627872 +9.579,0.31622776601683794,0.5011872336272722,0.001865983854968374,0.1,0.8134164858767834 +9.579,0.5011872336272722,0.7943282347242814,0.0007358633567411192,0.1,0.8706827596053548 +9.579,0.7943282347242814,1.2589254117941668,0.00029019269290702694,0.1,0.903105147449207 +9.579,1.2589254117941668,1.9952623149688788,0.00011443945162533516,0.1,0.9316411475681786 +9.579,1.9952623149688788,3.1622776601683795,4.512997194075903e-05,0.1,0.9390698016445345 +9.579,3.1622776601683795,5.01187233627272,1.77973097428125e-05,0.1,0.9568974012918506 +9.579,5.01187233627272,7.943282347242813,7.018489497343159e-06,0.1,0.9571810800605608 +9.579,7.943282347242813,12.589254117941662,2.7677888139363144e-06,0.1,0.9550247825644098 +9.579,12.589254117941662,19.952623149688787,1.0914962430948887e-06,0.1,0.9537668221860696 +9.579,19.952623149688787,31.622776601683793,4.3043892752637916e-07,0.1,0.9466113932901232 +9.579,31.622776601683793,50.118723362727195,1.6974650302479555e-07,0.1,0.9113429102140364 +9.579,50.118723362727195,79.43282347242814,6.854709739886736e-07,0.32,0.1 +16.08,0.012589254117941675,0.0199526231496888,0.01999735615095646,0.23587899012290253,0.117942388934193 +16.08,0.0199526231496888,0.03162277660168379,2.302646457107055,0.2514035018042374,0.20694299577277286 +16.08,0.03162277660168379,0.05011872336272725,4.623237882761817,0.23422839772906243,0.23553418601607756 +16.08,0.05011872336272725,0.07943282347242814,1.8006793896081452,0.18845985664607448,0.2899963008856928 +16.08,0.07943282347242814,0.12589254117941673,0.1713798825397484,0.1208416089315547,0.4954425122131275 +16.08,0.12589254117941673,0.19952623149688797,0.018924926393187574,0.1,0.6704718083416198 +16.08,0.19952623149688797,0.31622776601683794,0.005676468690662127,0.1,0.7539795241681223 +16.08,0.31622776601683794,0.5011872336272722,0.0022216219399217463,0.1,0.8072247840435847 +16.08,0.5011872336272722,0.7943282347242814,0.000811894918680494,0.1,0.870828522852448 +16.08,0.7943282347242814,1.2589254117941668,0.0002967081604363452,0.1,0.9032104519757597 +16.08,1.2589254117941668,1.9952623149688788,0.00010843242203387257,0.1,0.9306896639956173 +16.08,1.9952623149688788,3.1622776601683795,3.962678387692773e-05,0.1,0.9458204733612634 +16.08,3.1622776601683795,5.01187233627272,1.4481664902202484e-05,0.1,0.9534106981487875 +16.08,5.01187233627272,7.943282347242813,5.292345172170021e-06,0.1,0.9529725896325841 +16.08,7.943282347242813,12.589254117941662,1.9340951203153105e-06,0.1,0.9569630171986361 +16.08,12.589254117941662,19.952623149688787,7.068178307980026e-07,0.1,0.9665766973249578 +16.08,19.952623149688787,31.622776601683793,2.5830758822893196e-07,0.1,0.9487771937769978 +16.08,31.622776601683793,50.118723362727195,9.439887794188343e-08,0.1,0.9491118247571814 +16.08,50.118723362727195,79.43282347242814,3.532607915915502e-07,0.32,0.1 +23.16,0.012589254117941675,0.0199526231496888,0.0037997870728821463,0.32,0.1 +23.16,0.0199526231496888,0.03162277660168379,1.030311118455608,0.2578799183615753,0.18664904409216268 +23.16,0.03162277660168379,0.05011872336272725,4.067640247520111,0.23703638782089045,0.22665335685425472 +23.16,0.05011872336272725,0.07943282347242814,2.3476423010362364,0.20736984708463185,0.27127577525471624 +23.16,0.07943282347242814,0.12589254117941673,0.2974774476053728,0.13312368458559823,0.4363628788022292 +23.16,0.12589254117941673,0.19952623149688797,0.02446630431756125,0.1,0.651194465305772 +23.16,0.19952623149688797,0.31622776601683794,0.006437754261794234,0.1,0.746140744781241 +23.16,0.31622776601683794,0.5011872336272722,0.002394668248952174,0.1,0.8030298811164628 +23.16,0.5011872336272722,0.7943282347242814,0.0009582592202806372,0.1,0.8551003909714278 +23.16,0.7943282347242814,1.2589254117941668,0.0003834605205354247,0.1,0.8943986180625971 +23.16,1.2589254117941668,1.9952623149688788,0.0001534469668512409,0.1,0.9227229926235273 +23.16,1.9952623149688788,3.1622776601683795,6.140390046664688e-05,0.1,0.941335257607903 +23.16,3.1622776601683795,5.01187233627272,2.4571609787328856e-05,0.1,0.950615305688041 +23.16,5.01187233627272,7.943282347242813,9.832665399956234e-06,0.1,0.9609170501847444 +23.16,7.943282347242813,12.589254117941662,3.934675412164217e-06,0.1,0.9522570871256878 +23.16,12.589254117941662,19.952623149688787,1.574514129114834e-06,0.1,0.9572631795771868 +23.16,19.952623149688787,31.622776601683793,6.300633427392808e-07,0.1,0.9602385697244247 +23.16,31.622776601683793,50.118723362727195,2.521284557077757e-07,0.1,0.9462835716657525 +23.16,50.118723362727195,79.43282347242814,1.0331382779400924e-06,0.32,0.1 +30.39,0.012589254117941675,0.0199526231496888,0.00044520408537425323,0.32,0.1 +30.39,0.0199526231496888,0.03162277660168379,0.44202683877235854,0.2929976786720368,0.13606157250960138 +30.39,0.03162277660168379,0.05011872336272725,3.0676733282995214,0.23867266446185426,0.21709321056056893 +30.39,0.05011872336272725,0.07943282347242814,2.950966529004954,0.2187740774314825,0.25822369493169384 +30.39,0.07943282347242814,0.12589254117941673,0.5771193075206469,0.15255605129139066,0.37597735266542565 +30.39,0.12589254117941673,0.19952623149688797,0.04099737940259154,0.10197576489596409,0.6213120756069364 +30.39,0.19952623149688797,0.31622776601683794,0.008205667350469301,0.1,0.732327003039153 +30.39,0.31622776601683794,0.5011872336272722,0.0030436903367525143,0.1,0.7930150727145128 +30.39,0.5011872336272722,0.7943282347242814,0.0011845587191292868,0.1,0.8378616050145095 +30.39,0.7943282347242814,1.2589254117941668,0.0004610125222404682,0.1,0.8843225879167961 +30.39,1.2589254117941668,1.9952623149688788,0.00017941917292098515,0.1,0.9131757789367785 +30.39,1.9952623149688788,3.1622776601683795,6.982725643807814e-05,0.1,0.9355204655689834 +30.39,3.1622776601683795,5.01187233627272,2.717572298595095e-05,0.1,0.9384825559019611 +30.39,5.01187233627272,7.943282347242813,1.05763846022513e-05,0.1,0.9434221550172315 +30.39,7.943282347242813,12.589254117941662,4.116170572998801e-06,0.1,0.9443956027898699 +30.39,12.589254117941662,19.952623149688787,1.601951973495253e-06,0.1,0.9307896797127482 +30.39,19.952623149688787,31.622776601683793,6.23455729026243e-07,0.1,0.9449926556277641 +30.39,31.622776601683793,50.118723362727195,2.426396374465317e-07,0.1,0.9421297091509299 +30.39,50.118723362727195,79.43282347242814,9.443171490642629e-08,0.1,0.9321790193703262 +37.66,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +37.66,0.0199526231496888,0.03162277660168379,0.004903530549231421,0.32,0.1 +37.66,0.03162277660168379,0.05011872336272725,0.9298672248422799,0.24531484308734686,0.1849816184719492 +37.66,0.05011872336272725,0.07943282347242814,3.1541836942365538,0.23426394546475227,0.24395674538842407 +37.66,0.07943282347242814,0.12589254117941673,1.2907247948224125,0.18411193053790204,0.3126474920499498 +37.66,0.12589254117941673,0.19952623149688797,0.11924190967344299,0.11462927210397188,0.5245377824735699 +37.66,0.19952623149688797,0.31622776601683794,0.013637264534758217,0.1,0.6941997622809297 +37.66,0.31622776601683794,0.5011872336272722,0.004350075891330023,0.1,0.7573834462732065 +37.66,0.5011872336272722,0.7943282347242814,0.0015857260793042968,0.1,0.8143646888213829 +37.66,0.7943282347242814,1.2589254117941668,0.0006328318630421753,0.1,0.8601806857224763 +37.66,1.2589254117941668,1.9952623149688788,0.0002525506593529261,0.1,0.8847795693551539 +37.66,1.9952623149688788,3.1622776601683795,0.00010078796480471615,0.1,0.906062314172316 +37.66,3.1622776601683795,5.01187233627272,4.022248001847874e-05,0.1,0.9301432362329947 +37.66,5.01187233627272,7.943282347242813,1.6051994917960855e-05,0.1,0.9375599988109555 +37.66,7.943282347242813,12.589254117941662,6.406033161750984e-06,0.1,0.9425441734747769 +37.66,12.589254117941662,19.952623149688787,2.5565209233611213e-06,0.1,0.9426310939001388 +37.66,19.952623149688787,31.622776601683793,1.0202568526505638e-06,0.1,0.9454227432274414 +37.66,31.622776601683793,50.118723362727195,4.0716429733417003e-07,0.1,0.9476555141986623 +37.66,50.118723362727195,79.43282347242814,1.62491204634338e-07,0.1,0.9558454031429453 +44.92,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +44.92,0.0199526231496888,0.03162277660168379,0.0005034818183207038,0.32,0.1 +44.92,0.03162277660168379,0.05011872336272725,0.3481495825351589,0.31056843122863537,0.1 +44.92,0.05011872336272725,0.07943282347242814,1.7417044364784289,0.22870019947192644,0.20055015744422341 +44.92,0.07943282347242814,0.12589254117941673,2.0195832469545394,0.21188986854704878,0.26669948683652384 +44.92,0.12589254117941673,0.19952623149688797,0.3960863401859151,0.14211675071190638,0.3894695543197544 +44.92,0.19952623149688797,0.31622776601683794,0.03370043812298416,0.1,0.6319406502704973 +44.92,0.31622776601683794,0.5011872336272722,0.0070211845401037904,0.1,0.7401805275680142 +44.92,0.5011872336272722,0.7943282347242814,0.002246864365940871,0.1,0.7995393351331342 +44.92,0.7943282347242814,1.2589254117941668,0.0009679142605008971,0.1,0.8319093411166305 +44.92,1.2589254117941668,1.9952623149688788,0.00041696242545049693,0.1,0.8629713704659066 +44.92,1.9952623149688788,3.1622776601683795,0.00017962093475881803,0.1,0.879029958386691 +44.92,3.1622776601683795,5.01187233627272,7.737790801838572e-05,0.1,0.9032349091791307 +44.92,5.01187233627272,7.943282347242813,3.3333200594580627e-05,0.1,0.9167087730389729 +44.92,7.943282347242813,12.589254117941662,1.4359424935791001e-05,0.1,0.9147302906981755 +44.92,12.589254117941662,19.952623149688787,6.185817167528149e-06,0.1,0.9042515734100918 +44.92,19.952623149688787,31.622776601683793,2.664753929992821e-06,0.1,0.9145788316532099 +44.92,31.622776601683793,50.118723362727195,1.1479345921647578e-06,0.1,0.9097760484559024 +44.92,50.118723362727195,79.43282347242814,4.94512387450356e-07,0.1,0.9170924616126577 +52.16,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +52.16,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +52.16,0.03162277660168379,0.05011872336272725,0.0019554053555110216,0.32,0.1 +52.16,0.05011872336272725,0.07943282347242814,0.11696709028698668,0.2286017877369897,0.1 +52.16,0.07943282347242814,0.12589254117941673,1.5437040284343415,0.20882688124698462,0.195630797643331 +52.16,0.12589254117941673,0.19952623149688797,1.5350369090560867,0.18948531016313416,0.27962487191709734 +52.16,0.19952623149688797,0.31622776601683794,0.2103482496969563,0.1271152353053556,0.45899875297284454 +52.16,0.31622776601683794,0.5011872336272722,0.017626402866680133,0.1,0.6869338523202627 +52.16,0.5011872336272722,0.7943282347242814,0.004105712193174192,0.1,0.7640496741110404 +52.16,0.7943282347242814,1.2589254117941668,0.001638159631918383,0.1,0.809011634949459 +52.16,1.2589254117941668,1.9952623149688788,0.000743727554210192,0.1,0.8312713822693032 +52.16,1.9952623149688788,3.1622776601683795,0.0003376537085361608,0.1,0.8560218568749494 +52.16,3.1622776601683795,5.01187233627272,0.00015329541878987733,0.1,0.8610640097240986 +52.16,5.01187233627272,7.943282347242813,6.959640847376393e-05,0.1,0.8858931169655582 +52.16,7.943282347242813,12.589254117941662,3.159690035542566e-05,0.1,0.9025033051679945 +52.16,12.589254117941662,19.952623149688787,1.4345052193994404e-05,0.1,0.894107880719096 +52.16,19.952623149688787,31.622776601683793,6.512680678599796e-06,0.1,0.8862347667752004 +52.16,31.622776601683793,50.118723362727195,2.956769278202017e-06,0.1,0.8802914243877077 +52.16,50.118723362727195,79.43282347242814,1.3423788138801976e-06,0.1,0.893212561651399 +59.34,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +59.34,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +59.34,0.03162277660168379,0.05011872336272725,0.0,0.32,0.1 +59.34,0.05011872336272725,0.07943282347242814,0.002694531473847235,0.32,0.1 +59.34,0.07943282347242814,0.12589254117941673,0.5574034753482935,0.32,0.1 +59.34,0.12589254117941673,0.19952623149688797,0.8363832880051874,0.18896112393783765,0.16923710730878283 +59.34,0.19952623149688797,0.31622776601683794,0.9507893227399531,0.17708165358021902,0.2640702637641676 +59.34,0.31622776601683794,0.5011872336272722,0.16225121830657654,0.12301084561923045,0.4675715358796556 +59.34,0.5011872336272722,0.7943282347242814,0.02001906029749398,0.1,0.7001198865478234 +59.34,0.7943282347242814,1.2589254117941668,0.004804964695608002,0.1,0.7659993562236929 +59.34,1.2589254117941668,1.9952623149688788,0.0014958161500865784,0.1,0.7939874191858363 +59.34,1.9952623149688788,3.1622776601683795,0.000640994817720803,0.1,0.8175957633965458 +59.34,3.1622776601683795,5.01187233627272,0.0002746823908280065,0.1,0.8121226351902138 +59.34,5.01187233627272,7.943282347242813,0.00011770830862450671,0.1,0.8266512716247353 +59.34,7.943282347242813,12.589254117941662,5.0440968849428864e-05,0.1,0.8292790830485759 +59.34,12.589254117941662,19.952623149688787,2.1615222988085098e-05,0.1,0.8293737894005931 +59.34,19.952623149688787,31.622776601683793,9.262666350032505e-06,0.1,0.8354121746337494 +59.34,31.622776601683793,50.118723362727195,3.9692853485396905e-06,0.1,0.8409636718231623 +59.34,50.118723362727195,79.43282347242814,1.700938540021639e-06,0.1,0.8346364329941128 +66.44,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +66.44,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +66.44,0.03162277660168379,0.05011872336272725,0.0,0.32,0.1 +66.44,0.05011872336272725,0.07943282347242814,0.0,0.32,0.1 +66.44,0.07943282347242814,0.12589254117941673,0.0001791018992365193,0.32,0.1 +66.44,0.12589254117941673,0.19952623149688797,0.3756631456848193,0.32,0.1 +66.44,0.19952623149688797,0.31622776601683794,0.378698120932947,0.18430528156401446,0.1 +66.44,0.31622776601683794,0.5011872336272722,1.1540511175333352,0.17010809252661016,0.18907180317960973 +66.44,0.5011872336272722,0.7943282347242814,0.36330045784913123,0.13246925363741566,0.3423407488920865 +66.44,0.7943282347242814,1.2589254117941668,0.026877509326041583,0.1,0.6129610320422042 +66.44,1.2589254117941668,1.9952623149688788,0.0036020326669464198,0.1,0.7594603806998884 +66.44,1.9952623149688788,3.1622776601683795,0.0013513455207365842,0.1,0.7772482091153307 +66.44,3.1622776601683795,5.01187233627272,0.0006775027183780024,0.1,0.7653316236510899 +66.44,5.01187233627272,7.943282347242813,0.00033966881627682326,0.1,0.7497217406946608 +66.44,7.943282347242813,12.589254117941662,0.00017029437907956354,0.1,0.7345757077104096 +66.44,12.589254117941662,19.952623149688787,8.537779789140115e-05,0.1,0.7553842611535082 +66.44,19.952623149688787,31.622776601683793,4.280451540552178e-05,0.1,0.7705226017247755 +66.44,31.622776601683793,50.118723362727195,2.1460222497564374e-05,0.1,0.7916358634248051 +66.44,50.118723362727195,79.43282347242814,1.0759172140646597e-05,0.1,0.8106709918015099 diff --git a/lstchain/data/LST1_backg_irf_gheffi_0.70_theffi_0.70.csv b/lstchain/data/LST1_backg_irf_gheffi_0.70_theffi_0.70.csv new file mode 100644 index 0000000000..f3e2ba3c71 --- /dev/null +++ b/lstchain/data/LST1_backg_irf_gheffi_0.70_theffi_0.70.csv @@ -0,0 +1,191 @@ +ZD_deg,Ereco_min_TeV,Ereco_max_TeV,BckgRate_per_second,Theta_cut_deg,Gammaness_cut +6.0,0.012589254117941675,0.0199526231496888,0.15514275685330334,0.32,0.11281024702481211 +6.0,0.0199526231496888,0.03162277660168379,11.05469622088641,0.32,0.17608067775290787 +6.0,0.03162277660168379,0.05011872336272725,17.509024798429156,0.32,0.18909915303573652 +6.0,0.05011872336272725,0.07943282347242814,10.888996161065874,0.32,0.19761519995256535 +6.0,0.07943282347242814,0.12589254117941673,2.0898779630778845,0.2673516001927132,0.31688826851906404 +6.0,0.12589254117941673,0.19952623149688797,0.23999615011858383,0.21162990524807224,0.4483405640793372 +6.0,0.19952623149688797,0.31622776601683794,0.07211030101126466,0.1964156738356117,0.5117977755279097 +6.0,0.31622776601683794,0.5011872336272722,0.02296676261259157,0.1793650747665875,0.5752482874791648 +6.0,0.5011872336272722,0.7943282347242814,0.010091305532878039,0.1739656472302629,0.6283677225961126 +6.0,0.7943282347242814,1.2589254117941668,0.004284656279599664,0.17348059554416378,0.6651146560298868 +6.0,1.2589254117941668,1.9952623149688788,0.0015450244710438553,0.17624772620116916,0.724208825217094 +6.0,1.9952623149688788,3.1622776601683795,0.0005477990781407422,0.1775535739379539,0.7674645252753842 +6.0,3.1622776601683795,5.01187233627272,0.00018103057367059957,0.17268620694592587,0.7940713722648308 +6.0,5.01187233627272,7.943282347242813,6.58904018964151e-05,0.17626078400309486,0.8155952548324541 +6.0,7.943282347242813,12.589254117941662,2.607448110457876e-05,0.1875924714679997,0.8099803960631745 +6.0,12.589254117941662,19.952623149688787,7.258738446189413e-06,0.16745612152278178,0.8145827139310521 +6.0,19.952623149688787,31.622776601683793,3.7452352339724327e-06,0.2035040253875547,0.8166966955524747 +6.0,31.622776601683793,50.118723362727195,5.162808422191587e-07,0.1278319007462821,0.7245412020766999 +6.0,50.118723362727195,79.43282347242814,1.1302660241562094e-06,0.32,0.1 +9.579,0.012589254117941675,0.0199526231496888,0.13453279350789224,0.32,0.1062672796107913 +9.579,0.0199526231496888,0.03162277660168379,10.261078547984237,0.32,0.17180929867983782 +9.579,0.03162277660168379,0.05011872336272725,17.150348733802918,0.32,0.18601167944133717 +9.579,0.05011872336272725,0.07943282347242814,10.708411006873446,0.32,0.19501246351169796 +9.579,0.07943282347242814,0.12589254117941673,2.020196695528386,0.26249644100706687,0.31576995492959054 +9.579,0.12589254117941673,0.19952623149688797,0.2493115299814993,0.21260796591917902,0.43624271737327935 +9.579,0.19952623149688797,0.31622776601683794,0.07302200995517918,0.1959237053617897,0.503618439468562 +9.579,0.31622776601683794,0.5011872336272722,0.025647282916292518,0.18304408518215354,0.5621554621054391 +9.579,0.5011872336272722,0.7943282347242814,0.01032639074384704,0.1745120474721881,0.6247461567451729 +9.579,0.7943282347242814,1.2589254117941668,0.004707475782640064,0.17364226573653765,0.6546501202021747 +9.579,1.2589254117941668,1.9952623149688788,0.0017924382704006246,0.17598988293186216,0.7246773797513005 +9.579,1.9952623149688788,3.1622776601683795,0.0006797671057046483,0.17801224280917452,0.7737687181061742 +9.579,3.1622776601683795,5.01187233627272,0.00023227549530653382,0.17091322945173035,0.8160090457289869 +9.579,5.01187233627272,7.943282347242813,7.678671181623945e-05,0.16140656603972997,0.8276962270244919 +9.579,7.943282347242813,12.589254117941662,3.87924193756367e-05,0.18843251513484596,0.8094454349142098 +9.579,12.589254117941662,19.952623149688787,1.3803948604211396e-05,0.18462404760506107,0.8264931239980831 +9.579,19.952623149688787,31.622776601683793,3.5870040776471747e-06,0.1545812010486962,0.8204809696921489 +9.579,31.622776601683793,50.118723362727195,9.561222500539285e-07,0.1310845823395837,0.6623916168229874 +9.579,50.118723362727195,79.43282347242814,2.112033577641232e-06,0.32,0.1 +16.08,0.012589254117941675,0.0199526231496888,0.05525288706530069,0.32,0.1 +16.08,0.0199526231496888,0.03162277660168379,7.395023010988047,0.32,0.15871012759309444 +16.08,0.03162277660168379,0.05011872336272725,17.27044593826248,0.32,0.18455200504785835 +16.08,0.05011872336272725,0.07943282347242814,11.527288478989545,0.32,0.18949670544837852 +16.08,0.07943282347242814,0.12589254117941673,2.689946864017896,0.2724281634454505,0.2978261008792913 +16.08,0.12589254117941673,0.19952623149688797,0.29336840531118696,0.21098096755684217,0.43338077520875984 +16.08,0.19952623149688797,0.31622776601683794,0.07581084902102551,0.19115748735220828,0.5095112328716478 +16.08,0.31622776601683794,0.5011872336272722,0.0302904455773565,0.18628273923337543,0.5558435714740266 +16.08,0.5011872336272722,0.7943282347242814,0.011423725561458403,0.17555909293551242,0.6205931081610208 +16.08,0.7943282347242814,1.2589254117941668,0.0048865481631772964,0.17065888688454617,0.6548541467630675 +16.08,1.2589254117941668,1.9952623149688788,0.0021198674818263294,0.17544369021121614,0.719944269328782 +16.08,1.9952623149688788,3.1622776601683795,0.0009203929316937314,0.1804370073088594,0.7669775580574933 +16.08,3.1622776601683795,5.01187233627272,0.00038783237760301245,0.1828170030516919,0.7939089282972858 +16.08,5.01187233627272,7.943282347242813,0.0001574857493387026,0.18183218337620444,0.8206511143622254 +16.08,7.943282347242813,12.589254117941662,6.764499106605881e-05,0.18600453137256687,0.8325706439063199 +16.08,12.589254117941662,19.952623149688787,2.7833109367250163e-05,0.18622678135939866,0.8248789734517062 +16.08,19.952623149688787,31.622776601683793,7.848024094935826e-06,0.1543465581252543,0.8229067850806525 +16.08,31.622776601683793,50.118723362727195,1.99038240209587e-06,0.12132228128821079,0.813707871649978 +16.08,50.118723362727195,79.43282347242814,5.683858204287456e-06,0.32,0.1 +23.16,0.012589254117941675,0.0199526231496888,0.0041728407715988505,0.32,0.1 +23.16,0.0199526231496888,0.03162277660168379,3.1556705725517804,0.32,0.13287595749300685 +23.16,0.03162277660168379,0.05011872336272725,14.97632756139708,0.32,0.17674200189280798 +23.16,0.05011872336272725,0.07943282347242814,12.058594769952464,0.32,0.18683194617380808 +23.16,0.07943282347242814,0.12589254117941673,4.279203297691901,0.2980711569005066,0.27012102227819773 +23.16,0.12589254117941673,0.19952623149688797,0.41607744055349044,0.21999715435848038,0.41726404044314613 +23.16,0.19952623149688797,0.31622776601683794,0.08695029642188,0.19257474848592668,0.5016088022805569 +23.16,0.31622776601683794,0.5011872336272722,0.03184215274624511,0.18316071004588494,0.5505978790502144 +23.16,0.5011872336272722,0.7943282347242814,0.014351167983589378,0.17805703472453632,0.5826247541059284 +23.16,0.7943282347242814,1.2589254117941668,0.005869436091868877,0.17616453857941652,0.6470008848062724 +23.16,1.2589254117941668,1.9952623149688788,0.0025380630886366385,0.17228031578013675,0.7079716572463909 +23.16,1.9952623149688788,3.1622776601683795,0.0011228602955282321,0.1704164327776217,0.7491467671786793 +23.16,3.1622776601683795,5.01187233627272,0.0005682547621178823,0.18029526744332122,0.7793540495073453 +23.16,5.01187233627272,7.943282347242813,0.0002662169457879666,0.1835248164793628,0.8102841198456221 +23.16,7.943282347242813,12.589254117941662,0.00011751565826497445,0.18133806327636143,0.7984871443704252 +23.16,12.589254117941662,19.952623149688787,5.505480552198517e-05,0.1845877541112928,0.843535208980421 +23.16,19.952623149688787,31.622776601683793,1.909142995166195e-05,0.16165480321673856,0.8289874768781574 +23.16,31.622776601683793,50.118723362727195,6.300807295130527e-06,0.13811205978462782,0.8042892210088163 +23.16,50.118723362727195,79.43282347242814,1.5293388236766216e-05,0.32,0.1 +30.39,0.012589254117941675,0.0199526231496888,0.00044520408537425323,0.32,0.1 +30.39,0.0199526231496888,0.03162277660168379,0.9155897850387641,0.32,0.1 +30.39,0.03162277660168379,0.05011872336272725,11.10242449502949,0.32,0.1639943823130124 +30.39,0.05011872336272725,0.07943282347242814,13.120433772840112,0.32,0.18995715355916115 +30.39,0.07943282347242814,0.12589254117941673,6.768787991435292,0.32,0.23979310406872537 +30.39,0.12589254117941673,0.19952623149688797,0.6959854235423557,0.22402461427403417,0.40273486866668845 +30.39,0.19952623149688797,0.31622776601683794,0.10936131682439872,0.19288798399745657,0.5027879408643803 +30.39,0.31622776601683794,0.5011872336272722,0.03734158644717221,0.17883852724172133,0.5563226609339352 +30.39,0.5011872336272722,0.7943282347242814,0.01689205525244843,0.17493639469943909,0.5803931952439499 +30.39,0.7943282347242814,1.2589254117941668,0.007091993868415975,0.16971032286194984,0.6276184481944082 +30.39,1.2589254117941668,1.9952623149688788,0.0031242303490547625,0.16715020923254043,0.6812090477125519 +30.39,1.9952623149688788,3.1622776601683795,0.001462909122026292,0.17134919204561602,0.7369775023835263 +30.39,3.1622776601683795,5.01187233627272,0.0007940581277243386,0.18912003866642446,0.7489922589629406 +30.39,5.01187233627272,7.943282347242813,0.0003002067912910604,0.1742045944084464,0.775382140614782 +30.39,7.943282347242813,12.589254117941662,0.00014876597858109946,0.18371269022815082,0.7964947089057927 +30.39,12.589254117941662,19.952623149688787,8.368386133453983e-05,0.20641736395150476,0.7930503688994939 +30.39,19.952623149688787,31.622776601683793,3.0800929395405765e-05,0.18760535027670747,0.8230752954142034 +30.39,31.622776601683793,50.118723362727195,8.8033288982448e-06,0.15025352913817633,0.8310774813706188 +30.39,50.118723362727195,79.43282347242814,3.3560342374281397e-06,0.13898001973606883,0.9083706576565244 +37.66,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +37.66,0.0199526231496888,0.03162277660168379,0.005439726626267147,0.32,0.1 +37.66,0.03162277660168379,0.05011872336272725,3.3375735922542797,0.32,0.11722132626955144 +37.66,0.05011872336272725,0.07943282347242814,11.73090146393259,0.32,0.18230473862772834 +37.66,0.07943282347242814,0.12589254117941673,8.992373593228063,0.32,0.2089567583420639 +37.66,0.12589254117941673,0.19952623149688797,1.8131050615579232,0.2527672221001267,0.3395879542412909 +37.66,0.19952623149688797,0.31622776601683794,0.18577939629823192,0.19566267756040165,0.48358370727664507 +37.66,0.31622776601683794,0.5011872336272722,0.05656180643929607,0.1834213475404794,0.5241001139699583 +37.66,0.5011872336272722,0.7943282347242814,0.022124863563169547,0.17040517814542902,0.5451847511534684 +37.66,0.7943282347242814,1.2589254117941668,0.010311179893023055,0.16986375115070904,0.5692501405024707 +37.66,1.2589254117941668,1.9952623149688788,0.0053634770958229955,0.17073099613563647,0.5987840017675614 +37.66,1.9952623149688788,3.1622776601683795,0.002468614581690785,0.17116008788909587,0.6455422465080286 +37.66,3.1622776601683795,5.01187233627272,0.0011476689876657765,0.17245304908370357,0.6930177608297756 +37.66,5.01187233627272,7.943282347242813,0.0005680543914486575,0.17928511092161245,0.7351275510452387 +37.66,7.943282347242813,12.589254117941662,0.000242386379844833,0.1730573517272854,0.7672332972520108 +37.66,12.589254117941662,19.952623149688787,0.0001026359808893749,0.16640732885712728,0.7835020740940228 +37.66,19.952623149688787,31.622776601683793,5.369297807661835e-05,0.1778558255900631,0.7927161914703217 +37.66,31.622776601683793,50.118723362727195,2.4968087532488496e-05,0.179220913351681,0.7806032848847845 +37.66,50.118723362727195,79.43282347242814,6.854871039873389e-06,0.13876575937601177,0.7421657312459152 +44.92,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +44.92,0.0199526231496888,0.03162277660168379,0.0005034818183207038,0.32,0.1 +44.92,0.03162277660168379,0.05011872336272725,0.5086737255508019,0.32,0.1 +44.92,0.05011872336272725,0.07943282347242814,6.418129948766718,0.32,0.11897251766575202 +44.92,0.07943282347242814,0.12589254117941673,9.904341714888213,0.32,0.1837821792493808 +44.92,0.12589254117941673,0.19952623149688797,4.902337090843914,0.3024989789352129,0.2525458156292578 +44.92,0.19952623149688797,0.31622776601683794,0.5604466911148424,0.2121320161600014,0.4293445270705401 +44.92,0.31622776601683794,0.5011872336272722,0.08776583465231391,0.18041396202136545,0.5175175038429095 +44.92,0.5011872336272722,0.7943282347242814,0.026850582427664617,0.16526725365221795,0.5440768051602956 +44.92,0.7943282347242814,1.2589254117941668,0.012887162869378266,0.16097223287861506,0.5524953278466989 +44.92,1.2589254117941668,1.9952623149688788,0.007250183920382882,0.16371647071121967,0.5602967285967493 +44.92,1.9952623149688788,3.1622776601683795,0.0037801269469097875,0.16639429403474953,0.5875755949931851 +44.92,3.1622776601683795,5.01187233627272,0.001986945749602008,0.1681963148550561,0.6201994770798541 +44.92,5.01187233627272,7.943282347242813,0.0009890997311236382,0.16545570235079157,0.6666928181923851 +44.92,7.943282347242813,12.589254117941662,0.0005469293821485527,0.17154006394593865,0.6906162599772102 +44.92,12.589254117941662,19.952623149688787,0.0003444574615941554,0.1898042820273609,0.7281710124768749 +44.92,19.952623149688787,31.622776601683793,0.00019202748110055764,0.19758690018354005,0.7421912549100427 +44.92,31.622776601683793,50.118723362727195,9.134218318862706e-05,0.18999858388994198,0.7253371508623293 +44.92,50.118723362727195,79.43282347242814,3.657244211521825e-05,0.16762154445923916,0.7877360419394447 +52.16,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +52.16,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +52.16,0.03162277660168379,0.05011872336272725,0.0019554053555110216,0.32,0.1 +52.16,0.05011872336272725,0.07943282347242814,0.36170125871910935,0.32,0.1 +52.16,0.07943282347242814,0.12589254117941673,7.502925386172446,0.32,0.10652308487787197 +52.16,0.12589254117941673,0.19952623149688797,9.28456480272395,0.32,0.18473333371772352 +52.16,0.19952623149688797,0.31622776601683794,2.4957631418016617,0.2595946529302991,0.2959262029770815 +52.16,0.31622776601683794,0.5011872336272722,0.24780527925603688,0.1886623854150769,0.47686084897956843 +52.16,0.5011872336272722,0.7943282347242814,0.04808062092219167,0.1690366708729037,0.5263928545410059 +52.16,0.7943282347242814,1.2589254117941668,0.019024869246481526,0.15986004253069044,0.5290206247615351 +52.16,1.2589254117941668,1.9952623149688788,0.010147593472572753,0.15467617138729564,0.5075462889704585 +52.16,1.9952623149688788,3.1622776601683795,0.004946975878343615,0.1555130647391105,0.5203952668639182 +52.16,3.1622776601683795,5.01187233627272,0.003181292423022168,0.16069714049725675,0.5377443746126526 +52.16,5.01187233627272,7.943282347242813,0.0014446801650994203,0.15515195503284768,0.5759227203435799 +52.16,7.943282347242813,12.589254117941662,0.0007591774354520192,0.1611417616403506,0.6183355860929303 +52.16,12.589254117941662,19.952623149688787,0.0003964787770481321,0.16684434757697789,0.6528851954846248 +52.16,19.952623149688787,31.622776601683793,0.00022310549993776125,0.17931711428655356,0.6771689770598769 +52.16,31.622776601683793,50.118723362727195,0.00010457236571103917,0.17588958187093684,0.6855806258278904 +52.16,50.118723362727195,79.43282347242814,5.274901346304115e-05,0.17897974316562515,0.7148069911355314 +59.34,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +59.34,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +59.34,0.03162277660168379,0.05011872336272725,0.0,0.32,0.1 +59.34,0.05011872336272725,0.07943282347242814,0.002694531473847235,0.32,0.1 +59.34,0.07943282347242814,0.12589254117941673,0.7490480638793929,0.32,0.1 +59.34,0.12589254117941673,0.19952623149688797,4.574987635448585,0.3178846353573444,0.1 +59.34,0.19952623149688797,0.31622776601683794,6.20245097544997,0.30343081777545006,0.1613017852334765 +59.34,0.31622776601683794,0.5011872336272722,1.795317976820798,0.23905647977647984,0.2898554775852136 +59.34,0.5011872336272722,0.7943282347242814,0.22907801903668695,0.1757202337948652,0.4922966095686704 +59.34,0.7943282347242814,1.2589254117941668,0.05392950317835855,0.1587666456813387,0.5240522050270533 +59.34,1.2589254117941668,1.9952623149688788,0.016347751375240638,0.14965133240769898,0.49994863728429545 +59.34,1.9952623149688788,3.1622776601683795,0.007884754689762017,0.1441526636289405,0.4688845144612062 +59.34,3.1622776601683795,5.01187233627272,0.004671493509012499,0.14533876777071955,0.45309685949761846 +59.34,5.01187233627272,7.943282347242813,0.0025009838820515594,0.149213170331759,0.46354371638170516 +59.34,7.943282347242813,12.589254117941662,0.0015499149580913297,0.1586202192645271,0.4732790898387837 +59.34,12.589254117941662,19.952623149688787,0.0008498678511541781,0.15861096053452234,0.5008794110357272 +59.34,19.952623149688787,31.622776601683793,0.0005291796136841228,0.16900982400913736,0.5500547600738862 +59.34,31.622776601683793,50.118723362727195,0.00027617816908460685,0.16487619115421753,0.5862800757784927 +59.34,50.118723362727195,79.43282347242814,0.00016181676170701128,0.17042288768853997,0.5992219866891113 +66.44,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +66.44,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +66.44,0.03162277660168379,0.05011872336272725,0.0,0.32,0.1 +66.44,0.05011872336272725,0.07943282347242814,0.0,0.32,0.1 +66.44,0.07943282347242814,0.12589254117941673,0.0001791018992365193,0.32,0.1 +66.44,0.12589254117941673,0.19952623149688797,0.6138624957418299,0.32,0.1 +66.44,0.19952623149688797,0.31622776601683794,1.3995973259018928,0.30838869551613846,0.1 +66.44,0.31622776601683794,0.5011872336272722,6.545853525448173,0.27302051299491586,0.1 +66.44,0.5011872336272722,0.7943282347242814,2.8791051438911315,0.2385185727611427,0.20151533948234723 +66.44,0.7943282347242814,1.2589254117941668,0.28301727108012686,0.17951202616288997,0.42725626915203463 +66.44,1.2589254117941668,1.9952623149688788,0.029780751162947663,0.1494066867698645,0.5319012554914053 +66.44,1.9952623149688788,3.1622776601683795,0.01203016089856015,0.13832364567991648,0.49406518317601666 +66.44,3.1622776601683795,5.01187233627272,0.006788174228680906,0.13367638396697065,0.44175248770451625 +66.44,5.01187233627272,7.943282347242813,0.00402018645901959,0.1313429094712646,0.40625462802117973 +66.44,7.943282347242813,12.589254117941662,0.002493201006176852,0.13367578147823206,0.3873548159062976 +66.44,12.589254117941662,19.952623149688787,0.0012408610863292337,0.13483173498613846,0.39892274658665067 +66.44,19.952623149688787,31.622776601683793,0.0006402176564767995,0.13846844758081744,0.42660129373129496 +66.44,31.622776601683793,50.118723362727195,0.000291777326549409,0.13365007793852096,0.4744804124172203 +66.44,50.118723362727195,79.43282347242814,0.0001518523513920558,0.13785123669460694,0.5364545721486894 diff --git a/lstchain/data/LST1_backg_irf_gheffi_0.90_theffi_0.90.csv b/lstchain/data/LST1_backg_irf_gheffi_0.90_theffi_0.90.csv new file mode 100644 index 0000000000..44f63684e9 --- /dev/null +++ b/lstchain/data/LST1_backg_irf_gheffi_0.90_theffi_0.90.csv @@ -0,0 +1,191 @@ +ZD_deg,Ereco_min_TeV,Ereco_max_TeV,BckgRate_per_second,Theta_cut_deg,Gammaness_cut +6.0,0.012589254117941675,0.0199526231496888,0.17559648096477617,0.32,0.1 +6.0,0.0199526231496888,0.03162277660168379,15.736548406715594,0.32,0.11620018197905382 +6.0,0.03162277660168379,0.05011872336272725,24.719070277163965,0.32,0.11228635214620931 +6.0,0.05011872336272725,0.07943282347242814,16.14113612780422,0.32,0.10326804603068676 +6.0,0.07943282347242814,0.12589254117941673,5.524885886729925,0.32,0.13950141062486865 +6.0,0.12589254117941673,0.19952623149688797,1.3809775006953484,0.32,0.22828367421218862 +6.0,0.19952623149688797,0.31622776601683794,0.5388986696995435,0.32,0.2546930869373168 +6.0,0.31622776601683794,0.5011872336272722,0.24152775907106613,0.32,0.2816099714731492 +6.0,0.5011872336272722,0.7943282347242814,0.11803242158935899,0.32,0.30931832769270645 +6.0,0.7943282347242814,1.2589254117941668,0.05439237983742474,0.32,0.32972385843341534 +6.0,1.2589254117941668,1.9952623149688788,0.02740361620255927,0.32,0.33444928025707926 +6.0,1.9952623149688788,3.1622776601683795,0.011288435510305356,0.32,0.410147075635551 +6.0,3.1622776601683795,5.01187233627272,0.005975646622498844,0.32,0.45329619750767874 +6.0,5.01187233627272,7.943282347242813,0.003163268508234222,0.32,0.47244032671019254 +6.0,7.943282347242813,12.589254117941662,0.0016745079298216643,0.32,0.4666984680217185 +6.0,12.589254117941662,19.952623149688787,0.0008225144648288427,0.3082496059580866,0.49147566817020794 +6.0,19.952623149688787,31.622776601683793,0.00046923404012519956,0.32,0.5412176576412443 +6.0,31.622776601683793,50.118723362727195,6.667027008255354e-05,0.16578504144578401,0.3576959983091037 +6.0,50.118723362727195,79.43282347242814,0.00013148972333362836,0.32,0.1 +9.579,0.012589254117941675,0.0199526231496888,0.1481991677785433,0.32,0.1 +9.579,0.0199526231496888,0.03162277660168379,14.718800243714139,0.32,0.11015613325644188 +9.579,0.03162277660168379,0.05011872336272725,24.348141102914983,0.32,0.1093508969147617 +9.579,0.05011872336272725,0.07943282347242814,15.931940576847682,0.32,0.1 +9.579,0.07943282347242814,0.12589254117941673,5.534426659267126,0.32,0.1386379371118968 +9.579,0.12589254117941673,0.19952623149688797,1.368027052982405,0.32,0.21868047618084444 +9.579,0.19952623149688797,0.31622776601683794,0.5323497848356938,0.32,0.2501446318160868 +9.579,0.31622776601683794,0.5011872336272722,0.24433275473965724,0.32,0.27996738802813265 +9.579,0.5011872336272722,0.7943282347242814,0.11648402362321165,0.32,0.3035384465813813 +9.579,0.7943282347242814,1.2589254117941668,0.056219489786327466,0.32,0.3184095389054086 +9.579,1.2589254117941668,1.9952623149688788,0.02752388933883518,0.32,0.3652345161846954 +9.579,1.9952623149688788,3.1622776601683795,0.011307859534533645,0.32,0.4509771334446142 +9.579,3.1622776601683795,5.01187233627272,0.005469451122185254,0.32,0.48171932833846715 +9.579,5.01187233627272,7.943282347242813,0.0026454958594608396,0.32,0.48881721036185327 +9.579,7.943282347242813,12.589254117941662,0.001279588789821422,0.32,0.4749441813144524 +9.579,12.589254117941662,19.952623149688787,0.0006189189316555376,0.32,0.4745695200280388 +9.579,19.952623149688787,31.622776601683793,0.00019206303048174906,0.25631446180909623,0.5137827335639006 +9.579,31.622776601683793,50.118723362727195,6.59891033377573e-05,0.21602576938344725,0.2526272049317645 +9.579,50.118723362727195,79.43282347242814,7.00363958657308e-05,0.32,0.1 +16.08,0.012589254117941675,0.0199526231496888,0.05712969752402746,0.32,0.1 +16.08,0.0199526231496888,0.03162277660168379,10.510096452332292,0.32,0.10190888871381985 +16.08,0.03162277660168379,0.05011872336272725,24.75223918502653,0.32,0.10832150005806346 +16.08,0.05011872336272725,0.07943282347242814,16.911257020032323,0.32,0.1 +16.08,0.07943282347242814,0.12589254117941673,6.737544819392419,0.32,0.12289679410094916 +16.08,0.12589254117941673,0.19952623149688797,1.5764717901556917,0.32,0.21753616961917074 +16.08,0.19952623149688797,0.31622776601683794,0.5705973026151518,0.32,0.26065738337231786 +16.08,0.31622776601683794,0.5011872336272722,0.26773470580047337,0.32,0.27578920137471125 +16.08,0.5011872336272722,0.7943282347242814,0.12387952941388833,0.32,0.314850128704439 +16.08,0.7943282347242814,1.2589254117941668,0.06099758711348397,0.32,0.325454324523201 +16.08,1.2589254117941668,1.9952623149688788,0.028423235938300165,0.32,0.36894744146794856 +16.08,1.9952623149688788,3.1622776601683795,0.012697792669658139,0.32,0.4224723676467453 +16.08,3.1622776601683795,5.01187233627272,0.005633864740117978,0.32,0.46335199119040127 +16.08,5.01187233627272,7.943282347242813,0.0024996810654964925,0.32,0.5108895048274227 +16.08,7.943282347242813,12.589254117941662,0.0011090797733761056,0.32,0.5988211248295849 +16.08,12.589254117941662,19.952623149688787,0.0004920859547606634,0.32,0.48258021345782603 +16.08,19.952623149688787,31.622776601683793,0.00018343879420620128,0.2933161018614291,0.4862742934905485 +16.08,31.622776601683793,50.118723362727195,4.748483167002689e-05,0.22404150054737482,0.5646455548300677 +16.08,50.118723362727195,79.43282347242814,4.298090175087213e-05,0.32,0.1 +23.16,0.012589254117941675,0.0199526231496888,0.0041728407715988505,0.32,0.1 +23.16,0.0199526231496888,0.03162277660168379,4.085329145278117,0.32,0.1 +23.16,0.03162277660168379,0.05011872336272725,22.1583338576892,0.32,0.1002160464454597 +23.16,0.05011872336272725,0.07943282347242814,17.634969865505067,0.32,0.1 +23.16,0.07943282347242814,0.12589254117941673,8.560641907630382,0.32,0.11432499819699211 +23.16,0.12589254117941673,0.19952623149688797,1.9544998195484937,0.32,0.2146976687239026 +23.16,0.19952623149688797,0.31622776601683794,0.632303218717951,0.32,0.25873499020817864 +23.16,0.31622776601683794,0.5011872336272722,0.28143427610085064,0.32,0.2861409926057853 +23.16,0.5011872336272722,0.7943282347242814,0.14582926578850922,0.32,0.28607060600462253 +23.16,0.7943282347242814,1.2589254117941668,0.0702826964106501,0.32,0.3062777071535567 +23.16,1.2589254117941668,1.9952623149688788,0.03225750354416201,0.32,0.35014735631005917 +23.16,1.9952623149688788,3.1622776601683795,0.014861057358120013,0.32,0.40461859697891644 +23.16,3.1622776601683795,5.01187233627272,0.007831621051631366,0.32,0.4365758774608059 +23.16,5.01187233627272,7.943282347242813,0.004127181991047416,0.32,0.4968352080195375 +23.16,7.943282347242813,12.589254117941662,0.002174981536380379,0.32,0.48956967277584734 +23.16,12.589254117941662,19.952623149688787,0.001146192412609121,0.32,0.5883215945363233 +23.16,19.952623149688787,31.622776601683793,0.0006040313560128342,0.32,0.4256245535221811 +23.16,31.622776601683793,50.118723362727195,0.00015996630324149972,0.2268471004255254,0.5044199848424348 +23.16,50.118723362727195,79.43282347242814,0.00016775033394255662,0.32,0.1 +30.39,0.012589254117941675,0.0199526231496888,0.00044520408537425323,0.32,0.1 +30.39,0.0199526231496888,0.03162277660168379,1.050185033042877,0.32,0.1 +30.39,0.03162277660168379,0.05011872336272725,15.915030145752871,0.32,0.1 +30.39,0.05011872336272725,0.07943282347242814,19.306976247084496,0.32,0.1 +30.39,0.07943282347242814,0.12589254117941673,11.213696509608464,0.32,0.11019316618346073 +30.39,0.12589254117941673,0.19952623149688797,3.0889973115818297,0.32,0.19753178466885557 +30.39,0.19952623149688797,0.31622776601683794,0.7832783703342685,0.32,0.26763304986805486 +30.39,0.31622776601683794,0.5011872336272722,0.3412170249508419,0.32,0.2918270740287575 +30.39,0.5011872336272722,0.7943282347242814,0.17274182938029797,0.32,0.29551950640407904 +30.39,0.7943282347242814,1.2589254117941668,0.08500626806157821,0.32,0.3095920117078006 +30.39,1.2589254117941668,1.9952623149688788,0.04187648478036508,0.32,0.3525498890231738 +30.39,1.9952623149688788,3.1622776601683795,0.020873142253946236,0.32,0.37384226603178716 +30.39,3.1622776601683795,5.01187233627272,0.010118520466922519,0.32,0.405777918831999 +30.39,5.01187233627272,7.943282347242813,0.004905081141780332,0.32,0.4721346454914295 +30.39,7.943282347242813,12.589254117941662,0.0023778003005578423,0.32,0.5051540074317398 +30.39,12.589254117941662,19.952623149688787,0.0011526688562140343,0.32,0.5078703223644796 +30.39,19.952623149688787,31.622776601683793,0.0005587708487437167,0.32,0.5751575238410231 +30.39,31.622776601683793,50.118723362727195,0.00016794813524384343,0.2519740959959312,0.36249544819835267 +30.39,50.118723362727195,79.43282347242814,6.199561135441974e-05,0.2198791765944912,0.46335389409984495 +37.66,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +37.66,0.0199526231496888,0.03162277660168379,0.005439726626267147,0.32,0.1 +37.66,0.03162277660168379,0.05011872336272725,4.033267102327914,0.32,0.1 +37.66,0.05011872336272725,0.07943282347242814,17.523482783229554,0.32,0.1 +37.66,0.07943282347242814,0.12589254117941673,14.281379870888037,0.32,0.1 +37.66,0.12589254117941673,0.19952623149688797,5.678990057752521,0.32,0.15236231036817205 +37.66,0.19952623149688797,0.31622776601683794,1.2305712397647384,0.32,0.27694803767970055 +37.66,0.31622776601683794,0.5011872336272722,0.48448358873612357,0.32,0.29125177403665164 +37.66,0.5011872336272722,0.7943282347242814,0.2327857061679189,0.32,0.2933286427575119 +37.66,0.7943282347242814,1.2589254117941668,0.12593306045573732,0.32,0.2808656999339381 +37.66,1.2589254117941668,1.9952623149688788,0.06733221915550197,0.32,0.30008849367851587 +37.66,1.9952623149688788,3.1622776601683795,0.03382440263576279,0.32,0.32059189443595226 +37.66,3.1622776601683795,5.01187233627272,0.01641770122605414,0.32,0.3776824343131596 +37.66,5.01187233627272,7.943282347242813,0.008053027406759726,0.32,0.41732136463774966 +37.66,7.943282347242813,12.589254117941662,0.003950081044909462,0.32,0.4681958676935546 +37.66,12.589254117941662,19.952623149688787,0.0019375496286347814,0.32,0.4975242950021638 +37.66,19.952623149688787,31.622776601683793,0.0009503851998836204,0.32,0.4560346038262774 +37.66,31.622776601683793,50.118723362727195,0.00046617233169622414,0.32,0.3750940508406499 +37.66,50.118723362727195,79.43282347242814,8.019044891039412e-05,0.18950220722004738,0.41521133819475003 +44.92,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +44.92,0.0199526231496888,0.03162277660168379,0.0005034818183207038,0.32,0.1 +44.92,0.03162277660168379,0.05011872336272725,0.5484809263354772,0.32,0.1 +44.92,0.05011872336272725,0.07943282347242814,7.981817248040001,0.32,0.1 +44.92,0.07943282347242814,0.12589254117941673,14.655471398718488,0.32,0.1 +44.92,0.12589254117941673,0.19952623149688797,9.816223606160673,0.32,0.11136677915685111 +44.92,0.19952623149688797,0.31622776601683794,2.8933208884379975,0.32,0.22655455341978178 +44.92,0.31622776601683794,0.5011872336272722,0.7123626522119845,0.32,0.29012215229289173 +44.92,0.5011872336272722,0.7943282347242814,0.2842214873553521,0.32,0.28340725435932895 +44.92,0.7943282347242814,1.2589254117941668,0.1567624148375671,0.32,0.2737269917453996 +44.92,1.2589254117941668,1.9952623149688788,0.08875308840941676,0.32,0.2701209850027934 +44.92,1.9952623149688788,3.1622776601683795,0.04806543382272748,0.32,0.28067980473119364 +44.92,3.1622776601683795,5.01187233627272,0.026913273538198192,0.32,0.30028752922703805 +44.92,5.01187233627272,7.943282347242813,0.01209136211856364,0.32,0.3605503283144242 +44.92,7.943282347242813,12.589254117941662,0.005657276374754573,0.32,0.3926582316629572 +44.92,12.589254117941662,19.952623149688787,0.002646912371536695,0.32,0.4559466839610515 +44.92,19.952623149688787,31.622776601683793,0.0012384307639377004,0.32,0.47256881354857116 +44.92,31.622776601683793,50.118723362727195,0.0005794338994973615,0.32,0.4840381992257248 +44.92,50.118723362727195,79.43282347242814,0.00019247100719071677,0.26962765847382447,0.5393843603168154 +52.16,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +52.16,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +52.16,0.03162277660168379,0.05011872336272725,0.0019554053555110216,0.32,0.1 +52.16,0.05011872336272725,0.07943282347242814,0.3994934321970759,0.32,0.1 +52.16,0.07943282347242814,0.12589254117941673,8.500001960462598,0.32,0.1 +52.16,0.12589254117941673,0.19952623149688797,13.77121151419063,0.32,0.1 +52.16,0.19952623149688797,0.31622776601683794,7.350167883707211,0.32,0.12782242504578498 +52.16,0.31622776601683794,0.5011872336272722,1.7606005784971623,0.32,0.27287389189997224 +52.16,0.5011872336272722,0.7943282347242814,0.4662033361800171,0.32,0.29900323749359553 +52.16,0.7943282347242814,1.2589254117941668,0.22664777573667008,0.32,0.27923370595753566 +52.16,1.2589254117941668,1.9952623149688788,0.13279403949330562,0.32,0.2565722269327261 +52.16,1.9952623149688788,3.1622776601683795,0.0733972452119725,0.32,0.24227269207419672 +52.16,3.1622776601683795,5.01187233627272,0.043727538658971755,0.32,0.24930940071416305 +52.16,5.01187233627272,7.943282347242813,0.02150721255993863,0.32,0.26471029380820676 +52.16,7.943282347242813,12.589254117941662,0.011762951290919104,0.32,0.2937675351898848 +52.16,12.589254117941662,19.952623149688787,0.006433517253196768,0.32,0.340632595436056 +52.16,19.952623149688787,31.622776601683793,0.003518687038951978,0.32,0.3737626001490424 +52.16,31.622776601683793,50.118723362727195,0.0019244773878451216,0.32,0.42187429178046265 +52.16,50.118723362727195,79.43282347242814,0.00091399013595741,0.2981934958226351,0.46791605631265054 +59.34,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +59.34,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +59.34,0.03162277660168379,0.05011872336272725,0.0,0.32,0.1 +59.34,0.05011872336272725,0.07943282347242814,0.002694531473847235,0.32,0.1 +59.34,0.07943282347242814,0.12589254117941673,0.7654213051293,0.32,0.1 +59.34,0.12589254117941673,0.19952623149688797,5.273446261948501,0.32,0.1 +59.34,0.19952623149688797,0.31622776601683794,10.281973385352533,0.32,0.1 +59.34,0.31622776601683794,0.5011872336272722,6.131847616690863,0.32,0.12202484630493435 +59.34,0.5011872336272722,0.7943282347242814,1.8291917471893973,0.32,0.2924205960602483 +59.34,0.7943282347242814,1.2589254117941668,0.5972286277825254,0.32,0.3039130771250425 +59.34,1.2589254117941668,1.9952623149688788,0.21573662540572833,0.32,0.2655940062321337 +59.34,1.9952623149688788,3.1622776601683795,0.11928754625232638,0.32,0.22674592937409366 +59.34,3.1622776601683795,5.01187233627272,0.07339439344642261,0.32,0.20170444402723714 +59.34,5.01187233627272,7.943282347242813,0.04053636432479874,0.32,0.20614783961783145 +59.34,7.943282347242813,12.589254117941662,0.021372538870198188,0.32,0.21249691515716132 +59.34,12.589254117941662,19.952623149688787,0.010162926277315474,0.32,0.23612639342909444 +59.34,19.952623149688787,31.622776601683793,0.0051678215728923425,0.32,0.27812830854616394 +59.34,31.622776601683793,50.118723362727195,0.0026278238255907083,0.32,0.3079741076018832 +59.34,50.118723362727195,79.43282347242814,0.001336241578959414,0.32,0.34640349847763846 +66.44,0.012589254117941675,0.0199526231496888,0.0,0.32,0.1 +66.44,0.0199526231496888,0.03162277660168379,0.0,0.32,0.1 +66.44,0.03162277660168379,0.05011872336272725,0.0,0.32,0.1 +66.44,0.05011872336272725,0.07943282347242814,0.0,0.32,0.1 +66.44,0.07943282347242814,0.12589254117941673,0.0001791018992365193,0.32,0.1 +66.44,0.12589254117941673,0.19952623149688797,0.6138624957418299,0.32,0.1 +66.44,0.19952623149688797,0.31622776601683794,1.8041276536438466,0.32,0.1 +66.44,0.31622776601683794,0.5011872336272722,8.667526744433859,0.32,0.1 +66.44,0.5011872336272722,0.7943282347242814,8.771824605436732,0.32,0.1 +66.44,0.7943282347242814,1.2589254117941668,2.1405234715862367,0.32,0.23471231536544535 +66.44,1.2589254117941668,1.9952623149688788,0.2995577349270367,0.29488440944505034,0.3347833782701662 +66.44,1.9952623149688788,3.1622776601683795,0.135993795957133,0.2864069963983325,0.283200878397462 +66.44,3.1622776601683795,5.01187233627272,0.08799609209773075,0.2845632469458279,0.2299551460542766 +66.44,5.01187233627272,7.943282347242813,0.05817807877156834,0.2908403906472992,0.20391703600905553 +66.44,7.943282347242813,12.589254117941662,0.038460607554884396,0.31068923577798857,0.1862295796307054 +66.44,12.589254117941662,19.952623149688787,0.021556562587680985,0.31760805725788366,0.18385702585616567 +66.44,19.952623149688787,31.622776601683793,0.010124930615833161,0.31224688339635853,0.19600365005165704 +66.44,31.622776601683793,50.118723362727195,0.004644484386197977,0.3033689680537544,0.2193029044907234 +66.44,50.118723362727195,79.43282347242814,0.0020234883061995145,0.28724537451748516,0.27689025508229353 diff --git a/lstchain/data/LST1_gamma_irf_gheffi_0.40_theffi_0.40.csv b/lstchain/data/LST1_gamma_irf_gheffi_0.40_theffi_0.40.csv new file mode 100644 index 0000000000..d8a0e0ee4b --- /dev/null +++ b/lstchain/data/LST1_gamma_irf_gheffi_0.40_theffi_0.40.csv @@ -0,0 +1,351 @@ +ZD_deg,Etrue_min_TeV,Etrue_max_TeV,Aeff_m2,emig_mu_loc,emig_mu_scale,emig_mu_a,emig_model +6.0,0.011220184543019636,0.014125375446227547,601.8636022004555,2.1416143873932616,0.23758555043599847,,moyal +6.0,0.014125375446227547,0.01778279410038923,1229.8921003382998,1.7960288037431202,0.2188624433856196,,moyal +6.0,0.01778279410038923,0.0223872113856834,2391.6360577355,1.4318091256862098,0.20194091802666045,,moyal +6.0,0.0223872113856834,0.02818382931264455,3886.8925701871162,1.2220059711278857,0.18583359986525078,,moyal +6.0,0.02818382931264455,0.03548133892335755,5221.962832997907,1.0412027798004448,0.17579638320259053,,moyal +6.0,0.03548133892335755,0.04466835921509633,6785.535023650937,0.9145783919288173,0.17295748420269527,,moyal +6.0,0.04466835921509633,0.05623413251903491,9491.647193215018,0.6715237139146365,0.441722687262099,1.9999999999999911,skewnorm +6.0,0.05623413251903491,0.0707945784384138,13231.637162196841,0.6686668205737102,0.42765325970874074,1.9999999999999998,skewnorm +6.0,0.0707945784384138,0.08912509381337459,16841.09849510246,0.6430424073651757,0.43304714105038133,1.7291063553365082,skewnorm +6.0,0.08912509381337459,0.11220184543019636,21247.877896919395,0.9235413936237876,0.29686285173808485,-0.0003604688722291855,skewnorm +6.0,0.11220184543019636,0.14125375446227548,27214.702012542348,0.762158862205991,0.32153813948252724,0.8805451416623801,skewnorm +6.0,0.14125375446227548,0.17782794100389238,30417.867124438242,0.7866083804609031,0.29378826035323297,0.8634935147733573,skewnorm +6.0,0.17782794100389238,0.223872113856834,36680.96920919961,1.1280300221302992,0.29276998577519864,-1.1388125041256871,skewnorm +6.0,0.223872113856834,0.2818382931264455,41943.63176485503,0.80969922456544,0.2552237448878649,1.2629563994863944,skewnorm +6.0,0.2818382931264455,0.35481338923357547,52146.593884762355,0.8186525141647494,0.2456029394426815,1.4930686871419794,skewnorm +6.0,0.35481338923357547,0.44668359215096326,62614.48985238869,0.837853782342887,0.2257322368400467,1.2048265698927283,skewnorm +6.0,0.44668359215096326,0.5623413251903494,67483.29650785784,0.8352802589426394,0.24054882015993603,1.4958118976601982,skewnorm +6.0,0.5623413251903494,0.707945784384138,76941.67865242239,0.8264223819488132,0.22616353173889342,1.395257366218998,skewnorm +6.0,0.707945784384138,0.8912509381337459,88424.93981088648,0.8114273128889187,0.24488329260201955,1.9999999999990314,skewnorm +6.0,0.8912509381337459,1.1220184543019636,82451.00444782853,0.8377273903261895,0.22621715029817144,1.5901410902279516,skewnorm +6.0,1.1220184543019636,1.4125375446227548,91587.93946528672,0.8772816746913944,0.17731376751164787,1.0124588057373753,skewnorm +6.0,1.4125375446227548,1.7782794100389236,96634.37931229484,0.8575923119747371,0.19497171247089962,1.3613234131267542,skewnorm +6.0,1.7782794100389236,2.238721138568341,109674.25467072222,0.8818674605221172,0.2057950286935981,1.111704584170622,skewnorm +6.0,2.238721138568341,2.8183829312644564,93981.58154780796,1.0712835429815568,0.17026236532841052,-0.867978642086584,skewnorm +6.0,2.8183829312644564,3.548133892335755,89832.36761660763,1.129526145527036,0.18462666388786717,-1.2776976818785515,skewnorm +6.0,3.548133892335755,4.466835921509633,112172.8012039856,1.069634475163725,0.140474188989395,-0.5233687080926358,skewnorm +6.0,4.466835921509633,5.623413251903494,111121.72323763899,1.0987364224872207,0.1610931323622503,-1.9999999999999998,skewnorm +6.0,5.623413251903494,7.079457843841384,97634.32707788429,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,7.079457843841384,8.912509381337463,75216.23510884547,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,8.912509381337463,11.220184543019636,85453.42197674119,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,11.220184543019636,14.125375446227547,95949.26991601476,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,14.125375446227547,17.782794100389236,62226.68364806862,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,17.782794100389236,22.38721138568341,64514.26745336599,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,22.38721138568341,28.183829312644566,87019.98291463699,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +6.0,28.183829312644566,35.48133892335755,51124.111651727944,1.1017127009203993,0.18392012695263915,-1.418079505038764,skewnorm +9.579,0.011220184543019636,0.014125375446227547,577.6228719953638,2.1484868689560477,0.25723709884370793,,moyal +9.579,0.014125375446227547,0.01778279410038923,1260.5755037770095,1.7691650948065145,0.21476631343564587,,moyal +9.579,0.01778279410038923,0.0223872113856834,2002.8822859939603,1.455989576053346,0.19817625065276323,,moyal +9.579,0.0223872113856834,0.02818382931264455,3416.5753445368646,1.2251659615003505,0.17930714099302433,,moyal +9.579,0.02818382931264455,0.03548133892335755,5070.865599783437,1.0370176067044186,0.17560326084363387,,moyal +9.579,0.03548133892335755,0.04466835921509633,7195.150404019665,0.9046931452520488,0.16905598795059543,,moyal +9.579,0.04466835921509633,0.05623413251903491,9256.596706757788,0.8039770661842053,0.1658656824077297,,moyal +9.579,0.05623413251903491,0.0707945784384138,12799.597426198787,0.7436361986270849,0.17449273158318118,,moyal +9.579,0.0707945784384138,0.08912509381337459,17566.00995214436,0.6389437363807903,0.4308435768992372,1.9832041509643914,skewnorm +9.579,0.08912509381337459,0.11220184543019636,22382.13665760197,0.9091842052031807,0.3108944654564607,-0.0004828568071997037,skewnorm +9.579,0.11220184543019636,0.14125375446227548,27136.546600594393,0.7451959671961274,0.3312736931673839,0.8819114018855561,skewnorm +9.579,0.14125375446227548,0.17782794100389238,31043.748131082375,0.7612105687718028,0.3024940176762347,0.9549461641391844,skewnorm +9.579,0.17782794100389238,0.223872113856834,37292.645812300594,0.7909511031983759,0.288414944952662,1.097524962453271,skewnorm +9.579,0.223872113856834,0.2818382931264455,41069.10773741235,0.796510294496107,0.2717201804219316,1.187459896950261,skewnorm +9.579,0.2818382931264455,0.35481338923357547,52402.627076919634,0.7764751296514665,0.26837700075000526,1.7872807042674408,skewnorm +9.579,0.35481338923357547,0.44668359215096326,59526.88865181214,0.7868012234829366,0.2624607124831607,1.732439177464561,skewnorm +9.579,0.44668359215096326,0.5623413251903494,68986.74216511779,0.8348550918744085,0.21747059488346848,1.2364700154964643,skewnorm +9.579,0.5623413251903494,0.707945784384138,78031.98890139583,0.8089540099827887,0.24690792986187993,1.65285914474079,skewnorm +9.579,0.707945784384138,0.8912509381337459,82048.7142051187,0.8084770159795305,0.22521570741137542,1.8737945257405917,skewnorm +9.579,0.8912509381337459,1.1220184543019636,97936.40441749875,0.8630794350270525,0.20944655764187797,1.1729927604518828,skewnorm +9.579,1.1220184543019636,1.4125375446227548,97638.79003779501,0.8631757021614802,0.19256553927004416,1.2593411482529409,skewnorm +9.579,1.4125375446227548,1.7782794100389236,102617.85945827473,0.8774903778955321,0.17213836059009238,1.1308631070053208,skewnorm +9.579,1.7782794100389236,2.238721138568341,117105.8784386259,1.0878368001738883,0.18292238495479238,-1.127055815102826,skewnorm +9.579,2.238721138568341,2.8183829312644564,94482.34899703837,0.9072912469220058,0.16465334577867777,0.7130159243197379,skewnorm +9.579,2.8183829312644564,3.548133892335755,118209.72095355761,1.0608185060290445,0.14171946438339467,-0.8285934180315704,skewnorm +9.579,3.548133892335755,4.466835921509633,111265.2124332937,1.1437473099332087,0.20796228714413495,-1.9999999999999998,skewnorm +9.579,4.466835921509633,5.623413251903494,126650.78722360879,0.8987800290449995,0.17112836809195922,0.9976923689917807,skewnorm +9.579,5.623413251903494,7.079457843841384,89641.26785282741,0.9145144849616571,0.17263188905195076,1.5639729884506854,skewnorm +9.579,7.079457843841384,8.912509381337463,96201.42364523513,1.1098713684713082,0.20965331880829063,-1.359216213810098,skewnorm +9.579,8.912509381337463,11.220184543019636,107136.13800725206,1.0644614191669302,0.14591752986765005,-0.7964968614006599,skewnorm +9.579,11.220184543019636,14.125375446227547,112885.68818183517,1.1262250581952766,0.18062376285007736,-1.548317144434663,skewnorm +9.579,14.125375446227547,17.782794100389236,83054.02294284823,1.1262250581952766,0.18062376285007736,-1.548317144434663,skewnorm +9.579,17.782794100389236,22.38721138568341,83647.05602758989,1.1262250581952766,0.18062376285007736,-1.548317144434663,skewnorm +9.579,22.38721138568341,28.183829312644566,73128.75309368294,1.1262250581952766,0.18062376285007736,-1.548317144434663,skewnorm +9.579,28.183829312644566,35.48133892335755,81016.00813016412,1.1262250581952766,0.18062376285007736,-1.548317144434663,skewnorm +16.08,0.011220184543019636,0.014125375446227547,448.42003443324154,2.3052122669323105,0.25077700469375724,,moyal +16.08,0.014125375446227547,0.01778279410038923,1031.013583982261,1.8875594311278292,0.23594922430865098,,moyal +16.08,0.01778279410038923,0.0223872113856834,1791.6752771506478,1.5453666766991405,0.18145421395521016,,moyal +16.08,0.0223872113856834,0.02818382931264455,3187.894224338138,1.28344005338625,0.18028274836018995,,moyal +16.08,0.02818382931264455,0.03548133892335755,4739.000169650142,1.1000868419335221,0.1723110467213821,,moyal +16.08,0.03548133892335755,0.04466835921509633,6553.913962126864,0.9481639918666261,0.15423070751126713,,moyal +16.08,0.04466835921509633,0.05623413251903491,9005.970175378818,0.8384219442826707,0.1693551341700048,,moyal +16.08,0.05623413251903491,0.0707945784384138,12470.999652286637,0.7740204306933869,0.17018031549834833,,moyal +16.08,0.0707945784384138,0.08912509381337459,17166.24252163279,0.6220131523047413,0.44561299913913227,1.8630953616780925,skewnorm +16.08,0.08912509381337459,0.11220184543019636,22688.576666361372,0.7335495730912144,0.3473661482121987,0.9008429222054214,skewnorm +16.08,0.11220184543019636,0.14125375446227548,26427.892225421307,0.8982079972469815,0.27773783994007306,0.0005993180898318624,skewnorm +16.08,0.14125375446227548,0.17782794100389238,31879.19302619991,0.9145367333628674,0.24779963302362892,-0.00017113798408734785,skewnorm +16.08,0.17782794100389238,0.223872113856834,38122.86839714995,0.9410679367613531,0.237526215732425,-0.0003137074611441929,skewnorm +16.08,0.223872113856834,0.2818382931264455,43232.05610835921,0.7862549948311103,0.26868337310973833,1.2845108963335743,skewnorm +16.08,0.2818382931264455,0.35481338923357547,55246.26695777232,0.7752709161880613,0.2628085605908604,1.5312742886519275,skewnorm +16.08,0.35481338923357547,0.44668359215096326,62043.71884146762,0.9444000726799685,0.1784687377499227,0.00032646442506826634,skewnorm +16.08,0.44668359215096326,0.5623413251903494,70057.5470055841,0.7866181511877588,0.23899402242496853,1.6110344437283592,skewnorm +16.08,0.5623413251903494,0.707945784384138,81146.48436041152,0.7944562799414553,0.22539171294070515,1.6584921707449116,skewnorm +16.08,0.707945784384138,0.8912509381337459,91695.94920373558,0.8127190383277385,0.2189695232857249,1.3638892798813822,skewnorm +16.08,0.8912509381337459,1.1220184543019636,97379.29123293186,0.7982523869554485,0.22830975754398689,1.8078307679594792,skewnorm +16.08,1.1220184543019636,1.4125375446227548,110158.08689018288,0.8191719943805837,0.2173586580688261,1.71957974924571,skewnorm +16.08,1.4125375446227548,1.7782794100389236,107212.7622607319,1.089090580288389,0.18137527895990505,-1.062201391499608,skewnorm +16.08,1.7782794100389236,2.238721138568341,111880.16730651812,1.0339389891601434,0.15871027729220427,-0.7152039708407515,skewnorm +16.08,2.238721138568341,2.8183829312644564,111776.79434393164,0.8940627163001567,0.1558277388476233,0.8509667182684624,skewnorm +16.08,2.8183829312644564,3.548133892335755,123049.93782713456,0.8619523323319404,0.19199494526741517,1.1249651646165926,skewnorm +16.08,3.548133892335755,4.466835921509633,120353.6927590569,0.8890166143864605,0.17180277667038865,1.069278853872016,skewnorm +16.08,4.466835921509633,5.623413251903494,102010.9892158398,0.8734441696800795,0.18887307121684285,1.029756160583745,skewnorm +16.08,5.623413251903494,7.079457843841384,112056.43301903197,1.0011599425806093,0.12469502934666248,-0.00017696336358966075,skewnorm +16.08,7.079457843841384,8.912509381337463,94311.30471257173,1.0453950223275807,0.15889958706272453,-0.719596755251878,skewnorm +16.08,8.912509381337463,11.220184543019636,121724.114041165,1.116960551433944,0.20136075790811564,-1.6025532036022114,skewnorm +16.08,11.220184543019636,14.125375446227547,102998.43928158499,1.0752436443759719,0.1950215941310804,-1.9999999999999956,skewnorm +16.08,14.125375446227547,17.782794100389236,112272.95162996995,1.159792514008744,0.2574132998069626,-1.7726857826184697,skewnorm +16.08,17.782794100389236,22.38721138568341,95556.01478136779,1.159792514008744,0.2574132998069626,-1.7726857826184697,skewnorm +16.08,22.38721138568341,28.183829312644566,80198.59683869533,1.159792514008744,0.2574132998069626,-1.7726857826184697,skewnorm +16.08,28.183829312644566,35.48133892335755,82033.2918847564,1.159792514008744,0.2574132998069626,-1.7726857826184697,skewnorm +23.16,0.011220184543019636,0.014125375446227547,319.8570954587096,2.526474227719214,0.3375722296360266,,moyal +23.16,0.014125375446227547,0.01778279410038923,713.9889378563084,2.0935354606808327,0.2722177530066669,,moyal +23.16,0.01778279410038923,0.0223872113856834,1490.9980386849822,1.7057825291071593,0.23465590831117525,,moyal +23.16,0.0223872113856834,0.02818382931264455,2426.764377787693,1.4122882349989736,0.20185027571509415,,moyal +23.16,0.02818382931264455,0.03548133892335755,3848.4325875608006,1.2059644653198542,0.19229101538316581,,moyal +23.16,0.03548133892335755,0.04466835921509633,5932.864783999707,1.0301644179097256,0.16851563195932706,,moyal +23.16,0.04466835921509633,0.05623413251903491,8913.618227816893,0.8966016638205843,0.16369561539414665,,moyal +23.16,0.05623413251903491,0.0707945784384138,11625.719051380784,0.6862957286045551,0.465715280551969,1.9999999999999998,skewnorm +23.16,0.0707945784384138,0.08912509381337459,16153.774356680318,0.665996341447093,0.4269640278651275,1.9999999999999996,skewnorm +23.16,0.08912509381337459,0.11220184543019636,21916.33759098714,0.6432258344259812,0.41058644467601085,1.7901618488488857,skewnorm +23.16,0.11220184543019636,0.14125375446227548,27228.171235266633,0.7195282488918937,0.3601631311092426,1.0331112637724331,skewnorm +23.16,0.14125375446227548,0.17782794100389238,34548.94339105637,1.0500878592723308,0.3040410064494574,-0.6876072049579943,skewnorm +23.16,0.17782794100389238,0.223872113856834,38531.64931401624,1.1121215823818726,0.30295414061228393,-0.9181464375974147,skewnorm +23.16,0.223872113856834,0.2818382931264455,46765.49253990874,1.0449323792680734,0.23453584665759347,-0.5242852753487135,skewnorm +23.16,0.2818382931264455,0.35481338923357547,57962.79607756052,0.9496095211418712,0.1993589242074982,-4.5155971426998e-05,skewnorm +23.16,0.35481338923357547,0.44668359215096326,65406.009330770365,0.7839180207227946,0.26424231055443637,1.83699808492736,skewnorm +23.16,0.44668359215096326,0.5623413251903494,72550.21798871357,0.7931744840997517,0.25309732681972896,1.6112844731201978,skewnorm +23.16,0.5623413251903494,0.707945784384138,87159.21467257694,0.7927970413116902,0.2367972507502024,1.9025012129068453,skewnorm +23.16,0.707945784384138,0.8912509381337459,97007.44285586497,0.7869441519459668,0.2411748205510364,1.9999999999999998,skewnorm +23.16,0.8912509381337459,1.1220184543019636,102482.63072258966,0.8129873846433668,0.21994706487278654,1.69610014479073,skewnorm +23.16,1.1220184543019636,1.4125375446227548,117594.52038923147,0.8055008127953105,0.23335121433593883,1.9999999999999667,skewnorm +23.16,1.4125375446227548,1.7782794100389236,118264.99231316563,0.8339754176096734,0.21467928487528387,1.3492314746585512,skewnorm +23.16,1.7782794100389236,2.238721138568341,129716.82936559332,0.8449311554492258,0.19296465116652198,1.7278512537949964,skewnorm +23.16,2.238721138568341,2.8183829312644564,128176.06326878654,0.873961365048221,0.16777082468524412,1.0992228123233054,skewnorm +23.16,2.8183829312644564,3.548133892335755,129293.83166980627,1.0887161150280225,0.1736320459096489,-1.1544234895524454,skewnorm +23.16,3.548133892335755,4.466835921509633,140671.79394727087,1.0663595605286538,0.1618936823894502,-0.735818019242498,skewnorm +23.16,4.466835921509633,5.623413251903494,133222.7454217528,0.8629863185278506,0.182220832317389,1.8946711735316866,skewnorm +23.16,5.623413251903494,7.079457843841384,147510.57197290813,1.0787531985468883,0.15472835838953053,-0.8239873813005449,skewnorm +23.16,7.079457843841384,8.912509381337463,143306.44967343702,1.0801543325434826,0.16769323480721993,-1.2253800459335467,skewnorm +23.16,8.912509381337463,11.220184543019636,127035.76099987696,0.8374995166772008,0.1729119083177262,1.9999999999896927,skewnorm +23.16,11.220184543019636,14.125375446227547,130361.92546006244,0.9762206504865573,0.12732795420522935,-0.004280617554575218,skewnorm +23.16,14.125375446227547,17.782794100389236,111666.51634719495,1.136180060132352,0.2415803244508058,-1.999999999999998,skewnorm +23.16,17.782794100389236,22.38721138568341,110759.85430231286,1.0814513016631881,0.18803919203779465,-1.9999999999999998,skewnorm +23.16,22.38721138568341,28.183829312644566,117986.33438967797,1.0814513016631881,0.18803919203779465,-1.9999999999999998,skewnorm +23.16,28.183829312644566,35.48133892335755,124905.26819276328,1.0814513016631881,0.18803919203779465,-1.9999999999999998,skewnorm +30.39,0.011220184543019636,0.014125375446227547,190.33590524122098,3.1245423495517994,0.39404022139948836,,moyal +30.39,0.014125375446227547,0.01778279410038923,417.44894392547553,2.513158292721154,0.23211214648262693,,moyal +30.39,0.01778279410038923,0.0223872113856834,836.6853930786237,2.0279675594765036,0.23878211680292447,,moyal +30.39,0.0223872113856834,0.02818382931264455,1929.9006767417802,1.6710507610297642,0.20245693070179463,,moyal +30.39,0.02818382931264455,0.03548133892335755,3274.5488781705603,1.3710942260440329,0.18258791838573815,,moyal +30.39,0.03548133892335755,0.04466835921509633,4887.084406064533,1.150878797406265,0.16316299095099437,,moyal +30.39,0.04466835921509633,0.05623413251903491,7077.013943718682,0.986526545540296,0.16390805352764037,,moyal +30.39,0.05623413251903491,0.0707945784384138,10175.352086164288,0.8771409210120826,0.16460196418046122,,moyal +30.39,0.0707945784384138,0.08912509381337459,15023.059273574618,0.7195152994918327,0.41464338073216167,1.9999999999999998,skewnorm +30.39,0.08912509381337459,0.11220184543019636,20890.731849179374,0.6448006704652228,0.43565040213010775,1.9999999998088793,skewnorm +30.39,0.11220184543019636,0.14125375446227548,27999.728450265167,0.680696339133565,0.36995226521674285,1.3292744161463492,skewnorm +30.39,0.14125375446227548,0.17782794100389238,35144.077156469575,0.6895266374388364,0.36009329764755205,1.3724161275031859,skewnorm +30.39,0.17782794100389238,0.223872113856834,41663.5175328948,0.9274818646290719,0.26627914201665714,0.0008337319275586882,skewnorm +30.39,0.223872113856834,0.2818382931264455,48605.79055902879,0.8029996670982712,0.27151478455348904,0.933543932740076,skewnorm +30.39,0.2818382931264455,0.35481338923357547,57642.88521077106,0.774533292639352,0.27008994211224985,1.5554051957042336,skewnorm +30.39,0.35481338923357547,0.44668359215096326,72870.93324302262,0.7763736451204035,0.268617042364806,1.8111385679785286,skewnorm +30.39,0.44668359215096326,0.5623413251903494,83399.15421890686,0.7801375414001886,0.2545565118958336,1.6400057910752819,skewnorm +30.39,0.5623413251903494,0.707945784384138,94373.99075198393,0.7817087274120201,0.24342573230523903,1.9999999999999967,skewnorm +30.39,0.707945784384138,0.8912509381337459,108386.57645841916,0.7859419275628735,0.2356068348750221,1.9999999999999991,skewnorm +30.39,0.8912509381337459,1.1220184543019636,117623.28310166084,0.7978100888492264,0.23251937810805481,1.9999999999999998,skewnorm +30.39,1.1220184543019636,1.4125375446227548,130840.97133877648,0.8081080380445416,0.2170047401536409,1.9999999999992755,skewnorm +30.39,1.4125375446227548,1.7782794100389236,139099.41344520342,0.8172794891099574,0.20314538204661978,1.6307362638231135,skewnorm +30.39,1.7782794100389236,2.238721138568341,137550.3197238099,0.8140412523970304,0.21386621934769123,1.8439492611186932,skewnorm +30.39,2.238721138568341,2.8183829312644564,175315.23474362228,0.8461252847537504,0.17146086542286745,1.313173852828904,skewnorm +30.39,2.8183829312644564,3.548133892335755,152752.14234257577,0.8310271769407577,0.1908447013737304,1.999999999999999,skewnorm +30.39,3.548133892335755,4.466835921509633,144227.6652758006,0.8219458153943283,0.20286185007362173,1.8393741519210118,skewnorm +30.39,4.466835921509633,5.623413251903494,139157.44589669027,0.8539734997057823,0.17091617664347403,1.1281029417828707,skewnorm +30.39,5.623413251903494,7.079457843841384,171589.07409451588,1.025342414032965,0.15542114369914156,-0.6565513858406671,skewnorm +30.39,7.079457843841384,8.912509381337463,151816.73775708902,0.8398000902861535,0.17741678389327706,1.4696061016939104,skewnorm +30.39,8.912509381337463,11.220184543019636,162599.68803867447,0.8741450448172374,0.1474086945830828,1.0186814240478028,skewnorm +30.39,11.220184543019636,14.125375446227547,131678.92815430395,0.9832318782531893,0.1343707926988832,0.0015890668779381819,skewnorm +30.39,14.125375446227547,17.782794100389236,155224.69756073612,0.839937659284341,0.19332684101442743,1.191900760872225,skewnorm +30.39,17.782794100389236,22.38721138568341,134704.45103986902,1.0771388606571897,0.15732880884806244,-1.9999999999999998,skewnorm +30.39,22.38721138568341,28.183829312644566,133755.4924474502,0.9695428801350161,0.17601811889892854,0.0023519621456727625,skewnorm +30.39,28.183829312644566,35.48133892335755,138318.86905036395,0.9695428801350161,0.17601811889892854,0.0023519621456727625,skewnorm +37.66,0.011220184543019636,0.014125375446227547,64.92944788638971,3.921019318614137,0.6715515033519495,,moyal +37.66,0.014125375446227547,0.01778279410038923,158.3738305910112,3.0627383839823694,0.3848132034873237,,moyal +37.66,0.01778279410038923,0.0223872113856834,430.9198797713534,2.4334561252126004,0.2965817414264588,,moyal +37.66,0.0223872113856834,0.02818382931264455,999.9739463551232,2.0049441312290237,0.22684587324859545,,moyal +37.66,0.02818382931264455,0.03548133892335755,2010.66046766127,1.6494550138389177,0.212888171934335,,moyal +37.66,0.03548133892335755,0.04466835921509633,3512.98143873947,1.3911735610333238,0.1905086708658879,,moyal +37.66,0.04466835921509633,0.05623413251903491,5476.7311922861445,1.167940230268506,0.17927545035556897,,moyal +37.66,0.05623413251903491,0.0707945784384138,8542.225221370521,0.9965317264247259,0.1621505240826098,,moyal +37.66,0.0707945784384138,0.08912509381337459,12392.732885130868,0.7472911541713698,0.44651065527007644,1.9999999999999998,skewnorm +37.66,0.08912509381337459,0.11220184543019636,19066.74139756289,0.7304766976340002,0.3912175370310156,1.999999999999997,skewnorm +37.66,0.11220184543019636,0.14125375446227548,25606.55101019495,0.6486633913706725,0.4400584142454941,1.9999999999999998,skewnorm +37.66,0.14125375446227548,0.17782794100389238,34574.02898466665,0.7094133102652301,0.37850254757528046,1.2552992615053125,skewnorm +37.66,0.17782794100389238,0.223872113856834,42947.276077212184,0.9335677627575242,0.272675050333503,-0.00017515355032631563,skewnorm +37.66,0.223872113856834,0.2818382931264455,53419.66081844478,1.0767135090829711,0.2796413519653677,-0.6860243892597435,skewnorm +37.66,0.2818382931264455,0.35481338923357547,64550.10094532319,0.9475489843445613,0.21825389489574565,0.0038103418077241958,skewnorm +37.66,0.35481338923357547,0.44668359215096326,76836.44334046598,0.7854192794654812,0.27113583052862567,1.4569083157916947,skewnorm +37.66,0.44668359215096326,0.5623413251903494,92934.79700184656,0.798243628161538,0.2574585514815373,1.598361517192944,skewnorm +37.66,0.5623413251903494,0.707945784384138,110438.76893343353,0.7922596642801077,0.23965581311640283,1.7107528632763862,skewnorm +37.66,0.707945784384138,0.8912509381337459,128536.19645321678,0.7960030298379623,0.23647393555689536,1.9658882484815814,skewnorm +37.66,0.8912509381337459,1.1220184543019636,135626.9052075163,0.809430416391536,0.22856213563804767,1.790735899688733,skewnorm +37.66,1.1220184543019636,1.4125375446227548,154207.43873017558,0.8161293287947017,0.21771975945773728,1.8959124958252234,skewnorm +37.66,1.4125375446227548,1.7782794100389236,173955.77198786853,0.8340391284763842,0.20472950370551599,1.4959872851432297,skewnorm +37.66,1.7782794100389236,2.238721138568341,185713.60490145994,0.8256002495738853,0.20535471947760603,1.8059417360526213,skewnorm +37.66,2.238721138568341,2.8183829312644564,194731.76850073464,0.9669392333052053,0.13977657673421187,0.00034095227230425797,skewnorm +37.66,2.8183829312644564,3.548133892335755,181953.30467809268,0.8344839598923958,0.19842080001777887,1.9999999999997844,skewnorm +37.66,3.548133892335755,4.466835921509633,212382.98743794774,0.9761438380086012,0.1322035583670936,0.00321016166358821,skewnorm +37.66,4.466835921509633,5.623413251903494,203963.80915224907,0.881160798744469,0.15778190189679964,1.1683403249801991,skewnorm +37.66,5.623413251903494,7.079457843841384,207454.04109042717,0.908387872132375,0.1431609462311097,0.6945817928101907,skewnorm +37.66,7.079457843841384,8.912509381337463,227242.77759408377,0.8312351012494936,0.20364869818364326,1.9999999999999996,skewnorm +37.66,8.912509381337463,11.220184543019636,190989.75956138683,0.8481612410482346,0.19078793248955123,1.9999999999999998,skewnorm +37.66,11.220184543019636,14.125375446227547,190730.25316626992,0.8480452231810757,0.19927653852509647,1.3069263277555823,skewnorm +37.66,14.125375446227547,17.782794100389236,215848.20459581327,0.8796544543818551,0.191118641305355,1.3590634880430144,skewnorm +37.66,17.782794100389236,22.38721138568341,189733.37989932703,0.8523295655393635,0.22078863478105082,1.9886727552437686,skewnorm +37.66,22.38721138568341,28.183829312644566,155866.38435089984,0.8523632369518201,0.20217209926533652,1.8172292304056852,skewnorm +37.66,28.183829312644566,35.48133892335755,157998.6679277556,1.1079828256344317,0.1815993371847924,-0.804142546929144,skewnorm +44.92,0.011220184543019636,0.014125375446227547,11.213086235985653,,,,-- +44.92,0.014125375446227547,0.01778279410038923,34.33302891688256,,,,-- +44.92,0.01778279410038923,0.0223872113856834,109.40751649854745,3.316969324147977,0.38057560700315723,,moyal +44.92,0.0223872113856834,0.02818382931264455,386.0006163804636,2.7228143591302794,0.2710734470693519,,moyal +44.92,0.02818382931264455,0.03548133892335755,839.166634768259,2.1770419974526236,0.23136026539191032,,moyal +44.92,0.03548133892335755,0.04466835921509633,1972.7553145595766,1.7922539028230693,0.2048196543781041,,moyal +44.92,0.04466835921509633,0.05623413251903491,3179.0820130542825,1.470564100289258,0.17324027441557185,,moyal +44.92,0.05623413251903491,0.0707945784384138,6056.732202518578,1.216153541423012,0.1583433729010365,,moyal +44.92,0.0707945784384138,0.08912509381337459,9474.756229059174,1.0683421424748052,0.17278988058228748,,moyal +44.92,0.08912509381337459,0.11220184543019636,14459.83673286641,0.9121970528223324,0.15883390392409874,,moyal +44.92,0.11220184543019636,0.14125375446227548,21075.90621434399,0.8279461787014329,0.1619935161536369,,moyal +44.92,0.14125375446227548,0.17782794100389238,30191.60730378354,0.9904721595241381,0.27020857014952915,0.0008648734089735828,skewnorm +44.92,0.17782794100389238,0.223872113856834,41372.24975741374,0.9582612552548923,0.2519164055095203,-0.00012924623290141283,skewnorm +44.92,0.223872113856834,0.2818382931264455,55213.39213160205,0.7955982716128046,0.31333640783204547,0.7634347386764102,skewnorm +44.92,0.2818382931264455,0.35481338923357547,68653.24994137265,0.8100796514210622,0.2809055241022709,0.8200383321376793,skewnorm +44.92,0.35481338923357547,0.44668359215096326,83869.05107799724,0.7948998127318605,0.2693199118766278,1.2030184186796278,skewnorm +44.92,0.44668359215096326,0.5623413251903494,100359.92096707094,0.7847283302373297,0.26892347995242494,1.565413251114323,skewnorm +44.92,0.5623413251903494,0.707945784384138,120622.08561292324,0.8020964887517017,0.24903379459105868,1.4514453310607875,skewnorm +44.92,0.707945784384138,0.8912509381337459,138948.75196304705,0.7981304053865954,0.2351839091282725,1.5607627246917297,skewnorm +44.92,0.8912509381337459,1.1220184543019636,169036.57454104602,0.8062044931573847,0.23075258977878668,1.6433820648134825,skewnorm +44.92,1.1220184543019636,1.4125375446227548,189623.499097919,0.8113897054837546,0.2217601255444703,1.7900619718529436,skewnorm +44.92,1.4125375446227548,1.7782794100389236,209324.18567761805,0.8426123181458894,0.18901294303838212,1.2362590564868239,skewnorm +44.92,1.7782794100389236,2.238721138568341,222597.0212217358,0.8308077588074135,0.19226354547821392,1.6258384256918035,skewnorm +44.92,2.238721138568341,2.8183829312644564,237381.87685776548,0.9598120740939853,0.1345358869382871,-0.00013543041339252832,skewnorm +44.92,2.8183829312644564,3.548133892335755,258172.18406390908,0.8272589117182718,0.20747744094056925,1.9999999999999738,skewnorm +44.92,3.548133892335755,4.466835921509633,268154.58169218263,0.8360459051785001,0.18563911516266096,1.9558248665175046,skewnorm +44.92,4.466835921509633,5.623413251903494,279229.93561832106,0.8503904287153613,0.17399116302334247,1.6334656258507876,skewnorm +44.92,5.623413251903494,7.079457843841384,327183.1359047416,0.8390045408666583,0.18000149774512186,1.9999999999999971,skewnorm +44.92,7.079457843841384,8.912509381337463,311881.87360126706,0.8874008865651306,0.07557404084609906,,moyal +44.92,8.912509381337463,11.220184543019636,299892.6813042992,0.8398191671957083,0.15506260852390352,1.8283155844987502,skewnorm +44.92,11.220184543019636,14.125375446227547,281239.346976689,0.8301110671740033,0.1813188611337919,1.9999999999999998,skewnorm +44.92,14.125375446227547,17.782794100389236,261789.34549122906,0.9002914032046566,0.07589712356283235,,moyal +44.92,17.782794100389236,22.38721138568341,270142.01604579587,0.7822398416014035,0.21106206696985075,1.9999999990982784,skewnorm +44.92,22.38721138568341,28.183829312644566,304379.3406700574,0.8660053153372104,0.0867452990979236,,moyal +44.92,28.183829312644566,35.48133892335755,226917.50838121373,0.8075675653301224,0.22572278815378857,1.6258598662618187,skewnorm +52.16,0.011220184543019636,0.014125375446227547,0.0,,,,-- +52.16,0.014125375446227547,0.01778279410038923,0.0,,,,-- +52.16,0.01778279410038923,0.0223872113856834,10.985986853869262,,,,-- +52.16,0.0223872113856834,0.02818382931264455,49.78993688630152,,,,-- +52.16,0.02818382931264455,0.03548133892335755,167.15151146078102,3.209529610369136,0.3179905714056528,,moyal +52.16,0.03548133892335755,0.04466835921509633,501.9662953759566,2.5372151158829586,0.255043620241555,,moyal +52.16,0.04466835921509633,0.05623413251903491,1222.4829232534742,2.070442733968691,0.19111509845989538,,moyal +52.16,0.05623413251903491,0.0707945784384138,2643.77037510025,1.6737730755869502,0.17436096675283563,,moyal +52.16,0.0707945784384138,0.08912509381337459,5313.92417269014,1.3879635892914277,0.16417714819840085,,moyal +52.16,0.08912509381337459,0.11220184543019636,9382.285957365573,1.1580159983950304,0.15462117639712644,,moyal +52.16,0.11220184543019636,0.14125375446227548,14702.112563732522,1.003026434942104,0.15744552109930818,,moyal +52.16,0.14125375446227548,0.17782794100389238,23413.144617459187,0.7894690569382673,0.41196462339200823,1.9999999999999998,skewnorm +52.16,0.17782794100389238,0.223872113856834,33825.8535230635,0.7200645359543121,0.42215028995633996,1.9999999999999096,skewnorm +52.16,0.223872113856834,0.2818382931264455,47106.81250965088,0.9850370270646628,0.27010810733208324,0.0004064703871002305,skewnorm +52.16,0.2818382931264455,0.35481338923357547,67504.83436806753,0.7330032201992895,0.36110750458338353,1.4196361982197958,skewnorm +52.16,0.35481338923357547,0.44668359215096326,85005.47133047853,0.7964474835146724,0.30393910369168525,0.9283112775571716,skewnorm +52.16,0.44668359215096326,0.5623413251903494,109719.91202834771,0.7737958846925653,0.30163711511733904,1.3775836328762299,skewnorm +52.16,0.5623413251903494,0.707945784384138,139449.33358281778,0.8195941905559551,0.2539323583585262,1.2552897208013634,skewnorm +52.16,0.707945784384138,0.8912509381337459,160642.33322050102,0.7975274781903083,0.25829259842426455,1.5899273096197315,skewnorm +52.16,0.8912509381337459,1.1220184543019636,195629.47186924817,0.8223922822931042,0.22787934914235672,1.476499943470567,skewnorm +52.16,1.1220184543019636,1.4125375446227548,229230.95826621138,0.79796025808163,0.23619064969179815,1.9020954306189681,skewnorm +52.16,1.4125375446227548,1.7782794100389236,257082.0873762013,0.8434960166618213,0.19579002677163782,1.2982885560913746,skewnorm +52.16,1.7782794100389236,2.238721138568341,292446.9700499994,0.8240242356050972,0.20931421829978178,1.7454572699576023,skewnorm +52.16,2.238721138568341,2.8183829312644564,311602.02168010396,0.8852697305544086,0.07480775118351811,,moyal +52.16,2.8183829312644564,3.548133892335755,348928.78017438145,0.8252273605193622,0.19668456282077587,1.8404214678926494,skewnorm +52.16,3.548133892335755,4.466835921509633,363213.1665667945,0.8384525772705712,0.18904816984914094,1.6872402293079167,skewnorm +52.16,4.466835921509633,5.623413251903494,393236.60623841925,0.8476441818653986,0.17987993402423513,1.5701200071011177,skewnorm +52.16,5.623413251903494,7.079457843841384,420710.3711230494,0.865821991563782,0.1634599259142811,1.2868255094265313,skewnorm +52.16,7.079457843841384,8.912509381337463,437360.0142131806,0.8436135579156928,0.17285193779753855,1.6732345014874799,skewnorm +52.16,8.912509381337463,11.220184543019636,407997.2942727642,0.9754193759750935,0.13859324728326633,-0.00017653285169396142,skewnorm +52.16,11.220184543019636,14.125375446227547,456798.3111503883,0.9771709229908847,0.12214432690090965,1.8896006824992093e-05,skewnorm +52.16,14.125375446227547,17.782794100389236,394436.8754040832,0.8816577527100793,0.16212963207173353,1.0461764541753529,skewnorm +52.16,17.782794100389236,22.38721138568341,387805.33594158443,0.9112228326555228,0.15431741712795413,0.6684311195008507,skewnorm +52.16,22.38721138568341,28.183829312644566,396936.4412880153,0.876818713255631,0.16631945291524614,1.2273351076576502,skewnorm +52.16,28.183829312644566,35.48133892335755,330820.6997661303,1.1031060721224308,0.18712397117324234,-1.9999999999999998,skewnorm +59.34,0.011220184543019636,0.014125375446227547,0.0,,,,-- +59.34,0.014125375446227547,0.01778279410038923,0.0,,,,-- +59.34,0.01778279410038923,0.0223872113856834,0.0,,,,-- +59.34,0.0223872113856834,0.02818382931264455,0.0,,,,-- +59.34,0.02818382931264455,0.03548133892335755,6.594480616212384,,,,-- +59.34,0.03548133892335755,0.04466835921509633,34.86822874640203,,,,-- +59.34,0.04466835921509633,0.05623413251903491,117.05733128825993,3.340660967663978,0.427088026216017,,moyal +59.34,0.05623413251903491,0.0707945784384138,402.62619100583174,2.6361907019005133,0.23459930883315538,,moyal +59.34,0.0707945784384138,0.08912509381337459,1053.5076939407877,2.1366456164998087,0.20915251884764305,,moyal +59.34,0.08912509381337459,0.11220184543019636,2744.330961270939,1.7391474175522885,0.16718352252350757,,moyal +59.34,0.11220184543019636,0.14125375446227548,5686.421501999154,1.4188184278513696,0.16063985641117437,,moyal +59.34,0.14125375446227548,0.17782794100389238,11323.168097572157,1.192066867059588,0.14843041763573842,,moyal +59.34,0.17782794100389238,0.223872113856834,19597.537455586258,1.0251212100061018,0.15369548479653344,,moyal +59.34,0.223872113856834,0.2818382931264455,33105.31116777264,0.8356091705578744,0.36993131141815194,1.9999999999999998,skewnorm +59.34,0.2818382931264455,0.35481338923357547,48166.086420815285,0.8458505929064776,0.14696436447754113,,moyal +59.34,0.35481338923357547,0.44668359215096326,68790.03414111596,1.0107345373115948,0.24698944715322577,9.924687144840452e-05,skewnorm +59.34,0.44668359215096326,0.5623413251903494,97533.8406769678,0.7575384181285996,0.34093550940482165,1.4847349064204867,skewnorm +59.34,0.5623413251903494,0.707945784384138,133945.57596207122,0.8142846555296939,0.28963481408027186,0.9776766447228661,skewnorm +59.34,0.707945784384138,0.8912509381337459,172602.9901079089,0.8150466502655557,0.26729733303888603,1.1627505290071634,skewnorm +59.34,0.8912509381337459,1.1220184543019636,211163.33825098426,0.9748259786944541,0.19568131909237527,0.0001835253874288206,skewnorm +59.34,1.1220184543019636,1.4125375446227548,263108.57014328486,0.8130072671134978,0.2428241329439751,1.5908425151106045,skewnorm +59.34,1.4125375446227548,1.7782794100389236,304066.3937468299,0.8189107507077682,0.23075975469901566,1.55354056598251,skewnorm +59.34,1.7782794100389236,2.238721138568341,362658.46382779063,0.814295268885639,0.22861007452075677,1.631243645653721,skewnorm +59.34,2.238721138568341,2.8183829312644564,413606.862311286,0.8344162266525855,0.20362906016609164,1.4697238831634025,skewnorm +59.34,2.8183829312644564,3.548133892335755,451985.70143519825,0.9578690242381168,0.13973467590857497,0.00063115851350882,skewnorm +59.34,3.548133892335755,4.466835921509633,502434.5723172022,0.8279685673494052,0.19992521834319188,1.7008134564223027,skewnorm +59.34,4.466835921509633,5.623413251903494,528012.1764895428,0.8262625458369015,0.197935684742697,1.8479758819012846,skewnorm +59.34,5.623413251903494,7.079457843841384,606044.5214943996,0.8404535226211525,0.1907859258359077,1.6307660243829536,skewnorm +59.34,7.079457843841384,8.912509381337463,655957.6207555854,0.9590395261017436,0.1274320778736502,0.00036554988686013216,skewnorm +59.34,8.912509381337463,11.220184543019636,658973.0879647554,0.8436857441552686,0.1771195189900154,1.6786019973971422,skewnorm +59.34,11.220184543019636,14.125375446227547,693081.8451190108,0.8640194210911286,0.17604334383060455,1.4320478632538665,skewnorm +59.34,14.125375446227547,17.782794100389236,672912.1496397232,0.839847118199022,0.18613626190324123,1.714165628655833,skewnorm +59.34,17.782794100389236,22.38721138568341,641600.9077816121,0.8550484369776706,0.18811808318657436,1.9999999999999862,skewnorm +59.34,22.38721138568341,28.183829312644566,614962.5840342575,0.8431941735306215,0.18541260124695028,1.9694567483441339,skewnorm +59.34,28.183829312644566,35.48133892335755,684507.0879628478,0.8623239797088076,0.1708615193632962,1.5894207762028103,skewnorm +66.44,0.011220184543019636,0.014125375446227547,0.0,,,,-- +66.44,0.014125375446227547,0.01778279410038923,0.0,,,,-- +66.44,0.01778279410038923,0.0223872113856834,0.0,,,,-- +66.44,0.0223872113856834,0.02818382931264455,0.0,,,,-- +66.44,0.02818382931264455,0.03548133892335755,0.0,,,,-- +66.44,0.03548133892335755,0.04466835921509633,0.0,,,,-- +66.44,0.04466835921509633,0.05623413251903491,0.0,,,,-- +66.44,0.05623413251903491,0.0707945784384138,5.507685938437952,,,,-- +66.44,0.0707945784384138,0.08912509381337459,25.42380788963012,,,,-- +66.44,0.08912509381337459,0.11220184543019636,87.29093950044131,,,,-- +66.44,0.11220184543019636,0.14125375446227548,369.97236592019163,2.5868775763033183,0.44062169384326555,,moyal +66.44,0.14125375446227548,0.17782794100389238,1083.716723592641,1.977062496251662,0.6235073638051926,1.9999999999996878,skewnorm +66.44,0.17782794100389238,0.223872113856834,3344.031782850412,1.7319440134815496,0.20131059696746673,,moyal +66.44,0.223872113856834,0.2818382931264455,7856.99320047493,1.3402836820725519,0.433586530065737,1.9999999999999998,skewnorm +66.44,0.2818382931264455,0.35481338923357547,17408.808384654527,1.0934797479297846,0.387429510270976,1.9999999999999998,skewnorm +66.44,0.35481338923357547,0.44668359215096326,33222.09838614011,1.0642335987179057,0.1336177082901529,,moyal +66.44,0.44668359215096326,0.5623413251903494,53753.06638209194,0.9348010585178528,0.14532550676358008,,moyal +66.44,0.5623413251903494,0.707945784384138,79622.77971668489,0.7858701262471142,0.3846693364559444,1.9999999999999998,skewnorm +66.44,0.707945784384138,0.8912509381337459,113644.42126664656,0.7547877754928657,0.3723108778923076,1.9999999999998896,skewnorm +66.44,0.8912509381337459,1.1220184543019636,172079.53873520333,0.8004437519691939,0.3307593605626883,1.3537344557349742,skewnorm +66.44,1.1220184543019636,1.4125375446227548,240188.99042957384,0.7997299848990049,0.301140567335514,1.3431281270935835,skewnorm +66.44,1.4125375446227548,1.7782794100389236,315384.6243680885,0.8130771861500119,0.27344525322008156,1.3912272211057297,skewnorm +66.44,1.7782794100389236,2.238721138568341,388685.63864901825,0.9845484092849405,0.19293693464519746,0.00010052253876627766,skewnorm +66.44,2.238721138568341,2.8183829312644564,463745.3196001249,0.8296370821549627,0.2365744369839841,1.5054487882643852,skewnorm +66.44,2.8183829312644564,3.548133892335755,565326.2088547446,0.822616811143436,0.22973642285500084,1.5890805721422667,skewnorm +66.44,3.548133892335755,4.466835921509633,649614.8108419586,0.8269649019360225,0.22137280032568987,1.622304929376244,skewnorm +66.44,4.466835921509633,5.623413251903494,722881.036505778,0.8456594844721661,0.199628447836953,1.2789396069246413,skewnorm +66.44,5.623413251903494,7.079457843841384,853324.1413953209,0.853214291320395,0.18810262722378907,1.217705225319623,skewnorm +66.44,7.079457843841384,8.912509381337463,910865.699027567,0.832004117551249,0.19833356076427328,1.536654775217749,skewnorm +66.44,8.912509381337463,11.220184543019636,1011411.0190117816,0.849577655584695,0.18837428660185207,1.4114726756566975,skewnorm +66.44,11.220184543019636,14.125375446227547,1069989.3869830498,0.8588289829388677,0.1742588970321335,1.2699675521112113,skewnorm +66.44,14.125375446227547,17.782794100389236,1139977.7620089403,0.8440601747973817,0.17856494690510982,1.4449964542987168,skewnorm +66.44,17.782794100389236,22.38721138568341,1294651.1937077083,0.8437951061483908,0.1825837393344494,1.5404018277567963,skewnorm +66.44,22.38721138568341,28.183829312644566,1192070.3172069404,0.8486139612732978,0.18428997760921983,1.6722035275564837,skewnorm +66.44,28.183829312644566,35.48133892335755,1231130.3181114073,0.8418079971093185,0.19539161079854211,1.7343566461691893,skewnorm diff --git a/lstchain/data/LST1_gamma_irf_gheffi_0.70_theffi_0.70.csv b/lstchain/data/LST1_gamma_irf_gheffi_0.70_theffi_0.70.csv new file mode 100644 index 0000000000..880c44c6e2 --- /dev/null +++ b/lstchain/data/LST1_gamma_irf_gheffi_0.70_theffi_0.70.csv @@ -0,0 +1,351 @@ +ZD_deg,Etrue_min_TeV,Etrue_max_TeV,Aeff_m2,emig_mu_loc,emig_mu_scale,emig_mu_a,emig_model +6.0,0.011220184543019636,0.014125375446227547,1523.5581041209598,2.1610037007829788,0.24911300828762759,,moyal +6.0,0.014125375446227547,0.01778279410038923,3045.4471055995996,1.7884164513224,0.23669835172791936,,moyal +6.0,0.01778279410038923,0.0223872113856834,5419.198466082751,1.2022066181865934,0.6035684384735526,3.579431337301123,skewnorm +6.0,0.0223872113856834,0.02818382931264455,8858.63426071004,1.218007626486119,0.1921926525080171,,moyal +6.0,0.02818382931264455,0.03548133892335755,12992.827804060527,1.0563378246201154,0.18855286638193244,,moyal +6.0,0.03548133892335755,0.04466835921509633,18168.315998284892,0.6827788118420669,0.6123055842013111,6.380388591772884,skewnorm +6.0,0.04466835921509633,0.05623413251903491,24782.459317894332,0.6317298754630765,0.5573572953293215,2.058065840009622,skewnorm +6.0,0.05623413251903491,0.0707945784384138,34390.59879161294,0.6211980692549665,0.5022944499853085,2.053167793044295,skewnorm +6.0,0.0707945784384138,0.08912509381337459,44212.46990544327,0.6115496737979137,0.49810096284378047,1.9017488775179474,skewnorm +6.0,0.08912509381337459,0.11220184543019636,59401.676033561605,0.6846792221984536,0.40556258157561664,1.1669633146619578,skewnorm +6.0,0.11220184543019636,0.14125375446227548,73706.4846173022,0.6932053194233184,0.37473905225514026,1.218702327656692,skewnorm +6.0,0.14125375446227548,0.17782794100389238,87007.54531262298,0.7405223900162305,0.3256043583338413,0.9935886824790576,skewnorm +6.0,0.17782794100389238,0.223872113856834,97600.8703330211,0.7293930489117956,0.3325257530858587,1.3087090540533532,skewnorm +6.0,0.223872113856834,0.2818382931264455,109239.08521884098,0.7685678682783514,0.29071541892646463,1.1901186531318133,skewnorm +6.0,0.2818382931264455,0.35481338923357547,122332.69573806296,0.7814694646843571,0.27806185517810666,1.2386496799846154,skewnorm +6.0,0.35481338923357547,0.44668359215096326,145548.80680812232,0.7629713773545741,0.292986659742509,1.6086624829649936,skewnorm +6.0,0.44668359215096326,0.5623413251903494,149435.56739561658,0.8126906129907701,0.2648723725235948,1.208386336753635,skewnorm +6.0,0.5623413251903494,0.707945784384138,156943.5377057934,0.7884995425534035,0.25497853363248163,1.4943505315499954,skewnorm +6.0,0.707945784384138,0.8912509381337459,175565.69999796344,0.813781967651626,0.24576260507037234,1.4044317347918027,skewnorm +6.0,0.8912509381337459,1.1220184543019636,174833.08226052157,0.8193865981989024,0.23545159496477297,1.1686098185710727,skewnorm +6.0,1.1220184543019636,1.4125375446227548,191607.78446864747,0.8460610655238754,0.2108822667081557,1.054911011268091,skewnorm +6.0,1.4125375446227548,1.7782794100389236,204982.01672304966,0.8483509180118686,0.21816230238095358,1.1246555230516078,skewnorm +6.0,1.7782794100389236,2.238721138568341,212897.08259610785,1.1136998001698477,0.22467781288960492,-1.073466160930602,skewnorm +6.0,2.238721138568341,2.8183829312644564,190283.6959733396,1.1502348694342823,0.26737286420704415,-2.0855024409961285,skewnorm +6.0,2.8183829312644564,3.548133892335755,184046.80194622054,1.161278996767633,0.2836120603221115,-2.187563048557988,skewnorm +6.0,3.548133892335755,4.466835921509633,230781.7467393474,1.1365069136596004,0.21652212499355258,-1.865445521558706,skewnorm +6.0,4.466835921509633,5.623413251903494,232661.10802880663,1.1220558279927748,0.22589159751059873,-2.0642229588095224,skewnorm +6.0,5.623413251903494,7.079457843841384,243357.20331353246,1.148982284360027,0.26230018831682517,-2.0695106353204804,skewnorm +6.0,7.079457843841384,8.912509381337463,185288.77429252173,0.7141913612013501,0.3131577881344294,1.6568847423428024,skewnorm +6.0,8.912509381337463,11.220184543019636,205550.1231332423,1.1872527714174492,0.28647150131033045,-2.352973922437182,skewnorm +6.0,11.220184543019636,14.125375446227547,223881.62980403443,1.1091740288764065,0.26708535789685506,-1.77014986410527,skewnorm +6.0,14.125375446227547,17.782794100389236,212302.80303458704,0.7258627254629125,0.37391377341120946,1.3037848337501934,skewnorm +6.0,17.782794100389236,22.38721138568341,211975.45020391681,,,,-- +6.0,22.38721138568341,28.183829312644566,232053.2877723653,,,,-- +6.0,28.183829312644566,35.48133892335755,146068.8904335084,,,,-- +9.579,0.011220184543019636,0.014125375446227547,1385.4152640497937,2.1809344077359754,0.2734254843147913,,moyal +9.579,0.014125375446227547,0.01778279410038923,2877.360483731124,1.756367677326761,0.2234992186601322,,moyal +9.579,0.01778279410038923,0.0223872113856834,5065.293948337394,1.1957022847581742,0.6354066419082139,5.523925843179915,skewnorm +9.579,0.0223872113856834,0.02818382931264455,8225.522147977452,0.9679474126761177,0.6078808941622484,6.47848896843218,skewnorm +9.579,0.02818382931264455,0.03548133892335755,12616.401993360969,0.7913090598701766,0.5971001069254681,6.5356739962272545,skewnorm +9.579,0.03548133892335755,0.04466835921509633,17987.876010049164,0.6675862396061387,0.6135335999255693,6.2249832525051705,skewnorm +9.579,0.04466835921509633,0.05623413251903491,24729.003938545233,0.6319030533252338,0.5642371072583144,3.8442791530315628,skewnorm +9.579,0.05623413251903491,0.0707945784384138,33064.402076862534,0.6583613152116407,0.47668375770279237,1.3963771316607394,skewnorm +9.579,0.0707945784384138,0.08912509381337459,45649.425547041814,0.6290670626497145,0.4721007626326203,1.8460015692161895,skewnorm +9.579,0.08912509381337459,0.11220184543019636,60380.52995300011,0.6757443618843559,0.40557208131181355,1.1552367892102156,skewnorm +9.579,0.11220184543019636,0.14125375446227548,71733.72367191159,0.6646523804070077,0.3962344981671617,1.3426307304129808,skewnorm +9.579,0.14125375446227548,0.17782794100389238,87188.26764044777,0.7046288400087503,0.3521333055758652,1.1660688059769266,skewnorm +9.579,0.17782794100389238,0.223872113856834,101352.34955343003,0.7030334445106112,0.347714276206818,1.4219575674764349,skewnorm +9.579,0.223872113856834,0.2818382931264455,111067.95019868566,0.7503166762797937,0.2988376829097654,1.2033636367700091,skewnorm +9.579,0.2818382931264455,0.35481338923357547,124838.30343693434,0.7262577111841861,0.32096099625047164,2.0131284885323586,skewnorm +9.579,0.35481338923357547,0.44668359215096326,137644.33988102042,0.7763356616047778,0.2715142718726146,1.364038486505073,skewnorm +9.579,0.44668359215096326,0.5623413251903494,152447.86001293373,0.7745480719189156,0.2654833842869067,1.570425348578412,skewnorm +9.579,0.5623413251903494,0.707945784384138,165175.0574861938,0.7731756425072408,0.2665009969253811,1.7404334738700673,skewnorm +9.579,0.707945784384138,0.8912509381337459,172792.5570858644,0.7954127720862212,0.24951462268010488,1.5304333848320568,skewnorm +9.579,0.8912509381337459,1.1220184543019636,191680.52517384416,0.8285224324275979,0.23122778050740228,1.2296455886022302,skewnorm +9.579,1.1220184543019636,1.4125375446227548,196450.41839436235,0.8681627554751736,0.2014518908719833,0.8417472294392205,skewnorm +9.579,1.4125375446227548,1.7782794100389236,210588.08928397746,1.0816835407261822,0.21005597248816205,-1.0860445067409754,skewnorm +9.579,1.7782794100389236,2.238721138568341,217947.05153855373,1.1086323832700473,0.2230420564069321,-1.0952026728417374,skewnorm +9.579,2.238721138568341,2.8183829312644564,203005.41858806385,1.0961528907461568,0.22325047048115856,-1.4243791197814757,skewnorm +9.579,2.8183829312644564,3.548133892335755,239733.73314880376,0.8727787176820403,0.20500966497442114,0.9173817684741018,skewnorm +9.579,3.548133892335755,4.466835921509633,244783.46735324617,1.1692616813053232,0.26630692739994777,-2.751429097633909,skewnorm +9.579,4.466835921509633,5.623413251903494,259138.01625475715,1.1376757946839933,0.25005447997180513,-1.8645015112135765,skewnorm +9.579,5.623413251903494,7.079457843841384,199855.94144236934,1.1715862763494829,0.2903920461042626,-2.9725933370930395,skewnorm +9.579,7.079457843841384,8.912509381337463,210903.1210684001,1.1436400909594535,0.2744219704519533,-1.8137531279403298,skewnorm +9.579,8.912509381337463,11.220184543019636,230575.6013634338,1.1315755570684343,0.24928812328559113,-2.288859918274727,skewnorm +9.579,11.220184543019636,14.125375446227547,263888.6217237705,1.207014824965494,0.34237634315549614,-8.739473878013293,skewnorm +9.579,14.125375446227547,17.782794100389236,217786.1046056909,1.2373884469701888,0.4480182366212804,-5.31534230603821,skewnorm +9.579,17.782794100389236,22.38721138568341,232352.93340997197,1.178622245742845,0.3312996258660044,-4.143403837699235,skewnorm +9.579,22.38721138568341,28.183829312644566,225236.55952854347,1.1822944720223634,0.4694975395880055,-9.512508520457718,skewnorm +9.579,28.183829312644566,35.48133892335755,180444.74538082007,,,,-- +16.08,0.011220184543019636,0.014125375446227547,1087.7640050957623,2.301090239026353,0.25193534653277005,,moyal +16.08,0.014125375446227547,0.01778279410038923,2334.0123465610695,1.889348332742689,0.23035478605851228,,moyal +16.08,0.01778279410038923,0.0223872113856834,4489.141944416346,1.543850706468458,0.1971918812910741,,moyal +16.08,0.0223872113856834,0.02818382931264455,7701.571502665957,1.2958404868557518,0.18989059976985476,,moyal +16.08,0.02818382931264455,0.03548133892335755,11733.915866131078,1.105939581900147,0.18632101656754116,,moyal +16.08,0.03548133892335755,0.04466835921509633,16535.723529899475,0.9669904418815992,0.18462118050334209,,moyal +16.08,0.04466835921509633,0.05623413251903491,23542.53613863602,0.6252600498111353,0.5775442355887664,4.131468651427831,skewnorm +16.08,0.05623413251903491,0.0707945784384138,32622.271680035006,0.6090062088528094,0.5246975067215661,2.437670917123703,skewnorm +16.08,0.0707945784384138,0.08912509381337459,45364.53009871934,0.5949722992855373,0.49955542455969637,2.230093170555699,skewnorm +16.08,0.08912509381337459,0.11220184543019636,59574.97417539303,0.680763090862947,0.396457444021979,1.166174164959308,skewnorm +16.08,0.11220184543019636,0.14125375446227548,72010.98200016176,0.7032004690773048,0.34977951209604163,0.9861936881619721,skewnorm +16.08,0.14125375446227548,0.17782794100389238,89233.27690815777,0.6943980660850909,0.3464592178283736,1.1196833013703225,skewnorm +16.08,0.17782794100389238,0.223872113856834,101070.39313437612,0.7024787167599681,0.3403287760435061,1.2995769643812116,skewnorm +16.08,0.223872113856834,0.2818382931264455,114834.36584840689,0.6964128176737857,0.3347637597070631,1.7504083635475496,skewnorm +16.08,0.2818382931264455,0.35481338923357547,132515.31765999075,0.7316724287179386,0.3007703738173642,1.608724323006229,skewnorm +16.08,0.35481338923357547,0.44668359215096326,146013.25893805057,0.7487615330617691,0.27443651973775585,1.5254030404716388,skewnorm +16.08,0.44668359215096326,0.5623413251903494,158917.11947447987,0.7564377965891549,0.2619981303945415,1.525281403709265,skewnorm +16.08,0.5623413251903494,0.707945784384138,167895.79037177464,0.7699504387770432,0.25033621528632005,1.6614253782223634,skewnorm +16.08,0.707945784384138,0.8912509381337459,191237.9649339792,0.7755919984858726,0.25209228182889104,1.588616044392862,skewnorm +16.08,0.8912509381337459,1.1220184543019636,200745.01430395382,0.7871103156916746,0.24561917719447263,1.589150270246958,skewnorm +16.08,1.1220184543019636,1.4125375446227548,216924.76175524044,0.8658827361366912,0.19507454977199307,0.6795334591051769,skewnorm +16.08,1.4125375446227548,1.7782794100389236,216006.8336993507,0.9724877296074772,0.16826884556212207,0.0022496027355398303,skewnorm +16.08,1.7782794100389236,2.238721138568341,235107.61136832368,0.8475802105810858,0.20523396265969673,0.9195267397897997,skewnorm +16.08,2.238721138568341,2.8183829312644564,227062.2772995562,0.8446176355452749,0.2140457331215304,1.2236775537049644,skewnorm +16.08,2.8183829312644564,3.548133892335755,243260.26170441217,1.0902746246606216,0.2231207978098831,-1.0851610233905367,skewnorm +16.08,3.548133892335755,4.466835921509633,251431.9720015941,1.1164043150049032,0.22852106488511117,-1.3689830303425696,skewnorm +16.08,4.466835921509633,5.623413251903494,227524.51026081914,1.1203216832871103,0.2384800810007503,-1.37797374678181,skewnorm +16.08,5.623413251903494,7.079457843841384,226001.45760580045,1.1065709159325698,0.20734933124213287,-1.3671466643346013,skewnorm +16.08,7.079457843841384,8.912509381337463,239344.6556571148,1.1321385833537774,0.29170482084929816,-2.561210824594165,skewnorm +16.08,8.912509381337463,11.220184543019636,269389.43271405366,1.1324593860163694,0.25907319958671315,-1.339107795639018,skewnorm +16.08,11.220184543019636,14.125375446227547,257496.0982039625,1.1216069790076058,0.28003004207118576,-2.0841650544716117,skewnorm +16.08,14.125375446227547,17.782794100389236,253009.4684619041,1.1572939399106388,0.3539149405381489,-2.3996693919803143,skewnorm +16.08,17.782794100389236,22.38721138568341,234908.53633752913,1.0772198864040747,0.2145195484323544,-1.1633219046158063,skewnorm +16.08,22.38721138568341,28.183829312644566,248114.4089697137,1.1651268153607792,0.37017590867578953,-4.655746264617402,skewnorm +16.08,28.183829312644566,35.48133892335755,230324.24259950835,1.1499983537531582,0.41631616892025314,-284654.08744137455,skewnorm +23.16,0.011220184543019636,0.014125375446227547,748.5731183634505,2.5988447041573615,0.34595998404983125,,moyal +23.16,0.014125375446227547,0.01778279410038923,1806.96726452734,2.0965720028619192,0.2825598636352019,,moyal +23.16,0.01778279410038923,0.0223872113856834,3606.0852564195357,1.7179987850074505,0.24607261726608912,,moyal +23.16,0.0223872113856834,0.02818382931264455,6280.090798650582,1.4167470155925492,0.20623758491455402,,moyal +23.16,0.02818382931264455,0.03548133892335755,10110.575087495261,1.2078386517484736,0.20160129273502955,,moyal +23.16,0.03548133892335755,0.04466835921509633,15656.643169237048,1.036795181646683,0.18897873033448886,,moyal +23.16,0.04466835921509633,0.05623413251903491,22589.013300025767,0.921403396448888,0.19163090631876406,,moyal +23.16,0.05623413251903491,0.0707945784384138,30896.39240364059,0.5992390274505033,0.5935000277590965,4.465286366012604,skewnorm +23.16,0.0707945784384138,0.08912509381337459,42890.17884309136,0.5882465825833534,0.534586033778269,2.241995574554447,skewnorm +23.16,0.08912509381337459,0.11220184543019636,57838.634714061525,0.7227959347813532,0.381980899994428,0.9087715975697437,skewnorm +23.16,0.11220184543019636,0.14125375446227548,73513.37446046817,0.6715775994513239,0.3947108291464402,1.2965402501577916,skewnorm +23.16,0.14125375446227548,0.17782794100389238,94425.88298995377,0.6812700904652782,0.3701032982660915,1.2923679616348134,skewnorm +23.16,0.17782794100389238,0.223872113856834,105775.66085870904,0.7255935996742419,0.3395746717439297,1.042062417614452,skewnorm +23.16,0.223872113856834,0.2818382931264455,120480.13645745981,0.7697969413799999,0.2914639746927382,0.9279488976478885,skewnorm +23.16,0.2818382931264455,0.35481338923357547,139623.83493115334,0.7424165006689969,0.29609340013480007,1.3152901707077258,skewnorm +23.16,0.35481338923357547,0.44668359215096326,155036.46656182603,0.7412038546361884,0.29643556333585225,1.8012244275492846,skewnorm +23.16,0.44668359215096326,0.5623413251903494,165110.59934599564,0.7781934276178196,0.26834557775154083,1.3540957392348258,skewnorm +23.16,0.5623413251903494,0.707945784384138,185701.66526452443,0.7664444468590713,0.2654068255442937,1.8596581443628069,skewnorm +23.16,0.707945784384138,0.8912509381337459,202070.57371112428,0.768048001367066,0.2584320880206915,1.9067134961060965,skewnorm +23.16,0.8912509381337459,1.1220184543019636,214893.2662964302,0.8017305999229914,0.238676590515113,1.4605863616834542,skewnorm +23.16,1.1220184543019636,1.4125375446227548,240430.396544383,0.7943939143096558,0.24478576428769108,1.7566120965976957,skewnorm +23.16,1.4125375446227548,1.7782794100389236,238560.28492355297,0.8389341392719594,0.21404612587891456,0.9799415339197911,skewnorm +23.16,1.7782794100389236,2.238721138568341,250700.6702188889,0.8099830014709022,0.23044767077857428,1.662104584105471,skewnorm +23.16,2.238721138568341,2.8183829312644564,258765.48337736196,0.8234476154564766,0.21976999830071822,1.480258967282898,skewnorm +23.16,2.8183829312644564,3.548133892335755,253186.3544447903,0.8814651728666726,0.19088505874758996,0.9663107140194516,skewnorm +23.16,3.548133892335755,4.466835921509633,286868.461977063,0.9161694228853873,0.17539374961655121,0.57222238051741,skewnorm +23.16,4.466835921509633,5.623413251903494,277681.1440718462,0.8534741155577742,0.20815605404673213,1.5004434464907905,skewnorm +23.16,5.623413251903494,7.079457843841384,311186.6860798336,1.124203793761266,0.22847325376968644,-1.4862011352277373,skewnorm +23.16,7.079457843841384,8.912509381337463,321379.55281794455,1.1099871324227155,0.2217515963212223,-1.3615630947290076,skewnorm +23.16,8.912509381337463,11.220184543019636,279692.1796804014,0.8457740416151189,0.2008888572672537,0.891791593369608,skewnorm +23.16,11.220184543019636,14.125375446227547,294322.2853170482,1.12282620233425,0.2484538604172581,-1.9234458589907615,skewnorm +23.16,14.125375446227547,17.782794100389236,263939.0386388244,1.1622313487742792,0.30890102957407706,-4.6434549243412455,skewnorm +23.16,17.782794100389236,22.38721138568341,276899.63575578213,1.131839722717615,0.32150208774470945,-2.5074210097259875,skewnorm +23.16,22.38721138568341,28.183829312644566,305691.86637325655,1.137934777845102,0.2949213940214973,-2.189567970345534,skewnorm +23.16,28.183829312644566,35.48133892335755,300447.80727448466,1.1057091326108628,0.25759459007259666,-1.7510357897902955,skewnorm +30.39,0.011220184543019636,0.014125375446227547,429.7521382490461,3.125266714740322,0.38420691083531533,,moyal +30.39,0.014125375446227547,0.01778279410038923,1011.2210879927586,2.4897851344760094,0.2826554827626241,,moyal +30.39,0.01778279410038923,0.0223872113856834,2206.4968529488424,2.009516021461368,0.25231444882660253,,moyal +30.39,0.0223872113856834,0.02818382931264455,4683.830726597316,1.6376228269760484,0.21554513521871804,,moyal +30.39,0.02818382931264455,0.03548133892335755,8299.132143021807,1.3825134358372952,0.19479647927322982,,moyal +30.39,0.03548133892335755,0.04466835921509633,13207.620056358757,1.159334321009168,0.1802291809453713,,moyal +30.39,0.04466835921509633,0.05623413251903491,19543.99608295644,1.0073054214881514,0.184235588766871,,moyal +30.39,0.05623413251903491,0.0707945784384138,28156.207158236437,0.6464636654307794,0.5691155029224217,5.518320475549248,skewnorm +30.39,0.0707945784384138,0.08912509381337459,39917.983037125116,0.6215010296879273,0.5226448602544128,2.89346496669072,skewnorm +30.39,0.08912509381337459,0.11220184543019636,55008.14007624156,0.6338049458776338,0.46809897986108423,1.7578868720917564,skewnorm +30.39,0.11220184543019636,0.14125375446227548,75296.40528096104,0.6552017543376776,0.4166160136288835,1.3759897528104381,skewnorm +30.39,0.14125375446227548,0.17782794100389238,95485.79453833243,0.6979472586156862,0.35671390245837004,1.073305597135119,skewnorm +30.39,0.17782794100389238,0.223872113856834,115105.902036463,0.6888624930092364,0.35123942854239787,1.167466809409078,skewnorm +30.39,0.223872113856834,0.2818382931264455,128906.85584623017,0.7166413699183866,0.32622412414455837,1.217528493028233,skewnorm +30.39,0.2818382931264455,0.35481338923357547,144723.63407378254,0.724277964807403,0.3069382302991543,1.4015495043814616,skewnorm +30.39,0.35481338923357547,0.44668359215096326,172581.08293763123,0.750046172371181,0.28307355261942496,1.4206292945827386,skewnorm +30.39,0.44668359215096326,0.5623413251903494,186289.88219526107,0.742951930679233,0.2824394220040217,1.6944021847237816,skewnorm +30.39,0.5623413251903494,0.707945784384138,209386.66733911246,0.7547683475597206,0.2675035394439845,1.8654534789282975,skewnorm +30.39,0.707945784384138,0.8912509381337459,229688.9052334862,0.7522274964508803,0.26663258474551343,2.0669425673789408,skewnorm +30.39,0.8912509381337459,1.1220184543019636,239145.1552147753,0.7657456601586684,0.2582587573285203,2.0094444411100163,skewnorm +30.39,1.1220184543019636,1.4125375446227548,268864.4296677877,0.795386081178956,0.23159418146513042,1.6716881548384643,skewnorm +30.39,1.4125375446227548,1.7782794100389236,279103.04842959554,0.7782654144078354,0.23836351586143986,1.8761415791905256,skewnorm +30.39,1.7782794100389236,2.238721138568341,278895.13102620764,0.796153663583384,0.23536410431612298,1.637073703655159,skewnorm +30.39,2.238721138568341,2.8183829312644564,335821.8256805626,0.8852232693285348,0.17183639461176298,0.5388389180225934,skewnorm +30.39,2.8183829312644564,3.548133892335755,315427.1600735472,0.8300312580654607,0.21323516493917868,1.392391588482517,skewnorm +30.39,3.548133892335755,4.466835921509633,288833.88085416233,0.845829202776021,0.19956923691234138,1.0750078017507363,skewnorm +30.39,4.466835921509633,5.623413251903494,304526.05454789416,1.0548855902189687,0.19783361013780965,-0.9699554656350708,skewnorm +30.39,5.623413251903494,7.079457843841384,350377.68975943106,1.0622979994049104,0.20003226827185933,-0.9776364313716532,skewnorm +30.39,7.079457843841384,8.912509381337463,311941.8541974018,1.1063217583649567,0.2381754532074322,-1.4829677586618086,skewnorm +30.39,8.912509381337463,11.220184543019636,360381.76471729606,1.1302254494963873,0.2524449029893568,-1.6450799439317456,skewnorm +30.39,11.220184543019636,14.125375446227547,328000.2392207207,0.8715587140661766,0.19365159040587607,0.6869526756130713,skewnorm +30.39,14.125375446227547,17.782794100389236,358674.5438782057,1.1448489131672996,0.289673983513845,-2.229034052852358,skewnorm +30.39,17.782794100389236,22.38721138568341,295970.343129853,1.1494109697059067,0.2880764505006238,-1.8859346871048825,skewnorm +30.39,22.38721138568341,28.183829312644566,377381.5679767345,1.1789417704031293,0.3769909054471943,-2.4594166980942345,skewnorm +30.39,28.183829312644566,35.48133892335755,333769.4448824,1.1750124285084407,0.3340807169138664,-2.295024598772383,skewnorm +37.66,0.011220184543019636,0.014125375446227547,148.12030299082653,3.834178410738438,0.6149999939300179,,moyal +37.66,0.014125375446227547,0.01778279410038923,413.8154928345777,3.0568293476626214,0.35465013910451004,,moyal +37.66,0.01778279410038923,0.0223872113856834,1072.475969430943,2.1528084773675222,0.7892852231579126,3.2234142935907117,skewnorm +37.66,0.0223872113856834,0.02818382931264455,2552.5650735907093,2.010726495256655,0.2541945283609711,,moyal +37.66,0.02818382931264455,0.03548133892335755,5252.181525791987,1.381331636220864,0.6226578560255046,3.7932641133887253,skewnorm +37.66,0.03548133892335755,0.04466835921509633,9460.988367892873,1.3699261886359424,0.1980579638719429,,moyal +37.66,0.04466835921509633,0.05623413251903491,15707.232748378183,1.1638666753793172,0.19496074012433875,,moyal +37.66,0.05623413251903491,0.0707945784384138,23821.622120310058,0.7896842214408688,0.5442874040019661,3.788969916626774,skewnorm +37.66,0.0707945784384138,0.08912509381337459,34048.00940909664,0.6998607717971581,0.5116505874356773,3.0689018282216454,skewnorm +37.66,0.08912509381337459,0.11220184543019636,50672.726081097,0.6314034817554596,0.5163417662772076,3.020868679276829,skewnorm +37.66,0.11220184543019636,0.14125375446227548,69261.4594875598,0.6227201327968686,0.4817533753385742,2.0146075628136564,skewnorm +37.66,0.14125375446227548,0.17782794100389238,93300.06713446246,0.7141326741463496,0.38210790717465376,1.0834717192877859,skewnorm +37.66,0.17782794100389238,0.223872113856834,115914.23183849592,0.7072342580936424,0.35907357886032276,1.1327415271578885,skewnorm +37.66,0.223872113856834,0.2818382931264455,142546.89332455027,0.7356270499672558,0.327618670244314,1.0263786009378262,skewnorm +37.66,0.2818382931264455,0.35481338923357547,164114.7453959263,0.7621648290296418,0.2948585365612845,1.0097436555091108,skewnorm +37.66,0.35481338923357547,0.44668359215096326,191305.0988055109,0.7480374566783932,0.29420940357901876,1.3958284897541358,skewnorm +37.66,0.44668359215096326,0.5623413251903494,216040.08220898616,0.7618868465861806,0.27892211853538784,1.5103444516296436,skewnorm +37.66,0.5623413251903494,0.707945784384138,233589.18254152563,0.7685047162030696,0.2585243731403705,1.5995860174925758,skewnorm +37.66,0.707945784384138,0.8912509381337459,271667.14031247614,0.7711052867690248,0.2555847659972887,1.7015073208159466,skewnorm +37.66,0.8912509381337459,1.1220184543019636,282374.73312367033,0.7840727834354554,0.24631528780705905,1.7041450710764268,skewnorm +37.66,1.1220184543019636,1.4125375446227548,314704.91772434517,0.7985227055421198,0.23286447541237512,1.5805807970856942,skewnorm +37.66,1.4125375446227548,1.7782794100389236,336672.11083702015,0.8057399720578351,0.22331025020408873,1.4746730691670735,skewnorm +37.66,1.7782794100389236,2.238721138568341,364352.40580667375,0.7970711848389367,0.23270683947431628,1.776239784104557,skewnorm +37.66,2.238721138568341,2.8183829312644564,392499.89513817517,0.8214524431779179,0.22690963684140236,1.530014962003408,skewnorm +37.66,2.8183829312644564,3.548133892335755,356006.6759597976,0.8349318027086765,0.21140004798928164,1.3540872901708283,skewnorm +37.66,3.548133892335755,4.466835921509633,404875.12106750754,0.8386332040213564,0.21421277133144162,1.6041910256561012,skewnorm +37.66,4.466835921509633,5.623413251903494,413178.17180742737,0.9126839958706042,0.16804632833711491,0.6042579245727575,skewnorm +37.66,5.623413251903494,7.079457843841384,418975.80847674503,0.8786489616298514,0.1868907491707069,0.7636172789335297,skewnorm +37.66,7.079457843841384,8.912509381337463,455125.6756884326,0.8087929232324071,0.24094424349035512,1.7431709493629992,skewnorm +37.66,8.912509381337463,11.220184543019636,402126.1182326246,0.8154428176883052,0.2591780302364128,2.3565443463162232,skewnorm +37.66,11.220184543019636,14.125375446227547,408852.61715961056,0.8500033864449226,0.1942238117409088,0.9743240706789107,skewnorm +37.66,14.125375446227547,17.782794100389236,486616.36657399323,1.1020819912110995,0.20389803066416448,-0.9658114086714837,skewnorm +37.66,17.782794100389236,22.38721138568341,440567.339766234,1.144599678805072,0.28715868125730987,-1.449356096379902,skewnorm +37.66,22.38721138568341,28.183829312644566,386629.6027405437,0.8606650019848778,0.21657252816048167,1.1341363532456639,skewnorm +37.66,28.183829312644566,35.48133892335755,366964.6480902711,1.138977862998708,0.22848775955262576,-1.4106680471761177,skewnorm +44.92,0.011220184543019636,0.014125375446227547,22.426172471971306,,,,-- +44.92,0.014125375446227547,0.01778279410038923,86.90547944585899,4.1570031251761135,0.44794781670220374,,moyal +44.92,0.01778279410038923,0.0223872113856834,329.57325957587136,3.3254734530568943,0.3900583451060812,,moyal +44.92,0.0223872113856834,0.02818382931264455,957.3495463533085,2.689128384333837,0.29120762336984557,,moyal +44.92,0.02818382931264455,0.03548133892335755,2294.863858345851,2.139378358472836,0.24955993854033676,,moyal +44.92,0.03548133892335755,0.04466835921509633,5204.085399473419,1.7727754940079692,0.21877453501578417,,moyal +44.92,0.04466835921509633,0.05623413251903491,9825.63661665443,1.4569726572740416,0.19733345146448975,,moyal +44.92,0.05623413251903491,0.0707945784384138,18050.59963881771,0.9987886848783454,0.5293085988216514,5.1703179567518465,skewnorm +44.92,0.0707945784384138,0.08912509381337459,28397.382318763615,1.056371404816738,0.17909024992553657,,moyal +44.92,0.08912509381337459,0.11220184543019636,41924.049104232996,0.9391012639866951,0.1904073280027344,,moyal +44.92,0.11220184543019636,0.14125375446227548,59656.83117687341,0.6545931907042506,0.5043360554180882,2.2504596557845584,skewnorm +44.92,0.14125375446227548,0.17782794100389238,83997.90105946033,0.685511892305617,0.43211633561486196,1.5298869625102482,skewnorm +44.92,0.17782794100389238,0.223872113856834,116471.73021814518,0.6959964511309046,0.39172065912310106,1.2476451844461325,skewnorm +44.92,0.223872113856834,0.2818382931264455,147938.56222511176,0.7203186330217469,0.3533225152859798,1.1190896693291,skewnorm +44.92,0.2818382931264455,0.35481338923357547,183867.4037251168,0.7320953680664217,0.32004449592454315,1.1155031316576713,skewnorm +44.92,0.35481338923357547,0.44668359215096326,215493.87288549676,0.786558050305293,0.27512441118382824,0.806802126828349,skewnorm +44.92,0.44668359215096326,0.5623413251903494,236378.48863339526,0.7519814842794806,0.28430597739459934,1.2940715772641351,skewnorm +44.92,0.5623413251903494,0.707945784384138,268025.34958254016,0.7816110678508952,0.259217099393065,1.1769364650835819,skewnorm +44.92,0.707945784384138,0.8912509381337459,301019.78076205007,0.7764095373609852,0.2521735828124721,1.431686254726771,skewnorm +44.92,0.8912509381337459,1.1220184543019636,348227.52881022857,0.7733148158877775,0.24958089688405494,1.6333202520436048,skewnorm +44.92,1.1220184543019636,1.4125375446227548,378991.3260622228,0.784909524542993,0.24275983753388247,1.6546982123381055,skewnorm +44.92,1.4125375446227548,1.7782794100389236,411889.0562872249,0.7978991088265631,0.22622853363822,1.4336313736450719,skewnorm +44.92,1.7782794100389236,2.238721138568341,437089.7819620977,0.7996266356385815,0.21894194334663938,1.7212489529983648,skewnorm +44.92,2.238721138568341,2.8183829312644564,473063.31047156407,0.7907206875302821,0.23804627682817872,1.9031482853966704,skewnorm +44.92,2.8183829312644564,3.548133892335755,494722.98289526853,0.8059513700286657,0.22454801510131764,1.716926338253847,skewnorm +44.92,3.548133892335755,4.466835921509633,507472.4395239999,0.8210561548771442,0.20136880761581466,1.525350570185971,skewnorm +44.92,4.466835921509633,5.623413251903494,531995.794713885,0.8074642604106723,0.2280844805738095,2.1763388930642598,skewnorm +44.92,5.623413251903494,7.079457843841384,610371.6725951381,0.8118428742762924,0.2178470112145251,2.2157022496396475,skewnorm +44.92,7.079457843841384,8.912509381337463,590424.6503692953,0.8902164553743854,0.09519624440783009,,moyal +44.92,8.912509381337463,11.220184543019636,611970.6182823622,0.805062196161565,0.22633972173732167,2.706426721654225,skewnorm +44.92,11.220184543019636,14.125375446227547,557365.2512810747,0.7882886499593442,0.2740620606958969,2.9978055305853277,skewnorm +44.92,14.125375446227547,17.782794100389236,554692.9984383829,0.7921789966087043,0.2653493540623421,2.24579267983445,skewnorm +44.92,17.782794100389236,22.38721138568341,561895.3933752554,0.7732828924971754,0.23805028453203522,2.779898893830273,skewnorm +44.92,22.38721138568341,28.183829312644566,686979.0705625877,0.7939530169591409,0.24273066621133635,2.2983579218668013,skewnorm +44.92,28.183829312644566,35.48133892335755,612249.1263870483,1.0760648152586685,0.24417887078138084,-1.1306454411259086,skewnorm +52.16,0.011220184543019636,0.014125375446227547,,,,,-- +52.16,0.014125375446227547,0.01778279410038923,4.907019979208879,,,,-- +52.16,0.01778279410038923,0.0223872113856834,34.05655924699471,,,,-- +52.16,0.0223872113856834,0.02818382931264455,138.30538023972645,3.803294577586596,0.4050345095467987,,moyal +52.16,0.02818382931264455,0.03548133892335755,513.6426654263583,2.651552419497247,1.0225503082540943,6.991187195498465,skewnorm +52.16,0.03548133892335755,0.04466835921509633,1306.4275635112233,2.4861419838859913,0.26166882398643454,,moyal +52.16,0.04466835921509633,0.05623413251903491,3416.3292527941335,1.7889484948381744,0.5998243263116124,3.5139611964490984,skewnorm +52.16,0.05623413251903491,0.0707945784384138,7653.385198877596,1.6605116631625507,0.19527479679308307,,moyal +52.16,0.0707945784384138,0.08912509381337459,16186.694126029803,1.3774977396067722,0.18807844624254275,,moyal +52.16,0.08912509381337459,0.11220184543019636,29237.053071368075,1.1531864683111934,0.16683692716658538,,moyal +52.16,0.11220184543019636,0.14125375446227548,46490.838738780774,0.7824689643983869,0.5113630807525696,5.621751968306168,skewnorm +52.16,0.14125375446227548,0.17782794100389238,70143.44257738984,0.6779629452266029,0.5206436296204671,4.967531831526724,skewnorm +52.16,0.17782794100389238,0.223872113856834,100291.07398897262,0.6368821085789884,0.49602407495110384,3.3018864014166858,skewnorm +52.16,0.223872113856834,0.2818382931264455,134778.59304361357,0.7066729537595312,0.39705802245084393,1.2551155012039725,skewnorm +52.16,0.2818382931264455,0.35481338923357547,182264.79395536004,0.7030764508179057,0.3767409049984356,1.4046409157840893,skewnorm +52.16,0.35481338923357547,0.44668359215096326,221939.2463179719,0.7579360927193716,0.32050496635981557,0.97560087415231,skewnorm +52.16,0.44668359215096326,0.5623413251903494,271153.8872209619,0.7584431892826112,0.300262765659277,1.096729902116445,skewnorm +52.16,0.5623413251903494,0.707945784384138,321768.74131640716,0.792936648200749,0.26358415046484457,1.0472399296225923,skewnorm +52.16,0.707945784384138,0.8912509381337459,357060.7156036402,0.7792547986767511,0.266963352795656,1.2853987262073725,skewnorm +52.16,0.8912509381337459,1.1220184543019636,409704.16554435005,0.7915708050974571,0.24669679254539667,1.3749369289399007,skewnorm +52.16,1.1220184543019636,1.4125375446227548,469899.203533912,0.7737645713951308,0.2504197259821984,1.671628381077058,skewnorm +52.16,1.4125375446227548,1.7782794100389236,512418.87884353497,0.8312370863960871,0.20287726731096217,0.9633593006258938,skewnorm +52.16,1.7782794100389236,2.238721138568341,566986.7815281921,0.812680595368915,0.21052300810083027,1.2312417428251217,skewnorm +52.16,2.238721138568341,2.8183829312644564,611586.3914200709,0.8323993767887667,0.19378139813724454,1.051678259784677,skewnorm +52.16,2.8183829312644564,3.548133892335755,661119.0510589452,0.8326668616817205,0.19342062175225838,1.1387441445452013,skewnorm +52.16,3.548133892335755,4.466835921509633,691573.6514895816,0.8233286857520327,0.19437189182664463,1.2192805055076212,skewnorm +52.16,4.466835921509633,5.623413251903494,727694.6881759379,0.8284646093416488,0.1910978579760791,1.4219693261663227,skewnorm +52.16,5.623413251903494,7.079457843841384,781666.6680651207,0.8422339739233372,0.1781696160948594,1.183344372680759,skewnorm +52.16,7.079457843841384,8.912509381337463,814801.7064791555,0.834336218751661,0.18648959951405963,1.29246356797301,skewnorm +52.16,8.912509381337463,11.220184543019636,814342.7776375415,0.8771322064891557,0.1729018239403178,0.7015358668816518,skewnorm +52.16,11.220184543019636,14.125375446227547,859529.4473846457,0.8470937981117105,0.18955320536306433,1.1691594165643933,skewnorm +52.16,14.125375446227547,17.782794100389236,796727.5823980707,1.0575241497664376,0.18174342383905628,-1.0054380798977802,skewnorm +52.16,17.782794100389236,22.38721138568341,842625.1916917714,1.086793125599595,0.2154712786057985,-1.302796029180954,skewnorm +52.16,22.38721138568341,28.183829312644566,840896.7118575375,0.863587512633732,0.18419283074801518,0.9532980511310108,skewnorm +52.16,28.183829312644566,35.48133892335755,828792.9109930423,1.1011682323636653,0.24400602278927078,-1.937099005060706,skewnorm +59.34,0.011220184543019636,0.014125375446227547,,,,,-- +59.34,0.014125375446227547,0.01778279410038923,,,,,-- +59.34,0.01778279410038923,0.0223872113856834,,,,,-- +59.34,0.0223872113856834,0.02818382931264455,0.0,,,,-- +59.34,0.02818382931264455,0.03548133892335755,17.145649602152197,,,,-- +59.34,0.03548133892335755,0.04466835921509633,86.3403759434717,3.5650564572372634,1.9792715372644525,6.828474020689434,skewnorm +59.34,0.04466835921509633,0.05623413251903491,309.3658041189726,3.317024957345085,0.3553052028611334,,moyal +59.34,0.05623413251903491,0.0707945784384138,1121.0376298593746,2.6117527068163064,0.2807165474597177,,moyal +59.34,0.0707945784384138,0.08912509381337459,3034.6322253137155,2.116472566595116,0.22510363129564084,,moyal +59.34,0.08912509381337459,0.11220184543019636,8045.310675215261,1.7232281240727814,0.19045220760605402,,moyal +59.34,0.11220184543019636,0.14125375446227548,17951.869912590126,1.4229429218188494,0.1767408737655012,,moyal +59.34,0.14125375446227548,0.17782794100389238,35595.59848536256,1.201330645995971,0.17264681372091675,,moyal +59.34,0.17782794100389238,0.223872113856834,63511.00036561967,1.0304187912454827,0.16475485103028512,,moyal +59.34,0.223872113856834,0.2818382931264455,100824.52996159617,0.9209027442676935,0.17140308220158926,,moyal +59.34,0.2818382931264455,0.35481338923357547,145078.57355667252,0.6599057248913055,0.48962218524659223,3.078253708560259,skewnorm +59.34,0.35481338923357547,0.44668359215096326,203497.62453138237,0.7094132146173161,0.3962414377185636,1.5206956182590514,skewnorm +59.34,0.44668359215096326,0.5623413251903494,260055.40331379263,0.7225914714709202,0.36149371780425,1.424511999924881,skewnorm +59.34,0.5623413251903494,0.707945784384138,334758.67815589544,0.7739355230647932,0.3144187865554494,0.9929696724525231,skewnorm +59.34,0.707945784384138,0.8912509381337459,408085.1501246107,0.7904765082886667,0.2792270100690727,1.0425786877670382,skewnorm +59.34,0.8912509381337459,1.1220184543019636,463241.39797623584,0.7967988923996732,0.26621556165382987,1.1226710504883066,skewnorm +59.34,1.1220184543019636,1.4125375446227548,551682.6474746546,0.7931404329977362,0.2553873441334107,1.3561186918160184,skewnorm +59.34,1.4125375446227548,1.7782794100389236,624459.8308100656,0.7869496147956215,0.2526076495113524,1.4280769248998777,skewnorm +59.34,1.7782794100389236,2.238721138568341,712834.419722087,0.7921827617713478,0.2432832981795224,1.3888077048045804,skewnorm +59.34,2.238721138568341,2.8183829312644564,792432.1951678235,0.8249072318090924,0.20722256110201523,1.0593650988518517,skewnorm +59.34,2.8183829312644564,3.548133892335755,870075.7725030647,0.8235535859852416,0.20173695979130093,1.0560731194625608,skewnorm +59.34,3.548133892335755,4.466835921509633,940611.9802303207,0.815517911856439,0.20328998848012872,1.2299308524598265,skewnorm +59.34,4.466835921509633,5.623413251903494,976174.5305645941,0.8177177071880696,0.1966220174519741,1.2149206467796407,skewnorm +59.34,5.623413251903494,7.079457843841384,1116564.0055148667,0.8253614916626406,0.19744660354719645,1.269169376796332,skewnorm +59.34,7.079457843841384,8.912509381337463,1179729.8421770907,0.8117816716771501,0.2035003640017601,1.5338370564368875,skewnorm +59.34,8.912509381337463,11.220184543019636,1176558.9121193513,0.8338631050916945,0.1883537629310276,1.2832041662668077,skewnorm +59.34,11.220184543019636,14.125375446227547,1245972.135202585,0.8588166321483834,0.1855077599772143,0.9768109308006145,skewnorm +59.34,14.125375446227547,17.782794100389236,1314095.6321058641,0.8258022514978525,0.2064722642007338,1.4745860374143267,skewnorm +59.34,17.782794100389236,22.38721138568341,1248250.79334944,0.8379634837411614,0.20582406441369513,1.3868476211724134,skewnorm +59.34,22.38721138568341,28.183829312644566,1278116.44381907,0.8306363732488117,0.20094690599919762,1.4720254697560444,skewnorm +59.34,28.183829312644566,35.48133892335755,1365057.487555968,0.8414326166446661,0.19548183837125074,1.3718108661082955,skewnorm +66.44,0.011220184543019636,0.014125375446227547,,,,,-- +66.44,0.014125375446227547,0.01778279410038923,,,,,-- +66.44,0.01778279410038923,0.0223872113856834,,,,,-- +66.44,0.0223872113856834,0.02818382931264455,,,,,-- +66.44,0.02818382931264455,0.03548133892335755,,,,,-- +66.44,0.03548133892335755,0.04466835921509633,,,,,-- +66.44,0.04466835921509633,0.05623413251903491,0.0,,,,-- +66.44,0.05623413251903491,0.0707945784384138,11.015371876875903,,,,-- +66.44,0.0707945784384138,0.08912509381337459,76.27142366889035,,,,-- +66.44,0.08912509381337459,0.11220184543019636,261.87281850132393,3.316838380433461,0.37257804948534806,,moyal +66.44,0.11220184543019636,0.14125375446227548,948.7410175577193,2.676091034065133,0.29299793872509294,,moyal +66.44,0.14125375446227548,0.17782794100389238,3260.373291829776,2.158560421421704,0.23172906902656704,,moyal +66.44,0.17782794100389238,0.223872113856834,9724.397979643125,1.51445619094617,0.5761754230540027,3.048534389234322,skewnorm +66.44,0.223872113856834,0.2818382931264455,24776.93669731164,1.255656149730818,0.4960158393795277,3.989742435288241,skewnorm +66.44,0.2818382931264455,0.35481338923357547,54342.71792799664,1.2163877682702602,0.16428819911908038,,moyal +66.44,0.35481338923357547,0.44668359215096326,102712.81254878116,1.0459814859569239,0.15181665140079834,,moyal +66.44,0.44668359215096326,0.5623413251903494,165925.770291764,0.7358544188204161,0.4801607699963378,4.836162008698918,skewnorm +66.44,0.5623413251903494,0.707945784384138,241475.31049424864,0.692698481818057,0.45940415292542375,2.820896062831038,skewnorm +66.44,0.707945784384138,0.8912509381337459,326695.93138174684,0.7014230007052407,0.4257467259762937,2.258965135321032,skewnorm +66.44,0.8912509381337459,1.1220184543019636,446405.86460525694,0.8033532401173726,0.3184944594471217,0.9595825127377953,skewnorm +66.44,1.1220184543019636,1.4125375446227548,564958.792038328,0.780218626858482,0.31015866973533834,1.1929663897863534,skewnorm +66.44,1.4125375446227548,1.7782794100389236,687952.5992576611,0.782169115180376,0.2926440284624247,1.3180258428165705,skewnorm +66.44,1.7782794100389236,2.238721138568341,801406.5057372738,0.812400047409304,0.26280364355763886,1.0683618784538238,skewnorm +66.44,2.238721138568341,2.8183829312644564,926613.5794941502,0.8036621973268265,0.2573073970319791,1.302243912835922,skewnorm +66.44,2.8183829312644564,3.548133892335755,1082253.722094646,0.8043034969832995,0.24309389843697213,1.2585125635306107,skewnorm +66.44,3.548133892335755,4.466835921509633,1224167.1399746467,0.833490842001681,0.2149816917581346,0.9284002736698557,skewnorm +66.44,4.466835921509633,5.623413251903494,1333326.8744749501,0.8756647334504378,0.18875816643727553,0.4979272864481302,skewnorm +66.44,5.623413251903494,7.079457843841384,1530218.7432293457,0.8507776903174427,0.19125016140553325,0.7393104494580769,skewnorm +66.44,7.079457843841384,8.912509381337463,1634519.7217769485,0.8232243453248254,0.2012453511148324,1.030990367629061,skewnorm +66.44,8.912509381337463,11.220184543019636,1776079.6490356487,0.8396896489514701,0.19554835367239953,0.9460005021915273,skewnorm +66.44,11.220184543019636,14.125375446227547,1919826.900779926,0.8382335263868046,0.18847249180142178,0.9149128639136531,skewnorm +66.44,14.125375446227547,17.782794100389236,1964985.9400971255,0.8253194427621874,0.19531310162329427,1.0615078185117401,skewnorm +66.44,17.782794100389236,22.38721138568341,2261285.380937903,0.8220225716625172,0.19812362472797898,1.2558386488945499,skewnorm +66.44,22.38721138568341,28.183829312644566,2116637.424053525,0.8179481406729284,0.20730632102172328,1.5188887453269042,skewnorm +66.44,28.183829312644566,35.48133892335755,2213826.267097194,0.8195622290147621,0.2084128854679719,1.4030579675347268,skewnorm diff --git a/lstchain/data/LST1_gamma_irf_gheffi_0.90_theffi_0.90.csv b/lstchain/data/LST1_gamma_irf_gheffi_0.90_theffi_0.90.csv new file mode 100644 index 0000000000..c48e8300ec --- /dev/null +++ b/lstchain/data/LST1_gamma_irf_gheffi_0.90_theffi_0.90.csv @@ -0,0 +1,351 @@ +ZD_deg,Etrue_min_TeV,Etrue_max_TeV,Aeff_m2,emig_mu_loc,emig_mu_scale,emig_mu_a,emig_model +6.0,0.011220184543019636,0.014125375446227547,1857.9267720101018,2.1274856783497613,0.26131298420369603,,moyal +6.0,0.014125375446227547,0.01778279410038923,3806.8088819994987,1.7520223108737538,0.24344766700522993,,moyal +6.0,0.01778279410038923,0.0223872113856834,7147.259201440771,1.4383113947449457,0.2125539533767747,,moyal +6.0,0.0223872113856834,0.02818382931264455,11498.240409120694,1.2201479704198464,0.19368082063219103,,moyal +6.0,0.02818382931264455,0.03548133892335755,17528.266852020948,1.0461833764696975,0.18615801365655982,,moyal +6.0,0.03548133892335755,0.04466835921509633,24926.267546229934,0.9382498912667445,0.18798165113063858,,moyal +6.0,0.04466835921509633,0.05623413251903491,34042.602921030935,0.855972826635985,0.18818435294548144,,moyal +6.0,0.05623413251903491,0.0707945784384138,45538.39882364001,0.7959899748134717,0.19027471464771684,,moyal +6.0,0.0707945784384138,0.08912509381337459,59017.226425647714,0.6294864312941103,0.4909973655898151,1.9999999999999998,skewnorm +6.0,0.08912509381337459,0.11220184543019636,75799.49484531462,0.6149882969547472,0.5002991384185688,1.9999999999999998,skewnorm +6.0,0.11220184543019636,0.14125375446227548,96879.68707883665,0.6384814466513199,0.46323026947425533,1.9999999999999634,skewnorm +6.0,0.14125375446227548,0.17782794100389238,116510.31414811904,0.6663412075379233,0.4094472665687403,1.7577258961990063,skewnorm +6.0,0.17782794100389238,0.223872113856834,136816.3286207458,0.9413099457161127,0.26944692882896015,3.479851717739268e-05,skewnorm +6.0,0.223872113856834,0.2818382931264455,158260.34226075315,0.7041141397081735,0.3616335923558678,1.7260348738428462,skewnorm +6.0,0.2818382931264455,0.35481338923357547,180541.1485758159,0.7332637730651687,0.33959718999122795,1.5727524191763416,skewnorm +6.0,0.35481338923357547,0.44668359215096326,210645.80947404183,0.7523584969141519,0.3152603453216642,1.3473443983039544,skewnorm +6.0,0.44668359215096326,0.5623413251903494,226642.01468676783,0.9597951475178,0.2354369187489256,0.00042342652330608295,skewnorm +6.0,0.5623413251903494,0.707945784384138,235779.61374927923,0.9556856426171145,0.21445643515604376,0.0012141315811311768,skewnorm +6.0,0.707945784384138,0.8912509381337459,263623.73134490434,0.7838956562490718,0.2818001222055325,1.4067157605404452,skewnorm +6.0,0.8912509381337459,1.1220184543019636,264212.7425443021,1.0739840842279667,0.25104905271166167,-0.8174713671568862,skewnorm +6.0,1.1220184543019636,1.4125375446227548,294825.9384692087,1.133329117986968,0.29650258394819523,-1.5075959615912513,skewnorm +6.0,1.4125375446227548,1.7782794100389236,310035.3002936126,1.1226153413710402,0.27963887446891605,-1.2572054599591804,skewnorm +6.0,1.7782794100389236,2.238721138568341,323492.969659021,1.1469992191418523,0.2931783115784595,-1.4880779909149784,skewnorm +6.0,2.238721138568341,2.8183829312644564,291226.87615431857,1.1601796349615887,0.3212764582288072,-1.9999999997847593,skewnorm +6.0,2.8183829312644564,3.548133892335755,294328.81422351935,1.1692366438561232,0.34717009636332974,-1.9999999999988114,skewnorm +6.0,3.548133892335755,4.466835921509633,335598.95442176017,1.1559943279676326,0.2874838428610312,-1.9999999999999998,skewnorm +6.0,4.466835921509633,5.623413251903494,358830.5646215426,1.144880038249868,0.3453081918819603,-1.9999999999999998,skewnorm +6.0,5.623413251903494,7.079457843841384,389080.0795491807,1.1536422847190457,0.32487493305753556,-1.9999999999999865,skewnorm +6.0,7.079457843841384,8.912509381337463,284354.05955783045,1.0813550470550433,0.307950935234577,-1.7213223081129738,skewnorm +6.0,8.912509381337463,11.220184543019636,344123.2398522821,1.1450617657287379,0.3288691678194212,-1.9999999999999998,skewnorm +6.0,11.220184543019636,14.125375446227547,427410.38417133846,1.098219893686063,0.4211743149566496,-1.9999999999999998,skewnorm +6.0,14.125375446227547,17.782794100389236,340416.56348649296,1.184572032930787,0.45249748697836806,-1.8108691338903162,skewnorm +6.0,17.782794100389236,22.38721138568341,387085.60472019593,1.0239437995132856,0.49498536815817873,-1.9999999999999998,skewnorm +6.0,22.38721138568341,28.183829312644566,417695.9179902575,1.062093028439496,0.502968400408199,-1.9999999999997988,skewnorm +6.0,28.183829312644566,35.48133892335755,306744.6699103677,1.062093028439496,0.502968400408199,-1.9999999999997988,skewnorm +9.579,0.011220184543019636,0.014125375446227547,1721.140232798368,2.129525350832714,0.274896846981797,,moyal +9.579,0.014125375446227547,0.01778279410038923,3565.7860516796227,1.7451525139237785,0.2263296582451902,,moyal +9.579,0.01778279410038923,0.0223872113856834,6626.705660852407,1.4491014055788123,0.21742174255038177,,moyal +9.579,0.0223872113856834,0.02818382931264455,10840.60635860755,1.2336740772811394,0.19966419557383383,,moyal +9.579,0.02818382931264455,0.03548133892335755,16810.821687008996,1.0496703892327155,0.18805875236923877,,moyal +9.579,0.03548133892335755,0.04466835921509633,24686.96900863706,0.9327600824205626,0.19128627770139095,,moyal +9.579,0.04466835921509633,0.05623413251903491,33436.97511539431,0.8622608677261774,0.1911863300420297,,moyal +9.579,0.05623413251903491,0.0707945784384138,44651.63809357637,0.8018798986298421,0.1953476938627192,,moyal +9.579,0.0707945784384138,0.08912509381337459,59829.88539782502,0.6493673107839161,0.46378281626546825,1.9999999999999998,skewnorm +9.579,0.08912509381337459,0.11220184543019636,77499.02156937621,0.610028349908749,0.4936114377030635,1.9999999999999893,skewnorm +9.579,0.11220184543019636,0.14125375446227548,94296.20082929397,0.6275434626055414,0.4537976832305107,1.9999999999999998,skewnorm +9.579,0.14125375446227548,0.17782794100389238,117438.38844118737,0.9357512997041263,0.2825962795227593,-0.00017411192863433795,skewnorm +9.579,0.17782794100389238,0.223872113856834,141572.64232669628,0.6598621536439198,0.40845408981745485,1.924226026693,skewnorm +9.579,0.223872113856834,0.2818382931264455,160824.75380362757,0.6833301916178697,0.380252598175973,1.8759572814446401,skewnorm +9.579,0.2818382931264455,0.35481338923357547,185563.48407631626,0.7090662656376407,0.3546507763637783,1.76691135577822,skewnorm +9.579,0.35481338923357547,0.44668359215096326,205655.20098087122,0.7391047953234584,0.3220372822730919,1.6588825901968582,skewnorm +9.579,0.44668359215096326,0.5623413251903494,225695.2046975554,0.7486617813576477,0.30611459133992,1.501265976140073,skewnorm +9.579,0.5623413251903494,0.707945784384138,247542.15688211162,0.7588092847641694,0.2961421068237873,1.4631853044460204,skewnorm +9.579,0.707945784384138,0.8912509381337459,261686.37258881715,0.9510299492565193,0.2143099615792647,0.000784090775541258,skewnorm +9.579,0.8912509381337459,1.1220184543019636,283444.95642353385,1.0790825605889491,0.24778214357420222,-0.7869912627950525,skewnorm +9.579,1.1220184543019636,1.4125375446227548,298194.14254786045,1.1336142755074836,0.28023222362857886,-1.5302930388306835,skewnorm +9.579,1.4125375446227548,1.7782794100389236,319665.7060822515,1.1091185056061676,0.2662195091784872,-1.189237092382598,skewnorm +9.579,1.7782794100389236,2.238721138568341,329244.1066419303,1.1697370889747705,0.3203535824440817,-1.8158325747889512,skewnorm +9.579,2.238721138568341,2.8183829312644564,313283.5782533377,1.1379241902126742,0.2956712998609832,-1.6685306044487227,skewnorm +9.579,2.8183829312644564,3.548133892335755,360152.9815968204,1.158094888008935,0.3072445055703104,-1.9999999999992524,skewnorm +9.579,3.548133892335755,4.466835921509633,396382.3192936088,1.1633921981974489,0.33262918064329533,-1.9999999999999998,skewnorm +9.579,4.466835921509633,5.623413251903494,400963.55217796884,1.1536875129381847,0.33556288858000766,-1.9999999999999998,skewnorm +9.579,5.623413251903494,7.079457843841384,329174.4917874319,1.1474029187428267,0.37123504661261225,-1.9999999999999998,skewnorm +9.579,7.079457843841384,8.912509381337463,343180.0785805984,1.151765688509001,0.38088545817194974,-1.9999999999999998,skewnorm +9.579,8.912509381337463,11.220184543019636,377305.52950380073,1.1378921210258781,0.34058453207780665,-1.9999999999999998,skewnorm +9.579,11.220184543019636,14.125375446227547,416357.6031641713,1.1581747260553457,0.4134466591285519,-1.9999999999999998,skewnorm +9.579,14.125375446227547,17.782794100389236,394968.02021710045,1.103916952943355,0.45787758107973603,-1.9999999999999978,skewnorm +9.579,17.782794100389236,22.38721138568341,394999.9867969523,1.1185862330013951,0.32630110230404946,-1.9999999999999998,skewnorm +9.579,22.38721138568341,28.183829312644566,427071.91806710843,1.0507613526089035,0.4367499981541915,-1.9999999999999998,skewnorm +9.579,28.183829312644566,35.48133892335755,382984.7657062304,1.062764977416507,0.4134110132192726,-1.9999999999999998,skewnorm +16.08,0.011220184543019636,0.014125375446227547,1395.5032444126925,2.271774557388867,0.24948342182370342,,moyal +16.08,0.014125375446227547,0.01778279410038923,2950.722925936962,1.8706954270742069,0.235515366363515,,moyal +16.08,0.01778279410038923,0.0223872113856834,5705.490382570841,1.5396577463743173,0.2036096996492253,,moyal +16.08,0.0223872113856834,0.02818382931264455,10167.678355455837,1.2891562165657435,0.19305137544303572,,moyal +16.08,0.02818382931264455,0.03548133892335755,16005.957297360299,1.0992296510155857,0.18295011281445875,,moyal +16.08,0.03548133892335755,0.04466835921509633,22922.8105911722,0.9666011392971333,0.18352138632682397,,moyal +16.08,0.04466835921509633,0.05623413251903491,32398.490153403312,0.8684832933216083,0.18747615635776493,,moyal +16.08,0.05623413251903491,0.0707945784384138,43758.66662445453,0.8049148691995047,0.18982422976485858,,moyal +16.08,0.0707945784384138,0.08912509381337459,59661.80688035625,0.7621682065741856,0.18635252328824664,,moyal +16.08,0.08912509381337459,0.11220184543019636,76985.51343783832,0.6162239346067288,0.4744919254635792,1.999999999999997,skewnorm +16.08,0.11220184543019636,0.14125375446227548,94180.76801625901,0.6406394111929151,0.4283390320167125,1.9999999999999998,skewnorm +16.08,0.14125375446227548,0.17782794100389238,121096.65684257883,0.6360694222471361,0.4285668946145254,1.9207407964931842,skewnorm +16.08,0.17782794100389238,0.223872113856834,142199.29449652327,0.6598339279849006,0.39412116568624783,1.7587691415190383,skewnorm +16.08,0.223872113856834,0.2818382931264455,165033.67505712775,0.926207471201827,0.2580835463981576,0.00029844885990209407,skewnorm +16.08,0.2818382931264455,0.35481338923357547,195428.54228237682,0.6973342196323505,0.347863460819145,1.7417417724671047,skewnorm +16.08,0.35481338923357547,0.44668359215096326,217788.5469960096,0.7187340457834079,0.3133342684773533,1.5519441260373916,skewnorm +16.08,0.44668359215096326,0.5623413251903494,237875.6253430147,0.7328591111067667,0.3003732239985521,1.6110676099162526,skewnorm +16.08,0.5623413251903494,0.707945784384138,256093.01658506907,0.7483573597994009,0.29017107609582027,1.6432184527683333,skewnorm +16.08,0.707945784384138,0.8912509381337459,284598.23127970134,0.7633409576060431,0.28309786128416703,1.3693550676110542,skewnorm +16.08,0.8912509381337459,1.1220184543019636,304509.8328308484,0.7798433193712037,0.2695146883293678,1.166585452834914,skewnorm +16.08,1.1220184543019636,1.4125375446227548,332609.59416784986,1.0838396785530766,0.2505567724134734,-0.9492569083239024,skewnorm +16.08,1.4125375446227548,1.7782794100389236,325433.4288091243,1.1359861569878325,0.2767683348755537,-1.4852215404691282,skewnorm +16.08,1.7782794100389236,2.238721138568341,360325.8057380744,1.1337800836873286,0.3067413170992409,-1.6755889333509033,skewnorm +16.08,2.238721138568341,2.8183829312644564,338337.830413246,1.1373846892256854,0.2863539685699161,-1.4599366873525614,skewnorm +16.08,2.8183829312644564,3.548133892335755,364417.1235649754,1.1459588473382347,0.31211048233110056,-1.585479457926949,skewnorm +16.08,3.548133892335755,4.466835921509633,380127.0098033579,1.1506968141411085,0.3247763893179517,-1.999999999999997,skewnorm +16.08,4.466835921509633,5.623413251903494,351037.8158309781,1.160167252958713,0.3250130362597348,-1.9501529455301712,skewnorm +16.08,5.623413251903494,7.079457843841384,372052.5388440894,1.1561306270806064,0.34429367408305844,-1.9999999999999998,skewnorm +16.08,7.079457843841384,8.912509381337463,371697.49504366505,1.1398799317138533,0.38294212201225686,-1.999999999999863,skewnorm +16.08,8.912509381337463,11.220184543019636,432020.83098216756,1.1869139101524524,0.36621543877976237,-1.999999999999972,skewnorm +16.08,11.220184543019636,14.125375446227547,393152.5792089769,1.116693118482607,0.35167963003235153,-1.9999999999999998,skewnorm +16.08,14.125375446227547,17.782794100389236,401652.5311832727,1.1317963437029088,0.35022154064552385,-1.644784838580502,skewnorm +16.08,17.782794100389236,22.38721138568341,412085.31374464853,1.1185648657292417,0.37182420317742093,-1.9999999999999982,skewnorm +16.08,22.38721138568341,28.183829312644566,451117.1072176613,1.0726839583078185,0.37797359796024377,-1.9999999999999998,skewnorm +16.08,28.183829312644566,35.48133892335755,422786.96586759074,0.9943611608751742,0.36643197174397063,-1.9999999999999998,skewnorm +23.16,0.011220184543019636,0.014125375446227547,908.5016660928053,2.5794239284942555,0.32222018166297006,,moyal +23.16,0.014125375446227547,0.01778279410038923,2235.0222438582546,2.079305052899657,0.26540502874244004,,moyal +23.16,0.01778279410038923,0.0223872113856834,4628.483911517809,1.7009666518162136,0.23218396710516853,,moyal +23.16,0.0223872113856834,0.02818382931264455,8087.426920892466,1.4187390993625388,0.20287589276577347,,moyal +23.16,0.02818382931264455,0.03548133892335755,13688.942230314951,1.2000347407428265,0.1939175274359521,,moyal +23.16,0.03548133892335755,0.04466835921509633,21139.018220354254,1.0321471851225001,0.18111924342549865,,moyal +23.16,0.04466835921509633,0.05623413251903491,31021.103532342346,0.9200666697286745,0.18808430713712096,,moyal +23.16,0.05623413251903491,0.0707945784384138,42380.662961348724,0.8357616979656963,0.1890222030433546,,moyal +23.16,0.0707945784384138,0.08912509381337459,57763.52279144689,0.7849706679031384,0.19454886668353816,,moyal +23.16,0.08912509381337459,0.11220184543019636,74855.02152362485,0.6496770786571662,0.4545978604161516,1.9999999999999998,skewnorm +23.16,0.11220184543019636,0.14125375446227548,95527.06867733231,0.6185690482332782,0.4758933455769919,1.9999999999999853,skewnorm +23.16,0.14125375446227548,0.17782794100389238,122342.5120767525,0.6384925570467842,0.43916966640576804,1.9999999999999996,skewnorm +23.16,0.17782794100389238,0.223872113856834,145777.00823942907,0.6634808198235156,0.4142780935140345,1.7703367775043373,skewnorm +23.16,0.223872113856834,0.2818382931264455,174673.40504871876,0.6921949373671988,0.37109955596606153,1.6596542342986622,skewnorm +23.16,0.2818382931264455,0.35481338923357547,203595.5871542035,0.9278981066663361,0.24826100749623561,0.000359424248897094,skewnorm +23.16,0.35481338923357547,0.44668359215096326,227794.80832549004,0.7146998146362111,0.3319574398905043,1.7228352100003321,skewnorm +23.16,0.44668359215096326,0.5623413251903494,252641.68830583,0.7316435789248908,0.3273579883557315,1.6093229022547375,skewnorm +23.16,0.5623413251903494,0.707945784384138,281751.928110811,0.7490141638800676,0.30456585962184685,1.7038607072633436,skewnorm +23.16,0.707945784384138,0.8912509381337459,308914.4355978287,0.7521214020155008,0.2846461008362118,1.6260742528949022,skewnorm +23.16,0.8912509381337459,1.1220184543019636,325168.8470635501,0.7754203465925146,0.27890673228259927,1.3801300945101982,skewnorm +23.16,1.1220184543019636,1.4125375446227548,360444.00421019294,0.7865801502267257,0.2695007693664174,1.4412018057893456,skewnorm +23.16,1.4125375446227548,1.7782794100389236,364777.28673417034,1.0696125096216766,0.24594136266041428,-0.8494103669927177,skewnorm +23.16,1.7782794100389236,2.238721138568341,382334.49706279155,0.791963564171316,0.2597933249219881,1.3871353419188654,skewnorm +23.16,2.238721138568341,2.8183829312644564,409466.2104841779,1.0665918299339825,0.23017691976151075,-0.784025708336966,skewnorm +23.16,2.8183829312644564,3.548133892335755,402397.51266425336,1.1653706673892983,0.3013702037140571,-1.7202639008729896,skewnorm +23.16,3.548133892335755,4.466835921509633,438165.0137753362,1.1564474645599545,0.27980620773139525,-1.67085475655596,skewnorm +23.16,4.466835921509633,5.623413251903494,425349.72935860837,1.1620303060536372,0.29250155513473064,-1.7726123885916696,skewnorm +23.16,5.623413251903494,7.079457843841384,479577.7499758475,1.1710766570312587,0.33260499370440166,-1.9112463927843453,skewnorm +23.16,7.079457843841384,8.912509381337463,492668.9186998042,1.1692324651110242,0.3437179773931218,-1.9373924810409893,skewnorm +23.16,8.912509381337463,11.220184543019636,433416.12576428615,1.1338947375394892,0.34407427121962914,-1.9999999999999996,skewnorm +23.16,11.220184543019636,14.125375446227547,446187.2087911416,1.158177704592475,0.3313134848286514,-1.9999999999999998,skewnorm +23.16,14.125375446227547,17.782794100389236,407751.97635869676,1.1335835068639246,0.3624019077755211,-1.9999999999999998,skewnorm +23.16,17.782794100389236,22.38721138568341,451559.406001737,1.1192137054604614,0.36998905524295017,-1.9999999999999998,skewnorm +23.16,22.38721138568341,28.183829312644566,493397.3983568351,1.1218241434708074,0.3242496065641991,-1.9999999999999998,skewnorm +23.16,28.183829312644566,35.48133892335755,526627.6172451641,1.1095945214765408,0.3373592075027871,-1.9999999999999998,skewnorm +30.39,0.011220184543019636,0.014125375446227547,502.77408931643265,3.0697172480955075,0.33524280250684096,,moyal +30.39,0.014125375446227547,0.01778279410038923,1244.811652283187,2.441648407919946,0.2780653883734285,,moyal +30.39,0.01778279410038923,0.0223872113856834,2766.184360790552,1.9917107062982558,0.2550804605976299,,moyal +30.39,0.0223872113856834,0.02818382931264455,5863.745249258751,1.6278117089554691,0.2221814886997794,,moyal +30.39,0.02818382931264455,0.03548133892335755,10761.809398505451,1.3711753773209958,0.20063373234983026,,moyal +30.39,0.03548133892335755,0.04466835921509633,17697.226644734077,1.152911306354468,0.1846088840981204,,moyal +30.39,0.04466835921509633,0.05623413251903491,26978.43497332758,1.0048828976595559,0.1820304343605322,,moyal +30.39,0.05623413251903491,0.0707945784384138,38919.521805983335,0.8910413235418876,0.18224602934721473,,moyal +30.39,0.0707945784384138,0.08912509381337459,54480.30494735733,0.8231874192231426,0.1845300009847712,,moyal +30.39,0.08912509381337459,0.11220184543019636,74168.27875448299,0.7826087726678402,0.18218246073561442,,moyal +30.39,0.11220184543019636,0.14125375446227548,98555.69231767126,0.6313697788706518,0.45820176228301696,1.9999999999999998,skewnorm +30.39,0.14125375446227548,0.17782794100389238,125581.30143433146,0.9412271667198819,0.2620600163340633,0.00013434274457144662,skewnorm +30.39,0.17782794100389238,0.223872113856834,156408.94286939196,0.64914359525058,0.4110167222216276,1.91331361798301,skewnorm +30.39,0.223872113856834,0.2818382931264455,184510.9248493845,0.9230957316674361,0.2673073236585434,-0.00012925555950155896,skewnorm +30.39,0.2818382931264455,0.35481338923357547,214394.2470280635,0.6784455949678745,0.3655438585017867,1.7925553771146192,skewnorm +30.39,0.35481338923357547,0.44668359215096326,259571.94246618496,0.7104442180534197,0.33275161901535866,1.6250148132837055,skewnorm +30.39,0.44668359215096326,0.5623413251903494,289418.8934693836,0.7218556065509791,0.3214974849367867,1.500886369297606,skewnorm +30.39,0.5623413251903494,0.707945784384138,323439.4050501878,0.7271810738507521,0.31561601108672077,1.748594214135925,skewnorm +30.39,0.707945784384138,0.8912509381337459,346761.5139516393,0.7320633144868856,0.30668068407396837,1.9057448464571958,skewnorm +30.39,0.8912509381337459,1.1220184543019636,368749.4679613911,0.750777603990807,0.2866423202139152,1.7112308248740222,skewnorm +30.39,1.1220184543019636,1.4125375446227548,406409.0555307833,0.7796034908535391,0.26402190450009744,1.399286962207778,skewnorm +30.39,1.4125375446227548,1.7782794100389236,436437.5962484389,0.7690973831932351,0.2660123414800899,1.4797662742056608,skewnorm +30.39,1.7782794100389236,2.238721138568341,427828.9254857811,0.7588190699526662,0.28963888671194044,1.726794880556474,skewnorm +30.39,2.238721138568341,2.8183829312644564,497522.6620857837,1.0941747458178512,0.25080191216266406,-1.318744137812615,skewnorm +30.39,2.8183829312644564,3.548133892335755,495241.68983902026,1.122081139654597,0.2745560952797514,-1.3689247590756566,skewnorm +30.39,3.548133892335755,4.466835921509633,441389.6527863084,1.0999005224626839,0.252019093098836,-1.1406287094739345,skewnorm +30.39,4.466835921509633,5.623413251903494,477043.1621321472,1.1260952660085313,0.3125320395119026,-1.915429084679942,skewnorm +30.39,5.623413251903494,7.079457843841384,526166.4964366799,1.1424054712218517,0.3023336997733039,-1.926319502856402,skewnorm +30.39,7.079457843841384,8.912509381337463,502279.2567585283,1.1673310479974661,0.354064092159078,-1.9999999999999991,skewnorm +30.39,8.912509381337463,11.220184543019636,555311.215289976,1.1304090114068157,0.3039393334709209,-1.7314595107944246,skewnorm +30.39,11.220184543019636,14.125375446227547,511153.65747170715,1.1594839981461003,0.3565710858581665,-1.879867049579473,skewnorm +30.39,14.125375446227547,17.782794100389236,598293.2517632255,1.1493473893245407,0.37805236720317953,-1.9999999999999998,skewnorm +30.39,17.782794100389236,22.38721138568341,478105.9389020703,1.167253937202583,0.36193145027239937,-1.9999999999999998,skewnorm +30.39,22.38721138568341,28.183829312644566,573237.8247747866,1.1319415234006525,0.29961281511905913,-1.2558828227990926,skewnorm +30.39,28.183829312644566,35.48133892335755,550268.5442655784,1.160786394204726,0.4017170402608415,-1.9999999999999998,skewnorm +37.66,0.011220184543019636,0.014125375446227547,171.4543233249978,3.6784329710301598,0.5285902950595364,,moyal +37.66,0.014125375446227547,0.01778279410038923,471.28986683938007,3.008993033541901,0.3030062493652152,,moyal +37.66,0.01778279410038923,0.0223872113856834,1291.1517293149132,2.464345755291737,0.27276437882744264,,moyal +37.66,0.0223872113856834,0.02818382931264455,3078.8671506197215,2.0131009816628627,0.2426646544124947,,moyal +37.66,0.02818382931264455,0.03548133892335755,6623.202224907022,1.6378591389717494,0.20920412535558552,,moyal +37.66,0.03548133892335755,0.04466835921509633,12290.622732247406,1.3808229194927255,0.19513544801224963,,moyal +37.66,0.04466835921509633,0.05623413251903491,20860.852955868686,1.1706854947542487,0.18917825111414044,,moyal +37.66,0.05623413251903491,0.0707945784384138,32826.55120783815,1.015556302322552,0.18154245527186544,,moyal +37.66,0.0707945784384138,0.08912509381337459,48284.289334990775,0.9038845658989394,0.17982822612375354,,moyal +37.66,0.08912509381337459,0.11220184543019636,69554.11876684925,0.8317599030439307,0.1796733490614012,,moyal +37.66,0.11220184543019636,0.14125375446227548,94330.3135074456,0.6535737593189439,0.4811118956150625,1.9999999999999998,skewnorm +37.66,0.14125375446227548,0.17782794100389238,123020.70453650136,0.6509329729676441,0.45839996557007834,1.9999999999999751,skewnorm +37.66,0.17782794100389238,0.223872113856834,157864.60371623706,0.6471289085971138,0.4431603870709839,1.926827282799218,skewnorm +37.66,0.223872113856834,0.2818382931264455,197241.8245604115,0.6668255857841946,0.40533360550104064,1.9672169275501672,skewnorm +37.66,0.2818382931264455,0.35481338923357547,238323.15200973648,0.6936224634718618,0.3678478568035432,1.7122920411503608,skewnorm +37.66,0.35481338923357547,0.44668359215096326,284503.3735044894,0.7088501019196399,0.34856649115778,1.6973942970503613,skewnorm +37.66,0.44668359215096326,0.5623413251903494,328967.37139506306,0.732360997056354,0.320675137547888,1.5722931076191238,skewnorm +37.66,0.5623413251903494,0.707945784384138,357095.5222005081,0.7321490713907581,0.3130900973918259,1.6981836972462145,skewnorm +37.66,0.707945784384138,0.8912509381337459,414734.072121709,0.7479697485923943,0.2933644377914268,1.5490692744799286,skewnorm +37.66,0.8912509381337459,1.1220184543019636,438873.51500899217,0.770833630678645,0.2785176970589174,1.4600204233587808,skewnorm +37.66,1.1220184543019636,1.4125375446227548,480275.009834639,0.7704412742398293,0.2766180081492073,1.5049345626033075,skewnorm +37.66,1.4125375446227548,1.7782794100389236,517013.9243809778,0.7687057843633154,0.2710802753644303,1.559272736296874,skewnorm +37.66,1.7782794100389236,2.238721138568341,555854.4867050623,0.7751025673991642,0.27110896154819003,1.615805108694842,skewnorm +37.66,2.238721138568341,2.8183829312644564,602211.0304466586,0.811413990471518,0.2532925792839454,1.2235352891860358,skewnorm +37.66,2.8183829312644564,3.548133892335755,542037.365616671,1.0598231411458043,0.2200298125592479,-0.7524516536443839,skewnorm +37.66,3.548133892335755,4.466835921509633,622070.4118461943,1.1193376152098062,0.2552294909844531,-1.2892656537665688,skewnorm +37.66,4.466835921509633,5.623413251903494,628450.8654275239,1.1096897095543634,0.2387050179128118,-1.1598057122793364,skewnorm +37.66,5.623413251903494,7.079457843841384,658971.659934298,1.1478107328585592,0.300147251735972,-1.7459327091636392,skewnorm +37.66,7.079457843841384,8.912509381337463,688129.5377849016,1.1156725077264424,0.28015071674317255,-1.1464349054091538,skewnorm +37.66,8.912509381337463,11.220184543019636,650332.2192659881,1.1714157698994137,0.3404993947428983,-1.5932562941175814,skewnorm +37.66,11.220184543019636,14.125375446227547,665526.840835495,1.1678430042963155,0.3475811958820149,-1.9999999999999998,skewnorm +37.66,14.125375446227547,17.782794100389236,749721.2786848663,1.177241778122588,0.31643503800821854,-1.9999999999999643,skewnorm +37.66,17.782794100389236,22.38721138568341,689793.389633994,1.157003956252038,0.34111896372216954,-1.4804638734397546,skewnorm +37.66,22.38721138568341,28.183829312644566,667998.790075285,1.2006308673333859,0.41758170963497376,-1.9842402317025212,skewnorm +37.66,28.183829312644566,35.48133892335755,621801.2092640705,1.1481134632339047,0.32718679470557666,-1.9999999999999998,skewnorm +44.92,0.011220184543019636,0.014125375446227547,23.67207094263638,,,,-- +44.92,0.014125375446227547,0.01778279410038923,93.34292236777446,4.017816179546889,0.4909974090314828,,moyal +44.92,0.01778279410038923,0.0223872113856834,367.3931418222828,2.921743519396353,0.8907845246608525,1.9999999999999998,skewnorm +44.92,0.0223872113856834,0.02818382931264455,1050.8739247714825,2.6742582661572674,0.30455149588631275,,moyal +44.92,0.02818382931264455,0.03548133892335755,2658.788164240249,2.138062394373709,0.2538280996799505,,moyal +44.92,0.03548133892335755,0.04466835921509633,6036.846864226026,1.7749775029485855,0.22447367596417836,,moyal +44.92,0.04466835921509633,0.05623413251903491,11973.298211385872,1.4640002396402942,0.2005862299434763,,moyal +44.92,0.05623413251903491,0.0707945784384138,22381.718435259063,1.2322992573580234,0.18063125110463304,,moyal +44.92,0.0707945784384138,0.08912509381337459,37275.261169034166,1.0678756181823104,0.17933622679147646,,moyal +44.92,0.08912509381337459,0.11220184543019636,57169.15786940863,0.9416176444453241,0.17611172515891357,,moyal +44.92,0.11220184543019636,0.14125375446227548,82590.62156215431,0.8566215643867953,0.17819846540647902,,moyal +44.92,0.14125375446227548,0.17782794100389238,114050.03043326904,0.8132069541323005,0.17262512870018665,,moyal +44.92,0.17782794100389238,0.223872113856834,156493.26989532993,0.6690047636588552,0.4367123370359943,1.9999999999999998,skewnorm +44.92,0.223872113856834,0.2818382931264455,201485.51997762633,0.6664611424378484,0.4230472051482394,1.8539645961343114,skewnorm +44.92,0.2818382931264455,0.35481338923357547,259413.80816637157,0.6785279983665844,0.39184648740498573,1.907893630748769,skewnorm +44.92,0.35481338923357547,0.44668359215096326,319467.78003810387,0.6974142198423158,0.3592499899218503,1.6486302111725837,skewnorm +44.92,0.44668359215096326,0.5623413251903494,364423.90503965813,0.709183006831286,0.3413232764869068,1.6476136009269755,skewnorm +44.92,0.5623413251903494,0.707945784384138,418546.67737996986,0.7321647627279166,0.31902983509174143,1.513438692362349,skewnorm +44.92,0.707945784384138,0.8912509381337459,476802.8574521433,0.745235068526115,0.30164660356476536,1.5061136426091806,skewnorm +44.92,0.8912509381337459,1.1220184543019636,549419.6391570403,0.7513171274470741,0.288179904253016,1.4593762626174334,skewnorm +44.92,1.1220184543019636,1.4125375446227548,592221.885497276,0.7612968986708633,0.2783199856393788,1.5462084217665124,skewnorm +44.92,1.4125375446227548,1.7782794100389236,655117.1080202644,0.9414514440015271,0.19921359634336486,-0.00013974416874566221,skewnorm +44.92,1.7782794100389236,2.238721138568341,680082.525395291,0.7739717752643646,0.26391164705334125,1.6235236227572274,skewnorm +44.92,2.238721138568341,2.8183829312644564,741223.2100451286,0.7665813977280294,0.283593205923416,1.922584809190939,skewnorm +44.92,2.8183829312644564,3.548133892335755,766809.9198316105,0.7801384877673504,0.27438636119461757,1.8018616298192072,skewnorm +44.92,3.548133892335755,4.466835921509633,762690.9207928411,0.7901731209060258,0.25529852972903466,1.8708706028795545,skewnorm +44.92,4.466835921509633,5.623413251903494,816654.2588496949,0.7899001281400524,0.2713019333983877,1.8608652533277088,skewnorm +44.92,5.623413251903494,7.079457843841384,905947.0381905444,0.8921599302914988,0.10281490794312928,,moyal +44.92,7.079457843841384,8.912509381337463,925966.5281747963,0.8842967558860888,0.1163973327057996,,moyal +44.92,8.912509381337463,11.220184543019636,922694.6379633403,0.900303143725539,0.09909406145438183,,moyal +44.92,11.220184543019636,14.125375446227547,863319.571173897,0.8958194349751416,0.11311232882600297,,moyal +44.92,14.125375446227547,17.782794100389236,850815.3728464945,0.764979027141839,0.3205668740171627,1.849116444384866,skewnorm +44.92,17.782794100389236,22.38721138568341,854999.4807849439,0.8732409005370527,0.11316165047782017,,moyal +44.92,22.38721138568341,28.183829312644566,1055975.2545033835,1.1055834601524641,0.2642319441530294,-1.411551509800831,skewnorm +44.92,28.183829312644566,35.48133892335755,986877.0883371652,1.1458306223070578,0.3657048574108705,-1.6183866938253704,skewnorm +52.16,0.011220184543019636,0.014125375446227547,0.0,,,,-- +52.16,0.014125375446227547,0.01778279410038923,4.907019979208879,,,,-- +52.16,0.01778279410038923,0.0223872113856834,36.25375661776856,,,,-- +52.16,0.0223872113856834,0.02818382931264455,149.36981065890456,3.8247960166623334,0.4605133984927266,,moyal +52.16,0.02818382931264455,0.03548133892335755,548.4658969806877,3.0794108117620467,0.3048609309010585,,moyal +52.16,0.03548133892335755,0.04466835921509633,1394.1072657602986,2.497847579423691,0.2649158989520182,,moyal +52.16,0.04466835921509633,0.05623413251903491,3692.284765943902,2.042288477979181,0.2211344406438066,,moyal +52.16,0.05623413251903491,0.0707945784384138,8629.600015438922,1.6753718404531721,0.20407426468316253,,moyal +52.16,0.0707945784384138,0.08912509381337459,18740.876609034774,1.3968321580857777,0.18920883252302167,,moyal +52.16,0.08912509381337459,0.11220184543019636,35067.94557656182,1.1789275585501182,0.1776594255014463,,moyal +52.16,0.11220184543019636,0.14125375446227548,59078.78612951074,1.0249419343609414,0.17225064585948976,,moyal +52.16,0.14125375446227548,0.17782794100389238,92256.34174274263,0.9141768024759521,0.1701615803641573,,moyal +52.16,0.17782794100389238,0.223872113856834,135918.62935607068,0.8443420411773828,0.17063443031659745,,moyal +52.16,0.223872113856834,0.2818382931264455,188482.57219069943,0.6956858024611134,0.44000810181464467,1.9999999999999998,skewnorm +52.16,0.2818382931264455,0.35481338923357547,257813.7948124777,0.682238956099989,0.4194782007961768,1.9999999999999998,skewnorm +52.16,0.35481338923357547,0.44668359215096326,321762.5873285441,0.6923675676186579,0.399906607940871,1.7786979557859035,skewnorm +52.16,0.44668359215096326,0.5623413251903494,407034.38189590763,0.7062217722297798,0.3615181089950087,1.729527797110049,skewnorm +52.16,0.5623413251903494,0.707945784384138,495680.88977569615,0.7249324178514056,0.33909487660048643,1.6624237522433485,skewnorm +52.16,0.707945784384138,0.8912509381337459,565244.0823691139,0.7317636750744312,0.32607550007992664,1.6322733272178864,skewnorm +52.16,0.8912509381337459,1.1220184543019636,656154.3530159951,0.7650637084734404,0.28852362171045803,1.3857866402951273,skewnorm +52.16,1.1220184543019636,1.4125375446227548,749384.9077157578,0.7443420810993215,0.2975891391842063,1.6412970483321656,skewnorm +52.16,1.4125375446227548,1.7782794100389236,815140.4542366246,0.7863566669178169,0.2548367988844022,1.2309581035797028,skewnorm +52.16,1.7782794100389236,2.238721138568341,887118.4384499419,0.7852451355600747,0.2509263778449068,1.2530959216473558,skewnorm +52.16,2.238721138568341,2.8183829312644564,959977.6442439421,1.0271705325000962,0.21040982020792684,-0.6647724396733149,skewnorm +52.16,2.8183829312644564,3.548133892335755,1029722.957061523,0.7980418721951797,0.23953464782003672,1.2722331172356227,skewnorm +52.16,3.548133892335755,4.466835921509633,1085255.5145879297,1.0190266809042545,0.20715079012601423,-0.6563373912264494,skewnorm +52.16,4.466835921509633,5.623413251903494,1122587.0274932557,1.0525148217361422,0.2099769930990669,-0.9672122159912855,skewnorm +52.16,5.623413251903494,7.079457843841384,1206545.9280845176,1.0527558631058198,0.2078103471071983,-1.0806567918328849,skewnorm +52.16,7.079457843841384,8.912509381337463,1258284.7608913204,1.070010121814587,0.23040309754846047,-1.2007294548153287,skewnorm +52.16,8.912509381337463,11.220184543019636,1305481.2209456463,1.1103397665086414,0.2707164095768375,-1.6427407508797423,skewnorm +52.16,11.220184543019636,14.125375446227547,1320486.771990121,1.1163630608423438,0.2737190189817738,-1.7760238498641545,skewnorm +52.16,14.125375446227547,17.782794100389236,1288028.3807443073,1.1185581337755597,0.28611659001235923,-1.999999999916018,skewnorm +52.16,17.782794100389236,22.38721138568341,1304036.6395542796,1.1310891463378188,0.2980428193571546,-1.959115178349524,skewnorm +52.16,22.38721138568341,28.183829312644566,1363691.049163704,1.1249016629287278,0.2967848019933458,-1.8454010488169759,skewnorm +52.16,28.183829312644566,35.48133892335755,1328506.2837976706,1.1327489243568662,0.33699451376709555,-1.999999999999999,skewnorm +59.34,0.011220184543019636,0.014125375446227547,0.0,,,,-- +59.34,0.014125375446227547,0.01778279410038923,0.0,,,,-- +59.34,0.01778279410038923,0.0223872113856834,0.0,,,,-- +59.34,0.0223872113856834,0.02818382931264455,0.0,,,,-- +59.34,0.02818382931264455,0.03548133892335755,18.464545725394675,,,,-- +59.34,0.03548133892335755,0.04466835921509633,89.66115963360522,4.321315182278122,0.6002907329201759,,moyal +59.34,0.04466835921509633,0.05623413251903491,321.9076610427148,3.316100127353422,0.355869233385713,,moyal +59.34,0.05623413251903491,0.0707945784384138,1176.300048232724,2.6210711941664857,0.2818400902571692,,moyal +59.34,0.0707945784384138,0.08912509381337459,3339.4206147557043,2.1434389142287666,0.24042311116027476,,moyal +59.34,0.08912509381337459,0.11220184543019636,9000.404581189492,1.7511933486293878,0.2050014856943333,,moyal +59.34,0.11220184543019636,0.14125375446227548,20634.9367524069,1.450116613151804,0.18650608530514728,,moyal +59.34,0.14125375446227548,0.17782794100389238,42939.462908247944,1.2306994515096414,0.17466246797792415,,moyal +59.34,0.17782794100389238,0.223872113856834,79688.33064742846,1.0576248379502982,0.1695987073801281,,moyal +59.34,0.223872113856834,0.2818382931264455,133762.21930067122,0.942126615810544,0.16834467279855128,,moyal +59.34,0.2818382931264455,0.35481338923357547,203017.68025071456,0.8697672625624766,0.1695035875376232,,moyal +59.34,0.35481338923357547,0.44668359215096326,293092.3684911849,0.7237158526138695,0.4169383405696297,1.9999999999999998,skewnorm +59.34,0.44668359215096326,0.5623413251903494,386205.5808717653,0.7207849165243896,0.387883937201047,1.9999999999999998,skewnorm +59.34,0.5623413251903494,0.707945784384138,504519.5643113614,0.7030572437591043,0.399517228704227,1.867318267282071,skewnorm +59.34,0.707945784384138,0.8912509381337459,632005.2288461762,0.7237963972829706,0.3574024667823403,1.7100102024545447,skewnorm +59.34,0.8912509381337459,1.1220184543019636,748768.5979893187,0.7374512294340662,0.33724637661063955,1.592482869852303,skewnorm +59.34,1.1220184543019636,1.4125375446227548,893550.5181996451,0.7542483993071667,0.30947788864917575,1.496262923161751,skewnorm +59.34,1.4125375446227548,1.7782794100389236,1019019.0273915499,0.7597624481878245,0.29089591113817975,1.3948313679540356,skewnorm +59.34,1.7782794100389236,2.238721138568341,1160457.1542171962,0.7602929435628104,0.2903581475652521,1.4881355852231122,skewnorm +59.34,2.238721138568341,2.8183829312644564,1301897.7907655397,0.7836578917894391,0.2605697221027763,1.2511832319833105,skewnorm +59.34,2.8183829312644564,3.548133892335755,1427441.2741853371,0.7798888116034636,0.25726637325872287,1.3209038858121713,skewnorm +59.34,3.548133892335755,4.466835921509633,1511454.6965642737,0.7760337446424336,0.2573717188248899,1.3781487598251643,skewnorm +59.34,4.466835921509633,5.623413251903494,1560625.0632109765,0.9334088503069673,0.18279092210724818,1.1265992095845703e-05,skewnorm +59.34,5.623413251903494,7.079457843841384,1772871.0122916938,0.7879833199189711,0.24798842746132585,1.3327476568137475,skewnorm +59.34,7.079457843841384,8.912509381337463,1835024.8794773673,0.7843668388733914,0.2459463391913026,1.5085641181221547,skewnorm +59.34,8.912509381337463,11.220184543019636,1861807.5092877646,0.9407968219222826,0.18153000016567852,0.0005811142877605452,skewnorm +59.34,11.220184543019636,14.125375446227547,1952705.5924224248,0.8112200417020734,0.23593852344806548,1.187150758411722,skewnorm +59.34,14.125375446227547,17.782794100389236,2066329.4496795428,0.7980265211721244,0.252757725362054,1.4363395728233581,skewnorm +59.34,17.782794100389236,22.38721138568341,1988047.4302078749,1.059719683679896,0.21914485389911806,-0.7344085299489629,skewnorm +59.34,22.38721138568341,28.183829312644566,2004128.4893654762,1.0946242475310164,0.24485456995943652,-1.2296599781465691,skewnorm +59.34,28.183829312644566,35.48133892335755,2162989.6421176693,1.082807554808789,0.22885200346161932,-1.0894197584813705,skewnorm +66.44,0.011220184543019636,0.014125375446227547,0.0,,,,-- +66.44,0.014125375446227547,0.01778279410038923,0.0,,,,-- +66.44,0.01778279410038923,0.0223872113856834,0.0,,,,-- +66.44,0.0223872113856834,0.02818382931264455,0.0,,,,-- +66.44,0.02818382931264455,0.03548133892335755,0.0,,,,-- +66.44,0.03548133892335755,0.04466835921509633,0.0,,,,-- +66.44,0.04466835921509633,0.05623413251903491,3.0284276680418656,,,,-- +66.44,0.05623413251903491,0.0707945784384138,14.687162502501202,,,,-- +66.44,0.0707945784384138,0.08912509381337459,83.2051894569713,,,,-- +66.44,0.08912509381337459,0.11220184543019636,293.8794963181524,3.262364971790911,0.33318544331092664,,moyal +66.44,0.11220184543019636,0.14125375446227548,1062.2968922460948,2.664660230636009,0.2820942219906546,,moyal +66.44,0.14125375446227548,0.17782794100389238,3592.4056496964577,2.165591055675511,0.23903780313715925,,moyal +66.44,0.17782794100389238,0.223872113856834,10955.187455275569,1.7769203178492858,0.2048872238156751,,moyal +66.44,0.223872113856834,0.2818382931264455,28650.617065917886,1.482479553350125,0.18572856081387176,,moyal +66.44,0.2818382931264455,0.35481338923357547,65908.71800173383,1.2435467416198969,0.1766719108954862,,moyal +66.44,0.35481338923357547,0.44668359215096326,131324.59184228396,1.0720130408838837,0.16610307696329182,,moyal +66.44,0.44668359215096326,0.5623413251903494,227174.51657629633,0.958492320682369,0.16513300434776837,,moyal +66.44,0.5623413251903494,0.707945784384138,353171.1813257706,0.8833967890112088,0.16051224030041294,,moyal +66.44,0.707945784384138,0.8912509381337459,501796.63008341746,0.7418824015199696,0.41584229454395766,1.9999999999999998,skewnorm +66.44,0.8912509381337459,1.1220184543019636,690354.943529157,0.7371183051530632,0.39143177792235173,1.9999999999999998,skewnorm +66.44,1.1220184543019636,1.4125375446227548,892622.4369052699,0.7172564547703438,0.38820319161345745,1.9098952640891538,skewnorm +66.44,1.4125375446227548,1.7782794100389236,1107051.2198538284,0.7336739276189062,0.3574851419220619,1.72093366489278,skewnorm +66.44,1.7782794100389236,2.238721138568341,1336683.815422703,0.749814876498203,0.3337768093576935,1.5040192803690027,skewnorm +66.44,2.238721138568341,2.8183829312644564,1530030.657290625,0.7628820719180656,0.3125330679949241,1.3910902925840725,skewnorm +66.44,2.8183829312644564,3.548133892335755,1782378.579382366,0.7649252814948652,0.29614768805288944,1.3626240715430007,skewnorm +66.44,3.548133892335755,4.466835921509633,2053445.390833702,0.7674819217005542,0.28597747938007895,1.3365577850373231,skewnorm +66.44,4.466835921509633,5.623413251903494,2215746.312017105,0.9404642800086161,0.20727430969863225,0.00207322381085733,skewnorm +66.44,5.623413251903494,7.079457843841384,2505997.1019892707,0.7830428450863928,0.25946711019071744,1.1858109196040718,skewnorm +66.44,7.079457843841384,8.912509381337463,2695848.138405872,0.7734218949239624,0.25835036123201227,1.2844788007940127,skewnorm +66.44,8.912509381337463,11.220184543019636,2963818.365838322,0.77396157946228,0.2638255289741927,1.291773304711111,skewnorm +66.44,11.220184543019636,14.125375446227547,3181762.3469136492,0.7654908995131571,0.2641988949488634,1.4306831398400233,skewnorm +66.44,14.125375446227547,17.782794100389236,3300493.868405334,0.9310273637947536,0.19216863271465254,-0.00017125270443750374,skewnorm +66.44,17.782794100389236,22.38721138568341,3621539.9759411146,0.9341296501604117,0.1892115433229644,0.00033485180864913436,skewnorm +66.44,22.38721138568341,28.183829312644566,3483388.7993919547,0.7722802929786422,0.26542310146783166,1.4931862700122722,skewnorm +66.44,28.183829312644566,35.48133892335755,3691550.6997481054,0.7913556135004776,0.2510805767288358,1.2058824623523396,skewnorm diff --git a/notebooks/LST1_observation_simulator.ipynb b/notebooks/LST1_observation_simulator.ipynb new file mode 100644 index 0000000000..b073fa5c52 --- /dev/null +++ b/notebooks/LST1_observation_simulator.ipynb @@ -0,0 +1,984 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": null, + "id": "8d7fb8ab-7034-47be-a392-17612fd611b7", + "metadata": {}, + "outputs": [], + "source": [ + "import os\n", + "import pandas as pd\n", + "import numpy as np\n", + "import astropy.units as u\n", + "import matplotlib.pyplot as plt\n", + "import gammapy\n", + "from gammapy.modeling.models import EBLAbsorptionNormSpectralModel\n", + "import subprocess\n", + "from pyirf.spectral import CRAB_MAGIC_JHEAP2015, PowerLaw, LogParabola\n", + "from pyirf.statistics import li_ma_significance\n", + "from scipy.stats import moyal, norm, skewnorm\n", + "from pathlib import Path" + ] + }, + { + "cell_type": "markdown", + "id": "c30c3f8d", + "metadata": {}, + "source": [ + "## Set cut efficiency" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "ef3576a4", + "metadata": {}, + "outputs": [], + "source": [ + "#\n", + "# Cut efficiency: the same for the gammaness cut and the theta cut\n", + "# Options: 0.4, 0.7, 0.9\n", + "#\n", + "# Recommended: \n", + "# 0.7: standard cuts (safer for spectral analysis)\n", + "# 0.4: tight cuts, better for detection of weak sources\n", + "#\n", + "cut_efficiency = 0.7" + ] + }, + { + "cell_type": "markdown", + "id": "5dd2e6fd", + "metadata": {}, + "source": [ + "## Load files which contain the (approximate) instrument response function, and show table contents\n", + "They characterize the average performance of LST1 within 1 degree off-axis (computed from diffuse gamma MC and\n", + "real data for the cosmic ray rates)\n", + "The Ereco/Etrue distributions in each Etrue bin is parametrized with moyal or a skewnorm function (whatever fits better) " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "95ac1d89", + "metadata": {}, + "outputs": [], + "source": [ + "#\n", + "# Load files for the requested efficiency:\n", + "#\n", + "datadir = os.environ['LSTCHAIN'] + '/data/'" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "c59f8666", + "metadata": {}, + "outputs": [], + "source": [ + "input_filename = datadir + f'LST1_gamma_irf_gheffi_{cut_efficiency:.2f}_theffi_{cut_efficiency:.2f}.csv'\n", + "gamma_data = pd.read_csv(input_filename)\n", + "\n", + "input_filename = datadir + f'LST1_backg_irf_gheffi_{cut_efficiency:.2f}_theffi_{cut_efficiency:.2f}.csv'\n", + "background_data = pd.read_csv(input_filename)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "ad3a988e", + "metadata": { + "scrolled": false + }, + "outputs": [], + "source": [ + "gamma_data" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "1d2b1efe", + "metadata": { + "scrolled": false + }, + "outputs": [], + "source": [ + "background_data" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "90114613", + "metadata": {}, + "outputs": [], + "source": [ + "# CHECK that we have the same pointing zenith values in both tables:\n", + "assert np.alltrue(np.unique(gamma_data.ZD_deg) == np.unique(background_data.ZD_deg))\n", + "\n", + "# Available zeniths:\n", + "zenith = np.unique(gamma_data.ZD_deg)\n", + "print('Available zeniths:', zenith, '(degrees)')" + ] + }, + { + "cell_type": "markdown", + "id": "dc762a1a", + "metadata": {}, + "source": [ + "## Zenith distance bin selection" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "a5a9386a", + "metadata": {}, + "outputs": [], + "source": [ + "# Choose the bins among those above. Just set the bin number (from 0)\n", + "# Make sure you choose values which make sense for the declination of your source\n", + "\n", + "zd_bin = 3\n", + "\n", + "print('Selected ZD = ', zenith[zd_bin], 'degrees')\n", + "\n", + "# Cuts for tables:\n", + "zd_selection_gamma = abs(gamma_data.ZD_deg - np.unique(gamma_data.ZD_deg)[zd_bin])<0.01\n", + "zd_selection_backg = abs(background_data.ZD_deg - np.unique(background_data.ZD_deg)[zd_bin])<0.01" + ] + }, + { + "cell_type": "markdown", + "id": "8d1c2dcf", + "metadata": {}, + "source": [ + "## ON to OFF exposure" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "0c58505b", + "metadata": {}, + "outputs": [], + "source": [ + "# Ratio of ON to OFF exposure (Li & Ma's \"alpha\")\n", + "# For standard wobble offset (0.4 deg) reasonable values are alpha=0.333 (3 off regions) above 0.2 TeV, \n", + "# and alpha=1 below 0.2 TeV. For testing sensitivity with the standard definition, set alpha=0.2\n", + "\n", + "alpha = 0.333 # 1" + ] + }, + { + "cell_type": "markdown", + "id": "311a8237", + "metadata": {}, + "source": [ + "## Pulsar mode\n", + "(overrides the setting of alpha above)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "a589ca87", + "metadata": {}, + "outputs": [], + "source": [ + "pulsar_mode = False # Set to True to activate it\n", + "\n", + "on_phase_interval = 0.043 # Crab P1: [-0.017, 0.026] \n", + "off_phase_interval = 0.35 # [0.52 - 0.87]\n", + "\n", + "alpha = on_phase_interval / off_phase_interval\n", + "print(f'alpha = {alpha:.4f}')\n", + "\n", + "# The spectrum is interpreted as average flux in full period (i.e. not just in the on-phase)" + ] + }, + { + "cell_type": "markdown", + "id": "aaf2cba7", + "metadata": {}, + "source": [ + "## Observation time" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "7d8d3a2e", + "metadata": {}, + "outputs": [], + "source": [ + "# effective_obs_time = 103 * u.h # LST1 Crab pulsar paper\n", + "\n", + "# effective_obs_time = 2392 * u.s # GRB190114C, 62 - 2454 s, MAGIC\n", + "\n", + "# effective_obs_time = 8 * u.s # BOAT, 240 - 248 s, LHAASO\n", + "\n", + "effective_obs_time = 34 * u.h # Crab nebula. LST1 performance paper\n", + "\n", + "# effective_obs_time = 11.8 * u.h # 1ES 1011+496, February 2014 flare, MAGIC" + ] + }, + { + "cell_type": "markdown", + "id": "16d22a98", + "metadata": {}, + "source": [ + "## Source extension\n", + "For moderately extended sources, say below ~0.5 degrees radius\n", + "Set here the radius of the source within which (in 2d-gaussian approximation) the same fraction of the source is contained as the cut efficiency set above. That is, if you use 0.7 efficiency, set as source radius the angular distance within which 70% of the emission is expected to be contained.\n", + "We will simply increase the background by a factor (source_radius$^2$ + theta_cut$^2$) / (theta_cut$^2$), since we expect to integrate the same fraction of the signal within an angle sqrt(source_radius$^2$ + theta_cut$^2$) " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "4954c1de", + "metadata": { + "scrolled": false + }, + "outputs": [], + "source": [ + "source_radius = 0 # 0.3 * u.deg\n", + "theta_cut = background_data.Theta_cut_deg[zd_selection_backg].to_numpy()*u.deg\n", + "\n", + "if source_radius == 0:\n", + " print(\"Emission from a point-like source will be assumed!\")\n", + " background_increase_factor = np.ones_like(theta_cut.to_value())\n", + "else:\n", + " print(f'Source angular radius within which {cut_efficiency:.1%} of the emission is contained: {source_radius}')\n", + "\n", + "\n", + " background_increase_factor = (theta_cut**2 + \n", + " (source_radius.to_value(u.deg)*np.ones_like(theta_cut))**2) / (theta_cut**2)" + ] + }, + { + "cell_type": "markdown", + "id": "d662546b", + "metadata": {}, + "source": [ + "## Source redshift" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "ea9bca3c", + "metadata": {}, + "outputs": [], + "source": [ + "redshift = 0\n", + "\n", + "\n", + "# redshift = 0.151 # BOAT GRB\n", + "# redshift = 0.212 # 1ES 1011\n", + "# redshift = 0.42 # GRB190114C\n", + "\n", + "# We will apply the Dominguez EBL model to simulate the absorption\n", + "\n", + "\n", + "# Make sure we have the necessary EBL absorption data:\n", + "try:\n", + " os.environ['GAMMAPY_DATA']\n", + "except:\n", + " # WE SET HERE THE GAMMAPY_DATA ENV VARIABLE IN CASE IT IS NOT SET\n", + " gammapy_dir = Path(gammapy.__file__).parent\n", + " gammapy_dir\n", + "\n", + " ebl_file = subprocess.run(['find', str(gammapy_dir), '-name', 'ebl_dominguez11.fits.gz'], \n", + " stdout=subprocess.PIPE).stdout.decode()\n", + " gammapy_data = ebl_file[:ebl_file.find('/ebl/')]\n", + " os.environ['GAMMAPY_DATA'] = gammapy_data\n", + " print('Set GAMMAPY_DATA to', gammapy_data)\n", + "\n", + "\n", + "dominguez = EBLAbsorptionNormSpectralModel.read_builtin(\"dominguez\", redshift=redshift)" + ] + }, + { + "cell_type": "markdown", + "id": "e5766b4e", + "metadata": {}, + "source": [ + "## Source (intrinsic) spectrum " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "e24f72ce", + "metadata": {}, + "outputs": [], + "source": [ + "# Set here the simulated intrinsic spectrum intrinsic_dFdE\n", + "# (it must take as argument an astropy quantity with energy units!)\n", + "\n", + "# Crab Nebula:\n", + "def intrinsic_dFdE(E):\n", + " return CRAB_MAGIC_JHEAP2015(E) # Crab Nebula, MAGIC log-parabola\n", + "\n", + "# Crab pulsar P1, from LST1 paper (smoothly broken power-law):\n", + "# def intrinsic_dFdE(E):\n", + "# return PowerLaw(normalization=1.27e-4 / (u.TeV * u.cm**2 * u.s), \n", + "# index=-1.811, e_ref=1*u.GeV)(E) * (1+(E/(6.8*u.GeV))**((4.09-1.811)/3))**-3\n", + "\n", + "# GRB 190114C T0+62 s to T0+2454 s (set redshift above to 0.42)\n", + "# def intrinsic_dFdE(E):\n", + "# return PowerLaw(normalization=8.45e-9 / (u.TeV * u.cm**2 * u.s), \n", + "# index=-2.22, \n", + "# e_ref=0.46*u.TeV)(E)\n", + "\n", + "\n", + "# The BOAT (GRB 221009A) @ ~T0+240s: (set redshift above to 0.151)\n", + "# def intrinsic_dFdE(E):\n", + "# return PowerLaw(normalization=208e-8 / (u.TeV * u.cm**2 * u.s), \n", + "# index=-2.455, \n", + "# e_ref=1*u.TeV)(E)\n", + "\n", + "\n", + "# 1ES1011 February 2014 flare\n", + "# def intrinsic_dFdE(E):\n", + "# return PowerLaw(normalization=8.7e-10 / (u.TeV * u.cm**2 * u.s), \n", + "# index=-2.03, e_ref=0.25*u.TeV)(E)\n", + "\n", + "# A log-parabola spectrum:\n", + "# def intrinsic_dFdE(E):\n", + "# return LogParabola(normalization=5e-8 / (u.TeV * u.cm**2 * u.s), \n", + "# a=-2, b =-0.05, \n", + "# e_ref=0.1*u.TeV)(E)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "014476b0", + "metadata": {}, + "outputs": [], + "source": [ + "# After EBL absorption:\n", + "def dFdE(E):\n", + " return intrinsic_dFdE(E) * dominguez.evaluate(E, redshift, 1)" + ] + }, + { + "cell_type": "markdown", + "id": "8613146f", + "metadata": {}, + "source": [ + "## END of user settings\n", + "#" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "80627c34", + "metadata": {}, + "outputs": [], + "source": [ + "# Other settings (change only for tests)\n", + "\n", + "backg_systematics_uncertainty = 0.005 # 0.5% (relative) will be added in quadrature to statistical uncertainty of flux points \n", + "min_signi_in_flux_point = 2 # minimum significance to display a flux point\n", + "\n", + "integral_significance_threshold = 5 # \"Detection significance\"\n", + "integral_min_signal_to_backg_ratio = 0.05 # S/B must be larger or equal than this for detection\n", + "\n", + "if pulsar_mode: # unbiased background can be taken from off-phase\n", + " backg_systematics_uncertainty = 0\n", + " integral_min_signal_to_backg_ratio = 0\n", + "\n", + "\n", + "min_Aeff = 100 *u.m**2 # Minimum required Aeff in Etrue bins. \n", + " # Just to avoid Etrue bins with little MC stats, hence noisy!\n", + "\n", + "# To exclude Ereco bins with too strong deviation of the true energies that fall inside them:\n", + "max_ereco_to_etrue_deviation = 0.5 #  abs(mean(Etrue)/ereco_bin_center - 1) < max_ereco_to_etrue_deviation" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "24e8394c", + "metadata": {}, + "outputs": [], + "source": [ + "erecobins = background_data[zd_selection_backg].Ereco_min_TeV.to_numpy()\n", + "erecobins = np.append(erecobins, background_data[zd_selection_backg].Ereco_max_TeV.to_numpy()[-1])" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "056ea33b", + "metadata": {}, + "outputs": [], + "source": [ + "etruebins = gamma_data[zd_selection_gamma].Etrue_min_TeV.to_numpy()\n", + "etruebins = np.append(etruebins, gamma_data[zd_selection_gamma].Etrue_max_TeV.to_numpy()[-1])\n", + "\n", + "effective_area = gamma_data.Aeff_m2[zd_selection_gamma].to_numpy()*u.m**2\n", + "\n", + "# Model to characterize the energy migration matrix (skewnorm or moyal):\n", + "emig_model = gamma_data[zd_selection_gamma].emig_model.to_numpy()\n", + "\n", + "# Parameters to characterize the energy migration matrix:\n", + "loc = gamma_data[zd_selection_gamma].emig_mu_loc.to_numpy()\n", + "scale = gamma_data[zd_selection_gamma].emig_mu_scale.to_numpy()\n", + "a = gamma_data[zd_selection_gamma].emig_mu_a.to_numpy()\n", + "\n", + "\n", + "# Now we extrapolate Aeff to higher Etrue by using the same value of the highest available energy:\n", + "# (it is better than having zeros!) We also assume the same E-migration\n", + "factor = etruebins[-1]/etruebins[-2] # step in energy in each bin\n", + "while etruebins[-1] < 80: # extend to 80 TeV at least\n", + " etruebins = np.append(etruebins, etruebins[-1]*factor)\n", + " effective_area = np.append(effective_area, effective_area[-1])\n", + " emig_model = np.append(emig_model, emig_model[-1])\n", + " loc = np.append(loc, loc[-1])\n", + " scale = np.append(scale, scale[-1])\n", + " a = np.append(a, a[-1])" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "1b79d5b2", + "metadata": {}, + "outputs": [], + "source": [ + "etruebincenters = (etruebins[:-1]*etruebins[1:])**0.5\n", + "erecobincenters = (erecobins[:-1]*erecobins[1:])**0.5" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "2038bb5b", + "metadata": {}, + "outputs": [], + "source": [ + "# Exclude too low eff. areas (in general, unreliable due to low MC stats)\n", + "\n", + "effective_area[effective_area= etruebincenters[np.where(effective_area>0)[0][0]]" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "48a34458", + "metadata": {}, + "outputs": [], + "source": [ + "# Integrate dF/dE within a bin, assuming power-law approximation within it:\n", + "def integrate (dfde1, dfde2, e1, e2):\n", + " # We cannot let numpy deal with the units, sometimes rounding leads to wrong units in result!\n", + " # like TeV^(1e-15) :-D\n", + " if (dfde1 == 0) | (dfde2 == 0): # May happen, e.g because of EBL or a strong cutoff\n", + " return 0\n", + " \n", + " # In power-law approximation:\n", + " gamma = np.log(dfde2/dfde1) / np.log(e2/e1)\n", + " e1tev = e1.to_value(u.TeV)\n", + " e2tev = e2.to_value(u.TeV)\n", + " \n", + " integral = (dfde1.to_value(1/(u.TeV * u.cm**2 * u.s)) / \n", + " (gamma+1) * e1tev**(-gamma) * \n", + " (e2tev**(gamma+1) - e1tev**(gamma+1))\n", + " )\n", + " return integral # (1/u.s/u.cm**2)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "54b22527", + "metadata": {}, + "outputs": [], + "source": [ + "integrated_flux = []\n", + "for etruemin, etruemax in zip(etruebins[:-1], etruebins[1:]):\n", + " integrated_flux.append(integrate(dFdE(etruemin*u.TeV),\n", + " dFdE(etruemax*u.TeV),\n", + " etruemin*u.TeV, \n", + " etruemax*u.TeV))\n", + "integrated_flux = np.array(integrated_flux)\n", + "\n", + "# Too strong EBL absorption produces NaNs in high-E bins, jus trplace by 0's:\n", + "integrated_flux[np.isnan(integrated_flux)] = 0\n", + " \n", + "integrated_flux = np.array(integrated_flux) * 1/(u.s * u.cm**2)" + ] + }, + { + "cell_type": "markdown", + "id": "a0d668ca", + "metadata": {}, + "source": [ + "## Effective area" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "4eff20a0", + "metadata": {}, + "outputs": [], + "source": [ + "plt.figure(figsize=(8,4))\n", + "plt.plot(etruebincenters, effective_area)\n", + "plt.xscale('log')\n", + "plt.yscale('log')\n", + "plt.xlabel('Etrue (TeV)')\n", + "plt.ylabel('Aeff (m2)')\n", + "plt.grid()" + ] + }, + { + "cell_type": "markdown", + "id": "bb2fb8b7", + "metadata": {}, + "source": [ + "## Background rate" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "9c89f6e6", + "metadata": {}, + "outputs": [], + "source": [ + "plt.figure(figsize=(8,4))\n", + "\n", + "bgrate = background_data[zd_selection_backg].BckgRate_per_second.to_numpy()\n", + "\n", + "if pulsar_mode: # Scale background to the on-phase:\n", + " bgrate *= on_phase_interval\n", + "\n", + "bgrate *= background_increase_factor\n", + "\n", + "plt.plot(erecobincenters[ereco_mask], bgrate[ereco_mask])\n", + "plt.xlabel('Ereco (TeV)')\n", + "plt.ylabel('Background rate within theta cut\\n (events/s) in Ereco bins')\n", + "plt.yscale('log')\n", + "plt.xscale('log')\n", + "plt.grid()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "5164274a", + "metadata": { + "scrolled": true + }, + "outputs": [], + "source": [ + "total_gamma_rate = (integrated_flux*effective_area).to(1/u.s)\n", + "print(f'Total gamma rate after cuts: {total_gamma_rate.sum().to_value(1/u.s):.3f} events/s')" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "6a415916", + "metadata": {}, + "outputs": [], + "source": [ + "num_realizations_a = 100\n", + "num_realizations_b = 100\n", + "# Number of realizations (random numbers taken from skewnorm or moyal) for E-migration simulation\n", + "\n", + "num_realizations = num_realizations_a * num_realizations_b" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "25b8fe0c", + "metadata": {}, + "outputs": [], + "source": [ + "fine_etrue_binning = np.logspace(-2.5, 3.5, 601) # Just for calculation of mean Etrue within an Ereco bin\n", + "\n", + "total_bg_counts = effective_obs_time.to_value(u.s) * bgrate\n", + "\n", + "total_signal_counts_2d = np.zeros(shape=[len(total_bg_counts), len(fine_etrue_binning)-1])\n", + "\n", + "for ietrue in range(len(etruebincenters)):\n", + " \n", + " gamma_rate = total_gamma_rate[ietrue]\n", + "\n", + " \n", + " emin= etruebins[ietrue]\n", + " emax = etruebins[ietrue+1]\n", + "\n", + " etrue_values = np.exp(np.log(emin) + np.log(emax/emin) * np.random.uniform(0, 1, num_realizations_a))\n", + " \n", + " for etrue in etrue_values:\n", + " # print(f'{ietrue}: {etrue:.4f} TeV, {gamma_rate}, {loc[ietrue]:.4f}, {scale[ietrue]:.4f}, {emig_model[ietrue]}')\n", + " # Possible alternative: simulate according to the (interpolated) events vs Etrue graph,\n", + " # instead of all energies being equal to the bin center\n", + "\n", + " if emig_model[ietrue] == 'skewnorm':\n", + " ereco = etrue*skewnorm.rvs(a[ietrue], loc[ietrue], scale[ietrue], \n", + " num_realizations_b) \n", + " elif emig_model[ietrue] == 'moyal':\n", + " ereco = etrue*moyal.rvs(loc[ietrue], scale[ietrue], \n", + " num_realizations_b)\n", + " else:\n", + " continue\n", + "\n", + " # counts, _ = np.histogram(ereco, bins=erecobins)\n", + " # Histogram in very fine Etrue bins, for later calculation of mean true energy of\n", + " # gammas within each Ereco bin:\n", + " counts, _, _ = np.histogram2d(ereco, etrue*np.ones_like(ereco), bins=[erecobins, fine_etrue_binning]) \n", + "\n", + "\n", + " # Integrate in Ereco bins and add to total:\n", + " total_signal_counts_2d += (gamma_rate * effective_obs_time * counts / num_realizations)\n", + "\n", + " \n", + "total_signal_counts = np.sum(total_signal_counts_2d, axis=1) # Integrate in Ereco bins\n", + "# Now set to zero values below the \"reliable minimum energy\":\n", + "total_signal_counts[~ereco_mask] = 0\n", + "total_bg_counts[~ereco_mask] = 0\n", + " " + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "31415b87", + "metadata": {}, + "outputs": [], + "source": [ + "plt.errorbar(0.5*(erecobins[:-1]+erecobins[1:]), total_signal_counts, total_signal_counts**0.5,\n", + " label='gammas', fmt='o')\n", + "plt.errorbar(0.5*(erecobins[:-1]+erecobins[1:]), total_bg_counts, total_bg_counts**0.5,\n", + " label='Background')\n", + "\n", + "# PLOT GAMMAS VS. ETRUE FOR COMPARISON, correcting for the different bin width:\n", + "plt.plot(etruebincenters, total_gamma_rate*effective_obs_time.to_value(u.s)*\n", + " len(etruebincenters)/len(erecobincenters), '--', alpha=0.5, label='gammas vs. Etrue')\n", + "\n", + "plt.xscale('log')\n", + "plt.yscale('log')\n", + "plt.legend()\n", + "plt.xlabel('Ereco (TeV)')\n", + "plt.ylabel('Number of events after cuts in Ereco bins')\n", + "plt.grid()\n", + "plt.ylim(0.1, 2*max(np.max(total_bg_counts), np.max(total_signal_counts)))\n" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "1ce3c375", + "metadata": {}, + "outputs": [], + "source": [ + "total_bg_counts" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "a3769732", + "metadata": { + "scrolled": false + }, + "outputs": [], + "source": [ + "mean_etrue_vs_ereco = np.zeros_like(erecobincenters)\n", + "finebincenters = 0.5 * (fine_etrue_binning[1:]+fine_etrue_binning[:-1])\n", + "\n", + "for iereco in range(len(erecobincenters)):\n", + " if np.nansum(total_signal_counts_2d[iereco]) == 0:\n", + " continue\n", + " mean_etrue_vs_ereco[iereco] = (np.nansum(finebincenters * total_signal_counts_2d[iereco]) / \n", + " np.nansum(total_signal_counts_2d[iereco]))" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "e4d6ac4d", + "metadata": {}, + "outputs": [], + "source": [ + "fig = plt.figure(figsize=(12,4))\n", + "\n", + "fig.add_subplot(1, 2, 1)\n", + "plt.scatter(erecobincenters, mean_etrue_vs_ereco)\n", + "plt.xlabel('Ereco bin center (TeV)')\n", + "plt.ylabel('Mean Etrue in bin (TeV)')\n", + "plt.xscale('log')\n", + "plt.yscale('log')\n", + "plt.grid()\n", + "\n", + "fig.add_subplot(1, 2, 2)\n", + "plt.scatter(erecobincenters, mean_etrue_vs_ereco / erecobincenters)\n", + "plt.xlabel('Ereco bin center (TeV)')\n", + "plt.ylabel('Mean Etrue in bin / Ereco bin center')\n", + "plt.ylim(0, 1.1*np.nanmax(mean_etrue_vs_ereco/erecobincenters))\n", + "plt.xscale('log')\n", + "plt.grid()\n", + "\n", + "plt.show()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "80fcfbe2", + "metadata": {}, + "outputs": [], + "source": [ + "# Let's not use bind with too high bias in the spectrum:\n", + "not_too_high_bias = abs(mean_etrue_vs_ereco / erecobincenters - 1) < max_ereco_to_etrue_deviation" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "b22d2e36", + "metadata": {}, + "outputs": [], + "source": [ + "significance = li_ma_significance(total_signal_counts+total_bg_counts,\n", + " total_bg_counts/alpha, alpha)\n", + "\n", + "# integrating from each Ereco to max Ereco\n", + "integral_signal_counts = np.cumsum(total_signal_counts[::-1])[::-1]\n", + "integral_bg_counts = np.cumsum(total_bg_counts[::-1])[::-1]\n", + "\n", + "integral_significance = li_ma_significance(integral_signal_counts+integral_bg_counts,\n", + " integral_bg_counts/alpha, alpha)\n", + "integral_signal_to_background_ratio = integral_signal_counts / integral_bg_counts" + ] + }, + { + "cell_type": "markdown", + "id": "d27e4dfc", + "metadata": {}, + "source": [ + "## Signal to background ratio and significance in Ereco bins" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "2a3d41c1", + "metadata": {}, + "outputs": [], + "source": [ + "signal_to_background_ratio = total_signal_counts / total_bg_counts" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "9b5e3aa8", + "metadata": { + "scrolled": false + }, + "outputs": [], + "source": [ + "fig = plt.figure(figsize=(12,4))\n", + "fig.add_subplot(1, 2, 1)\n", + "plt.plot(erecobincenters, signal_to_background_ratio)\n", + "plt.scatter(erecobincenters, signal_to_background_ratio)\n", + "plt.grid()\n", + "plt.xlabel('Ereco (TeV)')\n", + "plt.ylabel('Signal / background ratio')\n", + "plt.yscale('log')\n", + "plt.xscale('log')\n", + "plt.ylim(1e-3, 2*np.nanmax(signal_to_background_ratio))\n", + "\n", + "fig.add_subplot(1, 2, 2)\n", + "\n", + "# Conditions to consider a bin reliable\n", + "reliable_points = ((signal_to_background_ratio >= backg_systematics_uncertainty) &\n", + " not_too_high_bias)\n", + "\n", + "plt.scatter(erecobincenters[reliable_points], significance[reliable_points], \n", + " label='Valid bins')\n", + "\n", + "plt.scatter(erecobincenters[~reliable_points], significance[~reliable_points], label='Unreliable bins', color='orange')\n", + "\n", + "plt.grid()\n", + "plt.legend()\n", + "plt.xlabel('Ereco (TeV)')\n", + "plt.ylabel('Li & Ma significance per Ereco bin')\n", + "plt.xscale('log')" + ] + }, + { + "cell_type": "markdown", + "id": "ca520d68", + "metadata": {}, + "source": [ + "## Integral significance (for Ereco>xx)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "0e1df16e", + "metadata": {}, + "outputs": [], + "source": [ + "plt.figure(figsize=(8,4))\n", + "\n", + "mask = (integral_significance >= integral_significance_threshold)\n", + "\n", + "plt.scatter(erecobins[:-1][mask & reliable_points], integral_significance[mask & reliable_points], \n", + " label=f'Significance $\\geq$ {integral_significance_threshold} sigma')\n", + "\n", + "if (mask & ~reliable_points).sum() > 0:\n", + " plt.scatter(erecobins[:-1][mask & ~reliable_points], integral_significance[mask & ~reliable_points], \n", + " label=f'Significance >= {integral_significance_threshold}, \\nUNRELIABLE points!')\n", + "\n", + "plt.scatter(erecobins[:-1][~mask], integral_significance[~mask],\n", + " label=f'Significance < {integral_significance_threshold} sigma')\n", + "plt.grid()\n", + "plt.xlabel('Minimum Ereco (TeV)')\n", + "plt.ylabel('Li & Ma integral significance')\n", + "plt.xscale('log')\n", + "\n", + "if (mask & reliable_points).sum() > 0:\n", + " print('*************************')\n", + " print('Detection successful! :-D')\n", + "else:\n", + " print('****************')\n", + " print('No detection :-C')\n", + "plt.legend()\n", + "plt.show()" + ] + }, + { + "cell_type": "markdown", + "id": "04cd5759", + "metadata": {}, + "source": [ + "## Simulated SED from observation (Asimov dataset)\n", + "The fluxes & uncertainties are computed in the Ereco bins, then placed at the mean Etrue of the gamma events falling within the bin, and at the \"expected\" flux level. \n", + "In reality proper energy unfolding will be needed... what is shown below is just a an estimate of what the spectrum will look like - but pretty accurate for the purpose of observation planning." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "64661959", + "metadata": {}, + "outputs": [], + "source": [ + "noff = total_bg_counts / alpha\n", + "non = total_signal_counts + total_bg_counts\n", + "# excess = non - alpha * noff (= total_signal_counts in Asimov dataset)\n", + "stat_excess_error = (non + alpha**2 * noff)**0.5\n", + "\n", + "np.seterr(divide='ignore')\n", + "relative_stat_excess_error = stat_excess_error/total_signal_counts\n", + "\n", + "syst_excess_error = backg_systematics_uncertainty * total_bg_counts\n", + "\n", + "relative_syst_excess_error = syst_excess_error / total_signal_counts\n", + "total_relative_excess_error = np.hypot(relative_stat_excess_error, relative_syst_excess_error)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "4d56445f", + "metadata": {}, + "outputs": [], + "source": [ + "displayed_points = reliable_points & (significance > min_signi_in_flux_point)\n", + "\n", + "if displayed_points.sum() > 0:\n", + "\n", + " SED = ((mean_etrue_vs_ereco[displayed_points]*u.TeV)**2 * \n", + " dFdE(mean_etrue_vs_ereco[displayed_points]*u.TeV))\n", + "\n", + " SED_stat_error = SED*relative_stat_excess_error[displayed_points]\n", + " SED_total_error = SED*total_relative_excess_error[displayed_points]\n", + "\n", + " plt.figure(figsize=(8,4))\n", + " intrinsic_SED = (finebincenters*u.TeV)**2 * intrinsic_dFdE(finebincenters*u.TeV)\n", + "\n", + " if redshift > 0:\n", + " plt.plot(finebincenters, intrinsic_SED, '--', color='lightgrey', label='intrinsic')\n", + "\n", + "\n", + " SED_fine = (fine_etrue_binning*u.TeV)**2 * dFdE(fine_etrue_binning*u.TeV)\n", + " plt.plot(fine_etrue_binning, SED_fine, '--')\n", + "\n", + "\n", + " if not pulsar_mode:\n", + " plt.errorbar(mean_etrue_vs_ereco[displayed_points], SED, yerr=SED_total_error, \n", + " fmt='o', markersize=2, \n", + " label='Observed, stat+syst')\n", + " plt.errorbar(mean_etrue_vs_ereco[displayed_points], SED, yerr=SED_stat_error, \n", + " fmt='o', markersize=2, \n", + " label='Observed, stat-only')\n", + "\n", + " plt.yscale('log')\n", + " plt.xscale('log')\n", + " plt.xlabel('Energy (TeV)')\n", + " plt.ylabel(f'SED ({SED.unit})')\n", + " plt.ylim(np.nanmin((SED-SED_total_error)).value*0.1, \n", + " np.nanmax((SED+SED_total_error)).value*5)\n", + "\n", + " plt.xlim(mean_etrue_vs_ereco[displayed_points][0]/3, mean_etrue_vs_ereco[displayed_points][-1]*3)\n", + " plt.legend()\n", + " plt.grid()\n", + " \n", + "else:\n", + " print(\"Nothing to display. No valid flux points!\")" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "5215017d", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.11.9" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}