void feynman_D0()
{

  TString outPlotName("FeynDiagrams/D0Decay.eps");

  TCanvas *c1 = new TCanvas("c1", "A canvas", 10,10, 600, 1200);
  c1->Range(0, 0, 100, 195);
  Int_t linsav = gStyle->GetLineWidth();
  gStyle->SetLineWidth(3);
  TLatex t;
  t.SetTextAlign(22);
  t.SetTextSize(0.06);
  TLine * l;
  TCurlyLine *gamma;

  // D0->K-Pi+
  l = new TLine(15, 188, 85, 188); l->Draw();
  l = new TLine(15, 180, 85, 180); l->Draw();
  t.DrawLatex(12,188,"#bar{u}");
  t.DrawLatex(12,180,"c");
  t.DrawLatex(5,184,"D^{0}");
  t.DrawLatex(88,188,"#bar{u}");
  t.DrawLatex(88,180,"s");
  t.DrawLatex(95,184,"K^{-}");
  
  gamma = new TCurlyLine(45, 180, 60, 164);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,161,"W");
  
  l = new TLine(60, 164, 85, 168); l->Draw();
  l = new TLine(60, 164, 85, 160); l->Draw();
  
  t.DrawLatex(88,168,"u");
  t.DrawLatex(88,160,"#bar{d}");
  t.DrawLatex(95,165,"#pi^{+}");
  
  // D0->K-Pi+Pi0
  l = new TLine(15, 143, 50, 143); l->Draw();
  l = new TLine(15, 135, 45, 135); l->Draw();
  t.DrawLatex(12,143,"#bar{u}");
  t.DrawLatex(12,135,"c");
  t.DrawLatex(5,139,"D^{0}");
  l = new TLine(50, 143, 85, 148 ); l->Draw();
  l = new TLine(60, 139, 85, 143); l->Draw();
  l = new TLine(60, 139, 85, 135); l->Draw();
  l = new TLine(45, 135, 85, 130); l->Draw();
  t.DrawLatex(88,148 ,"#bar{u}");
  t.DrawLatex(88,143,"u");
  t.DrawLatex(88,135,"#bar{u}");
  t.DrawLatex(88,130,"s");
  t.DrawLatex(95,146,"#pi^{0}");
  t.DrawLatex(95,132,"K^{-}");

  gamma = new TCurlyLine(45, 135, 60, 114);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,121,"W");
  
  l = new TLine(60, 114, 85, 118); l->Draw();
  l = new TLine(60, 114, 85, 110); l->Draw();
  
  t.DrawLatex(88,118,"u");
  t.DrawLatex(88,110,"#bar{d}");
  t.DrawLatex(95,115,"#pi^{+}");
  
  // D0->K-Pi+Pi+Pi-
  l = new TLine(15, 89, 30, 89); l->Draw();
  l = new TLine(15, 81, 30, 81); l->Draw();
  l = new TLine(30, 89, 50, 95); l->Draw();
  l = new TLine(30, 81, 45, 75); l->Draw();
  t.DrawLatex(12,89,"#bar{u}");
  t.DrawLatex(12,81,"c");
  t.DrawLatex(5,85,"D^{0}");
  l = new TLine(50, 95, 85, 100 ); l->Draw();
  l = new TLine(60, 91, 85, 95); l->Draw();
  l = new TLine(60, 91, 85, 87); l->Draw();
  l = new TLine(60, 79, 85, 83); l->Draw();
  l = new TLine(60, 79, 85, 75); l->Draw();
  l = new TLine(45, 75, 85, 70); l->Draw();
  t.DrawLatex(88,100 ,"#bar{u}");
  t.DrawLatex(88,95,"d");
  t.DrawLatex(95,98,"#pi^{-}");
  t.DrawLatex(88,87 ,"#bar{d}");
  t.DrawLatex(88,83,"u");
  t.DrawLatex(95,85,"#pi^{+}");
  t.DrawLatex(88,75,"#bar{u}");
  t.DrawLatex(88,70,"s");
  t.DrawLatex(95,72,"K^{-}");

  gamma = new TCurlyLine(45, 75, 60, 54);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,61,"W");
  
  l = new TLine(60, 54, 85, 58); l->Draw();
  l = new TLine(60, 54, 85, 50); l->Draw();
  
  t.DrawLatex(88,58,"u");
  t.DrawLatex(88,50,"#bar{d}");
  t.DrawLatex(95,55,"#pi^{+}");
  
  // D0->K0Pi+Pi-
  l = new TLine(15, 38, 50, 38); l->Draw();
  l = new TLine(15, 30, 45, 30); l->Draw();
  t.DrawLatex(12,38,"#bar{u}");
  t.DrawLatex(12,30,"c");
  t.DrawLatex(5,34,"D^{0}");
  l = new TLine(50, 38, 85, 43 ); l->Draw();
  l = new TLine(60, 34, 85, 38); l->Draw();
  l = new TLine(60, 34, 85, 30); l->Draw();
  l = new TLine(45, 30, 85, 25); l->Draw();
  t.DrawLatex(88,43 ,"#bar{u}");
  t.DrawLatex(88,38,"d");
  t.DrawLatex(88,30,"#bar{d}");
  t.DrawLatex(88,25,"s");
  t.DrawLatex(95,41,"#pi^{-}");
  t.DrawLatex(95,27,"#bar{K}^{0}");

  gamma = new TCurlyLine(45, 30, 60, 9);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,16,"W");
  
  l = new TLine(60, 9, 85, 13); l->Draw();
  l = new TLine(60, 9, 85, 5); l->Draw();
  
  t.DrawLatex(88,13,"u");
  t.DrawLatex(88,5,"#bar{d}");
  t.DrawLatex(95,10,"#pi^{+}");
  

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