{
 gStyle->SetOptStat(0);
 gStyle->SetOptTitle(0); 
 gStyle->SetPalette(1);

        TFile gg("temp19_gg.root", "OPEN");
        TFile fg("temp18_fg.root", "OPEN");
        TFile ff("temp17_ff.root", "OPEN");

        TH1F *hPtJets3 = (TH1F*)gg.Get("hPtJets3");
        TH1F *hPtJets5 = (TH1F*)gg.Get("hPtJets5");
        TH1F *hPtJets7 = (TH1F*)gg.Get("hPtJets7");
        TH1F *hPtJets3_eta1 = (TH1F*)gg.Get("hPtJets3_eta1");
        TH1F *hPtJets5_eta1 = (TH1F*)gg.Get("hPtJets5_eta1");
        TH1F *hPtJets7_eta1 = (TH1F*)gg.Get("hPtJets7_eta1");
        TH1F *hPtJets3_eta2 = (TH1F*)gg.Get("hPtJets3_eta2");
        TH1F *hPtJets5_eta2 = (TH1F*)gg.Get("hPtJets5_eta2");
        TH1F *hPtJets7_eta2 = (TH1F*)gg.Get("hPtJets7_eta2");
        TH1F *hPtJets3_eta3 = (TH1F*)gg.Get("hPtJets3_eta3");
        TH1F *hPtJets5_eta3 = (TH1F*)gg.Get("hPtJets5_eta3");
        TH1F *hPtJets7_eta3 = (TH1F*)gg.Get("hPtJets7_eta3");

  hPtJets3 -> Add( (TH1F*)fg.Get("hPtJets3"));
  hPtJets5 -> Add( (TH1F*)fg.Get("hPtJets5"));
  hPtJets7 -> Add( (TH1F*)fg.Get("hPtJets7"));

  hPtJets3 -> Add( (TH1F*)ff.Get("hPtJets3"));
  hPtJets5 -> Add( (TH1F*)ff.Get("hPtJets5"));
  hPtJets7 -> Add( (TH1F*)ff.Get("hPtJets7"));

  hPtJets3_eta1 -> Add( (TH1F*)fg.Get("hPtJets3_eta1"));
  hPtJets5_eta1 -> Add( (TH1F*)fg.Get("hPtJets5_eta1"));
  hPtJets7_eta1 -> Add( (TH1F*)fg.Get("hPtJets7_eta1"));
  hPtJets3_eta1 -> Add( (TH1F*)ff.Get("hPtJets3_eta1"));
  hPtJets5_eta1 -> Add( (TH1F*)ff.Get("hPtJets5_eta1"));
  hPtJets7_eta1 -> Add( (TH1F*)ff.Get("hPtJets7_eta1"));

  hPtJets3_eta2 -> Add( (TH1F*)fg.Get("hPtJets3_eta2"));
  hPtJets5_eta2 -> Add( (TH1F*)fg.Get("hPtJets5_eta2"));
  hPtJets7_eta2 -> Add( (TH1F*)fg.Get("hPtJets7_eta2"));
  hPtJets3_eta2 -> Add( (TH1F*)ff.Get("hPtJets3_eta2"));
  hPtJets5_eta2 -> Add( (TH1F*)ff.Get("hPtJets5_eta2"));
  hPtJets7_eta2 -> Add( (TH1F*)ff.Get("hPtJets7_eta2"));

  hPtJets3_eta3 -> Add( (TH1F*)fg.Get("hPtJets3_eta3"));
  hPtJets5_eta3 -> Add( (TH1F*)fg.Get("hPtJets5_eta3"));
  hPtJets7_eta3 -> Add( (TH1F*)fg.Get("hPtJets7_eta3"));
  hPtJets3_eta3 -> Add( (TH1F*)ff.Get("hPtJets3_eta3"));
  hPtJets5_eta3 -> Add( (TH1F*)ff.Get("hPtJets5_eta3"));
  hPtJets7_eta3 -> Add( (TH1F*)ff.Get("hPtJets7_eta3"));


  hPtJets3 -> SetLineColor(kBlue);
  hPtJets5 -> SetLineColor(kGreen);
  hPtJets7 -> SetLineColor(kRed);
  hPtJets3 -> SetXTitle("pt");
  hPtJets5 -> SetXTitle("pt");
  hPtJets7 -> SetXTitle("pt");
  hPtJets3 -> SetYTitle("Number of jets");
  hPtJets5 -> SetYTitle("Number of jets");
  hPtJets7 -> SetYTitle("Number of jets");

  hPtJets3_eta1 -> SetXTitle("pt");
  hPtJets5_eta1 -> SetXTitle("pt");
  hPtJets7_eta1 -> SetXTitle("pt");
  hPtJets3_eta1 -> SetYTitle("Number of jets");
  hPtJets5_eta1 -> SetYTitle("Number of jets");
  hPtJets7_eta1 -> SetYTitle("Number of jets");

  hPtJets3_eta2 -> SetXTitle("pt");
  hPtJets5_eta2 -> SetXTitle("pt");
  hPtJets7_eta2 -> SetXTitle("pt");
  hPtJets3_eta2 -> SetYTitle("Number of jets");
  hPtJets5_eta2 -> SetYTitle("Number of jets");
  hPtJets7_eta2 -> SetYTitle("Number of jets");
  hPtJets3_eta3 -> SetXTitle("pt");
  hPtJets5_eta3 -> SetXTitle("pt");
  hPtJets7_eta3 -> SetXTitle("pt");
  hPtJets3_eta3 -> SetYTitle("Number of jets");
  hPtJets5_eta3 -> SetYTitle("Number of jets");
  hPtJets7_eta3 -> SetYTitle("Number of jets");
  
  hPtJets3_eta1 -> SetLineColor(kBlue+1);
  hPtJets5_eta1 -> SetLineColor(kGreen+1);
  hPtJets7_eta1 -> SetLineColor(kRed+1);
  hPtJets3_eta2 -> SetLineColor(kBlue+2);
  hPtJets5_eta2 -> SetLineColor(kGreen+2);
  hPtJets7_eta2 -> SetLineColor(kRed+2);

  hPtJets3_eta3 -> SetLineColor(kBlue-1);
  hPtJets5_eta3 -> SetLineColor(kGreen-1);
  hPtJets7_eta3 -> SetLineColor(kRed-1);



  TH1F *ratio35 = new TH1F("ratio35", "ratio37", 50, 0, 150);
  TH1F *ratio37 = new TH1F("ratio37", "ratio37", 50, 0, 150);
  TH1F *ratio57 = new TH1F("ratio57", "ratio57", 50, 0, 150);

  TH1F *ratio35_eta1 = new TH1F("ratio35_eta1", "ratio37_eta1", 50, 0, 150);
  TH1F *ratio37_eta1 = new TH1F("ratio37_eta1", "ratio37_eta1", 50, 0, 150);
  TH1F *ratio57_eta1 = new TH1F("ratio57_eta1", "ratio57_eta1", 50, 0, 150);
  TH1F *ratio35_eta2 = new TH1F("ratio35_eta2", "ratio35_eta2", 50, 0, 150);
  TH1F *ratio37_eta2 = new TH1F("ratio37_eta2", "ratio37_eta2", 50, 0, 150);
  TH1F *ratio57_eta2 = new TH1F("ratio57_eta2", "ratio57_eta2", 50, 0, 150);
  TH1F *ratio35_eta3 = new TH1F("ratio35_eta3", "ratio37_eta3", 50, 0, 150);
  TH1F *ratio37_eta3 = new TH1F("ratio37_eta3", "ratio37_eta3", 50, 0, 150);
  TH1F *ratio57_eta3 = new TH1F("ratio57_eta3", "ratio57_eta3", 50, 0, 150);
 
 ratio35 -> SetLineColor(kBlue);  
 ratio37 -> SetLineColor(kBlue+1);  
 ratio57 -> SetLineColor(kBlue+2);  
 ratio35_eta1 -> SetLineColor(kGreen);  
 ratio37_eta1 -> SetLineColor(kGreen+1);  
 ratio57_eta1 -> SetLineColor(kGreen+2);  
 ratio35_eta2 -> SetLineColor(kRed);  
 ratio37_eta2 -> SetLineColor(kRed+1);  
 ratio57_eta2 -> SetLineColor(kRed+2);  

 ratio35_eta3 -> SetLineColor(kBlue);  
 ratio37_eta3 -> SetLineColor(kBlue+1);  
 ratio57_eta3 -> SetLineColor(kBlue+2);  

 ratio35 -> SetXTitle("pt");  
 ratio37 -> SetXTitle("pt");  
 ratio57 -> SetXTitle("pt");  
 ratio35 -> SetYTitle("ratio");  
 ratio37 -> SetYTitle("ratio");  
 ratio57 -> SetYTitle("ratio");  
 ratio35_eta1 -> SetXTitle("pt");  
 ratio37_eta1 -> SetXTitle("pt");  
 ratio57_eta1 -> SetXTitle("pt");  
 ratio35_eta1 -> SetYTitle("ratio = #frac{#sigma_{zeus}}{#sigma_{cteq}}");  
 ratio37_eta1 -> SetYTitle("ratio = #frac{#sigma_{zeus}}{#sigma_{mrst}}");  
 ratio57_eta1 -> SetYTitle("ratio = #frac{#sigma_{cteq}}{#sigma_{mrst}}");  

 ratio35_eta2 -> SetXTitle("pt");  
 ratio37_eta2 -> SetXTitle("pt");  
 ratio57_eta2 -> SetXTitle("pt");  
 ratio35_eta2 -> SetYTitle("ratio");  
 ratio37_eta2 -> SetYTitle("ratio");  
 ratio57_eta2 -> SetYTitle("ratio");  

 ratio35_eta3 -> SetXTitle("pt");  
 ratio37_eta3 -> SetXTitle("pt");  
 ratio57_eta3 -> SetXTitle("pt");  
 ratio35_eta3 -> SetYTitle("ratio");  
 ratio37_eta3 -> SetYTitle("ratio");  
 ratio57_eta3 -> SetYTitle("ratio");  


  TCanvas *Sum = new TCanvas("Sum", "  ... ", 10, 10, 1200, 400);
  Sum -> Divide(3,1);
  Sum -> SetHighLightColor(45);


  Sum_1 -> SetLogy();
  Sum_2 -> SetLogy();
  Sum_3 -> SetLogy();
  Sum_1 -> SetGrid();
  Sum_2 -> SetGrid();
  Sum_3 -> SetGrid();

  Sum -> cd(1);
  hPtJets3 -> Draw("hist");
  hPtJets3_eta1 -> Draw("hist same");
  hPtJets3_eta2 -> Draw("hist same");
  hPtJets3_eta3 -> Draw("hist same");
  leg = new TLegend(0.5,0.65,0.95,0.95);
  leg->AddEntry(hPtJets3, "all #eta-region", "l");
  leg->AddEntry(hPtJets3_eta1, "0 < |#eta| < 1.3", "l");
  leg->AddEntry(hPtJets3_eta2, "1.3 < |#eta| < 3.0", "l");
  leg->AddEntry(hPtJets3_eta3, "3.0 < |#eta| < 5.0", "l");
  leg->SetHeader("   ZEUS2005");
  leg->Draw();
  
  Sum -> cd(2);
  hPtJets5 -> Draw("hist");
  hPtJets5_eta1 -> Draw("hist same");
  hPtJets5_eta2 -> Draw("hist same");
  hPtJets5_eta3 -> Draw("hist same");
  leg = new TLegend(0.5,0.65,0.95,0.95);
  leg->AddEntry(hPtJets5, "all #eta-region", "l");
  leg->AddEntry(hPtJets5_eta1, "0 < |#eta| < 1.3", "l");
  leg->AddEntry(hPtJets5_eta2, "1.3 < |#eta| < 3.0", "l");
  leg->AddEntry(hPtJets5_eta3, "3.0 < |#eta| < 5.0", "l");
  leg->SetHeader("   CTEQ6M");
  leg->Draw();
  
  Sum -> cd(3);
  hPtJets7 -> Draw("hist");
  hPtJets7_eta1 -> Draw("hist same");
  hPtJets7_eta2 -> Draw("hist same");
  hPtJets7_eta3 -> Draw("hist same");
  leg = new TLegend(0.5,0.65,0.95,0.95);
  leg->AddEntry(hPtJets7, "all #eta-region", "l");
  leg->AddEntry(hPtJets7_eta1, "0 < |#eta| < 1.3", "l");
  leg->AddEntry(hPtJets7_eta2, "1.3 < |#eta| < 3.0", "l");
  leg->AddEntry(hPtJets7_eta3, "3.0 < |#eta| < 5.0", "l");
  leg->SetHeader("  MRST2006");
  leg->Draw();

  Sum -> cd();
  Sum -> SaveAs("~andrey/public/html/pythia/h057a.png");

//////////  RATIOS

  ratio35 ->  Divide(hPtJets3, hPtJets5);
  ratio37 ->  Divide(hPtJets3, hPtJets7);
  ratio57 ->  Divide(hPtJets5, hPtJets7);
  ratio35 -> SetAxisRange(0, 70);
  ratio37 -> SetAxisRange(0, 70);
  ratio57 -> SetAxisRange(0, 70);

  ratio35_eta1 ->  Divide(hPtJets3_eta1, hPtJets5_eta1);
  ratio37_eta1 ->  Divide(hPtJets3_eta1, hPtJets7_eta1);
  ratio57_eta1 ->  Divide(hPtJets5_eta1, hPtJets7_eta1);
  ratio35_eta1 -> SetAxisRange(0, 100);
  ratio37_eta1 -> SetAxisRange(0, 100);
  ratio57_eta1 -> SetAxisRange(0, 100);

  ratio35_eta2 ->  Divide(hPtJets3_eta2, hPtJets5_eta2);
  ratio37_eta2 ->  Divide(hPtJets3_eta2, hPtJets7_eta2);
  ratio57_eta2 ->  Divide(hPtJets5_eta2, hPtJets7_eta2);
  ratio35_eta2 -> SetAxisRange(0, 100);
  ratio37_eta2 -> SetAxisRange(0, 100);
  ratio57_eta2 -> SetAxisRange(0, 100);

  ratio35_eta3 ->  Divide(hPtJets3_eta3, hPtJets5_eta3);
  ratio37_eta3 ->  Divide(hPtJets3_eta3, hPtJets7_eta3);
  ratio57_eta3 ->  Divide(hPtJets5_eta3, hPtJets7_eta3);
  ratio35_eta3 -> SetAxisRange(0, 100);
  ratio37_eta3 -> SetAxisRange(0, 100);
  ratio57_eta3 -> SetAxisRange(0, 100);

  ratio35_eta1 -> SetMinimum(0.4);
  ratio37_eta1 -> SetMinimum(0.4);
  ratio57_eta1 -> SetMinimum(0.4);
  ratio35_eta1 -> SetMaximum(1.5);
  ratio37_eta1 -> SetMaximum(1.5);
  ratio57_eta1 -> SetMaximum(1.5);

  Sum_1 -> SetLogy(0);
  Sum_2 -> SetLogy(0);
  Sum_3 -> SetLogy(0);


//////////////////

    TF1 *f1 = new TF1("f1","[0]+[1]*x+[2]*x^2", 10, 100);
    TF1 *f2 = new TF1("f2","[0]+[1]*x+[2]*x^2", 10, 100);
    TF1 *f3 = new TF1("f3","[0]+[1]*x+[2]*log(x)", 10, 100);
//    f1 -> SetLineWidth(1);
//    f2 -> SetLineWidth(1);
//    f3 -> SetLineWidth(1);
    f1 -> SetLineColor(kGreen);
    f2 -> SetLineColor(kRed);
    f3 -> SetLineColor(kBlue);

  Sum -> cd(1);
  ratio35_eta1 -> Fit("f1","WR");
  ratio35_eta2 -> Fit("f2","WR+", "same");
  ratio35_eta3 -> Fit("f3","WR+","same");
//  ratio35 -> Draw("same");
  leg = new TLegend(0.12,0.8,0.5,0.97);
  leg->AddEntry(ratio35_eta1, "0 < |#eta| < 1.3", "l");
  leg->AddEntry(ratio35_eta2, "1.3 < |#eta| < 3.0", "l");
  leg->AddEntry(ratio35_eta3, "3.0 < |#eta| < 5.0", "l");
  leg->SetHeader(" ZEUS/CTEQ");
  leg->Draw();
  Sum -> cd(2);
//  ratio37_eta1 -> Draw("e1p");
//  ratio37_eta2 -> Draw("same");
//  ratio37_eta3 -> Draw("same");
  ratio37_eta1 -> Fit("f1","WR");
  ratio37_eta2 -> Fit("f2","WR+","same");
  ratio37_eta3 -> Fit("f3","WR+","same");

//  ratio37 -> Draw("same");
  leg = new TLegend(0.12,0.8,0.5,0.97);
  leg->AddEntry(ratio37_eta1, "0 < |#eta| < 1.3", "l");
  leg->AddEntry(ratio37_eta2, "1.3 < |#eta| < 3.0", "l");
  leg->AddEntry(ratio37_eta3, "3.0 < |#eta| < 5.0", "l");
  leg->SetHeader("  ZEUS/MRST");
  leg->Draw();
  Sum -> cd(3);
//  ratio57_eta1 -> Draw("e1p");
//  ratio57_eta2 -> Draw("same");
//  ratio57_eta3 -> Draw("same");
  ratio57_eta1 -> Fit("f1","WR");
  ratio57_eta2 -> Fit("f2","WR","same");
  ratio57_eta3 -> Fit("f3","WR","same");

//  ratio57 -> Draw("same");
  leg = new TLegend(0.12,0.8,0.5,0.97);
  leg->AddEntry(ratio57_eta1, "0 < |#eta| < 1.3", "l");
  leg->AddEntry(ratio57_eta2, "1.3 < |#eta| < 3.0", "l");
  leg->AddEntry(ratio57_eta3, "3.0 < |#eta| < 5.0", "l");
  leg->SetHeader("  CTEQ/MRST");
  leg->Draw();
  Sum -> SaveAs("~andrey/public/html/pythia/h058_ratio.png");

}