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

        TFile f("temp5.root", "OPEN");
        TH1F *hPtJets3 = (TH1F*)f.Get("hPtJets3");
        TH1F *hPtJets5 = (TH1F*)f.Get("hPtJets5");
        TH1F *hPtJets7 = (TH1F*)f.Get("hPtJets7");
        TH1F *hPtJets3_eta1 = (TH1F*)f.Get("hPtJets3_eta1");
        TH1F *hPtJets5_eta1 = (TH1F*)f.Get("hPtJets5_eta1");
        TH1F *hPtJets7_eta1 = (TH1F*)f.Get("hPtJets7_eta1");
        TH1F *hPtJets3_eta2 = (TH1F*)f.Get("hPtJets3_eta2");
        TH1F *hPtJets5_eta2 = (TH1F*)f.Get("hPtJets5_eta2");
/*
  hPtJets3_eta1 -> Sumw2();
  hPtJets5_eta1 -> Sumw2();
  hPtJets7_eta1 -> Sumw2();
  hPtJets3_eta2 -> Sumw2();
  hPtJets5_eta2 -> Sumw2();
  hPtJets7_eta2 -> Sumw2();
*/
  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_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);


        TH1F *x1_3 = (TH1F*)f.Get("x1_3");
        TH1F *x1_5 = (TH1F*)f.Get("x1_5");
        TH1F *x1_7 = (TH1F*)f.Get("x1_7");
        TH1F *x2_3 = (TH1F*)f.Get("x2_3");
        TH1F *x2_5 = (TH1F*)f.Get("x2_5");
        TH1F *x2_7 = (TH1F*)f.Get("x2_7");
        TH1F *xf_3 = (TH1F*)f.Get("xf_3");
        TH1F *xf_5 = (TH1F*)f.Get("xf_5");
        TH1F *xf_7 = (TH1F*)f.Get("xf_7");
        TH1F *Q2_3 = (TH1F*)f.Get("Q2_3");
        TH1F *Q2_5 = (TH1F*)f.Get("Q2_5");
        TH1F *Q2_7 = (TH1F*)f.Get("Q2_7");
/*
        TH1F *x1_7_eta1 = (TH1F*)f.Get("x1_7_eta1");
        TH1F *x1_7_eta2 = (TH1F*)f.Get("x1_7_eta2");
        TH1F *x1_5_eta1 = (TH1F*)f.Get("x1_5_eta1");
        TH1F *x1_5_eta2 = (TH1F*)f.Get("x1_5_eta2");
        TH1F *x1_3_eta1 = (TH1F*)f.Get("x1_3_eta1");
        TH1F *x1_3_eta2 = (TH1F*)f.Get("x1_3_eta2");
*/
  x1_3 -> SetLineColor(kBlue);
  x2_3 -> SetLineColor(kBlue+1);
  xf_3 -> SetLineColor(kBlue+2);
  Q2_3 -> SetLineColor(kBlue+3);
  x1_5 -> SetLineColor(kGreen);
  x2_5 -> SetLineColor(kGreen+1);
  xf_5 -> SetLineColor(kGreen+2);
  Q2_5 -> SetLineColor(kGreen+3);
  x1_7 -> SetLineColor(kRed);
  x2_7 -> SetLineColor(kRed+1);
  xf_7 -> SetLineColor(kRed+2);
  Q2_7 -> SetLineColor(kRed+3);
  x1_3 -> SetXTitle("log10 (x)");
  x1_5 -> SetXTitle("log10 (x)");
  x1_7 -> SetXTitle("log10 (x)");
  x2_3 -> SetXTitle("log10 (x)");
  x2_5 -> SetXTitle("log10 (x)");
  x2_7 -> SetXTitle("log10 (x)");
  xf_3 -> SetXTitle("log10 (x)");
  xf_5 -> SetXTitle("log10 (x)");
  xf_7 -> SetXTitle("log10 (x)");
  Q2_3 -> SetXTitle("GeV");
  Q2_5 -> SetXTitle("GeV");
  Q2_7 -> SetXTitle("GeV");
