.. note::
    :class: sphx-glr-download-link-note

    Click :ref:`here <sphx_glr_download_auto_examples_offline_analysis_k40summary.py>` to download the full example code
.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_offline_analysis_k40summary.py:


=======================
K40 Calibration Summary
=======================

Combine k40calib results into a single CSV file.

.. code-block:: bash

    Usage:
        k40summary.py CALIB_FOLDER
        k40summary.py (-h | --help)

    Options:
        -h --help   Show this screen.




.. code-block:: python

    from __future__ import absolute_import, print_function, division

    from glob import glob
    import os
    import re
    import pickle
    import numpy as np


    def write_header(fobj):
        """Add the header to the CSV file"""
        fobj.write("# K40 calibration results\n")
        fobj.write("det_id\trun_id\tdom_id")
        for param in ['t0', 'qe']:
            for i in range(31):
                fobj.write("\t{}_ch{}".format(param, i))


    def main():
        from docopt import docopt
        args = docopt(__doc__)

        file_pattern = os.path.join(args['CALIB_FOLDER'], '*.p')
        files = glob(file_pattern)

        with open('k40calib_summary.csv', 'w') as csv_file:
            write_header(csv_file)

            for fn in files:
                det_id, run_id = [
                    int(x) for x in re.search("_(\\d{8})" * 2, fn).groups()
                ]
                with open(fn, 'rb') as f:
                    data = pickle.load(f)
                if not data:
                    print("Empty dataset found for '{}'".format(fn))
                else:
                    for dom_id in data.keys():
                        t0s = data[dom_id]['opt_t0s'].x
                        qes = data[dom_id]['opt_qes'].x
                        cols = np.concatenate([[det_id, run_id, dom_id], t0s, qes])
                        line = '\n' + '\t'.join(str(c) for c in cols)
                        csv_file.write(line)


    if __name__ == '__main__':
        main()

**Total running time of the script:** ( 0 minutes  0.000 seconds)


.. _sphx_glr_download_auto_examples_offline_analysis_k40summary.py:


.. only :: html

 .. container:: sphx-glr-footer
    :class: sphx-glr-footer-example



  .. container:: sphx-glr-download

     :download:`Download Python source code: k40summary.py <k40summary.py>`



  .. container:: sphx-glr-download

     :download:`Download Jupyter notebook: k40summary.ipynb <k40summary.ipynb>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.readthedocs.io>`_