{
    TGraph kinematics("a(p,p)a_kinematics_7MeV_p.txt","%*lf %*lf %lf %lf");
    kinematics.Scale(4.0);
    kinematics.SetLineColor(kBlue);
    kinematics.SetLineWidth(2.);
    //TFile fin("4He(pp)_candidate_kinematic_curve.root");
    //TH2F *h2 = (TH2F*)(fin.Get("Ef_thetaf_AnodeQQQR_TC1_PC1_pidlow0_1"));

//    TFile fin("../results_run18.root");
//    TFile fin("out_zminus5_catima.root");
    TFile fin("out.root");
    //TH2F *h2 = (TH2F*)(fin.Get("Ef_thetaf_AnodeQQQR_TC1_PC1_pidlow0_1"));
//    TFile fin("run15_sx3_qqq_alphas_kinematic.root");
    TH2F *h2 = (TH2F*)(fin.Get("a(p,p)/ap_Ef_vs_theta_qqq_a1c2"));
    h2->Draw("box");
    h2->RebinY(8);
    kinematics.Draw("L SAME");
    gPad->Modified();
    gPad->Update();
    while(gPad->WaitPrimitive());

    TGraph kinematicsp("a(p,p)a_kinematics_7MeV_p.txt","%lf %lf %*lf %*lf");
    kinematicsp.SetLineWidth(2.);
    TH2F *h2p = (TH2F*)(fin.Get("a(p,p)/ap_Ef_vs_theta_sx3_a1c2"));
    h2p->RebinY(8);
    h2p->SetLineColor(kBlack);
    //h2p->SetLineColor(kRed);
    h2p->Draw("box same");
    kinematicsp.Draw("L SAME");

    gPad->SaveAs("c1.png");

}