MPPC 読み出しによる無機シンチレーターを用いた γ 線検出
井上裕貴

奥野あや
平成 28 年 3 月

日野純子

概要

MPPC(Multi Pixel Photon Counter) は、低バイアス電圧で高い増倍率を誇り、優れた時間分解能を持つ小型のフォトンカウンティングデバイスであり、光電子増倍管に代わる検出器として様々な分野で期待されている。また、GAGG は発光量が大きく潮解性と自己放射性を持たない無機シンチレーターである。放射性物質を利用し、放射線を光に変換するシンチレーターと、その光を検出する受光素子から構成される γ 線検出器には様々な応用例があり、医療機器分野における典型的なものの一つが PET(Positoron Emission
Tomography) 検査装置である。PET をはじめとして、GAGG シンチレーターと MPPCを組み合わせた γ 線検出器の特性評価とその考察を行った結果について記す。

目次
第1章
1.1
1.2
1.3
1.4
1.5
1.6

原子核崩壊から生じる γ 線とその検出
γ線 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
実験に用いた γ 線源 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
光電効果 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
コンプトン効果 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
電子対生成 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
シンチレーション光 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1
1
1
2
2
3
3

第 2 章 MPPC(Multi Pixel Photon Counter)
2.1 MPPC . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.1.1 APD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.1.2 プロポーショナルモードとガイガーモード . . . . . . . . . . . . . .
2.2 MPPC の動作原理 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.3 仕様 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.4 基本接続回路 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.5 動作に関して留意すべき点 . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.5.1 ダークパルス . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.5.2 クロストーク . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.5.3 アフターパルス . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4
4
5
5
6
7
8
9
9
9
9

第 3 章 シンチレーター
3.1 有機シンチレーター
3.2 無機シンチレーター

10. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

第 4 章 エレクトロニクス
4.1 ＮＩＭモジュール . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1.1 High Voltage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1.2 Clock generator . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1.3 Attenuator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1.4 AMP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

12
12
12
12
12
12

第5章
5.1
5.2
5.3
5.4

14
14
14
16
17

MPPC の動作確認
原理 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
セットアップ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
実験方法 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
生データのフォーマット . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1

5.6

オフライン解析 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.5.1 段階１：信号電荷への変換 . . . . . . . . . . . . . . . . . . . . . . .
5.5.2 段階２：波高分布のフィット . . . . . . . . . . . . . . . . . . . . . .
5.5.3 段階３：増倍率の算出 . . . . . . . . . . . . . . . . . . . . . . . . .
結果・考察 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

18
18
19
20
21

第6章
6.1
6.2
6.3
6.4

GAGG シンチレーターとの組み合わせによる γ 線検出
セットアップ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
実験方法 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
波高分布 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
結果・考察 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

23
23
24
24
25

5.5

第 7 章 まとめ

29

謝辞

30

付 録 A 増倍率測定時に用いたデータ収集プログラムのソースコード

31

付 録 B 増倍率データ解析プログラム

33

付 録 C 検出光子数のシミュレーション

35

参考文献

36

2

図目次
1.1
1.2

137

Cs の壊変図 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
コンプトン散乱 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1
2

2.1
2.2
2.3
2.4
2.5
2.6
2.7
2.8

浜松ホトニクス社製 S10362-11-050C 型 MPPC の外観 . . . . . . . . . . . .
空乏化 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
MPPC の 1 ピクセル . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
MPPC の構造 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
１ピクセルの光子入射後の動作サイクル . . . . . . . . . . . . . . . . . . .
MPPC 外寸図 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
基本接続回路 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
ダークパルスをオシロスコープで観測した波形 . . . . . . . . . . . . . . . .

4
5
6
6
6
7
8
9

3.1
3.2

有機シンチレーターの発光原理 . . . . . . . . . . . . . . . . . . . . . . . . 10
無機シンチレーターの発光原理 . . . . . . . . . . . . . . . . . . . . . . . . 11

4.1

アンプの増幅率 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

5.1
5.2
5.3
5.4
5.5
5.6

15
15
15
15
15

5.9

LED 点灯実験の Block Diagram . . . . . . . . . . . . . . . . . . . . . . . .
デジタルオシロスコープ . . . . . . . . . . . . . . . . . . . . . . . . . . . .
アルミシャーシ内の様子 . . . . . . . . . . . . . . . . . . . . . . . . . . . .
恒温槽 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
NIM モジュール . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
サンプル 616、温度 25 ℃、逆電圧 70V、LED 点灯をトリガーした時の波形
スナップショット . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
サンプル 616、温度 25 ℃、逆電圧 70V の時の信号電荷の分布 . . . . . . . .
サンプル 616、温度 25 ℃、逆電圧 70V の時の信号電荷の分布における個々のピークをガウス分布でフィットしたもの . . . . . . . . . . . . . . . . . .
温度 25 ℃時の増倍率の印加電圧による変化 . . . . . . . . . . . . . . . . .

6.1
6.2
6.3
6.4
6.5
6.6
6.7
6.8
6.9

Block Diagram . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
上から見た図 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
アルミシャーシ内のセットアップ . . . . . . . . . . . . . . . . . . . . . . .
γ 線入射時の波高分布 . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
反射材のテフロンシートをかぶせた場合の想定図 . . . . . . . . . . . . . .
テフロンを被せていないとき . . . . . . . . . . . . . . . . . . . . . . . . .
テフロンを被せたとき . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
MPPC の検出効率 (浜松ホトニクス社のカタログより抜すい) . . . . . . .
各ピクセルに入射した光子数のヒストグラム . . . . . . . . . . . . . . . .

23
24
24
24
25
25
25
26
28

5.7
5.8

3

17
19
19
22

6.10 光子を検出できたピクセル数のヒストグラム

4. . . . . . . . . . . . . . . . 28

表目次
2.1

MPPC の仕様 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

3.1
3.2

有機シンチレーターの例 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
無機シンチレーターの例 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

5.1
5.2
5.3
5.4