/*
  x1_3_eta1 -> SetLineColor(kBlue+1);
  x1_3_eta2 -> SetLineColor(kBlue+2);
  x1_5_eta1 -> SetLineColor(kBlue+1);
  x1_5_eta2 -> SetLineColor(kBlue+2);
  x1_7_eta1 -> SetLineColor(kBlue+1);
  x1_7_eta2 -> SetLineColor(kBlue+2);
*/
//        TH1F *x1Fx_Q100_3 = (TH1F*)f.Get("x1Fx_Q100_3");
//        TH1F *x2Fx_Q100_3 = (TH1F*)f.Get("x2Fx_Q100_3");
//        TH1F *x1Fx_Q100_5 = (TH1F*)f.Get("x1Fx_Q100_5");
//        TH1F *x2Fx_Q100_5 = (TH1F*)f.Get("x2Fx_Q100_5");
//        TH1F *x1Fx_Q100_7 = (TH1F*)f.Get("x1Fx_Q100_7");
//        TH1F *x2Fx_Q100_7 = (TH1F*)f.Get("x2Fx_Q100_7");
//  x1Fx_Q100_3 -> SetYTitle("x*F(x)");
//  x1Fx_Q100_5 -> SetYTitle("x*F(x)");
//  x1Fx_Q100_7 -> SetYTitle("x*F(x)");

  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);
 
 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 -> 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");  
 ratio37_eta1 -> SetYTitle("ratio");  
 ratio57_eta1 -> SetYTitle("ratio");  

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

  TCanvas *Sum = new TCanvas("Sum", "CTEQ3D, GRV94M, CTEQ5L difference", 10, 10, 1200, 400);
  Sum -> Divide(3,1);
  Sum -> SetHighLightColor(45);
/*
  Sum -> cd(1);
  x1Fx_Q100_3 -> Draw("e1p");
//  x2Fx_Q100_3 -> Draw("same");
  leg = new TLegend(0.75,0.85,0.95,0.95);
  leg->AddEntry(x1Fx_Q100_3, "x1*F(x1)", "l");
//  leg->AddEntry(x2Fx_Q100_3, "x2*F(x2)", "l");
  leg->SetHeader("x*F(x)");
  leg->Draw();
  Sum -> cd(2);
  x1Fx_Q100_5 -> Draw("e1p");
//  x2Fx_Q100_5 -> Draw("same");
  Sum -> cd(3);
  x1Fx_Q100_7 -> Draw("e1p");
//  x2Fx_Q100_7 -> Draw("same");
  Sum -> SaveAs("~andrey/public/html/pythia/h000_xFx.png");
*/

Double_t cs3, cs5, cs7;
  Double_t cs_c3d[30] = {67.29,67.31,67.41,67.35,67.32,67.34,67.42,67.37,67.30,67.30,
                         67.31,67.35,67.37,67.32,67.31,67.43,67.36,67.28,67.28,67.35,
                         67.38,67.31,67.38,67.29,67.30,67.42,67.27,67.47,67.30,67.38};
  Double_t cs_grv[30] = {80.60,80.34,80.36,80.30,80.34,80.34,80.40,80.31,80.33,80.24,
                         80.35,80.39,80.17,80.44,80.14,80.25,80.36,80.45,80.44,80.25,
                         80.50,80.32,80.28,80.57,80.39,80.40,80.57,80.26,80.35,80.33};
  Double_t cs_c5l[30] = {89.34,89.42,89.21,89.31,89.20,89.04,89.39,89.31,89.28,89.18,
                         89.21,89.13,89.24,89.40,89.24,89.52,89.24,89.25,89.44,89.35,
                         89.37,89.16,89.26,89.00,89.14,89.40,89.29,89.21,89.41,89.28};
