31:Tdc8,viewof params = { // create individual inputs const energy = Inputs.range([60, 300], { value: 300, step: 10, label: "energy [kV]" }); const semiangle = Inputs.range([5, 30], { value: 20, step: 1, label: "angle [mrad]" }); const C10 = Inputs.range([-500, 500], { value: 100, step: 1, label: "C10 [Å]" }); const C30 = Inputs.range([-1000, 1000], { value: 0, step: 1, label: "C30 [um]" }); const stig = xySlider({ min: -100, max: 100, value: [0, 0], width: "125px", height: "125px" }); const coma = xySlider({ min: -10000, max: 10000, value: [0, 0], width: "125px", height: "125px" }); const resetPosition = Inputs.button("reset position", { value: null, reduce: () => { mutable shift = [0, 0]; } }); const resetAberrations = Inputs.button("reset aberrations", { value: null, reduce: () => { C10.value = 100; C10.dispatchEvent(new Event("input", { bubbles: true })); C30.value = 0; C30.dispatchEvent(new Event("input", { bubbles: true })); stig.value = [0, 0]; stig.dispatchEvent(new Event("input", { bubbles: true })); coma.value = [0, 0]; coma.dispatchEvent(new Event("input", { bubbles: true })); } }); // magnitude labels for 2d sliders const stigLabel = htl.html`

C12: 0.0 Å

`; const comaLabel = htl.html`

C21: 0.0 Å

`; // 2. Build the responsive Grid const form = htl.html`

${energy} ${semiangle} ${C10} ${C30}

${resetPosition} ${resetAberrations}

${stigLabel} ${stig}

${comaLabel} ${coma}

`; const update = () => { const [sx, sy] = stig.value; const [cx, cy] = coma.value; const stigMag = Math.sqrt(sx ** 2 + sy ** 2).toFixed(1); const comaMag = Math.sqrt(cx ** 2 + cy ** 2).toFixed(1); stigLabel.textContent = `C12: ${stigMag} Å`; comaLabel.textContent = `C21: ${comaMag} Å`; form.value = { energy: energy.value, semiangle: semiangle.value, C10: C10.value, C30: C30.value, stig: stig.value, coma: coma.value }; form.dispatchEvent(new Event("input", { bubbles: true })); }; energy.addEventListener("input", update); semiangle.addEventListener("input", update); C10.addEventListener("input", update); C30.addEventListener("input", update); stig.addEventListener("input", update); coma.addEventListener("input", update); update(); return form; }32:[] 2f:{"title":"STEM Probes PSF / CTF","cells":[{"id":0,"value":"# STEM Probes PSF / CTF","pinned":false,"hidden":false,"mode":"md"},{"id":338,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":157,"value":"visualization_output = hstack(\n { gap: gap, minWidth: 272.5, maxWidth: half_panel_width },\n show2d({\n data: probeIntensity,\n width: gpts[0],\n height: gpts[1],\n plotWidth: Math.max(half_panel_width, 272.5),\n plotHeight: Math.max(half_panel_width, 272.5),\n label: \"Real-space probe intensity\",\n showTooltip: false,\n scalebar: { sampling: sampling[1], units: \"Å\", loc: \"lower right\" },\n onHover: ({ ix, iy }) => {\n mutable shift = [iy - gpts[0] / 2, ix - gpts[1] / 2];\n }\n }),\n show2d({\n data: probeCentered,\n width: gpts[0],\n height: gpts[1],\n plotWidth: Math.max(half_panel_width, 272.5),\n plotHeight: Math.max(half_panel_width, 272.5),\n complex: true,\n showTooltip: false,\n scalebar: {\n sampling: angularSampling[1],\n units: \"mrad\",\n loc: \"lower right\"\n },\n label: \"Fourier-space complex probe\"\n })\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":307,"value":"---\n### Inputs","pinned":false,"hidden":false,"mode":"md"},{"id":437,"value":"gap = 4","pinned":false,"hidden":false,"mode":"ojs"},{"id":441,"value":"half_panel_width = (viz_width - gap) / 2","pinned":false,"hidden":false,"mode":"ojs"},{"id":409,"value":"// viewof viz_width = Inputs.range([400, 800], { value: 680 })\nviz_width = 677.5","pinned":false,"hidden":false,"mode":"ojs"},{"id":11,"value":"gpts = [256, 256]","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"sampling = [0.2, 0.2]","pinned":false,"hidden":false,"mode":"ojs"},{"id":53,"value":"aberrationCoefs = ({\n C10: params.C10,\n C12: Math.sqrt(params.stig[0] ** 2 + params.stig[1] ** 2),\n phi12: Math.atan2(params.stig[1], params.stig[0]) / 2,\n C21: Math.sqrt(params.coma[0] ** 2 + params.coma[1] ** 2),\n phi21: Math.atan2(params.coma[1], params.coma[0]) / 1,\n C30: params.C30 * 1e4\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":8,"value":"wavelength = electronWavelength(params.energy * 1e3)","pinned":false,"hidden":false,"mode":"ojs"},{"id":253,"value":"mutable shift = [0, 0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":318,"value":"---\n### Imports / Functions","pinned":false,"hidden":false,"mode":"md"},{"id":175,"value":"import {\n electronWavelength,\n spatialFrequencies,\n polarCoordinates,\n scatteringAngles,\n probeFourierSpace,\n createFFTUtils,\n fftshift2d,\n complexToIntensity,\n show2d,\n hstack\n} from \"@gvarnavi/em-utilities\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":344,"value":"import { xySlider } from \"@mootari/2d-slider\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":15,"value":"webfft = import(\"https://cdn.skypack.dev/webfft@1.0.3?min\").then(\n (m) => m.default\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":37,"value":"fftInstance = {\n const instance = new webfft(gpts[1]);\n instance.profile();\n // invalidation.then(() => instance.dispose());\n return instance;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":177,"value":"fft = createFFTUtils(fftInstance)","pinned":false,"hidden":false,"mode":"ojs"},{"id":57,"value":"cartesianFrequencies = spatialFrequencies(gpts, sampling)","pinned":false,"hidden":false,"mode":"ojs"},{"id":204,"value":"polarFrequencies = polarCoordinates(\n cartesianFrequencies.kxa,\n cartesianFrequencies.kya\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":208,"value":"angularFrequencies = scatteringAngles(\n polarFrequencies.k,\n polarFrequencies.phi,\n wavelength\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":63,"value":"angularSampling = [\n (1 / sampling[0] / gpts[0]) * wavelength * 1e3,\n (1 / sampling[1] / gpts[1]) * wavelength * 1e3\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":224,"value":"probeFS = probeFourierSpace(\n angularFrequencies.alpha,\n angularFrequencies.phi,\n wavelength,\n params.semiangle,\n angularSampling,\n aberrationCoefs\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":282,"value":"probeFS_shifted = fft.fourierShift(\n probeFS,\n gpts[0],\n gpts[1],\n shift[0],\n shift[1],\n { inputFS: true, outputFS: true }\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":82,"value":"probeReal = {\n const [nx, ny] = gpts;\n const N = nx * ny;\n const result = fft.ifft2d(probeFS_shifted, nx, ny);\n // normalize\n for (let i = 0; i < 2 * N; i++) result[i] /= N;\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":88,"value":"probeIntensity = complexToIntensity(fftshift2d(probeReal, gpts[0], gpts[1]))","pinned":false,"hidden":false,"mode":"ojs"},{"id":111,"value":"probeCentered = fftshift2d(probeFS_shifted, gpts[0], gpts[1])","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W32","stdlib":"1","titling":"auto"}