{

TFile *F = new TFile("4lepMassPseudoNoPhzz.root");
TFile *f = new TFile("4lepMassPseudozz.root");

TFile *G = new TFile("4lepMassNCNoPhzz.root");
TFile *g = new TFile("4lepMassNCzz.root");

TFile *z = new TFile("NewFit4lepMassNCzz.root");
TFile *Z = new TFile("NewFit4lepMassPseudozz.root");

TH1F *NChist = z->Get("Fit Stats");
TH1F *Pseudhist = Z->Get("Fit Stats");

TF1 *NCfunc = NChist->GetFunction("fitf");

TF1 *Pseudfunc = Pseudhist->GetFunction("fitf");

TCanvas *c1 = new TCanvas("c1","Ratio",1000,800);

TH1F *main2 = (TH1F*)F->Get("main117");
TH1F *main = (TH1F*)f->Get("main113");

TH1F *main3 = (TH1F*)G->Get("main125");
TH1F *main4 = (TH1F*)g->Get("main121");

added1 = new TH1F("added1","4 lep Mass Pseudo ZZ",200,0,600);
added2 = new TH1F("added2","4 lep Mass NC ZZ",200,0,600);


added1->Add(main);
added1->Add(main2);

added2->Add(main3);
added2->Add(main4);

ratio = new TH1F("efficiency","ZZ 4 lep Mass Cut Efficiency (14TeV)",200,0,600);

ratio->Divide(added1,added2);

ratio->GetXaxis()->SetTitle("Mass (GeV)");
ratio->GetYaxis()->SetTitle("Efficiency of pT/Eta cuts");
ratio->GetYaxis()->SetTitleOffset(1.3);

TFile *g = new TFile("4lepMassNCzz.root");
TFile *e = new TFile("4lepMassNC+tauzz.root");

TF1 *func2 = new TF1("func2","((.5+.5*TMath::Erf((x-[0])/[1]))*([3]/(1+TMath::Power(TMath::E(),(x-[0])/[2])))+(.5+.5*TMath::Erf((x-[4])/[5]))*(([7]/(1+TMath::Power(TMath::E(),(x-[4])/[6]))+([9]/(1+TMath::Power(TMath::E(),(x-[4])/[8]))))))/((.5+.5*TMath::Erf((x-[10])/[11]))*([13]/(1+TMath::Power(TMath::E(),(x-[10])/[12])))+(.5+.5*TMath::Erf((x-[14])/[15]))*(([17]/(1+TMath::Power(TMath::E(),(x-[14])/[16]))+([19]/(1+TMath::Power(TMath::E(),(x-[14])/[18]))))))",85,600);

func2->SetParameter(0,Pseudfunc->GetParameter(0));
func2->SetParameter(1,Pseudfunc->GetParameter(1));
func2->SetParameter(2,Pseudfunc->GetParameter(2));
func2->SetParameter(3,Pseudfunc->GetParameter(3));
func2->SetParameter(4,Pseudfunc->GetParameter(4));
func2->SetParameter(5,Pseudfunc->GetParameter(5));
func2->SetParameter(6,Pseudfunc->GetParameter(6));
func2->SetParameter(7,Pseudfunc->GetParameter(7));
func2->SetParameter(8,Pseudfunc->GetParameter(8));
func2->SetParameter(9,Pseudfunc->GetParameter(9));
func2->SetParameter(10,NCfunc->GetParameter(0));
func2->SetParameter(11,NCfunc->GetParameter(1));
func2->SetParameter(12,NCfunc->GetParameter(2));
func2->SetParameter(13,NCfunc->GetParameter(3));
func2->SetParameter(14,NCfunc->GetParameter(4));
func2->SetParameter(15,NCfunc->GetParameter(5));
func2->SetParameter(16,NCfunc->GetParameter(6));
func2->SetParameter(17,NCfunc->GetParameter(7));
func2->SetParameter(18,NCfunc->GetParameter(8));
func2->SetParameter(19,NCfunc->GetParameter(9));


gStyle->SetOptFit(1);
gStyle->SetOptStat("ne");

ratio->Draw();
func2->Draw("same");
c1->Update();
/*
c1->SaveAs("/home/msnowball/zPrograms/Fits14/4lepMassEfficiency+tauzz.jpg");

TCanvas *c2 = new TCanvas("c2","Ratio",454,347);
ratio->SetMarkerSize(.15);
ratio->Draw("HIST C");
func2->Draw("same");
c2->Update();
c2->SaveAs("/home/msnowball/zPrograms/Fits14/4lepMassEfficiency+tauzzSm.jpg");
*/
}