void eeqqpp()
{
   //
   // e+e- -> qqbar with double ISR
   //
   TCanvas *c1 = new TCanvas("eeqqpp", "e+e- -> qqpp", 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");

   // photon radiations
   lc = new TCurlyLine(16.5,16, 30.8, 16); l->Draw();
   lc->SetWavy();
   lc->Draw();
   lc = new TCurlyLine(16.5,44, 30.8, 44); l->Draw();
   lc->SetWavy();
   lc->Draw();

   t.DrawLatex(33.5,16,"#gamma");
   t.DrawLatex(33.5,44,"#gamma");

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