A simple Photon tracing program

originally made by Sekiguchi

Parameters used in this program are

#define Lx 10.0          /*X length(mm)*/

#define Ly 45.0          /*Y length(mm)*/

#define Lz 3.0           /*Thickness(mm)*/

#define R 0.5            /*Fiber radius(mm)*/

#define Fx 25.0          /*Xposition of fiber*/

#define Fz 1.5           /*Zposition of fiber*/

#define atts 400.0       /*Attenuation length of scintillator (mm)*/

#define attf 3500.0      /*Attenuation length of fiber (mm)*/

#define nf 1.49          /*Refractive index of fiber*/

#define ns 1.50          /*Refractive index of scintillator*/         

#define ef 0.01338       /*Absorptivity of fiber*/

#define cf 200.0         /*density of fiber (??) */

#define Lf 150.0         /*fiber length from PMT*/

#define rrp 0.95         /*refective of 3M radiant mirror film (for up and bottom) */

#define rrp2 0.80       /*reflective of 3M sheet (for side of scinti)*/

 この状態でmappingしたら、ｘ−ｙ面にｚは光子数

あまりにがたがたしてないかい？

（ここでファイバー中の発光はどうなっているか、未チェック）

xy軸へ射影すると

あるいは

次のプロットはファイバー中に入って光センサーに到達した光子数分布

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

#define rrp 0.95         /*refective of 3M radiant mirror film (for up and bottom) */

#define rrp2 0.95      /*reflective of 3M sheet (for side of scinti)*/

2番目のside reflectionを0.8 (上の図）から0.95に変えると 

光量は３０％程度増えた

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

 check  光子数分布のなぞ、

　    aa=Lx*rn1/rm; // aa : souce position in x random

// bb=Ly*rn2/rm; // aa : souce position in y random

    bb=20.0;

yを一定にしてみた：やはり依存性は変だ。つまり一定yから来ているのではない。

ーーーーーーーーーーー

次は一定xである

//    aa=Lx*rn1/rm; // aa : souce position in x random

bb=Ly*rn2/rm; // aa : souce position in y random

    aa=3.0;

としてx=3.0でいろいろyでスキャンすると、つまりファイバーからの距離依存な訳だ。当然と言えば当然の依存性、結局光の直接のファイバーの立体角が一番効いている。

こいつですね。x=3.0　のプロット

ーーーーーーーーーーーーーー

x＝４．５では

同じだ？？いや、よく似ている。

-------------------------

x=0.5もガウス形

ーーーーーーーーーーーーー

x=5.0はどうなっているか？

---------------------------------------

これより前のバージョンでは　ファイバーもシンチレーションを発していた。

with fiber correction where no light emission in the fibre

ref=0.95 for both up/down and sides

----------------------------

x dependence

x=3.0  Np=470

x=4.0 Np=485

----

x=2.0 Np=415

----

x=1.0 Np=350

-----

x=0.5  Np= 290

--------

------------------------

x=2.0

rrp=rrp2=0.95 abov  Np= 415

rrp=rrp2=0.90  Np=345

rrp=rrp2=0.80 Np=285

rrp=rrp2=0.70 Np=250

rrp=rrp2=0.60 Np=225

---------------------