void bhabha()
{
   //
   // Draws the feynman diagrams for D0->mu+mu- via Rp violating SUSY.
   //
   TCanvas *c1 = new TCanvas("c1", "A canvas", 10, 10, 900, 300);
   c1->Range(0, 0, 150, 50);
   Int_t linsav = gStyle->GetLineWidth();
   gStyle->SetLineWidth(2);
   TLatex t;
   t.SetTextAlign(22);
   t.SetTextSize(0.1);
   TLine * l;
   //
   // The ee --> ee diagram
   //
   l = new TLine(10, 10, 50, 10); l->Draw();
   l = new TLine(10, 40, 50, 40); l->Draw();
   g = new TCurlyLine(30, 10, 30, 40, 0.08, 0.03); g->SetWavy(); g->Draw();
   //   g = new TCurlyArc(30, 10, 30, 40, 0.05, 0.02); g->Draw();
   t.DrawLatex(7,10,"e^{-}");
   t.DrawLatex(7,40,"e^{+}");
   t.DrawLatex(55,10,"e^{-}");
   t.DrawLatex(55,40,"e^{+}");
   t.DrawLatex(35,25,"#gamma");
   t.DrawLatex(30,7,"#sqrt{#alpha}");
   t.DrawLatex(30,43,"#sqrt{#alpha}");
   //
   // put a plus sign between the two
   //
   t.DrawLatex(65,25,"+");
   //
   // The crossed diagram
   //
   l = new TLine(80, 10, 95, 25); l->Draw();
   //   l->SetLineStyle(1); // dashed line
   l = new TLine(80, 40, 95, 25); l->Draw();
   g = new TCurlyLine(95, 25, 120, 25, 0.08, 0.03); g->SetWavy(); g->Draw();
   l = new TLine(120, 25, 135, 10); l->Draw();
   l = new TLine(120, 25, 135, 40); l->Draw();
   t.DrawLatex(77,10,"e^{-}");
   t.DrawLatex(77,40,"e^{+}");
   t.DrawLatex(140,10,"e^{-}");
   t.DrawLatex(140,40,"e^{+}");
   t.DrawLatex(107,20,"#gamma");
   t.DrawLatex(95,20,"#alpha");
   t.DrawLatex(120,20,"#alpha");
   
   c1->Update();
   gStyle->SetLineWidth(linsav);
}