サンプル 606 の温度 25 ℃時の電圧変化による増倍率 . . . . . . . . . . . . .
サンプル 608 の温度 25 ℃時の電圧変化による増倍率 . . . . . . . . . . . . .
サンプル 616 の温度 25 ℃時の電圧変化による増倍率 . . . . . . . . . . . . .
図 5.9 における各 MPPC の近似式。p0：切片、p1：傾き . . . . . . . . . .

5

7

21
21
21
22

第 1 章 原子核崩壊から生じる γ 線とその
検出
1.1

γ線

γ 線は波長の極めて短い電磁波で、ある種の放射性同位元素の原子核がエネルギー準位間の遷移を起こすと放出される。波長が短いことからγ線は２つの特徴を示す。１つ目
は、粒子としての性質を強く持つことである。２つ目は、透過力が強いことである。その透過力の強さから、γ 線は人体の深くまで到達できるため、癌の放射線治療などにも使用される。また、透過性が強い事、荷電粒子ではないことから、γ線は直接的に物質を電離することはほとんどないが、以下に述べるようなで示す様な過程を経ることで、物質中か
ら電子（二次電子）をたたき出すため、この電子が物質を電離する作用を持つ。

1.2

実験に用いた γ 線源

Cs は半減期 30.2 年で、137 Cs は β − 崩壊により、94 ％が励起状態の 137 Ba、残りの 6％が基底状態の同位体 137 Ba へと変化する。励起状態の 137 Ba はエネルギー 662keV のγ線を放出して、基底状態の 137 Ba となる。
137

図 1.1: 137 Cs の壊変図

1

1.3

光電効果

物質に光が入射すると、その光子のエネルギーが物質中の電子に全て与えられ、その電子（光電子）が放出される現象である。このとき、光子はある一定の値 ϕ を上回るエネ
ルギーを持つ必要があり、ϕ の値は物質によって異なる。光電子の最大の運動エネルギーは、照射される光の振動数によってのみ決まり、単位時間当たりに飛び出す光電子の個数は、照射する光の強度による。
1 個の光子が 1 個の電子に当たり、その電子が放出されるとき、その電子のエネルギー
E は、入射する光子のエネルギーを Ephoton とすると

E = Ephoton − ϕで与えられる。

1.4

コンプトン効果

光子が束縛の緩い軌道電子や自由電子に衝突すると、その電子を散乱することでエネル
ギーを失い、散乱光子の波長が長くなる現象をコンプトン効果という。
プランク定数を h、散乱前の光子のエネルギーを hν0 、散乱された光子のエネルギーをhν1 とすると、前後のエネルギーの差E = hν0 − hν1は、散乱された電子 (コンプトン電子）の運動エネルギーとなる。

図 1.2: コンプトン散乱

2

1.5

電子対生成

1 個の光子（γ 線) が物質中で消えて、電子とその反粒子である陽電子が生成する過程
を電子対生成という。γ 線が原子核近傍などの電磁場が強い自由空間でないところに入射（衝突）すると、γ 線のエネルギー hν が電子と陽電子に与えられるので、電子の質量を me
として 2me c2 =1.02MeV のエネルギーが必要である。
この過程で作られた電子は、物質中を進んでいく間に次第に運動エネルギーを失い、最終的には止まってしまう。一方、陽電子は物質中に存在するの電子と対消滅する。対消滅が起こる直前には、陽電子が電子と互いに相手の周りをまわる束縛状態を作る。このような系のことをポジトロニウムという。

1.6

シンチレーション光

物質中を荷電粒子が通過すると、物質中の原子（あるいは分子）が励起、つまり高いエ
ネルギー状態に移る。これらの励起状態の原子は、すぐに低いエネルギー状態に移ろうとする。その際、余分になったエネルギーは光としてそして放出される。この光がシンチ
レーション光である。 シンチレーション光を効率よく出す素材をシンチレーターと呼ぶ。
本研究で用いたシンチレーターについては第３章で詳しく述べる。

3

第 2 章 MPPC(Multi Pixel Photon
Counter)
2.1

MPPC

MPPC とは、Silicon Photomultiplier とも呼ばれるデバイスの一種で、Multi Pixel Photon Counter の略称である。浜松ホトニクス社が開発した光検出半導体素子であり、優れたフォトンカウンティング能力をもち、微弱光を検出する様々な用途に利用できる。特徴は 100V 以下の低電圧で動作し、ガイガーモードで動作する APD ピクセルを持つことにより 105 〜106 におよぶ高い増幅率を持つことである。また、半導体であるため磁場の影
響を受けず、個々のピクセルの静電容量が小さいため良好な時間分解能を示す。本実験で使用した MPPC（浜松ホトニクス社製 S10362-11-050C 型）の画像を図 2.1 に記す。

図 2.1: 浜松ホトニクス社製 S10362-11-050C 型 MPPC の外観

4

2.1.1APD

APD（Avalanche Photodiode）を含め、PN 接合のダイオードでは P 層側が N 層側より低い電位になるように電圧を印加すると、電子-正孔対の電子は N 層に、正孔は P 層に集まる。このためダイオードの内部で電子または正孔のいずれのキャリアも存在しない空乏化した領域が形成され、素子内の電流が止まる。（図 2.2）この空乏層に光子が入射し、光電効果を起こして電子-正孔対が生成されると、素子内部の電位差により電子が N 層に、正孔が P 層に向かって流れて電流が発生するので、これを信号パルスとして取り出す。これがフォトダイオードの基本的な原理である。
APD では、不純物濃度や構造を工夫することにより、素子の内部に局所的に強い電場を作るようにする。すると、入射した光が光電効果を起こして生成した電子が局所的に強い電場が作られている場所を通るとき、半導体内の電子に衝突、新たな電子-正孔対を生成するに十分なエネルギーを得る。このため、電子-正孔対の連鎖的な生成が可能になる結果、電子雪崩が形成される。APD は、こうして電子雪崩を利用して信号を増幅する機能を持つフォトダイオードである。

図 2.2: 空乏化

2.1.2

プロポーショナルモードとガイガーモード

