Turn-off of truncated Gaussian via standard deviations (closes #11) (#12

)

* Turn-off of truncated Gaussian via standard deviations (closes #11)

* Version 0.1.1

Co-authored-by: Stefanos Carlström <[email protected]>

Project.toml

@@ -1,7 +1,7 @@
 name = "ElectricFields"
 uuid = "2f84ce32-9bb1-11e8-0d9f-3dce90a4beca"
 authors = ["Stefanos Carlström <[email protected]>"]
-version = "0.1.0"
+version = "0.1.1"
 
 [deps]
 FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"

src/envelopes.jl

@@ -250,10 +250,15 @@ show(io::IO, env::TruncatedGaussianEnvelope) =
 function TruncatedGaussianEnvelope(field_params::Dict{Symbol,Any}, carrier)
     test_field_parameters(field_params, [:T]) # Period time required to round time window up
     test_field_parameters(field_params, [:τ, :σ])
-    test_field_parameters(field_params, [:toff])
+    test_field_parameters(field_params, [:toff, :σoff])
     test_field_parameters(field_params, [:σmax, :tmax, :Tmax])
 
     gaussian_common!(field_params, carrier, Tmax_rounder=NoUnits)
+    @namespace!(field_params) do
+        if σoff
+            toff = σ*σoff
+        end
+    end
 
     @unpack τ, σ, α, toff, tmax, Tmax = field_params
     toff < tmax || throw(ArgumentError("Hard turn-off must occur before end of pulse"))

test/gaussian_pulses.jl

@@ -99,4 +99,19 @@ import ElectricFields: time_integral
 
         @test time_integral(F) ≈ σ*√(2π)
     end
+
+    @testset "Truncated Gaussian turn-off via standard deviations" begin
+        @field(F) do
+            I₀ = 1.0
+            T = 2.0
+            τ = 3.0
+            σoff = 4.0
+            σmax = 5.0
+            env = :trunc_gauss
+        end
+
+        σ = params(F)[:σ]
+        @test params(F)[:toff] ≈ 4σ
+        @test params(F)[:tmax] ≈ 5σ
+    end
 end