{

//OPEN ROOT FILE WITH MC DATA
TFile *fMC = new TFile("D:\\\EEE\\ICD\\2019\\ROOT\\ALTA-01from2018-10-02to2018-10-03-MC.root");
TTree *eeeMC = (TTree*)fMC->Get("eee");

//CREATE HISTOS FOR MC THETA AND PHI DISTRIBUTIONS
TH1F *hPhiMC = new TH1F("hPhiMC", "MC Phi Distribution", 72, 0, 360);
TH1F *hThetaMC = new TH1F("hThetaMC", "MC Theta Distribution", 45, 0, 90);

//FILL AND DRAW MC PHI DISTRIBUTION
TCanvas *cPhiMC = new TCanvas();
eeeMC->Draw("(Phi+180) >> hPhiMC");
hPhiMC->Draw();
cPhiMC->SaveAs("D:\\\EEE\\ICD\\2019\\ROOT\\MCPhi.png");

//FILL AND DRAW MC THETA DISTRIBUTION
TCanvas *cThetaMC = new TCanvas();
eeeMC->Draw("Theta >> hThetaMC");
hThetaMC->Draw();
cThetaMC->SaveAs("D:\\\EEE\\ICD\\2019\\ROOT\\MCTheta.png");


//OPEN ROOT FILE WITH EXPERIMENTAL DATA
TFile *fDATA = new TFile("D:\\\EEE\\ICD\\2019\\ROOT\\ALTA-01from2018-10-02to2018-10-03.root");
TTree *eeeDATA = (TTree*)fDATA->Get("eee");

//CREATE HISTOS FOR EXPERIMENTAL THETA AND PHI DISTRIBUTIONS
TH1F *hPhiDATA = new TH1F("hPhiDATA", "Experimental Phi Distribution", 72, 0, 360);
TH1F *hThetaDATA = new TH1F("hThetaDATA", "Experimental Theta Distribution", 45, 0, 90);

//FILL AND DRAW EXPERIMENTAL PHI DISTRIBUTION
TCanvas *cPhiDATA = new TCanvas();
eeeDATA->Draw("(Phi+180) >> hPhiDATA");
hPhiDATA->Draw();
cPhiDATA->SaveAs("D:\\\EEE\\ICD\\2019\\ROOT\\DATAPhi.png");

//FILL AND DRAW EXPERIMENTAL THETA DISTRIBUTION
TCanvas *cThetaDATA = new TCanvas();
eeeDATA->Draw("Theta >> hThetaDATA");
hThetaDATA->Draw();
cThetaDATA->SaveAs("D:\\\EEE\\ICD\\2019\\ROOT\\DATATheta.png");


// CREATE HISTOS FOR THE RATIOS
TH1F *hPhiRATIO = new TH1F("hPhiRATIO", "Ratio between exp and MC Phi Distribution", 72, 0, 360);
TH1F *hThetaRATIO = new TH1F("hThetaRATIO", "Ratio between exp and MC Theta Distribution", 45, 0, 90);

hPhiRATIO->Divide(hPhiDATA,hPhiMC);
TCanvas *cPhiRATIO = new TCanvas();
hPhiRATIO->Draw();
cPhiRATIO->SaveAs("D:\\\EEE\\ICD\\2019\\ROOT\\RATIOPhi.png");

hThetaRATIO->Divide(hThetaDATA,hThetaMC);
TCanvas *cThetaRATIO = new TCanvas();
hThetaRATIO->Draw();
cThetaRATIO->SaveAs("D:\\\EEE\\ICD\\2019\\ROOT\\RATIOTheta.png");


}