Ops changes

Author: ocerri

DQM_SiPM.py

@@ -111,6 +111,10 @@ def define_range_around_peak(h, perc = [0.4, 0.4], Range=[0.0, 99999.0]):
         if len(args.input_file) > 1:
             aux = re.search(r'Run[0-9]+', args.input_file[-1])
             flag += '_'+aux.group(0)[3:]
+    elif len(args.runs_interval)==1:
+        N = args.runs_interval[0]
+        args.input_file = [args.input_file[0].replace('XXX', str(N))]
+        flag = '_'+str(N)
     elif len(args.runs_interval)==2:
         deduced_file_list = []
         for i in range(args.runs_interval[0], args.runs_interval[1]+1):
@@ -202,7 +206,7 @@ def define_range_around_peak(h, perc = [0.4, 0.4], Range=[0.0, 99999.0]):
 
             amp_aux = np.concatenate(list(tree2array( chain, 'amp['+str(k)+']')))
 
-            h = rt.TH1D(name, title, 80, np.percentile(amp_aux, 1), min(990., np.percentile(amp_aux, 99.9)))
+            h = rt.TH1D(name, title, 80, np.percentile(amp_aux, 1), min(990., np.percentile(amp_aux, 99.99)))
             h.SetXTitle('Peak amplitude [mV]')
             h.SetYTitle('Events / {:.1f} mV'.format(h.GetBinWidth(1)))
             chain.Project(name, 'amp['+str(k)+']')
@@ -427,6 +431,8 @@ def define_range_around_peak(h, perc = [0.4, 0.4], Range=[0.0, 99999.0]):
                 canvas['space_corr'][k].Divide(2)
 
                 selection += ' && ntracks == 1 && chi2 < 8'
+                if chain.GetEntries(selection) == 0:
+                    continue
                 delta_t = np.concatenate(list(tree2array(chain, var_dT, selection)))
 
                 add_sel = ''

config/FNAL_TB_1811/VME_v3.txt

@@ -4,7 +4,7 @@
 ch_num idx_time idx_dut idx_ref var_ref amp_min amp_max
 0 0 0 -1 gaus_mean 20 240
 1 0 0 0 LP2_20 100 600
-2 0 0 0 LP2_20 100 600
+2 0 0 0 LP2_20 400 1200
 3 0 0 0 LP2_20 100 600
 4 0 0 0 LP2_20 400 1200
 5 0 0 0 LP2_20 100 800
@@ -16,5 +16,5 @@ ch_num idx_time idx_dut idx_ref var_ref amp_min amp_max
 12 2 0 -1 LP2 20 200
 13 2 0 -1 LP2 20 200
 14 2 0 -1 None 20 200
-15 2 0 -1 None 20 200
-16 2 0 -1 None 20 200
+15 2 0 -1 None 400 1200
+16 2 0 -1 None 100 850

config/FNAL_TB_1811/VME_v4.txt

@@ -0,0 +1,37 @@
+#-->Print: Amp PosRaw PosWeight TimeResRaw TimeCorrected
+-->Print: Amp PosRaw PosWeight WaveColor TimeResRaw TimeCorrected
+-->XYcenter-width: 20 20 30
+ch_num idx_time idx_dut idx_ref var_ref amp_min amp_max
+0 0 0 -1 gaus_mean 20 240
+1 0 0 0 LP2_20 300 800
+2 0 0 0 LP2_20 40 1200
+3 0 0 0 LP2_20 40 600
+4 0 0 0 LP2_20 80 1200
+5 0 0 0 LP2_20 80 800
+6 0 0 -1 LP2_20 800 1200
+7 0 0 -1 LP2_20 800 1200
+9 2 0 -1 gaus_mean 20 200
+10 2 0 -1 None 20 300
+11 2 0 -1 None 20 300
+12 2 0 -1 None 20 200
+13 2 0 -1 None 20 200
+14 2 0 -1 None 20 200
+15 2 0 -1 None 20 1200
+16 2 0 -1 None 20 850
+18 2 0 -1 gaus_mean 20 200
+19 2 0 -1 None 20 300
+20 2 0 18 LP2_20 20 300
+21 2 0 18 LP2_20 20 200
+22 2 0 18 LP2_20 20 200
+23 2 0 -1 None 20 200
+24 2 0 -1 None 20 1200
+25 2 0 18 LP2_20 20 200
+26 2 0 18 LP2_20 20 200
+27 2 0 18 LP2_20 20 200
+28 2 0 18 LP2_20 20 200
+29 2 0 18 LP2_20 20 200
+30 2 0 18 LP2_20 20 200
+31 2 0 18 LP2_20 20 200
+32 2 0 18 LP2_20 20 200
+33 2 0 18 LP2_20 20 200
+34 2 0 -1 None 200 1200

lib/histo_utilities.py

@@ -19,6 +19,8 @@ def create_TH1D(x, name='h', title=None, binning=[None, None, None], weights=Non
             bin_w = 4*(np.percentile(x,75) - np.percentile(x,25))/(len(x))**(1./3.)
             if bin_w == 0:
                 bin_w = 0.5*np.std(x)
+            if bin_w == 0:
+                bin_w = 1
             binning[0] = int((binning[2] - binning[1])/bin_w) + 5
 
         h = rt.TH1D(name, title, binning[0], binning[1], binning[2])

runDQM.sh

@@ -0,0 +1,3 @@
+source ~/.bash_profile
+cd /Users/fqnet1/FNAL_TB/18_11/DataQualityMonitor
+python DQM_SiPM.py -C config/FNAL_TB_1811/VME_v3.txt -S ~/cernbox/ocerri/www/FNAL_TB_1811/ -i ../data/VME/RECO/v3/DataVMETiming_RunXXX.root -N $1