Czas martwy, pile-up i gubienie impulsów w torach zliczających

Streszczenie

Czas martwy, pile-up i utracone impulsy są praktycznymi ograniczeniami torów zliczających. Najważniejsze jest zrozumienie, dlaczego przy dużej częstości zdarzeń licznik może zaniżać wynik albo deformować widmo.

Każdy detektor promieniowania zamienia pojedyncze zdarzenie fizyczne w impuls elektryczny. Ten impuls ma skończony czas trwania, musi zostać wzmocniony, ukształtowany, rozpoznany przez elektronikę i zapisany przez licznik albo analizator. Jeśli w tym czasie przyjdzie następne zdarzenie, tor pomiarowy może go nie zauważyć, połączyć oba impulsy w jeden albo nadać im błędną amplitudę.^1,2

Dlatego czas martwy nie jest drobnym detalem aparaturowym. Jest granicą, która decyduje, czy liczba zliczeń nadal opisuje aktywność próbki, czy już głównie ograniczenia elektroniki. W spektrometrii dochodzi drugi problem: pile-up, czyli nakładanie się impulsów, które może tworzyć fałszywe zdarzenia energetyczne i zniekształcać piki widmowe.^1,3

Rozszerzenie tematu

Impuls nie jest punktem w czasie

Najprostszy opis promieniowania mówi: cząstka albo foton trafia w detektor, a licznik dodaje jedno zliczenie. Dydaktycznie to dobry początek, ale technicznie jest to duże uproszczenie. W prawdziwym torze pomiarowym zdarzenie tworzy ładunek, światło scyntylacyjne albo lawinę elektronową. Potem układ wejściowy zamienia to na impuls napięciowy, wzmacniacz nadaje mu użyteczny poziom, układ kształtujący poprawia stosunek sygnału do szumu, a dyskryminator albo analizator wielokanałowy decyduje, czy impuls zostanie zaliczony.

Każdy z tych etapów trwa skończony czas. Impuls może mieć strome narastanie, długi ogon albo kształt specjalnie dobrany przez elektronikę analogową lub cyfrową, aby poprawić rozdzielczość energetyczną. W tym czasie tor nie jest idealnie gotowy na następne zdarzenie. Jeżeli druga cząstka przyjdzie za szybko, układ może być nadal zajęty poprzednią. To właśnie intuicyjny sens czasu martwego.

W liczniku Geigera-Müllera problem jest bardzo widoczny, bo po wyładowaniu w rurze potrzeba chwili, aby gaz i pole elektryczne wróciły do stanu umożliwiającego kolejne poprawne zliczenie. W licznikach scyntylacyjnych i półprzewodnikowych problem jest inny: detektor może być szybszy, ale elektronika kształtująca, dyskryminator, ADC albo analizator widma nadal mają skończony czas przetwarzania. Dlatego czas martwy jest pojęciem całego toru, nie tylko samego detektora.^1,2

Dwa modele czasu martwego

W praktyce używa się dwóch prostych modeli: nieparalizowalnego i paralizowalnego. Model nieparalizowalny zakłada, że po zarejestrowaniu impulsu układ jest martwy przez czas tau, a zdarzenia przychodzące w tym okresie są po prostu tracone. Nie wydłużają jednak martwego okresu. Dla prawdziwej częstości zdarzeń n i obserwowanej częstości m typowa zależność ma postać:

m = n / (1 + n tau)

Po odwróceniu można oszacować poprawioną częstość:

n = m / (1 - m tau)

Ten wzór od razu pokazuje niebezpieczną granicę. Jeżeli m tau zbliża się do jedności, poprawka przestaje być stabilna. Mały błąd pomiaru czasu martwego albo częstości zliczeń daje wtedy ogromny błąd aktywności.

Model paralizowalny zakłada coś bardziej surowego: zdarzenie, które przychodzi w czasie martwym, nie zostaje zliczone, ale ponownie „paraliżuje” układ i wydłuża okres niewrażliwości. W uproszczeniu obserwowana częstość może być opisana zależnością:

m = n exp(-n tau)

Ten model ma ważną konsekwencję jakościową. Przy bardzo dużej aktywności licznik może zacząć pokazywać mniej zliczeń, a nie więcej. To szczególnie zdradliwe, bo niedoświadczony użytkownik może uznać niski odczyt za niski poziom promieniowania, podczas gdy układ został przeciążony.

Czas martwy a pile-up

Czas martwy dotyczy przede wszystkim pytania: ile zdarzeń zostało utraconych? Pile-up dotyczy pytania: co elektronika zrobiła z impulsami, które przyszły zbyt blisko siebie? W prostym liczniku dwa bliskie impulsy mogą zostać policzone jako jeden. W spektrometrze gamma problem jest subtelniejszy, bo amplituda impulsu niesie informację o energii fotonu. Jeśli dwa impulsy częściowo się nałożą, analizator może zobaczyć jeden impuls o błędnej amplitudzie.

W skrajnym przypadku suma dwóch zdarzeń może trafić do kanału energetycznego, w którym fizycznie nie było pojedynczego fotonu o takiej energii. Widmo zaczyna wtedy zawierać artefakty aparaturowe. Przy dużych aktywnościach rośnie tło ciągłe, poszerzają się piki, pojawiają się zdarzenia sumacyjne, a automatyczna identyfikacja radionuklidów staje się mniej wiarygodna.

To tłumaczy, dlaczego spektrometria nie polega na ustawieniu próbki jak najbliżej detektora. Więcej zliczeń poprawia statystykę tylko do pewnej granicy. Po jej przekroczeniu rośnie udział efektów aparaturowych: czasu martwego, pile-upu, niestabilności wzmocnienia i przeciążenia toru.

