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

        TFile gg("temp23_gg.root", "OPEN");
        TFile fg("temp22_fg.root", "OPEN");
        TFile ff("temp21_ff.root", "OPEN");

        TH1F *x1_3 = (TH1F*)gg.Get("x1_5");
        TH1F *x1_5 = (TH1F*)fg.Get("x1_5");
        TH1F *x1_7 = (TH1F*)ff.Get("x1_5");

        TH1F *x2_3 = (TH1F*)gg.Get("x2_5");
        TH1F *x2_5 = (TH1F*)fg.Get("x2_5");
        TH1F *x2_7 = (TH1F*)ff.Get("x2_5");

        TH1F *xf_3 = (TH1F*)gg.Get("xf_5");
        TH1F *xf_5 = (TH1F*)fg.Get("xf_5");
        TH1F *xf_7 = (TH1F*)ff.Get("xf_5");

        TH1F *XX_3 = (TH1F*)gg.Get("XX_5");
        TH1F *XX_5 = (TH1F*)fg.Get("XX_5");
        TH1F *XX_7 = (TH1F*)ff.Get("XX_5");

        TH1F *xx_3 = (TH1F*)gg.Get("xx_5");
        TH1F *xx_5 = (TH1F*)fg.Get("xx_5");
        TH1F *xx_7 = (TH1F*)ff.Get("xx_5");

        TH1F *Q2_3 = (TH1F*)gg.Get("Q2_5");
        TH1F *Q2_5 = (TH1F*)fg.Get("Q2_5");
        TH1F *Q2_7 = (TH1F*)ff.Get("Q2_5");

  x1_3 -> SetLineColor(kBlue);
  x2_3 -> SetLineColor(kBlue+1);
  xf_3 -> SetLineColor(kBlue);
  XX_3 -> SetLineColor(kOrange);
  xx_3 -> SetLineColor(kSpring);
  XX_3 -> SetFillColor(kOrange);
  xx_3 -> SetFillColor(kSpring);
  XX_3 -> SetFillStyle(3010);
  xx_3 -> SetFillStyle(3012);

  Q2_3 -> SetLineColor(kBlue+3);
  x1_5 -> SetLineColor(kGreen);
  x2_5 -> SetLineColor(kGreen+1);
  xf_5 -> SetLineColor(kBlue);
  XX_5 -> SetLineColor(kOrange);
  xx_5 -> SetLineColor(kSpring);
  XX_5 -> SetFillColor(kOrange);
  xx_5 -> SetFillColor(kSpring);
  XX_5 -> SetFillStyle(3010);
  xx_5 -> SetFillStyle(3012);

  Q2_5 -> SetLineColor(kGreen+3);
  x1_7 -> SetLineColor(kRed);
  x2_7 -> SetLineColor(kRed+1);
  xf_7 -> SetLineColor(kBlue);
  XX_7 -> SetLineColor(kOrange);
  xx_7 -> SetLineColor(kSpring);
  XX_7 -> SetFillColor(kOrange);
  xx_7 -> SetFillColor(kSpring);
  XX_7 -> SetFillStyle(3010);
  xx_7 -> SetFillStyle(3012);
  xf_7 -> SetFillStyle(3020);

  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)");

  XX_3 -> SetXTitle("log10 (x)");
  XX_5 -> SetXTitle("log10 (x)");
  XX_7 -> SetXTitle("log10 (x)");
  xx_3 -> SetXTitle("log10 (x)");
  xx_5 -> SetXTitle("log10 (x)");
  xx_7 -> SetXTitle("log10 (x)");
  Q2_3 -> SetXTitle("GeV^{2}");
  Q2_5 -> SetXTitle("GeV^{2}");
  Q2_7 -> SetXTitle("GeV^{2}");

        TH2F *x1_Q2_3 = (TH2F*)gg.Get("x1_Q2_5");
        TH2F *x1_Q2_5 = (TH2F*)fg.Get("x1_Q2_5");
        TH2F *x1_Q2_7 = (TH2F*)ff.Get("x1_Q2_5");
        TH2F *x2_Q2_3 = (TH2F*)gg.Get("x2_Q2_5");
        TH2F *x2_Q2_5 = (TH2F*)fg.Get("x2_Q2_5");
        TH2F *x2_Q2_7 = (TH2F*)ff.Get("x2_Q2_5");
