Details | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
25 | f9daq | 1 | { |
2 | TVector3 dum; |
||
3 | gROOT->ProcessLine(".L src/RTUtil.cpp+"); |
||
4 | gROOT->ProcessLine(".L src/guide.cpp+"); |
||
5 | gROOT->ProcessLine(".L src/raySimulator.cpp+"); |
||
6 | gROOT->ProcessLine(".L src/userFunctions.cpp+"); |
||
7 | |||
8 | // Set the parameters of the light guide: |
||
9 | // SiPM0: light guide exit surface linear dimension (surface = SiPM0^2) |
||
10 | // M0: Magnification ratio = b / SiPM |
||
11 | // b: Light guide entry surface linear dimension = M0*SiPM (surface = (M0*SiPM0)^2) |
||
12 | // d0: length of light array |
||
13 | // n10: refractive index outside of the light guide - entry surface and lateral sides |
||
14 | // n20: refractive index of light guide material and the glass plate |
||
15 | // n30: refractive index outside of the light guide at the exit surface (ie optical cement, detector) - it is used only by Fresenl eq, not for the refraction |
||
16 | //gROOT->ProcessLine("SetLG(3.0, 1.666, 3.0, 1.0, 1.48, 1.48)"); |
||
17 | |||
18 | // Set light guide surface types |
||
19 | // |
||
20 | // void SetLGType(int in = 1, int side = 1, int out = 1) |
||
21 | // 0 = dummy; 1 = refractor; 2 = reflector; 3 = total reflection (oboslete, total reflection now included with refractor type) |
||
22 | // |
||
23 | //gROOT->ProcessLine("SetLGType(1, 1, 1)"); |
||
24 | |||
25 | // Gap between light guide exit surface and sensitive area (SiPM) |
||
26 | // supposed it is epoxy of the same refractive index as light guide |
||
27 | gROOT->ProcessLine("SetGap(0.3, 0, 0)"); |
||
28 | |||
29 | // Glass surface of the detector on/off; width |
||
30 | gROOT->ProcessLine("SetGlass(0, 0)"); |
||
31 | |||
32 | // Glass plate above the light guide on/off |
||
33 | gROOT->ProcessLine("SetPlate(1, 1)"); |
||
34 | |||
35 | // Draw the detector and propagate certain number of rays |
||
36 | gROOT->ProcessLine("Show3D(1)"); |
||
37 | |||
38 | // Draw the statistics |
||
39 | gROOT->ProcessLine("ShowData(1)"); |
||
40 | |||
41 | // Propagate 1e4 rays, isotropicaly distributed in 0-30 degree cone (relative to entry surface normal) |
||
42 | // draws first 20 rays (slow) |
||
43 | gROOT->ProcessLine("LGI(10000,30,20)"); |
||
44 | //gROOT->ProcessLine("LGI_ad()"); |
||
45 | //gROOT->ProcessLine("LGS()"); |
||
46 | } |
||
47 |