void feynman_DstarCh()
{

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

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

  // D*->DPi+
  l = new TLine(15, 72, 50, 72); l->Draw();
  l = new TLine(15, 64, 50, 64); l->Draw();
  t.DrawLatex(12,72,"#bar{d}");
  t.DrawLatex(12,64,"c");
  t.DrawLatex(5,68,"D^{*+}");
  l = new TLine(50, 72, 85, 77); l->Draw();
  l = new TLine(60, 68, 85, 72); l->Draw();
  l = new TLine(60, 68, 85, 64); l->Draw();
  l = new TLine(50, 64, 85, 59); l->Draw();
  t.DrawLatex(88,77,"#bar{d}");
  t.DrawLatex(88,72,"u");
  t.DrawLatex(88,64,"#bar{u}");
  t.DrawLatex(88,59,"c");
  t.DrawLatex(95,76,"#pi^{+}");
  t.DrawLatex(95,63,"D^{0}");

  
  // D*->DPi0
  l = new TLine(15, 44, 50, 44); l->Draw();
  l = new TLine(15, 36, 50, 36); l->Draw();
  t.DrawLatex(12,44,"#bar{d}");
  t.DrawLatex(12,36,"c");
  t.DrawLatex(5,40,"D^{*+}");
  l = new TLine(50, 44, 85, 49); l->Draw();
  l = new TLine(60, 40, 85, 44); l->Draw();
  l = new TLine(60, 40, 85, 36); l->Draw();
  l = new TLine(50, 36, 85, 31); l->Draw();
  t.DrawLatex(88,49,"#bar{d}");
  t.DrawLatex(88,44,"d");
  t.DrawLatex(88,36,"#bar{d}");
  t.DrawLatex(88,31,"c");
  t.DrawLatex(95,48,"#pi^{0}");
  t.DrawLatex(95,35,"D^{+}");

  // D*->Dgamma
  l = new TLine(15, 21, 85, 21); l->Draw();
  l = new TLine(15, 13, 85, 13); l->Draw();
  t.DrawLatex(12,21,"#bar{d}");
  t.DrawLatex(12,13,"c");
  t.DrawLatex(5,17,"D^{*+}");
  t.DrawLatex(88,21,"#bar{d}");
  t.DrawLatex(88,13,"c");
  t.DrawLatex(95,17,"D^{+}");
  
  gamma = new TCurlyLine(55, 13, 85, 5);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(88,5,"#gamma");
  
  

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