Koincydencje, veto i bramkowanie: logika czasowa w pomiarach jądrowych

Streszczenie

Nie każde zdarzenie zarejestrowane przez detektor jest zdarzeniem, którego szukamy. Czasem interesuje nas tylko sytuacja, w której dwa detektory zareagują prawie jednocześnie. Czasem przeciwnie: chcemy zaakceptować impuls tylko wtedy, gdy drugi detektor milczy.

Koincydencja, antykoincydencja, veto i bramkowanie są językiem takich decyzji czasowych. W klasycznej aparaturze jądrowej realizowały je moduły NIM i CAMAC; dziś często robi to FPGA albo cyfrowy system akwizycji. Fizyka decyzji pozostaje jednak ta sama: trzeba określić okno czasowe, opóźnienia, progi i konsekwencje fałszywych zbieżności.¹

Rozszerzenie tematu

Koincydencja odpowiada na pytanie: czy dwa sygnały pojawiły się w zadanym oknie czasu? Jeżeli dwa detektory rejestrują fotony anihilacyjne w PET, koincydencja pomaga odróżnić pary pochodzące z jednego zdarzenia od tła. Jeżeli dwa detektory gamma patrzą na kaskadę jądrową, koincydencja gamma-gamma pozwala badać relacje między poziomami energetycznymi.

Antykoincydencja działa odwrotnie. Zdarzenie w głównym detektorze jest odrzucane, jeśli detektor osłonowy także zarejestrował impuls. Tak działa np. supresja Comptona albo osłony antykosmiczne. Veto nie mówi "to jest dobry sygnał", lecz "ten sygnał jest podejrzany, bo towarzyszy mu coś, czego nie chcemy".

Bramkowanie jest ogólniejsze. Gate może otworzyć analizator tylko w określonym czasie, zsynchronizować pomiar z impulsem akceleratora albo dopuścić do zliczania tylko zdarzenia spełniające warunek logiczny. W NIM standardowe piny reset/gate istniały dla wolniejszych sygnałów, natomiast szybkie sygnały prowadzono kablami front-panel, żeby uniknąć przesłuchów na magistrali.²

Najtrudniejsze jest okno czasowe. Jeśli jest zbyt wąskie, tracimy prawdziwe zdarzenia przez jitter, różne czasy propagacji i rozrzut odpowiedzi detektorów. Jeśli jest zbyt szerokie, rośnie liczba koincydencji przypadkowych. Dlatego ćwiczenia z koincydencji są dobrym wprowadzeniem do pojęć: rozdzielczość czasowa, opóźnienie, częstość tła i prawdopodobieństwo przypadkowego zbiegu.

W aparaturze NIM logika czasowa była widoczna fizycznie: kabel opóźniający, moduł coincidence, moduł gate, licznik. W cyfrowym systemie DPP te same decyzje mogą być w konfiguracji programu. To wygodniejsze, ale mniej przejrzyste. Student powinien rozumieć, że kliknięcie "coincidence window = 100 ns" jest odpowiednikiem konkretnej logiki impulsowej.

W kontekście broni jądrowej i testów jądrowych bezpieczne zastosowanie tego tematu to diagnostyka i monitoring: rejestracja sygnałów z wielu czujników, odrzucanie tła, synchronizacja kanałów pomiarowych, CTBT, bhangmetry i pomiary laboratoryjne. Nie jest potrzebne opisywanie sekwencji uzbrajania ani inicjacji, aby nauczyć studentów sensu logiki czasowej.

Fizyczne podstawy techniki koincydencji

Technika koincydencji opiera się na probabilistycznym argumencie: jeśli dwa niezależne detektory mają szanse rejestracji $\varepsilon_1$ i $\varepsilon_2$ z częstościami $n_1$ i $n_2$, to prawdziwe koincydencje (pochodzące z tego samego zdarzenia fizycznego) pojawiają się z częstością:

$$n_{true} = n_{source} \cdot \varepsilon_1 \cdot \varepsilon_2$$

podczas gdy przypadkowe koincydencje (dwa niezależne zdarzenia mieszczące się w oknie $\Delta t$):

$$n_{acc} = 2 \cdot \Delta t \cdot n_1 \cdot n_2$$

Stosunek sygnału do tła:

$$\frac{n_{true}}{n_{acc}} = \frac{n_{source} \cdot \varepsilon_1 \cdot \varepsilon_2}{2 \cdot \Delta t \cdot n_1 \cdot n_2}$$

Jeśli $n_1 \approx n_{source} \cdot \varepsilon_1$ i $n_2 \approx n_{source} \cdot \varepsilon_2$ (zakładając, że detektor mierzy głównie interesujący sygnał), to stosunek sygnał/tło rośnie gdy $\Delta t$ maleje. Stąd motywacja do minimalizacji okna koincydencji — przy zachowaniu wystarczającego marginesu dla jitter detektorów.

Przykład liczbowy: Dwa detektory beta, każdy o $n_1 = n_2 = 1000$ cps, badające źródło o aktywności $A_{true} = 100$ cps. Prawdziwa koincydencja: $n_{true} = 100 \times 0.9 \times 0.9 = 81$ cps. Przypadkowe przy $\Delta t = 100$ ns: $n_{acc} = 2 \times 10^{-7} \times 1000 \times 1000 = 0.2$ cps. SNR = 81/0.2 = 400. Jeśli okno rozszerzymy do $\Delta t = 10$ μs: $n_{acc} = 20$ cps, SNR = 4. Widać jak okno koincydencji jest krytyczne.

Rozdzielczość czasowa i FWHM okna koincydencji

Rozdzielczość czasowa układu koincydencji (CRT — Coincidence Resolving Time) jest parametrem opisującym jak dokładnie system może określić jednoczesność dwóch zdarzeń. Mierzy się ją jako FWHM rozkładu coincidence time spectrum.

Czynniki wpływające na CRT:

