Subversion Repositories f9daq

Rev

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