APD 素子中の構造および印加電圧により、電子雪崩による増幅が数十倍か 100 倍程度であるとき、形成された電子な誰は素子中（あるいはピクセル中）の限定された領域に形成される。そのため素子受光面に複数の光子が入射して電子‐正孔対を作ると、独立に電子雪崩が発達する。したがって、最終的な信号電荷の大きさが入射光子の数に比例する。
これをプロポーショナルモードという。
一方、降伏電圧を超えた電圧のときにガイガーモードで動作する。ガイガーモードでは、形成された電子雪崩がピクセルの全面に広り、これをガイガー放電と呼ぶ。これにより、出力される信号パルスは入射して光電効果を起こした光子数に関係なく、ほぼ一定の出力となる。したがって、信号電荷の大きさは入射光子数の情報を失っている代わり、105
倍程度に達する大きな増幅率を得ることが可能である。

5

2.2

MPPC の動作原理

MPPC は、個々のピクセルが独立なガイガーモードの APD としてはたらき、その各々にクエンチング抵抗が直列に接続されたものをピクセルの数だけ並列接続した構造をしている（図 2.3 および図 2.4）。光子が入射したピクセルで、ガイガーモードの電子雪崩を形成し信号パルスを作る。このときの電流によるクエンチング抵抗の電圧降下によって逆バイアスが降伏電圧 VBR 以下に下がり、電子雪崩が停止する。その後再充電され、逆バイ
アスがもともとの印加電圧 VR まで戻ると再び入射する光子を検出可能な状態になる（図2.5）。受光部に入る粒子の数がピクセルの数より十分少ない条件下で使用すれば、信号を出したピクセルの数を検出した光子の数とみなしてよい。ガイガーモードで信号を出した
ピクセル 1 個あたり 105 個の電子がパルス出力されるので、ピクセルに一光子入射した信号の大きさは十分雑音から分離するのでフォトンカウンティングが可能である。このとき出力信号の電荷はとびとびの値を取るという特徴がある。

図 2.4: MPPC の構造
図 2.3: MPPC の 1 ピクセル

図 2.5: １ピクセルの光子入射後の動作サイクル

6

2.3

仕様

今回使用した MPPC は、浜松ホトニクス社製 S10362-11-050C である。表 2.1 に仕様を、
図 2.6 に外寸図を示す [1]。
表 2.1: MPPC の仕様
項目
有効受光面サイズ
ピクセル数
ピクセルサイズ
動作温度
開口率
感度波長範囲最大感度波長検出効率
動作電圧範囲ダークカウント
端子間容量
時間分解能
逆電圧の温度計数
増倍率

1× 1
400
50 × 50-20〜+40
61.5
320〜9004405070 ± 10400
35200〜300
56
7.5 × 105

図 2.6: MPPC 外寸図

7

単位
mm2
‐
µm2
℃
％
nm
nm
％
V
kcps
pF
ps
mV/℃
-

2.4

基本接続回路MPPC の基本接続回路を図 2.7 に示す。この回路において、47kΩ の抵抗と 0.1µF のコンデンサーによるローパスフィルター２段を配置して、バイアス電源のリップルやノイズを遮断する。さらに信号の読み出しは 0.47µF の結合コンデンサー接続して、バイアス電源の直流成分をカットし、AC 結合で信号を読み出す。使用したコンデンサはいずれも耐圧 100V のセラミック製で、TDK 社の製品である。

図 2.7: 基本接続回路

8

2.5

動作に関して留意すべき点

2.5.1

ダークパルスMPPC では、電子-正孔対は入射光子による光電効果だけでなく、熱励起によっても発生する。熱励起によって発生した電子-正孔対もガイガーモードによって光子が入射した
場合と同様の信号パルスが形成される。これをダークパルスと呼ぶ。図 2.8 に外部からの光を遮断し、MPPC に 70V の電圧を印加し、オシロスコープによって観測された波形を示す。2.5GSample/sec のレートで光子の検出を行っている。

図 2.8: ダークパルスをオシロスコープで観測した波形

2.5.2

クロストークAPD ピクセルにおいて、電子雪崩の過程で二次的な光子が放出される。放出された光子が別のピクセルで光電効果を起こし、信号パルスを出力することがある。光電効果を起こした二次光子が一つなら、まるで 2 つの光子が同時に 2 つのピクセルに入射したかのよ
うに、2 倍の波高が検出される。それより確率は下がるが 2 つの二次光子が他のピクセルで相互作用したときは 3 倍の波高になる。こうして、信号を出力したピクセルの数が実際に検出した光子数より大きくなることがある。この現象をクロストークと呼ぶ。

2.5.3

アフターパルス

APD ピクセル内で電子雪崩を生成した際に、発生した電子の一部が半導体内の格子欠
陥にトラップし、遅れて解放されて電子雪崩を作ってパルスを出す。これをアフターパル
スと呼ぶ。温度が低くなるにつれ電子が格子欠陥にトラップされる確率が高くなり、アフ
ターパルスの発生頻度は増加すると考えられている。

9

第 3 章 シンチレーター
シンチレーターとは、その内部で入射した粒子がエネルギー損失を起こしたとき、それにより蛍光を発する物質のことである。シンチレーターには、大別して有機シンチレー
ターと無機シンチレーターの 2 種が存在し、それぞれの発光の仕組みは大きく異なっている。本研究では、MeV 程度の γ 線を検出する目的に向く無機シンチレーターを用いた。

3.1

有機シンチレーター

有機シンチレーターは、π 電子を持つ有機分子の状態遷移を発光の源とするものが大多数である。図 3.1 は有機分子中の π 電子エネルギー準位図である。S は電子のスピンが 0の シングレット状態、T はスピンが 1 の トリブレット状態であり、添字の 0 は基底準位、
1,2,3,… は励起準化を表し、実線はその最低の準位を表す。入射した放射線により基底準
位 (S0) にある有機分子がそれぞれの励起準位（S1,S2,S3）に励起され、全ての励起準位から第１励起状態（S1）の最低準位に非常に短い時間で遷移する。そして、基底状態 (S0)に遷移して蛍光が発する。蛍光減衰時間が比較的短いため、良好な時間分解能を得ることができる。表 3.1 に有機シンチレーターの例をあげる。[2]図 3.1: 有機シンチレーターの発光原理

