Changeset 3232 for trunk

Timestamp:

12/15/11 17:03:38 (5 months ago)

Author:

erkn

Message:

Perfect crystal now works for any hkl

added support for arbitrary hkls. Still only single material diamond
structures, though.
Can now handle hits from below and backwards.
Fix to detector out data print routine to allow values closely spaced.
There was an implicit hard limit at 4 digits precision.

Files:

• trunk/xlib/samples/Perfect_crystal.comp (modified) (8 diffs)

trunk/xlib/samples/Perfect_crystal.comp

@@ -38 +38 @@
 * width [m]                     width of the crystal (along x-axis)
 * Material                      Si, Ge (maybe also GaAs?)       
-<<<<<<< HEAD
-* V [Ã
-^3]                     unit cell volume
-* h,k,l                         indices of reflection (111)(220)(400)
-=======
 * V [Ã
 ^3]                     unit cell volume160.1826
@@ -48 +43 @@
 * k []                          Miller index of reflection
 * l []                          Miller index of reflection
->>>>>>> Perfect crystal now works for any hkl
 * alpha [rad]                   asymmetry angle (alpha=0 for symmetric reflection, ie the Bragg planes are parallel to the crystal surface)
 * R0 []                         Reflectivity. Overrides the computed Darwin reflectivity. Probably only useful for debugging. 
@@ -86 +80 @@
     theta = asin(lambda/(2*d));                 /* kinematical bragg angle (rad) */
     b = sin(theta - alpha)/sin(theta + alpha);  /* asymmetry factor */
-    
-<<<<<<< HEAD
-    *Theta0 = Thetain + alpha;                  /* (rad) angle between Bragg planes and incident ray */ 
-    *Thetah = b*(*Theta0 - theta) + theta;      /* (rad) Angle betweeb Bragg planes and reflected ray */
-=======
-    b = sin(theta + alpha)/sin(theta - alpha);  /* asymmetry factor */
 
     *Theta0 = Thetain - alpha;                  /* (rad) angle between Bragg planes and incident ray */ 
@@ -100 +88 @@
       *DeltaTheta0 = -1; /*to mark it irrelevant*/
     }
->>>>>>> Perfect crystal now works for any hkl
     
     /* Define polarization factor: */
@@ -248 +235 @@
         PROP_DL(dist);                  /* now the photon is on the mirror surface, ready to be reflected... */ 
         SCATTER;
-<<<<<<< HEAD
-
-        Thetain = atan(kyu/kzu);        /* (rad )The angle of incidence*/
-        
-=======
         /*Check which quadrant the k vector is in to determine sense of alpha. This to allow for hitting the crystal from behind.*/
         int quadrant;
         Thetain=fabs(atan(ky/kz));
->>>>>>> Perfect crystal now works for any hkl
         /* for a start, just define the value of atomic form factors for Si(111) at E = 8.04778keV */
         double d=cbrt(V)/(sqrt(h*h+k*k+l*l));/*this is valid only for cubic structures*/
@@ -264 +245 @@
         fpp = Table_Value(prms.m_t,E,2);
 //        printf("%g %g %g %g\n",f00,f0h,fp,fpp);    
-<<<<<<< HEAD
-/*      printf("Thetain for this ray is: %g: \n",Thetain);*/
-        DarwinReflectivity(&R, &Thetah, &Theta0, &DeltaTheta0, f00, f0h, fp, fpp, V, alpha, h, k, l, M, E, Thetain,pol);
-        Thetaout = Thetah + alpha;      /* (rad) the angle between the crystal surface and the reflected ray */
-/*      printf("Thetaout for this ray is: %g: \n", Thetaout);*/
-
-=======
         //printf("Thetain for this ray is: %g: \n",Thetain);
         if (ky<0 && kz>0){
@@ -301 +275 @@
           kz=kz-2*vnz;
         }while(0);
->>>>>>> Perfect crystal now works for any hkl
         /* The direction of the reflected ray: */
        // printf("kout: %g %g %g\n",kx,ky,kz);
 
-<<<<<<< HEAD
-        NORM(kxout,kyout,kzout); 
-        kx=K*kxout;
-        ky=K*kyout;
-        kz=K*kzout;
-                  
-        /* apply reflectivity */
-        p*=R; 
-
-=======
         /* apply Darwin reflectivity if not is supplied from outside*/
         if (!R0){
@@ -323 +286 @@
         /*catch dead rays*/
         if (p==0) ABSORB;
->>>>>>> Perfect crystal now works for any hkl
-        } else {
-          RESTORE_XRAY(INDEX_CURRENT_COMP, x, y, z, kx, ky, kz, phi, t, Ex, Ey, Ez, p);
-        } 
+    } else {
+      RESTORE_XRAY(INDEX_CURRENT_COMP, x, y, z, kx, ky, kz, phi, t, Ex, Ey, Ez, p);
+    } 
   }