void feynman_Bch()
{

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

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

  // B->Dlnu
  l = new TLine(15, 145, 85, 145); l->Draw();
  l = new TLine(15, 137, 85, 137); l->Draw();
  t.DrawLatex(12,145,"u");
  t.DrawLatex(12,137,"#bar{b}");
  t.DrawLatex(5,141,"B^{+}");
  t.DrawLatex(88,145,"u");
  t.DrawLatex(88,137,"#bar{c}");
  t.DrawLatex(95,141,"#bar{D}^{0}");
  
  gamma = new TCurlyLine(45, 137, 60, 120);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,128,"W");
  
  l = new TLine(60, 120, 85, 125); l->Draw();
  l = new TLine(60, 120, 85, 115); l->Draw();
  
  t.DrawLatex(88,125,"l^{+}");
  t.DrawLatex(88,115,"#nu");
  
  // B->DPi0
  l = new TLine(15, 95, 50, 95); l->Draw();
  l = new TLine(15, 87, 45, 87); l->Draw();
  t.DrawLatex(12,95,"u");
  t.DrawLatex(12,87,"#bar{b}");
  t.DrawLatex(5,91,"B^{+}");
  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(45, 87, 85, 82); l->Draw();
  t.DrawLatex(88,100,"u");
  t.DrawLatex(88,95,"#bar{u}");
  t.DrawLatex(88,87,"u");
  t.DrawLatex(88,82,"#bar{c}");
  t.DrawLatex(95,99,"#pi^{0}");
  t.DrawLatex(95,86,"#bar{D}^{0}");

  gamma = new TCurlyLine(45, 87, 60, 65);
  gamma->SetWavy();
  gamma->Draw();
  t.DrawLatex(45,73,"W");
  
  l = new TLine(60, 65, 85, 70); l->Draw();
  l = new TLine(60, 65, 85, 60); l->Draw();
  
  t.DrawLatex(88,70,"l^{+}");
  t.DrawLatex(88,60,"#nu");
  
  // B->DPi-
  l = new TLine(15, 40, 50, 40); l->Draw();
  l = new TLine(15, 32, 45, 32); l->Draw();
  t.DrawLatex(12,40,"u");
  t.DrawLatex(12,32,"#bar{b}");
  t.DrawLatex(5,36,"B^{+}");
  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,"u");
  t.DrawLatex(88,40,"#bar{d}");
  t.DrawLatex(88,32,"d");
  t.DrawLatex(88,27,"#bar{c}");
  t.DrawLatex(95,44,"#pi^{+}");
  t.DrawLatex(95,31,"D^{-}");

  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);
}