表 3.1: 有機シンチレーターの例
Anthracene Plastic(NE 102A)
3
密度 [g/cm ]
1.25
1.32
屈折率
1.62
1.58
蛍光減衰時間 [ns]322.4発光波長 [nm]
447
423

10

Liquid(NE213)
0.87
1.51
3.7
425

3.2

無機シンチレーター

図 3.2 に示すように、入射した放射線のエネルギー損失により励起した結晶中の電子が脱励起する際に光子を放出する。この過程をシンチレーションと呼ぶ。無機シンチレー
ターは有機シンチレーターに比べ、原子番号の大きな元素を含んで密度の高いものが多い
ため、γ 線の検出効率を高くできる。また、発光量が豊富なものが多く、高いエネルギー
分解能を実現可能であるという利点がある。一方、無機シンチレーターは一般的に有機シンチレーターに比べ、蛍光減衰時間が長く、時間応答性が悪くなるため高速の測定に不向きであるというデメリットがある。温度が上昇すると、熱エネルギーが増加するため励起状態の分布が空間的、エネルギー的に広がる。これにより、励起された電子が基底状態に戻る時間が短くなるので、蛍光減衰時間が短くなるものも多い。表 3.2 に無機シンチレータの例を挙げる。[3] 本実験では以下に挙げる 5mm × 5mm 角のシンチレーター GAGG を用いた。LSO や LYSO、LuAG といった無機シンチレーターはルテチウム（Lu) を含む。
Lu は放射性同位元素の存在比が高く、シンチレーター中にバックグラウンド源を含有することになってしまう.。そこで、Lu を使わない酸化物シンチレーターとして Gd、Ga、Al の酸化物をホスト物質とし、発光する添加物として Ce をドープした GAGG が開発された。

図 3.2: 無機シンチレーターの発光原理

表 3.2: 無機シンチレーターの例
GAGG LSO BSO
3
密度 [g/cm ]6.63
7.40 6.80
発光量 [photon/MeV] 60000 26000 900
蛍光減衰時間 [ns]
88
40
100
発光波長 [nm]
520
420
480

11

第 4 章 エレクトロニクス
本実験において、ＭＰＰＣへの逆バイアス電圧印加、トリガー論理信号の生成にはＮＩ
Ｍ規格のエレクトロニクスを使用した。それについて、本章で記す。

4.1

ＮＩＭモジュール

ＮＩＭモジュールとはＮｕｃｌｅａｒ Instrument Modules の略称である。1960 年代に米国原子力委員会（ＡＥＣ）において制定された「放射線測定モジュール標準規格 Ｔ
ＩＤ−２０８９３」のことをＮＩＭ規格と言い、それに準拠して製作されたモジュールのことをＮＩＭモジュールという。この規格を適用したものであれば、国やメーカーやによらず物理的、電気的に完全な互換性を持つ。このＮＩＭモジュールを装着する箱状の標準ビンのことをＮＩＭビンといい、ＮＩＭビン電源は装着したＮＩＭモジュールに電源を供
給する。

4.1.1

High Voltage

High Voltage は高電圧を印加する電源モジュールである。本実験においては、半導体検出器向けの最大＋ 100V まで出力できるサトウ電子工業製のものを使用した。

4.1.2

Clock generator

Clock generator は、出力パルス幅を可変でき、Fast NIM 信号と TTL 信号を出力するモジュールである。本実験においては、LED のパルスを点灯する際に使用した。

4.1.3

Attenuator

Attenuator は、入力信号を減衰させるモジュールである。本実験においては、アンプの増幅率を測定する際に使用した。

4.1.4

AMP

AMP は、光電子増倍管や MPPC のように応答が速い光検出器の入力信号を線形増幅するモジュールである。本実験において、MPPC は増倍率が十分に高いため光電子増倍管の入力信号を線形増幅する既製の PM AMP を用いて目的を達成する。
AMP の増幅率は Clock generator を接続したときの、AMP 使用時と未使用時のオシロス
コープ画面の比較により求めた。
12

図 4.1: アンプの増幅率
図 4.1 のように、100 番目までのデータをそれぞれ基準面とし、基準面からのパルスの高さが AMP 使用時は平均 80.3 であるのに対し、AMP 未使用時は平均 8.8 であることから

80.3 ÷ 8.8 = 9.1となり、AMP の増幅率を 9.1 と求めた。

13

第 5 章 MPPC の動作確認
MPPC の動作確認のために 3 つのサンプルについて、それぞれの増倍率を測定した。十分に減光した LED のパルスから１パルス当たり数個の光子が MPPC の受光部に届く状態にし、MPPC の出力波高分布を得ることにより測定した。信号の数値化にはデジタル
オシロスコープを使用した。

5.1

原理

増倍率は、１ピクセルあたり１フォトンの信号を検知した時の電荷量を、１電子あたりの電荷量で割った値である。MPPC の出力パルスは検出した光子の数に応じて離散的な波高になるため、１フォトン検出時の信号の電荷量は隣り合う２つの波高分布のピークの差から求められる。つまり、増倍率=

２つのピークの電荷量の差
１電子あたりの電荷量(5.1)

で求まる。MPPC には温度特性がある。温度が上がると結晶の格子振動が激しくなり、
光電効果で発生したキャリアのエネルギーが十分大きくならないうちに、結晶中のフォノンと衝突する確率が高まる。これが、電子雪崩の発生と発達を抑制する向きにはたらくた
め、同じ逆電圧では高温になると増倍率が小さくなる。一定の温度なら、逆電圧を大きくすれば増倍率は大きくなる。したがって、一定の増倍率で動作させるためには、素子の温度を一定に保つ必要がある。

5.2

セットアップ

図 5.1 に LED 点灯実験の Block Diagram、図 5.2 から図 5.5 に実験に使用したハードウェ
アを示す。
LED 回路と MPPC 回路はアルミシャーシの中に設置し、アルミシャーシは恒温槽に入れて温度を一定に保つ。２つの回路の間には、LED を点灯させるために Clock generatorから送ったパルスが MPPC の読み出し系に対して雑音とならないよう、両社の間を静電遮蔽するために段ボール表面にアルミテープを貼った壁を置き、LED の光が MPPC に届
くように壁に 2mm 程度の穴を開けた。LED 回路は Clock generator と LED の間に直列に100 Ωの抵抗を入れ、TTL パルスを送って LED をパルス点灯させた。MPPC の逆バイア
ス電圧印加と信号読み出しには、図 2.7 基本接続回路を使用した。

