void feynman_BToDPilnu()
{

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

  TCanvas *c1 = new TCanvas("c1", "A canvas", 10,10, 600, 300);
  c1->Range(0, 0, 100, 50);
  Int_t linsav = gStyle->GetLineWidth();
  gStyle->SetLineWidth(3);
  TLatex t;
  t.SetTextAlign(22);
  t.SetTextSize(0.1);
  TLine * l;
  l = new TLine(15, 40, 50, 40); l->Draw();
  l = new TLine(15, 32, 45, 32); l->Draw();
  t.DrawLatex(12,40,"d");
  t.DrawLatex(12,32,"#bar{b}");
  t.DrawLatex(5,36,"B^{0}");
  l = new TLine(50, 40, 85, 45); l->Draw();
  l = new TLine(60, 36, 85, 40); l->Draw();
  l = new TLine(60, 36, 85, 32); l->Draw();
  l = new TLine(45, 32, 85, 27); l->Draw();
  t.DrawLatex(88,45,"d");
  t.DrawLatex(88,40,"#bar{d}");
  t.DrawLatex(88,32,"d");
  t.DrawLatex(88,27,"#bar{c}");
  t.DrawLatex(95,44,"#pi^{0}");
  t.DrawLatex(95,31,"D^{-}");
  
  TCurlyLine *gamma = new TCurlyLine(45, 32, 60, 10);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,18,"W");
  
  l = new TLine(60, 10, 85, 15); l->Draw();
  l = new TLine(60, 10, 85, 5); l->Draw();
  
  t.DrawLatex(88,15,"l^{+}");
  t.DrawLatex(88,5,"#nu");
  
  c1->Update();
  gStyle->SetLineWidth(linsav);
  
  c1->Print(outPlotName);
}