;$Id: pST_kpq_transfer.pro,v 1.5 2014/07/20 08:23:08 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/pST_kpq_transfer.ps
;
!p.charsize=1.7
!x.margin=[6.0,0.5]
!y.margin=[3.2,0.5]
xx=2.3 & yy=19.
col=255
siz=1.8
;
;  read
;
restore,'SM_kpq.sav'
restore,'TM_kpq.sav'
;
m=24
ki=1+findgen(m)
print,ki
;
sk=1./60.
!p.multi=[0,2,2]
!x.range=sk*[1,24]
;
Pi_SS=fltarr(m)
Pi_LS=fltarr(m)
Pi_SL=fltarr(m)
Pi_LL=fltarr(m)
for k=0,m-1 do begin
for p=0,m-1 do begin
for q=0,m-1 do begin
  Transf=TM_kpq(k,p,q)-TM_kpq(p,k,q)-SM_kpq(p,k,q)
  if p lt k and q lt k then Pi_LL(k)=Pi_LL(k)+Transf
  if p lt k and q ge k then Pi_LS(k)=Pi_LS(k)+Transf
  if p ge k and q le k then Pi_SL(k)=Pi_SL(k)+Transf
  if p ge k and q ge k then Pi_SS(k)=Pi_SS(k)+Transf
endfor
endfor
endfor
;
;plot,ki,Pi_LL,ytit='LL'
;plot,ki,Pi_LS,ytit='LS'
;plot,ki,Pi_SL,ytit='SL'
;plot,ki,Pi_SS,ytit='SS'
;
!p.multi=0
!x.title='!8k!6/!8k!6!d0!n'
Pi_tot=Pi_LL+Pi_LS+Pi_SL+Pi_SS
s=1./max(abs(Pi_tot))
plot,sk*ki,s*Pi_LL,yr=[-.4,.5]
oplot,sk*ki,s*Pi_LS,li=1
oplot,sk*ki,s*Pi_SL,li=3
oplot,sk*ki,s*Pi_SS,li=2
;oplot,sk*ki,s*(Pi_LS+Pi_SL),li=2,col=122
;
xyouts,siz=siz,sk*17,.37,'!7P!6!s!d!8p!6<!8q!6!n!r!u!8q!6<!8k!6!n'
xyouts,siz=siz,sk*17,.04,'!7P!6!s!d!8p!6>!8q!6!n!r!u!8q!6>!8k!6!n'
xyouts,siz=siz,sk*15,-.14,'!7P!6!s!d!8p!6<!8q!6!n!r!u!8q!6>!8k!6!n'
xyouts,siz=siz,sk*17.2,-.3,'!7P!6!s!d!8p!6>!8q!6!n!r!u!8q!6<!8k!6!n'
END