Jitter dyskryminatora: Zmienność czasu wyzwolenia dyskryminatora (SLD lub CFD) spowodowana szumem elektronicznym. Dla sygnałów o amplitudzie $V$ i czasie narastania $\tau_r$: $\sigma_{disc} \approx V_{noise} / (dV/dt)|_{threshold} = V_{noise} \cdot \tau_r / V_{signal}$. CFD redukuje jitter walk (zależność od amplitudy), ale nie usuwa jitter szumowy.
Statystyczny jitter scyntylatora: Dla scyntylatora o czasie zaniku $\tau_{dec}$ i liczbie fotoelektronów N_pe, jitter: $\sigma_{scint} \approx \tau_{dec} / \sqrt{N_{pe}}$. Dla NaI(Tl): $\tau_{dec} = 250$ ns, $N_{pe} = 10^4$ przy 1 MeV → $\sigma_{scint} \approx 2.5$ ns. Dla LSO/LYSO (używany w PET): $\tau_{dec} = 40$ ns, $N_{pe} = 3000$ → $\sigma_{scint} \approx 0.7$ ns — dlatego LSO/LYSO jest preferowany w PET.
Czas tranzytu PMT i jego rozrzut (TTS — Transit Time Spread): Każdy foton wędruje różną drogą przez PMT, powodując rozrzut czasowy. TTS typowe: 200-800 ps dla klasycznych PMT bialkali; 50-100 ps dla MCP-PMT (Microchannel Plate PMT).
Opóźnienia kablowe i ich różnice: Nawet krótkie różnice długości kabli (~1 ns/20 cm) powodują przesunięcie centrum okna koincydencji. Kompensacja: kable wyrównujące długości.

Łączna CRT dla dwóch detektorów:

$$CRT_{FWHM} = 2.35 \sqrt{\sigma_{det1}^2 + \sigma_{det2}^2}$$

Współczesne systemy PET z kryształami LSO/LYSO i SiPM osiągają CRT < 200 ps FWHM. Klasyczne NaI(Tl) + PMT: CRT ~ 5-10 ns.

Opóźnienia kablem i elektryczne: kalibracja toru koincydencji

Kalibracja toru koincydencji polega na wyznaczeniu opóźnień każdego kanału i ustawieniu ich tak, by prawdziwe koincydencje (zdarzenia jednoczesne) trafiały w centrum okna koincydencji:

Pomiar opóźnienia kablowego: Prędkość propagacji sygnału w kablu RG-58: $v_p = 0.66c = 2 \times 10^8$ m/s. Opóźnienie: $t_{delay} = l/v_p = l / (0.20 \text{ m/ns})$. Dla kabla 50 cm: $t_{delay} = 2.5$ ns. Dla kabla 10 m: $t_{delay} = 50$ ns.

Pomiar całkowitego opóźnienia kanału: Używając generatora impulsów jako wspólnego źródła dla obu kanałów, zmierz opóźnienie od wejścia do wyjścia CFD → wyjście koincydencji. Różnica opóźnień między kanałami = korekcja kablowa do zastosowania.

Kalibracja CRT: Umieść źródło ⁶⁰Co (emituje dwa fotony gamma 1173 i 1332 keV — ale te nie są koincydentne; lepiej użyć ²²Na emitującego 511+511 keV w koincydencji lub ¹⁷⁴Lu kaskadowe gamma) między dwoma detektorami. Zarejestruj spectrum różnic czasu (time difference spectrum). Gaussian peaked at zero = właściwe opóźnienia. Shift od zera = dodaj kabel wyrównujący.

Realizacja hardwarowa koincydencji NIM: szczegóły implementacyjne

W klasycznym systemie NIM koincydencja realizowana jest przez moduł AND gate (bramka koincydencyjna). Oto szczegóły implementacyjne:

Moduł coincidence NIM: Typowy moduł (np. CAEN N89, LeCroy 622) ma 4 wejścia NIM fast logic i 1 wyjście. Wyjście = "1" gdy $n$ z 4 wejść jest aktywne (konfigurowalne, np. 2-fold, 3-fold, 4-fold). Czas odpowiedzi: 5-20 ns od zmiany wejść do wyjścia. Wewnętrzna realizacja: bramki ECL szybkie.

Opóźniacze kablowe (delay cables): Standardowy sposób wyrównywania opóźnień: długość kabla dobierana by wyrównać różnice opóźnień między kanałami. Kable dedykowane: 1 ns, 2 ns, 5 ns, 10 ns, 20 ns (często w zestawach kalibracyjnych). Alternatywa: Delay Line Unit NIM (programowalny delay, np. CAEN N108A: 0-128 ns krokiem 1 ns).

Look-At-Me timing w systemach z CAMAC: W CAMAC, moduł ADC zgłasza LAM (Look-At-Me) gdy konwersja jest gotowa. Kontroler musi odczytać dane przed utratą (następnym zdarzeniem). Okno czasu między LAM a odczytem jest kluczowe dla systemów o wysokich rate. Przy koincydencji z CAMAC: gate do modułu ADC CAMAC musi przychodzić z opóźnieniem odpowiednim do czasu konwersji ADC (typowo 2-10 μs dla Wilkinson ADC).

Układy veto: realizacja i zastosowania

Veto to sygnał NIM logic "odrzuć bieżące zdarzenie". Realizacja jako sygnał antykoincydencji: zdarzenie głównego detektora akceptowane tylko gdy veto = 0:

$$Wyjście = A \cdot \overline{B}$$

gdzie A = sygnał głównego detektora, B = sygnał detektora veto.

Realizacja NIM: moduł ANTICOINCIDENCE (AND-NOT gate). Wejście A podawane przez Delay Unit (by pewnie "zakryć" czas trwania sygnału veto B). Typowy schemat:

Detektor A (main) → CFD → Delay (10-50 ns) → wejście A modułu veto
Detektor B (veto) → CFD → Pulse Width Stretcher (rozszerzenie do 50-200 ns) → wejście B modułu veto
Wyjście = A · B̄ → Gate do ADC lub MCA

Pulse Width Stretcher (rozszerzacz impulsów): moduł NIM który odbiera krótki impuls logiczny i generuje dłuższy impuls wyjściowy (np. 20 ns wejście → 100 ns wyjście). Ważne: veto musi trwać wystarczająco długo, by pokryć całe okno potencjalnego zdarzenia w detektorze A.

Zastosowania veto:

Supresja Comptona: Spektrometr antykoincydencyjny Comptona (Compton Suppression Spectrometer): HPGe w centrum, scyntylatory NaI(Tl) otaczające jako "guard ring". Zdarzenia w HPGe rejestrowane tylko gdy żaden NaI nie dał sygnału. Efekt: eliminacja zdarzeń, gdzie foton oddał tylko część energii przez Compton i uciekł do otoczenia. FWHM w widmie pełnej absorpcji (full energy peak) pozostaje, ale tło Comptona jest znacznie zredukowane (factor 10-100×).
Antykosmiczne veto: Eksperyment liczy rzadkie zdarzenia (np. poszukiwanie ciemnej materii). Wyładowania kosmiczne (miony) przechodzą przez detektor główny i powodują fałszywy sygnał. Otaczające plastikowe scyntylatory (muon veto) rejestrują mion i generują sygnał veto. Przy wyzwoleniu veto, detektor główny nie zapisuje zdarzenia przez ~1 μs.
Veto pulsera: Generator testowy pulsów jest wyzwalany periodycznie (np. 10 Hz). Gdy pulsator wyzwala, veto blokuje ADC — zapobiegając zarejestrowaniu sztucznego zdarzenia z pulsatora w widmie rzeczywistym. Jednocześnie scaler zlicza impulsy pulsatora dla wyznaczenia live time.

Bramkowanie (gating) i jego zastosowania

Bramkowanie (gating) polega na aktywowaniu lub dezaktywowaniu modułu pomiarowego sygnałem logicznym. Różne typy:

Gate do ADC: Analogowy gate sygnał do modułu ADC kontroluje okno integracji sygnału. Dla scyntylatora NaI(Tl): sygnał z PMT jest stosunkowo powolny (μs). ADC integruje ładunek przez cały czas świecenia scyntylatora. Gate musi być: szeroki (pokrywać cały impuls), zsynchronizowany (zaczynać się przed impulsem), zakończony przed następnym zdarzeniem.

Trigger z akceleratorem: Eksperymenty przy akceleratorach mają pulsy wiązki o ustalonej częstotliwości (np. CERN SPS: 40 MHz pakiety protonów, LHC: 40 MHz). Trigger zewnętrzny zsynchronizowany z pakietem wiązki otwiera gate detektorów tylko podczas pakietu → eliminacja zdarzeń tła między pakietami.

Bramkowanie MCA: Analogowy MCA (Multichannel Analyzer) ma wejście GATE: MCA akumuluje widmo tylko gdy GATE = 1. Zastosowania: zbieranie widma tylko podczas naświetlania aktywacyjnego (gate = czas naświetlania), lub wyłącznie w wybranych interwałach czasu.

Digital gate (soft gate) w DPP: W cyfrowym DPP, "gate" jest warunkiem programowym: zdarzenie jest zapisywane tylko gdy spełnione są warunki (timestamp w oknie koincydencji, amplituda w zakresie, brak veto). Jest to pełnowartościowy odpowiednik hardware gate, ale realizowany w FPGA lub software.

Techniki redukcji tła przez koincydencje: przykłady eksperymentalne

Koincydencje gamma-gamma: Badanie kaskad jądrowych. Dwa detektory gamma mierzą fotony z kaskady wzbudzonego jądra (np. ²⁴⁴Am: α→²⁴⁰Np* → ²⁴⁰Np + γ kaskada). Macierz koincydencji E₁×E₂ pokazuje linie korelacji dla par gamma z tej samej kaskady. Tło (nieskorelowane) tworzy płaski kontynuum. Gated spectrum: przekrój pionowy lub poziomy macierzy = widmo gamma w koincydencji z wybraną linią.

β-γ koincydencja: Detektor beta (licznik GM, scyntylator plastikowy) + detektor gamma (NaI, HPGe) w koincydencji. Odrzuca zdarzenia gamma niespowodowane rozpadem beta (np. tło kosmiczne). Ważne dla pomiaru aktywności izotopów beta-gamma (np. ⁶⁰Co).

Coincidence summing (efekt sumowania koincydencji): Problem, nie zaleta: jeśli dwa fotony gamma z kaskady zostaną zarejestrowane przez jeden detektor prawie jednocześnie, ADC mierzy sumę ich energii jako jedno zdarzenie. Artefakt: "suma peak" w widmie (np. dla ⁶⁰Co: 1173+1332 = 2505 keV). Efekt ten jest silniejszy przy krótszych odległościach źródło-detektor i większej efektywności detekcji. Korekcja: pomiar geometria-zależnych współczynników sumowania lub numeryczna korekcja w oprogramowaniu.

Triple coincidence: Trzy detektory w koincydencji dla jeszcze lepszego odrzucenia tła. Używane w pomiarach pół-rozpadu podwójnego beta (2νββ i 0νββ): detektor bolometryczny w koincydencji z wewnętrznym detektorem Czerenkowa i zewnętrznym muon veto. Przykład: CUORE (CUORE-0), AMoRE, NEXUS, LEGEND.

PET (Positron Emission Tomography): koincydencje w zastosowaniach medycznych

PET jest jedynym zastosowaniem koincydencji jądrowej, które niemal każdy student spotyka w kontekście medycznym. Warto omówić je szczegółowo:

Zasada PET: Pacjent otrzymuje znacznik radioaktywny emitujący pozytrony (np. ¹⁸F-FDG). Pozytron anihiluje z elektronem, emitując dwa fotony gamma 511 keV w przeciwnych kierunkach. Pierścień detektorów LOR (Lines of Response) identyfikuje parę fotonów w koincydencji (okno ~3-5 ns) jako LOR — linię prostą przechodząca przez miejsce anihilacji. Po zebraniu milionów LOR, algorytm (MLEM, FBP) rekonstruuje obraz 3D rozkładu znacznika.

Realizacja koincydencji w PET: Każdy kryształ detektora (LSO, LYSO, BGO) podłączony do PMT lub SiPM. Szybki trigger hardware dla par kryształów ustawionych naprzeciwlegle. Okno koincydencji: 3-6 ns dla LYSO+SiPM (CRT~200 ps), 8-12 ns dla BGO+PMT (CRT~5 ns).

Time-of-Flight PET (TOF-PET): Gdy CRT < 600 ps, różnica czasu dotarcia dwóch fotonów gamma 511 keV (max. 3 ns dla 45 cm pierścienia) pozwala lokalizować miejsce anihilacji wzdłuż LOR z precyzją ~60-90 mm. TOF-PET poprawia SNR obrazu przez factor $\sqrt{D/(\Delta x)}$, gdzie D = średnica FOV, Δx = localization uncertainty. Najnowsze kliniczne PET (Siemens Biograph Vision, GE Omni Legend) z LYSO+SiPM osiągają CRT~210 ps FWHM.