//  x1_Q2_3 -> SetFillColor(kViolet);
//  x1_Q2_5 -> SetFillColor(kViolet);
//  x1_Q2_7 -> SetFillColor(kViolet);
  x1_Q2_3 -> SetXTitle("log10 (x)");
  x1_Q2_5 -> SetXTitle("log10 (x)");
  x1_Q2_7 -> SetXTitle("log10 (x)");
  x1_Q2_3 -> SetYTitle("GeV^{2}");
  x1_Q2_5 -> SetYTitle("GeV^{2}");
  x1_Q2_7 -> SetYTitle("GeV^{2}");

  x1_Q2_3 -> SetTitle("Q^{2} vs logx");
  x1_Q2_5 -> SetTitle("Q^{2} vs logx");
  x1_Q2_7 -> SetTitle("Q^{2} vs logx");

        TH2F *XX_pt_3 = (TH2F*)gg.Get("XX_pt_5");
        TH2F *xx_pt_3 = (TH2F*)gg.Get("xx_pt_5");
        TH2F *XX_pt_5 = (TH2F*)fg.Get("XX_pt_5");
        TH2F *xx_pt_5 = (TH2F*)fg.Get("xx_pt_5");
        TH2F *XX_pt_7 = (TH2F*)ff.Get("XX_pt_5");
        TH2F *xx_pt_7 = (TH2F*)ff.Get("xx_pt_5");

  xx_pt_3 -> SetTitle("x1, x2  vs  Pt jet 1,2 at 3<|eta|<5 ");
  xx_pt_5 -> SetTitle("x1, x2  vs  Pt jet 1,2 at 3<|eta|<5 ");
  xx_pt_7 -> SetTitle("x1, x2  vs  Pt jet 1,2 at 3<|eta|<5 ");
  xx_pt_3 -> SetXTitle("log10(x1, x2)");
  xx_pt_5 -> SetXTitle("log10(x1, x2)");
  xx_pt_7 -> SetXTitle("log10(x1, x2)");
  xx_pt_3 -> SetYTitle("Pt of jet; GeV");
  xx_pt_5 -> SetYTitle("Pt of jet; GeV");
  xx_pt_7 -> SetYTitle("Pt of jet; GeV");


    TH2F *XX_xx_1_1_3 = (TH2F*)gg.Get("XX_xx_1_1_5");
    TH2F *XX_xx_2_1_3 = (TH2F*)gg.Get("XX_xx_2_1_5");
    TH2F *XX_xx_3_1_3 = (TH2F*)gg.Get("XX_xx_3_1_5");
    TH2F *XX_xx_1_1_5 = (TH2F*)fg.Get("XX_xx_1_1_5");
    TH2F *XX_xx_2_1_5 = (TH2F*)fg.Get("XX_xx_2_1_5");
    TH2F *XX_xx_3_1_5 = (TH2F*)fg.Get("XX_xx_3_1_5");
    TH2F *XX_xx_1_1_7 = (TH2F*)ff.Get("XX_xx_1_1_5");
    TH2F *XX_xx_2_1_7 = (TH2F*)ff.Get("XX_xx_2_1_5");
    TH2F *XX_xx_3_1_7 = (TH2F*)ff.Get("XX_xx_3_1_5");

    TH2F *XX_xx_1_2_3 = (TH2F*)gg.Get("XX_xx_1_2_5");
    TH2F *XX_xx_2_2_3 = (TH2F*)gg.Get("XX_xx_2_2_5");
    TH2F *XX_xx_3_2_3 = (TH2F*)gg.Get("XX_xx_3_2_5");
    TH2F *XX_xx_1_2_5 = (TH2F*)fg.Get("XX_xx_1_2_5");
    TH2F *XX_xx_2_2_5 = (TH2F*)fg.Get("XX_xx_2_2_5");
    TH2F *XX_xx_3_2_5 = (TH2F*)fg.Get("XX_xx_3_2_5");
    TH2F *XX_xx_1_2_7 = (TH2F*)ff.Get("XX_xx_1_2_5");
    TH2F *XX_xx_2_2_7 = (TH2F*)ff.Get("XX_xx_2_2_5");
    TH2F *XX_xx_3_2_7 = (TH2F*)ff.Get("XX_xx_3_2_5");

    TH2F *XX_xx_1_3_3 = (TH2F*)gg.Get("XX_xx_1_3_5");
    TH2F *XX_xx_2_3_3 = (TH2F*)gg.Get("XX_xx_2_3_5");
    TH2F *XX_xx_3_3_3 = (TH2F*)gg.Get("XX_xx_3_3_5");

    TH2F *XX_xx_1_3_5 = (TH2F*)fg.Get("XX_xx_1_3_5");
    TH2F *XX_xx_2_3_5 = (TH2F*)fg.Get("XX_xx_2_3_5");
    TH2F *XX_xx_3_3_5 = (TH2F*)fg.Get("XX_xx_3_3_5");


    TH2F *jets_Q2_3 = (TH2F*)gg.Get("jets_Q2_5");
    TH2F *jets_Q2_5 = (TH2F*)fg.Get("jets_Q2_5");
    TH2F *jets_Q2_7 = (TH2F*)ff.Get("jets_Q2_5");


  TCanvas *Sum = new TCanvas("Sum", "CTEQ, gg->jets; fg->jets; ff->jets", 10, 10, 1200, 400);
  Sum -> Divide(3,1);
  Sum -> SetHighLightColor(45);
  Sum_1 -> SetGrid();
  Sum_2 -> SetGrid();
  Sum_3 -> SetGrid();


  Sum -> cd(1);
  XX_3 -> Draw("");
  xx_3 -> Draw("same");
  xf_3 -> Draw("same");
  leg = new TLegend(0.05,0.75,0.35,0.97);
  leg->AddEntry(XX_3, "Xmax", "l");
  leg->AddEntry(xx_3, "Xmin", "l");
  leg->AddEntry(xf_3, "xf=|x1-x2|", "l");
  leg->SetHeader("gg CTEQ");
  leg->Draw();

  Sum -> cd(2);
  XX_5 -> Draw("");
  xx_5 -> Draw("same");
  xf_5 -> Draw("same");
  leg = new TLegend(0.05,0.75,0.35,0.97);
  leg->AddEntry(XX_5, "Xmax", "l");
  leg->AddEntry(xx_5, "Xmin", "l");
  leg->AddEntry(xf_5, "xf=|x1-x2|", "l");
  leg->SetHeader("fg CTEQ");
  leg->Draw();

  Sum -> cd(3);
  XX_7 -> Draw("");
  xx_7 -> Draw("same");
  xf_7 -> Draw("same");
  leg = new TLegend(0.05,0.75,0.35,0.97);
  leg->AddEntry(XX_7, "Xmax", "l");
  leg->AddEntry(xx_7, "Xmin", "l");
  leg->AddEntry(xf_7, "xf=|x1-x2|", "l");
  leg->SetHeader("ff CTEQ");
  leg->Draw();
  Sum -> SaveAs("~andrey/public/html/pythia/h052_xplots.png");

 gStyle->SetOptTitle(100); 