14

図 5.1: LED 点灯実験の Block Diagram

図 5.2: デジタルオシロスコープ
図 5.3: アルミシャーシ内の様子

図 5.4: 恒温槽
図 5.5: NIM モジュール15

• PC 環境： CPU:AMD Athlon(tm)• LED：日亜化学工業株式会社製

X2 220 processor OS:LinuxNSPB320BS• デジタルオシロスコープ：Tektronix 製• 恒温槽：日本ブロアー株式会社製• Bias：サトウ電子工業製

DPO 3034 Dihital Phosphor Oscilloscope

LS-5 バイオチェンバー

HV-07WS Dual High Voltage Power Supply• Clock generator：株式会社テクノランドコーポレーション製
Clock Generator• AMP：株式会社海津製作所

5.3

N-TM 203 100MHz

KM2107 12ch PMT AMP

実験方法

Clock generator から LED 回路に周波数 10Hz で TTL パルスを送ることで LED を点灯させる。LED に拡散キャップをつけて光を拡散させ、薄い青と白の紙製の箱を LED 全体を覆うように被せることで十分に減光した。また、段ボールにアルミテープを張り２ mm程度の穴を開けた仕切りを作成し、LED 回路と MPPC 回路の間に置いて更に光量を絞る。
アルミボックス、LED 回路、MPPC 回路、仕切りは一点接地して静電遮蔽した。
デジタルオシロスコープは LAN インターフェイスを装備しており、これを介して PCとの命令やデータの送受信が可能である。 LXI と呼ばれる通信プロトコルに基づいてコ
マンドやデータの送受信を行なうソフトウェアパッケージとして VXI11 がある。 VXI11がサポートする関数をコールするデータ収集用プログラムを作成、実行することで、デジタルオシロスコープよりダイナミックレンジ 8bit、最大サンプリング周波数 2.5GHz の波形データを取得できる。このデータをテキスト形式で PC のディスク上のファイルに生データとして書き出して、オフライン解析処理を進める。データ収集用プログラムのソー
スコードの一例を付録 A に示す。
MPPC は温度を 25 ℃に保ち、逆電圧を 70V から 0.5V 刻みで 71.5V まで変化させて増倍率の変化を測定した。

16

5.4

生データのフォーマット

以下に、例として製造番号 616 の MPPC サンプル（以後サンプル 616 と略記）を、温度 25 ℃、逆電圧 70V の条件で動作させた際の生データの一例、図 5.6 にオシロスコープ
画面上のスナップショットを示す。

• 1 行目は、ラン開始レコード (begin run record) にあたり、デジタルオシロスコープ
が 2.50GSample/s、すなわち時間間隔が1s= 0.4ns
2.50GSample(5.2)

で収集されたデータであることを示す。

• 2 行目は、最初のイベントレコード (event record) で、イベント番号、サンプリング
数を表し、そのあと 8bit のデータが 1000 個分続く。
• その後、イベントレコードがファイルの末尾まで繰り返される。

図 5.6: サンプル 616、温度 25 ℃、逆電圧 70V、LED 点灯をトリガーした時の波形スナッ
プショット

17

5.5

オフライン解析

5.5.1

段階１：信号電荷への変換

5.4 で得た生データを、付録 B 増倍率データ解析用プログラムで処理をする。以下に、例として MPPC サンプル 616、温度 25 ℃、逆電圧 70V の時の生データを、プログラムで処理したデータを示す。

• １列目は、イベント番号を示す。
• ２列目は、ゼロ点の値を示す。ここでゼロ点は、triger position を左から 10 ％、つまり 100 番目に設定しているため、100 番目までの平均と定義する。
∑100
ゼロ点 =

i=1 xi(5.3)

100

ここで、i はサンプリング数、xi は i 番目のデータである。

• ３列目は、ゼロ点の分散を示す。
∑100
分散 =

i=1 (ゼロ点 − xi )

2

100

(5.4)

triger position より前でパルスが出ている時は分散が大きくなるため、指標として分散を計算する。
• ４列目は、ゼロ点からピークまでのパルスハイトを示す。
• ５列目は、信号パルスの積分値（信号電荷）を示す。積分値は、信号パルスが出ていると思われる 110 番目から 300 番目のデータを、ゼロ点からそれぞれ引いたものを積算する。

積分値 =

300
∑
i=110

18

(ゼロ点 − xi )(5.5)

5.5.2段階２：波高分布のフィット

1. 段階１で得た信号電荷、すなわち波高データを読み込み、得た分布の一例を図 5.7に示す。検出した光電子の個数（photoelectron=p.e.）に応じて等間隔にピークが立つ分布を示す。
2. 各 p.e. ピークの中央値を Maximum Likelihood 法によりガウス分布でフィットして求める。（図 5.8）
3. 隣り合う２つのピークの間隔の平均と分散を求める。

図 5.7: サンプル 616、温度 25 ℃、逆電圧 70V の時の信号電荷の分布

図 5.8: サンプル 616、温度 25 ℃、逆電圧 70V の時の信号電荷の分布における個々のピー
クをガウス分布でフィットしたもの

19

5.5.3

段階３：増倍率の算出

式（5.1）より MPPC の増倍率を求める。ここで、２つのピークの電荷量の差 =

２つのピークの間隔の平均× LSB 電圧×データの時間間隔
デジタルオシロスコープの抵抗値×アンプの増倍率(5.6)

1div 当たりの電圧× 8[div]28
データの時間間隔 = 0.4[ns]LSB 電圧 =アンプの増幅率 = 9.1
デジタルオシロスコープの抵抗値 = 50[Ω]
１電子当たりの電荷量 = 1.602 × 1019 [C]

(5.7)
(5.8)
(5.9)
(5.10)(5.11)である。これらの値を式（5.1）に代入することにより、電圧変化による MPPC の増倍率を計算する。

20

5.6

結果・考察

