Programa de procesamiento

El procesamiento de los datos contenidos en los archivos de salida es realizado por un segundo programa, escrito en C que realiza las siguientes tareas:

  • Captura desde el archivo de datos especificado por linea de comandos (opcion -f) las trazas del osciloscopio en una matriz de 4x500xn (n=número de eventos registrados) y los timestamps en un vector de n elementos.
  • Calcula para cada evento el tiempo desde el evento anterior.
  • Calcula para cada traza la linea de base (promedio de las 50 primeras muestras) y se la resta a cada muestra de la traza.
  • Encuentra el máximo de cada traza.
  • Calcula el área bajo cada traza (ver normalización).
  • Realiza una salida en múltiples columnas a la standard output.
El programa permite definir un threshold de nivel (opcion -t en la linea de comando). Solo envía a la salida los eventos que en algunos de sus canales superaron dicho nivel. La diferencia de tiempo se calcula desde el evento anterior que superó el threshold. El área se calcula como la suma de todas las muestras que superaron el threshold.

Usage:

$analisis [-f FILE] [-t THRESHOLD] [-p PERIOD]
-f FILE                 set intput file
-t THRESHOLD    set threshold level (0 to 128)
-p PERIOD           set stats period in seconds

Example:

veiga@beba:~/SCOPE$ ./analisis -f 2008-11-07-10001.dat -t 20 -p 10000
# 1226056272
# File starts at Fri Nov 7 08:11:12 2008
# It contains 216 events
# Baseline compensation is performed by trace (first: -2 -2 -1 0)
# The following events present at least one strip above a threshold of 20:
#
# |---------------|-------|-------|-------|-------|-------|-------|-------|-------|-------|
# |Timestamp |Diff |Peak0 |Peak1 |Peak2 |Peak3 |Area0 |Area1 |Area2 |Area3
# |---------------|-------|-------|-------|-------|-------|-------|-------|-------|-------|
4 08:23:14 722 22 8 8 12 22 0 0 0
11 08:41:29 1095 34 22 20 25 256 44 0 72
13 08:46:26 297 22 16 15 41 44 0 0 293
14 08:56:40 614 40 23 24 43 532 178 177 512
41 09:53:49 3429 32 15 12 19 257 0 0 0
43 10:05:40 711 35 17 20 31 275 0 0 206
54 10:21:56 976 9 6 10 25 0 0 0 184
55 10:22:45 49 22 18 10 9 65 0 0 0
58 10:27:08 263 42 19 14 19 441 0 0 0
64 10:36:58 590 24 11 8 9 70 0 0 0
67 10:44:42 464 93 62 60 64 1279 781 696 943
# |---------------|-------|-------|-------|-------|-------|-------|-------|-------|-------| 11 events, mean diff 837 s
75 11:06:04 1282 128 128 127 127 7916 6932 8690 7086
85 12:20:58 4494 25 13 11 14 94 0 0 0
98 12:45:38 1480 12 9 10 25 0 0 0 91
112 13:15:10 1772 23 9 9 12 44 0 0 0
117 13:26:53 703 11 8 9 20 0 0 0 0
127 13:40:54 841 20 12 14 27 0 0 0 163
128 13:43:43 169 24 19 16 14 68 0 0 0
# |---------------|-------|-------|-------|-------|-------|-------|-------|-------|-------| 7 events, mean diff 1534 s
129 13:45:36 113 24 9 8 11 67 0 0 0
131 13:49:39 243 20 13 12 27 0 0 0 180
147 14:21:24 1905 12 9 10 20 0 0 0 0
151 14:27:33 369 20 11 16 32 0 0 0 308
152 14:27:57 24 21 8 6 9 21 0 0 0
154 14:28:35 38 26 14 15 21 141 0 0 21
157 14:32:41 246 26 14 12 18 163 0 0 0
166 14:41:58 557 11 8 10 25 0 0 0 141
177 14:56:58 900 18 11 14 27 0 0 0 96
179 14:58:43 105 51 52 83 122 444 737 1227 1981
182 15:00:08 85 25 16 19 26 181 0 0 274
184 15:06:14 366 28 10 11 14 123 0 0 0
197 15:22:46 992 21 14 14 35 21 0 0 279
200 15:29:56 430 23 17 19 27 23 0 0 75
206 15:36:24 388 25 16 22 67 48 0 65 667
208 15:40:34 250 128 128 127 127 9074 8845 6856 6998
# |---------------|-------|-------|-------|-------|-------|-------|-------|-------|-------| (incomp) 16 events ->[8], mean diff 438 s
# |Timestamp |Diff |Peak0 |Peak1 |Peak2 |Peak3 |Area0 |Area1 |Area2 |Area3
# |---------------|-------|-------|-------|-------|-------|-------|-------|-------|-------| 34 events, mean diff 793 s
# 34 events matched the selection criteria
# Mean arrival time is 793 seconds



Sobre esta entrada

Esta página contiene una sola entrada realizada por Alejandro Veiga y publicada el 4 de Noviembre 2008 9:10 PM.

Sexta sesión de medición de muones: ganancia es la entrada anterior en este blog.

Apertura del detector es la entrada siguiente en este blog.

Encontrará los contenidos recientes en la página principal. Consulte los archivos para ver todos los contenidos.