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