Ders 4 / 20Orta~40 dk kaldıGlenn F. Knoll (4th ed.)

Radyasyon Dedektörlerinin Genel Özellikleri

Bu bölüm, kullandığı fiziksel mekanizmadan bağımsız olarak her radyasyon dedektörüne uygulanan evrensel çerçeveyi geliştirir. Tek bir radyasyon kuantumunun enerji bıraktığı, Q = E·q/W kadar yük serbest bıraktığı ve zaman integrali enerji bilgisini taşıyan bir akım darbesi ürettiği basitleştirilmiş dedektör modeliyle başlar. Bu modelden üç çalışma modu türetilir: darbe modu (her olay tek tek kaydedilir, spektroskopiyi mümkün kılar), akım modu (doza veya akı hızına orantılı zaman-ortalamalı akım) ve ortalama-kare-gerilim (Campbelling) modu (olayları yüklerinin karesiyle ağırlıklandırır ve gama arka planında nötron dedeksiyonu gibi karışık radyasyon alanları için paha biçilmezdir). Bölüm daha sonra dedektör fiziğinin merkezi tanı aracı olan diferansiyel darbe yüksekliği spektrumunu tanımlar ve bir diskriminatör eşiği veya çalışma gerilimi taranırken sayım eğrileri ile platoların nasıl ortaya çıktığını gösterir. İki temel başarım ölçütünü nicelendirir: Poisson istatistiği, Fano faktörü ve gürültü kaynaklarının karesel toplamından türetilen enerji çözünürlüğü; ve mutlak, içsel, toplam ve fotopik kategorilerine ayrılan dedeksiyon verimi. Son olarak ölü zamanı felç olabilen (paralyzable) ve felç olmayan (nonparalyzable) modellerle ele alır, sayım hızı düzeltme formüllerini türetir ve iki-kaynak ile bozunan-kaynak yöntemlerini sunar. Bu bölüme hakim olmak, öğrenciye kitabın ilerleyen kısımlarındaki herhangi bir dedektör türünden gelen verileri yorumlamak ve iyi anlaşılmış belirsizliklerle ölçümler tasarlamak için gereken kavramsal ve nicel araçları kazandırır.

Basitleştirilmiş Dedektör Modeli

Tüm dedektörler ortak bir soyutlamayı paylaşır: etkileşen bir radyasyon kuantumu aktif hacme E enerjisi bırakır ve bu enerji ölçülebilir bir Q elektrik yüküne dönüştürülür. Bilgi taşıyıcılarının (elektron-iyon çiftleri, elektron-boşluk çiftleri veya fotoelektronlar) sayısı N = E/W olup, W taşıyıcı başına harcanan ortalama enerjidir; toplanan yük ise Q = N·q'dur. Yük toplanması ani değildir; dedektör, toplama süresi boyunca zaman integrali Q'ya eşit olan geçici bir i(t) akım darbesi üretir. i(t) biçimi taşıyıcı sürüklenme hızlarına ve geometriye bağlıdır, ancak enerji ölçümleri için yalnızca toplam integralli yük önemlidir. Bu model gaz, sintilasyon ve yarıiletken dedektörlerin ayrıntılı fiziğini soyutlayarak aynı istatistiksel ve sinyal işleme kavramlarının tekdüze uygulanmasını sağlar. Ayrıca Q'nun büyüklüğünün, dolayısıyla tipik sinyalin, dedektör teknolojileri arasında büyüklük mertebeleriyle değiştiğini açıkça ortaya koyar; bu da çalışma modu ve elektronik seçimini belirler.

information carriersW-valuecharge collectioncurrent pulse i(t)active volumeenergy deposition

Dedektör Çalışma Modları

Basitleştirilmiş modelden gelen aynı akım darbesi dizisi üç farklı şekilde işlenebilir. Darbe modunda her etkileşim tek tek kaydedilir; olay başına integralli yük korunur, böylece enerji spektrumu yeniden oluşturulabilir. Bu, tüm radyasyon spektroskopisinin temelidir ancak olay hızları darbeler örtüşecek kadar yüksek olduğunda başarısız olur. Akım modunda T tepki süreli bir pikoampermetre akımı birçok olay üzerinden ortalar ve etkileşim hızı ile olay başına ortalama yükün çarpımına orantılı I0 = r·Q kararlı sinyalini verir; göreli istatistiksel dalgalanma 1/sqrt(2·r·T) ile azalır. Akım modu, bireysel olayların ayırt edilmesi gerekmeyen dozimetri ve yüksek-akı izlemede kullanılır. Ortalama-kare-gerilim (MSV veya Campbelling) modu yalnızca akımın dalgalanan bileşenini geçirir ve ortalama karesini oluşturur; r·Q^2'ye orantılı bir sinyal üretir. Ağırlıklandırma yükün karesiyle gittiğinden MSV modu büyük yük bırakan olayları güçlü biçimde kayırır ve dedektörlerin darbe modunun imkansız olduğu yoğun karışık alanlarda nötronları gama ışınlarından ayırt etmesini sağlar.

pulse modecurrent modemean square voltage (MSV) modeCampbelling moderesponse time Tfluctuating current component

Darbe Yüksekliği Spektrumları

Bir dedektör darbe modunda çalıştığında, ardışık darbelerin genlikleri olaydan olaya, bırakılan enerjideki farklar ve istatistiksel dalgalanmalar nedeniyle değişir. Bu dağılımı temsil etmenin en bilgilendirici yolu, darbe yüksekliği H'ye karşı çizilen ve pratikte darbeleri genlik kutularına ayıran bir çok-kanallı analizör (MCA) ile kaydedilen diferansiyel darbe yüksekliği spektrumu dN/dH'dir. Genliği H1 ile H2 arasında olan darbelerin sayısı, bu sınırlar arasındaki diferansiyel spektrumun altındaki alana eşittir ve kaydedilen toplam olay sayısı toplam alandır. Spektrumdaki tepeler tek-enerjili enerji bırakımlarına karşılık gelir; süreklilikler kısmi enerji bırakımından veya sürekli kaynak spektrumlarından doğar. İntegral darbe yüksekliği spektrumu, belirli bir eşiği aşan darbe sayısı, tamamlayıcı temsildir; diferansiyel spektrum integral spektrumun negatif eğimidir. Bu spektrumları okumak ve yorumlamak radyasyon spektroskopisindeki merkezi beceri ve çözünürlük, verim ile dedektör tepki fonksiyonlarını anlamanın temelidir.