//  Sum_1 -> SetLogy();
//  Sum_2 -> SetLogy();
Sum_3 -> SetLogy(0);
  Sum -> cd(1);
  x1_Q2_3 -> Draw("col");
  Sum -> cd(2);
  x1_Q2_5 -> Draw("colz");
  Sum -> cd(3);
  x1_Q2_7->Draw("col");
  Sum -> SaveAs("~andrey/public/html/pythia/h053_xQ2.png");

  Sum_1 -> Clear();
  Sum_2 -> Clear();
  Sum_3 -> Clear();

Sum -> cd(1);
xx_pt_3 -> Draw("col");
//XX_pt_3->SetPalette(3);
XX_pt_3 -> Draw("same col");
Sum -> cd(2);
xx_pt_5 -> Draw("col");
XX_pt_5 -> Draw("same colz");
Sum -> cd(3);
xx_pt_7 -> Draw("col");
XX_pt_7 -> Draw("same col");
Sum -> SaveAs("~andrey/public/html/pythia/h054_xpt.png");


  Sum -> cd(1);
  XX_xx_1_1_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_1_1_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_1_1_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_1_1.png");

  Sum -> cd(1);
  XX_xx_2_1_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_2_1_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_2_1_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_2_1.png");

  Sum -> cd(1);
  XX_xx_3_1_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_3_1_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_3_1_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_3_1.png");

  Sum -> cd(1);
  XX_xx_1_2_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_1_2_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_1_2_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_1_2.png");

  Sum -> cd(1);
  XX_xx_2_2_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_2_2_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_2_2_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_2_2.png");

  Sum -> cd(1);
  XX_xx_3_2_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_3_2_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_3_2_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_3_2.png");


  Sum -> cd(1);
  XX_xx_1_3_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_1_3_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_1_3_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_1_3.png");

  Sum -> cd(1);
  XX_xx_2_3_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_2_3_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_2_3_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_2_3.png");

  Sum -> cd(1);
  XX_xx_3_3_3 -> Draw("col");
  Sum -> cd(2);
  XX_xx_3_3_5 -> Draw("colz");
  Sum -> cd(3);
  XX_xx_3_3_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h055_Xx_3_3.png");

  Sum -> cd(1);
  jets_Q2_3 -> Draw("col");
  Sum -> cd(2);
  jets_Q2_5 -> Draw("colz");
  Sum -> cd(3);
  jets_Q2_7 -> Draw("col");
    Sum -> SaveAs("~andrey/public/html/pythia/h056_qjets.png");
}