{
#include "TStyle.h"

    
    gROOT->Reset();
    gROOT->SetStyle("Plain");
    gStyle->SetOptStat(kFALSE);
    
    // Create a histograms and stack
    THStack hs1("hs1","");
    TH1F* JetPt_50PU_TopJets;

    
    TFile *bkg1 = TFile::Open("SingleS_P+MET_50PU_TopJets_mt2w_his.root");
    JetPt_50PU_TopJets= (TH1F*) bkg1->Get("hJetPt");


    
    
TFile* file = new TFile ("stacked_plots_300/stack_JetPt_50PU.root","RECREATE");
    TCanvas* c1 = new TCanvas("c1","stack_JetPt",800,20,1000,800) ;
    c1->cd(1);
    c1->SetLogy();
    file->Write();
	
	hs1.Draw();    
    
	TLatex T1;
	T1.DrawLatexNDC(.05,.95, "Delphes samples #sqrt{s}=14 TeV, #int Ldt=300 fb^{-1}");


    JetPt_50PU_TopJets->SetFillColor(kAzure+1);
    hs1->Add(JetPt_50PU_TopJets);



    

    hs1->GetXaxis()->SetRangeUser(0, 600);
    hs1->GetHistogram()->GetXaxis()->SetTitle("P_{T} [GeV]");
    hs1->GetHistogram()->GetYaxis()->SetTitle("Events/20 GeV");
    
    TLegend *pl = new TLegend(0.7, 0.6, 0.9, 0.8);
    pl->SetTextSize(0.02);
    pl->SetFillColor(0);
    pl->AddEntry(JetPt_50PU_TopJets, "TopJets JetPt 50PU",  "lpf");
    pl->SetFillColor(0);
    pl->SetLineColor(2);
    pl->Draw("same");
    c1->SaveAs("stacked_plots_300/stack_50PU_JetPt.png");
    c1->SaveAs("stacked_plots_300/stack_50PU_JetPt.pdf");
    c1->SaveAs("stacked_plots_300/stack_50PU_JetPt.root");  

}