

.. _sphx_glr_auto_examples_monitoring_k40_calibration.py:


=========================
K40 Intra-DOM Calibration
=========================

The following script calculates the PMT time offsets using K40 coincidences




.. code-block:: python

    # Author: Jonas Reubeld <jreubelt@km3net.de> and Tamas Gal <tgal@km3net.de>
    # License: MIT
    import km3pipe as kp
    from km3pipe.logger import logging
    from km3modules import k40
    from km3modules.common import StatusBar, MemoryObserver
    import matplotlib.pylab as plt
    import numpy as np
    import pandas as pd
    import multiprocessing as mp
    from datetime import datetime
    import time
    import km3pipe.style

    db = kp.db.DBManager()


    class IntraDOMCalibrationPlotter(kp.Module):
        def process(self, blob):
            calibration = blob["IntraDOMCalibration"]
            for process in (self.create_plot, self.save_hdf5):
                proc = mp.Process(target=process, args=(calibration,))
                proc.daemon = True
                proc.start()
                proc.join()

        def create_plot(self, calibration):
            print("Creating plot...")
            fig, axes = plt.subplots(8, 4, figsize=(16, 20),
                                     sharex=True, sharey=True)
            for ax, (dom_id, calib) in zip(axes.flatten(), calibration.items()):
                ax.plot(np.cos(calib['angles']), calib["means"], '.')
                ax.plot(np.cos(calib['angles']), calib["corrected_means"], '.')
                ax.set_title("{0} - {1}".format(db.doms.via_dom_id(dom_id), dom_id))
                ax.set_ylim((-20, 20))
            plt.suptitle("{0} UTC".format(datetime.utcnow().strftime("%c")))
            plt.savefig("www/plots/intradom.png", bbox_inches='tight')
            plt.close('all')

            fig, axes = plt.subplots(8, 4, figsize=(16, 20),
                                     sharex=True, sharey=True)
            for ax, (dom_id, calib) in zip(axes.flatten(), calibration.items()):
                ax.plot(np.cos(calib['angles']), calib["rates"], '.')
                ax.plot(np.cos(calib['angles']), calib["corrected_rates"], '.')
                ax.set_title("{0} - {1}".format(db.doms.via_dom_id(dom_id), dom_id))
                ax.set_ylim((0, 15))
            plt.suptitle("{0} UTC".format(datetime.utcnow().strftime("%c")))
            plt.savefig("www/plots/angular_k40rate_distribution.png",
                        bbox_inches='tight')
            plt.close('all')

        def save_hdf5(self, calibration):
            print("Saving calibration information...")
            store = pd.HDFStore('data/k40calib.h5', mode='r')
            now = int(time.time())
            timestamps = (now,) * 31
            for dom_id, calib in calibration.items():
                tdc_channels = range(31)
                t0s = calib['opt_t0s'].x
                dom_ids = (dom_id,) * 31
                df = pd.DataFrame({
                    'timestamp': timestamps,
                    'dom_id': dom_ids,
                    'tdc_channel': tdc_channels,
                    't0s': t0s
                    })
                store.append('t0s', df, format='table', data_columns=True)
            store.close()


    pipe = kp.Pipeline()
    pipe.attach(kp.io.ch.CHPump,
                host='127.0.0.1',
                port=5553,
                tags='IO_TSL',
                timeout=7*60*60*24,
                max_queue=42)
    pipe.attach(StatusBar, every=100)
    pipe.attach(MemoryObserver, every=500)
    pipe.attach(kp.io.daq.TimesliceParser)
    pipe.attach(k40.CoincidenceFinder, accumulate=10*60*30)
    pipe.attach(k40.IntraDOMCalibrator)
    pipe.attach(IntraDOMCalibrationPlotter)
    pipe.drain()

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



.. container:: sphx-glr-footer


  .. container:: sphx-glr-download

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



  .. container:: sphx-glr-download

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

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

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