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

Timestamp:

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

Author:

delerue

Message:

Added double gaussian to GFWnew

File:

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

CLIO/CTR+gfw/CTR.m

@@ -1 +1 @@
   function [nu,CTR_spec_1,CTR_spec_N,E_total_1]=CTR(var)
 % clear all; close all; clc
-%    var='PrSelect=1; P_file=''../p20.mat''; PlotFunc=[1 1 1 1 1]; ';
+%    
 %PlotFunc:
 %   1 -- Form factor
@@ -26 +26 @@
 basefig=0;%base number of figures
 icolor='r';%color of lines on plot
+PlotArgTxt='MarkerSize'; %one manual argument for the plots
+PlotArgVal=6; %one manual argument for the plots
 if (~isempty(var))
 disp(var);
@@ -89 +91 @@
 % imin =find(nu>=3e+11,1);%set max freq
 imin=1;
-disp(['Emited energy (from 1um to 1mm) equal=' num2str(sum(CTR_spec_N(imin:end))*(omega(2)-omega(1))) ' J'])
+disp(['Emitted energy (from 1um to 1mm) equal=' num2str(sum(CTR_spec_N(imin:end))*(omega(2)-omega(1))) ' J'])
 E_total_1=sum(CTR_spec_N(imin:end))*(omega(2)-omega(1));
 if PlotFunc(1)==1
 figure(basefig+1)
-loglog(3e+8./nu*1e+3,FFT_spec_om,icolor,'linew',2)
+loglog(3e+8./nu*1e+3,FFT_spec_om,icolor,'linew',2,PlotArgTxt,PlotArgVal)
 hold on
 % loglog(3e+8./(nu2(1:round(length(nu))))*1e+3,FFT_spec_om,'r','linew',2)
@@ -108 +110 @@
 if PlotFunc(2)==1
 figure(basefig+2)
-loglog(c./nu*1e+3,CTR_spec_1*1e+18,icolor,'linew',2)
+loglog(c./nu*1e+3,CTR_spec_1*1e+18,icolor,'linew',2,PlotArgTxt,PlotArgVal)
 hold on
 grid on
@@ -119 +121 @@
 if PlotFunc(3)==1
 figure(basefig+3)
-semilogy(c./nu*1e+3,CTR_spec_N*1e+18,icolor,'linew',2)
+semilogy(c./nu*1e+3,CTR_spec_N*1e+18,icolor,'linew',2,PlotArgTxt,PlotArgVal)
 set(gca,'fontsize',16)
 ylabel('dI/d\omega, [\muJ/THz]')
@@ -137 +139 @@
 if PlotFunc(4)==1
 figure(basefig+4)
-plot(time*1e+3,Profile./max(Profile),icolor,'linew',2)
+plot(time*1e+3,Profile./max(Profile),icolor,'linew',2,PlotArgTxt,PlotArgVal)
 set(gca,'fontsize',16)
 title(['Profile, fwhm=' num2str(pulse_len_ps) ' ps'])