Skip to content
Snippets Groups Projects

Add acoustics monitoring

Merged Carlo Guidi requested to merge cguidi/km3mon:undefined into master
1 unresolved thread
Compare and Show latest version
1 file
+ 19
17
Compare changes
  • Side-by-side
  • Inline
+ 19
17
@@ -46,12 +46,14 @@ table=db.run_table(detid)
ACOUSTIC_BEACONS = [12, 14, 16]
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
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:
@@ -67,13 +69,13 @@ while check==True:
print(now)
COL=[]
for i in DU:
DUx=[]
for j in AB:
for k in DOM:
for du in DUS:
N_Pulses_Indicator=[] # Array indicating for each DU how many pulses each piezo reveals.
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"]
@@ -153,20 +155,20 @@ while check==True:
print(NUM)
if (NUM>7):
DUx.append(1.5)
N_Pulses_Indicator.append(1.5)
elif (NUM<8 and NUM>3):
DUx.append(0.5)
N_Pulses_Indicator.append(0.5)
elif (NUM<4 and NUM>0):
DUx.append(-0.5)
N_Pulses_Indicator.append(-0.5)
elif (NUM==0):
DUx.append(-1.5)
N_Pulses_Indicator.append(-1.5)
except:
stop=[]
DUx.append(-1.5)
N_Pulses_Indicator.append(-1.5)
COL.append(DUx)
COL.append(N_Pulses_Indicator)
Loading