6.1.3. Example of masks

6.1.3.1. Creating an instance

[1]:

from diffractio import np, plt
from diffractio import degrees, mm, um
from diffractio.scalar_masks_X import Scalar_mask_X

from numpy import loadtxt

6.1.3.2. Procedures to generate masks

6.1.3.2.1. mask_from_function

[2]:

num_data = 1024
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

f1 = "R1-h1+np.sqrt(R1**2-(self.x-x0)**2)"
f2 = "R2-h2+np.sqrt(R2**2-(self.x-x0)**2)"
v_globals = {
    'R1': 5 * mm,
    'R2': 1 * mm,
    'x0': 0 * um,
    'y0': 0 * um,
    'h1': 1 * mm,
    'h2': -1 * mm
}
index = 1.5
print(v_globals)

t1 = Scalar_mask_X(x, wavelength)
t1.mask_from_function(x0=0 * um,
                      index=index,
                      f1=f1,
                      f2=f2,
                      v_globals=v_globals,
                      radius=100 * um)
t1.draw(kind='phase')

{'R1': 5000.0, 'R2': 1000.0, 'x0': 0.0, 'y0': 0.0, 'h1': 1000.0, 'h2': -1000.0}

../../../_images/source_tutorial_scalar_X_masks_x_5_1.png

6.1.3.2.2. mask_from_array

[3]:

num_data = 1024 * 8
x = np.linspace(-1 * mm, 1 * mm, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)

profile1 = loadtxt('profile1.txt')
profile2 = loadtxt('profile2.txt')
t1.mask_from_array(x0=0,
                   index=1.25,
                   array1=profile1 * mm,
                   array2=profile2 * mm,
                   interp_kind='quadratic',
                   radius=1.5 * mm)

t1.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_7_0.png

6.1.3.3. slit

[4]:

num_data = 128
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.slit(x0=0, size=100 * um)
t1.draw()

../../../_images/source_tutorial_scalar_X_masks_x_9_0.png

6.1.3.4. double_slit

[5]:

num_data = 128
length = 25 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.double_slit(x0=0, size=5 * um, separation=15 * um)
t1.draw()

../../../_images/source_tutorial_scalar_X_masks_x_11_0.png

6.1.3.5. two levels

[6]:

num_data = 128
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.two_levels(level1=0, level2=.5, x_edge=0)
t1.draw()

../../../_images/source_tutorial_scalar_X_masks_x_13_0.png

6.1.3.6. gray scale

[7]:

x = np.linspace(0, 250 * um, 1024)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.gray_scale(num_levels=16, levelMin=2, levelMax=3)
t1.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_15_0.png

6.1.3.7. prism[9]:
num_data = 1024
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.prism(x0=0, n=1.5, angle=1 * degrees)
t1.draw(kind='phase')

6.1.3.8. Fresnel biprism

[10]:

num_data = 1024
length = 500 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.biprism_fresnel(angle=1 * degrees, x0=0 * um, radius=125)
t1.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_19_0.png

6.1.3.9. Fresnel biprism from material properties n and h

[11]:

num_data = 1024
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.biprism_fresnel_nh(x0=0, width=100 * um, height=5 * um, n=1.5)

t1.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_21_0.png

6.1.3.10. standard lens

[12]:

num_data = 2048
length = 500 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.lens(x0=0 * um, radius=250 * um, focal=5 * mm)
t1.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_23_0.png

6.1.3.11. Fresnel lens: amplitude and phase

[13]:

num_data = 1024
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

[14]:

t1 = Scalar_mask_X(x, wavelength)
t1.fresnel_lens(x0=0 * um,
                radius=100 * um,
                focal=.25 * mm,
                kind='amplitude',
                phase=np.pi)
t1.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_26_0.png

[15]:

t1.fresnel_lens(x0=0 * um,
                radius=100 * um,
                focal=.25 * mm,
                kind='phase',
                phase=np.pi)
t1.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_27_0.png

6.1.3.12. rough surface

[16]:

num_data = 128
length = 250 * um
x = np.linspace(-length / 2, length / 2, num_data)
wavelength = 0.6328 * um

t1 = Scalar_mask_X(x, wavelength)
t1.roughness(t=15 * um, s=.2 * um)
t1.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_29_0.png

6.1.3.13. Dust particles with different sizas

TODO: does not work properly

[17]:

wavelength = 0.6328 * um
x = np.linspace(-500 * um, 500 * um, 2048)

t1 = Scalar_mask_X(x, wavelength)
t1.dust_different_sizes(percentage=0.2, size=20 * um, std=5 * um)
t1.draw()

