1 | \section{Performance}
|
---|
2 |
|
---|
3 | \subsection{Calibration}
|
---|
4 |
|
---|
5 | In this section, we describe the tests performed using light pulses of different colour,
|
---|
6 | pulse shapes and intensities with the MAGIC calibration pulser box.
|
---|
7 | \par
|
---|
8 | The LED pulser system is able to provide fast light pulses of 3--4\,ns FWHM
|
---|
9 | with intensities ranging from 3--4 photo-electrons to more than 500 in one inner pixel of the
|
---|
10 | camera. These pulses can be produced in three colours $green$, $blue$ and $UV$.
|
---|
11 |
|
---|
12 | \begin{table}[htp]
|
---|
13 | \centering
|
---|
14 | \begin{tabular}{|c|c|c|c|c|c|}
|
---|
15 | \hline
|
---|
16 | \hline
|
---|
17 | \multicolumn{6}{|c|}{The possible pulsed light colours} \\
|
---|
18 | \hline
|
---|
19 | \hline
|
---|
20 | Colour & Wavelength & Spectral Width & Min. Nr. & Max. Nr. & Secondary \\
|
---|
21 | & [nm] & [nm] & Phe's & Phe's & Pulses \\
|
---|
22 | \hline
|
---|
23 | Green & 520 & 40 & 6 & 120 & yes \\
|
---|
24 | \hline
|
---|
25 | Blue & 460 & 30 & 6 & 500 & yes \\
|
---|
26 | \hline
|
---|
27 | UV & 375 & 12 & 3 & 50 & no \\
|
---|
28 | \hline
|
---|
29 | \hline
|
---|
30 | \end{tabular}
|
---|
31 | \caption{The pulser colours available from the calibration system}
|
---|
32 | \label{tab:pulsercolours}
|
---|
33 | \end{table}
|
---|
34 |
|
---|
35 | Table~\ref{tab:pulsercolours} lists the available colours and intensities and
|
---|
36 | figures~\ref{fig:pulseexample1leduv} and~\ref{fig:pulseexample23ledblue} show exemplary pulses
|
---|
37 | as registered by the FADCs.
|
---|
38 | Whereas the UV-pulse is very stable, the green and blue pulses show sometimes smaller secondary
|
---|
39 | pulses after about 10--40\,ns from the main pulse.
|
---|
40 | One can see that the very stable UV-pulses are unfortunately only available in such intensities as to
|
---|
41 | not saturate the high-gain readout channel. However, the brightest combination of light pulses easily
|
---|
42 | saturates all channels in the camera, but does not reach a saturation of the low-gain channel.
|
---|
43 | \par
|
---|
44 | Our tests can be classified into three subsections:
|
---|
45 |
|
---|
46 | \begin{enumerate}
|
---|
47 | \item Un-calibrated pixels and events: These tests measure the percentage of failures of the extractor
|
---|
48 | resulting either in a pixel declared as un-calibrated or in an event which produces a signal ouside
|
---|
49 | of the expected Gaussian distribution.
|
---|
50 | \item Number of photo-electrons: These tests measure the reconstructed numbers of photo-electrons, their
|
---|
51 | spread over the camera and the ratio of the obtained mean value for outer and inner pixels.
|
---|
52 | \item Linearity tests: These test the linearity of the extractor with respect to pulses of different intensity
|
---|
53 | and colour.
|
---|
54 | \item Time resolution: These tests show the time resolution and stability obtained with different
|
---|
55 | intensities and colours.
|
---|
56 | \end{enumerate}
|
---|
57 |
|
---|
58 | \begin{figure}[htp]
|
---|
59 | \centering
|
---|
60 | \includegraphics[height=0.25\textheight]{1LedUV_Pulse_Inner.eps}
|
---|
61 | \includegraphics[height=0.25\textheight]{1LedUV_Pulse_Outer.eps}
|
---|
62 | \caption{Example of a calibration pulse from the lowest available intensity (1\,Led UV).
|
---|
63 | The left plot shows the signal obtained in an inner pixel, the right one the signal in an outer pixel.
|
---|
64 | Note that the pulse height fluctuates much more than suggested from these pictures. Especially, a
|
---|
65 | zero-pulse is also possible.}
|
---|
66 | \label{fig:pulseexample1leduv}
|
---|
67 | \end{figure}
|
---|
68 |
|
---|
69 | \begin{figure}[htp]
|
---|
70 | \centering
|
---|
71 | \includegraphics[height=0.25\textheight]{23LedsBlue_Pulse_Inner.eps}
|
---|
72 | \includegraphics[height=0.25\textheight]{23LedsBlue_Pulse_Outer.eps}
|
---|
73 | \caption{Example of a calibration pulse from the highest available intensity (23\,Leds Blue).
|
---|
74 | The left plot shows the signal obtained in an inner pixel, the right one the signal in an outer pixel.
|
---|
75 | One the left side, the (saturated) high-gain channel is visible, on the right side, the delayed low-gain
|
---|
76 | pulse appears. Note that on the left side, there is a secondary pulses visible in the tail of the
|
---|
77 | high-gain pulse. }
|
---|
78 | \label{fig:pulseexample23ledblue}
|
---|
79 | \end{figure}
|
---|
80 |
|
---|
81 | We used data taken on the 7$^{th}$ of June, 2004 with different pulser LED combinations, each taken with
|
---|
82 | 16384 events. The corresponding run numbers range from nr. 31741 to 31772. This data was taken before the
|
---|
83 | latest camera repair access which replaced about 2\% of the pixels known to be mal-functionning at that time.
|
---|
84 | Thus, there is a lower limit to the number of un-calibrated pixels of about 1.5--2\%.
|
---|
85 | \par
|
---|
86 | Although, we had looked at and tested all colour and extractor combinations resulting from these data,
|
---|
87 | we refrain ourselves to show here only exemplary behaviour and results of extractors.
|
---|
88 | All plots, including those which are not displayed in this TDAS, can be retrieved from the following
|
---|
89 | locations:
|
---|
90 |
|
---|
91 | \begin{verbatim}
|
---|
92 | http://www.magic.ifae.es/~markus/pheplots/
|
---|
93 | http://www.magic.ifae.es/~markus/timeplots/
|
---|
94 | \end{verbatim}
|
---|
95 |
|
---|
96 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
---|
97 |
|
---|
98 | \subsubsection{Un-calibrated pixels and events}
|
---|
99 |
|
---|
100 | The MAGIC calibration software incorporates a series of checks to sort out mal-functionning pixels.
|
---|
101 | Except for the software bug searching criteria, the following exclusion reasons can apply:
|
---|
102 |
|
---|
103 | \begin{enumerate}
|
---|
104 | \item The reconstructed mean signal is less than 2.5 times the extractor resolution $R$ from zero.
|
---|
105 | (2.5 Pedestal RMS in the case of the simple fixed window extractors). This criterium cuts out
|
---|
106 | dead pixels.
|
---|
107 | \item The reconstructed mean signal error is smaller than its value. This criterium cuts out
|
---|
108 | signal distributions which fluctuate so much that their RMS is bigger than its mean value. This
|
---|
109 | criterium cuts out ``ringing'' pixels or mal-functionning extractors.
|
---|
110 | \item The reconstructed mean number of photo-electrons lies 4.5 sigma outside
|
---|
111 | the distribution of photo-electrons obtained with the inner or outer pixels in the camera.
|
---|
112 | \item All reconstructed negative mean signal, signal sigma's and mean numbers of photo-electrons
|
---|
113 | smaller than one.
|
---|
114 | \end{enumerate}
|
---|
115 |
|
---|
116 | Moreover, the number of events are counted which have been reconstructed outside a 5 sigma region
|
---|
117 | from the mean signal. These events are called ``outliers''. Figure~\ref{fig:outlier} shows a typical
|
---|
118 | outlier obtained with the digital filter.
|
---|
119 |
|
---|
120 | \begin{figure}[htp]
|
---|
121 | \centering
|
---|
122 | \includegraphics[width=0.95\linewidth]{Outlier.eps}
|
---|
123 | \caption{Example of an event classified as ``un-calibrated''. The histogram has been obtained
|
---|
124 | using the digital filter (extractor \#32) applied to a high-intensity blue pulse (run 31772).
|
---|
125 | The event marked as ``outlier'' clearly has been mis-reconstructed. It lies outside the 5 sigma
|
---|
126 | region from the fitted mean.}
|
---|
127 | \label{fig:outlier}
|
---|
128 | \end{figure}
|
---|
129 |
|
---|
130 | The following figures~\ref{fig:unsuited:5ledsuv},~\ref{fig:unsuited:1leduv},~\ref{fig:unsuited:2ledsgreen}
|
---|
131 | and~\ref{fig:unsuited:23ledsblue} show the resulting numbers of un-calibrated pixels and events for
|
---|
132 | different colours and intensities.
|
---|
133 |
|
---|
134 | \par
|
---|
135 |
|
---|
136 | \begin{figure}[htp]
|
---|
137 | \centering
|
---|
138 | \includegraphics[height=0.95\textheight]{UnsuitVsExtractor-5LedsUV-Colour-13.eps}
|
---|
139 | \caption{Uncalibrated pixels and pixels outside of the Gaussian distribution for a typical calibration
|
---|
140 | pulse of UV-light which does not saturate the high-gain readout.}
|
---|
141 | \label{fig:unsuited:5ledsuv}
|
---|
142 | \end{figure}
|
---|
143 |
|
---|
144 | \begin{figure}[htp]
|
---|
145 | \centering
|
---|
146 | \includegraphics[height=0.95\textheight]{UnsuitVsExtractor-1LedUV-Colour-04.eps}
|
---|
147 | \caption{Uncalibrated pixels and pixels outside of the Gaussian distribution for a very low
|
---|
148 | intensity pulse.}
|
---|
149 | \label{fig:unsuited:1leduv}
|
---|
150 | \end{figure}
|
---|
151 |
|
---|
152 | \begin{figure}[htp]
|
---|
153 | \centering
|
---|
154 | \includegraphics[height=0.95\textheight]{UnsuitVsExtractor-2LedsGreen-Colour-02.eps}
|
---|
155 | \caption{Uncalibrated pixels and pixels outside of the Gaussian distribution for a typical green pulse.}
|
---|
156 | \label{fig:unsuited:2ledsgreen}
|
---|
157 | \end{figure}
|
---|
158 |
|
---|
159 | \begin{figure}[htp]
|
---|
160 | \centering
|
---|
161 | \includegraphics[height=0.95\textheight]{UnsuitVsExtractor-23LedsBlue-Colour-00.eps}
|
---|
162 | \caption{Uncalibrated pixels and pixels outside of the Gaussian distribution for a high-intensity blue pulse.}
|
---|
163 | \label{fig:unsuited:23ledsblue}
|
---|
164 | \end{figure}
|
---|
165 |
|
---|
166 | One can see that in general, big extraction windows raise the
|
---|
167 | number of un-calibrated pixels and are thus less stable. Especially for the very low-intensity
|
---|
168 | $1LedUV$-pulse, the big extraction windows summing 8 or more slices, cannot calibrate more than 50\%
|
---|
169 | of the inner pixels (fig.~\ref{fig:unsuited:1leduv}). This is an expected behavior since big windows
|
---|
170 | add up more noise which in turn makes the for the small signal more difficult.
|
---|
171 | \par
|
---|
172 | In general, one can also say that all ``sliding window''-algorithms (extractors \#17-32) discard
|
---|
173 | less pixels than the ``fixed window''-ones (extractors \#1--16). The digital filter with
|
---|
174 | the correct weights (extractor \#32) discards the least number of pixels, but is also robust against
|
---|
175 | slight modifications of its weights (extractors \#28--31). Also the ``spline'' algorithms on small
|
---|
176 | windows (extractors \#23--25) discard less pixels than the previous extractors, although slightly more
|
---|
177 | then the digital filter.
|
---|
178 | \par
|
---|
179 | In the low-gain, there is one extractor discarding a too high amount of events which is the
|
---|
180 | MExtractFixedWindowPeakSearch. The reason becomes clear when one keeps in mind that this extractor
|
---|
181 | defines its extraction window by searching for the highest signal found in a sliding peak search window
|
---|
182 | looping only over {\textit non-saturating pixels}. In the case of an intense calibration pulse, only
|
---|
183 | the dead pixels match this requirement and define thus an alleatory window fluctuating like the noise
|
---|
184 | does in these channels. It is clear that one cannot use this extractor for the intense calibration pulses.
|
---|
185 | \par
|
---|
186 | It seems also that the spline algorithm extracting the amplitude of the signal produces an over-proportional
|
---|
187 | number of excluded pixels in the low-gain. The same, however in a less significant manner, holds for
|
---|
188 | the digital filter with high-low-gain inverted weights. The limit of stability with respect to
|
---|
189 | changes in the pulse form seems to be reached, there.
|
---|
190 | \par
|
---|
191 | Concerning the numbers of outliers, one can conclude that in general, the numbers are very low never exceeding
|
---|
192 | 0.25\%. There seems to be the opposite trend of larger windows producing less
|
---|
193 | outliers. However, one has to take into account that already more ``unsuited'' pixels have
|
---|
194 | been excluded thus cleaning up the sample somewhat. It seems that the ``digital filter'' and a
|
---|
195 | medium-sized ``spline'' (extractors \#25--26) yield the best result except for the outer pixels
|
---|
196 | in fig~\ref{fig:unsuited:5ledsuv} where the digital filter produces a worse result than the rest
|
---|
197 | of the extractors.
|
---|
198 | \par
|
---|
199 | In conclusion, one can say that this test excludes all extractors with too big window sizes because
|
---|
200 | they are not able to extract small signals produced by about 4 photo-electrons. The excluded extractors
|
---|
201 | are:
|
---|
202 | \begin{itemize}
|
---|
203 | \item: MExtractFixedWindow Nr. 3--5
|
---|
204 | \item: MExtractFixedWindowSpline Nr. 6--11
|
---|
205 | \item: MExtractFixedWindowPeakSearch Nr. 14--16
|
---|
206 | \item: MExtractTimeAndChargeSlidingWindow Nr. 21--22
|
---|
207 | \item: MExtractTimeAndChargeSpline Nr. 27
|
---|
208 | \end{itemize}
|
---|
209 |
|
---|
210 | The best extractors after this test are:
|
---|
211 | \begin{itemize}
|
---|
212 | \item: MExtractFixedWindow Nr. 1--2
|
---|
213 | \item: MExtractFixedWindowPeakSearch Nr. 13
|
---|
214 | \item: MExtractTimeAndChargeSlidingWindow Nr. 17--19
|
---|
215 | \item: MExtractTimeAndChargeSpline Nr. 24--25
|
---|
216 | \item: MExtractTimeAndChargeDigitalFilter Nr. 28--32
|
---|
217 | \end{itemize}
|
---|
218 |
|
---|
219 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
---|
220 |
|
---|
221 | \subsubsection{Number of photo-electrons}
|
---|
222 |
|
---|
223 | Assuming that the readout chain is clean and adds only negligible noise with respect to the one
|
---|
224 | introduced by the photo-multiplier itself, one can make the assumption that variance of the
|
---|
225 | true (non-extracted) signal $ST$ is the amplified Poisson variance on the number of photo-electrons,
|
---|
226 | multiplied with the excess noise of the photo-multiplier, characterized by the excess-noise factor $F$.
|
---|
227 |
|
---|
228 | \begin{equation}
|
---|
229 | Var(ST) = F^2 \cdot Var(N_{phe}) \cdot \frac{<ST>^2}{<N_{phe}>^2}
|
---|
230 | \label{eq:excessnoise}
|
---|
231 | \end{equation}
|
---|
232 |
|
---|
233 | After introducing the effect of the night-sky background (eq.~\ref{eq:rmssubtraction})
|
---|
234 | in formula~\ref{eq:excessnoise} and assuming that the number of photo-electrons per event follows a
|
---|
235 | Poisson distribution, one can
|
---|
236 | get an expression to retrieve the mean number of photo-electrons impinging on the pixel from the
|
---|
237 | mean extracted signal $<SE>$, its variance $Var(SE)$ and the RMS of the extracted signal obtained from
|
---|
238 | pure pedestal runs $R$ (see section~\ref{sec:determiner}):
|
---|
239 |
|
---|
240 | \begin{equation}
|
---|
241 | <N_{phe}> \approx F^2 \cdot \frac{Var(SE) - R^2}{<SE>^2}
|
---|
242 | \label{eq:pheffactor}
|
---|
243 | \end{equation}
|
---|
244 |
|
---|
245 | Equation~\ref{eq:pheffactor} must not depend on the extractor! Effectively, we will use it to test the
|
---|
246 | quality of our extractors by requiring that a valid extractor yields the same number of photo-electrons
|
---|
247 | for all pixels of a same type and does not deviate from the number obtained with other extractors.
|
---|
248 | As the camera is flat-fielded, but the number of photo-electrons impinging on an inner and an outer pixel is
|
---|
249 | different, we also use the ratio of the mean numbers of photo-electrons from the outer pixels to the one
|
---|
250 | obtained from the inner pixels as a test variable. In the ideal case, it should always yield its central
|
---|
251 | value of about 2.4--2.8.
|
---|
252 | \par
|
---|
253 | In our case, there is an additional complication due to the fact that the green and blue coloured pulses
|
---|
254 | show secondary pulses which destroy the Poisson behaviour of the number of photo-electrons. We will thus
|
---|
255 | have to split our sample of extractors into those being affected by the secondary pulses and those without
|
---|
256 | showing any effect.
|
---|
257 | \par
|
---|
258 | Figures~\ref{fig:phe:5ledsuv},~\ref{fig:phe:1leduv},~\ref{fig:phe:23ledsblue}~and~\ref{fig:phe:2ledsgreen} show
|
---|
259 | some of the obtained results. Although one can see an amazing stability for the standard 5Leds UV pulse,
|
---|
260 | there is a considerable difference for all shown non-standard pulses. Especially the pulses from green
|
---|
261 | and blue LEDs
|
---|
262 | show a clear dependency on the extraction window of the number of photo-electrons. Only the largest
|
---|
263 | extraction windows seem to catch the entire range of (jittering) secondary pulses and get also the ratio
|
---|
264 | of outer vs. inner pixels right.
|
---|
265 | \par
|
---|
266 | The strongest discrepancy is observed in the low-gain extraction (fig.~\ref{fig:phe:23ledsblue}) where all
|
---|
267 | fixed window extractors
|
---|
268 |
|
---|
269 |
|
---|
270 | \begin{figure}[htp]
|
---|
271 | \centering
|
---|
272 | \includegraphics[height=0.92\textheight]{PheVsExtractor-5LedsUV-Colour-13.eps}
|
---|
273 | \caption{Number of photo-electrons from a typical, not saturating calibration pulse of colour UV,
|
---|
274 | reconstructed with each of the tested signal extractors.
|
---|
275 | The first plots shows the number of photo-electrons obtained for the inner pixels, the second one
|
---|
276 | for the outer pixels and the third shows the ratio of the mean number of photo-electrons for the
|
---|
277 | outer pixels divided by the mean number of photo-electrons for the inner pixels. Points
|
---|
278 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
279 | \label{fig:phe:5ledsuv}
|
---|
280 | \end{figure}
|
---|
281 |
|
---|
282 | \begin{figure}[htp]
|
---|
283 | \centering
|
---|
284 | \includegraphics[height=0.92\textheight]{PheVsExtractor-1LedUV-Colour-04.eps}
|
---|
285 | \caption{Number of photo-electrons from a typical, very low-intensity calibration pulse of colour UV,
|
---|
286 | reconstructed with each of the tested signal extractors.
|
---|
287 | The first plots shows the number of photo-electrons obtained for the inner pixels, the second one
|
---|
288 | for the outer pixels and the third shows the ratio of the mean number of photo-electrons for the
|
---|
289 | outer pixels divided by the mean number of photo-electrons for the inner pixels. Points
|
---|
290 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
291 | \label{fig:phe:1leduv}
|
---|
292 | \end{figure}
|
---|
293 |
|
---|
294 | \begin{figure}[htp]
|
---|
295 | \centering
|
---|
296 | \includegraphics[height=0.92\textheight]{PheVsExtractor-2LedsGreen-Colour-02.eps}
|
---|
297 | \caption{Number of photo-electrons from a typical, not saturating calibration pulse of colour green,
|
---|
298 | reconstructed with each of the tested signal extractors.
|
---|
299 | The first plots shows the number of photo-electrons obtained for the inner pixels, the second one
|
---|
300 | for the outer pixels and the third shows the ratio of the mean number of photo-electrons for the
|
---|
301 | outer pixels divided by the mean number of photo-electrons for the inner pixels. Points
|
---|
302 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
303 | \label{fig:phe:2ledsgreen}
|
---|
304 | \end{figure}
|
---|
305 |
|
---|
306 |
|
---|
307 | \begin{figure}[htp]
|
---|
308 | \centering
|
---|
309 | \includegraphics[height=0.92\textheight]{PheVsExtractor-23LedsBlue-Colour-00.eps}
|
---|
310 | \caption{Number of photo-electrons from a typical, high-gain saturating calibration pulse of colour blue,
|
---|
311 | reconstructed with each of the tested signal extractors.
|
---|
312 | The first plots shows the number of photo-electrons obtained for the inner pixels, the second one
|
---|
313 | for the outer pixels and the third shows the ratio of the mean number of photo-electrons for the
|
---|
314 | outer pixels divided by the mean number of photo-electrons for the inner pixels. Points
|
---|
315 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
316 | \label{fig:phe:23ledsblue}
|
---|
317 | \end{figure}
|
---|
318 |
|
---|
319 | One can see that all extractor using a large window belong to the class of extractors being affected
|
---|
320 | by the secondary pulses. The only exception to this rule is the digital filter which - despite of its
|
---|
321 | 6 slices extraction window - seems to filter out all the secondary pulses.
|
---|
322 | \par
|
---|
323 | Moreover, one can see in fig.~\ref{fig:phe:1leduv} that all peak searching extractors show the influence of
|
---|
324 | the bias at low numbers of photo-electrons.
|
---|
325 | \par
|
---|
326 | The extractor MExtractFixedWindowPeakSearch at low extraction windows apparently yields chronically low
|
---|
327 | numbers of photo-electrons. This is due to the fact that the decision to fix the extraction window is
|
---|
328 | made sometimes by an inner pixel and sometimes by an outer one since the camera is flat-fielded and the
|
---|
329 | pixel carrying the largest non-saturated peak-search window is more or found by a random signal
|
---|
330 | fluctuation. However, inner and outer pixels have a systematic offset of about 0.5 to 1 FADC slices.
|
---|
331 | Thus, the extraction fluctuates artificially for one given channel which results in a systematically
|
---|
332 | large variance and thus in a systematically low reconstructed number of photo-electrons. This test thus
|
---|
333 | excludes the extractors \#11--13.
|
---|
334 | \par
|
---|
335 | Moreover, one can see that the extractors applying a small fixed window do not get the ratio of
|
---|
336 | photo-electrons from outer to inner pixels correctly for the green and blue pulses.
|
---|
337 | \par
|
---|
338 | The extractor MExtractTimeAndChargeDigitalFilter seems to be stable against modifications in the
|
---|
339 | exact form of the weights in the high-gain readout channel since all applied weights yield about
|
---|
340 | the same number of photo-electrons and the same ratio of outer vs. inner pixels. This statement does not
|
---|
341 | hold any more for the low-gain, as can be seen in figure~\ref{fig:phe:23ledsblue}. There, the application
|
---|
342 | of high-gain weights to the low-gain signal (extractors \#30--31) produces a too low number of photo-electrons
|
---|
343 | and also a too low ratio of outer per inner pixels.
|
---|
344 | \par
|
---|
345 | All sliding window and spline algorithms yield a stable ratio of outer vs. inner pixels in the low-gain,
|
---|
346 | however the effect of raising the number of photo-electrons with the extraction window is very pronounced.
|
---|
347 | Note that in figure~\ref{fig:phe:23ledsblue}, the number of photo-electrons raises by about a factor 1.4,
|
---|
348 | which is slightly higher than in the case of the high-gain channel (figure~\ref{fig:phe:2ledsgreen}).
|
---|
349 | \par
|
---|
350 | Concluding, there is now fixed window extractor yielding the correct number of photo-electrons
|
---|
351 | for the low-gain, except for the largest extraction window of 10 low-gain slices.
|
---|
352 | Either the number of photo-electrons itself is wrong or the ratio of outer vs. inner pixels is
|
---|
353 | not correct. All sliding window algorithms seem to reproduce the correct numbers if one takes into
|
---|
354 | account the after-pulse behaviour of the light pulser itself. The digital filter seems to be not
|
---|
355 | stable against exchanging the pulse form to match the slimmer high-gain pulses, though.
|
---|
356 |
|
---|
357 |
|
---|
358 | \subsubsection{Linearity tests}
|
---|
359 |
|
---|
360 | In this section, we test the lineary of the extractors. As the photo-multiplier is a linear device over a
|
---|
361 | wide dynamic range, the number of photo-electrons per charge has to remain constant over the tested
|
---|
362 | linearity region. We will show here only examples of extractors which were not already excluded in the
|
---|
363 | previous section.
|
---|
364 | \par
|
---|
365 | A first test concerns the stability of the conversion factor photo-electrons per FADC counts over the
|
---|
366 | tested intensity region.
|
---|
367 |
|
---|
368 |
|
---|
369 | \begin{figure}[htp]
|
---|
370 | \centering
|
---|
371 | \includegraphics[width=0.95\linewidth]{PheVsCharge-3.eps}
|
---|
372 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
373 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
374 | A fixed window extractor on a window size of 6 high-gain and 6 low-gain slices has been used (extractor \#3). }
|
---|
375 | \label{fig:linear:phevscharge3}
|
---|
376 | \end{figure}
|
---|
377 |
|
---|
378 | \begin{figure}[htp]
|
---|
379 | \centering
|
---|
380 | \includegraphics[width=0.95\linewidth]{PheVsCharge-8.eps}
|
---|
381 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
382 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
383 | A fixed window spline extractor on a window size of 6 high-gain and 6 low-gain slices has been used
|
---|
384 | (extractor \#8). }
|
---|
385 | \label{fig:linear:phevscharge8}
|
---|
386 | \end{figure}
|
---|
387 |
|
---|
388 | \begin{figure}[htp]
|
---|
389 | \centering
|
---|
390 | \includegraphics[width=0.95\linewidth]{PheVsCharge-14.eps}
|
---|
391 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
392 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
393 | A fixed window peak search extractor on a window size of 6 high-gain and 6 low-gain slices has been used
|
---|
394 | (extractor \#14). }
|
---|
395 | \label{fig:linear:phevscharge14}
|
---|
396 | \end{figure}
|
---|
397 |
|
---|
398 | \begin{figure}[htp]
|
---|
399 | \centering
|
---|
400 | \includegraphics[width=0.95\linewidth]{PheVsCharge-20.eps}
|
---|
401 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
402 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
403 | A sliding window extractor on a window size of 6 high-gain and 6 low-gain slices has been used
|
---|
404 | (extractor \#20). }
|
---|
405 | \label{fig:linear:phevscharge20}
|
---|
406 | \end{figure}
|
---|
407 |
|
---|
408 | \begin{figure}[htp]
|
---|
409 | \centering
|
---|
410 | \includegraphics[width=0.95\linewidth]{PheVsCharge-25.eps}
|
---|
411 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
412 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
413 | An integrating spline extractor on a sliding window and a window size of 2 high-gain and 3 low-gain slices
|
---|
414 | has been used (extractor \#25). }
|
---|
415 | \label{fig:linear:phevscharge25}
|
---|
416 | \end{figure}
|
---|
417 |
|
---|
418 | \begin{figure}[htp]
|
---|
419 | \centering
|
---|
420 | \includegraphics[width=0.95\linewidth]{PheVsCharge-27.eps}
|
---|
421 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
422 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
423 | An integrating spline extractor on a sliding window and a window size of 6 high-gain and 7 low-gain slices
|
---|
424 | has been used (extractor \#27). }
|
---|
425 | \label{fig:linear:phevscharge27}
|
---|
426 | \end{figure}
|
---|
427 |
|
---|
428 | \begin{figure}[htp]
|
---|
429 | \centering
|
---|
430 | \includegraphics[width=0.95\linewidth]{PheVsCharge-32.eps}
|
---|
431 | \caption{Example of a the development of the conversion factor FADC counts to photo-electrons for two
|
---|
432 | exemplary inner pixels (upper plots) and two exemplary outer ones (lower plots).
|
---|
433 | A digital filter extractor on a window size of 6 high-gain and 6 low-gain slices has been used
|
---|
434 | (extractor \#32). }
|
---|
435 | \label{fig:linear:phevscharge32}
|
---|
436 | \end{figure}
|
---|
437 |
|
---|
438 |
|
---|
439 |
|
---|
440 | \subsubsection{Time resolution}
|
---|
441 |
|
---|
442 | The extractors \#17--32 are able to extract also the arrival time of each pulse. In the calibration,
|
---|
443 | we have a fast-rising pulse, uniform over camera also in time. We estimate the time-uniformity to better
|
---|
444 | than 300\,ps, a limit due to the different travel times of the light between inner and outer parts of the
|
---|
445 | camera. Since the calibraion does not have an absolute measurement of the arrival time, we measure
|
---|
446 | the relative arrival time, i.e.
|
---|
447 |
|
---|
448 | \begin{equation}
|
---|
449 | \delta t_i = t_i - t_1
|
---|
450 | \end{equation}
|
---|
451 |
|
---|
452 | where $t_i$ denotes the reconstructed arrival time of pixel number $i$ and $t_1$ the reconstructed
|
---|
453 | arrival time of pixel number 1 (software numbering). For one calibration run, one can then fill
|
---|
454 | histograms of $\delta t_i$ for each pixel which yields then a mean $<\delta t_i>$, comparable to
|
---|
455 | systematic offsets in the signal delay and a sigma $\sigma(\delta t_i)$ which is a measure of the
|
---|
456 | combined time resolutions of pixel $i$ and pixel 1. Assuming that the PMTs and readout channels are
|
---|
457 | of a same kind, we obtain an approximate absolute time resolution of pixel $i$ by:
|
---|
458 |
|
---|
459 | \begin{equation}
|
---|
460 | tres_i \approx \sigma(\delta t_i)/sqrt(2)
|
---|
461 | \end{equation}
|
---|
462 |
|
---|
463 | Figures~\ref{fig:reltimesinner10leduv} and~\ref{fig:reltimesouter10leduv} show distributions of $<\delta t_i>$
|
---|
464 | for
|
---|
465 | one typical inner pixel and one typical outer pixel and a non-saturating calibration pulse of UV-light,
|
---|
466 | obtained with three different extractors. One can see that the first two yield a Gaussian distribution
|
---|
467 | to a good approximation, whereas the third extractor shows a three-peak structure and cannot be fitted.
|
---|
468 | We discarded that particular extractor for this reason.
|
---|
469 |
|
---|
470 | \begin{figure}[htp]
|
---|
471 | \centering
|
---|
472 | \includegraphics[width=0.3\linewidth]{RelArrTime_Pixel97_10LedUV_Extractor32.eps}
|
---|
473 | \includegraphics[width=0.32\linewidth]{RelArrTime_Pixel97_10LedUV_Extractor23.eps}
|
---|
474 | \includegraphics[width=0.32\linewidth]{RelArrTime_Pixel97_10LedUV_Extractor17.eps}
|
---|
475 | \caption{Example of a two distributions of relative arrival times of an inner pixel with respect to
|
---|
476 | the arrival time of the reference pixel Nr. 1. The left plot shows the result using the digital filter
|
---|
477 | (extractor \#32), the central plot shows the result obtained with the half-maximum of the spline and the
|
---|
478 | right plot the result of the sliding window with a window size of 2 FADC slices (extractor \#17). A
|
---|
479 | medium sized UV-pulse (10Leds UV) has been used which does not saturate the high-gain readout channel.}
|
---|
480 | \label{fig:reltimesinner10leduv}
|
---|
481 | \end{figure}
|
---|
482 |
|
---|
483 | \begin{figure}[htp]
|
---|
484 | \centering
|
---|
485 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel400_10LedUV_Extractor32.eps}
|
---|
486 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel400_10LedUV_Extractor23.eps}
|
---|
487 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel400_10LedUV_Extractor17.eps}
|
---|
488 | \caption{Example of a two distributions of relative arrival times of an outer pixel with respect to
|
---|
489 | the arrival time of the reference pixel Nr. 1. The left plot shows the result using the digital filter
|
---|
490 | (extractor \#32), the central plot shows the result obtained with the half-maximum of the spline and the
|
---|
491 | right plot the result of the sliding window with a window size of 2 FADC slices (extractor \#17). A
|
---|
492 | medium sized UV-pulse (10Leds UV) has been used which does not saturate the high-gain readout channel.}
|
---|
493 | \label{fig:reltimesouter10leduv}
|
---|
494 | \end{figure}
|
---|
495 |
|
---|
496 | Figures~\ref{fig:reltimesinner10ledsblue} and~\ref{fig:reltimesouter10ledsblue} show distributions of
|
---|
497 | $<\delta t_i>$ for
|
---|
498 | one typical inner and one typical outer pixel and a high-gain-saturating calibration pulse of blue-light,
|
---|
499 | obtained with two different extractors. One can see that the first (extractor \#23) yields a Gaussian
|
---|
500 | distribution to a good approximation, whereas the second (extractor \#32) shows a two-peak structure
|
---|
501 | and cannot be fitted.
|
---|
502 | \par
|
---|
503 | \ldots {\it Unfortunately, this happens for all digital filter extractors in the low-gain.
|
---|
504 | The reason is not yet understood, and has to be found by Hendrik... } \ldots
|
---|
505 | \par
|
---|
506 |
|
---|
507 | \begin{figure}[htp]
|
---|
508 | \centering
|
---|
509 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel97_10LedBlue_Extractor23.eps}
|
---|
510 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel97_10LedBlue_Extractor32.eps}
|
---|
511 | \caption{Example of a two distributions of relative arrival times of an inner pixel with respect to
|
---|
512 | the arrival time of the reference pixel Nr. 1. The left plot shows the result using the half-maximum of the spline (extractor \#23), the right plot shows the result obtained with the digital filter
|
---|
513 | (extractor \#32). A
|
---|
514 | medium sized Blue-pulse (10Leds Blue) has been used which saturates the high-gain readout channel.}
|
---|
515 | \label{fig:reltimesinner10ledsblue}
|
---|
516 | \end{figure}
|
---|
517 |
|
---|
518 |
|
---|
519 |
|
---|
520 | \begin{figure}[htp]
|
---|
521 | \centering
|
---|
522 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel400_10LedBlue_Extractor23.eps}
|
---|
523 | \includegraphics[width=0.31\linewidth]{RelArrTime_Pixel400_10LedBlue_Extractor32.eps}
|
---|
524 | \caption{Example of a two distributions of relative arrival times of an outer pixel with respect to
|
---|
525 | the arrival time of the reference pixel Nr. 1. The left plot shows the result using the half-maximum of the spline (extractor \#23), the right plot shows the result obtained with the digital filter
|
---|
526 | (extractor \#32). A
|
---|
527 | medium sized Blue-pulse (10Leds Blue) has been used which saturates the high-gain readout channel.}
|
---|
528 | \label{fig:reltimesouter10ledsblue}
|
---|
529 | \end{figure}
|
---|
530 |
|
---|
531 | \begin{figure}[htp]
|
---|
532 | \centering
|
---|
533 | \includegraphics[width=0.95\linewidth]{TimeResExtractor-5LedsUV-Colour-12.eps}
|
---|
534 | \caption{Reconstructed arrival time resolutions from a typical, not saturating calibration pulse
|
---|
535 | of colour UV, reconstructed with each of the tested arrival time extractors.
|
---|
536 | The first plots shows the time resolutions obtained for the inner pixels, the second one
|
---|
537 | for the outer pixels. Points
|
---|
538 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
539 | \label{fig:time:5ledsuv}
|
---|
540 | \end{figure}
|
---|
541 |
|
---|
542 | \begin{figure}[htp]
|
---|
543 | \centering
|
---|
544 | \includegraphics[width=0.95\linewidth]{TimeResExtractor-1LedUV-Colour-04.eps}
|
---|
545 | \caption{Reconstructed arrival time resolutions from the lowest intensity calibration pulse
|
---|
546 | of colour UV (carrying a mean number of 4 photo-electrons),
|
---|
547 | reconstructed with each of the tested arrival time extractors.
|
---|
548 | The first plots shows the time resolutions obtained for the inner pixels, the second one
|
---|
549 | for the outer pixels. Points
|
---|
550 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
551 | \label{fig:time:1leduv}
|
---|
552 | \end{figure}
|
---|
553 |
|
---|
554 | \begin{figure}[htp]
|
---|
555 | \centering
|
---|
556 | \includegraphics[width=0.95\linewidth]{TimeResExtractor-2LedsGreen-Colour-02.eps}
|
---|
557 | \caption{Reconstructed arrival time resolutions from a typical, not saturating calibration pulse
|
---|
558 | of colour Green, reconstructed with each of the tested arrival time extractors.
|
---|
559 | The first plots shows the time resolutions obtained for the inner pixels, the second one
|
---|
560 | for the outer pixels. Points
|
---|
561 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
562 | \label{fig:time:2ledsgreen}
|
---|
563 | \end{figure}
|
---|
564 |
|
---|
565 | \begin{figure}[htp]
|
---|
566 | \centering
|
---|
567 | \includegraphics[width=0.95\linewidth]{TimeResExtractor-23LedsBlue-Colour-00.eps}
|
---|
568 | \caption{Reconstructed arrival time resolutions from the highest intensity calibration pulse
|
---|
569 | of colour blue, reconstructed with each of the tested arrival time extractors.
|
---|
570 | The first plots shows the time resolutions obtained for the inner pixels, the second one
|
---|
571 | for the outer pixels. Points
|
---|
572 | denote the mean of all not-excluded pixels, the error bars their RMS.}
|
---|
573 | \label{fig:time:23ledsblue}
|
---|
574 | \end{figure}
|
---|
575 |
|
---|
576 |
|
---|
577 | \begin{figure}[htp]
|
---|
578 | \centering
|
---|
579 | \includegraphics[width=0.95\linewidth]{TimeResVsCharge-Area-21.eps}
|
---|
580 | \caption{Reconstructed mean arrival time resolutions as a function of the extracted mean number of
|
---|
581 | photo-electrons for the weighted sliding window with a window size of 8 FADC slices (extractor \#21).
|
---|
582 | Error bars denote the
|
---|
583 | spread (RMS) of the time resolutions over the investigated channels.
|
---|
584 | The marker colours show the applied
|
---|
585 | pulser colour, except for the last (green) point where all three colours were used.}
|
---|
586 | \label{fig:time:dep20}
|
---|
587 | \end{figure}
|
---|
588 |
|
---|
589 | \begin{figure}[htp]
|
---|
590 | \centering
|
---|
591 | \includegraphics[width=0.95\linewidth]{TimeResVsCharge-Area-23.eps}
|
---|
592 | \caption{Reconstructed mean arrival time resolutions as a function of the extracted mean number of
|
---|
593 | photo-electrons for the half-maximum searching spline (extractor \#23). Error bars denote the
|
---|
594 | spread (RMS) of the time resolutions over the investigated channels.
|
---|
595 | The marker colours show the applied
|
---|
596 | pulser colour, except for the last (green) point where all three colours were used.}
|
---|
597 | \label{fig:time:dep23}
|
---|
598 | \end{figure}
|
---|
599 |
|
---|
600 |
|
---|
601 | \begin{figure}[htp]
|
---|
602 | \centering
|
---|
603 | \includegraphics[width=0.95\linewidth]{TimeResVsCharge-Area-30.eps}
|
---|
604 | \caption{Reconstructed mean arrival time resolutions as a function of the extracted signal
|
---|
605 | for the digital filter with inverted high- and low-gain weights (extractor \#30). Error bars denote the
|
---|
606 | spread (RMS) of the time resolutions over the investigated channels.
|
---|
607 | The marker colours show the applied
|
---|
608 | pulser colour, except for the last (green) point where all three colours were used.}
|
---|
609 | \label{fig:time:dep30}
|
---|
610 | \end{figure}
|
---|
611 |
|
---|
612 |
|
---|
613 | \begin{figure}[htp]
|
---|
614 | \centering
|
---|
615 | \includegraphics[width=0.95\linewidth]{TimeResVsCharge-Area-32.eps}
|
---|
616 | \caption{Reconstructed mean arrival time resolutions as a function of the extracted signal
|
---|
617 | for the digital filter (extractor \#32). Error bars denote the
|
---|
618 | spread (RMS) of the time resolutions over the investigated channels.
|
---|
619 | The marker colours show the applied
|
---|
620 | pulser colour, except for the last (green) point where all three colours were used.}
|
---|
621 | \label{fig:time:dep32}
|
---|
622 | \end{figure}
|
---|
623 |
|
---|
624 |
|
---|
625 |
|
---|
626 |
|
---|
627 |
|
---|
628 | \clearpage
|
---|
629 |
|
---|
630 | \subsection{Pulpo Pulses}
|
---|
631 | \subsection{MC Data}
|
---|
632 | \subsection{Cosmics Data?}
|
---|
633 | The results of this subsection are based on the following runs taken
|
---|
634 | on the 21st of September 2004.
|
---|
635 | \begin{itemize}
|
---|
636 | \item{Run 39000}: OffCrab11 at 19.1 degrees zenith angle and 106.2
|
---|
637 | azimuth.
|
---|
638 | \item{Run 39182}: CrabNebula at 19.0 degrees zenith angle and 106.0 azimuth.
|
---|
639 | \end{itemize}
|
---|
640 |
|
---|
641 | \subsection{Pedestals}
|
---|
642 |
|
---|
643 |
|
---|
644 | %%% Local Variables:
|
---|
645 | %%% mode: latex
|
---|
646 | %%% TeX-master: "MAGIC_signal_reco"
|
---|
647 | %%% TeX-master: "MAGIC_signal_reco."
|
---|
648 | %%% TeX-master: "MAGIC_signal_reco"
|
---|
649 | %%% TeX-master: "MAGIC_signal_reco"
|
---|
650 | %%% TeX-master: "MAGIC_signal_reco"
|
---|
651 | %%% TeX-master: "MAGIC_signal_reco"
|
---|
652 | %%% End:
|
---|