Subversion Repositories f9daq

Rev

Blame | Last modification | View Log | RSS feed

  1. // own include
  2. #include "include/RTUtil.h"
  3. #include "include/guide.h"
  4.  
  5. // general include
  6. #include "TROOT.h"
  7. #include "TSystem.h"
  8. #include "TStyle.h"
  9. #include "TCanvas.h"
  10. #include "TMath.h"
  11. #include "TView3D.h"
  12. #include "TH1F.h"
  13. #include "TH2F.h"
  14. #include "TBenchmark.h"
  15. #include "TPolyMarker.h"
  16. #include "TGraph.h"
  17. #include "TF1.h"
  18.  
  19. //=================================================================================
  20. TCanvas *c3dview;
  21. TCanvas *cacc;
  22. int show_3d = 0, show_data = 0;
  23. int draw_width = 2;
  24.  
  25.  // calls default constructor CVodnik.cpp
  26.  
  27. //void SetCenter(double x = -2.0, double y = 0.0, double z = 0.0)
  28.         //{center.SetXYZ(x,y,z); detector = new CDetector(center);}
  29.  
  30.  
  31. //void SetLGType(int in = 1, int side = 1, int out = 0)
  32.         //{detector->SetLGType(in, side, out);}
  33.  
  34. //void SetLG(double SiPM0 = 3.0, double M0 = 1.666, double d0 = 3.0, double n10 = 1.0, double n20 = 1.48, double n30 = 1.48, double R0 = 0.96)
  35.         //{detector->SetLG(SiPM0, M0, d0, n10, n20, n30, R0);}
  36.        
  37. //void SetR(double R0) {detector->SetR(R0);}
  38.         /*
  39. void SetGlass(int glass_on0 = 0, double glass_d0 = 0.3)
  40.         {detector->SetGlass(glass_on0, glass_d0);}
  41.  
  42. void SetGap(double x_gap0 = 0.3, double y_gap0 = 0.0, double z_gap0 = 0.0)
  43.         {detector->SetGap(x_gap0 , y_gap0, z_gap0);}
  44.        
  45. void SetRCol(int in = 2, int lg = 8, int out = 4, int gla = 6)
  46.         {detector->SetRCol(in, lg, out, gla);};
  47. void SetDCol(int LG0 = 1, int glass0 = 2, int active0 = 3)
  48.         {detector->SetDCol(LG0, glass0, active0);};
  49.        
  50.  
  51. void SetDWidth(int w = 1) {draw_width = w;}
  52.  
  53. void SetFresnel(int b = 1) {detector->SetFresnel(b);}
  54. void SetAbsorb(int b = 0, double A0 = 1e4) {detector->SetAbsorption(b, A0);}
  55.  
  56. void SetGuideOn(int b = 0) {detector->SetGuideOn(b);}
  57.  
  58. void SetTypeU() {SetGlass(1, 0.3);SetGap(0.7);}
  59. void SetTypeC() {SetGlass(0, 0.3);SetGap(1.03);}
  60. */
  61. //-----------------------------------------------------------------------------
  62. void Init(double);
  63. //-----------------------------------------------------------------------------
  64. void DrawData(CDetector*, DetectorParameters&, double, double, int);
  65. //-----------------------------------------------------------------------------
  66. void DrawAcc(int, double*, double*, char*,
  67.                          double, double, double, double);
  68. //-----------------------------------------------------------------------------
  69. void PrintGlassStat(CDetector*);
  70. //-----------------------------------------------------------------------------
  71. void PrintGuideHead();
  72. //-----------------------------------------------------------------------------
  73. void PrintGuideStat(CDetector*, double);
  74. //-----------------------------------------------------------------------------
  75.        
  76. //-----------------------------------------------------------------------------
  77. // en zarek
  78. double Single(CDetector*, DetectorParameters&, TVector3, double, double);
  79. //-----------------------------------------------------------------------------
  80. // zarki, razporejeni v mrezi
  81. double Grid(CDetector*, DetectorParameters&, int, double);
  82. //-----------------------------------------------------------------------------
  83. // zarki z nakljucnim polozajem vpada (na vstopni pov. vodnika)
  84. // vsi pod kotom (theta, phi)
  85. double RandYZ(CDetector*, DetectorParameters&, int, double, double, int);
  86. //-----------------------------------------------------------------------------
  87. // zarki, izotropno porazdeljeni znotraj kota theta
  88. // = nakljucni vstopni polozaj in kot
  89. double RandIso(CDetector*, DetectorParameters&, int, double, int, int);
  90.  
  91.