for(Int_t i =0; i<30; i++)
{ 
cs3 +=cs_c3d[i];
cs5 +=cs_grv[i];
cs7 +=cs_c5l[i];
}
cs3 = cs3/30;
cs5 = cs5/30;
cs7 = cs7/30;

cout<<"cteq3d: "<<cs3<<"  grv94: "<<cs5<<" cteq5l: "<<cs7<<endl;
cout<<"cteq3d: "<<log10(cs3)<<"  grv94: "<<log10(cs5)<<" cteq5l: "<<log10(cs7)<<endl;

  //ftemp.Close();
//  Double_t xbins[14] = {0.0, 0.001, 0.002, 0.003, 0.005, 0.01, 0.02, 0.03, 0.05, 0.1, 0.2, 0.3, 0.5, 1.0};
//Double_t abins[50];


  Sum_1 -> SetLogy();
  Sum_2 -> SetLogy();
  Sum_3 -> SetLogy();
  Sum -> cd(1);
  hPtJets3 -> Draw("hist");
  hPtJets3_eta1 -> Draw("hist same");
  hPtJets3_eta2 -> 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, "2 < |#eta| < 4", "l");
  leg->AddEntry(hPtJets3_eta2, "4 < |#eta| < 5", "l");
  leg->SetHeader("   CTEQ3D");
  leg->Draw();
  
  Sum -> cd(2);
  hPtJets5 -> Draw("hist");
  hPtJets5_eta1 -> Draw("hist same");
  hPtJets5_eta2 -> 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, "2 < |#eta| < 4", "l");
  leg->AddEntry(hPtJets5_eta2, "4 < |#eta| < 5", "l");
  leg->SetHeader("    GRV94M");
  leg->Draw();
  
  Sum -> cd(3);
  hPtJets7 -> Draw("hist");
  hPtJets7_eta1 -> Draw("hist same");
  hPtJets7_eta2 -> 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, "2 < |#eta| < 4", "l");
  leg->AddEntry(hPtJets7_eta2, "4 < |#eta| < 5", "l");
  leg->SetHeader("   CTEQ5L");
  leg->Draw();

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

  Sum -> cd(1);
  hPtJets3 -> Draw("hist");
  hPtJets5 -> Draw("hist same");
  hPtJets7 -> Draw("hist same");
  leg = new TLegend(0.5,0.65,0.95,0.95);
  leg->AddEntry(hPtJets3, "CTEQ3D", "l");
  leg->AddEntry(hPtJets5, "GRV94M", "l");
  leg->AddEntry(hPtJets7, "CTEQ5L", "l");
  leg->SetHeader("  All #eta-region");
  leg->Draw();
  Sum -> cd(2);
  hPtJets3_eta1 -> Draw("hist");
  hPtJets5_eta1 -> Draw("hist same");
  hPtJets7_eta1 -> Draw("hist same");
  leg = new TLegend(0.5,0.65,0.95,0.95);
  leg->AddEntry(hPtJets3_eta1, "CTEQ3D", "l");
  leg->AddEntry(hPtJets5_eta1, "GRV94M", "l");
  leg->AddEntry(hPtJets7_eta1, "CTEQ5L", "l");
  leg->SetHeader("  2 < |#eta| <4");
  leg->Draw();
  Sum -> cd(3);
  hPtJets3_eta2 -> Draw("hist");
  hPtJets5_eta2 -> Draw("hist same");
  hPtJets7_eta2 -> Draw("hist same");
  leg = new TLegend(0.5,0.65,0.95,0.95);
  leg->AddEntry(hPtJets3_eta2, "CTEQ3D", "l");
  leg->AddEntry(hPtJets5_eta2, "GRV94M", "l");
  leg->AddEntry(hPtJets7_eta2, "CTEQ5L", "l");
  leg->SetHeader("  4 < |#eta| < 5");
  leg->Draw();
  
  Sum -> cd();
  Sum -> SaveAs("~andrey/public/html/pythia/h033.png");

  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, 70);
  ratio37_eta1 -> SetAxisRange(0, 70);
  ratio57_eta1 -> SetAxisRange(0, 70);

  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, 70);
  ratio37_eta2 -> SetAxisRange(0, 70);
  ratio57_eta2 -> SetAxisRange(0, 70);

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

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

  Sum -> cd(1);
  ratio35_eta1 -> Draw("e1p");
  ratio35_eta2 -> Draw("same");
