Subversion Repositories f9daq

Rev

Blame | Last modification | View Log | RSS feed

  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.  
  55.