differential pulse height spectrumintegral pulse height spectrummultichannel analyzer (MCA)channel binningfull-energy peakcontinuum

Sayım Eğrileri ve Platolar

Birçok basit sayım uygulamasında tek bir integral diskriminatör, ayarlanabilir bir eşiğin üzerindeki tüm darbeleri kaydeder (veya dedektör besleme gerilimi değiştirilir). Kaydedilen sayım hızının bu kontrol parametresine karşı çizilmesi bir sayım eğrisi üretir. Darbeleri kusursuzca tekdüze olan ideal tek-enerjili bir kaynak keskin bir basamak verirdi; sonlu genişlikli ve elektronik gürültülü gerçek spektrumlar bu basamağı eğimli bir bölgeye yayar. İyi tasarlanmış bir dedektör, sayım hızının diskriminatör ayarından veya uygulanan gerilimden neredeyse bağımsız olduğu bir plato sergiler; çünkü gerçek darbelerin büyük kısmı gürültünün epey üzerinde ama herhangi bir sahte çoğalmanın altındadır. Platonun varlığı, uzunluğu ve eğimi pratik kalite göstergeleridir: uzun, düz bir plato ölçümün eşik veya yüksek gerilimdeki küçük kaymalara duyarsız olduğu, kararlı ve tekrarlanabilir sonuçlar verdiği anlamına gelir. Çalışma noktası platonun merkezine yakın seçilir. Geiger-Mueller tüpleri ve orantılı sayaçlar rutin olarak plato eğrileriyle karakterize edilir ve plato eğimi (her 100 V için yüzde) standart bir spesifikasyondur.

counting curveintegral discriminatorplateauplateau slopeoperating pointthreshold setting

Enerji Çözünürlüğü

