void DrawMassError(TString nSt="") {

  TStopwatch timer;
  timer.Start();

//  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
//  basiclibs();

TStyle *plain = new TStyle("Plain","Plain Style(no colors/fill areas)");
//plain->SetCanvasBorderMode(0);
//plain->SetPadBorderMode(0);
plain->SetPadColor(0);
plain->SetCanvasColor(0);
plain->SetTitleColor(0);
plain->SetStatColor(0);
plain->SetOptStat(110010);
plain->SetOptFit(0011);
plain->SetStatW(0.225);
plain->SetStatH(0.06);
//plain->SetStatX(0.7);
plain->SetOptTitle(0);
plain->cd();

  TString name = "out.root";
  TFile *dirOld = new TFile("out_old.root","read");
  TFile *dirNew = new TFile("out_new.root","read");

  TH1D *hMassErrOld, *hMassErrNew;
  TH1D *hEErrOld, *hEErrNew;
    
  hMassErrOld = (TH1D*) dirOld->Get("SigmaMass");
  hMassErrOld->SetLineColor(1);
  hEErrOld = (TH1D*) dirOld->Get("SigmaE");
  hEErrOld->SetLineColor(1);

  hMassErrNew = (TH1D*) dirNew->Get("SigmaMass");
  hMassErrNew->SetLineColor(1);
  hEErrNew = (TH1D*) dirNew->Get("SigmaE");
  hEErrNew->SetLineColor(1);
  
  TCanvas *c1;
  c1 = new TCanvas("c1","position",0,0,1200,800);
//  c1 -> Divide(NParameters-2,2);
  c1 -> Divide(2,2);

  c1->UseCurrentStyle();

  c1->cd(1);
  hMassErrOld->GetXaxis()->SetTitle("Old Mass Error D^{0} [GeV/c^{2}]");
  hMassErrOld->GetXaxis()->SetTitleOffset(1);
  hMassErrOld->Draw();
  
  c1->cd(2);
  hEErrOld->GetXaxis()->SetTitle("Old E Error D^{0} [GeV/c^{2}]");
  hEErrOld->GetXaxis()->SetTitleOffset(1);
  hEErrOld->Draw();

  c1->cd(3);
  hMassErrNew->GetXaxis()->SetTitle("New Mass Error D^{0} [GeV/c^{2}]");
  hMassErrNew->GetXaxis()->SetTitleOffset(1);
  hMassErrNew->Draw();
  
  c1->cd(4);
  hEErrNew->GetXaxis()->SetTitle("New E Error D^{0} [GeV/c^{2}]");
  hEErrNew->GetXaxis()->SetTitleOffset(1);
  hEErrNew->Draw();

  c1->SaveAs("MassErrorComparison.pdf");
 
  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);

}