// Macro to make plot showing where low and hi mA come from 
// in Eglu1 vs Eglu2 space.
// Requires that you have already run A_to_tautau.C analysis and made
// a pointer t=a->GetHist()->HiggsTree
{
  TH2F *Elt100=new TH2F("Elt100","E_{g2} vs E_{g1} for (tan#beta=40, m_{A}=110 GeV) in Pythia v6.155",50,0.0,300.0,50,0.0,300.0);
  TH2F *Egt100=new TH2F("Egt100","E_{g2} vs E_{g1} for (tan#beta=40, m_{A}=110 GeV) in Pythia v6.155",50,0.0,300.0,50,0.0,300.0);
  Elt100->SetStats(0);
  Egt100->SetStats(0);

  Elt100->SetMarkerStyle(20);
  Elt100->SetMarkerSize(0.4);
  Elt100->SetMarkerColor(2);

  Egt100->SetMarkerStyle(20);
  Egt100->SetMarkerSize(0.4);
  Egt100->SetMarkerColor(1);

  t->Project("Elt100","fg1.E():fg2.E()","fA.M()<100");
  t->Project("Egt100","fg1.E():fg2.E()","fA.M()>=100");

  Egt100->GetXaxis()->SetTitle("E_{g1} (GeV)");
  Egt100->GetYaxis()->SetTitle("E_{g2} (GeV)");

  Egt100->Draw();
  Elt100->Draw("same");

  TLegend *leg=new TLegend(0.66,0.75,0.89,0.89);
  leg->AddEntry(Elt100,"m_{A}<100 GeV","p");
  leg->AddEntry(Egt100,"m_{A}>100 GeV","p");
  leg->Draw();
}