../../../_images/source_tutorial_scalar_X_masks_x_31_0.png

6.1.3.14. dust particles

[18]:

wavelength = 0.6328 * um
x = np.linspace(-500 * um, 500 * um, 2048 * 8)

t1 = Scalar_mask_X(x, wavelength)
t1.dust(percentage=0.9, size=20 * um)
t1.draw()

../../../_images/source_tutorial_scalar_X_masks_x_33_0.png

6.1.3.15. bit structure

[19]:

num_data = 2048 * 4

wavelength = 0.85 * um
code = [1, 1, 0, 0, 1, 0, 1, 1, 0, 1]
bit_width = 80 * um
x = np.linspace(0, bit_width * len(code), num_data)

t1 = Scalar_mask_X(x, wavelength)
t1.binary_code(kind='normal', code=code, bit_width=bit_width, x0=0 * um)
t1.draw()

../../../_images/source_tutorial_scalar_X_masks_x_35_0.png

6.1.3.16. Gratings

6.1.3.16.1. sine grating

[20]:

x = np.linspace(-250 * um, 250 * um, 512)
wavelength = 0.6238 * um
period = 40 * um
red = Scalar_mask_X(x, wavelength)
red.sine_grating(period=period, amp_min=0, amp_max=1, x0=0 * um)
red.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_37_0.png

6.1.3.16.2. Ronchi grating

[21]:

x = np.linspace(-250 * um, 250 * um, 512)
wavelength = 0.6238 * um
red = Scalar_mask_X(x, wavelength)
red.ronchi_grating(period=50 * um, x0=0 * um, fill_factor=0.5)
red.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_39_0.png

6.1.3.16.3. binary grating: amplitude and phase

[22]:

x = np.linspace(-250 * um, 250 * um, 512)
wavelength = 0.6238 * um
red = Scalar_mask_X(x, wavelength)
red.binary_grating(period=50 * um,
                   amin=0.25,
                   amax=0.5,
                   phase=np.pi,
                   x0=25 * um,
                   fill_factor=0.55)
red.draw(kind='amplitude')
plt.ylim(0, 1)

../../../_images/source_tutorial_scalar_X_masks_x_41_0.png

[23]:

x = np.linspace(-250 * um, 250 * um, 512)
wavelength = 0.6238 * um
period = 40 * um
red = Scalar_mask_X(x, wavelength)
red.binary_grating(period=period,
                   amin=1,
                   amax=1,
                   phase=np.pi / 2,
                   x0=0,
                   fill_factor=0.5)
red.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_42_0.png

6.1.3.16.4. blazed grating

[24]:

x = np.linspace(-250 * um, 250 * um, 512)
wavelength = 0.6238 * um
period = 50 * um
red = Scalar_mask_X(x, wavelength)
red.blazed_grating(x0=period/2, period=period, phase_max=2*np.pi)
red.draw(kind='phase')

../../../_images/source_tutorial_scalar_X_masks_x_44_0.png

6.1.3.16.5. chirped grating: general

[25]:

x = np.linspace(-250 * um, 250 * um, 512)
wavelength = 0.6238 * um

fx = '10+20*(self.x/length)**2'

red = Scalar_mask_X(x, wavelength)
red.chirped_grating(kind='amplitude_binary',
                    p_x=fx,
                    length=500 * um,
                    x0=0 * um,
                    amp_max=1,
                    amp_min=0,
                    delta_x=0,
                    phase_max=np.pi)
red.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_46_0.png

6.1.3.16.6. chirped grating: linear p

[26]:

x = np.linspace(0 * um, 100 * um, 4096 * 4)
wavelength = 0.6238 * um
red = Scalar_mask_X(x, wavelength)
conds = {
    'kind': 'amplitude_binary',
    'p0': 20 * um,
    'p1': 4 * um,
    'amp_min': 0,
    'amp_max': 1,
    'delta_x': 0,
    'phase_max': np.pi,
    'length': 0,  # equal to length of x
}
red.chirped_grating_p(**conds)
red.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_48_0.png

6.1.3.16.7. chirped grating: linear q

[27]:

x = np.linspace(-100 * um, 100 * um, 4096 * 4)
wavelength = 0.6238 * um
red = Scalar_mask_X(x, wavelength)
conds = {
    'kind': 'amplitude_binary',
    'p0': 20 * um,
    'p1': 4 * um,
    'amp_min': 0,
    'amp_max': 1,
    'delta_x': 0,
    'phase_max': np.pi,
    'length': 0,
}
red.chirped_grating_q(**conds)
red.draw(kind='amplitude')

../../../_images/source_tutorial_scalar_X_masks_x_50_0.png