void eeqqgg()
{
   //
   // e+e- -> qqbar with double gluon radiation
   //
   TCanvas *c1 = new TCanvas("eeqqgg", "e+e- -> qqgg", 740,40, 400, 300);
   c1->Range(0, 0, 90, 60);
   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(7,6,"e^{-}");
   t.DrawLatex(7,55,"e^{+}");

   // Z propagator
   lc = new TCurlyLine(30, 30, 55, 30);
   lc->SetWavy();
   lc->Draw();
   t.DrawLatex(42.5, 25.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(80,6,"#bar{q}");
   t.DrawLatex(80,55,"q");

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

   t.DrawLatex(80,24,"g");
   t.DrawLatex(80,36,"g");

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