Changeset 737 in ETALON for CLIO/CTR+gfw/GFWnew/GFWnew.m

Timestamp:

Sep 1, 2017, 6:03:25 PM (7 years ago)

Author:

delerue

Message:

Added double gaussian to GFWnew

File:

• CLIO/CTR+gfw/GFWnew/GFWnew.m (modified) (3 diffs)

CLIO/CTR+gfw/GFWnew/GFWnew.m

@@ -1 +1 @@
   function  [nu_Hz,SP_spec_1,SP_spec_N,lam,E_total,sum_tr,thet]=GFWnew(var)
+global Profile
+global time
 
 %  single_elec
@@ -46 +48 @@
 ThetaDiv=10;%steps in theta
 phiDiv=3;%steps in phi
+double_gaussian_offset=0; % Offset of the second gaussian with respect to the first one
+double_gaussian_amplitude=1; % Amplitude of the second gaussian with respect to the first one
+double_gaussian_sigma_factor=1; % Sigma ration of the second gaussian with respect to the first one
 
 if (~isempty(var))
@@ -74 +79 @@
 sigma_L= sigma_t*epsl;
 sigma_H= sigma_t;
+double_gaussian_sigma=sigma_t*double_gaussian_sigma_factor;
 time=(1:(length(nu)*2))*(1/(length(nu)*2*(nu(2)-nu(1))));
 xc=mean(time);
 G1=exp(-(time-xc).^2./(2*sigma_L^2));
 G2=exp(-(time-xc).^2./(2*sigma_H^2));
+G3=double_gaussian_amplitude*exp(-(time-xc+double_gaussian_offset).^2./(2*double_gaussian_sigma^2));
 G1(time>xc)=0;
 G2(time<xc)=0;
-Profile=(G1+G2)./sum(G1+G2);
+Profile=(G1+G2+G3)./sum(G1+G2+G3);
 end