{
    // ================================================================

    gROOT->Reset();
    gStyle->SetOptStat(00000);

    TCanvas *c1 = new TCanvas("c1", "c1",200,0,1200,1200);
    c1->Range(0,0,1,1);
    c1->SetFillColor(0);
    c1->SetBorderMode(0);
    c1->SetBorderSize(2);
    c1->SetFrameBorderMode(0);


    // get the RGA data and the SIDIS sim results

    TFile *inputfile1 = new TFile("data.root");
    TFile *inputfile2 = new TFile("sim10p6sidiscut45na.root");

    // SIDIS sim backgnd only, peak only.
    
    TH1D *my_mmBkgd10   = (TH1D*)inputfile2->Get("expMM_bg_Pmmbin20");
    TH1D *my_mmPeak10   = (TH1D*)inputfile2->Get("expMM_sig_Pmmbin20");

    gStyle->SetStatW(0.4);
    gStyle->SetStatFontSize(0.05);

    TSpline3 *splineBkgd10 = new TSpline3(my_mmBkgd10);
    TSpline3 *splinePeak10 = new TSpline3(my_mmPeak10);
    gROOT->ProcessLine(".L mySplineFit.C");
    hist10 = (TH1F*)inputfile1->Get("expMM_Pmmbin20");
    // hist1->GetXaxis()->SetRangeUser(0, 1.5);
    myfuncTest10();
    myfit10(hist10);


}