void eeqqggg()
{
   //
   // e+e- -> qqbar with single gluon radiation plus splitting
   //
   TCanvas *c1 = new TCanvas("eeqqggg", "e+e- -> qqggg", 740,40, 400, 300);
   c1->SetFillColor(0);
   c1->Range(0, 0, 100, 60);
   c1->Range(5, 5,  95, 55);
   c1->SetFrameLineColor(0);
   c1->SetFillStyle(0);
   Int_t linsav = gStyle->GetLineWidth();
   gStyle->SetLineWidth(3);
   TLatex t;
   t.SetTextAlign(22);
   t.SetTextSize(0.1);
   TLine * l;
   TCurlyLine * lc; 

   // Incoming e+e-
   l = new TLine(10, 10, 30, 30); l->Draw();
   l = new TLine(10, 50, 30, 30); l->Draw();

   t.DrawLatex(23.0,44,"e^{-}");
   t.DrawLatex(23.0,16,"e^{+}");

   // Z propagator
   lc = new TCurlyLine(30, 30, 55, 30);
   lc->SetWavy();
   lc->Draw();
   t.DrawLatex(42.5, 24.0,"Z^{0}/#gamma");

   // Outgoing qqbar
   l = new TLine(55, 30, 75, 10); l->Draw();
   l = new TLine(55, 30, 75, 50); l->Draw();

   t.DrawLatex(63.,44.,"#bar{q}");
   t.DrawLatex(63.,16.,"q");

   // gluon radiations
   lc = new TCurlyLine(60.8,36, 75, 36); l->Draw();
   lc->SetCurly();
   lc->Draw();
   lc = new TCurlyLine(75, 36, 86, 47); l->Draw();
   lc->SetCurly();
   lc->Draw();
   lc = new TCurlyLine(75, 36, 86, 25); l->Draw();
   lc->SetCurly();
   lc->Draw();

   t.DrawLatex(69,32,"g");
   t.DrawLatex(86,41,"g");
   t.DrawLatex(86,32,"g");

   c1->Update();
   gStyle->SetLineWidth(linsav);
}