void wwh()
{
   //
   // e+e- -> vvWWh
   //
   TCanvas *c1 = new TCanvas("wwh", "WW fusion", 100,10, 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;

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

   t.DrawLatex(7,6,"e^{-}");
   t.DrawLatex(7,55,"e^{+}");

   // WW fusion propagators
   l = new TLine(45, 15, 50, 30); l->Draw();
   l = new TLine(45, 45, 50, 30); l->Draw();
   t.DrawLatex(42.5, 37.5,"W^{+}");
   t.DrawLatex(42.5, 22.5,"W^{-}");

   // Outgoing vv and h
   l = new TLine(45, 15, 75, 10); l->Draw();
   l = new TLine(45, 45, 75, 50); l->Draw();
   l = new TLine(50, 30, 75, 30); l->Draw();

   t.DrawLatex(80,55,"#bar{#nu_{e}}");
   t.DrawLatex(80,6,"#nu_{e}");
   t.DrawLatex(80,30,"H^{0}");

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