#include "TFile.h"
#include "TTree.h"
#include "TH1F.h"


void teste(){


 TH1::SetDefaultSumw2(true);

 TFile *MC = new TFile("mc_weighted.root");
 TFile *Dados= new TFile("NtupleWithSweights.root");
 TTree *Dados_Tree = (TTree*)Dados->Get("sTree");
 TTree *MC_Tree = (TTree*)MC->Get("weightedTree");


 TH1F *Histo_Dados_D_IPCHI2_OWNPV = new TH1F("Histo_Dados_D_IPCHI2_OWNPV ","D_IPCHI2_OWNPV",100,0,12);
 TH1F *Histo_MC_D_IPCHI2_OWNPV = new TH1F("Histo_MC_D_IPCHI2_OWNPV","D_IPCHI2_OWNPV",100,0,12);
 TH1F *Histo_MC_CIN_D_IPCHI2_OWNPV = new TH1F("Histo_MC_CIN_D_IPCHI2_OWNPV","D_IPCHI2_OWNPV",100,0,12);

 TCanvas *Canvas_D_IPCHI2_OWNPV = new TCanvas("Canvas_D_IPCHI2_OWNPV","D_IPCHI2_OWNPV",700,400);
 Canvas_D_IPCHI2_OWNPV->Divide(2);
 Canvas_D_IPCHI2_OWNPV->cd(1);

 Dados_Tree->Draw("D_IPCHI2_OWNPV>>Histo_Dados_D_IPCHI2_OWNPV","((D_IPCHI2_OWNPV>0)&&(D_IPCHI2_OWNPV<12))*(signalweight)");
 MC_Tree->Draw("D_IPCHI2_OWNPV>>Histo_MC_D_IPCHI2_OWNPV","((D_IPCHI2_OWNPV>0)&&(D_IPCHI2_OWNPV<12))*(D_PIDEff)*(wdata)*(wnTracks)");
 MC_Tree->Draw("D_IPCHI2_OWNPV>>Histo_MC_CIN_D_IPCHI2_OWNPV","((D_IPCHI2_OWNPV>0)&&(D_IPCHI2_OWNPV<12))*(D_PIDEff)*(wdata)*(wnTracks)*(weight)");

 Histo_MC_D_IPCHI2_OWNPV->Scale(1./(Histo_MC_D_IPCHI2_OWNPV->Integral()));
 Histo_MC_CIN_D_IPCHI2_OWNPV->Scale(1./(Histo_MC_CIN_D_IPCHI2_OWNPV->Integral()));
 Histo_Dados_D_IPCHI2_OWNPV->Scale(1./(Histo_Dados_D_IPCHI2_OWNPV->Integral())); 
 
 Histo_Dados_D_IPCHI2_OWNPV->Draw("P");
 Histo_MC_D_IPCHI2_OWNPV->Draw("sameHIST");
 Histo_MC_CIN_D_IPCHI2_OWNPV->Draw("sameHIST");

 Canvas_D_IPCHI2_OWNPV->cd(2);

 TH1F *Histo_Dados_D_IPCHI2_OWNPV_copy = (TH1F*)Histo_Dados_D_IPCHI2_OWNPV->Clone("");
 TH1F *Histo_MC_CIN_D_IPCHI2_OWNPV_copy = (TH1F*)Histo_MC_CIN_D_IPCHI2_OWNPV->Clone("");

 Histo_MC_D_IPCHI2_OWNPV->SetLineColor(kBlue);
 Histo_MC_CIN_D_IPCHI2_OWNPV->SetLineColor(kRed);
 

 Histo_Dados_D_IPCHI2_OWNPV->Divide(Histo_MC_CIN_D_IPCHI2_OWNPV_copy);

 Histo_Dados_D_IPCHI2_OWNPV->Draw("E");

//#####################################################################################################################
/*
TH1F *Histo_Dados_D_IP_OWNPV = new TH1F("Histo_Dados_D_IP_OWNPV ","D_IP_OWNPV",100,0,0.07);
 TH1F *Histo_MC_D_IP_OWNPV = new TH1F("Histo_MC_D_IP_OWNPV","D_IP_OWNPV",100,0,0.07);
 TH1F *Histo_MC_CIN_D_IP_OWNPV = new TH1F("Histo_MC_CIN_D_IP_OWNPV","D_IP_OWNPV",100,0,0.07);

 TCanvas *Canvas_D_IP_OWNPV = new TCanvas("Canvas_D_IP_OWNPV","D_IP_OWNPV",700,400);
 Canvas_D_IP_OWNPV->Divide(2);
 Canvas_D_IP_OWNPV->cd(1);

 MC_Tree->Draw("D_IP_OWNPV>>Histo_MC_D_IP_OWNPV","((D_IP_OWNPV>0)&&(D_IP_OWNPV<0.07))*(D_PIDEff)*(wdata)*(wnTracks)");
 MC_Tree->Draw("D_IP_OWNPV>>Histo_MC_CIN_D_IP_OWNPV","((D_IP_OWNPV>0)&&(D_IP_OWNPV<0.07))*(D_PIDEff)*(wdata)*(wnTracks)*(weight)");
 Dados_Tree->Draw("D_IP_OWNPV>>Histo_Dados_D_IP_OWNPV","((D_IP_OWNPV>0)&&(D_IP_OWNPV<0.07))*(signalweight)");

 
 Histo_MC_D_IP_OWNPV->Scale(1./(Histo_MC_D_IP_OWNPV->Integral()));
 Histo_MC_CIN_D_IP_OWNPV->Scale(1./(Histo_MC_CIN_D_IP_OWNPV->Integral()));
 Histo_Dados_D_IP_OWNPV->Scale(1./(Histo_Dados_D_IP_OWNPV->Integral())); 
 
 Histo_Dados_D_IP_OWNPV->Draw("P");
 Histo_MC_D_IP_OWNPV->Draw("sameHIST");
 Histo_MC_CIN_D_IP_OWNPV->Draw("sameHIST");

 Canvas_D_IP_OWNPV->cd(2);

 TH1F *Histo_Dados_D_IP_OWNPV_copy = (TH1F*)Histo_Dados_D_IP_OWNPV->Clone("");
 TH1F *Histo_MC_CIN_D_IP_OWNPV_copy = (TH1F*)Histo_MC_CIN_D_IP_OWNPV->Clone("");

 Histo_MC_D_IP_OWNPV->SetLineColor(kBlue);
 Histo_MC_CIN_D_IP_OWNPV->SetLineColor(kRed);
 

 Histo_Dados_D_IP_OWNPV->Divide(Histo_MC_CIN_D_IP_OWNPV_copy);

 Histo_Dados_D_IP_OWNPV->Draw("E");
*/
}