void feynman_Dch()
{

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

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

  // D+->K-Pi+Pi+
  l = new TLine(15, 78, 50, 78); l->Draw();
  l = new TLine(15, 70, 45, 70); l->Draw();
  t.DrawLatex(12,78,"#bar{d}");
  t.DrawLatex(12,70,"c");
  t.DrawLatex(5,74,"D^{+}");
  l = new TLine(50, 78, 85, 83 ); l->Draw();
  l = new TLine(60, 74, 85, 78); l->Draw();
  l = new TLine(60, 74, 85, 70); l->Draw();
  l = new TLine(45, 70, 85, 65); l->Draw();
  t.DrawLatex(88,83 ,"#bar{d}");
  t.DrawLatex(88,78,"u");
  t.DrawLatex(88,70,"#bar{u}");
  t.DrawLatex(88,65,"c");
  t.DrawLatex(95,81,"#pi^{+}");
  t.DrawLatex(95,67,"K^{-}");

  gamma = new TCurlyLine(45, 70, 60, 49);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,56,"W");
  
  l = new TLine(60, 49, 85, 53); l->Draw();
  l = new TLine(60, 49, 85, 45); l->Draw();
  
  t.DrawLatex(88,53,"u");
  t.DrawLatex(88,45,"#bar{d}");
  t.DrawLatex(95,50,"#pi^{+}");
  
  // D+->K0Pi+
  l = new TLine(15, 35, 85, 35); l->Draw();
  l = new TLine(15, 27, 85, 27); l->Draw();
  t.DrawLatex(12,35,"#bar{d}");
  t.DrawLatex(12,27,"c");
  t.DrawLatex(5,31,"D^{+}");
  t.DrawLatex(88,35,"#bar{d}");
  t.DrawLatex(88,27,"s");
  t.DrawLatex(95,31,"#bar{K}^{0}");
  
  gamma = new TCurlyLine(45, 27, 60, 11);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,18,"W");
  
  l = new TLine(60, 11, 85, 15); l->Draw();
  l = new TLine(60, 11, 85, 7); l->Draw();
  
  t.DrawLatex(88,15,"u");
  t.DrawLatex(88,7,"#bar{d}");
  t.DrawLatex(95,12,"#pi^{+}");
  

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