#include "TFile.h"
#include "TDirectory.h"
#include "TH1.h"
#include "TCanvas.h"
void checksum(){
  TFile *f=new TFile("fitLRWW160.root");
  TDirectory * dir=dynamic_cast<TDirectory *>(f->Get("Templates"));
  dir->ls();
  TH1* h=dynamic_cast<TH1*>(dir->Get("All_Higgs"));
  TH1* h1=dynamic_cast<TH1*>(dir->Get("All_WH"));
  TH1* h2=dynamic_cast<TH1*>(dir->Get("All_ZH"));
  TH1* h3=dynamic_cast<TH1*>(dir->Get("All_VBF"));
  TH1* h4=dynamic_cast<TH1*>(h->Clone("Sum"));
  h4->Add(h1);
  h4->Add(h2);
  h4->Add(h3);
  TCanvas* c=new TCanvas();
  c->Divide(2,3);
  int i=1;
  c->cd(i++);
  h->Draw();
  c->cd(i++);
  h1->Draw();
  c->cd(i++);
  h2->Draw();
  c->cd(i++);
  h3->Draw();
  c->cd(i++);
  h4->Draw();
  h4->SetLineColor(2);
  TFile *f2=new TFile("../../run21030/templatesBl97_v100/fit/out_1/fitLRWW160.root");
  TDirectory * dir2=dynamic_cast<TDirectory *>(f2->Get("Templates"));
  TH1* g=dynamic_cast<TH1*>(dir2->Get("HiSB_Higgs"));
  TH1* g2=dynamic_cast<TH1*>(dir2->Get("LoSB_Higgs"));
  g2->Add(g);
  TCanvas* c2=new TCanvas();
  c2->Divide(2,2);
  i=1;
  c2->cd(i++);
  g2->SetLineColor(3);
  g2->SetMarkerSize(1);
  g2->Draw();
  c2->cd(i++);
  h4->Draw();
  c2->cd(i++);
  g2->Draw();
  h4->Draw("SAME*");
    
}