Jak widać problem na oscyloskopie

Oscyloskop pozwala zobaczyć to, co licznik ukrywa w jednej liczbie. Przy małej częstości impulsy są od siebie dobrze oddzielone. Można łatwo rozpoznać czas narastania, amplitudę, czas opadania i poziom dyskryminacji. Gdy częstość rośnie, odstępy między impulsami stają się porównywalne z czasem ich trwania. Wtedy pojawia się nakładanie, zmiana linii bazowej i niestabilność progu.

W torze z dyskryminatorem problem można ująć tak: dyskryminator widzi tylko przekroczenie progu. Jeżeli drugi impuls zaczyna się na ogonie pierwszego, może przekroczyć próg wcześniej, później albo wcale. W torze spektrometrycznym dochodzi jeszcze pomiar amplitudy. Układ może próbować zmierzyć wysokość impulsu, który nie zaczął się od prawidłowej linii bazowej.

Dlatego artykuł o oscyloskopach w pomiarach nanosekundowych powinien być czytany razem z tym tekstem. Czas martwy i pile-up są zjawiskami licznikowymi, ale ich źródło widać właśnie w kształcie przebiegu, paśmie toru, terminacji i czasie odpowiedzi elektroniki.

Metoda dwóch źródeł jako idea dydaktyczna

Klasyczny sposób wyznaczania czasu martwego licznika GM wykorzystuje porównanie zliczeń od dwóch źródeł mierzonych osobno i razem. Sens jest prosty: gdyby układ nie miał czasu martwego, częstość dla dwóch źródeł razem byłaby sumą częstości dla każdego z nich, po poprawce na tło. Jeżeli wynik wspólny jest mniejszy od sumy, część impulsów została utracona.²

W serwisie najlepiej przerobić tę metodę na ćwiczenie z danymi syntetycznymi. Student dostaje cztery liczby: tło, źródło A, źródło B oraz A+B. Następnie oblicza częstości netto, porównuje sumę z pomiarem wspólnym i wyznacza przybliżony czas martwy. Nie trzeba do tego żadnego realnego źródła promieniotwórczego, a cała wartość dydaktyczna zostaje zachowana.

Taka forma ma jeszcze jedną zaletę. Można wygenerować kilka zestawów danych: jeden prawie idealny, jeden z małym czasem martwym i jeden z przeciążeniem. Wtedy ćwiczenie nie jest podstawieniem do wzoru, ale rozpoznaniem sytuacji pomiarowej.

Jak ograniczać błąd

Najprostsza zasada brzmi: nie pracować przy zbyt dużej częstości zliczeń. Jeżeli tor pokazuje duży udział czasu martwego, lepiej zwiększyć odległość od źródła, użyć kolimatora, skrócić geometrię widzenia albo dobrać mniej aktywną próbkę. W spektrometrii gamma często lepiej mierzyć dłużej przy mniejszej częstości niż krótko przy przeciążonym torze.

Druga zasada to zapisywanie informacji o czasie martwym razem z widmem. Współczesne analizatory często pokazują live time i real time, czyli czas efektywnie dostępny dla zliczeń oraz rzeczywisty czas pomiaru. Różnica między nimi jest pierwszym ostrzeżeniem, że wynik wymaga korekty.

Trzecia zasada dotyczy doboru kształtowania impulsu. Dłuższy czas kształtowania może poprawiać rozdzielczość energetyczną, ale zwiększa ryzyko pile-upu przy dużej częstości. Krótsze kształtowanie bywa lepsze dla dużej liczby zdarzeń, ale może pogarszać rozdzielczość. Nie ma jednego ustawienia najlepszego dla wszystkich pomiarów.

Dlaczego to jest temat elektroniki

Czas martwy bywa opisywany w podręcznikach radiometrii, ale w praktyce jest to także temat elektroniki impulsowej. Decyduje o nim detektor, pojemności wejściowe, przedwzmacniacz, układ kształtujący, dyskryminator, przetwornik ADC, bufor danych i oprogramowanie akwizycji. Każdy z tych elementów może stać się wąskim gardłem.

To dobry przykład, dlaczego serwis potrzebuje osobnej kategorii Elektronika. Bez zrozumienia czasu martwego łatwo traktować licznik promieniowania jak zwykły licznik zdarzeń. Tymczasem przyrząd pomiarowy jest dynamicznym układem, który sam wprowadza ograniczenia do obserwowanego zjawiska.

Najkrótszy wniosek jest taki: liczba zliczeń jest wynikiem współpracy promieniowania z elektroniką. Im większa częstość zdarzeń, tym mniej wolno udawać, że elektronika jest przezroczysta. Każdy tor pomiarowy ma swoją specyficzną charakterystykę czasu martwego — model fizyczny (paralizowalny lub nieparalizowalny lub hybryda jako kombinacja obu), wartość τ w μs, odporność na pile-up i sposób korekty przy rejestrowanym live time. Poznanie tych parametrów dla używanego systemu jest równie ważne jak poznanie charakterystyki samego detektora (rozdzielczość energetyczna, wydajność geometryczna, tło własne). Bez tej wiedzy wyniki spektrometrii można mieć, ale nie można ich w pełni rozumieć, prawidłowo dokumentować ani prawidłowo oceniać ich niepewności pomiaru i wynikających z nich granic wykrywalności.

Historia pomiaru czasu martwego i zliczania impulsów

