From 25444667e0ddf08f45ff4fcfec8c4cf4de0b2523 Mon Sep 17 00:00:00 2001
From: Michiel Cottaar <michiel.cottaar@ndcn.ox.ac.uk>
Date: Thu, 1 Feb 2024 14:31:07 +0000
Subject: [PATCH] Switch slice_thickness unit from um to mm

---
 src/overlapping/trapezoid_gradients.jl | 4 ++--
 src/variables.jl                       | 4 ++--
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/overlapping/trapezoid_gradients.jl b/src/overlapping/trapezoid_gradients.jl
index 04b3d04..39fb32e 100644
--- a/src/overlapping/trapezoid_gradients.jl
+++ b/src/overlapping/trapezoid_gradients.jl
@@ -37,7 +37,7 @@ If not set, they will be determined during the sequence optimisation.
 - [`qval`](@ref)/[`qvec`](@ref): Spatial scale on which spins will be dephased due to this pulsed gradient in rad/um (given by `δ` * `gradient_strength`).
 ### Pulse variables
 Any variables defined for the specific pulse added. Also:
-- [`slice_thickness`](@ref): vector with the slice thickness
+- [`slice_thickness`](@ref): vector with the slice thickness in mm.
 
 The [`bvalue`](@ref) can be constrained for multiple gradient pulses.
 """
@@ -130,7 +130,7 @@ slew_rate(g::TrapezoidGradient) = g.slew_rate
 δ(g::TrapezoidGradient) = rise_time(g) + flat_time(g)
 duration(g::TrapezoidGradient) = 2 * rise_time(g) + flat_time(g)
 qvec(g::TrapezoidGradient, ::Nothing, ::Nothing) = δ(g) .* gradient_strength(g) .* 2π
-inverse_slice_thickness(g::TrapezoidGradient) = isnothing(g.pulse) ? nothing : inverse_bandwidth(g) .* gradient_strength(g)
+inverse_slice_thickness(g::TrapezoidGradient) = isnothing(g.pulse) ? nothing : inverse_bandwidth(g) .* gradient_strength(g) .* 1000
 
 function variables(tg::TrapezoidGradient) 
     list = [qvec, δ, gradient_strength, duration, rise_time, flat_time]
diff --git a/src/variables.jl b/src/variables.jl
index 92025d4..7539c5b 100644
--- a/src/variables.jl
+++ b/src/variables.jl
@@ -19,8 +19,8 @@ all_variables_symbols = [
         :inverse_bandwidth => "Inverse of the [`bandwidth`](@ref) of the RF pulse in ms",
         :N_left => "The number of zero crossings of the RF pulse before the main peak",
         :N_right => "The number of zero crossings of the RF pulse after the main peak",
-        :slice_thickness => "Slice thickness of an RF pulse that is active during a gradient.",
-        :inverse_slice_thickness => "Inverse of the [`slice_thickness`](@ref) in 1/um.",
+        :slice_thickness => "Slice thickness of an RF pulse that is active during a gradient in mm.",
+        :inverse_slice_thickness => "Inverse of the [`slice_thickness`](@ref) in 1/mm.",
     ],
 
     # gradients
-- 
GitLab