Add acoustics monitoring

I do not know if I'm doing the right thing. I tried to upload the script "acoustic.py" and I am asked to send a merge request. Is this the correct procedure? @tgal Thank you

scripts/acoustics.py

+#!/usr/bin/env python
+# coding=utf-8
+"""
+Online Acoustic Monitoring
+
+Usage:
+    acoustics.py [options]
+    acoustics.py (-h | --help)
+
+Options:
+    -d DET_ID       Detector ID 
+    -o PLOT_DIR     The directory to save the plot 
+    -h --help       Show this screen.
+
+"""
 import numpy as np
 import matplotlib.pyplot as plt
 import matplotlib
@@ -10,7 +25,7 @@ import argparse
 
 parser = argparse.ArgumentParser(description='Online_monitoring')
 parser.add_argument('-d', dest='det_id', type=str, nargs=1, required=True,
-                    help='The detector ID (e.g. D0ARCA001)')
+                    help='The detector ID (e.g. D_ORCA005)')
 parser.add_argument('-o', dest='plot_dir', type=str, nargs=1, required=True,
                     help='Directory in which the plot is saved')
 
@@ -22,7 +37,6 @@ def diff(first, second):
     second = set(second)
     return [item for item in first if item not in second]
 
-#detid='D_ORCA005'
 #print(time.time())
 
 db = kp.db.DBManager()
@@ -31,13 +45,15 @@ table=db.run_table(detid)
 
 
 
-AB=[12,14,16] # Acoustic Beacons
-DOM=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18] # DOMs
-DU=[1,2,3,4,5] # DUs
+ACOUSTIC_BEACONS = [12, 14, 16]
+N_DOMS = 18
+DOMS = range(N_DOMS + 1)
+N_DUS = 5
+DUS = range(1, N_DUS + 1)
 
-#AB=[16] # Acoustic Beacons
-#DOM=[7] #
-#DU=[3] # DUs
+#ACOUSTIC_BEACONS=[16] # Acoustic Beacons
+#DOMS=[7] # DOMs
+#DUS=[3] # DUs
 
 check=True
 while check==True:
@@ -53,13 +69,13 @@ while check==True:
     print(now)
 
     COL=[]
-    for i in DU:
+    for du in DUS:
         DUx=[]
-        for j in AB:
-            for k in DOM:
+        for ab in ACOUSTIC_BEACONS:
+            for dom in DOMS:
                 try:
-                    macaddress = db.doms.via_omkey((i,k), detid).dom_id
-                    toas_all = sds.toashort(detid=detid, minrun=minrun, maxrun=maxrun, domid=macaddress, emitterid=j) # Prendere i dati basandosi sul tempo e non sul run????
+                    macaddress = db.doms.via_omkey((du,dom), detid).dom_id
+                    toas_all = sds.toashort(detid=detid, minrun=minrun, maxrun=maxrun, domid=macaddress, emitterid=ab) # Prendere i dati basandosi sul tempo e non sul run????
                 
                     QF_abdom=toas_all["QUALITYFACTOR"]
                     UTB_abdom=toas_all["UNIXTIMEBASE"]