Bachelorarbeit - Desy

26
8 break ;
9 }
10 }
11 f o r ( i = alpha ; i < alpha +100; i ++){
12 int d i f f 3 = buf [ i +1] − buf [ i ] ;
13 // Wenn d i e D i f f e r e n z z w i s c h e n einem Punkt und dem Nächsten k l e i n e r a l s 100 i s t ,
14 // s e t z t e d i e s e n Punkt a l s beta und beende d i e S c h l e i f e .
15 // Dadurch wird der Endpunkt der A n s t i e g s f l a n k e e r k a n n t .
16 i f ( d i f f 3 < 100){
17 beta = i ;
18 break ;
19 }
20 }
Das Verwerfen von Signalen mit Pile-up kann nun durch eine einfache Methode rea-
lisiert werden. Das Eintreffen eines plötzlichen Anstiegs der Datenpunkte in einer Trace
wird durch die Punkte alpha und beta gekennzeichnet. Durch die Forderung, dass nach
Zeitpunkt beta die Differenz zweier Punkte einen vom Benutzer festgelegten Schwel-
lenwert nicht überschreitet, kann Pile-up erkannt werden. Um Traces, in denen Pile-up
auftritt, zu verwerfen, kann von einer Booleschen Variablen 3 Gebrauch gemacht werden.
Diese wird pileup genannt und kann die Werte true und false annehmen. Zu Beginn
des Programms wird sie auf false gesetzt und die Auswertung der Daten einer Trace nur
fortgesetzt, wenn pileup==false erfüllt ist. Wird durch die beschriebene Routine Pile-up
erkannt, wird pileup=true gesetzt, diese Trace also nicht weiter bearbeitet:
1 bool p i l e u p = f a l s e ;
2 int d i f f , int n=4000; // n = Länge der Trace
3 double beta ; // Endpunkt der A n s t i e g s f l a n k e
4 f o r ( i = beta ; i < n ; i ++) {
5 d i f f = ( buf [ i +1] − buf [ i ] ) ;
6 // Wenn h i n t e r der e r s t e n A n s t i e g s f l a n k e e i n e w e i t e r e e i n t r i t t , s e t z e p i l e u p=t r u e
7 i f ( d i f f > 100){
8 p i l e u p = true ;
9 break ;
10 }
11 }
12 // Nur wenn p i l e u p==f a l s e , f a h r e mit der V e r a r b e i t u n g f o r t
13 i f ( p i l e u p == f a l s e ){
14 // H i e r f o l g t d i e V e r a r b e i t u n g der Traces
15 }
Neben dem Auftreten von zwei Signalen in einer Trace kommt es vor, dass in einer
Trace ein Signal beginnt, obwohl das vorige Signal noch nicht auf die Baseline abgefallen
ist. Dies macht sich über eine nicht-konstante Baseline bemerkbar. Um dies zu verhindern,
kann gefordert werden, dass die Baseline, welche durch Mittelwertbildung bestimmt wur-
de, eine gewisse Standardabweichung des Mittelwerts nicht überschreiten darf. Des Wei-
teren kommt es gelegentlich vor, dass ein Signal zunächst abfällt und anschließend wieder
leicht ansteigt. Da solche Signale keinen exponentiellen Abfall über die gesamte Abstiegs-
flanke zeigen, werden solche Traces ebenfalls verworfen. Die letzte Art unerwünschter
Signale sind solche, die unter die Baseline abfallen. Diese Verwerfungsmethoden werden
folgendermaßen durchgeführt:
1 double sd = 0 ; // s t a n d a r d d e v i a t i o n
2 f o r ( i = 0 ; i