void myGraph(void){ 
  gROOT->Reset();//Reset session memory
  gROOT->SetStyle("Plain");//Make style plain
  gStyle->SetOptStat(0);//Remove the statistics box from histogram

  //Make a canvas
  TCanvas * myCanvas = new TCanvas("myCanvas","");

  //For this example, I will plot 5 points
  const Int_t nPoints = 5;
  //Create vectors that have the x and y coordinates 
  Double_t c1 = 1281.9, c1f = 1284.29, c1fErr = 0.213808, s1f = 8.6965, s1fErr = 0.321985; 
  Double_t c2 = 1426.3, c2f = 1428.54, c2fErr = 0.277389, s2f = 9.93108, s2fErr = 0.627892; 
  Double_t c3 = 1518.0, c3f = 1516.98, c3fErr = 0.340637, s3f = 9.49579, s3fErr = 0.107779; 
  Double_t c4 = 1413.9, c4f = 1416.19, c4fErr = 0.267824, s4f = 9.92736, s4fErr = 0.537668; 
  Double_t c5 = 1409.0, c5f = 1413.18, c5fErr = 0.410183, s5f = 10.6733, s5fErr = 1.14787; 
  Double_t xVec[nPoints] = {c1,c2,c3,c4,c5};
  Double_t yVec[nPoints] = {s1f,s2f,s3f,s4f,s5f};
  //Create vectors that have the x and y errors
  Double_t xErrVec[nPoints] = {0,0,0,0,0}; //Common for x-error = 0
  Double_t yErrVec[nPoints] = {s1fErr,s2fErr,s3fErr,s4fErr,s5fErr};
  //Create the TGraphErrors object
  TGraphErrors * myGraph = new TGraphErrors(nPoints,xVec,yVec,xErrVec,yErrVec);
  //Remove the title box
  myGraph->SetTitle("");
  //Set the marker style and color. 
  //See: https://root.cern.ch/doc/master/classTAttMarker.html
  myGraph->SetMarkerStyle(21);
  myGraph->SetMarkerColor(2);
  //Set the marker size
  myGraph->SetMarkerSize(1.0);
  //Draw the graph
  myGraph->Draw("AP");//A->Draw axis, L->Polyline, P->Draw marker
  //Set log-scale on x-axis
  //myCanvas->SetLogx();
  //Set axis titles (Must define AFTER SetLogx)
  myGraph->GetXaxis()->SetTitle("Mass generated (MeV)");
  myGraph->GetYaxis()->SetTitle("Gaussian #sigma (MeV)");


  return;
}