//  ratio35 -> Draw("same");
  leg = new TLegend(0.12,0.83,0.5,0.97);
  leg->AddEntry(ratio35_eta1, "2 < |#eta| < 4", "l");
  leg->AddEntry(ratio35_eta2, "4 < |#eta| < 5", "l");
  leg->SetHeader("  CTEQ3D/GRV94M");
  leg->Draw();
  Sum -> cd(2);
  ratio37_eta1 -> Draw("e1p");
  ratio37_eta2 -> Draw("same");
//  ratio37 -> Draw("same");
  leg = new TLegend(0.12,0.83,0.5,0.97);
  leg->AddEntry(ratio37_eta1, "2 < |#eta| < 4", "l");
  leg->AddEntry(ratio37_eta2, "4 < |#eta| < 5", "l");
  leg->SetHeader("  CTEQ3D/CTEQ5L");
  leg->Draw();
  Sum -> cd(3);
  ratio57_eta1 -> Draw("e1p");
  ratio57_eta2 -> Draw("same");
//  ratio57 -> Draw("same");
  leg = new TLegend(0.12,0.83,0.5,0.97);
  leg->AddEntry(ratio57_eta1, "2 < |#eta| < 4", "l");
  leg->AddEntry(ratio57_eta2, "4 < |#eta| < 5", "l");
  leg->SetHeader("  GRV94M/CTEQ5L");
  leg->Draw();
  Sum -> SaveAs("~andrey/public/html/pythia/h034_ratio.png");


  Sum -> cd(1);
  x1_7 -> Draw("");
  x1_5 -> Draw("same");
  x1_3 -> Draw("same");
  leg = new TLegend(0.65,0.75,0.97,0.97);
  leg->AddEntry(x1_7, "x1 CTEQ5L", "l");
  leg->AddEntry(x1_5, "x1 GRV94 ", "l");
  leg->AddEntry(x1_3, "x1 CTEQ3D", "l");
  leg->SetHeader("x-variable");
  leg->Draw();
  Sum -> cd(2);
  xf_7 -> Draw("");
//  xf_5 -> Draw("same");
  xf_3 -> Draw("same");
  leg = new TLegend(0.30,0.75,0.60,0.97);
  leg->AddEntry(xf_7, "xf CTEQ5L", "l");
//  leg->AddEntry(xf_5, "xf GRV94 ", "l");
  leg->AddEntry(xf_3, "xf CTEQ3D", "l");
  leg->SetHeader("xf = |x1-x2|");
  leg->Draw();
  Sum -> cd(3);
  Sum_3 -> SetLogy();
  Q2_7->Draw("");
  Q2_5->Draw("same");
  Q2_3->Draw("same");
  leg = new TLegend(0.65,0.75,0.97,0.97);
  leg->AddEntry(Q2_7, "Q^{2} CTEQ5L", "l");
  leg->AddEntry(Q2_5, "Q^{2} GRV94 ", "l");
  leg->AddEntry(Q2_3, "Q^{2} CTEQ3D", "l");
  leg->SetHeader("Q^{2}");
  leg->Draw();
  Sum -> SaveAs("~andrey/public/html/pythia/h035_xplots.png");

//  Sum_1 -> SetLogy();
//  Sum_2 -> SetLogy();
  Sum -> cd(1);
  xf_3 -> Draw();
  xf_7 -> Draw("same");
  Sum -> cd(2);
//  xf_5 -> Draw();
  Sum -> cd(3);
  Sum -> SaveAs("~andrey/public/html/pythia/h036.png");

}