Pierwsze elektroniczne układy zliczające pojawiły się razem z licznikiem Geigera-Müllera w 1928 roku. Oryginalna rura GM miała naturalny czas martwy rzędu 100–300 μs wynikający z jonizacji gazu i quenchingu. Liczniki mechaniczne (rejestratory elektromagnetyczne) wprowadzały dodatkowe opóźnienie — każde zliczenie wymagało przesunięcia kółka licznikowego, co zajmowało kilka milisekund. Przy typowych źródłach laboratoryjnych epoki (Rn, Ra) nie był to problem.

Dopiero w latach 40. XX w., przy pomiarach w reaktorach i akceleratorach, gdzie częstości zdarzeń przekraczały 10⁶ s⁻¹, elektronika pomiarowa stała się wąskim gardłem. Bell Telephone Laboratories i MIT Radiation Laboratory opracowały w latach 1942–1945 pierwsze skalery binarne (scaling circuits) — układy, które zliczały co 2, 4, 8 lub 16 zdarzeń, redukując wymaganą szybkość układu zliczającego. Szczegóły opisano w „Radiation Laboratory Series" — wielotomowej pracy zebranych przez MIT w 1947 roku, stanowiącej podstawę całej powojennej elektroniki jądrowej.

W latach 50. pojawiły się pierwsze analizatory wielokanałowe (Multichannel Analyzer, MCA). Każde zdarzenie było digitalizowane (mierzono amplitudę impulsu proporcjonalną do energii fotonu) i dodawane do odpowiedniego kanału histogramu. Wymagało to ADC i logiki zapisu, co tworzyło nowy rodzaj czasu martwego — nie fizyczny czas wyładowania detektora, lecz czas przetwarzania elektronicznego. Pierwsze MCA miały czas martwego na zdarzenie rzędu 10–50 μs, co ograniczało częstości mierzalnych zdarzeń do kilkudziesięciu tysięcy na sekundę.

Postęp w technice półprzewodnikowej skrócił czasy martwego do mikrosekund i poniżej. Nowoczesne cyfroweprzetworniki impulsów (DPP, Digital Pulse Processor) mogą przetwarzać miliony zdarzeń na sekundę. Jednocześnie pojawiły się nowe metody kompensacji czasu martwego, opisane poniżej.

Formalny opis modeli: obliczenia i przykłady

Model nieparalizowalny (Non-Extendable Dead Time)

Zakładamy stały czas martwy τ po każdym zaliczonym zdarzeniu. Prawdziwa częstość zdarzeń n (źródłowa) jest związana z obserwowaną częstością m wzorem:

m = n / (1 + n·τ)

Przykład 1: n = 10 000 s⁻¹, τ = 10 μs (10⁻⁵ s)

• n·τ = 0,1
• m = 10 000 / 1,1 = 9 091 s⁻¹
• Utrata: ~9%, mieszcząca się w granicy korekcji

Przykład 2: n = 100 000 s⁻¹, τ = 10 μs

• n·τ = 1,0
• m = 100 000 / 2 = 50 000 s⁻¹
• Utrata: 50% — wynik wymaga korekty, ale korekcja jest niestabilna

Odwrócona formuła: n = m / (1 - m·τ). Dla m = 50 000 s⁻¹ i τ = 10 μs: n = 50 000 / (1 - 0,5) = 100 000 s⁻¹ — wniosek poprawny. Ale przy m = 90 000 s⁻¹ i τ = 10 μs: n = 90 000 / 0,1 = 900 000 s⁻¹ — mały błąd pomiaru m lub τ daje ogromny błąd n. Dlatego model nieparalizowalny daje stabilne korekty tylko dla m·τ < 0,2.

Model paralizowalny (Extendable Dead Time)

Tu każde zdarzenie, nawet niezaliczone, wydłuża czas martwy:
m = n · exp(-n·τ)

Przy n = 100 000 s⁻¹ i τ = 10 μs: m = 100 000 · e⁻¹ = 36 788 s⁻¹. To mniej niż w modelu nieparalizowalnym przy tej samej aktywności.

Krytyczna własność: maksymalne m_max = 1/(e·τ) — przy τ = 10 μs wynosi 36 788 s⁻¹. Powyżej tego maksimum m zaczyna maleć przy rosnącym n. Układ paralizowalny przy n → ∞ daje m → 0. To niebezpieczna własność liczników GM przy dużych aktywności — „wyłączają się" bez ostrzeżenia.

Który model obowiązuje?

Większość liczników GM i scyntylacyjnych przy prostym dyskryminatorze zachowuje się bliżej modelu paralizowalnego. Liczniki z aktywnym reszetem (ang. forced reset po czasie τ) zachowują się nieparalizowalnie. Cyfrowe przetworniki sygnałów (DPP) mogą implementować oba typy, a nawet hybrydę. W praktyce zawsze warto sprawdzić dane producenta albo zmierzyć eksperymentalnie metodą dwóch źródeł.

Elektronika kształtowania impulsów a czas martwy i pile-up

Kształt impulsu wychodzącego z przedwzmacniacza determinuje zachowanie całego toru. Istnieją trzy główne podejścia:

Kształtowanie gaussowskie (CR-RC):
Układ CR-RC lub CR-RC^n tworzy impuls o symetrycznym kształcie dzwonowym. Czas kształtowania (peaking time) to czas do osiągnięcia maksimum amplitudy — typowo 0,5–10 μs dla spektrometrii HPGe. Dłuższy czas kształtowania: lepsza rozdzielczość energetyczna (szum szeregowy jest niższy), ale więcej pile-upu przy dużych częstościach. Krótszy czas: gorsza rozdzielczość, ale zdolność do pracy przy wyższych zliczeniach.

