{
  TF1 *erf = new TF1("erf", "(TMath::Erf(0.4*(x-6))+1)/2", 0, 20);
  TF1 *fgaus = new TF1("fgaus","TMath::Gaus(x, 8, 1.5)",0,20);
  fgaus->SetNpx(10000);
  fgaus->SetFillColor(19);
  fgaus->SetFillStyle(0);
  fgaus->SetLineColor(4);
  fgaus->SetLineWidth(3);
  fgaus->SetTitle("I am not working");
  fgaus->GetXaxis()->SetTitle("Jet p_{T} (GeV)");
  fgaus->Draw();

  TF1 * gaus_erf = new TF1("gaus_erf", "fgaus*erf", 0, 20);
  gaus_erf->Draw();
}