- Timestamp:
- 12/13/11 12:30:03 (5 months ago)
- Files:
-
- 1 modified
-
branches/mcxtrace-1.0/xlib/samples/PowderN.comp (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/mcxtrace-1.0/xlib/samples/PowderN.comp
r3164 r3217 742 742 p_inc=0; /*is turned off till we fix the Compton scattering code*/ 743 743 /* How to handle this one? Transmit (1) / Incoherent (2) / Coherent (3) ? */ 744 printf("xi,yi,zi=%g %g %g\n",x,y,z);745 744 if (type < p_transmit) { 746 745 type = 1; 747 746 l = l_full; /* Passing through, full length */ 748 printf("passing through %g\n",l3);749 747 PROP_DL(l3); 750 748 // } else if (type >= p_transmit && type < (p_transmit + p_inc)) { … … 757 755 PROP_DL(dl+l0); /* Point of scattering */ 758 756 SCATTER; 759 printf("scattering at %g\n",dl+l0);760 757 } else { 761 758 exit(fprintf(stderr,"PowderN %s: DEAD - this shouldn't happen!\n", NAME_CURRENT_COMP)); 762 759 } 763 printf("xf,yf,zf=%g %g %g\n",x,y,z);764 760 765 761 if (type == 3) { /* Make coherent scattering event */ … … 820 816 /* tmp_k = rotate k_out by alpha0 around 'k' (Debye-Scherrer cone) */ 821 817 rotate(tmp_kx,tmp_ky,tmp_kz, kout_x,kout_y,kout_z, alpha0, kx, ky, kz); 822 printf("ki=%g %g %g\n",kx,ky,kz);823 818 kx = tmp_kx; 824 819 ky = tmp_ky; 825 820 kz = tmp_kz; 826 printf("kf=%g %g %g tth=%g\n",kx,ky,kz,theta*RAD2DEG);827 821 /* Since now scattered and new direction given, calculate path to exit */ 828 822 if (line_info.shape == 1) {