Kształtowanie trapezoidalne (Digital Trapezoidal):
W cyfrowych przetwarznikach sygnałów (DPP-PHA, ang. Pulse Height Analysis) stosuje się kształtowanie trapezoidalne: impuls prostokątny z krótkimi rampami narastania i opadania. Parametry: rise time, flat top, fall time — dobierane niezależnie. Pozwala na krótki flat top (pomiar szczytowy) i krótki fall time, skracając efektywny czas martwego bez rezygnacji z dobrego pomiaru amplitudy.

Kształtowanie bipolarne:
Impuls bipolarny (dodatni + ujemny) skraca czas powrotu do linii bazowej — odbywa się samoistnie przez zerowanie całki impulsu. Stosowany w licznikach szybkich (np. w fizyce cząstek). Wadą jest większy szum i mniejsza liniowość pomiaru amplitudy w spektrometrii energetycznej.

Układy odrzucania pile-upu (PUR)

Pile-Up Rejection (PUR) to elektroniczny system, który rozpoznaje nakładające się impulsy i odrzuca je zamiast rejestrować z błędną energią. Działanie podstawowe:

1. Układ szybki (fast channel) — dyskryminator o niskim progu, rejestruje każde przekroczenie progu, niezależnie od toru spektrometrycznego. Jego czas martwy jest bardzo krótki (ns).
2. Tor wolny (slow channel, główny spektrometr) — mierzy amplitudę impulsu ze starannym kształtowaniem, ale jest wolniejszy.
3. Jeżeli układ szybki wykryje drugie zdarzenie podczas przetwarzania pierwszego przez tor wolny, sygnał PUR odrzuca bieżący impuls i nie zapisuje go do histogramu.

Wynik: widmo zawiera mniej zdarzeń (zdarzenia odrzucone przez PUR wchodzą do czasu martwego), ale te, które są zapisane, mają prawidłowe energie. Alternatywnie, systemy z VPG (Virtual Pulse Generator) rejestrują odrzucone zdarzenia jako znany rodzaj korekty statystycznej.

Nowoczesne cyfrowe systemy DPP integrują PUR bezpośrednio w logice FPGA. Czas martwy na zdarzenie wynosi od 100 ns do kilku mikrosekund w zależności od konfiguracji, a odrzucanie pile-upu odbywa się cyfrowo przez analizę kształtu impulsu (pulse shape discrimination, PSD).

Straty bezstratne (Loss-Free Counting, LFC) i metoda VPG

W pomiarach środowiskowych i bezpieczeństwa jądrowego nieraz mierzymy próbki o szybko zmieniającej się aktywności (wczesne opady, monitoring reaktora). Standardowe podejście (korekta z live time) ma dwa ograniczenia: działa tylko dla stacjonarnych aktywności i traci dokładność przy dużym czasie martwym.

Metoda VPG (Virtual Pulse Generator) rozwiązuje oba problemy:

1. Obok normalnych impulsów z detektora do analizatora podawane są regularne, znane impulsy z generatora pulsów (o stałej amplitudzie i częstości).
2. Każde zdarzenie generatora, które zostało zagubione przez czas martwy, jest policzalne — porównujemy wliczone pulsy generatora z oczekiwaną liczbą.
3. Dla każdego kanału histogramu stosujemy proporcjonalną korektę: N_korekt = N_zaobs × (f_gen × T) / N_gen_zaobs.

Gdzie f_gen × T to oczekiwana liczba impulsów generatora przy częstości f_gen przez czas T, a N_gen_zaobs to faktycznie zliczone impulsy generatora.

VPG automatycznie koryguje czas martwy niezależnie od tego, czy zmienia się w czasie czy nie. W polskim rynku aparatury jądrowej metodę VPG oferuje m.in. oprogramowanie GENIE 2000 (Mirion Technologies) i code GammaVision (ORTEC/AMETEK). IFJ PAN Kraków stosuje ją standardowo w pomiarach środowiskowych o szybko zmieniającym się tłe.

Sumowanie kaskadowe jako błąd podobny do pile-upu

Sumowanie kaskadowe (true coincidence summing, TCS) to efekt różny od pile-upu, ale o podobnych konsekwencjach widmowych. Dotyczy nuklidów, które emitują dwie lub więcej kwant gamma w kaskadzie z krótkim czasem między emisjami (< 1 ns).

Przykład: ⁶⁰Co emituje dwa fotony kaskadowo: 1173 keV i 1332 keV. Oba mogą jednocześnie wejść do detektora (szczególnie przy bliskiej geometrii), tworząc jedno zdarzenie o energii sumarycznej 2505 keV. W widmie pojawia się wówczas nieprawidłowy pik sumacyjny, a piki 1173 i 1332 keV tracą zliczenia.

Korekcja TCS zależy od geometrii (efektywność pełnego pochłaniania, total efficiency) i nie zależy od aktywności ani czasu martwego. Przy wydłużaniu czasu kształtowania (większy okienko czasowe pomiaru) TCS rośnie — co jest efektem odwrotnym do przypadku pile-upu. Dlatego:

• Pile-up jest błędem statystycznym: rośnie z aktywnością, spada przy większej odległości geometrycznej.
• TCS jest błędem geometrycznym: rośnie z bliskością geometryczną, niezależnie od aktywności.

Rozróżnienie jest kluczowe w kalibracji wydajności detektora dla nuklidów kaskadowych (⁶⁰Co, ¹³³Ba, ¹⁵²Eu).

Czas martwy w różnych typach detektorów