温度を 25 ℃に保った恒温槽において、製造番号 606,608,616 の 3 つの MPPC サンプルについて、逆電圧を 70.0V から 71.5V まで 0.5V 刻みで変化させて測定した結果を以下に示す。表 5.1〜表 5.3 に波高分布のフィットの結果より求めた平均ピーク間隔と、5.5.3 より求めた増倍率を示す。図 5.9 に横軸に印加電圧、縦軸に増倍率をとったグラフを示す。
表 5.1: サンプル 606 の温度 25 ℃時の電圧変化による増倍率

HV[V]平均ピーク間隔
増倍率 (× 105 )

70.0
853.2 ± 7.2
7.3 ± 0.2

70.5
1072.2 ± 12.3
9.2 ± 0.1

71.0
71.5
1387.2 ± 11.5 1418.2 ± 18.1
11.9 ± 0.9
12.2 ± 0.6

表 5.2: サンプル 608 の温度 25 ℃時の電圧変化による増倍率

HV[V]平均ピーク間隔
増倍率 (× 105 )

70.0
70.5
71.0
71.5
818.6 ± 5.3 986.2 ± 7.8 1249.6 ± 12.3 1294.6 ± 16.0
7.0 ± 0.2
8.5 ± 0.1
10.7 ± 0.7
11.1 ± 0.4

表 5.3: サンプル 616 の温度 25 ℃時の電圧変化による増倍率

HV[V]平均ピーク間隔
増倍率 (× 105 )

70.0
70.5
71.0
748.9 ± 4.0 1021.5 ± 9.1 1270.0 ± 13.5
6.5 ± 0.3
8.8 ± 0.2
10.9 ± 0.6

21

71.5
1355.3 ± 7.4
11.6 ± 0.5

図 5.9: 温度 25 ℃時の増倍率の印加電圧による変化

表 5.4: 図 5.9 における各 MPPC の近似式。p0：切片、p1：傾き
サンプル
606
608
616

χ2 /ndfp0
1.152/2 (-236.06 ± 48.03) × 105
0.615/2 (-195.85 ± 35.09) × 105
0.712/2 (-236.76 ± 37.76) × 105

p1(3.48 ± 0.68) × 105(2.90 ± 0.50) × 105(3.48 ± 0.53) × 105

図 5.9 より、同じ電圧を印加したときの３つの MPPC の増倍率個体差は大きくても相対値で 10 ％程度であることがわかった。印加電圧による増倍率の変化は、測定した 3 つのサンプルについては誤差の範囲で一致した。

22

第 6 章 GAGG シンチレーターとの組み合わせによる γ 線検出
Cs 線源を用いて、GAGG シンチレーターを光学グリースを用いて MPPC 受光部にとりつけ、137 Cs が発する 662keV の γ 線を検出する実験を行った。
137

6.1

セットアップ

下記に本実験のセットアップを示す。Bias 電源、MPPC が出力する信号の読み出しは前章で説明したものと同じである。ただし、LED 点灯時には LED を光らせる Clock generator
でデジタルオシロスコープをトリガーしていたが、この測定では MPPC が出力する信号
でトリガーするセルフトリガーのセットアップで行った点が異なる。線源が発する γ 線が
シンチレーターに届くよう、アルミシャーシ天板の対応する位置に直径 6mm の穴をあけた。図 6.2 はアルミシャーシを上から見た写真であり、図 6.3 はアルミシャーシ内のセットアップである。これらを恒温槽内に入れ、温度を 25 ℃に保って測定した。

図 6.1: Block Diagram23図 6.2: 上から見た図

6.2

図 6.3: アルミシャーシ内のセットアップ

実験方法

前小節で説明したセットアップを用いて、MPPC サンプル 616 を使用し、逆電圧は 70.0V、恒温槽で 25 ℃に保った。信号電荷すなわち波高を求める際に、トリガー位置は波形デー
タ 1000 点中の 100 番目付近にしていることと、GAGG シンチレーターの発光減衰時間が約 90ns であるので、パルス終了までの時間がダークパルスや LED 点灯時より遅いことを考慮して、積分するデータサンプルを波形データの終端である 1000 番目までとり入れた。

6.3

波高分布

Cs 線源の有無に対応して、出力信号パルスの有無をオシロスコープで確認した上で、Cs 線源の γ 線が入射した時の波高分布を図 6.4 に示す。

137
137

図 6.4: γ 線入射時の波高分布
波高分布の主成分はガウス分布に近い形をしている。そこでこの分布のピークをガウス分布でフィットして平均値を求めると 2.88 × 104 となっていた。LED 点灯により求めた

24

単一光電子検出時の波高をこのデータとオシロスコープの使用レンジの違いを考慮して換算した値 3.80 × 102 と比較すると2.88 × 104= 75.83.80 × 102となり、平均約 76 個の光子が検出できたということになる。

6.4

結果・考察

前小節で得た結果をさらに吟味・考察するために行った追加測定と見積もりについて述べる。
ここまでの測定ではシンチレーターを裸の状態で MPPC にとりつけていた。この場合、シンチレーション光のうち MPPC 受光部以外へ向かったものはそのまま外部へ逃げて検出されない可能性がある。そこで、反射材のテフロンシートでシンチレーターを包み、シンチレーションが反射してから MPPC の受光面に届く寄与を検討することにした。図 6.5にこの際の想定図をあげる。

図 6.5: 反射材のテフロンシートをかぶせた場合の想定図

図 6.6: テフロンを被せていないとき

図 6.7: テフロンを被せたとき

結果の波高分布を、図 6.6 と図 6.7 に示す。波高分布のピークをガウス分布でフィット
した平均値を比較すると、テフロンを被せていないとき 2.63 × 104 に対して、テフロンを被せたとき 2.75 × 104 であり、検出できる光子は平均して約 3 個の変化しかみられなかっ25

た。そこで、直接 MPPC の受光面に向かったシンチレーション光子の数を検討するため、
以下のようにモンテカルロシミュレーションを行った。
5mm 角のシンチレーターの中央で光電効果を起こしたと仮定し、この点に対し MPPCの受光面が占める立体角 ω は

立体角 =

受光面積
距離2

