Add IC Irene interactive monitor and config support

.gitignore

@@ -1,5 +1,6 @@
 # Byte-compiled / optimized
 __pycache__/
+*/__pycache__
 *.py[cod]
 
 # C extensions

invisible_cities/config/irene.conf

@@ -1,4 +1,4 @@
-files_in = '$ICDIR/database/test_data/electrons_40keV_z25_RWF.h5'
+files_in = '/analysis/16057/hdf5/data/ldc1/run_16057_0000_ldc1_trg0.waveforms.h5'
 
 # REPLACE /tmp with your output directory
 file_out = '/tmp/electrons_40keV_z25_PMP.h5'
@@ -7,8 +7,8 @@ file_out = '/tmp/electrons_40keV_z25_PMP.h5'
 compression = 'ZLIB4'
 
 # run number 0 is for MC
-run_number = 0
-detector_db = 'new'
+run_number = 16057
+detector_db = 'next100'
 
 # How frequently to print events
 print_mod = 1
@@ -20,35 +20,40 @@ n_baseline =   28000 # for a window of 800 mus
 
 # Set MAW for calibrated sum
 n_maw   = 100
-thr_maw =   3 * adc
+thr_maw =   3.0 * adc
 
 # Set thresholds for calibrated sum
-thr_csum_s1 = 0.5 * pes
-thr_csum_s2 = 1.0 * pes
+thr_csum_s1 = 3.0 * pes
+thr_csum_s2 = 3.5 * pes
 
 # Set parameters to search for S1
 # Notice that in MC file S1 is in t=100 mus
-s1_tmin       =  99 * mus # position of S1 in MC files at 100 mus
-s1_tmax       = 101 * mus # change tmin and tmax if S1 not at 100 mus
+s1_tmin       = 99.0 * mus # position of S1 in MC files at 100 mus
+s1_tmax       = 1500.0 * mus # change tmin and tmax if S1 not at 100 mus
 s1_stride     =   4       # minimum number of 25 ns bins in S1 searches
 s1_lmin       =   8       # 8 x 25 = 200 ns
-s1_lmax       =  20       # 20 x 25 = 500 ns
+s1_lmax       =   20       # 20 x 25 = 500 ns
 s1_rebin_stride = 1       # Do not rebin S1 by default
 
 # Set parameters to search for S2
-s2_tmin     =    101 * mus # assumes S1 at 100 mus, change if S1 not at 100 mus
-s2_tmax     =   1199 * mus # end of the window
+s2_tmin     =    1550.0 * mus # assumes S1 at 100 mus, change if S1 not at 100 mus
+s2_tmax     =   1650.0 * mus # end of the window
 s2_stride   =     40       #  40 x 25 = 1   mus
 s2_lmin     =    100       # 100 x 25 = 2.5 mus
 s2_lmax     = 100000       # maximum value of S2 width
 s2_rebin_stride = 40       # Rebin by default, 40 25 ns time bins to make one 1us time bin
 
-pmt_samp_wid  = 25 * ns
-sipm_samp_wid = 1 * mus
+# Set S2Si parameters
+thr_sipm      = 12.0 * pes
+thr_sipm_s2   = 4.0 * pes  # Threshold for the full sipm waveform
+thr_sipm_type = common
 
-cutting_function = threshold 
+pmt_samp_wid  = 25.0 * ns
+sipm_samp_wid = 1.0 * mus
+
+cutting_function = threshold
 cutting_params   = dict(  thr_sipm_type = common 
-                        , thr_sipm      = 3.5 * pes # Threshold for each SiPM time bin
-                        , thr_sipm_s2   = 10 * pes  # Threshold for the full sipm waveform
+                        , thr_sipm      = thr_sipm
+                        , thr_sipm_s2   = thr_sipm_s2
                         , detector_db   = detector_db 
                         , run_number    = run_number)