{
    TGraph kinematics("a(p,p)a_kinematics_7MeV_p.txt","%*lf %lf %*lf %lf");
    kinematics.Scale(4.0);
    TGraph kin2(kinematics.GetN(), kinematics.GetY(),kinematics.GetX());
    //TFile fin("4He(pp)_candidate_kinematic_curve.root");
    //TH2F *h2 = (TH2F*)(fin.Get("Ef_thetaf_AnodeQQQR_TC1_PC1_pidlow0_1"));

//    TFile fin("../results_run15.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_qqqEf_sx3Ef_matrix_a1c2"));
    h2->Draw();
    kin2.Draw("L SAME");
    gPad->Modified();
    gPad->Update();
    while(gPad->WaitPrimitive());
    gPad->Clear();

    TGraph angles("a(p,p)a_kinematics_7MeV_p.txt","%lf %*lf %lf %*lf");
    TGraph angles2(angles.GetN(),angles.GetY(),angles.GetX());
    h2 = (TH2F*)(fin.Get("a(p,p)/ap_theta_vs_theta_qqq_sx3_a1c2"));
    h2->Draw();
    angles2.Draw("L SAME");
    gPad->Modified();
    gPad->Update();
    while(gPad->WaitPrimitive());
    gPad->SaveAs("c1.png");
}