//Scaled Histogrammakes

{
  //gROOT->LoadMacro("~/Desktop/atlasstyle/AtlasStyle.C");
  //SetAtlasStyle();
  TFile* file_sig = new TFile("inputs/susymcoutput.root","OPEN");
  TFile* file_bkg = new TFile("inputs/ttbarmcoutput.root","OPEN");

  TCanvas* canvas = new TCanvas("canvas", "canvas", 10, 10, 500, 500);


  double Lumi_data = 10; //fb-1
  double Lumi_mc_s = 1.04649; //fb-1 signal
  double Lumi_mc_b = 6.84981; //fb-1 background

  //Make some scaled histograms

  double scaleFactor_s = Lumi_data/Lumi_mc_s;
  double scaleFactor_b = Lumi_data/Lumi_mc_b;
  TH1D* emiss_s = (TH1D*)file_sig->Get("et_miss");
  TH1D* emiss_b = (TH1D*)file_bkg->Get("et_miss");
  emiss_s.SetLineStyle(2);
  emiss_s.SetLineColor(2);
  //Scale the histograms
  //emiss_s.Sumw2();
  //emiss_b.Sumw2();
  emiss_s.Scale(scaleFactor_s);
  emiss_b.Scale(scaleFactor_b);
  THStack emiss_stack("emiss_stack","Emiss Comparison");
  
  emiss_stack.Add(emiss_s);
  emiss_stack.Add(emiss_b);
 
  TLegend* leg = new TLegend(0.5, 0.75, 0.9, 0.9);
  leg->AddEntry(emiss_s, "signal", "l");  // l for line, f for box, p for point
  leg->AddEntry(emiss_b, "background", "l");
     
  emiss_stack.Draw("nostackhist");

  emiss_stack->GetXaxis()->SetTitle("Emiss_{T}");
  emiss_stack->GetYaxis()->SetTitle("Number of scaled events");

  leg->Draw(); 

  canvas->Print("./outputs/stack_emiss.ps");
}