WebSVN – f9daq – Diff – //lightguide/trunk/include/guide.h

 class CRay
+{
 public:
   CRay() :
     r(TVector3(0,0,0)),
-    n(TVector3(1,0,0)),
+    k(TVector3(1,0,0)),
     //p(TVector3(0,1,0)),
     p(TVector3(0,0,1)),
     color(kBlack)
   {};
-  CRay(TVector3 r0, TVector3 n0) :
+  CRay(TVector3 r0, TVector3 k0) :
     r(r0),
-    n(n0.Unit()),
+    k(k0.Unit()),
     //p(TVector3(0,1,0)),
     p(TVector3(0,0,1)),
     color(kBlack)
   {};
   CRay(double x0, double y0, double z0, double l0, double m0, double n0) :
     r(TVector3(x0,y0,z0)),
-    n(TVector3(l0,m0,n0).Unit()),
+    k(TVector3(l0,m0,n0).Unit()),
     //p(TVector3(0,1,0)),
     p(TVector3(0,0,1)),
     color(kBlack)
   {};
-  void Set(TVector3 r0, TVector3 n0);
+  void Set(TVector3 r0, TVector3 k0);
   //void Set(double x0, double y0, double z0, double l0, double m0, double n0);
   void SetColor(int c){color = c;};
-  void SetPolarization(TVector3 p0) {p = p0.Unit();};
+  void setPolarization(TVector3 p0) {
+    p = p0.Unit();
+    if (p.Dot(k)>1e-3) printf("*** ERROR in CRay: E has component || with k\n");
+  };
   //inline CRay & operator = (const CRay &);
   TVector3 GetR() const {return r;};
-  TVector3 GetN() const {return n;};
+  TVector3 GetK() const {return k;};
   TVector3 GetP() const {return p;};
   void Print();
   void Draw();
   void Draw(double x_from, double x_to);
   void DrawS(double x_from, double t);
+  //r = position, k = unit wave vector, p = polarization
 private:
   TVector3 r;
-  TVector3 n;
+  TVector3 k;
-  TVector3 p; //r = point on line, n = normal, p = polarization
+  TVector3 p;
   int color;
 };
-//=================================================================================
 //=================================================================================
 // ravnina, definirana in omejena s 4 tockami na njej (stirikotni poligon)
 // zanasam se na taksno definicijo: (normala kaze noter)
 private:
   TVector3 r[4], n;
   double A, B, C, D;
   TVector3 edge[4]; // vektorji stranic
   double angle_r[4]; // koti ob posameznem vogalu
 };
 //=================================================================================
 // ravnina - krog
 class CPlaneR
+{
   void Set(TVector3 c, TVector3 n0, double R0)
   {center = c; n = n0; _r = R0;};
   int TestIntersection(TVector3 *vec, CRay in);
   void Draw(int color = 1, int NN = 32);
 private:
   TVector3 n, center;
   double _r;
 };
-//=================================================================================
 //=================================================================================
 // ravna opticna povrsina: refractor, zrcalo ali povrsina s totalnim odbojem
 #define SURF_DUMMY 0
 #define SURF_REFRA 1
 #define SURF_REFLE 2
   {type = type0; CPlane4::Set(r1, r2, r3, r4); SetIndex(n10, n20); reflection = reflectivity;};
   int PropagateRay(CRay in, CRay *out, TVector3 *intersection);
   double N1_N2(int sign) { return ((sign > 0) ? n1/n2 : n2/n1);};
   void Print(){printf("Type = %d\n", type); CPlane4::Print();};
 private:
   int type; //0 = dummy; 1 = refractor; 2 = reflector; 3 = total reflection
   double n1, n2, n1_n2; //index of refraction, n1 @ +normal, n2 @ -normal
   double reflection; // odbojnost stranic
   TRandom rand; // za racunanje verjetnosti odboja od zrcala
   double cosTtotal; //cosinus mejnega kota totalnega odboja za dana n1 in n2
   int fresnel; // ali naj uposteva Fresnelove enacbe
 };
-//=================================================================================
 //=================================================================================
 //    +-------   n1
 //    |       -------
   int getPlateOn() { return _plateOn; };
   double getOffsetY() { return offsetY; };
   double getOffsetZ() { return offsetZ; };
   // bad coupling in the case the small amount of grease was applied
   bool badCoupling() { return _coupling; };
 private:
   double _a;
   double _b;
   double _d;
   double _active;
   double _n1;
   TH2F* GetHOut() const {return hout;};
   int GetExitHits() {return (int)hfate->GetBinContent(5);};
   int GetEnteranceHits() {return (int)hfate->GetBinContent(6);};
   void GetVFate(int *out);
   int GetMAXODB() const {return MAX_REFLECTIONS;};
   void Draw(int color = 1, int width = 1);
   void DrawSkel(int color = 1, int width = 1);
 private:
   Fate fate;
   double _d; // parameters
   double _n1, _n2, _n3; // refractive index: n1 above entry surface, n2 inside, n3 after exit
   double _r;
   TVector3 center, vodnik_edge[8];
   TH1F *hfate, *hnodb_all, *hnodb_exit;
   TH2F *hin, *hout;
 };
-//=================================================================================
-//=============================================================================================================================== <<<<<<<<
 class Plate
+{
 public:
 private:
   TVector3 plate_edge[8];
   CSurface *sides[6];
 };
-// ================================================================================
 class CDetector
+{
 public:
   CDetector(TVector3 center0, DetectorParameters& parameters);
   double offsetZ;
 };
-//=================================================================================
 #endif

Subversion Repositories f9daq

(root)//lightguide/trunk/include/guide.h – Rev 72 → 73