280225 commit

Author: stefanocovino

Project.toml

@@ -1,7 +1,7 @@
 name = "FittingFunction"
 uuid = "7c8a6357-aefd-4570-aae0-a8cd991d38f9"
 authors = ["Stefano Covino <stefano.covino@inaf.it> and contributors"]
-version = "0.14.0"
+version = "0.15.0"
 
 [deps]
 AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"

src/FittingFunction.jl

@@ -24,6 +24,7 @@ export Gaussian
 export GetAtomicData
 export Jy2Ecm2sA
 export Mag2Counts
+export NorrisPulse
 export PL
 export Pol2Stokes
 export SBPL
@@ -221,7 +222,7 @@ Compute the discrete Fourier periodogram for the inpout signal.
 
 - `signal` array of input data.
 - `fs` sampling in frequency of the input data (1/dt).
-- 'zerofreq' is true (false) to (not) include the zero frequency in the output. 
+- 'zerofreq' is true (false) to (not) include the zero frequency in the output.
 
 Outputs are two arrays: the frequencies and the powers.
 
@@ -420,6 +421,64 @@ end
 
 
 
+
+"""
+    NorrisPulse(x;pulsNorm,tmax,σ_rise,σ_decay,pulsSharpness,base)
+
+Compute an asymmetric pulse shape according to the recipe reported in [Norris et al. (1996).](https://ui.adsabs.harvard.edu/abs/1996ApJ...459..393N/abstract).
+
+# Arguments
+
+- `x` input vector.
+- `pulsNorm` the nornalization of the pulse.
+- `tmax` the maximum of the pulse.
+- `σ_rise` the time-scale of the rising phase.
+- `σ_decay` the time-scale of the decaying phase.
+- `pulsSharpness` pulse sharpness parameter.
+- `base` signal level without the pulse.
+
+
+# Examples
+```jldoctest
+
+NorrisPulse([1.,2.,3.,4.,5.],pulsNorm=1.5,tmax=3.,σ_rise=1.,σ_decay=2.,pulsSharpness=1.1,base=0.3)
+
+# output
+
+5-element Vector{Float64}:
+ 0.47585740398559895
+ 0.8518191617571635
+ 1.8
+ 1.2407748943132098
+ 0.8518191617571635
+```
+"""
+function NorrisPulse(x;pulsNorm,tmax,σ_rise,σ_decay,pulsSharpness,base)
+    #
+    f1(x) = base .+ pulsNorm .* exp.(-((abs.(x .- tmax)) ./ σ_rise).^pulsSharpness)
+    f2(x) = base .+ pulsNorm .* exp.(-((abs.(x .- tmax)) ./ σ_decay).^pulsSharpness)
+    #
+    fv = [f1,f2]
+    #
+    return piecewise(x,[tmax,],fv)
+end
+
+
+
+function piecewise(x::Float64, breakpts::Vector{Float64}, f::Vector{Function})
+    @assert(issorted(breakpts))
+    @assert(length(f) == length(breakpts)+1)
+    b = searchsortedfirst(breakpts, x)
+    return f[b](x)
+end
+
+piecewise(X::Vector{Float64}, bpts, f) = [piecewise(x,bpts,f) for x in X ]
+
+
+
+
+
+
 """
     PL(E,N,α;E0=1.)
 

test/runtests.jl

@@ -50,4 +50,6 @@ using Test
     #
     @test Z2N([1.,0.5,0.25], [1.,2.,2.5,3.5,5.]) == [0.4, 0.4000000000000002, 0.537258300203048]
     #
+    @test NorrisPulse([1.,2.,3.,4.,5.],pulsNorm=1.5,tmax=3.,σ_rise=1.,σ_decay=2.,pulsSharpness=1.1,base=0.3) == [0.47585740398559895, 0.8518191617571635, 1.8, 1.2407748943132098, 0.8518191617571635]
+    #
 end