PMT Measurements of PEN and Friends at MPP Oliver Schulz oschulz@mpp.mpg.de PEN Meeting, TU-Dortmund, LKT, March 12 2018
First results with new PMT setup at MPP (see previous talk) Source: Cs-137 DAQ: Struck SIS3316, 14 bit, 250 Msps Analysis pipeline: Sample dithering (triangular) Baseline estimation (Hamming window) Baseline subtraction (per event) Energy reco via simple integrals (sum of sample values): Short integral: 20 ns (5 samples), suitable for BC408 Long integral: 160 ns (40 samples), suitable for PEN
BC408 Pulse Shapes In [18]: waveforms_bc408 = av(parent(wfanalysis_bc408[:waveforms])[250:(250 + 12), :]); plot_pulse_hist(waveforms_bc408, -20:1:1500) Signal Shapes Out[18]: 1500 5.0 4.5 4.0 1000 3.5 3.0 Counts 2.5 2.0 500 1.5 1.0 0.5 0 2.5 5.0 7.5 10.0 12.5 Sample [4ns]
PEN Pulse Shapes In [9]: waveforms_pen = av(parent(wfanalysis_pen[:waveforms])[245:(245 + 40), :]); plot_pulse_hist(waveforms_pen, -20:1:300) Signal Shapes Out[9]: 300 5.0 4.5 4.0 200 3.5 3.0 Counts 2.5 100 2.0 1.5 1.0 0 0.5 10 20 30 40 Sample [4ns]
BC408 Pulse Shape Examples (same Integral) In [19]: waveforms_bc408_sel = wfanalysis_bc408[entrysel(E_dep -> 2000 < E_dep < 2010, wfan alysis_bc408, :peak_integral), :][:waveforms]; @manipulate for i in eachindex(waveforms_bc408_sel) plot(waveforms_bc408_sel[i][25 2:(252 + 5)]) end Out[19]: 800 y1 600 400 200 1 2 3 4 5 6
PEN Pulse Shape Examples (same Integral) In [11]: waveforms_pen_sel = wfanalysis_pen[entrysel(E_dep -> 600 < E_dep < 601, wfanalysis _pen, :peak_integral), :][:waveforms]; @manipulate for i in eachindex(waveforms_pen_sel) plot(waveforms_pen_sel[i][245:(2 45 + 40)]) end Out[11]: y1 80 60 40 20 0 10 20 30 40
BC408 In [12]: plots_bc408 = plot_wfanalysis(wfanalysis_bc408) Signal Shapes Peak Integral Out[12]: 100 5.0 10 4 4.5 80 Long window 4.0 Short window 10 3 60 3.5 Counts Counts 3.0 40 10 2 2.5 2.0 20 10 1 1.5 0 1.0 10 0 0.5 -20 100 200 300 400 500 0 2500 5000 7500 10000 Sample [4ns] ADC value Post-Pulse Baseline Level PSA vs. Peak Integral 2250 10 4 2000 PSA short/long integral 1.2 1750 10 3 0.9 1500 Counts 1250 0.6 10 2 1000 750 0.3 10 1 500 0.0 250 10 0 -2.5 0.0 2.5 5.0 0 2500 5000 7500 10000 Integral Peak integral (long)
PEN, LKT Exp. 4 In [13]: plots_pen = plot_wfanalysis(wfanalysis_pen) Signal Shapes Peak Integral Out[13]: 100 5.0 10 4 4.5 80 Long window 4.0 Short window 10 3 60 3.5 Counts Counts 3.0 40 10 2 2.5 2.0 20 10 1 1.5 0 1.0 10 0 0.5 -20 100 200 300 400 500 0 1000 2000 3000 Sample [4ns] ADC value Post-Pulse Baseline Level PSA vs. Peak Integral 1250 10 4 PSA short/long integral 1.2 1000 10 3 0.9 750 Counts 0.6 10 2 500 0.3 10 1 250 0.0 10 0 -4 -2 0 2 4 6 500 1000 1500 2000 2500 3000 3500 Integral Peak integral (long)
PET-G, 3D-Printed In [14]: plots_petg = plot_wfanalysis(wfanalysis_petg) Signal Shapes Peak Integral Out[14]: 100 4.0 80 10 3 3.5 Long window Short window 60 3.0 Counts Counts 10 2 2.5 40 2.0 20 10 1 1.5 1.0 0 10 0 0.5 -20 100 200 300 400 500 0 500 1000 1500 2000 Sample [4ns] ADC value Post-Pulse Baseline Level PSA vs. Peak Integral 10 4 250 PSA short/long integral 1.2 10 3 200 0.9 Counts 150 10 2 0.6 100 10 1 0.3 50 0.0 10 0 -4 -2 0 2 4 6 500 1000 1500 2000 Integral Peak integral (long)
PMMA Blank, same size at PET-G In [15]: plots_pmma = plot_wfanalysis(wfanalysis_pmma) Signal Shapes Peak Integral Out[15]: 100 4.0 10 3 3.5 80 Long window Short window 3.0 60 10 2 2.5 Counts Counts 40 2.0 10 1 20 1.5 1.0 0 10 0 0.5 -20 100 200 300 400 500 0 500 1000 1500 Sample [4ns] ADC value Post-Pulse Baseline Level PSA vs. Peak Integral 100 PSA short/long integral 1.2 10 3 0.9 75 Counts 10 2 0.6 50 10 1 0.3 25 0.0 10 0 -2 0 2 4 6 250 500 750 1000 1250 1500 1750 Integral Peak integral (long)
Energy Spectra (uncalibrated) In [17]: plot_wfanalyses_spectra(wfanalyses) Out[17]: BC408 PEN 0.8 PET-G PMMA Count rate [1 / s / energy-unit] 0.6 0.4 0.2 0.0 0 1000 2000 3000 E [arb. unit]
PET-G, 3D-Printed, UV-Excitation Spectra
Recommend
More recommend