Detektor	Typowy czas martwy	Główny mechanizm	Uwagi
Licznik GM	100–300 μs	Wyładowanie gazu	Model paralizowalny
Detektor proporcjonalny	1–10 μs	Kolekcja jonów	Quenching elektroniczny
Scyntylator NaI(Tl)	0,5–2 μs	ADC + zapis MCA	Zależy od toru cyfrowego
HPGe (analog)	5–20 μs	Kształtowanie CR-RC	Zależy od shaping time
HPGe (DPP cyfrowy)	0,5–3 μs	Trapezoidalny filtr FPGA	Nowoczesne systemy
SiPM + DPP	0,1–1 μs	Czas powrotu lawinowego	Krótki, ale pile-up SiPM specyficzny
Detektor foliowy (dozymetria)	Stały odczyt	Nie dotyczy zliczania	Prąd, nie impulsy

Licznik GM przy dużych dawkach może „wyłączyć się" (utrata zliczeń jak model paralizowalny) — klasyczny błąd ratownictwa radiologicznego, gdy licznik przy bardzo dużym polu zwraca niski odczyt, co jest błędnie interpretowane jako brak zagrożenia. Standard IAEA Safety Reports Series nr 64 zaleca używanie wielozakresowych mierników z automatyczną sygnalizacją saturacji.

Czas martwy w systemach wielokanałowych i akceleratorowych

W dużych systemach fizyki cząstek czas martwy nabiera wymiaru systemowego. Eksperyment w CERN (LHC, CMS, ATLAS) ma wiele warstw elektroniki: detektor → FE electronics → readout → trigger → DAQ. Każda warstwa ma własny czas martwy. Trigger (ang. wyzwolenie) decyduje, które zdarzenia są w ogóle zapisywane. Przy kolizjach LHC (40 MHz częstości zderzań) tylko ~1 000 zdarzeń na sekundę przechodzi przez trigger hardware i jest zapisywanych — to stosunek 1:40 000.

W dozymetrii reaktorowej czas martwy staje się krytyyczny przy pomiarach mocy neutronowej wysokim strumieniu. Detektory oparte na jonizacji gazu (komory jonizacyjne kompensowane) pracują w trybie prądowym, a nie impulsowym — w ten sposób omijają problem czasu martwego impulsowego. Komory rekombinacyjne (artykuł o spektrometrii neutronowej w serwisie) pracują podobnie.

Cyfrowe przetwarzanie sygnałów (DPP) jako nowoczesne rozwiązanie

Tradycyjne tory analogowe (detektor → wzmacniacz → kształtujący → MCA) zastępowane są przez DPP (Digital Pulse Processor). W systemie DPP sygnał z przedwzmacniacza ładunkowego (lub napięciowego) jest próbkowany przez szybki ADC (14–16 bitów, 125–500 MHz), a następnie przetwarzany numerycznie przez FPGA lub DSP.

Zalety DPP:

• Konfigurowalny shaping time bez wymiany układów analogowych
• Zintegrowany PUR (analiza kształtu impulsu, wykrywanie nakładania)
• Precyzyjna korekcja czasu martwego (live time rejestrowany z rozdzielczością ns)
• Możliwość pracy w trybach PHA (energii), PSD (dyskryminacji kształtu), MCS (zliczanie czasowe)
• Niższy szum niż analogowe łańcuchy przy tej samej rozdzielczości energetycznej

Producenci: CAEN (Włochy, systemy DT5782, N6780A), ORTEC (AUSA DPP-PHA), Mirion/Canberra (DSA-LX), Else Nuclear, XIA LLC (PIXIE-16). Każdy producent implementuje własny algorytm trapezoidalny, choć zasada jest ta sama.

Ograniczenia DPP: szybki ADC jest drogi, przetwarzanie cyfrowe przy bardzo krótkich czasach kształtowania (< 0,5 μs) wymaga dużego pasma próbkowania, a noise roll-off ADC może pogarszać stosunek sygnału do szumu przy krótkich czasach kształtowania.

Polska perspektywa: doświadczenia laboratoryjne i dydaktyka

Czas martwy i pile-up są tematami obowiązkowymi w laboratoriach promieniowania jonizującego polskich uczelni i instytutów:

IFJ PAN (Kraków): laboratorium HpGe prowadzi pomiary środowiskowe w geometrii bliskości, gdzie TCS i pile-up są istotne. Stosują systemy DPP firmy CAEN z analogowym wzmacniaczem ORTEC 672 jako porównaniem.

AGH (Kraków, Wydział Fizyki i Informatyki Stosowanej): ćwiczenie laboratoryjne z czasem martwym licznika GM (metoda dwóch źródeł ¹³⁷Cs), obowiązkowe dla studentów fizyki jądrowej i techniki jądrowej. Szczegółowe instrukcje bazują na podręczniku Dziunikowski-Kalita.

CLOR (Warszawa): rutynowe pomiary środowiskowe wymagają sprawdzenia live time/real time ratio dla każdego widma. Przy czasie martwym > 20% wymagane jest pomierzenie przy innej geometrii. Kodeks CLOR-T-2018 zaleca sprawdzenie PUR przez porównanie widma przy 10 cm i 50 cm od detektora.

WAT (Warszawa, Instytut Optoelektroniki): badania nad impulsowym promieniowaniem gamma (symulacja EMP i flash) wymagają detektorów o czasie martwym < 1 ns — stosuje się detektory diamentowe i diamondowe (CVD diamond detectors) właśnie dla skrajnie małego czasu martwego.

Politechnika Gdańska: katedra Fizyki Medycznej prowadzi badania nad OKD (Odrzucanie Kaskadowych Dodatkowych) dla kamer gamma w diagnostyce medycznej (SPECT, PET), gdzie pile-up od aktywności terapeutycznych jest kluczowym problemem klinicznym.

Statystyczny wpływ czasu martwego na niepewność pomiaru