で表され、受光面積は 1mm2 , 距離は 2.5mm である。
MPPC の検出効率は、開口率 (ピクセルの境界付近など不感部分をのぞいた割合)61.5 ％と量子効率および電子雪崩形成確率の積で決まり、それを図 6.8 に示す。(検出効率にはクロストークとアフターパルスによる影響が含まれている) これによると、GAGG の発光波長
520[nm] のとき光子検出効率は 40 ％である。[1]図 6.8: MPPC の検出効率 (浜松ホトニクス社のカタログより抜すい)
すると、全立体角が 4π[sr] であるので、受光部に届いて検出される光子の割合は

11ωεdet =
×× 0.40 = 0.0051
24π
2.5 4πと見積もられる。

26

シンチレーター GAGG の発光量は、製造者である古河機械金属の仕様表（表 3.2 参照)によれば 60,000[photon/MeV] である。137 Cs の γ 線エネルギー 0.662MeV がすべてシンチ
レーター内で失われたと仮定すると MPPC で検出される光子の期待値は、

60,000[photon/MeV] × 0.662[MeV] × 0.51[％]=202[photons]と見積もられる。ここで、MPPC は同じピクセルに 2 個以上の光子が入った場合には 1 個の光子が入った時と同じ信号が出ることを考慮し、MPPC によって検出できる光子数のシミュレーションを以下のようなアルゴリズムで行った。
√
検出される光子数の期待値が 202 であるから、個々の事象では ± 202 個の標準偏差を持
つガウス分布するものとした。表 2.1 にあるように本実験で用いた MPPC は 400 のピクセルを持つため、光子は 400 のピクセルに一様分布で入射し、同一ピクセルに複数個の光子が入ったときはそれを１光子とカウントした。
図 6.9 には各ピクセルに入射した光子数の分布、図 6.10 は光子を検出したピクセル数の分布を示す。

27

図 6.9: 各ピクセルに入射した光子数のヒストグラム

図 6.10: 光子を検出できたピクセル数のヒストグラム
これらにより、平均 202 個の光子が MPPC の有感領域に達するが、同一ピクセルに複数個の光子が入る場合があるため、結果として、平均して 159 個程度の光子が検出されるという見積もりになる。これは実際に実験データから得た 76 個と比較すると76 ÷ 159=0.477となり、実験値は期待値の 48 ％程度であるという結論に達する。この約 2 倍の減少が生じている理由は現時点では未解明である。

28

第 7 章 まとめ
LED 点灯を MPPC で受け、波高分布からその増倍率を測定し、105 オーダーの増倍率を示し、フォトンカウンティンぐ能力を持つことを確認した。GAGG シンチレーターと組み合わせて、1 37Cs が発する 662keV の γ 線を検出できることを確認した。ただし、シンチレーター内で 662keV の全てが失われ、発生したシンチレーション光のうち MPPC 受
光部に直接入ったもののみが検出されるという仮定に基づいた期待値と比べ、実験データから得た測定値は半分程度であった。この約 2 倍の減少が生じている理由は現時点では未解明である。

29

謝辞
本研究を進めるにあたり、始終ご指導ご鞭撻を頂きました指導教官の宮林謙吉先生に心より感謝いたします。また、たくさんのアドバイスを頂きました林井久樹先生、下村真弥
先生に深く感謝いたします。そして、たびたびご助言、ご協力いただいた先輩方、1 年間
高エネルギー物理学研究室で共に研究を進めてきた同回生にも大変お世話になりありがとうございました。
たくさんの方々のご支援とご協力により、卒業論文を完成させることができました。皆様へ心から感謝の気持ちと御礼を申し上げたく、謝辞にかえさせていただきます。

30

付 録A

増倍率測定時に用いたデータ収集プログラムのソースコード//header file#include <unistd.h>#include <stdio.h>#include "vxi11_user.h"
int main(int argc , char *argv[]){
FILE *fp;
fp=fopen("Cs03032.txt", "w");
int nevt = 5000;//event number
CLINK *clink;
clink = new CLINK;
static char* serverIP = "10.0.1.108";
int A;//-- open scopeif (vxi11_open_device(serverIP , clink) != 0){
printf ("Couldn ’t open scope.\n");exit(1);}//-- buffer to receive data block.
static char wf1[10000],wf2[10000];char buffer[40];int ret;ret = vxi11_send(clink , "DAT:SOU CH2");//Digitize ch=2 signal.
ret = vxi11_send(clink , "HOR:SCA 40E-09"); //Horizontal 40 ns/div
ret = vxi11_send(clink , "CH2:SCA 10E-03");//ch=2 vertical 10mV/div
ret = vxi11_send(clink , "TRIG:A:EDGE:SOU CH2");//triger type = edge.ch=2
ret = vxi11_send(clink , "TRIG:A:EDGE:SLO FALL");//slope control=fall
ret = vxi11_send(clink , "TRIG:A:LEV:CH2 -37.6E-03");//triger level = -37.6mV
ret = vxi11_send(clink , "TRIG:A:MODE NORM");//trigered data get
ret = vxi11_send(clink , "CH1:POS 3.0"); //vertical position
ret = vxi11_send(clink , "HOR:DEL:MOD OFF"); //delay mode off
ret = vxi11_send(clink , "HOR:POS 20");//trig pos 20% from left
ret = vxi11_send(clink , "HOR:RECO 1000");//record length
ret = vxi11_send(clink , "ACQ:MOD SAM");//acquisition mode = 8bit
double xinc = vxi11_obtain_double_value(clink , "WFMO:XIN?"); //get horizontal range
printf("%10.3e\n",xinc);fprintf(fp,"%10.3e\n",xinc );for (int iev = 0 ; iev < nevt + 1; iev++){
ret = vxi11_send(clink , "DAT:ENC FAS");// fastest encording
ret = vxi11_send(clink , "WFMO:BYT_N 1"); //renge = -128~127
ret = vxi11_send(clink , "DAT:STAR 1");
ret = vxi11_send(clink , "DAT:STOP 1000");
ret = vxi11_send(clink , "CURV?");//get data
long bytes_returned=vxi11_receive_data_block(clink , wf1, 10000, 1000);
if (iev == 0){
for (int j = 0 ; j < 200; j++){
wf2[j]=wf1[j];31continue;}}for (int k = 0; k < 200; k++){A = wf2[k]-wf1[k];if (A==0){
continue;}break;}
printf(" %5d %5d ", iev, bytes_returned);
fprintf(fp," %5d %5d ",iev,bytes_returned);for (int i = 0 ; i < 1000; i++){fprintf(fp,"%d ",(char)wf1[i]);}printf("\n");fprintf(fp,"\n");}fclose(fp);printf("ending... \n");vxi11_close_device(serverIP , clink);
}

