Changeset 6559


Ignore:
Timestamp:
02/16/05 20:38:27 (20 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/TDAS-Extractor
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/TDAS-Extractor/Criteria.tex

    r6512 r6559  
    1 \section{Criteria for the Optimal Signal Extraction (STILL TO DO!!)}
     1\section{Criteria for the Optimal Signal Extraction}
    22
    3 The goal for the optimal signal reconstruction algorithm is to compute an unbiased estimate of the strenght and arrival time of the Cherenkov signal with the highest possible resolution for all signal intensities. The algorithm shall be stable with respect to changes in observation conditions and background levels. Also the needed computing time is of concern.
    4 
     3The goal for the optimal signal reconstruction algorithm is to compute an unbiased estimate of the strength and arrival time of the
     4Cherenkov signal with the highest possible resolution for all signal intensities. The MAGIC telescope has been optimized to
     5lower the energy treshold of observation in any respect. Particularly the choice for an FADC system has been made with an eye on the
     6possibility to extract the smallest possible signals from air showers. It would be inconsequent not to continue the optimization procedure
     7in the signal extraction algorithms and the subsequent image cleaning.
     8\par
     9In the image analysis, one hake the decision whether the extracted signal of a certain pixel is considered as signal or background.
     10Those considered as signal are further used to compute the image parameters while the background ones are simply rejected. The calculation
     11of the second moments of the image ``ellipse'' usually fails when applied to un-cleaned images, therefore the decision is yes or no. Moreover,
     12already low contributions of mis-estimated background can degrade the resolution of the image parameters considerably. If one wants to
     13lower the threshold for signal recognition, it is therefore mandatory to increase the efficiency with which the background is recognized as
     14such. If the background resolution is bad, the signal threshold goes up and vice versa.
     15\par
     16The algorithm must be stable with respect to changes
     17in observation conditions and background levels and between signals induced from gamma or hadronic showers or from muons.
    518The reconstructed signal shall be proportional to the total integrated charge in the FADCs due to the PMT pulse from the Cherenkov signal.
    619
    7 Discussion about the signal to noise for the image cleaning a la Maxim.
     20Also the needed computing time is of concern.
    821
     22\subsection{Bias and Mean-squared Error}
    923
     24Consider a large number of same signals $S$. By applying a signal extractor
     25we obtain a distribution of estimated signals $\widehat{S}$ (for fixed $S$ and
     26fixed background fluctuations $BG$). The distribution of the quantity
    1027
    11 \subsection{Resolution and Bias}
    12 \ldots {\textit{The jitter to identical input pulses is measured, for times, amplitudes,
    13 high-gain and low-gain pulses and different signal and background levels }}
     28\begin{equation}
     29X = \widehat{S}-S
     30\end{equation}
     31
     32has the mean $B$ and the Variance $MSE$ defined as:
     33
     34\begin{eqnarray}
     35   B   \ \ \ \  = \ \ \ \ \ \ <X> \ \ \ \ \  &=& \ \ <\widehat{S}> \ -\ S\\
     36   R   \ \ \ \  = \ <(X-B)^2> &=& \ Var[\widehat{S}]\\
     37   MSE \      = \ \ \ \ \ <X^2> \ \ \ \  &=& \ Var[\widehat{S}] +\ B^2
     38\end{eqnarray}
     39
     40The parameter $B$ is also called the {\textit{\bf BIAS}} of the estimator and $MSE$
     41the {\textit{\bf MEAN-SQUARED ERROR}} which combines the variance of $\widehat{S}$ and
     42the bias. Both depend generally on the size of $S$ and the background fluctuations $BG$,
     43thus: $B = B(S,BG)$ and $MSE = MSE(S,BG)$.
     44
     45\par
     46Usually, one measures easily the parameter $R$, but needs the $MSE$ for statistical analysis (e.g.
     47in the image cleaning).
     48However, only in case of a vanishing bias $B$, the two numbers are equal. Otherwise,
     49the bias $B$ has to be known beforehand. Note that every sliding window extractor has a
     50bias, especially at low or vanishing signals $S$.
    1451
    1552\subsection{Linearity}
  • trunk/MagicSoft/TDAS-Extractor/Pedestal.tex

    r6498 r6559  
    1111By definition, the $\boldsymbol{B}$ and thus the ``pedestal RMS''
    1212is independent from the signal extractor.
    13 
    14 \subsection{Bias and Mean-squared Error}
    15 
    16 Consider a large number of same signals $S$. By applying a signal extractor
    17 we obtain a distribution of estimated signals $\widehat{S}$ (for fixed $S$ and
    18 fixed background fluctuations $BG$). The distribution of the quantity
    19 
    20 \begin{equation}
    21 X = \widehat{S}-S
    22 \end{equation}
    23 
    24 has the mean $B$ and the Variance $MSE$ defined as:
    25 
    26 \begin{eqnarray}
    27    B   \ \ \ \  = \ \ \ \ \ \ <X> \ \ \ \ \  &=& \ \ <\widehat{S}> \ -\ S\\
    28    R   \ \ \ \  = \ <(X-B)^2> &=& \ Var[\widehat{S}]\\
    29    MSE \      = \ \ \ \ \ <X^2> \ \ \ \  &=& \ Var[\widehat{S}] +\ B^2
    30 \end{eqnarray}
    31 
    32 The parameter $B$ is also called the {\textit{\bf BIAS}} of the estimator and $MSE$
    33 the {\textit{\bf MEAN-SQUARED ERROR}} which combines the variance of $\widehat{S}$ and
    34 the bias. Both depend generally on the size of $S$ and the background fluctuations $BG$,
    35 thus: $B = B(S,BG)$ and $MSE = MSE(S,BG)$.
    36 
    37 \par
    38 Usually, one measures easily the parameter $R$, but needs the $MSE$ for statistical analysis (e.g.
    39 in the image cleaning).
    40 However, only in case of a vanishing bias $B$, the two numbers are equal. Otherwise,
    41 the bias $B$ has to be known beforehand. Note that every sliding window extractor has a
    42 bias, especially at low or vanishing signals $S$.
    4313
    4414\subsection{Pedestal Fluctuations as Contribution to the Signal Fluctuations \label{sec:ffactor}}
Note: See TracChangeset for help on using the changeset viewer.