From e857b01ae94875b445f4980e993a561eedf32f63 Mon Sep 17 00:00:00 2001
From: Ludwig Neste <31670556+The-Ludwig@users.noreply.github.com>
Date: Mon, 18 Nov 2024 21:54:04 +0100
Subject: [PATCH] Fix pdgid (#104)

* fix #103

* add tests for bug #103

* update corsika version

* update corsika version in tests

* update corsika compare files
---
 .github/workflows/ci.yml      |   2 +-
 .github/workflows/release.yml |   2 +-
 panama/read.py                |   6 ++++++
 pdm.lock                      |  37 ++++++++++++++++++++++++++++++++--
 tests/files/compare/DAT000000 | Bin 137640 -> 137640 bytes
 tests/files/noEHIST/DAT101001 | Bin 0 -> 22940 bytes
 tests/test_cli.py             |   4 ++--
 tests/test_read.py            |   8 ++++++++
 tests/test_run.py             |   4 ++--
 9 files changed, 55 insertions(+), 8 deletions(-)
 create mode 100644 tests/files/noEHIST/DAT101001

diff --git a/.github/workflows/ci.yml b/.github/workflows/ci.yml
index 559e1b9..84ae294 100644
--- a/.github/workflows/ci.yml
+++ b/.github/workflows/ci.yml
@@ -15,7 +15,7 @@ jobs:
     strategy:
       matrix:
         python-version: ["3.8"]
-        corsika-version: ["77500"]
+        corsika-version: ["77550"]
         os: [ubuntu-latest]
     runs-on: ${{ matrix.os }}
     steps:
diff --git a/.github/workflows/release.yml b/.github/workflows/release.yml
index 4766d9a..b1a394d 100644
--- a/.github/workflows/release.yml
+++ b/.github/workflows/release.yml
@@ -15,7 +15,7 @@ jobs:
       matrix:
         python-version: ["3.11"]
         poetry-version: ["1.3"]
-        corsika-version: ["77500"]
+        corsika-version: ["77550"]
         os: [ubuntu-latest]
     runs-on: ${{ matrix.os }}
     permissions:
diff --git a/panama/read.py b/panama/read.py
index 3ab62eb..16fb5d1 100644
--- a/panama/read.py
+++ b/panama/read.py
@@ -386,6 +386,8 @@ def add_mother_columns(
         else False
         for pdgid in pdgids
     }
+    # explicitly force an invalid PDGID to not be charm
+    has_charm[PDGID_ERROR_VAL] = False
 
     # this follows the MCEq definition
     lifetimes = {
@@ -395,6 +397,8 @@ def add_mother_columns(
     for pdgid in lifetimes:
         if lifetimes[pdgid] is None:
             lifetimes[pdgid] = 0
+    # explicitly force invalid PDGID to have infinite lifetime
+    lifetimes[PDGID_ERROR_VAL] = inf
 
     is_resonance = {
         pdgid: "*" in Particle.from_pdgid(pdgid).name
@@ -402,6 +406,8 @@ def add_mother_columns(
         else False
         for pdgid in pdgids
     }
+    # explicitly force an invalid PDGID to not be a resonance
+    is_resonance[PDGID_ERROR_VAL] = False
 
     df_particles["mother_lifetimes"] = df_particles["mother_pdgid"].map(
         lifetimes, na_action=None
diff --git a/pdm.lock b/pdm.lock
index feb87e0..0e10def 100644
--- a/pdm.lock
+++ b/pdm.lock
@@ -4,8 +4,11 @@
 [metadata]
 groups = ["default", "dev", "hdf"]
 strategy = ["cross_platform"]
-lock_version = "4.4.1"
-content_hash = "sha256:841616f19ebdf6cfc93ad4b8292b98d3747ecd22d145988992325f2b0321245a"
+lock_version = "4.5.0"
+content_hash = "sha256:e2bcaec8f930c4c515ce93e885c5b0128ee92cb8db45f7272db99c43c742ac96"
+
+[[metadata.targets]]
+requires_python = ">=3.8,<3.12"
 
 [[package]]
 name = "accessible-pygments"
@@ -61,6 +64,7 @@ requires_python = ">=3.7"
 summary = "Argon2 for Python"
 dependencies = [
     "argon2-cffi-bindings",
+    "typing-extensions; python_version < \"3.8\"",
 ]
 files = [
     {file = "argon2_cffi-23.1.0-py3-none-any.whl", hash = "sha256:c670642b78ba29641818ab2e68bd4e6a78ba53b7eff7b4c3815ae16abf91c7ea"},
@@ -119,6 +123,7 @@ version = "2.4.1"
 summary = "Annotate AST trees with source code positions"
 dependencies = [
     "six>=1.12.0",
+    "typing; python_version < \"3.5\"",
 ]
 files = [
     {file = "asttokens-2.4.1-py2.py3-none-any.whl", hash = "sha256:051ed49c3dcae8913ea7cd08e46a606dba30b79993209636c4875bc1d637bc24"},
@@ -143,6 +148,9 @@ name = "attrs"
 version = "23.1.0"
 requires_python = ">=3.7"
 summary = "Classes Without Boilerplate"
+dependencies = [
+    "importlib-metadata; python_version < \"3.8\"",
+]
 files = [
     {file = "attrs-23.1.0-py3-none-any.whl", hash = "sha256:1f28b4522cdc2fb4256ac1a020c78acf9cba2c6b461ccd2c126f3aa8e8335d04"},
     {file = "attrs-23.1.0.tar.gz", hash = "sha256:6279836d581513a26f1bf235f9acd333bc9115683f14f7e8fae46c98fc50e015"},
@@ -155,6 +163,7 @@ requires_python = ">=3.7"
 summary = "Internationalization utilities"
 dependencies = [
     "pytz>=2015.7; python_version < \"3.9\"",
+    "setuptools; python_version >= \"3.12\"",
 ]
 files = [
     {file = "Babel-2.13.1-py3-none-any.whl", hash = "sha256:7077a4984b02b6727ac10f1f7294484f737443d7e2e66c5e4380e41a3ae0b4ed"},
@@ -194,6 +203,7 @@ dependencies = [
     "pathspec>=0.9.0",
     "platformdirs>=2",
     "tomli>=1.1.0; python_full_version < \"3.11.0a7\"",
+    "typed-ast>=1.4.2; python_version < \"3.8\" and implementation_name == \"cpython\"",
     "typing-extensions>=3.10.0.0; python_version < \"3.10\"",
 ]
 files = [
@@ -405,6 +415,7 @@ requires_python = ">=3.7"
 summary = "Composable command line interface toolkit"
 dependencies = [
     "colorama; platform_system == \"Windows\"",
+    "importlib-metadata; python_version < \"3.8\"",
 ]
 files = [
     {file = "click-8.1.7-py3-none-any.whl", hash = "sha256:ae74fb96c20a0277a1d615f1e4d73c8414f5a98db8b799a7931d1582f3390c28"},
@@ -451,6 +462,7 @@ requires_python = ">=3.8"
 summary = "Python library for calculating contours of 2D quadrilateral grids"
 dependencies = [
     "numpy<2.0,>=1.16; python_version <= \"3.11\"",
+    "numpy<2.0,>=1.26.0rc1; python_version >= \"3.12\"",
 ]
 files = [
     {file = "contourpy-1.1.1-cp310-cp310-macosx_10_9_x86_64.whl", hash = "sha256:46e24f5412c948d81736509377e255f6040e94216bf1a9b5ea1eaa9d29f6ec1b"},
@@ -790,6 +802,9 @@ name = "fqdn"
 version = "1.5.1"
 requires_python = ">=2.7, !=3.0, !=3.1, !=3.2, !=3.3, !=3.4, <4"
 summary = "Validates fully-qualified domain names against RFC 1123, so that they are acceptable to modern bowsers"
+dependencies = [
+    "cached-property>=1.3.0; python_version < \"3.8\"",
+]
 files = [
     {file = "fqdn-1.5.1-py3-none-any.whl", hash = "sha256:3a179af3761e4df6eb2e026ff9e1a3033d3587bf980a0b1b2e1e5d08d7358014"},
     {file = "fqdn-1.5.1.tar.gz", hash = "sha256:105ed3677e767fb5ca086a0c1f4bb66ebc3c100be518f0e0d755d9eae164d89f"},
@@ -841,6 +856,7 @@ version = "6.8.0"
 requires_python = ">=3.8"
 summary = "Read metadata from Python packages"
 dependencies = [
+    "typing-extensions>=3.6.4; python_version < \"3.8\"",
     "zipp>=0.5",
 ]
 files = [
@@ -1205,6 +1221,9 @@ name = "kiwisolver"
 version = "1.4.5"
 requires_python = ">=3.7"
 summary = "A fast implementation of the Cassowary constraint solver"
+dependencies = [
+    "typing-extensions; python_version < \"3.8\"",
+]
 files = [
     {file = "kiwisolver-1.4.5-cp310-cp310-macosx_10_9_universal2.whl", hash = "sha256:05703cf211d585109fcd72207a31bb170a0f22144d68298dc5e61b3c946518af"},
     {file = "kiwisolver-1.4.5-cp310-cp310-macosx_10_9_x86_64.whl", hash = "sha256:146d14bebb7f1dc4d5fbf74f8a6cb15ac42baadee8912eb84ac0b3b2a3dc6ac3"},
@@ -1777,6 +1796,9 @@ name = "overrides"
 version = "7.4.0"
 requires_python = ">=3.6"
 summary = "A decorator to automatically detect mismatch when overriding a method."
+dependencies = [
+    "typing; python_version < \"3.5\"",
+]
 files = [
     {file = "overrides-7.4.0-py3-none-any.whl", hash = "sha256:3ad24583f86d6d7a49049695efe9933e67ba62f0c7625d53c59fa832ce4b8b7d"},
     {file = "overrides-7.4.0.tar.gz", hash = "sha256:9502a3cca51f4fac40b5feca985b6703a5c1f6ad815588a7ca9e285b9dca6757"},
@@ -1895,6 +1917,9 @@ files = [
 name = "pickleshare"
 version = "0.7.5"
 summary = "Tiny 'shelve'-like database with concurrency support"
+dependencies = [
+    "pathlib2; python_version in \"2.6 2.7 3.2 3.3\"",
+]
 files = [
     {file = "pickleshare-0.7.5-py2.py3-none-any.whl", hash = "sha256:9649af414d74d4df115d5d718f82acb59c9d418196b7b4290ed47a12ce62df56"},
     {file = "pickleshare-0.7.5.tar.gz", hash = "sha256:87683d47965c1da65cdacaf31c8441d12b8044cdec9aca500cd78fc2c683afca"},
@@ -1968,6 +1993,9 @@ name = "platformdirs"
 version = "3.11.0"
 requires_python = ">=3.7"
 summary = "A small Python package for determining appropriate platform-specific dirs, e.g. a \"user data dir\"."
+dependencies = [
+    "typing-extensions>=4.7.1; python_version < \"3.8\"",
+]
 files = [
     {file = "platformdirs-3.11.0-py3-none-any.whl", hash = "sha256:e9d171d00af68be50e9202731309c4e658fd8bc76f55c11c7dd760d023bda68e"},
     {file = "platformdirs-3.11.0.tar.gz", hash = "sha256:cf8ee52a3afdb965072dcc652433e0c7e3e40cf5ea1477cd4b3b1d2eb75495b3"},
@@ -2123,6 +2151,7 @@ summary = "pytest: simple powerful testing with Python"
 dependencies = [
     "colorama; sys_platform == \"win32\"",
     "exceptiongroup>=1.0.0rc8; python_version < \"3.11\"",
+    "importlib-metadata>=0.12; python_version < \"3.8\"",
     "iniconfig",
     "packaging",
     "pluggy<2.0,>=0.12",
@@ -2980,6 +3009,7 @@ summary = "Virtual Python Environment builder"
 dependencies = [
     "distlib<1,>=0.3.7",
     "filelock<4,>=3.12.2",
+    "importlib-metadata>=6.6; python_version < \"3.8\"",
     "platformdirs<4,>=3.9.1",
 ]
 files = [
@@ -2991,6 +3021,9 @@ files = [
 name = "wcwidth"
 version = "0.2.9"
 summary = "Measures the displayed width of unicode strings in a terminal"
+dependencies = [
+    "backports-functools-lru-cache>=1.2.1; python_version < \"3.2\"",
+]
 files = [
     {file = "wcwidth-0.2.9-py2.py3-none-any.whl", hash = "sha256:9a929bd8380f6cd9571a968a9c8f4353ca58d7cd812a4822bba831f8d685b223"},
     {file = "wcwidth-0.2.9.tar.gz", hash = "sha256:a675d1a4a2d24ef67096a04b85b02deeecd8e226f57b5e3a72dbb9ed99d27da8"},
diff --git a/tests/files/compare/DAT000000 b/tests/files/compare/DAT000000
index 7a968c97bd466e213812b187fbca7ddc5c219404..60e74b4bad3f8382854a7c2de977e81a659f0551 100644
GIT binary patch
delta 75
zcmZ3nnPbIf4xuTL3=Bb`ejW@!&`_T3@l9j$KNg#feydqn!Mx4?SdR0+nSIQgBSag|
QLio-19&Nw(i1Cye03=c$C;$Ke

delta 75
zcmZ3nnPbIf4xuTL3=Bb`ejW@!(6Av9NKgL9Vzbe2H47`4xA`B-aUM9ck9l*1XyaK3
Pzxm#y?e`uro-zXfufH2?

diff --git a/tests/files/noEHIST/DAT101001 b/tests/files/noEHIST/DAT101001
new file mode 100644
index 0000000000000000000000000000000000000000..9dff7a0fecaa8df51f677f5ca75cf6e27a2e0e16
GIT binary patch
literal 22940
zcmeIa2UOHt+cruSQ0Xe7GzAeA5J4c>$xxGwfT-B9V8MnRdqJs+A}Sz7sVbnTAc)9J
zGT2~}L9AdeD4+<~uz>}XKV#za^ZI@7`#)!`bIv;7Sx?rQHQD>R?t9<+-h1D>WG0FJ
z3JMdaj2oliwgpqzxNMAq!dEv1g)mozg^FGZf8)A-&H9ppLZt`yVdYra>(|B#dfqoj
z4McvB{Ik5tV4OmTo|l5f=f8AnWFZ2@Qa#tw1=9rz3TAHgTU6W}9F*O!U%SSpxl79G
zGT35*Dze{NvVP!q{`G79JT6MiV5x0C_f>kzNZqVw6XfN(`So<IXsdKhx-<cdQqM<!
z^SSTC?*E$p_bC;_oh~sm7qN9wP$>El<y!GE#>H7f!}H%u{9i@ne&Vj~_5Gd4*2RUr
zXE|=a{3*X@ed_YVYK7~8`u(n{hx&TTKkg%gT~2<O=n^dG>GE|_ATOh!aJqD`>!azT
zTobO&an;MUakZav%VpbaBR7KsL4TF|zsr)}vGz=;5dLePg2J`gS6!MqU{|rpde`r7
zmjQVj{mJQ}XW=TaR{0mbOX6)@n>zkoS$}HyC;k712kHzpudr=AM!{koI}X=*hb8y-
z3ET%NqVJtApY_t{SNnB6hERpO2QgLS{g{SBHa4%98Yb@>$HZ&A4lS7K^<msVuk0w)
z%WB!masO)b7*BlaWHxrl`?dTdj|befxv=K1G2cC(Mcp_qPyChs>;C`QoBtbP-Jdb;
z-x$osw6pG;TTWj3-@cLi=aSpM^4<PO|10|MQ@{QF?-~DdkTADjK5|PZ{<F=8J$~cO
z!>oON-~aRd{{s)O<8nmgZ8GR(PfFr&3_h(p#B)`ssxV&3R48AcIBk3YIzlZW{X69#
zaYgSj0Hnw^x0?OJ+^V07hPxEOvbGN3UVM=rG`E3{`abX%wzw-6C{&CFHqSR9Rp&n8
z`T6a{ioQm_u=)f`%4~Hm_=;Jh*0?|BTG&Nw(dV$k=glc&Q!V7tfuqRmJs;@#yC(_@
zzHwNZ@ig&FWQ2IZs&~-LKo!|}>tI>tE6$c_Pa$gT!sq0N!>onRWklbQeA!!Dt`2t#
zbH!~N#=#3z7sBnSH*u4$So*>s4trX-lw_Q?!C25^*p|G*#Eo*tW?bdU8|n%}GaB_M
zzm*w)%*tchUmT#33Jx>djDfqpGvKvZ2a$_@<-~md0(|cv4*L>ejC@Eb0XB{ai0X)R
zVaM7WVttMAuexh%RETvZ9D;*hj6}|K*<!aR8RCl{a6Std)C7Mvy9%E<)Dw{9`*HnU
zx%gU`(;1m`j!exQA^ujG3x?H<qgzGKJbhkqa^tOFs_5bj*zkx6a_gc$t^czM>-Uz!
z4!y5(a!`yfDL%R%N^u>Iiyn$SE$^|I*>|jG$PCb6=Pl{kJ568{ppToCm3bb$%gGUf
zV(}!ur%-xE5rn-JGsC6jw00xsGuIUd!T1jd^)~-0vS6yNFtw$Kxiycovvw$sgg$kJ
zVU=W5wdDs9Wuk?z8^K{^D>5a=J9QB|tyE;pMZL<=RGDym*pOd-*=|`RnrLuJ6gajQ
zoHS@Z&h!mquH<uahv-G(@+Y|<_dy^s>(WDdEr=$bx^kF)dl<NN=otA#RR&@<9~CY?
z{(_0}=j=2rbE2$YG1SAsGf_k2EuFk6jG;es^^2={Pnn!`2HM8a$f;ebve^ZSL{KfK
zQ~K}|{M=MRDK)yeEPI?@v1|1JI`<AI=VH7K2$#JE-v^Z-8VTmYskSLXx`}JsaHU%6
zVPG1$VT-mPaPa^}pn%e@wVa$$*h4DI?<WO+xrmNb3S|1ci)FljaM(OkOUgMXlA0s_
zhJ=^zV%nbHVGa!CY<acuICv?TOx1rJj^wYrjeE|Y%LGs1u-cmS)byZFq~Y#f=slOS
zjPDJc8LGx%5x2vD+m*%e(|eJ~<_L9M(#XzR@3``ML=x0ds4toKA{IH}8pr@Y7uxm<
zhaEqBKzvyxS){YF4)n9wiw}6PgPHiA!={E!fcB<-f=iztfYK0c=}Qp6XszeUo3i+)
zc$cmx)wD_nz4E-XV)l?s<_X^iXIwcZ9woRzW_N`lJwkk>_HEhB1YYOCA|=U<rM=5;
z4yqPt9QBYwDw)hzUhc)pG{=3lv7+JyJ&~Z!5X}4fP^mTF4xvA*9ePyL=i;xzaC!f2
zSaAi9#UZ%K9TDZcv`3l9tgq{mw;|Y|YdpsG`Ro!BX+J$sGK}u)W)ry`GmEwOOXqj1
zo}@xhBz$II7l?mE(y^noWXl$D^?Uv_j~Y2?Gv(fv364F8Wo%y`kVU`cVvk9MAG!TO
z71_Mc3GF(Xi*3v6!z3zmaz+OSQqS`8sour|1^06L&}Z{hnQIR@e~!_@V36M6960a-
zoHwk8>{jp_X;A^!mwe_uraFT5OHLo&g9h}{mz`Z7CX0X1>2#WZntU|u0ySaV2xO(F
zzRd8{E}5$(Cl~6!n|iMwLEX<#LPKr($`YJHWjY@@xx3|RqH85yqSH#xARB`|vhRwa
zve&%MVwJ(b{Y)6B5~(7|ke<x!W{6Cwjgvz(+Nci`zLH1M0+5q>ePtC5+hv`fS&a2{
zTf=0a73~N@1iMgGO?{bn$adK)UZ=9JuB0krgXFPopx|t5Ad~GPVHkdXFw8hiRV7ZB
zEVN8T8}1*Z>D!~3yF7NVa7*ddph@J^!kv(9;6~c<QYK@#m8;*j6CKph8@{BbPYlu^
zj$w)=Wzw}?oSo$_kCU=F&*2WoN64kuo9UDDHSxkHoUcDMdy#Wh!^OK#AAohPsWWe{
zD-(U}IJt1GOz{Y-^XHCUvw&YBnat0`>)7^44oe8RN0}}>3m&Lg2uh2m)3b(n<0JWa
zEuLm5*15PJ)*jjn2WBhC5=#w<3SZ7nhm^<U^tuBgEAOM|jQAYs7lka(Q+(eYJF6MK
zKSP^(=JF0;Q+wi1_j}R(`M9)o-XU?sMPsKW-V``9WF(yw>W=li&(-0orUBJ-^f>%T
ztp;)NT}%(ozKwn1`%*_%x`+wzcFq_u7QQ}vCB0+td#vmcC#UpGhjcO2Cp){oK;}WV
zjH1pk=B*mv&rBbPd&D)0Q-`a=en6esp?yJ`%KOC%8bLm*S^{slo(HbXE1^q+9q<eC
zIfrfcwN25I%l;Y2;;3ZQk}9B;oIP>hH(Yt&C+m|$$Y!XaF$oUqf1l2%8-zb^=djD}
zDrL*ecEc973WC+<S7|1HKNix=VKt}cmG#k@i@5GQA;=IPre~F1!$O~N7-l#Lp3<%c
zw(nCDoGf(09^Jo;ZQ%Wxmsktyc_o8~DJOw_gDw{Pb1N?U$jKe->>!`vD&jz;=}7U0
z+KOVu7Ho@rj$`YPfCwd4-}jOF{?@2VkGhIf$CubOe*QbLUzyxQ_`$w2+))>hS<!X-
z8kWMBrxXweZ|Ii@n$FJ^gj{O#6pydRGI%UEu#8HrJr9Rw?Lt0eo~-Eo=`PmAx6$P4
zzNB7NFIc433Y|4lp`R@Z$6HLeeg<V%Ir+3VIn5??z(j}9*lCqgEQHs&W|uK6Rvij|
zDD*<Njc~!{#E7tP9=p?#;e5s|27+d<MGhq@WA`+!Vmo;37`YdO?R*O8^V%*6p32y^
z%u;NP9AoRAZBT|3%`idle-gRa-23KPJF^xO^K!lBdr}8Z5t3~08R$*NB*vnjuk1k-
zXJ_535739;(c%fudVp}%T;{vfQKq+%n{S^F-z^^LY7BSOmm!aLSW8`UJn<d8KXv;g
z0Q)VysK;AQAsTlt;=R3NJUw4?ZM28c2IrnyL0?W#;7;m&JgYGSKhF0T59^`iljKt%
z^S&wa*<1<#@y&z~_2g{1Zd(L1U)A8{YSyU#K68A=f<Cx|d~L(7TNm&D0=7O~3bv^X
z7A%RKkI()H@MeA-9oOiANT1FFcP^O=oXuC^Z>GVx8$YI#Ci{y%?DI#)_#Z;z@}2O4
zQF(;H8Lqs(vCTm9Z~+|B`UF`xY5?)>PBf;&$E+`%32?^?HQ+q)5m-~FN*EXw;`e)T
za_<BpXw!ERo_@p!r0o2R+oXNK(ZL)>pNWG?roVAM^!hLyF;t0Y4c>sSQ|GXT;v%u@
z+^3Lv+XpxbIfHk8OUA8qIqVhf1A3*INM<)5MruFh3P1m7mF@NB=A#(x81eRo`|$Vp
zPUuFQIT17Hs!)1h&@aCP1Y40t8sVVykq@d-XF)IU+(8KW`PNvooI0hbMa_pHv~0m6
zW@eKyb7vr52P&FWJgZN}yIg<^+H{Df{6);iER$b4>r2*wgZ5*njm~qC-ZLdcuHPc2
zXbvaWXD|upTcnWY9YYZ&cqLJBtCK#Y$YJhRPf$;a?oi5WJCJc#uF&3=nv7X*4jWw>
z35e?k;LhZ?fZp<h5T1FDS6||=y*9zr6I}~f>~#)&xu+;28N>1(AzVC@j12=ixfm(%
zRYODfcoFK~%&~@-T%Ry`X9&)$Z-sY$bwz%D4I%W4x8nVHTa*`FrL^X2kZIFqA<31i
zn5dSX^p{Cod7}OK)Jx(DInn(ilH*fL7r)#_m-F!@HA_*FXSklq{bqp<e)5JMRCbh(
z{mPZ6?wJkw=?)RcUr>QJ0!L|Z)-yWCn$!7G`+{ghq?M>@iVtiWO3{bg3h3*fI612V
ztW?K4l6)E$1>DzW)5Gt7poj762}kt1M3it|QgV8{U~>Ky`bow)dgf<N&aW&2r0Qpa
z+P(##UrsGimYYuO;OC?Dw{jsg+mM|1xhFhj%t^vvZwzC{uUpTlSwjVW0QvyU1%{Ke
z@x#VN*uvj?lHpk(b4)4d+QOhk7jyB}gOyl4kF8cc0c=;F1`WFpqF)+GyzatBtcjm<
z>@$j}{wf6U+3A7gG>7Bot`5crwy<`xu_vHcDRQ7I9Bf;Hqj%hI;y$}tur+*~Hy#uR
z?|HWfPKkFwrWaS^$By>IFY;^X(jNnzGkiZf-8=Ib3j28#-!fhqZ<5dNtj>eY@o3<a
zIcUys$fcCNiCd3V!?*BqFS=UcR(v=$?&e`YwtS;=)fu928P^Uy7ns0MO$m8C&>qbl
z2H_*V&7wto{pPwkQ+^<y^b@AT%)3p(>f(Dc%ie>3^?9vFmdH_)Xk=z)4?+LEiiFmo
z;Y1L>z7Ae_3MhMAhjPD034D!{aO%Wo&tg6f9*eyIIpPxMr9sMwMQIOWslW!;e$2I>
zo*oJn*}S8m)XgrH({Rtg^2<_xJ|;W<%y9}`8%A!uR1Kd#f0LQ-G)0zunHzr_{iEQN
zun_9lGc`oP%-cgn!9u2@#Eplp*R~?;z9ad%BpN&lJuMYV<}>T~wf<L=9I|msICX2$
zCa}AFp{(SKi7XuAV!6-OA>gJU1Kw|119~ipC8ll~#Z)cfe5NjoB=z0|f%m=BVZ+Up
zOnT2CVoo|&UelvQYVyidXWJHaG(K}JlbT}7bn&sg^JtDZO0-6NFTE6;v~gyVI^3B=
zKK4xh>_hq#9Hbt2&PQvE2Qr=wOPLscj6J?-3@a-417A{R!%VVr<=oZJ3FAF%tYh0<
zZ$cL}V^<&$nD0ZfvQ;Xp=TsBDrg3plzwHQY>9~Om$Xy5*r^X7el|;&*SZ=>rdN+x(
zjg6x8k9-3iF&46>nj~gv5a-X~#&)EFyDj|sMmw-Fsi(V!#1jF$Ez^|u!VSsBi0$(f
z^wrKI(y3TBVZ+z&YE&xOZfyZ3UT{Qu6&xdA8(i4b#QE%at_xBNJw=Xd8VP!u#uGE9
z*a$V}a@bS!Fjd&nm-4=;h34m1GlRdO%sBomp>dTL&@CGbBz14WauwFcI^M*WdQR?!
z{Q&r2>k<mvCPLQz3?s1py_ik>zG}t>ZSwj`1&~r01(c?x5@Sy$GI^gkxpzYjiH|MG
zfXW@sAk#BaX4r`1IL@C>Xl&92_kNbcFBfk^r+D1OPnfL0bNKa3i|T3WRm4kp<7xw!
zs`PlzcYOdhp6}P0S?8(LfLKb)I0^M_>%&YO^MyXmpHFBX9RV+2>jOd?vcarjmGtcs
zhw1NnoGnVW4X}&qIQV0xH}cc=2>tO`Ha4M&8xI!gv&5d_DUPpkEm-BoOL`tS>^V-p
zR%F}#gvKt(7DowKT4sqBZfd5Br}f07e0d9fZ<eZTdIfo3W9P=Gp|qJ@QiU1654e5o
z8CZQlii}-B3x<~#O0N~@<F0jFd5hL=7T2GTa~`R71hN;b#Mix=K$!Ap7a`)A5I(uy
z`DdR)P*`dK5!=U_zGA`k`De4|JD0w)5-azlVW55*3wy3ed;Q+y?KKAY@<yolyCiVr
ztQqli;Z1D3C0F;^b-J*v@pST0bRd#HwGtC*<VzpO*C=fLUQGH4e0LWC+vR5k_ih(s
zqsul*H}LD1qSec&T#bQ}!%^=9(?0c;UVEP|t>@!JSkht9ECoXnc~gO$oOHnbozrQl
z0pCWx@A`>V<`<EcM#|`n(jS-KRH)F+JhpoE3b<y#8?r58KU!9Wd1h%|k}l)144($L
z@oNHAqZ}siej4Y28=aC0c`Rn?59h4x7f?~3qj3J{IavR_wb(#;8?k%s%pRs?F1ve`
z4L)*4&@y?9^nA^5I+^!NQuRQ}^3zdDr!5Y+N(agkmY$J~oWiY9EJg*v>kv)oS<)^z
zxpgXGn{-m>rOIJVK@?;?j0BT#KV<KtX+*z;`-SjLt{rHrQIezWVv%cjCK}4bGWU)o
z;zRiPed_8cYSWb?VA+Yi$kf*zj5~gvSSRM%A*Ht`)&HY_+-NZY8L9u3Sr$};&s)f1
zY&+Bu@l<+R3N;|n6B%WHlSU5&;|@=`{!)bvk%$gyf`=;|QCZ*?rgKa>-p;SjV82kZ
zbEPNTW3we(`Q#j<okY+EM>spRW;p=EQzW_g*=I0TJ(Q@{+Jd))a<*ir^#ZlKOsL?+
z+Y!$zS;XO*&y48`4r{eL1Yfj_gPW!M(RKMjSX0am&)xidyKm+NFw$lZx$RQ~I=TO3
z;z#dfdJ7*5#%8NQy*iX&S=d>iTv#l19c4o#*K>8aP&ZCwc^8FVR2YJ2Q;Z2p&?fHj
z*p`9QM7^h?FtkGp(YAGA);!iFv>G_MnTZ~v9~bt5S#U6txdmm4QpOXvzH!(ajaH}6
zhw9-8hR;ER8Nw{MgA;-BH45v`f_*Om@!1Zop<B@51jGzn<4=6%+vvzK8Zvy*1W9de
zku|sW(CrI05VKUdx?hUfDq-qv#MTNBdU+irJ0Ay_X=xmWhD@h~sWf!CegzVsZ^?|g
zYezWo>sG_d{oti0sp9dsze6gX{bVWwLx|Xq+#L62bt?6o%%gI4#31Veon=L76J+yt
za&2^Hwx>ioTPShJI40;4xHGk;=8Vd7u6_wa#iD?7+r<+GRzgeMm1HF`J%sA6oKE{n
z1M%{xJa}bRFW`5thQ_Xo3A>$~ouZymRG+poP#SBCtjM&Iy*!~OOiSj<+tm4x`d;Eo
zg_I^C)+2<n^sCEpQ4A*+-e>~G7UWZ+@v4Y^=xQc&b}Hi@#mS9yT}3KiT}(aLFc6K7
zxWzbp*u=ziadw`+RX}aixCE{5Dh0hJ>&f=TL&U{sPUk=;KrRbSgLTV}gD<)-nJ_<`
z7$l!_*fz3n)svVB+d=%)uZXwROXllrZ{p1mPR{g`3DvF}27l@rg|3?9#5iX>X3Y40
zcFbijJfdnXIcQZbNFLEnPc8k(%vWIL*jTW8>QN*kJ)z7eu%}>oeJB0x^B2a2Z(IMA
zG1OYG%}DJ9b%9x~BU3(fDC5V^0p`68kc;g*sKGx~1)rEk!nw5{bGehN-++Bp@cPNy
z$cM8DG~;v;e#>AT{`Pm=oi*1fHcMObv+w~TJSfEio(3^LQLYY2r!=V|vL~4Nq7fKa
zWHM2Q-Vx8oa+vtYedt)TpZK8gA=G@(n+aHcgE)MP!!+mALE#f0k=0Gtk(L}^eDgkY
znTZOADI7PVqW66P_xG721)cuPBtK)}kQrQ`7{Btg^HPUnq_t5K2n&5)QL8Xh=IzAE
z-KiT$4$f937g?vkV>%Ei8qvu_KH@O^O$CTNZB7L=RsxOh*1~n-SlNnI+#1CzxSUED
z^MEWAk%(jO-L(J5(X!#gIXREP17Pv{bMWF1r{UYwMWTQCTjBx2^_Tpmy{K9G5G?dR
z3B>h-WoG)mL>@nno!GvSs<j9JOX5AiS)W+h8-NHccXI1pBxpLN=KK~|zgvxTDNJI*
z-4+lIZ@BVabh-nDj7Vy$R2?a-c!-A<=in0go{Q~Ed+d(F1_?(2)e(uDth|L&L@~am
zk!wE{r(ghDn&5TRBA|I;9kFHqNjzScv(wmnC>eH!-570e1Xjb(&_hddnDsBY@=9z*
zfV*@b>SNtO)S##yzu%*kzRu6@s;kFPH?J$g7AAMVh`v^`Ac>J|<W)}RV2cm1O|3n^
zOyiJ2R(qIN@8`<Sui|vR*2|=hv?NlkeuZeHiaVp*cNSyDpJ%zdhJhU#%f)3@6T!8<
zvBF8mlV~cQYum6r+axpLOe$gZ4b*+nc-c(9W@b9S_cf`zM-8(}0_OdDBA*jxS8kO}
zWOP?q{@Op}U9J#o&UT|Fof(OkKbVQ{bUMno@Mj6<1M<nrj_u;%Pjt{N3vSSRHO!gm
z{2FCVyaXt%+D}F~9R@^MtPt%qmK~kQ^@+z+uJ~c#Te!6~3bEbRMfZQ{kB9JnIqo_c
zj=N>*JhgfQ+AseR)+bHksl}g_v_C%tm+krpw>VUR%cj=En{Prw-I?>tnf}(0ul0A(
zyI2P`vR+$pq5c+;#jlC-&RSAY{R`kw$1LQ?<|yW^q64A*`^;md63CNo0IF}-qC*@?
z==1>*gb_bKG(U+zwhsSB4GdQnxaq7CrcZgulr7-<(uY(y>G)E(&?p|gy=D{Mw!Q)%
z!rSTKGXxYp8X-F6a29>IlEPDh!f?Q^iS9L>g^qr>2!VqoNSgr0w{Q4>-Qd?mr}plJ
zv*xA2w=b>$i8>L4Y;qW}m~Z!EU#F70AGyG`KW@RRumXIRXr^Z~-wtM`$?()9h}7#9
z4boFRh<<OQ=_~x%)2G9EaC+WUFff4zA7_ptp6K-_1QWPAEMKAm8LNySS1z=HQH4>2
zLUb=i*PX-i6X${8^o8VwQ9Tfi2%MOSrqk7i9F}!S5gIzlhkSL2gpAf;#Fh)=nV1_K
zmeV>7?rij-Qnsdo4||-5J$<&x&MoEI9hvZ$df)FA)wtvVBDpwGW;@PC_A`j<^P*lj
zb?9-0IQ3^cs<dnk9k&%?4EeqGl(;A4sG*MVBo9wSZ$~8a_)HqDq`{T<n7Z#=8B++~
zIS~cZXYa~l-|c6FKe;uS_jpT6-&~ti(0mKmsxOiCuUH|>b>#M3?j6<=i@Y*2J2n&@
zbY6owH+dVqjNcdQwT=~8Rz^4*eJO_4q69IvQGlKQz{LW~*EiwIK!@~cm;(D2h7pT4
zn8_A?<LZEPJOJql8Im0T@#w;HxrDva6j?aRVP~>}#TH6^#G+UOIO?O4?8qyG@!!t*
zrT^4c5TuetRd;kE2bcQbJ8rII7W4JHJtz#R6HP}x%@U!RWtPHAy$pmi^SHXZ+uabm
zEx$po$&LmSg)^Bfe22^h;>x4@8X{ReGbkoI7_GQbO(gDG!`P4H<kEKsfUVQD5ZeY#
zblLfAVyBV0usMsvG@trYb8mHu2Y%=!aF0l#JBm}KllXpZHtV!=Y}eCLVp<NI-t>*0
zPHV_|26N?Q`(i-qzXBkcOOQhi9mLGHn~2+foX!dT8bvu{kBV-Ln};xIq;$r<i_+iE
z7!qFQLiTlW;(cjfp>uO~3dfz#7EYVR#h#O&0wtpoDG;*zIyz(Jc4p9l<M@02S)}>G
zgXG(hyP)|R+F<?4MCOHF1@6&@EAPdZr&LGZBjCD|0jlPn$k^}wj34Fut^|1m<1JoP
zZ!#5(%y}#AI2FM-baHZcraMvz%YrDA`omzF%28(dcsKmkL=H>#uY!i0)nLz6pTn7f
zSIP^P_ZMztxc0*Whry}NmSCIC5K!2!AF=cFW<rO5Mr+_>hiH9l1vXJ;Xy8_X@X&Zo
zVR=8+XY8C4FTM!oJKrTQh;)#CG3V*3pmZ#OU!%P2u?LzQ-V2UWc7WWM8xk2-4=~^B
zoGoA14->mCi6-G!DezFk{>;WCMd>U)mb+#7!!L#ba`dHQFyZB5=|ss{>3aTLQ8z-v
z*+RHWyxa<d90LZ?dI59k<>stT)|QB@sgnABGbAH(B<R>T1~gV3BUR_uk%L0hs6bgd
zyv(x>4S#lzwpXjcKJxJ_RYmFCCygejnWK8bb1qBi@C_sJ4E|hdqHt%a$G%g}n6C~z
zCAf;dk{ylf$@e5|c@63psijlWkkXhc^ttL~TKUur`~W`(1YT-@#|4dqhus_u{Fo&~
zr?Wn>LmppPxvtqdlzV<2`An|`_Idc6Uea8J9pvSVM8UAuf}vpfq%Yu|MgwiI^a1vj
zx23@4KBY1v6bx#PL%&)pGGXGuxX$l6XG|&?>I1>4lU$M7(Mn9{ElYg49AnE1ELm5!
zw1+F&xwe-}fSxjwqHguqv$_`TDyM34h;zuKG<c8SNJ7yupE$tptDc&Ug9DeWaPCw+
z4}03|BnAvSLY(KH|Lr)gCf*Y=6k0s+05n!JlAg9&6SL*}XK8k>cv-TMc*w-B@I9;1
z^zD7wo^A4TH+F6M{qsl38ha~Ajm~I+r_h0p8+gp~6@MmoMsYSd?1UNlY6ODJtsl%(
zb#>Ak__gWe)*4Cv=vv9_Z+b5C9^9v2J$_DG^0D<w-V<uj{`Zt?KOGc%Q^rglc!$}(
zjrBEK-m|tOpdFnGUtAWCIO^&X;+}JfwLdw(IO|>n)~W!B+dLYr>#a}dM9d{V@MCkx
z><W0kkr3IZ`x?z3I0v`#*+z_%@2lA79lH)Cfe(-OQSCEJ!MF4Kh1(0?GP5ym94$X`
z8-8b}L}|qgLPp##z|ZQ=U?%f&v5OBt6Phc;K2IQ^p^=E^Wldys_?SFAXrW|4b#E|(
z?tw}^+0cudM=%PVoSjca>#0l4KO_T+-lMM<Mq@cnvzX`n^PGiqA|zAF-vM)CEBfi@
z2pIwP5tj<MF;?*N7<|>d5VRjGN5{Cfd(5o)jPK|B%iM>FWdES8$hrGZ(5hTxI%@oW
zLY?2stex)<x3T9X^z6CF<%Bi(wGE%~8h$TYuyQ+GxPBP?n3WwnYArsj{xjadWAn_8
zKwG~ni$fm+P;&kte)vQU{)JzE-n)Di)V?A~wYr69gVP~=ar9mMHIIc(+zb1^za`FB
zb3}erq~M7C9sCB5U2GFVCyU;RoNg$BTZdC{9lIL5QI4_w`ukHe^wByE$?!`Bf-Nbj
zxW~M^_%J?B?9e+-1}#jWBHw2L_tbriCM9BOvN^wmES^WgZg0i8`T9s_;|AvbOnc^(
z{CNPYb19xHX)x3kJF4fR8UA6+(0m8xJAdYwyap<}xnu#=w0x$ZaatIodcuKe=dqS4
zrf_1Dt~k{`4qEoah8b*|P0akxty{5eaa7f%W8&0o0kT7JD`UL|Fl+d|%w~ZRy#Hd2
zIQ(ZLtae*PwhosvE&SQVwL7umMD;wV+;kF@@O28a%HN!}9M8>Fbt2YhpH_i6f)XU)
z?hfqGqnEgVpKq0Sod)CmX2F+6wu8BMdkBw?+QVGq_c<w5?Ai6LACO&50a!d!O}5w1
zp16?6`K4^l5jgXWkpywMik_+t!Bn)a(c1j^i-Y}Skd?cXe12>v!o;Z&-mjAA5BxI}
z%Yutez6Z=j?`;~P2O)*hh&7|=-tw3wk3Ai#a3-2gO%J?*JQ;RQnkaUrPx1Rkx4HYM
z!kd=P>bo}!Zn}m`mmc3K{r%i`;9Dv6ZMPojVtP`rQWPN_JRx41%C9k7KG}gOtv+Nw
z6-;ot_a#qRWQKG%k3AgHQ~b4;K>Xou8u%fDqzCnT(|h^zv!(sZMf2M9q3{i9(DYqL
z=|M$ee5HKuV8^sy^(#@RLXrquTLuTJRY|R@wdqm({9sZ)lq_F0lf32U4h1|dllGsc
zO=En_(ywu*a{3&mn$Fb%n{Rie8#-<19eiA>Sgj>Ge!Ic>+0OuYRp$w50MR8qD_@7O
zei>=Gm3s9p2)<EdhxR>iLW-Fw(&zX!O1pkPao6D9Vw=cxkmDLJ^}VJ{58|;;Gpb5Q
zNQI>ara|zd;&|y@CxyQ-?G!U|N4Y=jpq!1AL?%cdm9<Kp`0^$M=Tmd$T|<UjWeQR>
z<D~tUy_S0MI&07MpcWS<kmo<Qfr1K>IjMp&ZUMYs4ivypNR{~GvSwIvZ51(9`hceS
z^W^WYR-pa$J=jug3D&iw5eHwKmyRmp<X-FMg1KrL)Yv`?&~QgB;#FWTW*h(fq4UOE
zX!adHk^MRwxMlHe?8=f>2JmY>-~HyK(!5e<)z{D9+ZK>)+mAMS<4tZHUC8n!RbG^n
z*xe~$)6ZJkK2JgBUBi{<YZorsx%j&%BIg1!@=>nG*~VIF9N#B)Jy#aN=Wmo&?3)dz
zTw9J$u4t3)<LA`qt`bP$riXZPS0?*gtHZd)<YBa@eE-L`QR3sHNa1HlqMo@HRqV*b
zzggMQTKsw4#rpS<Nkk*OeOX^L!RsU@jl7J(e1EZvI9>X#Cj$9mX`sK>CC_Hj3Jm3a
zy}Efnq|kVQoZz_xq!(sjk-^@ckNM}YW~a44XUP*N6dw;Zo_~vvt0~4eJ970~SLFo9
zAM8)9{b3B=ZPmo8mwMrK1{~I9{2D^Ee4wEAVmRbt1XlCqDxTB9)`9KUIYLF}r~8&U
zh4*~~*{hIPOj0>6mhT(c`qAb+VDtRhz&FnX(f&bUb5kquZM-dQ!U`m!XM<!_WuD91
z1`%dB%$Z)qk0a|fyI}K&cgU_eih^l_!m#qq#`IU-*P7NA@ckvK$iZ_5kRGPzG40Fu
z@!kD6JFS+jf%093QEx+fpx@Uo#Y$vL>94$Bkob>cGj@L%><z(3^xWvKV1yaBo0D6a
z^cspjCUG9_cmm$2dk#PA+lQ`e=h{e1O#-)Q?E_0%E0OB66?j0s4n1CplN+&dFnss=
zKC%gY1jTAz!Z$Ph8C1k!(FI=4I$}3xP^bZxTUKC#700D(`T7l06gyW1On?@dZgh^?
zph&1aIDwl#=j3MFmN~uF)f0`VJOkN{n}ol-u0wy5$1`@VvE8r+sT5>M;_Gn1)at4D
zs$E)iAO6f(_0U<ca+R)Rh}v|Qm|d%|CG7h2HSe=9A6<$1>lCV{b~I|G*jv`zlFx+u
zbLG8qV$V!N!l<KZ#RBV+OLW!MU>fIrwqUgfaw#~KO8*##uClciD(pPSm<DljpH?)0
zvHO#$mRZ55lu{I8f_&yS%weH}4uDI75=mp-Bp1oM<rO1ZdNMcp@vzf4033Q11ilV`
z0m$WjiCc;jg&X+2#mny&60A&7;u%nYyvXsWET7dOGjQgf@m+ng9ez_20qmxDAmcsm
z5KdUK@L?i1R~=X~39{+ZrOqtrg?!OSrroCGFz@)e%Cys#{JuGps#v)OwdlE-Zciea
z{rs9}WAHNgr<(!1Bxg8y>)S3Y304v|4(HCu;^V9-uX+<IAZ`TWgPSvMgL=xAEamEF
z`$iEq9G^-#)}2Q$T##W~=mX5}-@i=F-ACR0sx7hUsf$h}LYYa|_A(m$IFbw?$ppUu
zIBoU;!N9FO7>9^LdN1#nIrn=|y`NnG(eBX#c%=qYH+>(i%g<xZK1tBzt~f|weG_TR
zIES~d8%s3v>ofaht3jcuHbO>LBRhH}2=CvN5#P^qd#7;+rXzU?n_ydzhop%Tgx-&7
z;_zk8XQi|%<)(fC?9#;1Td#c?yJOlycm7#chvQRFt*b?SuHA*~39uH5huPvq-?+MO
ze`E`n{D^VB@<0LZc&){JNj!ur+i+`BuM2uqO~p0%n$Zlf|6UYxX1ucSrzuz7MK4Pz
zv;C@=vY!gC8D&H?rkuhGy*ZsP(lF%o5N%0H{UUU%-&A2%v%1XVBZt91G622#1Tc3g
zLndu#B$Q07h_XNqGcX$fr+7!fs?p(Ku7*hHZZ}A_d9C@cxHP#h4ByLn4z|qbkH*g1
z;ki3;fN(n>>(tHFB<)}Bf=h2=&=(J7%tPf?;(-^Zvt7T0Jk@%OGOaa3GX0C0`fX`4
z^MTy_FlQ@N_Oxy@;%u=&aE-E&eS5WlFr3cGkr#`oDL0$QmHOk6rCSwb)7IzEbMA3i
z<h?%BM)iK+YRd^U=cNnXRQLr?<=Y){H>Q-fr~;i39c0tGaHjETG~w*Z$!&Dj6Tg}{
zyma38L}aaGScUZnXZ$HYrbiC!E7p}9f_```1!v#%ma1$I!^`<`JK4J0DaU=VvuoXI
z$bChMr{cf_ypr!rT8~?R>eiLW-HXxa+9CTr-!!M-ZSokx#_PDBcaSyP-XY~G!vq6w
zC3-eYO2J?7@nzzZA+VCqDKRlI3k1JO#(Lzb<Bsw=uyJXRMt`&{Axh#Blqy(hnu59Z
z(ZEOZ^;>Xa7*$=Aj1&rPyVNfBEmu3@iP!OUm_KO%6{yL;E=5Pc&!;nFhqIm$vWr}Q
z8L@evC?n>@Iho@q_^#b7rnT}cUC;0P6JGBD`yJH5vzpz=7V#}&+~9-6HU7EGnx$V!
z&FJqS`Is82Ut`D^oplmE-p=_mvC06}Tj&VB9N2(BZrOxZw6^f;WDe_;?iY*v&xv+b
zSwauD<l{4`2pWCDosmr&VFfdNN5CWBFM^t0xe;xniit6+I6IdPSWGoucW_j8|AIW)
z|B;y?x{PJ(a5`tN6G3w~g^+`y`-0B1eHo}B7w@Rzu;$#ORBU(@Y_c{EDZxL|TK-*_
z$#4$4{^=1MuU!F>?r%g!pK=!(9-k}>^5Djp3Tz{GeYQ+|p0t2pcd0TNbpk?@pZ}J)
zeh?e@KZkC7*pD<_P-MdCt9WmI?-aSdoJ0j8asZ(Rj7FShv^%P4rhzN(VW<x|r)vS}
zZ|jWo`}~YPEOBKP@N20ryRL#0YkEiw28<Sb{-}Wyhk7u__&uKb%`OmxAE7D>El^4A
zP9g@ZXI5EpI*D5dnY=E5yx3F%sZ1@zlxG;ow)k@6s4P=WJYipe^98T<pijaDY4G&{
zQhR=GJpXhgIWe0&_hXV6Za_QfV)aMTH2xeqVNYMUagHyz{Q3)ME8blp|1HaLz?F9~
zy%9QA7Y4Po=?Qk24#t`rPnM73`x&)ORpO)c23}oci!%PRn0vS>Uc&z_PdjQqb*oND
z4ekKQn-ghFkV7bwlg{a!64nP4e^i3A=3GVRj}E}+C}~R9@OG-fx5$(6>tW*cP}D>v
z51ZrWE`7m|zi|fk&=e=JD5W|D^tX1xs~`2C+vFHKRy>2U5wD*IsSSFU1Tl*tTr{T#
zt;4VDqMpA&?#&GZv3v7e=Ioq-UxJ;a^Z4=(A6$+^W?!a4Z0k|$$*YJ9=WA#sehs!Y
zcMLMix(_8@X@trORD@;$&2-c*t`3Tu!$C*^4s2E(MDu6ez|DRv#_#cauA#0spv@^F
z=jLt7(CKMXOq?ExL41FKkT~*9t)qCSixxE7VY@7?J(8Hw!Oe}9kEX$6@;u1q*;_!O
zks%TM=_h?an=7xiVlg#Za)9#NIUI2-)Rck|Ud**|9Hzfc6OCm`B}Yp>2t-jE@z?!p
zr0#rt*^OKmee?V5yyHqAWc^8%3g?_0Y#rZzO`nX(fDaqsWh8<2s~%fnx~&q+<}sVy
z-{9)@P|72FB05+-6LZWarGxp{)3iGq`o8)c6uj;bwBS@S(Xm=txF?zGgX0pviPNSj
zKx-R2k@0o6Dqc43#A11!zFWdbk6bZ)SH&CEYSNX?E53j|=dt^{<Eg!O{E)g`I|Vmw
zdq~eLsK=c7HcD8ZM=BVu0DBjvq4q2mcJLYYyI)WC{Q?iyE``r&X}}|sRfx@>J8&)e
zoXW;s`xaMtt(Asz!M7Elrt)aH&jt&89RC|7vl)w>L(Ys6ZJhWH99n(3;+VyG%#ttf
z=g)giFN%gb>5Ha9p{M6pEb7l-`{Wqg&!(PF0-LVe!P!pF(YvM*6}EdSu#IwzE$`^9
zDBvp2LS7CXA`rxcSFGK62eaVi3htD^cZYkzlhU<NsY!T+UFS`#FJE4dlkb4j35*Jv
zJVRjoF03N(@l8yN#|$gOVXNB&utN1+1}L2I>=?HJbK`wRZC@hJP_}dKmvRsirp^`i
z8`Q+i;P+fN*H}=kALqhR8~u>mF;|(++uN9Pkz9Y7e_0oq59+CXNk)Rpx<B#$l?n7j
z{@GdP9(~Eo13Hqu%g>^<QDt;SsTCvU$Ag{4SMprNT`2j24{`|ZqUUQ0=^Xxd>Ei}E
zQ@xdK#EVXpz!MgPGTM**Y4;;sd8<e42g<wW!#~S>Kv&34B7Im4{ydh$zRuD>;_|d5
zr|q8$CWx<N71pHGgRlE=o2xLH_<^K1t^}^;TExaF3z!M(IJwlyebBa4HL)b<8`M^Y
z<MBZ&3EjQi8s+kd*<jr4R!M`=JV6%|ff+e>(69J&sShuGAn)gqWPflGtQu-89htI)
zj^@Ya@zpy4vDO+C6?-ALjG`xfw#M(Ya&4>lMO~t@ca%t_U4$4v*JKQW)$mdBGsl0f
ze?0Jy2mU|lfwAMqvE%dif6>|h{y)gW<>|f4Pyh8_?f?3Jc{BEX9CDfe(<Ogz$KT<-
z|C>BB)~3IsZ_MB1t^S+5*uTk>|G)k3*jK{x+}t?sf8n$Fzpr2KVqkrv;KqM`1MGL0
zJdeMgu=*zOJh|TP*WqsI_bg9^{e-c3-Sg#q9l1PzDY5TM!TtuXH_J<Ad6(Q26zbSd
z_v@EEd`3T(ca7!A`%U+JIV|p0t{KbIWIu9F_k1}|hn4?*4QDwX>}Qu;wp+d&m)rgO
z8tthN7RmB{y+5W~z8s&>E$zZ`l6aonKi#j%-O^)N&KdUefGzh=p1d#pz9zDqa`t1y
z<}K-#FUQNfrDw667WQMs=5;Se&b!c!x5!iBQV;fh9@cE$f^PY8-mg9)zhSS{o(gUT
zEYE`Ftowr}_y6x}3#)Gs%ai+bU$=ZYey>|Pg7u%ASI5@(&wAwl#ryALh=SXHjv>*k
zt%j^_dEMv!rHivr+AST+`p1IhnX!3)*1fwg<5_*9cwSUDdAY_r-O|Zy-C>pov3YUb
z^5wiY-O_1nx$Z2_j?GK$mM`ah?3O;vmg~v#<n=h%E&oq_xoo*(S)Mw}3HyWBy}l>e
za(!7|Z<e#^4_^2B7P9&XmS@d!V*lVZbW4}8_KjzGa{v8lUvoDe$+oqer^4p#>Xt9B
zM^(4<1-9JjEN?WM*WGS8Pn)&<_f^4`JD25QEbnMHc{#4sjaSL?7PGwREGO{~o<cX?
zb@qNOudllw<ofiw@oL!nRV+{5ZaIJG>+F`k&(>o#%agZTX}5g2KGkmep0Ivd$nq>$
zPTn89@7>b%EN>ReQ(*H}cgvUS6L(9$WXttnd6n##lH1t*is_bq&D!9^^4_yNxv#rl
z@|gSkde7>!VtMkIG@)C*9M9~QZe#sCkhgDVw|qH|>6ZS^>g&UoE015@ukQABvGwiA
z^5k(^<qzK6ZfPY<A*`LPOA*VH_v!A}g1^%LeGGB?&oLxSh1LCz)$R7*=tlq2t)L)(
z7t$-1XUFo?yXDJm`O~-E?fd8c9}oQFf&VEUVBcf*YwyC1^={+ZO#$1}`<q<=p<;hR
zugM-7l<n}oA38clLp1tTjgjNv0UI1o_?J3j{z^{rx_YtSVm3XE&9}U5<oKTDZ8kth
z%UHPt{U-a*tez{H-|TL)I>XufI#%aQ7Sm@r$5?I$tJ98^J<008d7Wd}?{jQ=9Gh>%
z>U_rX=Cb^JRxXg$d6(6*lhwI})j5T|Ph)kCWU(gJ&PbNq$m%p>Wed9L^k#LEY}%jA
zw_<g+vAj@L=M7dagVou=>UqfO{K4u>V(+_HovT>PlI5IZxyq~z#Oow@JB93bF{^V7
zo8OnkZ?n7!EI*Bv6S6uhSUq!Co&Qt%#QzjKgI}n7wAQzHi1NBTdY;JeSgh5{!)^R^
zcE0-8%S80mCF|xCSB>?tu9<U>yB<DX=IXJ%{J-qQ{cDGLn&5Gs<8|6RmzYn$?lDDJ
zhj~AIh(;)$ob~{p_|=(EVBgmkZdc?wNH5=Y!^>RPD+@AQPYyfa3Icb#9*c@}UE&<&
zI@fd$ue*|+H|6^ucKxituGtj6dsbXy=ifiBWj|zIHE-g*G#wUrJ-OlFHGf8LumA3S
zq*omldWue#dsa>vg5k?!u|?7Eu()JzypK*ku5GSI3_~~lAKoSRN5TIwDWChh#e@IA
lW#6ZY%YR$KuF?OCPyd71KhQrO_{RhPc;J@@6501M|1T-7v6cV;

literal 0
HcmV?d00001

diff --git a/tests/test_cli.py b/tests/test_cli.py
index 8e0af46..f0bf483 100644
--- a/tests/test_cli.py
+++ b/tests/test_cli.py
@@ -5,8 +5,8 @@
 from panama import read_DAT
 import pytest
 
-CORSIKA_VERSION = "corsika-77500"
-CORSIKA_EXECUTABLE = "corsika77500Linux_SIBYLL_urqmd"
+CORSIKA_VERSION = "corsika-77550"
+CORSIKA_EXECUTABLE = "corsika77550Linux_SIBYLL_urqmd"
 
 def test_cli_missing_executable(
     test_file_path=Path(__file__).parent / "files" / "example_corsika.template",
diff --git a/tests/test_read.py b/tests/test_read.py
index aa9bb1e..7a75f6f 100644
--- a/tests/test_read.py
+++ b/tests/test_read.py
@@ -15,6 +15,7 @@
 
 SINGLE_TEST_FILE = Path(__file__).parent / "files" / "DAT000000"
 GLOB_TEST_FILE = Path(__file__).parent / "files" / "DAT*"
+NOEHIST_TEST_FILE = Path(__file__).parent / "files" / "noEHIST" / "DAT101001"
 
 
 def test_noparse(test_file_path=SINGLE_TEST_FILE):
@@ -61,6 +62,13 @@ def check_eq(file, df_run, df_event, particles, skip_mother=False):
                 num += 1
 
 
+def test_noehist(test_file_path=NOEHIST_TEST_FILE):
+    # This used to fail (see issue #103)
+    df_run, df_event, df = panama.read_DAT(glob=test_file_path, mother_columns=True)
+
+    check_eq(test_file_path, df_run, df_event, df, skip_mother=False)
+
+
 def test_noadd(test_file_path=SINGLE_TEST_FILE):
     with pytest.raises(ValueError, match="requires"):
         df_run, df_event, particles = panama.read_DAT(
diff --git a/tests/test_run.py b/tests/test_run.py
index 2ac234d..2680a8d 100644
--- a/tests/test_run.py
+++ b/tests/test_run.py
@@ -9,8 +9,8 @@
 import pytest
 import numpy as np
 
-CORSIKA_VERSION = "corsika-77500"
-CORSIKA_EXECUTABLE = "corsika77500Linux_SIBYLL_urqmd"
+CORSIKA_VERSION = "corsika-77550"
+CORSIKA_EXECUTABLE = "corsika77550Linux_SIBYLL_urqmd"
 
 
 def test_corsika_runner_cleanup(