32

付 録B

増倍率データ解析プログラム#include <stdio.h>
int main(){
int i, n, j;
float xdiv; /* horizonta axises range. */int iev, nsample; /* Event number and number of sampling. */int siny[10000];
FILE *fp;
FILE *fpout;fp = fopen("Cs03032.txt", "r"); //file open for readif(fp == NULL){
printf("can’t open file \n");return 0;}
fpout = fopen("Cs03032i.txt","w"); //file open for write
if(fp == NULL){
printf("can’t open file \n");return 0;}/* Top line is x/div. */fscanf(fp, "%e\n", &xdiv );printf("xdiv read done.\n");/* Read event number unless EOF. */while(fscanf(fp,"%d",&iev)!=EOF){
printf("event number=%d",iev);/* Check number of samplings. */fscanf(fp, "%d", &nsample );/* Read sample and hold data.*/const int nstot =1000;//taking 1000 sample is normal.
for(i=0; i<nstot; i++){fscanf(fp, "%d", &(siny[i]) );//get data}
printf("\n");if(nsample == nstot){/* Do needed instructions for the read data. */double sum = 0.0;
double sum2 = 0.0;
double sum3 = 0.0;
double dev = 0.0;
double ave = 0.0;
double delta = 0.0;
double max = 0.0;
const int nped = 100;/* sum */for(i=0;i<nped;i++){sum = sum + (double)siny[i];}/* average */ave = sum / nped;/* deviation*/33

for(i=0;i<nped;i++){sum2 = sum2 + (ave - (double)siny[i])*(ave - (double)siny[i]);}dev = sum2 / nped;/* delta */for(i=nped+10;i<300;i++){delta = ave - (double)siny[i];/*pulse hight*/
if(max < delta){
max = delta;}/*integral*/sum3 = sum3 + delta;}/* select data *///fprintf(fpout,"%d %f %f %f %f\n" ,iev,ave,dev,max,sum3);
printf("%d %f %f %f %f\n",iev,ave,dev,max,sum3);
fprintf(fpout,"%f\n" ,sum3); //only integral}}/* Event loop end. */fclose(fp);fclose(fpout);
printf("ending ...\n");
}

34

付 録C

検出光子数のシミュレーション#include<stdio.h>#include<stdlib.h>#include<math.h>//#include<random.h>/*ramdom number of gaus generation*/
int main(int argc, char* argv[]){
int ipix[400];//ピクセル数
int nran[8];FILE *fp;
FILE *fpp;
fp=fopen("np2.txt","w");fpp=fopen("ne2.txt","w");/* expected value */const double mu=202.0;
const double sigma=sqrt(202.0);
srand(10);for(int i=0;i<10000;i++){double r=((double)rand())/((double)RAND_MAX);
double rr=((double)rand())/((double)RAND_MAX);
double z1=sqrt( -2.0*log(r) )*cos( 2.0*M_PI*rr);
double z2=sqrt( -2.0*log(r) )*sin( 2.0*M_PI*rr);
double rand_normal=0.0;
rand_normal= (double)mu +(double)sigma*(double)z1;
printf("event number =%d rand_normal =%f\n",i,rand_normal);
int npix=0;
for(int i=0;i<400;i++){ipix[i]=0.0;}
for(int t=0;t<rand_normal;t++){int iran=rand()%400;
ipix[iran]++;}/* number of photon in pixel
for(int it=0;it<400;it++){if(ipix[it]!=0){
npix++;}}*/

fprintf(fpp,"%d\n",npix);for(int jj=0;jj<400;jj++){fprintf(fp,"%3d",ipix[jj]);}
fprintf(fp,"\n");35}fclose(fp);fclose(fpp);
}

36

参考文献
[1] 浜 松 ホ ト ニ ク ス http://akizukidenshi.com/download/ds/hamamatsu/s1036211 series kapd1022j05.pdf[2] William R.Leo.Techniques for Nuclear and Particle Physics Experiments(1994)p160[3] FURUKAWA シンチレータ結晶物性比較表
http://www.furukawakk.co.jp/pdf/others/hikaku.pdf[4] 長坂憲子、横山紗依 卒業論文 「デジタルオシロスコープによる高速サンプリン
グ記録を用いた MPPC 読み出しシンチレーションカウンターの特性評価 」奈良女子大学（2014 年度）[5] 下雅意美紀、谷川祥子 卒業論文 「シンチレーションカウンターのデジタルオシロ
スコープ読み出し 」奈良女子大学（2013 年度）[6] 瓜生奈都美、福井千尋 卒業論文 「デジタルオシロスコープによるシンチレーショ
ンカウンターの高速サンプリング読み出しの研究」奈良女子大学（2012 年度）[7] 栗林和加、立川晶絵、峰村さつき 卒業論文 「ポジトロニウム消滅ガンマ線対の検出と時間分解能の研究」奈良女子大学（2011 年度）[8] 木原理美、脇田紗弥佳 卒業論文 「MPPC 読み出しによるシンチレーションカウ
ンターを用いた荷電粒子検出」奈良女子大学（2010 年度）[9] 小池 博子、米山 知佐子 卒業論文 「シンチレーションカウンターの MPPC の読み出し」奈良女子大学（2008 年度）[10] 辻 知佳、宮田 香織 卒業論文 「MPPC の基本特性およびシンチレーションカウン
ターへの応用」奈良女子大学（2007 年度）[11] 猪木慶治、川合光/著 「基礎量子力学」 講談社[12] 富永洋、野口正安/著 「放射線応用計測-基礎から実用まで-」 日刊工業新聞社

37