PET we fizyce jądrowej: Pozytronowe emittery (²²Na, ⁵⁸Co, ⁶⁸Ga) używane są jako źródła kalibracyjne dla systemów koincydencji. Podstawy rozumienia PET są bezpośrednim zastosowaniem techniki koincydencji omówionej w tym artykule.

Ćwiczenie obliczeniowe: optymalizacja okna koincydencji

Dla konkretnego systemu: dwa detektory NaI(Tl) (5"×5") z PMT, mierzące ⁶⁰Co (2 fotony gamma nieskalowalnie z kaskadą, emitowane w tym samym czasie, więc nie tworzą "prawdziwych" koincydencji, ale źródło ²²Na jest lepszym wyborem dla ćwiczeń koincydencji).

Ćwiczenie z ²²Na:

²²Na emituje pozyton → anihilacja → 2 × 511 keV (koincydentne, ~180°) + 1274 keV (w koincydencji z 511 keV przez kaskadę, ale z innym kątem)
Ustaw dwa detektory naprzeciwlegle (180°) od źródła ²²Na
Ustaw SCA na obu detektorach: okno energetyczne 450-550 keV (wybór tylko 511 keV)
Zmierz: n₁ (singles det.1), n₂ (singles det.2), n_coinc (koincydencje)
Oblicz przypadkowe: n_acc = 2·Δt·n₁·n₂
Zmieniaj Δt (przez różne kable delay lub oprogramowanie DPP)
Wyznacz optymalny Δt minimalizujący n_acc/n_coinc

Wynik: optymalne okno zależy od CRT systemu. Δt ≈ 2-3 × CRT (w ns) dla zachowania ~99% prawdziwych koincydencji przy minimalizacji przypadkowych.

Python kod do analizy:

import numpy as np

# Parametry przykładowe
n1 = 5000  # singles det. 1 [cps]
n2 = 4800  # singles det. 2 [cps]
n_coinc_true = 200  # prawdziwe koincydencje [cps] (mierzone przy małym Δt)

# Oblicz tło przypadkowe dla różnych okien
dt_values = np.logspace(-9, -6, 100)  # 1 ns do 1 μs
n_acc = 2 * dt_values * n1 * n2
ratio = n_acc / n_coinc_true

import matplotlib.pyplot as plt
plt.loglog(dt_values * 1e9, ratio)
plt.xlabel('Okno koincydencji Δt [ns]')
plt.ylabel('n_accidental / n_true')
plt.grid(True, alpha=0.3)
plt.axhline(0.01, color='r', linestyle='--', label='1% contamination')
plt.legend()
plt.title('Optymalizacja okna koincydencji')
plt.savefig('coincidence_optimization.png', dpi=150)

# Optimum: n_acc/n_true = 0.01 → dt_opt
dt_opt = 0.01 * n_coinc_true / (2 * n1 * n2)
print(f"Optymalne okno dla 1% tła: {dt_opt*1e9:.1f} ns")

Takie ćwiczenie obliczeniowe + eksperymentalne weryfikacja jest idealnym projektem laboratoryjnym łączącym teorię i praktykę koincydencji.

Koincydencje w detektorach neutronowych: problem n/γ dyskryminacji

Osobnym zagadnieniem techniki koincydencji jest rozróżnianie neutronów od fotonów gamma metodą PSD (Pulse Shape Discrimination) połączonej z koincydencją:

PSD + koincydencja: Detektor ciekłego scyntylatora (NE-213, BC-501) wykazuje różne kształty impulsów dla neutronów i fotonów gamma (różne składowe powolna/szybka). Parametr PSD: $Q_{slow}/Q_{total}$ (stosunek ładunku "wolnej" składowej do całkowitego ładunku). Neurony mają wyższy $Q_{slow}/Q_{total}$ niż gamma.

