{
  gSystem->Load("kkp/libkkp.so");
  gROOT->LoadMacro("kkp_func.C");
 
  const Float_t Eparton = 20;

  TCanvas *c1 = new TCanvas("c1","c1: Frag f",500,800);
  gPad->SetLogy();
  h1=gPad->DrawFrame(0,1e-5,1,1e3);
  h1->SetXTitle("z");
  h1->SetYTitle("dN/dz");

  TF1 *frag_f = new TF1("frag_f",kkp_func,0,1,2);
  frag_f->SetParameter(0,Eparton); // Q for fragmentation
  frag_f->SetParameter(1,1); // parton flavour 1=u

  frag_f->SetLineWidth(1);  
  frag_f->SetLineStyle(2);
  frag_f->Draw("same");

  TF1 *frag_f_glue = new TF1("frag_f_glue",kkp_func,0,1,2);
  frag_f_glue->SetParameter(0,Eparton); // Q for fragmentation
  frag_f_glue->SetParameter(1,0); // parton flavour 0=g

  frag_f_glue->SetLineStyle(2);
  frag_f_glue->SetLineWidth(1);
  frag_f_glue->Draw("same");
}