Skip to content

Commit

Permalink
Initial footwork toward #1736
Browse files Browse the repository at this point in the history
  • Loading branch information
willend committed Oct 11, 2024
1 parent 1004777 commit 9865c3f
Showing 1 changed file with 120 additions and 0 deletions.
120 changes: 120 additions & 0 deletions mcxtrace-comps/optics/Collimator_linear.comp
Original file line number Diff line number Diff line change
@@ -0,0 +1,120 @@
/*******************************************************************************
*
* McXtrace, x-ray tracing package
* Copyright, All rights reserved
* DTU Physics, Kgs. Lyngby, Denmark
* Synchrotron SOLEIL, Saint-Aubin, France
*
* Component: Collimator_linear
*
* %I
* Written by: Peter Willendrup
* Date: October 2024
* Origin: DTU
* Modified by: Peter Willendrup from McStas Collimator_linear from Krisitan Nielsen
*
* A simple analytical Soller collimator (with triangular transmission).
*
* %D
* Soller collimator with rectangular opening and specified length. The
* transmission function is an average and does not utilize knowledge of the
* actual neutron trajectory. A zero divergence disables collimation (then the
* component works as a double slit).
*
* Example: Collimator_linear(xmin=-0.1, xmax=0.1, ymin=-0.1, ymax=0.1, length=0.25, divergence=40,transmission=0.7)
*
* %P
* INPUT PARAMETERS:
*
* xmin: [m] Lower x bound on slits
* xmax: [m] Upper x bound on slits
* ymin: [m] Lower y bound on slits
* ymax: [m] Upper y bound on slits
* xwidth: [m] Width of slits
* yheight: [m] Height of slits
* length: [m] Distance between input and output slits
* divergence: [minutes of arc] Divergence horizontal angle (calculated as atan(d/length), where d is the blade spacing)
* transmission: [1] Transmission of Soller (0<=t<=1)
* divergenceV: [minutes of arc] Divergence vertical angle
*
* %E
*******************************************************************************/
DEFINE COMPONENT Collimator_linear
DEFINITION PARAMETERS ()
SETTING PARAMETERS (xmin=-0.01, xmax=0.01, ymin=-0.025, ymax=0.025, xwidth=0, yheight=0, length=0.1, divergence=10, transmission=1, divergenceV=0)
OUTPUT PARAMETERS ()
/* X-ray parameters: (x,y,z,kx,ky,kz,phi,t,Ex,Ey,Ez,p) */

DECLARE
%{
double slope;
double slopeV;
%}
INITIALIZE
%{
slope = tan(MIN2RAD*divergence);
slopeV= tan(MIN2RAD*divergenceV);
if (xwidth > 0) { xmax = xwidth/2; xmin = -xmax; }
if (yheight > 0) { ymax = yheight/2; ymin = -ymax; }

if ((xmin >= xmax) || (ymin >= ymax)) {
printf("Collimator_linear: %s: Null slit opening area !\n"
"ERROR (xwidth,yheight,xmin,xmax,ymin,ymax). Exiting",
NAME_CURRENT_COMP);
exit(0);
}
%}
TRACE
%{
double Phi, dl;

PROP_Z0;
if (x<xmin || x>xmax || y<ymin || y>ymax)
ABSORB;

double dirx,diry,dirz;

/* Use direction vector z-component to decide needed propagation length */
dirx=kx;diry=ky;dirz=kz;
NORM(dirx,diry,dirz);

PROP_DL(length/dirz);
if (x<xmin || x>xmax || y<ymin || y>ymax)
ABSORB;

if(slope > 0.0)
{
Phi = fabs(kx/kz);
if (Phi > slope)
ABSORB;
else
p *= transmission*(1.0 - Phi/slope);
SCATTER;
}
if (slopeV > 0) {
Phi = fabs(ky/kz);
if (Phi > slopeV)
ABSORB;
else
p *= transmission*(1.0 - Phi/slopeV);
SCATTER;
}
%}

MCDISPLAY
%{
double x;
int i;


for(x = xmin, i = 0; i <= 10; i++, x += (xmax - xmin)/10.0)
multiline(5, x, (double)ymin, 0.0, x, (double)ymax, 0.0,
x, (double)ymax, (double)length, x, (double)ymin, (double)length,
x, (double)ymin, 0.0);
line(xmin, ymin, 0, xmax, ymin, 0);
line(xmin, ymax, 0, xmax, ymax, 0);
line(xmin, ymin, length, xmax, ymin, length);
line(xmin, ymax, length, xmax, ymax, length);
%}

END

0 comments on commit 9865c3f

Please sign in to comment.