W koincydencji podwójnej n-γ: detektor neutronów (NE-213 daje n-γ PSD) w koincydencji z detektorem gamma (NaI lub HPGe). Koincydencja identyfikuje pary n-γ z tej samej reakcji jądrowej (np. reakcja rozszczepienia: n + ²³⁵U → n' + γ + fragmenty). Macierz E_n × E_γ z bramką PSD pozwala na pełną charakteryzację kinetyki reakcji.

Zastosowania: pomiary przekrojów czynnych reakcji neutronowych, badanie struktury jąder przez (n,γ) reakcje, diagnostyka reaktorów, detekcja materiałów rozszczepialnych (NDA — Non-Destructive Assay).

Cyfrowe systemy koincydencji: timestamping i offline matching

Nowoczesne systemy DPP realizują koincydencje inaczej niż klasyczny NIM: zamiast hardware AND gate, używają timestampów i offline matching.

Zasada timestampowania: każde zdarzenie w każdym kanale digitizera otrzymuje timestamp z rozdzielczością 1-10 ns. Zdarzenia są zapisywane do bufora list-mode: (kanał, timestamp, energia). Offline coincidence matching: algorytm przegląda pary zdarzeń z różnych kanałów i sprawdza czy timestampy są bliskie (różnica < Δt). Złożoność: O(N·k), gdzie k to średnia liczba zdarzeń w oknie.

Zalety timestamp matching: elastyczne okno Δt (zmiana w software bez modyfikacji hardware), możliwość retrospektywnej zmiany Δt po zebraniu danych, wielokrotne koincydencje (3-fold, 4-fold) bez dedykowanego hardware, łatwa eliminacja przypadkowych przez shifted window method. Wady: wymaga dokładnej synchronizacji zegarów między kanałami, drift zegarów = systematyczne błędy koincydencji.

Metoda shifted window do pomiaru tła przypadkowego

Elegancka technika pomiaru tła przypadkowego: pomiar ten sam co koincydencji, ale z przesuniętym oknem czasowym.

Zasada: koincydencje prawdziwe — zdarzenia w kanale A i B z różnicą czasu $|t_A - t_B| < \Delta t$. Koincydencje przypadkowe — zdarzenia z różnicą czasu $T < |t_A - t_B| < T + \Delta t$, gdzie $T \gg \Delta t$ (np. $T = 10\,\mu s$). Ponieważ T jest dużo większe niż jakiekolwiek fizyczne korelacje, przesunięte koincydencje zawierają wyłącznie tło przypadkowe. Stanowią więc empiryczną miarę $n_{acc}$ bez założeń.

Implementacja w Python: wczytaj zdarzenia z obu kanałów (list-mode), przesuń timestampy jednego kanału o T, oblicz koincydencje w przesuniętym oknie. Wynik: empiryczna miara n_acc. Metoda shifted window jest standardem w fizyce cząstek i medycynie nuklearnej (PET) — w klinicznym badaniu PET tło przypadkowe może stanowić 20-50% wszystkich koincydencji, stąd korekcja jest niezbędna dla poprawnego obrazowania.

Standardy czasu w fizyce jądrowej: GPS i White Rabbit

Dla eksperymentów wymagających absolutnej synchronizacji czasu (sieci monitoringu jak IMS CTBTO):

GPS timing: odbiornik GPS z wyjściem 1 PPS (Pulse Per Second) — precyzyjny sygnał sekundy ze standardem UTC. Dokładność: ±20-100 ns wobec UTC. Używany w stacjach IMS, laboratoriach rozproszonych geograficznie.

IEEE 1588 PTP (Precision Time Protocol): sieciowy protokół synchronizacji przez Ethernet. Z hardware timestamping + Boundary Clocks: do ~10-100 ns.

White Rabbit (WR): rozszerzenie IEEE 1588 opracowane w CERN. Łączy synchronous Ethernet (SyncE) z IEEE 1588 hardware timestamping i kalibrację opóźnień asymetrycznych. Wynik: sub-nanosekundowa precyzja synchronizacji przez Ethernet. WR jest open hardware (OHWR.org). Użycie w Polsce: NCBJ Świerk.

Porównanie dla zastosowań: koincydencje <1 μs, 1 detektor — kabel koaksjalny; ~100 ns, 2-4 detektory — wspólny zegar FPGA; ~10 ns, wielomodułowy — VME Trigger Bus lub White Rabbit; między laboratoriami — GPS 1PPS + UTC timestamping.

Koincydencje w fizyce reaktorów: metoda Rossi-alfa

W fizyce reaktorów, technika koincydencji stosowana jest w metodzie Rossi-alfa: pomiarze reaktywności reaktora w stanie podkrytycznym.

Metoda: mierzy rozkład $C(\tau)$ — częstość koincydencji par zdarzeń w detektorze neutronów jako funkcja czasu $\tau$ między nimi. Dla reaktora podkrytycznego:

$$C(\tau) = B + A \cdot e^{-|\alpha| \cdot \tau}$$

gdzie $\alpha = (k_{eff}-1)/\Lambda < 0$, $\Lambda$ — czas życia neutronów, B — tło przypadkowe. Zmierzona korelacja pochodzi od neutronów z tych samych lawin rozszczepiennych. Fit wykładniczy do danych daje $|\alpha|$ i stąd reaktywność. Realizacja tradycyjna: detektor neutronów (He-3 lub BF3) + moduł koincydencji NIM z różnymi opóźnieniami τ + zestaw skalerów. Realizacja nowoczesna: digitizer z list-mode + offline autocorrelation.

Historia techniki koincydencji

1929 — Walther Bothe i Hans Geiger: pierwsze elektroniczne układy koincydencji z lampami próżniowymi, do weryfikacji jednoczesności impulsów w detektorach Geigera. Bothe — Nobel z fizyki 1954.

1930-1950: Koincydencja stosowana do badania promieniowania kosmicznego, weryfikacji kwantów gamma w parze, badań kaskad jądrowych.

1960-1980: NIM standard umożliwia standaryzację układów koincydencji. Moduły COINCIDENCE jako standardowe produkty (LeCroy, Ortec). Powszechne zastosowanie w spektroskopii gamma-gamma i fizyce cząstek.

1975: PET (Positron Emission Tomography) — pierwsze kliniczne zastosowanie koincydencji do obrazowania medycznego.

2000-2015: Cyfrowe DPP z listmode umożliwiają offline koincydencje. TOF-PET z LSO i SiPM osiąga sub-nanosekund CRT. Poszukiwania ciemnej materii wymagają ultra-niskich tła — veto kosmiczne + ultraczyste materiały + koincydencje.

2020-teraz: Klasyczna fizyka koincydencji w LZ, XENONnT (ciemna materia) i SNO+, DUNE (neutrina). Kwantowe sieci sensorów ze splątanymi parami fotonów — nowe zastosowania.

Błędy w implementacji koincydencji: typowe pułapki

Z doświadczeń polskich i zagranicznych laboratoriów:

Pułapka 1 — niezrównoważone opóźnienia: efficiency koincydencji jest niska mimo dobrego ustawienia progów. Przyczyna: różnica opóźnień propagacji między kanałami. Diagnostyka: połącz oba wejścia do tego samego generatora impulsów, wyznacz efficiency w funkcji delay. Rozwiązanie: kabel wyrównujący lub Delay Unit NIM.

Pułapka 2 — zbyt szerokie impulsy po CFD: przypadkowe koincydencje wyższe niż obliczone $2\Delta t \cdot n_1 \cdot n_2$. Przyczyna: CFD generuje impulsy szersze niż okno koincydencji. Rozwiązanie: sprawdź szerokość impulsu CFD na oscyloskopie, użyj "fast coincidence" (ECL-based).

Pułapka 3 — saturacja ADC przy koincydencji: widmo energetyczne ścięte lub zniekształcone. Przyczyna: gate do ADC przybywa za wcześnie lub za późno. Rozwiązanie: sprawdź timing gate względem impulsu analogowego.

Pułapka 4 — koincydencje z własnego szumu elektronicznego: koincydencje rejestrowane bez źródła. Przyczyna: szum przesłuchuje się przez wspólne zasilanie. Diagnostyka: wyłącz detektory — czy koincydencje znikają? Rozwiązanie: separacja galwaniczna, filtry na zasilaniu.

Praktyczne wskazówki projektowania układu koincydencji

Dla studenta planującego swój pierwszy układ koincydencji:

Punkt startowy: zacznij od najprostszego układu (2 detektory, 1 moduł koincydencji, 2 skalery). Rozumiej każdy element przed dodaniem następnego.

Wybór okna Δt: zmierz CRT systemu empirycznie (time difference spectrum z ²²Na). Ustaw okno = 2-3 × CRT (FWHM), by zachować >98% prawdziwych koincydencji.

Kalibracja opóźnień: używaj generatora impulsów do wyrównania opóźnień kanałów przed pierwszym pomiarem. Zapisz długości kabli i zmierzone opóźnienia.

Weryfikacja: sprawdź przypadkowe koincydencje (shifted window lub obliczenie z singles rate). Proporcja przypadkowych do prawdziwych powinna być <5%.

Monitoring stabilności: mierz efficiency koincydencji periodycznie podczas długich pomiarów. Zmiany mogą sygnalizować dryfty HV, temperatury lub niestabilności elektryczne.

Te wskazówki, oparte na fizycznych zasadach omówionych w tym artykule, pozwalają zbudować wiarygodny układ koincydencji od podstaw i systematycznie diagnozować ewentualne problemy.

Rozbudowane zastosowania koincydencji: eksperymenty naukowe nowej generacji

Nowoczesne eksperymenty naukowe wychodzą daleko poza prostą koincydencję 2-fold. Przykłady:

KM3NeT (podwodny teleskop neutrin Morza Śródziemnego): Sieć ok. 6000 optycznych modułów (OM) na dnie Morza Śródziemnego (2500-5000 m głębokości). Każdy OM zawiera 31 małych PMT. Rejestracja neutrina: mion Czerenkowa ze zderzenia neutrina z materią dna morza wytwarzający stożek świetlny. Koincydencja: kilka OM musi zarejestrować fotony w krótkim oknie (<200 ns), by odróżnić mion od tła (izolowane fotony z rozpadów ⁴⁰K w wodzie morskiej, blumina bioluminescencji). Timestamp White Rabbit: wspólna synchronizacja wszystkich OM z <1 ns precyzją. Przepływ danych: ~10 GB/s z całej sieci do centrum brzegowego.

XENONnT (poszukiwania ciemnej materii): 5800 kg ciekłego ksenonu jako detektor. Ciemna materia rejestrowana przez pary scyntylacja+jonizacja. Koincydencja: S1 (scyntylacja bezpośrednia, <100 fotonów γ) + S2 (scyntylacja opóźnionej rekombinacji, po dryftie elektronów). Okno koincydencji S1-S2: 1-5 ms (czas dryfu). Veto antykosmiczne: zewnętrzny detektor wody rejestruje miony kosmiczne, weto wyłącza detektor na kilka ms. Antykoincydencja neutronowa: PMT w wodzie + kryształy NaI w zewnętrznych warstwach eliminują rozproszone neutrony.

CUORE (poszukiwania bezmasowego podwójnego beta): 988 kryształów TeO₂ po 750 g każdy, schłodzonych do 10 mK. Rejestracja: wzrost temperatury ~100 μK od pochłoniętego fotonu lub elektronu — odczyt przez SQUID. Koincydencja: zdarzenia 0νββ powinny dawać energię 2528 keV w jednym krysztale. Zdarzenia tła często dają koincydencje w sąsiednich kryształach (wielokrotne rozproszenia). Antykoincydencja sąsiednich kryształów redukuje tło o czynnik ~10. Lista zdarzeń z 988 kanałów, offline multiplet analysis.

Belle II (B-Factory w KEK, Japonia): Badanie naruszeń symetrii CP w rozpadach mezonów B. Koincydencja jest fundamentem całego detektora: system triggera wyselekcjonuje 1 spośród 40 000 zderzeń e+e- na sekundę. Trigger L1 (hardware, <5 μs latencja): koincydencja sygnałów z CDC (drift chamber) + ECL (kalorimetr). Trigger L2/L3 (software): bardziej złożone kryteria na GPU farmie. Polscy fizycy z AGH, UJ, NCBJ uczestniczą w Belle II.

Zastosowania koincydencji w medycynie nuklearnej i radioterapii

Poza PET, technika koincydencji pojawia się w innych zastosowaniach medycznych:

Dozymetria in vivo w radioterapii protonowej: Promieniowanie gamma z anihilacji pozytronów (emitowanych przez aktywowane fragmenty tkanki) monitoruje rozkład dawki. Dwa detektory gamma w koincydencji 511 keV+511 keV identyfikują miejsce interakcji wiązki protonowej z tkanką. Metoda: PGS (Prompt Gamma Spectroscopy) + PET inline. Prowadzone badania w NCBJ i IFJ PAN.

SPECT/CT z technetowym ⁹⁹ᵐTc: Nie jest to koincydencja, lecz kolimacja. Ale logika akwizycji jest podobna: bramkowanie EKG (rejestracja obrazu SPECT serca zsynchronizowana z fazą cyklu sercowego) dla redukcji artefaktów ruchu.

Radiochirurgia stereotaktyczna (Gamma Knife, CyberKnife): Wiele wiązek gamma przecinających się w punkcie docelowym. Koincydencja geometryczna: promieniowanie różnych koliatorów zbiegające się w jednym punkcie. Zasada analogiczna do fizycznej koincydencji: tylko gdzie "wszystkie wiązki są jednocześnie" następuje ablacja tkanki nowotworowej.

Koincydencje a pomiary bezwzględnej aktywności

Jednym z ważnych zastosowań koincydencji jest wyznaczanie bezwzględnej aktywności źródeł promieniotwórczych bez konieczności kalibrowanego detektora (metoda 4π-β-γ):

Zasada: izotop emituje zarówno cząstkę beta jak i kwant gamma przy każdym rozpadzie (np. ⁶⁰Co: beta + 2 gamma; ¹³⁷Cs: beta + gamma 662 keV). Detektor beta (4π — otaczający źródło) rejestruje każdy rozpad z wydajnością ε_β. Detektor gamma rejestruje gamma z wydajnością ε_γ.

Zliczenia: $n_\beta = A \cdot \varepsilon_\beta$, $n_\gamma = A \cdot \varepsilon_\gamma$, $n_{coinc} = A \cdot \varepsilon_\beta \cdot \varepsilon_\gamma$ (przy założeniu korelacji beta-gamma).

Aktywność:

$$A = \frac{n_\beta \cdot n_\gamma}{n_{coinc}}$$

Niezwykła właściwość: wynik nie zależy od efektywności żadnego z detektorów! Aktywność bezwzględna wyznaczona tylko ze stosunku zliczeń. To jest zasada metody 4πβ-γ koincydencji, która jest podstawą wzorcowania aktywności przez krajowe laboratoria metrologii radiacyjnej (PTB w Niemczech, NPL w UK, CCRI w BIPM, NCBJ w Polsce).

Techniczne wymagania: detektor beta musi mieć efektywność ε_β jak największą (możliwie blisko 4π steradianów), a zależność A(ε_β, ε_γ) musi być liniowa. W praktyce dokonuje się ekstrapolacji do ε_β = 1.

Eksperyment dydaktyczny: wyznaczanie półczasu przez koincydencje

Dla laboratoriów z dostępem do krótkożyjącego izotopu (np. Ag-108m, T₁/₂ = 418 min; lub aktywacja neutronowa srebra: Ag-108 T₁/₂ = 2.4 min, Ag-110 T₁/₂ = 24.6 s):

Schemat ćwiczenia:

Napromieniuj próbkę srebrną przez 5 minut wiązką neutronów termicznych (jeśli dostępne) lub użyj gotowego źródła ¹¹¹In (T₁/₂=2.8 d, emituje koincydentne fotony gamma).
Ustaw układ koincydencji: 2 detektory NaI po obu stronach próbki, SCA na 245+171 keV piki ¹¹¹In, moduł koincydencji, scaler.
Mierz koincydencje i singles każde 5 minut przez kilka godzin.
Dopasuj wykładniczy zanik: $n(t) = n_0 \cdot e^{-\lambda t}$, wyznacz $\lambda$ i $T_{1/2} = \ln 2 / \lambda$.
Porównaj z wartością tabelaryczną: 2.8047 d dla ¹¹¹In.

Ćwiczenie łączy: technikę koincydencji (budowanie układu, kalibracja), statystykę pomiarową (błąd Poissona) i kinetykę promieniotwórczą (dopasowanie wykładnicze). Dla ¹¹¹In mierzenie przez kilka godzin daje ~50% zmiany zliczalności — dostatecznie dużo dla dokładnego wyznaczenia T₁/₂ bez wielodniowego pomiaru.

Literatura zalecana do pogłębienia wiedzy

Kluczowe pozycje dla tematyki koincydencji:

Książki:

G.F. Knoll, "Radiation Detection and Measurement" (4. wyd., 2010) — rozdz. 17 "Timing methods and systems": wyczerpujące omówienie CRT, CFD, PSD, PET. Klasyczny podręcznik fizyki eksperymentalnej.
K. Debertin, R.G. Helmer, "Gamma- and X-Ray Spectrometry with Semiconductor Detectors" (1988): szczegółowe omówienie koincydencji sumowania i korekcji.
H. Spieler, "Semiconductor Detector Systems" (2005): elektronika detektorów, timing, szum i koincydencje z perspektywy projektowania układów.

Artykuły przeglądowe:

C.L. Melcher, J.S. Schweitzer, "Cerium-doped Lutetium Oxyorthosilicate: A Fast, Efficient New Scintillator", IEEE Trans. Nucl. Sci. 39(4), 1992 — artykuł opisujący LSO używany w TOF-PET.
M. Conti, B. Bendriem, "The new opportunities for high time resolution clinical TOF PET", Clin. Transl. Imaging, 2019 — przegląd TOF-PET i CRT.

Online:

CERN Online Technical Notes: dokumentacja systemów trigger/DAQ dużych eksperymentów.
XIA LLC Application Notes: praktyczne wskazówki dla systemów Pixie z DPP koincydencji.

Znajomość tych zasobów umożliwia doktorantowi samodzielne pogłębienie wiedzy zarówno w aspekcie teoretycznym (CRT, PSD, PET rekonstrukcja) jak i praktycznym (implementacja hardware, oprogramowanie DAQ, kalibracja).

Głęboka analiza przypadkowych koincydencji i ich konsekwencje

Przypadkowe koincydencje (ang. random coincidences, accidentals) są nie tylko szumem tła — ich właściwa charakterystyka jest niezbędna do poprawnej interpretacji wyników. Kilka aspektów wartych pogłębienia:

Rozkład statystyczny przypadkowych: Przy stałych rate n₁ i n₂ niezależnych detektorów, liczba przypadkowych koincydencji w czasie T ma rozkład Poissona ze średnią $\mu = 2\Delta t \cdot n_1 \cdot n_2 \cdot T$. Odchylenie standardowe: $\sigma = \sqrt{\mu}$. Dla małych rate (μ << 1 na interwał Δt), rozkład Bernoulliego z prawdopodobieństwem $p = 2\Delta t \cdot n_1 \cdot n_2$.

Korelacje między singletons a koincydencjami: W eksperymencie fizycznym, singles rate i coincidence rate nie są niezależne — oba zmieniają się z aktywnością źródła, tłem kosmicznym, zmianami geometrii. Przy analizie danych należy to uwzględnić: nie zakładaj stałego tła przypadkowego jeśli singles rate się zmienia.

Trzecia kategoria: "scattered coincidences": W PET i spektrometrii, obok prawdziwych i przypadkowych, istnieją "scattered coincidences" — fotony, które uległy rozproszeniu Comptona i zmieniły kierunek, ale nadal zostały zarejestrowane w oknie energetycznym. W PET: błędna lokalizacja LOR. Korekcja: energia-based scatter correction, model-based scatter estimation.

Pomijaność przypadkowych a dynamika pomiarów: Przy bardzo niskich rate zdarzeń (np. poszukiwania ciemnej materii, rzadkich rozpadów), przypadkowe koincydencje mogą być dominującym tłem. Projekty takie jak LZ (ciemna materia, 10 t ciekłego ksenonu) osiągają rate zdarzeń kandydatów poniżej 1 zdarzenia na 10 t·rok — przy takiej skali nawet $n_{acc} = 0.01$ zdarzenia/rok jest znaczące.

Gating w akceleratorowych eksperymentach jądrowych

Specyficznym zastosowaniem bramkowania są eksperymenty przy akceleratorach, gdzie wiązka cząstek jest dostarczana w pulsach o określonej strukturze czasowej:

CERN PS/SPS puls wiązki: Protonowa wiązka SPS (Super Proton Synchrotron) dostarcza paczki protonów przez ~10 μs co kilkanaście sekund (spill). Eksperymenty NA61/SHINE i NA62 rejestrują zdarzenia tylko podczas spilla. Bramkowanie "beam gate": zewnętrzny sygnał z systemu kontroli akceleratora (BCM — Beam Current Monitor) otwiera DAQ podczas spilla, zamyka między spillami. Cel: eliminacja tła kosmicznego i szumu elektronicznego pomiędzy spillami.

LHC bunch structure: LHC dostarcza wiązkę w paczkach (bunch) co 25 ns (40 MHz). Każde zdarzenie w detektorach ATLAS/CMS/LHCb/ALICE jest skojarzone z numerem cross-section (BC# — Bunch Crossing Number). Trigger L1 Hardware (latencja <2.5 μs) taguje zdarzenie BC# i decyduje o zapisaniu w 0.25 μs. Bramkowanie: każdy trigger decision zawiera w metadanych BC# dla korekty pile-up i identyfikacji tła.

Tandemowa bramkowanie (multiple gating): Eksperymenty z laserowymi wiązkami pulsowymi (laser pump-probe, photonics research) mają ultra-krótkie pulsy (femtosekund-y). Detektor otwiera się tylko podczas pulsu (okno pikosekund). Po pulsie: detektor zamknięty by uniknąć tła fluorescencji tła. Realizacja: Gating przez Pockels Cell lub elektro-optyczny modulator z synchronizacją do repetycji lasera.

Zastosowanie koincydencji w badaniach materiałowych

Technika koincydencji ma zastosowania nie tylko w fizyce cząstek i medycynie, ale też w badaniach materiałowych:

Pozytronowa spektroskopia anihilacyjna (PAS): Pozytrony z radioaktywnego źródła (²²Na) wnikają do badanego materiału i anihilują z elektronami. Dwa fotony 511 keV rejestrowane w koincydencji. Parametry mierzone:

PALS (Positron Annihilation Lifetime Spectroscopy): czas między emisją pozytonu a anihilacją. Krótki czas → gęsty materiał. Długi czas → wakancje, pory, defekty krysztalu. Koincydencja: foton gamma towarzyszący emisji pozytonu z ²²Na (1274 keV) + foton gamma z anihilacji (511 keV). Różnica czasu = czas życia pozytronu w materiale (100-5000 ps).
ACAR (Angular Correlation of Annihilation Radiation): fotony 511 keV nie są dokładnie antyparallelne — kąt odchylenia koreluje z pędem elektronu w punkcie anihilacji. Mierzony przez dwa detektory obracające się wokół próbki. Rozkład kąta Δθ = informacja o strukturze elektronowej.
CDB (Coincidence Doppler Broadening): rozszerzenie linii 511 keV spowodowane pędem elektronu (efekt Dopplera). Koincydencja dwóch HPGe mierzących energie obu fotonów 511 keV. Suma energii: zawsze blisko 1022 keV. Różnica energii: związana z pędem elektronu w materiale.

PAS jest stosowana do badania: defektów w metalach i stopach (po napromieniowaniu lub deformacji plastycznej), wolnych objętości w polimerach i materiałach porowatych, starzenia materiałów reaktorowych (embryttlement — kruchość pod wpływem neutronów). W Polsce: PAS laboratorium w IFJ PAN Kraków i na Wydziale Fizyki UW.

Koincydencje w detekcji materiałów jądrowych: NDA

Non-Destructive Assay (NDA) — badanie składu radioaktywnego bez naruszania próbki — intensywnie używa technik koincydencji. Kluczowe metody:

Multiplicity Counting: Jądra roszczepialne (²³⁵U, ²³⁹Pu) wykazują łańcuchowe reakcje rozszczepienia — jedna neutron może spowodować emisję 2-3 kolejnych. Statystyka wielokrotności emisji neutronów zależy od masy ²³⁹Pu. Pomiar przez detektor He-3 multiplicity counter (wielokrotnie bardziej niż 2-fold): analiza singletów, dubletów i tripletów (zdarzenia z 1, 2, 3 neutronami koincydentnie).

PDET (Passive Differential Die-Away Technique): Mierzy właściwości mnożące masy rozszczepialnej przez korelację neutronową. Dla kadetów nuklearnych kontroli proliferacji.

Prompt Fission Neutron Correlation: Przy każdym rozszczepieniu spontanicznym Pu-240 lub U-238 emitowane jest 2-3 szybkich neutronów w koincydencji (w ciągu <1 μs). Mierząc koincydencje neutronowe, IAEA weryfi

kuje deklarowaną masę plutonu w paliwie jądrowym. Technologia używana przez inspektorów IAEA do weryfikacji klauzul NPT.

Granica koincydencji i kwantowe zjawiska nielokalne

Na koniec wątek fundamentalny: w kwantowej mechanice koincydencja ma głębsze znaczenie niż klasyczna zbieżność czasowa.

Nierówności Bella i eksperymenty Aspecta: Alain Aspect (Nobel z fizyki 2022) wykazał eksperymentalnie korelacje par splątanych fotonów przekraczające granicę klasyczną (nierówność Bella). Koincydencja: dwa fotony ze splątanej pary rejestrowane w detektorach polaryzacyjnych. Korelacja mierzona przy różnych ustawieniach polaryzatorów. Wynik: korelacje silniejsze niż możliwe klasycznie — dowód splątania kwantowego.

Klasyczna fizyka koincydencji (jak opisana w tym artykule) zakłada lokalne ukryte zmienne — model klasyczny. Eksperymenty Aspecta obalają ten model. Jednak dla wszystkich praktycznych zastosowań (PET, spektroskopia, monitoring reaktorów) klasyczna koincydencja działa doskonale, bo mierzone są makroskopowe zdarzenia a nie stany kwantowe.

Zrozumienie tej granicy między klasycznym a kwantowym podejściem do koincydencji jest ważne dla doktoranta fizyki — uświadamia że technika pomiarowa (NIM, DPP, timestamp matching) jest narzędziem, które mierzy pewien aspekt rzeczywistości, a nie samą rzeczywistość.

Ćwiczenia kontrolne

Oblicz jakościowo, co stanie się z liczbą koincydencji przypadkowych po dwukrotnym zwiększeniu szerokości okna.
Narysuj układ dwóch detektorów z bramką koincydencji i osobny układ z veto.
Wyjaśnij różnicę między opóźnieniem kabla a jitterem detektora.

Przejdź do ćwiczenia interaktywnego

Powiązane materiały

Przypisy

Wikipedia, Nuclear electronics, https://en.wikipedia.org/wiki/Nuclear_electronics. Zakres wykorzystania: elementy torów: fast logic, coincidence i veto gates.
↩
U.S. DOE, Standard NIM Instrumentation System, 1990, https://inis.iaea.org/records/5sq70-myv38/files/21080298.pdf?download=1. Zakres wykorzystania: ograniczenia sygnałów na busach, fast signals i standardowe poziomy logiczne.
↩