Standardowe wzory statystyki Poissona zakładają, że zliczenia są niezależnymi zdarzeniami. Czas martwy narusza to założenie — jeżeli zliczenie Ni wyklucza możliwość zliczenia N{i+1} przez czas τ, to kolejne zliczenia nie są niezależne. Ma to konsekwencje dla oceny niepewności.

Dla modelu nieparalizowalnego:
Niepewność standardowa obserwowanej częstości m wynosi √m/T przy braku czasu martwego. Przy czasie martwym niepewność m jest nieco niższa (bo impulsy są „rozrzedzone" regularniej niż czysta statystyka Poissona), ale niepewność skorygowanego n jest większa ze względu na propagację niepewności przez wzór n = m/(1 - m·τ).

Pochodna: dn/dm = 1/(1-m·τ)². Dla m·τ = 0,1: dn/dm = 1/(0,9)² ≈ 1,23. Niepewność n jest ~23% większa niż przy braku korekcji. Dla m·τ = 0,3: dn/dm = 1/(0,7)² ≈ 2,04 — dwa razy większa.

W praktyce: przy pomiarach wymagających niepewności < 1% nie wolno dopuszczać m·τ > 0,05 bez specjalnej analizy propagacji niepewności.

Korelacja statystyczna przy metodzie dwóch źródeł: jeżeli używamy metody dwóch źródeł do wyznaczenia τ (mierzymy A, B i A+B), wszystkie cztery liczby (A, B, A+B, tło) są pomierzane z tym samym torem i tym samym τ, który wyznaczamy. Prowadzi to do korelacji między mierzonymi wielkościami. Analiza propagacji błędów w metodzie dwóch źródeł jest bardziej złożona niż proste dodawanie kwadratów niepewności. IAEA Technical Reports Series nr 323 podaje poprawne wyrażenia na niepewność τ z tej metody.

Wpływ na granicę wykrywalności (MDC): przy dużym czasie martwym wzrasta niepewność pomiaru tła, co podnosi granicę wykrywalności. Dlatego laboratoria środowiskowe często wolą mierzyć z mniejszą aktywnością preparatu (lepsza geometria, dłuższy czas) niż z dużą aktywnością i korekcją czasu martwego.

Czas martwy a pomiary impulsów nanosekundowych — powiązanie z szybką elektroniką

Czas martwy w elektronice jądrowej i czas martwy w szybkiej elektronice impulsowej to dwa aspekty tego samego zjawiska, ale napotykane w różnych kontekstach:

Pomiary detonacji i fali uderzeniowej (artykuł o tensometrach i czujnikach cisnienia): czujniki ciśnienia piezoelektryczne mogą generować impulsy nanosekundowej szerokości z amplitudami rzędu kilkunastu woltów. Przetwornik (charge amplifier lub IEPE) musi się zresetować po każdym impulsie. Dla pomiarów blisko ładunku (< 1 m), gdzie ciśnienie szczytowe może przekroczyć 1 GPa, ochrona wejścia i czas resetu wzmacniacza są krytyczne — i są formą czasu martwego.

Tory detektorów neutronowych w reaktorach: w pomiarach przejściowych mocy reaktora (np. testy SCRAM lub reakcji na zaburzenie) odpowiedź detektora musi być szybka. Komory szybkich neutronów (fission chambers) produkują impulsy < 10 ns, a ich czas martwy decyduje o wiarygodności pomiaru przy szybkich wzniosach mocy. WAT prowadzi badania nad takimi systemami w kontekście energetyki jądrowej.

Pomiary z szybkimi kamerami: w badaniach implosji lub wybuchów konwencjonalnych szybka kamera może rejestrować 10⁶–10⁹ klatek na sekundę. Układ sterowania (trigger, delay generator) musi działać z rozdzielczością nanosekundową. Jeżeli generacja jednego sygnału wyzwalającego „blokuje" generator na czas niedostosowany do kadencji kamery, pojawia się analogia do czasu martwego: okna ekspozycji są tracone. Artykuł o oscyloskopach i artykuł o szybkich kamerach omawiają sprzężony problem szybkości i wyłączności zasobu pomiarowego.

Numeryczne przykłady korekcji czasu martwego

Przykład A: pomiar aktywności ²⁴¹Am w glebie

Dane: Pomiar trwa T = 3600 s, zarejestrowano M = 120 000 zliczeń przy energii 59,5 keV. Live time = 3420 s. Real time = 3600 s. Czas martwy = (3600 - 3420)/3600 = 5%.

Częstość obserwowana: m = 120 000 / 3420 = 35,09 s⁻¹

Szacowana prawdziwa częstość (model nieparalizowalny, τ obliczony z live time): n ≈ m / (1 - czas_martwy) = 35,09 / 0,95 = 36,94 s⁻¹

Korekta jest mała (5%), wynik jest wiarygodny. Aktywność przeliczana z uwzględnieniem wydajności i geometrii.

Przykład B: spektrometria ⁶⁰Co przy przeciążeniu

Próbka aktywna. Pomiar T = 600 s, M = 3 500 000 zliczeń. Live time = 240 s. Czas martwy = (600 - 240)/600 = 60%. m = 3 500 000 / 240 = 14 583 s⁻¹.

Model nieparalizowalny: n = 14 583 / (1 - 0,6) = 36 458 s⁻¹. To bardzo duża korekta — przy błędzie pomiaru czasu martwego o 2% wartość bezwzględna n zmieni się o ~15%. Wniosek: wynik jest mało wiarygodny. Należy odsunąć próbkę lub skrócić czas pomiaru w nowej geometrii.

Przykład C: licznik GM przy bardzo dużej dawce

Miernik dozymetryczny GM wskazuje 0,5 mSv/h. Fizycznie źródło emituje ~10 mSv/h, ale licznik jest w głębokim nasyceniu modelu paralizowalnego. Odczyt 0,5 mSv/h jest poniżej prawdziwego — błąd w stronę mniejszej dawki. Standard PN-EN 60846-1 wymaga, aby mierniki dawki miały wskazaną wartość maksymalną i automatyczną flagę przy saturacji.

Diagnostyka widmowa i rozpoznawanie pile-upu

W praktycznym analizowaniu widm gamma istnieje kilka charakterystycznych sygnałów obecności pile-upu:

Podwyższone tło ciągłe powyżej piku: jeżeli pile-up tworzy przypadkowe zdarzenia sumacyjne, tło widma rośnie powyżej najwyższej energii rzeczywistego emitowanego fotonu. Dla ⁶⁰Co (max 1332 keV) nieoczekiwane zdarzenia przy 2000 keV sugerują pile-up 1173+662 (z ⁶⁰Co+Cs-137 w tle) albo 1332+1332 keV.

Poszerzenie piku: przy dużym pile-upie lewy ogon piku energetycznego (ang. low-energy tail) rośnie, bo część nakładających się impulsów jest zliczona z niedoborem amplitudy. Prawy ogon (powyżej piku) może rosnąć przez zdarzenia sumacyjne.

Asymetria piku: normalny pik HPGe jest prawie symetryczny gaussian. Przy pile-upie pik może być asymetryczny — bardziej rozciągnięty w prawą stronę (nadmiar energii z sumy impulsów).

Test diagnostyczny: pomiar tej samej próbki przy dwóch odległościach (np. 10 cm i 30 cm od detektora). Przy braku pile-upu pik powinien mieć tę samą energię i rozdzielczość (FWHM), a aktywność różnić się tylko przez geometrię (efektywność). Przy obecnym pile-upie pik przy 10 cm jest szerszy, tło jest wyższe i mogą pojawić się artefakty.

Test koincydencji: podłączenie generatora pulsów o znanych impulsach (np. 1000 s⁻¹ przy energii 500 keV wzorcowej) do wejścia sumacyjnego wzmacniacza. W widmie pojawi się pik 500 keV i pile-upowy pik „500+X keV" — jego intensywność obrazuje stopień pile-upu rzeczywistych zdarzeń detekcji.

Zasady projektowania systemu zliczającego z uwzględnieniem czasu martwego

Przy projektowaniu toru pomiarowego promieniowania warto stosować następujące zasady:

1. Określ spodziewaną maksymalną częstość zdarzeń: na podstawie szacowanej aktywności próbki, geometrii i wydajności detektora. Jeżeli spodziewasz się > 50 000 s⁻¹, zaplanuj korekcję czasu martwego już na etapie projektu.

2. Dobierz shaping time do wymagań: dla spektrometrii energetycznej minimalizuj szum kosztem czasu kształtowania tylko do granicy, przy której pile-up nie przekracza 5–10%. Dla liczenia bez spektrometrii (np. monitoring radioekologiczny) możesz wybrać krótki shaping time kosztem trochę gorszej rozdzielczości energetycznej.

3. Stosuj PUR jeśli dostępny: każdy nowoczesny detektor HPGe powinien pracować z aktywnym PUR. Koszt: zmniejszona live time. Zysk: wiarygodność widma.

4. Rejestruj live time i real time zawsze: nawet jeśli spodziewasz się małego czasu martwego. Zarejestrowany stosunek live/real w każdym pomiarze pozwala ex post ocenić wiarygodność wyników.

5. Sprawdzaj regularnie przez test stabilności wzmocnienia: przy dużej częstości impulsów wzmacniacz i zasilacz wysokiego napięcia mogą zmieniać punkt pracy. Kalibracyjne piki referencyjne (np. z wbudowanego generatora pulsów lub zewnętrznego źródła wzorcowego) pozwalają śledzić dryf wzmocnienia w czasie pomiaru.

6. Dla zmiennych aktywności — VPG lub LFC: jeżeli aktywność zmienia się w czasie (np. pomiary wczesnych opadów, monitoring reaktora), stosuj metodę wirtualnego generatora pulsów lub algorytmy LFC.

Błędy interpretacyjne i pułapki dydaktyczne

Błąd 1: więcej aktywności = lepszy wynik. Przy dużej aktywności pile-up i czas martwy mogą sprawić, że widmo jest gorszej jakości niż przy mniejszej aktywności mierzonej dłużej. Optymalna geometria to taka, gdzie zliczenia na sekundę pozwalają na czas martwy < 10–15%.

Błąd 2: live time to prawdziwy czas pomiaru. Live time (czas efektywny) jest podawany przez analizator jako czas pomiaru pomniejszony o czas martwy. Ale jeżeli analizator stosuje model uproszczony, live time może być niedokładny przy dużym pile-upie lub zmiennej aktywności.

Błąd 3: pile-up widać zawsze jako piki sumacyjne. Pile-up przy impulsach tego samego detektora może tworzyć ciągłe tło o podwyższonej amplitudzie, a nie wyraźne piki. Szczególnie przy ¹³⁷Cs (monoenergetyczna gamma 662 keV) nakładanie n impulsów rozciąga widmo do n × 662 keV, tworząc ciągłe rozszerzenie, nie ostrą linię.

Błąd 4: PUR usuwa problem pile-upu całkowicie. PUR odrzuca rozpoznane nakładki, ale nie te, których oba impulsy przyszły zbyt blisko siebie (poniżej rozdzielczości czasowej systemu szybkiego). Przy ekstremalnie krótkich czasach między impulsami (<< czasu narastania impulsu) PUR nie może ich rozróżnić.

Błąd 5: czas martwy ADC dominuje. W nowoczesnych DPP czas martwy ADC (konwersja) jest pomijalny (< 100 ns). Dominuje czas kształtowania i logiki zapisu. Przy starych analogowych MCA czas ADC był istotny — do 10 μs dla 12-bitowego sukcesywnego aproksymowania.

Błąd 6: wszystkie systemy DPP mają identyczny czas martwy. Czas martwy w DPP zależy od ustawień kształtowania, od modelu sprzętu i od konfiguracji oprogramowania. Systemy DPP różnych producentów mogą mieć czas martwy różniący się o rząd wielkości przy tych samych ustawieniach czasu kształtowania. Warto zmierzyć czas martwy własnego systemu metodą dwóch źródeł i porównać z deklarowanymi danymi producenta.

Błąd 7: czas martwy jest stały dla danego toru. Przy zmiennej aktywności (np. krótkotrwałe tarczki aktywacyjne, monitoring opadu) czas martwy może zmieniać się w czasie rzeczywistego pomiaru. Analizator zapisuje wtedy sumaryczny live time, który jest uśredniony — i może nie opisywać prawidłowo rozkładu zdarzeń w czasie. W takich przypadkach niezbędne są systemy rejestracji zdarzenie-po-zdarzeniu (ang. list mode acquisition), gdzie każde zdarzenie ma znacznik czasu. List mode pozwala na rekonstrukcję widma w dowolnych oknach czasowych po zakończeniu pomiaru.

Błąd 8: sumowanie kaskadowe i pile-up mają tę samą zależność od odległości. Pile-up maleje z odległością (liniowo przez aktywność źródła × efektywność). Sumowanie kaskadowe maleje z odległością szybciej — proporcjonalnie do kwadratu efektywności całkowitej detektora. Dlatego przy małych geometriach TCS dominuje nad pile-upem nawet przy nieduże aktywności, a przy dużych geometriach pile-up może dominować przy wysokich aktywnościach. Zrozumienie tej różnicy jest kluczowe przy kalibracji wydajnościowej detektora HPGe.

Znajomość tych pułapek odróżnia użytkownika, który tylko odczytuje wynik oprogramowania, od użytkownika, który rozumie granice fizyczne systemu pomiarowego.

Dodatkowe materiały multimedialne

Powiązane kalkulatory i narzędzia

• Wizualizacja: Tor impulsowy — Impulsy detektora przesuwają się po osi czasu, a próg i czas martwy decydują, które zdarzenia zostaną zapisane.
• Wizualizacja: Koincydencje — Dwa tory impulsowe pokazują, jak szerokość okna czasowego zmienia udział koincydencji prawdziwych i przypadkowych.
• Kalkulator: Czas martwy — Porównanie modelu nieparalizowalnego i paralizowalnego oraz klasyczna metoda dwóch źródeł.
• Kalkulator: Koincydencje — Korekta przypadkowych zbiegów w dwukanałowym układzie koincydencyjnym.

Powiązane artykuły

• Tor impulsowy NIM w praktyce dydaktycznej
• Od analogowego NIM do cyfrowego DPP
• Koincydencje, veto i bramkowanie
• Spektrometria gamma w praktyce

Ćwiczenia praktyczne

Pierwsze ćwiczenie powinno być oparte wyłącznie na danych syntetycznych. Dane wejściowe:

Pomiar	Czas pomiaru	Liczba zliczeń
Tło	60 s	1 200
Źródło A + tło	60 s	24 600
Źródło B + tło	60 s	18 900
Źródła A+B + tło	60 s	40 200

Zadanie:

1. obliczyć częstości zliczeń brutto,
2. odjąć tło,
3. porównać sumę częstości netto A+B z pomiarem wspólnym,
4. wyjaśnić, dlaczego pomiar wspólny jest niższy niż suma,
5. oszacować, czy pomiar wymaga korekty czasu martwego.

Drugie ćwiczenie powinno dotyczyć interpretacji widma. Student dostaje dwa syntetyczne widma tej samej próbki: jedno zmierzone przy małej częstości, drugie przy przeciążonym torze. Należy wskazać, w którym widmie widać objawy pile-upu, dlaczego wzrost liczby zliczeń nie zawsze poprawia jakość wyniku i jak należałoby zmienić geometrię albo czas pomiaru.

Szybka lista kontrolna przed pomiarem

Przy każdym pomiarze, w którym liczy się dokładność aktywności lub widmo energetyczne, warto sprawdzić kilka punktów związanych z czasem martwym i pile-upem:

1. Sprawdź czas martwy: przy rate wyższym niż 5000 zliczeń/s (cps) czas martwy powyżej 1% zaczyna być istotny. Jeśli MCA pokazuje >10% dead time — aktywność próbki lub geometria wymagają zmiany.

2. Zmierz rate bez okna energetycznego: jeśli MCA pozwala podejrzeć całkowite zliczenia (total count rate przed ADC), porównaj je z ratą w fotopiku — różnica to głównie Compton i tło.

3. Sprawdź pile-up rejector: jeśli moduł ma PUR, upewnij się, że jest włączony. Bez PUR widmo przy dużym rate ma artefakty sumacyjne.

4. Zmierz live time, nie real time: przy obliczaniu aktywności z formuły A = N/(ε·I_γ·t) użyj t = live time (t_live), nie t_real. Błąd jest równy frakcji czasu martwego.

5. Weryfikacja przez kalibrację rate: zmierz znane źródło przy różnych odległościach (różnych rate) — jeśli aktywność obliczona z pomiaru jest stała niezależnie od rate, korekcja czasu martwego działa poprawnie.

Przejdź do ćwiczenia interaktywnego