#set term postscript color enhanced
#xsset output 'kicker_pulse.ps'
cc              =0
aa              = -2.84411      
bb              = 3.3351        
dd              = 1.27428       
ee              = 0.0970149     
ff              = -3.92405        
set xlabel 'Time [ns]
set ylabel 'Amplitude [normalized]'
gg(x) = bb*tanh(aa*(x-dd))+cc +bb*tanh(-ff*(x+ee))
hh(x) = bb*tanh(aa*(x-dd+0.5))+cc +bb*tanh(-ff*(x+ee+0.5))
normalization = hh(0)
#set yrange [-0.5:]
#set xrange [-100:]
#fit gg(x) 'pulse_127cm.dat' via aa,bb,dd,ee,ff
#set xrange [180:800]
a=0.02
k1=0.12
l1=4.2
p1=1.
k2=-0.025
l2=4.2
p2=1.
q1=-0.07
q2=-0.01
#set xrange [1.8:8.]
t(x) = a+k1*sin(l1*x+p1)*exp(q1*x) + k2*sin(l2*x+p2)*exp(q2*x)
#fit t(x) 'pulse_tail_127cm.dat' via k1,l1,p1,q1 #,k2,l2,p2,q2
#set xrange [180:800]
#plot gg((x/100))/hh(0) lw 2 t 'functional form', 'pulse_127cm.dat'  u ($1*100):($2/normalization) ps 1.5 pt 5 t 'data'
plot t((x/100)) lw 2 t 'functional form', 'pulse_tail_127cm.dat'  u ($1*100):($2/normalization) ps 1.5 pt 5 t 'data'


#, hh(x)/hh(0)
#fit g(x) 'pulse.dat' via aa,bb,cc,dd,ee,a,b,c,d,e,f,g,h,j
#fit g(x) 'pulse.dat' via a,b,c,d,e,f,g,h,j