Skip to main content
This is a pre-production deployment of Warehouse. Changes made here affect the production instance of PyPI (pypi.python.org).
Help us improve Python packaging - Donate today!

Python read-only implementation of the EventIO file format

Project Description
pyeventio |Build Status|
========================

A Python (read-only) implementation of the EventIO data format invented
by Konrad Bernloehr as used for example by the IACT extension for
CORSIKA: https://www.ikp.kit.edu/corsika

Reading the data of the following Objects is currently supported:

+--------+-------------------------------+
| Code | Description |
+========+===============================+
| 1200 | CORSIKA Run Header |
+--------+-------------------------------+
| 1201 | CORSIKA Telescope Positions |
+--------+-------------------------------+
| 1202 | CORSIKA Event Header |
+--------+-------------------------------+
| 1203 | CORSIKA Array Offsets |
+--------+-------------------------------+
| 1204 | CORSIKA Telescope Data |
+--------+-------------------------------+
| 1205 | IACT Photons |
+--------+-------------------------------+
| 1209 | CORSIKA Event End Block |
+--------+-------------------------------+
| 1210 | CORSIKA Run End Block |
+--------+-------------------------------+
| 1211 | CORSIKA Longitudinal Block |
+--------+-------------------------------+
| 1212 | CORSIKA Input Card |
+--------+-------------------------------+

The following blocks are known, but reading their data is not (yet)
implemented:

+--------+-----------------------+
| Code | Description |
+========+=======================+
| 1206 | IACT Layout |
+--------+-----------------------+
| 1207 | IACT Trigger Time |
+--------+-----------------------+
| 1208 | IACT PhotoElectrons |
+--------+-----------------------+

install with
------------

::

pip install eventio

Open a file produced by the IACT CORSIKA extension
--------------------------------------------------

First Example
~~~~~~~~~~~~~

One may iterate over an instance of ``IACTFile`` class in order to
retrieve events. Events have a small number of fields. The most
important one is the ``photon_bunches`` field, which is dictionary
mapping telescope ids to a simple structured ``np.array``, containing
the typical parameters Cherekov photons bunches in Corsika have, like:

- ``x``, ``y`` coordinate in the observation plane (in cm)
- direction cosine ``cx``, ``cy`` in x and y direction of the incident
angle of the photon
- wavelength ``lambda`` of the photon (in nm)
- number of ``photons`` associated with this bunch
- the ``time`` since the first interaction (in ns, I believe)
- the production height of the photon bunch (called ``zem``)
- a bool flag, whether the photon was scattered in the atmosphere

An event has the following attributes: \* ``header``: a dictionary
containing the Corsika Event Header data \* ``end_block``: a numpy array
containing the Corsika Event End data \* ``time_offset``, ``x_offset``,
``y_offsett``, the offset of the array

This prints the number of photons for the first telescope in every
event:

.. code:: {python}

import eventio

with eventio.IACTFile('data/telescope.dat') as f:
for event in f:
print(event.photon_buches[0]['photons'].sum())

Second Example
~~~~~~~~~~~~~~

If you like to plot the origin of the Cherenkov photons of the first
event in file ``data/telescope.dat`` for the first telescope, have a
look into `this
example <https://github.com/fact-project/pyeventio/blob/new_api/examples/plot_production_3d.py>`__

It might look similar to this picture:

.. figure:: https://raw.githubusercontent.com/fact-project/pyeventio/master/a_shower.png
:alt: an example shower

an example shower

Low level access
----------------

For more low level access to the items of an ``EventIO`` file (or to
implement a higher level abstraction like ``IACTFile``) one can use the
``EventIOFile`` class which gives access to the ``Objects`` and
``subitems`` in ``EventIO`` files.

This is how our test file looks like in the low level view:

::

In [1]: import eventio

In [2]: eventio.EventIOFile('eventio/resources/one_shower.dat')
Out[2]:
EventIOFile(path=eventio/resources/one_shower.dat, objects=[
CorsikaRunHeader(first=0, length=1096)
CorsikaInputCard(first=1112, length=448)
CorsikaTelescopeDefinition(first=1576, length=20)
CorsikaEventHeader(first=1612, length=1096)
CorsikaArrayOffsets(first=2724, length=16)
CorsikaTelescopeData(first=2756, length=6136, subitems=1)
CorsikaEventEndBlock(first=8908, length=1096)
CorsikaRunEndBlock(first=10020, length=16)
])

And this is how a ``sim_telarray`` file looks like (sim\_telarray
objects are not implemted yet):

::

In [3]: eventio.EventIOFile('gamma_test.simtel')
Out[3]:
EventIOFile(path=../../CTA/ctapipe/ctapipe-extra/datasets/gamma_test.simtel, objects=[
UnknownObject[70](first=0, length=11960, subitems=131)
UnknownObject[70](first=11976, length=1744732, subitems=21526)
UnknownObject[70](first=1756724, length=838000, subitems=11186)
UnknownObject[2000](first=2594740, length=1876)
...
UnknownObject[2010](first=50007852, length=1782080, subitems=19)
UnknownObject[2021](first=51789948, length=12)
UnknownObject[2026](first=51789976, length=3536)
UnknownObject[2010](first=51793528, length=1560656, subitems=9)
])

.. |Build Status| image:: https://travis-ci.org/fact-project/pyeventio.svg?branch=master
:target: https://travis-ci.org/fact-project/pyeventio
Release History

Release History

This version
History Node

0.3.0

History Node

0.2.3

Download Files

Download Files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

File Name & Checksum SHA256 Checksum Help Version File Type Upload Date
eventio-0.3.0.tar.gz (421.7 kB) Copy SHA256 Checksum SHA256 Source Nov 25, 2016

Supported By

WebFaction WebFaction Technical Writing Elastic Elastic Search Pingdom Pingdom Monitoring Dyn Dyn DNS Sentry Sentry Error Logging CloudAMQP CloudAMQP RabbitMQ Heroku Heroku PaaS Kabu Creative Kabu Creative UX & Design Fastly Fastly CDN DigiCert DigiCert EV Certificate Rackspace Rackspace Cloud Servers DreamHost DreamHost Log Hosting