Enerji çözünürlüğü, bir dedektörün birbirine yakın iki enerjiyi ayırt etme yeteneğini ölçer ve R = FWHM/H0 olarak tanımlanır; bir tam-enerji tepesinin yarı yükseklikteki tam genişliğinin tepe merkezine bölümüdür; daha küçük R daha iyidir. Gauss tepesinde FWHM = 2.35·sigma'dır. Temel istatistiksel sınır, yük taşıyıcı sayısı N = E/W'nin dalgalanmasından doğar. Taşıyıcı üretimi bir Poisson süreci olsaydı sigma = sqrt(N) ve R_Poisson = 2.35/sqrt(N) olurdu; dolayısıyla R < %1 elde etmek N > 55.000 taşıyıcı gerektirir. Birçok gerçek dedektör bundan daha iyidir çünkü taşıyıcı üretimi bağımsız değil ilişkilidir; bu sapma Fano faktörü F = (N'deki gözlenen varyans)/N ile nicelenir ve R_istatistiksel = 2.35·sqrt(F/N) verir. F, yarıiletkenler ve orantılı sayaçlar için birden epey küçük, sintilatörler için ise bire yakındır. İstatistiğin ötesinde elektronik gürültü, yük toplama değişimleri ve kayma tepeyi genişletir; her kaynak bağımsız ve simetrikse toplam tepki, FWHM'leri karesel olarak birleşen bir Gauss'a yaklaşır: (FWHM_toplam)^2 = (FWHM_ist)^2 + (FWHM_gürültü)^2 + (FWHM_kayma)^2 + ... Bu ayrıştırma deneycinin baskın genişleme mekanizmasını izole edip ona müdahale etmesini sağlar.

FWHMenergy resolution RFano factorPoisson limitGaussian response functionquadrature noise additionelectronic noise

Dedeksiyon Verimi

Dedeksiyon verimi, kaydedilen darbe sayısını radyasyon kuantumlarının sayısına bağlar. Yüklü parçacıklar girişte hemen iyonlaştığından dedektörler çoğu zaman aktif hacme giren her parçacığı görebilir (yaklaşık %100 sayım verimi). Yüksüz radyasyonlar (gama ışınları, nötronlar) önce etkileşmelidir ve büyük mesafeler kat edebilir, bu yüzden verim genellikle %100'ün epey altındadır ve sayımları kaynak aktivitesine dönüştürmek için bilinmelidir. İki tamamlayıcı tanım kullanılır: kaynak-dedektör geometrisine bağlı olan mutlak verim eps_abs = (kaydedilen darbeler)/(kaynak tarafından yayılan kuantumlar); ve esas olarak dedektör malzemesi, kalınlığı ve radyasyon enerjisine bağlı olan içsel verim eps_int = (kaydedilen darbeler)/(dedektöre gelen kuantumlar). İzotropik bir nokta kaynak için bunlar geometrik katı-açı faktörüyle ilişkilidir: eps_abs = eps_int · (Omega/4*pi). Verimler ayrıca neyin sayıldığına göre sınıflandırılır: toplam verim darbe üreten herhangi bir etkileşimi sayar, tepe (veya fotopik) verimi ise yalnızca tam enerjiyi bırakan ve tam-enerji tepesini dolduran olayları sayar. Tepe-toplam oranı etkileşimlerin tepeye düşen kesrini yakalar. Bu ayrımlar nicel spektroskopi ve kaynak-aktivitesi belirlemede temeldir.

absolute efficiencyintrinsic efficiencytotal efficiencypeak (photopeak) efficiencysolid anglepeak-to-total ratiocounting efficiency

Ölü Zaman

Her sayım sistemi, iki olayın ayrı kaydedilmesi için aralarında minimum bir tau süresi gerektirir; dedektör veya elektronik tarafından belirlenen bu ölü zaman, yüksek hızlarda sayımların kaybedilmesine neden olur ve düzeltilmesi gerekir. İki idealleştirilmiş model gerçek davranışı sınırlar. Felç olmayan (nonparalyzable) modelde her kaydedilen olayı, sonraki olayların basitçe kaybolduğu ve başka etkisinin olmadığı sabit bir tau ölü periyodu izler; bu m = n/(1 + n·tau) ve açık düzeltme n = m/(1 - m·tau) verir; gözlenen hız 1/tau'da doyar. Felç olabilen (paralyzable) modelde her gerçek olay, ölü periyot sırasında gerçekleşen bile olsa, ölü periyodu yeniden başlatır; böylece m = n·exp(-n·tau) olur; bu kapalı biçimde tersine çevrilemez, iteratif çözülmelidir ve gözlenen hız bir tepe yapıp düşer, dolayısıyla belirli bir m iki farklı gerçek hıza n karşılık gelebilir. Her iki model de düşük kayıpta birinci mertebeden uyuşur: m ~ n(1 - n·tau). İki deneysel yöntem tau'yu belirler: doğrusal olmamayı kullanan, iki kaynağı ayrı ayrı ve birlikte sayan iki-kaynak yöntemi ve gözlenen hızın bilinen üstel bozunumdan sapmasını uyduran bozunan-kaynak yöntemi. %30-40 üzeri ölü zaman kayıpları, düzeltmeler model ve parametre belirsizliklerine çok duyarlı hale geldiğinden kaçınılmalıdır.

dead time tauparalyzable modelnonparalyzable modelcount-rate correctiontwo-source methoddecaying-source methodrate saturation

Olay Başına Üretilen Yük

E enerjisi bırakan tek bir radyasyon etkileşiminin serbest bıraktığı toplam yük; burada N = E/W bilgi taşıyıcısı oluşur. Basitleştirilmiş dedektör modelinin temeli; bırakılan enerjiyi tüm dedektör türlerinde ölçülebilir sinyale bağlar.

Temel Denklemler (23)

Bu denklemi çalıştır: Ölü Zaman Modelleri: Felç Olabilen vs Felç Olmayan →

Olay Başına Üretilen Yük

(4.1)

Q = N q = \frac{E}{W} q

sembolanlambirim

Q

olay başına toplanan toplam yükC

N

bilgi taşıyıcılarının sayısı–

E

aktif hacimde bırakılan enerjiJ

W

bir taşıyıcı oluşturmak için gereken ortalama enerjiJ

q

temel yük, 1.602e-19 CC

Akım Darbesinin İntegrali Olarak Yük

(4.2)

Q = \int_{0}^{t_c} i(t)\, dt

Toplam yük, geçici akım darbesinin yük toplama süresi t_c boyunca zaman integraline eşittir. Enerji bilgisinin darbe biçiminden bağımsız olarak integralli yükte korunduğunu gösterir; bu da yük-duyarlı önyükselteçleri gerekçelendirir.

sembolanlambirim

i(t)

anlık dedektör akımıA

t_c

yük toplama süresis

Q

toplanan toplam yükC

Akım Modu Ortalama Sinyali

(4.3)

I_0 = r\,Q = r\,\frac{E}{W} q

Akım modunda kararlı durum akımı, olay hızı ile olay başına ortalama yükün çarpımına eşittir. Bireysel olayların ayırt edilmediği dozimetri ve yüksek-akı izlemede kullanılır.

sembolanlambirim

I_0

ortalama (DC) akımA

r

olay (etkileşim) hızı1/s

Q

olay başına ortalama yükC

Akım Modu Zaman Ortalaması

(4.4)

I(t) = \frac{1}{T}\int_{t-T}^{t} i(t')\, dt'

Kaydedilen akım, anlık dedektör akımının ölçüm aletinin tepki süresi T üzerinden hareketli zaman ortalamasıdır. Uzun T dalgalanmaları bastırır ama radyasyon şiddetindeki değişimlere tepkiyi yavaşlatır.

sembolanlambirim

I(t)

kaydedilen ortalama akımA

T

ölçüm aleti tepki (integrasyon) süresis

i(t')

anlık akımA

Akım Modunda Kesirli İstatistiksel Dalgalanma

(4.5)

\frac{\sigma_I}{I_0} = \frac{1}{\sqrt{2\,r\,T}}

Eşit darbe yükü varsayımıyla rastgele darbe varış zamanlarından kaynaklanan akım-modu sinyalinin göreli standart sapması. Etkin olay sayısı 2rT'nin kareköküyle iyileşir; integrasyon süresi seçimini yönlendirir.

sembolanlambirim

\sigma_I

ölçülen akımın standart sapmasıA

I_0

ortalama akımA

r

olay hızı1/s

T

tepki süresis

Ortalama Kare Gerilim (Campbelling) Sinyali

(4.6)

\overline{\sigma_i^2} \propto r\,Q^2

MSV modunda dalgalanan akım bileşeninin ortalama karesi, hız ile olay başına yükün karesinin çarpımına orantılıdır. Q^2 ağırlıklandırması nedeniyle MSV modu büyük-yüklü olayları (örn. nötronlar) küçük-yüklü olaylara (örn. gama ışınları) göre güçlendirir ve karışık alanlarda ayrımı mümkün kılar.

sembolanlambirim

\overline{\sigma_i^2}

dalgalanan akımın ortalama karesiA^2

r

olay hızı1/s

Q

olay başına yükC

Enerji Çözünürlüğü Tanımı

(4.7)

R = \frac{\mathrm{FWHM}}{H_0}

Enerji çözünürlüğü, bir tam-enerji tepesinin yarı yükseklikteki tam genişliğinin tepe merkezine bölümüdür; boyutsuz bir kesir (çoğunlukla yüzde). Daha küçük R, komşu enerjileri ayırma yeteneğinin daha iyi olduğu anlamına gelir.

sembolanlambirim

R

enerji çözünürlüğü–

\mathrm{FWHM}

tepenin yarı yükseklikteki tam genişliği–

H_0

tepenin merkezi (ortalama darbe yüksekliği)–

Gauss FWHM

(4.8)

\mathrm{FWHM} = 2.355\,\sigma

Bir Gauss tepesinin yarı yükseklikteki tam genişliği ile standart sapması arasındaki ilişki. Birçok bağımsız dalgalanmanın baskın olduğu dedektör tepki fonksiyonları Gauss biçimine yaklaşır; bu da bu dönüşümü her yerde geçerli kılar.

sembolanlambirim

\mathrm{FWHM}

yarı yükseklikteki tam genişlik–

\sigma

Gauss standart sapması–

Poisson İstatistiksel Çözünürlük Sınırı

(4.9)

R_{\mathrm{Poisson}} = \frac{2.355\,K\sqrt{N}}{K N} = \frac{2.355}{\sqrt{N}}

Yük taşıyıcı üretimi saf Poisson istatistiğine uyuyorsa elde edilebilecek en iyi çözünürlük; yalnızca taşıyıcı sayısı N'ye bağlıdır. R < %1 gereksinimi N > 55.000 anlamına gelir; bu da yüksek-taşıyıcı-sayılı yarıiletken dedektörleri motive eder.

sembolanlambirim

R_{\mathrm{Poisson}}

Poisson-sınırlı çözünürlük–

N

olay başına yük taşıyıcı sayısı–

K

doğrusal orantı sabiti (H_0 = KN)–

Fano Faktörü

(4.10)

F = \frac{\sigma_N^2}{N}

Taşıyıcı sayısındaki gözlenen varyansın Poisson-öngörülü varyans N'ye oranı. F < 1 ilişkili taşıyıcı üretimini gösterir (yarıiletkenler, orantılı sayaçlar); sintilatörler için F ~ 1. Poisson istatistiğinden sapmayı nicelendirir.

sembolanlambirim

F

Fano faktörü–

\sigma_N^2

taşıyıcı sayısındaki gözlenen varyans–

N

ortalama taşıyıcı sayısı–

Fano Faktörlü İstatistiksel Çözünürlük

(4.11)

R_{\mathrm{stat}} = 2.355\,\sqrt{\frac{F}{N}}

Fano faktörünü içeren çözünürlük üzerindeki istatistiksel sınır. F < 1 için elde edilebilir çözünürlük naif Poisson tahmininden daha iyidir; bu da gerçek yarıiletken dedektörlerin neden sqrt(N) ölçeklenmesinden daha iyi olduğunu açıklar.

sembolanlambirim

R_{\mathrm{stat}}

istatistiksel çözünürlük sınırı–

F

Fano faktörü–

N

yük taşıyıcı sayısı–

FWHM Katkılarının Karesel Birleşimi

(4.12)

(\mathrm{FWHM})_{\mathrm{total}}^2 = (\mathrm{FWHM})_{\mathrm{stat}}^2 + (\mathrm{FWHM})_{\mathrm{noise}}^2 + (\mathrm{FWHM})_{\mathrm{drift}}^2 + \cdots

Birden çok bağımsız, simetrik genişleme kaynağı varken genel tepki Gauss'a yaklaşır ve toplam FWHM karesel toplanır. Deneycilerin her katkıyı izole etmesini sağlar (örn. bir darbe üreteci yalnızca elektronik gürültüyü ölçer).

sembolanlambirim

(\mathrm{FWHM})_{\mathrm{total}}

genel tepe genişliği–

(\mathrm{FWHM})_{\mathrm{stat}}

istatistiksel (taşıyıcı) katkı–

(\mathrm{FWHM})_{\mathrm{noise}}

elektronik gürültü katkısı–

(\mathrm{FWHM})_{\mathrm{drift}}

kayma/kararsızlık katkısı–

Gauss Tepki Fonksiyonu

(4.13)

G(H) = \frac{A}{\sigma\sqrt{2\pi}}\,\exp\!\left[-\frac{(H-H_0)^2}{2\sigma^2}\right]

Bir darbe yüksekliği spektrumundaki tam-enerji tepesinin biçimi için standart model; alan A, merkez H_0 ve genişlik sigma ile parametrelenir. Spektroskopide tepe uydurma, çözünürlük çıkarımı ve net-alan belirlemede kullanılır.

sembolanlambirim

G(H)

H darbe yüksekliğindeki diferansiyel sayım–

A

tepe alanı (toplam sayım)–

H_0

merkez darbe yüksekliği–

\sigma

standart sapma (genişlik)–

Mutlak Verim

(4.14)

\varepsilon_{\mathrm{abs}} = \frac{\text{number of pulses recorded}}{\text{number of quanta emitted by source}}

Kaynak tarafından yayılan tüm radyasyon kuantumlarından kaydedilen darbe üretenlerin kesri. Hem dedektör özelliklerine hem de kaynak-dedektör geometrisine bağlıdır; sayımları kaynak aktivitesine dönüştürmek için doğrudan kullanılır.

sembolanlambirim

\varepsilon_{\mathrm{abs}}

mutlak verim–

İçsel Verim

(4.15)

\varepsilon_{\mathrm{int}} = \frac{\text{number of pulses recorded}}{\text{number of quanta incident on detector}}

Dedektöre fiilen çarpan kuantumlardan kaydedilen darbe üretenlerin kesri. Esas olarak dedektör malzemesine, kalınlığına ve radyasyon enerjisine bağlıdır, geometriye değil; bir dedektör özelliği olarak tablolanır.

sembolanlambirim

\varepsilon_{\mathrm{int}}

içsel verim–

Mutlak-İçsel Verim İlişkisi

(4.16)

\varepsilon_{\mathrm{abs}} = \varepsilon_{\mathrm{int}} \cdot \frac{\Omega}{4\pi}

İzotropik bir nokta kaynak için mutlak verim, dedektörün kapsadığı kesirli katı açı Omega/4*pi ile ölçeklenmiş içsel verimdir. Verim hesaplarında geometrik ve fizik faktörlerini ayırır.

sembolanlambirim

\varepsilon_{\mathrm{abs}}

mutlak verim–

\varepsilon_{\mathrm{int}}

içsel verim–

\Omega

dedektörün kapsadığı katı açısr

Eşeksenli Diskin Katı Açısı

(4.17)

\Omega = 2\pi\left(1 - \frac{d}{\sqrt{d^2 + a^2}}\right)

d mesafesindeki, a yarıçaplı dairesel bir dedektörün ekseni üzerindeki bir nokta kaynakta kapsanan katı açı. Geometrik verim faktörünü verir; d >> a için Omega ~ pi*a^2/d^2 (ters-kare yasası) haline gelir.

sembolanlambirim

\Omega

katı açısr

a

dedektör yüzü yarıçapım

d

kaynak-dedektör mesafesim

Felç Olmayan Gözlenen Hız

(4.18)

m = \frac{n}{1 + n\tau}

Ölü zamanı tau olan felç olmayan bir sistem için gözlenen hız m ile gerçek hız n arasındaki ilişki. Gözlenen hız n büyüdükçe 1/tau'da doyar. Eşdeğer ifade: n - m = n*m*tau.

sembolanlambirim

m

gözlenen (kaydedilen) sayım hızı1/s

n

gerçek etkileşim hızı1/s

\tau

sistem ölü zamanıs

Felç Olmayan Gerçek-Hız Düzeltmesi

(4.19)

n = \frac{m}{1 - m\tau}

Felç olmayan bir sistem için ölçülen hızdan gerçek hızı geri kazanmaya yönelik kapalı formdaki düzeltme. Kayıplar orta düzeydeyken sayım verilerine uygulanan standart ölü-zaman düzeltmesi.

sembolanlambirim

n

gerçek etkileşim hızı1/s

m

gözlenen sayım hızı1/s

\tau

ölü zamans

Felç Olabilen Gözlenen Hız

(4.20)

m = n\,e^{-n\tau}

Her olayın ölü periyodu uzattığı felç olabilen bir sistem için gözlenen hız. Analitik olarak tersine çevrilemez (iteratif çözülmelidir) ve n = 1/tau'da bir maksimumdan geçer; böylece bir gözlenen hız iki olası gerçek hıza karşılık gelir.

sembolanlambirim

m

gözlenen sayım hızı1/s

n

gerçek etkileşim hızı1/s

\tau

ölü zamans

Her İki Modelin Düşük-Kayıp Sınırı

(4.21)

m \approx n\,(1 - n\tau), \quad n\tau \ll 1

Her iki ölü-zaman modeli de kayıplar küçükken birinci mertebeden uyuşur. Felç olabilen davranışı felç olmayandan ayırt etmenin yüksek sayım hızlarında ölçüm gerektirdiğini doğrular.

sembolanlambirim

m

gözlenen hız1/s

n

gerçek hız1/s

\tau

ölü zamans

İki-Kaynak Ölü Zaman Denklemi

(4.22)

\frac{m_{12}}{1 - m_{12}\tau} + \frac{m_b}{1 - m_b\tau} = \frac{m_1}{1 - m_1\tau} + \frac{m_2}{1 - m_2\tau}

İki-kaynak yönteminin temeli (felç olmayan model): düzeltilmiş gerçek hızlar toplandığından, birleşik gözlenen hız bireysel hızların toplamından küçüktür ve tau bu farktan çözülebilir. En iyi sonuçlar kesirli ölü zaman yaklaşık %20'yi aştığında elde edilir.

sembolanlambirim

m_1

kaynak 1'in gözlenen hızı (arka plan dahil)1/s

m_2

kaynak 2'nin gözlenen hızı1/s

m_{12}

iki kaynak birlikte gözlenen hız1/s

m_b

gözlenen arka plan hızı1/s

\tau

ölü zamans

Bozunan-Kaynak Gerçek Hızı

(4.23)

n = n_0\,e^{-\lambda t}

Bozunan-kaynak yönteminde kullanılan gerçek hızın bilinen üstel bozunumu. Ölçülen m'nin bu eğriden sapması bir ölü-zaman modeline yerleştirilerek grafiksel veya uydurmalı analizle tau verir.

sembolanlambirim

n

t anındaki gerçek hız1/s

n_0

ölçüm başındaki gerçek hız1/s

\lambda

izotopun bozunum sabiti1/s

t

geçen süres

Simülasyon tasarımları

SİMÜLASYON TASARIMIYÜKSEK SADAKAT

Ölü Zaman Modelleri: Felç Olabilen vs Felç Olmayan

Öğrenci gerçek etkileşim hızı n'yi tarar ve gözlenen hız m'nin her iki model eğrisini çizmesini izler. Felç olmayan eğri 1/tau'da doyarken, felç olabilen eğri n=1/tau'da bir tepeye yükselip geri katlanır ve iki-gerçek-hız belirsizliğini canlı biçimde gösterir. Sürüklenebilir bir gözlenen-hız çizgisi, karşılık gelebileceği bir veya iki gerçek hızı gösterir; canlı bir olay-treni animasyonu rastgele olayları ölü periyotlarıyla çizerek öğrencinin tek tek sayımların kaybedildiğini görmesini sağlar.

GİRDİLER

Dead time tau (0.1 - 100 us, log slider) · True rate n sweep range (1e2 - 1e7 1/s, log) · Model selector (paralyzable / nonparalyzable / both) · Live simulation true rate n (1e3 - 5e6 1/s) for the event-train animation

ÇIKTILAR

m-vs-n curve plot for both models with 1/tau and 1/(e*tau) reference lines · Fractional-loss (1 - m/n) readout vs n · Animated event train showing live periods, dead periods, recorded vs lost counts · Empirical m from the Monte Carlo overlaid on the analytic curves

YÖNTEMAnalytic curves m=n/(1+n*tau) and m=n*exp(-n*tau) evaluated on a log-spaced n grid; Monte Carlo event generation with exponentially distributed inter-arrival times (inverse-CDF sampling -ln(U)/n) processed through both a nonparalyzable and a paralyzable state machine to count recorded events and validate the analytic curves.

KOD KANCASIÖğrenci felç olabilen ve felç olmayan durum makinelerini uygular (gelen bir olayın kaydedilip kaydedilmeyeceği ve bir sonraki canlı zamanı nasıl güncellediği kuralı). Kodları Monte Carlo sayımını yönlendirir; doğru bir uygulama simüle edilmiş noktaların analitik eğrilerin tam üzerine düşmesini sağlar, hatalar ise özellikle geri-katlanma bölgesinde gözle görülür şekilde sapar.

Konsola Gir: dead-time →

SİMÜLASYON TASARIMIYÜKSEK SADAKAT

Enerji Çözünürlüğü ve Fano Faktörü

Öğrenci bir dedektör enerjisi ve dedektör türü (gaz, sintilatör, Si, Ge) seçer; bu W-değerini ve Fano faktörünü sabitler. Simülatör N=E/W'yi hesaplar, taşıyıcı-sayısı varyansı F*N olan bir Monte Carlo darbe yüksekliği spektrumu üretir ve FWHM ile R'yi çıkarmak için bir Gauss uydurur. Kaydırıcılar öğrencinin elektronik-gürültü ve kayma genişlemesi eklemesine izin verir; gösterilen toplam FWHM karesel toplam olarak güncellenir ve istatistik ile gürültünün ne zaman baskın olduğunu tam olarak gösterir. İki yakın tepe etkinleştirilerek ne zaman ayırt edilebilir hale geldikleri gösterilebilir.

GİRDİLER

Deposited energy E (10 keV - 3 MeV, log slider) · Detector type preset (gas / NaI / Si / Ge) setting W and F · Electronic noise FWHM (0 - 5 keV) · Drift FWHM (0 - 5 keV) · Two-peak separation toggle and spacing (0 - 50 keV)

ÇIKTILAR

Histogram pulse-height spectrum with fitted Gaussian overlay · Numeric N, sigma, FWHM, and R (percent) · Quadrature-sum bar chart of stat/noise/drift contributions · Resolved/unresolved indicator for the two-peak case

YÖNTEMCarrier count sampled from a Gaussian (or compound Poisson) with mean N and variance F*N; pulse height = K*N plus Gaussian noise and drift terms; histogram binned and least-squares fit to G(H)=A*exp(-(H-H0)^2/(2*sigma^2)); FWHM=2.355*sigma; R=FWHM/H0.

KOD KANCASIÖğrenci taşıyıcı-istatistiği örnekleyicisini ve karesel-birleşim fonksiyonunu yazar. F<1 seçmek tepeyi gerçek zamanlı daraltır; gürültü terimlerini yanlış birleştirmek (örn. karesel yerine doğrusal toplamak) gözle görülür şekilde yanlış tepe genişliği üretir ve karesel toplam kuralını doğrudan geri bildirimle öğretir.

SİMÜLASYON TASARIMIORTA SADAKAT

Dedeksiyon Verimi ve Katı Açı

Bir 3D sahne, yarıçapını ve mesafesini öğrencinin kontrol ettiği dairesel bir dedektöre doğru ışınlar yayan izotropik bir nokta kaynağı gösterir. Dedektöre çarpan ışınlar renklendirilir; yayılana karşı isabet sayısı geometrik (mutlak) verimi verir ve öğrenci bunu analitik katı-açı formülüyle karşılaştırır. Ardından bir içsel-verim kaydırıcısı (malzeme/kalınlık/enerji) isabetleri kaydedilen darbelere dönüştürür, geometriyi etkileşim fiziğinden ayırır ve kaynak geri çekildikçe ters-kare azalımı ortaya çıkar.

GİRDİLER

Source-to-detector distance d (1 - 100 cm) · Detector face radius a (0.5 - 10 cm) · Number of emitted rays (1e3 - 1e6) · Intrinsic efficiency (0 - 1, or material/energy preset)

ÇIKTILAR

3D ray-hit visualization (hits vs misses) · Monte Carlo geometric efficiency vs analytic Omega/4*pi · Absolute efficiency = intrinsic * geometric readout · Efficiency-vs-distance curve showing inverse-square limit

YÖNTEMIsotropic direction sampling (uniform in cos(theta) and phi); ray-disk intersection test to count hits; geometric efficiency = hits/emitted compared to Omega/4*pi with Omega=2*pi*(1 - d/sqrt(d^2+a^2)); absolute efficiency obtained by multiplying by the intrinsic factor.

KOD KANCASIÖğrenci izotropik yön örnekleyicisini ve ışın-disk kesişim testini uygular. Doğru bir küre-üzerinde-tekdüze örnekleyici Monte Carlo verimini analitik katı-açı değerine yakınsatır; yaygın bir hata (cos(theta) yerine theta'da tekdüze) sonucu saptırır ve öğrenci eğrinin analitik çizgiyi ıskaladığını görür.

Sayısal veriler

Typical gas-detector W-value (air)~34 eV per ion pairGazda bir elektron-iyon çifti oluşturmak için gereken ortalama enerji; iyon odaları ve orantılı sayaçlar için N'yi belirler.

Silicon W-value (e-h pair)3.62 eVSi'de elektron-boşluk çifti başına enerji; gazdan ~10 kat düşük olması ~10 kat fazla taşıyıcı ve çok daha iyi çözünürlük sağlar.

Germanium W-value (e-h pair)2.96 eVGe'de e-b çifti başına enerji; HPGe dedektörleri büyük N ve düşük Fano faktörü sayesinde keV-altı gama çözünürlüğüne ulaşır.

Fano factor, Ge and Si~0.06 - 0.12Birden epey düşük; yarıiletken çözünürlüğünün naif Poisson sqrt(N) ölçeklenmesini yenmesinin nedeni.

Fano factor, scintillators~1Sintilatörler Poisson'a yakın istatistik gösterir; çözünürlükleri düşük fotoelektron verimiyle sınırlıdır.

N for 1% resolution (Poisson)> 55,000 carriersR = 2.35/sqrt(N)'den; yüksek-taşıyıcı-sayılı dedektörleri motive eder.

Typical NaI(Tl) energy resolution at 662 keV6 - 8% FWHMSintilasyon gama spektroskopisi referans değeri (Cs-137 fotopiki).

Typical HPGe energy resolution at 1332 keV~0.13% (≈1.8 keV) FWHMCo-60 çizgisi; büyük N ve düşük F nedeniyle NaI(Tl)'den büyük mertebede daha iyi.

Recommended max dead-time loss< 30 - 40%Bunun ötesinde düzeltmeler model ve tau belirsizliklerine çok duyarlı hale gelir.

Two-source method target dead-time fractionm12*tau > ~20%Yöntem iki büyük ve neredeyse eşit sayımı farkettiğinden güvenilir tau için gereklidir.

Nonparalyzable saturation ratem_max = 1/tauGözlenen hızın asimptotik tavanı; tau = 1 us için bu 1e6 sayım/s'dir.

Paralyzable peak ratem_max = 1/(e*tau) at n = 1/tauGözlenen hız maksimumu; ötesinde sistem geri katlanır ve iki-değerli belirsizliği oluşturur.

YAZILIM NOTU

Bu bölüm, öğrencilerin dedektör fiziğini ilk kez gerçek hesaplamayla bağladığı yerdir. Çekirdek algoritmalar istatistikseldir: ters-CDF ile üstel dağılımlı olaylar-arası süreleri örnekleyerek (t = -ln(U)/n) sentetik bir olay akışı oluşturmak, ardından bu akışı felç olabilen ve felç olmayan ölü-zaman durum makinelerinden geçirerek m-vs-n eğrilerini yeniden üretmek ve analitik sayım-hızı düzeltmelerini doğrulamak. Çözünürlük için öğrenciler Fano-değiştirilmiş varyansla (F*N) taşıyıcı-sayısı örneklemeyi, Gauss tepe üretimini, histogram kutulamayı ve FWHM ile R'yi çıkarmak için en küçük kareler Gauss uydurmayı kodlar, sonra bağımsız genişleme kaynaklarını karesel olarak birleştirir. Verim için izotropik yön örneklemeyi (cos-theta'da tekdüze), ışın-geometri kesişim testlerini ve katı açının Monte Carlo tahminini uygular, kapalı-form disk formülüne karşı doğrularlar. Bunlar profesyonel araçlara doğrudan eşlenir: MCA spektrum toplama ve tepe-uydurma paketleri, sayım laboratuvarlarındaki ölü-zaman düzeltme rutinleri ve MCNP, GEANT4 ve PENELOPE gibi Monte Carlo taşıma kodlarının geometri/verim motorları; burada aynı izotropik örnekleme ve sayaç-verimi kavramları ölçeklenir. Sıcak döngüleri (olay üretimi, histogramlama, uydurma) WASM'a derlenen Rust'ta uygulamak tarayıcıda etkileşimli, neredeyse-yerel başarım sağlarken, TypeScript arayüzü, çizimi ve 3D görselleştirmeyi yönetir. Öğrenciler doğru örnekleme dağılımlarının ve doğru hata-birleştirme kurallarının isteğe bağlı ayrıntılar değil, güvenilir ile yanıltıcı nükleer ölçümler arasındaki fark olduğunu öğrenir.

ALGORİTMALAR

Inverse-CDF exponential inter-arrival samplingParalyzable/nonparalyzable dead-time state machinesMonte Carlo solid-angle estimationIsotropic direction sampling (uniform in cos-theta)Ray-disk intersection testCompound Poisson / Fano-modified carrier samplingHistogram binningLeast-squares Gaussian peak fittingQuadrature error combinationNewton/iterative inversion of m = n*exp(-n*tau)

DİLLER

Rust/WASMTypeScript

Pratiğe dökün

İlgili simülasyonu aç

← Önceki

Sayım İstatistiği ve Hata Öngörüsü

İyonizasyon Odaları

Ders 4 / 20Orta~40 dk kaldıGlenn F. Knoll (4th ed.)

Radyasyon Dedektörlerinin Genel Özellikleri

Basitleştirilmiş Dedektör Modeli

information carriersW-valuecharge collectioncurrent pulse i(t)active volumeenergy deposition

Dedektör Çalışma Modları

pulse modecurrent modemean square voltage (MSV) modeCampbelling moderesponse time Tfluctuating current component

Darbe Yüksekliği Spektrumları

differential pulse height spectrumintegral pulse height spectrummultichannel analyzer (MCA)channel binningfull-energy peakcontinuum

Sayım Eğrileri ve Platolar

counting curveintegral discriminatorplateauplateau slopeoperating pointthreshold setting

Enerji Çözünürlüğü

FWHMenergy resolution RFano factorPoisson limitGaussian response functionquadrature noise additionelectronic noise

Dedeksiyon Verimi

absolute efficiencyintrinsic efficiencytotal efficiencypeak (photopeak) efficiencysolid anglepeak-to-total ratiocounting efficiency

Ölü Zaman

dead time tauparalyzable modelnonparalyzable modelcount-rate correctiontwo-source methoddecaying-source methodrate saturation

Olay Başına Üretilen Yük

Temel Denklemler (23)

Bu denklemi çalıştır: Ölü Zaman Modelleri: Felç Olabilen vs Felç Olmayan →

Olay Başına Üretilen Yük

(4.1)

Q = N q = \frac{E}{W} q

sembolanlambirim

Q

olay başına toplanan toplam yükC

N

bilgi taşıyıcılarının sayısı–

E

aktif hacimde bırakılan enerjiJ

W

bir taşıyıcı oluşturmak için gereken ortalama enerjiJ

q

temel yük, 1.602e-19 CC

Akım Darbesinin İntegrali Olarak Yük

(4.2)

Q = \int_{0}^{t_c} i(t)\, dt

sembolanlambirim

i(t)

anlık dedektör akımıA

t_c

yük toplama süresis

Q

toplanan toplam yükC

Akım Modu Ortalama Sinyali

(4.3)

I_0 = r\,Q = r\,\frac{E}{W} q

Akım modunda kararlı durum akımı, olay hızı ile olay başına ortalama yükün çarpımına eşittir. Bireysel olayların ayırt edilmediği dozimetri ve yüksek-akı izlemede kullanılır.

sembolanlambirim

I_0

ortalama (DC) akımA

r

olay (etkileşim) hızı1/s

Q

olay başına ortalama yükC

Akım Modu Zaman Ortalaması

(4.4)

I(t) = \frac{1}{T}\int_{t-T}^{t} i(t')\, dt'

sembolanlambirim

I(t)

kaydedilen ortalama akımA

T

ölçüm aleti tepki (integrasyon) süresis

i(t')

anlık akımA

Akım Modunda Kesirli İstatistiksel Dalgalanma

(4.5)

\frac{\sigma_I}{I_0} = \frac{1}{\sqrt{2\,r\,T}}

sembolanlambirim

\sigma_I

ölçülen akımın standart sapmasıA

I_0

ortalama akımA

r

olay hızı1/s

T

tepki süresis

Ortalama Kare Gerilim (Campbelling) Sinyali

(4.6)

\overline{\sigma_i^2} \propto r\,Q^2

sembolanlambirim

\overline{\sigma_i^2}

dalgalanan akımın ortalama karesiA^2

r

olay hızı1/s

Q

olay başına yükC

Enerji Çözünürlüğü Tanımı

(4.7)

R = \frac{\mathrm{FWHM}}{H_0}

sembolanlambirim

R

enerji çözünürlüğü–

\mathrm{FWHM}

tepenin yarı yükseklikteki tam genişliği–

H_0

tepenin merkezi (ortalama darbe yüksekliği)–

Gauss FWHM

(4.8)

\mathrm{FWHM} = 2.355\,\sigma

sembolanlambirim

\mathrm{FWHM}

yarı yükseklikteki tam genişlik–

\sigma

Gauss standart sapması–

Poisson İstatistiksel Çözünürlük Sınırı

(4.9)

R_{\mathrm{Poisson}} = \frac{2.355\,K\sqrt{N}}{K N} = \frac{2.355}{\sqrt{N}}

sembolanlambirim

R_{\mathrm{Poisson}}

Poisson-sınırlı çözünürlük–

N

olay başına yük taşıyıcı sayısı–

K

doğrusal orantı sabiti (H_0 = KN)–

Fano Faktörü

(4.10)

F = \frac{\sigma_N^2}{N}

sembolanlambirim

F

Fano faktörü–

\sigma_N^2

taşıyıcı sayısındaki gözlenen varyans–

N

ortalama taşıyıcı sayısı–

Fano Faktörlü İstatistiksel Çözünürlük

(4.11)

R_{\mathrm{stat}} = 2.355\,\sqrt{\frac{F}{N}}

sembolanlambirim

R_{\mathrm{stat}}

istatistiksel çözünürlük sınırı–

F

Fano faktörü–

N

yük taşıyıcı sayısı–

FWHM Katkılarının Karesel Birleşimi

(4.12)

(\mathrm{FWHM})_{\mathrm{total}}^2 = (\mathrm{FWHM})_{\mathrm{stat}}^2 + (\mathrm{FWHM})_{\mathrm{noise}}^2 + (\mathrm{FWHM})_{\mathrm{drift}}^2 + \cdots

sembolanlambirim

(\mathrm{FWHM})_{\mathrm{total}}

genel tepe genişliği–

(\mathrm{FWHM})_{\mathrm{stat}}

istatistiksel (taşıyıcı) katkı–

(\mathrm{FWHM})_{\mathrm{noise}}

elektronik gürültü katkısı–

(\mathrm{FWHM})_{\mathrm{drift}}

kayma/kararsızlık katkısı–

Gauss Tepki Fonksiyonu

(4.13)

G(H) = \frac{A}{\sigma\sqrt{2\pi}}\,\exp\!\left[-\frac{(H-H_0)^2}{2\sigma^2}\right]

sembolanlambirim

G(H)

H darbe yüksekliğindeki diferansiyel sayım–

A

tepe alanı (toplam sayım)–

H_0

merkez darbe yüksekliği–

\sigma

standart sapma (genişlik)–

Mutlak Verim

(4.14)

\varepsilon_{\mathrm{abs}} = \frac{\text{number of pulses recorded}}{\text{number of quanta emitted by source}}

sembolanlambirim

\varepsilon_{\mathrm{abs}}

mutlak verim–

İçsel Verim

(4.15)

\varepsilon_{\mathrm{int}} = \frac{\text{number of pulses recorded}}{\text{number of quanta incident on detector}}

sembolanlambirim

\varepsilon_{\mathrm{int}}

içsel verim–

Mutlak-İçsel Verim İlişkisi

(4.16)

\varepsilon_{\mathrm{abs}} = \varepsilon_{\mathrm{int}} \cdot \frac{\Omega}{4\pi}

sembolanlambirim

\varepsilon_{\mathrm{abs}}

mutlak verim–

\varepsilon_{\mathrm{int}}

içsel verim–

\Omega

dedektörün kapsadığı katı açısr

Eşeksenli Diskin Katı Açısı

(4.17)

\Omega = 2\pi\left(1 - \frac{d}{\sqrt{d^2 + a^2}}\right)

sembolanlambirim

\Omega

katı açısr

a

dedektör yüzü yarıçapım

d

kaynak-dedektör mesafesim

Felç Olmayan Gözlenen Hız

(4.18)

m = \frac{n}{1 + n\tau}

Ölü zamanı tau olan felç olmayan bir sistem için gözlenen hız m ile gerçek hız n arasındaki ilişki. Gözlenen hız n büyüdükçe 1/tau'da doyar. Eşdeğer ifade: n - m = n*m*tau.

sembolanlambirim

m

gözlenen (kaydedilen) sayım hızı1/s

n

gerçek etkileşim hızı1/s

\tau

sistem ölü zamanıs

Felç Olmayan Gerçek-Hız Düzeltmesi

(4.19)

n = \frac{m}{1 - m\tau}

sembolanlambirim

n

gerçek etkileşim hızı1/s

m

gözlenen sayım hızı1/s

\tau

ölü zamans

Felç Olabilen Gözlenen Hız

(4.20)

m = n\,e^{-n\tau}

sembolanlambirim

m

gözlenen sayım hızı1/s

n

gerçek etkileşim hızı1/s

\tau

ölü zamans

Her İki Modelin Düşük-Kayıp Sınırı

(4.21)

m \approx n\,(1 - n\tau), \quad n\tau \ll 1

sembolanlambirim

m

gözlenen hız1/s

n

gerçek hız1/s

\tau

ölü zamans

İki-Kaynak Ölü Zaman Denklemi

(4.22)

\frac{m_{12}}{1 - m_{12}\tau} + \frac{m_b}{1 - m_b\tau} = \frac{m_1}{1 - m_1\tau} + \frac{m_2}{1 - m_2\tau}

sembolanlambirim

m_1

kaynak 1'in gözlenen hızı (arka plan dahil)1/s

m_2

kaynak 2'nin gözlenen hızı1/s

m_{12}

iki kaynak birlikte gözlenen hız1/s

m_b

gözlenen arka plan hızı1/s

\tau

ölü zamans

Bozunan-Kaynak Gerçek Hızı

(4.23)

n = n_0\,e^{-\lambda t}

sembolanlambirim

n

t anındaki gerçek hız1/s

n_0

ölçüm başındaki gerçek hız1/s

\lambda

izotopun bozunum sabiti1/s

t