Changeset 5925 for trunk/MagicSoft/TDAS-Extractor
- Timestamp:
- 01/21/05 11:12:16 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/TDAS-Extractor/Algorithms.tex
r5919 r5925 325 325 the computed spline value at the arrival time position. Can be adjusted for the half-maximum time extraction 326 326 type. 327 \item[Lo GainStretch:\xspace] Can be adjusted to account for the bigger rise and fall timein the327 \item[Low Gain Stretch:\xspace] Can be adjusted to account for the larger rise and fall times in the 328 328 low-gain as compared to the high gain pulses (default: 1.5) 329 329 \end{description} … … 333 333 \includegraphics[width=0.49\linewidth]{MExtractTimeAndChargeSpline_23Led_Blue.eps} 334 334 \caption[Sketch calculated arrival times MExtractTimeAndChargeSpline]{% 335 Sketch of the calculated arrival times for the extractor {\textit{ MExtractTimeAndChargeSpline}}335 Sketch of the calculated arrival times for the extractor {\textit{\bf MExtractTimeAndChargeSpline}} 336 336 for two typical calibration pulses (pedestals have been subtracted) and a typical inner pixel. 337 337 The extraction window sizes modify the position of the (amplitude-weighted) mean FADC-slices slightly. 338 The pulse would be shifted half a slice to the right for an outer pixel s. }338 The pulse would be shifted half a slice to the right for an outer pixel. } 339 339 \label{fig:splinesketch} 340 340 \end{figure} … … 342 342 \subsubsection{Digital Filter} 343 343 344 This extractor is implemented in the MARS-class {\textit{MExtractTimeAndChargeDigitalFilter}}. 345 346 347 The goal of the digital filtering method \cite{OF94,OF77} is to optimally reconstruct the amplitude and time origin of a signal with a known signal shape from discrete measurements of the signal. Thereby the noise contribution to the amplitude reconstruction is minimized. 348 349 For the digital filtering method two assumptions have to be made: 344 This extractor is implemented in the MARS-class {\textit{\bf MExtractTimeAndChargeDigitalFilter}}. 345 346 347 The goal of the digital filtering method \cite{OF94,OF77} is to optimally reconstruct the amplitude and time origin of a signal with a known signal shape 348 from discrete measurements of the signal. Thereby, the noise contribution to the amplitude reconstruction is minimized. 349 350 For the digital filtering method, three assumptions have to be made: 350 351 351 352 \begin{itemize} … … 355 356 \end{itemize} 356 357 357 Let $g(t)$ be the normalized signal shape, $E$ the signal amplitude and $\tau$ the time shift of the physical signal from the predicted signal shape. Then the time dependence of the signal, $y(t)$, is given by: 358 \par 359 \ldots {\textit{\bf IS THIS TRUE FOR MAGIC???? }} \ldots 360 \par 361 362 Let $g(t)$ be the normalized signal shape, $E$ the signal amplitude and $\tau$ the time shift 363 of the physical signal from the predicted signal shape. Then the time dependence of the signal, $y(t)$, is given by: 358 364 359 365 \begin{equation} … … 361 367 \end{equation} 362 368 363 where $b(t)$ is the time dependent noise distribution. For small time shifts $\tau$ (usually smaller than369 where $b(t)$ is the time-dependent noise contribution. For small time shifts $\tau$ (usually smaller than 364 370 one FADC slice width), 365 371 the time dependence can be linearized by the use of a Taylor expansion: … … 376 382 \end{equation} 377 383 378 The correlation of the noise contributions at times $t_i$ and $t_j$ is expressed in the noise autocorrelation matrix $\boldsymbol{B}$: 384 The correlation of the noise contributions at times $t_i$ and $t_j$ can be expressed in the 385 noise autocorrelation matrix $\boldsymbol{B}$: 379 386 380 387 \begin{equation} 381 \boldsymbol{B_{ij}} = \langle b_i b_j \rangle - \langle b_i \rangle \langle b_j388 B_{ij} = \langle b_i b_j \rangle - \langle b_i \rangle \langle b_j 382 389 \rangle \ . 383 390 \label{eq:autocorr} … … 385 392 %\equiv \langle b_i b_j \rangle with $\langle b_i \rangle = 0$. 386 393 387 The signal amplitude, $E$, and the product of amplitude and time shift, $E \tau$, can be estimated from the given set of 388 measurements $\boldsymbol{y} = (y_1, ... ,y_n)$ by minimizing: 394 The signal amplitude $E$, and the product of amplitude and time shift $E \tau$, can be estimated from the given set of 395 measurements $\boldsymbol{y} = (y_1, ... ,y_n)$ by minimizing the excess noise contribution with respect to the known noise 396 auto-correlation: 389 397 390 398 \begin{eqnarray} 391 \chi^2(E, E\tau) &=& \sum_{i,j}(y_i-E g_i-E\tau \dot{g}_i) \boldsymbol{B}^{-1}_{ij} (y_j - E g_j-E\tau \dot{g}_j) \\399 \chi^2(E, E\tau) &=& \sum_{i,j}(y_i-E g_i-E\tau \dot{g}_i) B^{-1}_{ij} (y_j - E g_j-E\tau \dot{g}_j) \\ 392 400 &=& (\boldsymbol{y} - E 393 401 \boldsymbol{g} - E\tau \dot{\boldsymbol{g}})^T \boldsymbol{B}^{-1} (\boldsymbol{y} - E \boldsymbol{g}- E\tau \dot{\boldsymbol{g}}) \ , 394 402 \end{eqnarray} 395 403 396 where the last expression is matricial. The minimum is obtained for: 404 where the last expression is matricial. $\chi^2$ is a continuous function of $\tau$ and will have to be discretized itself for a 405 desired resolution. 406 The minimum of $\chi^2$ is obtained for: 397 407 398 408 \begin{equation} … … 410 420 \end{eqnarray} 411 421 412 Solving these equations one gets the solutions:422 Solving these equations one gets the following solutions: 413 423 414 424 \begin{equation} … … 422 432 423 433 Thus $\overline{E}$ and $\overline{E\tau}$ are given by a weighted sum of the discrete measurements $y_i$ 424 with the digital filtering weights for the amplitude, $w_{\text{amp},i}$, and time shift, $w_{\text{time},i}$. 425 434 with the digital filtering weights for the amplitude, $w_{\text{amp}}(t)$, and time shift, $w_{\text{time}}(t)$. 435 Note the remaining time dependence of the two weights which follow from the dependency of $\boldsymbol{g}$ and 436 $\dot{\boldsymbol{g}}$ on the position of the pulse with respect to the FADC bin positions. 437 \par 426 438 Because of the truncation of the Taylor series in equation (\ref{shape_taylor_approx}) the above results are 427 439 only valid for vanishing time offsets $\tau$. For non-zero time offsets one has to iterate the problem using … … 449 461 450 462 \begin{equation}\label{of_noise} 451 \sigma_E \approx 8.3 \ \mathrm{FADC\ counts} \qquad \sigma_{\tau} \approx \frac{6.5\ \Delta T_{\mathrm{FADC}}}{ E\ /\ \mathrm{FADC\ counts}} \ ,463 \sigma_E \approx 8.3 \ \mathrm{FADC\ counts} \qquad \sigma_{\tau} \approx \frac{6.5\ \Delta T_{\mathrm{FADC}}}{(E\ /\ \mathrm{FADC\ counts})} \ , 452 464 \end{equation} 453 465 … … 717 729 %%% TeX-master: "MAGIC_signal_reco" 718 730 %%% TeX-master: "MAGIC_signal_reco" 731 %%% TeX-master: "MAGIC_signal_reco" 719 732 %%% End:
Note:
See TracChangeset
for help on using the changeset viewer.