#include <TAttFill.h>
#include <TEllipse.h>
#include <TGaxis.h>
#include <TStyle.h>
#include <TLatex.h>
{
//  gROOT->Reset();
  gROOT->SetStyle("Plain");

  gStyle->SetOptStat(1111);
  gStyle->SetPalette(1);

  c1 = new TCanvas ("c1","",200,10,800,400);
  c1->Divide(2,1);

  TFile fb1a("purgmag_b1_400keV.root");
  TFile fb1b("purgmag_b1_600keV.root");
  TFile fb1c("purgmag_b1_800keV.root");
  TFile fa1a("purgmag_a1_400keV.root");
  TFile fa1b("purgmag_a1_600keV.root");
  TFile fa1c("purgmag_a1_800keV.root");

  TNtuple* ntupleB1a;
  TNtuple* ntupleB1b;
  TNtuple* ntupleB1c;
  TNtuple* ntupleA1a;
  TNtuple* ntupleA1b;
  TNtuple* ntupleA1c;

  ntupleB1a = (TNtuple*)fb1a.Get("n101");
  ntupleB1b = (TNtuple*)fb1b.Get("n101");
  ntupleB1c = (TNtuple*)fb1c.Get("n101");
  ntupleA1a = (TNtuple*)fa1a.Get("n101");
  ntupleA1b = (TNtuple*)fa1b.Get("n101");
  ntupleA1c = (TNtuple*)fa1c.Get("n101");

  // Plot #1 
  c1->cd(1);
  ntupleB1c->SetMarkerSize(.2);
  ntupleB1c->SetMarkerColor(4);
  ntupleB1c->SetMarkerStyle(20);
  ntupleB1c->Draw("ey:ex","");
  ntupleB1b->SetMarkerSize(.2);
  ntupleB1b->SetMarkerColor(3);
  ntupleB1b->SetMarkerStyle(20);
  ntupleB1b->Draw("ey:ex","","same");
  ntupleB1a->SetMarkerSize(.2);
  ntupleB1a->SetMarkerColor(2);
  ntupleB1a->SetMarkerStyle(20);
  ntupleB1a->Draw("ey:ex","","same");
  TLatex *title1 = new TLatex(-10.0,gPad->GetUymax()+0.5,"Hit Pos [cm] Through B1");
  title1->SetTextSize(0.05);
  title1->Draw();

  // Draw outline of incident beam... 
  TEllipse *circle2 = new TEllipse(-1.4,-1.4,1.27,1.27);
  circle2->SetFillColorAlpha(kWhite, 0.05);
//  circle2->Draw("same");
  circle2->Draw();

  // Plot #2 
  c1->cd(2);
  ntupleA1c->SetMarkerSize(.2);
  ntupleA1c->SetMarkerColor(4);
  ntupleA1c->SetMarkerStyle(20);
  ntupleA1c->Draw("ey:ex","");
  ntupleA1b->SetMarkerSize(.2);
  ntupleA1b->SetMarkerColor(3);
  ntupleA1b->SetMarkerStyle(20);
  ntupleA1b->Draw("ey:ex","","same");
  ntupleA1a->SetMarkerSize(.2);
  ntupleA1a->SetMarkerColor(2);
  ntupleA1a->SetMarkerStyle(20);
  ntupleA1a->Draw("ey:ex","","same");
  TLatex *title2 = new TLatex(-9.0,gPad->GetUymax()+0.5,"Hit Pos [cm] Through A1");
  title2->SetTextSize(0.05);
  title2->Draw();

  // Draw outline of incident beam... 
  TEllipse *circle3 = new TEllipse(1.4,1.4,1.27,1.27);
  circle3->SetFillColor(0);
  circle3->Draw();

}