|
5 | 5 | import pandas as pd |
6 | 6 | import seaborn as sns |
7 | 7 | import tables |
| 8 | + |
| 9 | +# Imports from ctapipe |
8 | 10 | from ctapipe.coordinates import EngineeringCameraFrame |
9 | 11 | from ctapipe.instrument import CameraGeometry |
10 | 12 | from ctapipe.visualization import CameraDisplay |
| 13 | +from ctapipe_io_nectarcam.constants import N_PIXELS |
11 | 14 | from scipy import stats |
12 | 15 |
|
13 | 16 | sns.set(style="whitegrid") |
|
91 | 94 | pixel_display = [100, 144, 240, 723, 816, 1034, 1516] |
92 | 95 | fill_value = np.nan |
93 | 96 |
|
94 | | -pixel_process = np.arange(1855) |
| 97 | +pixel_process = np.arange(N_PIXELS) |
95 | 98 | # pixel_process = pixel_display |
96 | 99 |
|
97 | | -slopes_width_hg = np.full([len(pedfiles), 1855], fill_value=fill_value) |
98 | | -slopes_width_lg = np.full([len(pedfiles), 1855], fill_value=fill_value) |
| 100 | +slopes_width_hg = np.full([len(pedfiles), N_PIXELS], fill_value=fill_value) |
| 101 | +slopes_width_lg = np.full([len(pedfiles), N_PIXELS], fill_value=fill_value) |
99 | 102 |
|
100 | 103 | for itemp, uberset in enumerate(pedfiles): |
101 | 104 | temp = uberset["T"] |
|
106 | 109 | if not os.path.exists(outdir): |
107 | 110 | os.makedirs(outdir) |
108 | 111 |
|
109 | | - slopes_ped_hg = np.full(1855, fill_value=fill_value) |
110 | | - slopes_rms_hg = np.full(1855, fill_value=fill_value) |
111 | | - slopes_ped_lg = np.full(1855, fill_value=fill_value) |
112 | | - slopes_rms_lg = np.full(1855, fill_value=fill_value) |
| 112 | + slopes_ped_hg = np.full(N_PIXELS, fill_value=fill_value) |
| 113 | + slopes_rms_hg = np.full(N_PIXELS, fill_value=fill_value) |
| 114 | + slopes_ped_lg = np.full(N_PIXELS, fill_value=fill_value) |
| 115 | + slopes_rms_lg = np.full(N_PIXELS, fill_value=fill_value) |
113 | 116 |
|
114 | 117 | removed_pixels = 0 |
115 | 118 | flagged_pixels = 0 |
|
297 | 300 | if not os.path.exists(outdir): |
298 | 301 | os.makedirs(outdir) |
299 | 302 |
|
300 | | -slopes_hg = np.full(1855, fill_value=fill_value) |
301 | | -slopes_lg = np.full(1855, fill_value=fill_value) |
| 303 | +slopes_hg = np.full(N_PIXELS, fill_value=fill_value) |
| 304 | +slopes_lg = np.full(N_PIXELS, fill_value=fill_value) |
302 | 305 | x = [data["T"] for data in pedfiles] |
303 | 306 |
|
304 | 307 | print("Working on temperature variations") |
|
0 commit comments