!
! Gluon exchange diagrams for parton-parton scattering
!
exec lecture_fig_setup fragfunc 10 10 2 1 1 1.5
set font -1001
set font title -1004
wind 2 1 2 1
set y marg =+0.5
set tsize 0.5
!set title off
!
!exec lecture_fig_copy 'Fragmentation Functions' 2
!
!set tick all off
!set x grid on all 3
!set y grid on all 3
!
set x label 'z' =
!
fun del 0
fun add comis frag_lq.f 0 n
2
Light Quark Fragmentation Function
a
b
3.5
2.5
fun add frag 1
1
0.016
!
set function bin_width off
set fit area on
set x lim 0 1
set y mode log 0.00001 7
fun plot 1 &1 -1/red
0.02 1
title 98765&1 'Light Quarks'
comment 98765&1 
  new 'a/z (1-z)^b!' 6.4 11   = = l
  new 'a=3.5, b=2.5' 6.4 10.4 = = l
end
!
set y mode int
fun plot 2 &2 -1/red
0 1
title 98765&2 'Heavy Quarks'
!
dep p2(2) = 0.15
fun over 2 &3 -1/blue
!
comment 98765&2 
  new '[e]?b! = 0.016' 17.5 11   = = l = = red
  new '[e]?c! = 0.15'  17.5 10.4 = = l = = blue
end
!
redraw
!
@9 hard epost
@9 close
!
! Heavy Quarks
!
set hard peterson.eps
set y mode int
!
dep p2(2) = 0.016
fun plot 2 &3 -1/red
0 1
title 98765&3 'Heavy Quarks'
!
dep p2(2) = 0.15
fun over 2 &4 -1/blue
!
comment 98765&3 
  new '[e]?b! = 0.016' 6.5 11 = = l = = red
  new '[e]?c! = 0.15'  6.5 10.4 = = l = = blue
end
!
dep p2(2) = 0.0023
fun plot 2 &5 -1/red
0 1
title 98765&5 'Heavy Quarks'
!
dep p2(2) = 0.022
fun over 2 &6 -1/blue
!
comment 98765&5 
  new '[e]?b! = 0.0023' 17.5 11 = = l = = red
  new '[e]?c! = 0.022'  17.5 10.4 = = l = = blue
end
!
redraw
!
@9 hard epost
@9 close