~vpec/maus/tof_calib_read

4 by tunnell
moving stuff from tucs
1
# Author: Christopher Tunnell <tunnell@hep.uchicago.edu>
2
#
3
# This is a script intended to provide summary status information about the
4
# CIS system to somebody like the DQ leader.
5
#
6
# March 04, 2009
7
#
8
9
execfile('src/load.py', globals()) # don't remove this!
10
11
12
"""
13
Step 1:
14
You need to set the 'runs' variable to reflect which runs you would like to look at.  There are various examples below: the set of runs used for the Feb. 2009 TileWeek, the set of runs of Feb. 27th 2009, and the last example uses the last month of CIS runs.
15
"""
16
17
#runs = [100489,100699,100814,101152,101441,102009] # Runs as of Feb. 2009 TileWeek
18
runs = '-1 month' # Use the last month of runs
19
#runs = [101184]  # Just use one run, useful for debugging
20
#runs = [104131]
21
#runs = [103347, 104675, 104835, 104934, 102208]
22
#runs = [104467]
23
runs = '-2 week'
24
#runs = [110934]
25
        
26
"""
27
Step 2:
28
Do you want to look at the CIS scans for the last month?  If you set 'scans' equal to true, then you will have plots of the CIS scans for all channels that failed to be calibrated. However, be mindful that this will load in CIS scans for each channel for each run.  This can be slow (~30 minutes).
29
"""
30
31
useScans = True  #  Read in the CIS scan graphs in order to plot them?
32
33
"""
34
Step 3:
35
Which regions in the detector do you want to look at?
36
"""
37
38
selected_region = '' # Look at who detector
39
#selected_region = 'LBA' # Look at just LBA
40
#selected_region = 'LBA_m42_c0_lowgain' # Look at LBA module 42 channel 0 lowgain
41
#selected_region = 'lowgain' # Look at low-gain only
42
#selected_region = 'EBC_m06'
43
44
#
45
#  Experts only below.
46
#
47
48
# Define a plot of qflag and DB problems shown over whole detector
49
p = PlotModuleChannel(runType2='CIS', parameter2='isCalibrated_qflag',   text2='1% DB Variation', invert2=False,\
50
                      runType1='CIS', parameter1='isCalibrated', text1='DB Variation and QFlag Failure', invert1=True)
51
52
x = ReadCIS(getScans=useScans, getScansRMS=True)
53
54
#for partition in ['EBA', 'LBC', 'EBC', 'LBC']:
55
# This is the worker processing loop.
56
Go([\
57
    Use(run=runs, runType='CIS', region=selected_region),       # Tell Tucs which runs to use                       \
58
    x,                                                          # Tell Tucs to load in the CIS information          \
59
    ReadDB(),
60
    ScanRMSvsScan(),])