int graphs ()
{
//PG prepare the graphs
//PG ---- ---- ---- ---- ---- ---- ---- ----
TGraph ** tg_sig_par = new TGraph * [7] ; // [parameter][mass]
for (int k = 0 ; k < 7 ; ++k) tg_sig_par[k] = new TGraph (5) ;
TGraph ** tg_sAi_par = new TGraph * [9] ;
for (int k = 0 ; k < 9 ; ++k) tg_sAi_par[k] = new TGraph (5) ;
//PG fill the graphs from the input files
//PG ---- ---- ---- ---- ---- ---- ---- ----
int i = 0 ;
intCont ic_350 ("results_interference.350.1.root", 350) ;
ic_350.setsigGraphs (tg_sig_par, i) ;
ic_350.setsAiGraphs (tg_sAi_par, i++) ;
intCont ic_500 ("results_interference.500.1.root", 500) ;
ic_500.setsigGraphs (tg_sig_par, i) ;
ic_500.setsAiGraphs (tg_sAi_par, i++) ;
intCont ic_650 ("results_interference.650.1.root", 650) ;
ic_650.setsigGraphs (tg_sig_par, i) ;
ic_650.setsAiGraphs (tg_sAi_par, i++) ;
intCont ic_800 ("results_interference.800.1.root", 800) ;
ic_800.setsigGraphs (tg_sig_par, i) ;
ic_800.setsAiGraphs (tg_sAi_par, i++) ;
intCont ic_1000 ("results_interference.1000.1.root", 1000) ;
ic_1000.setsigGraphs (tg_sig_par, i) ;
ic_1000.setsAiGraphs (tg_sAi_par, i++) ;
//PG plot the interpolation for the signal
//PG ---- ---- ---- ---- ---- ---- ---- ----
compareFits (tg_sig_par, 7, tg_sAi_par, 9) ;
drawGraphs (tg_sig_par, "sig", 7) ;
drawGraphs (tg_sAi_par, "sAi", 9) ;
showInterpolation (tg_sig_par, "sig", 7) ;
showInterpolation (tg_sAi_par, "sAi", 9) ;
testInterpolation (tg_sig_par, "sig", 7) ;
testInterpolation (tg_sAi_par, "sAi", 9) ;
return 0 ;
}
int main(int argc, char *argv[])
{
testSequence();
testUnivariateStats();
testBivariateStats();
testUtility();
testVectorAlgebra();
testMatrixDeterminant();
testMatrixInverse();
testLUFactorisation();
testLUSolver();
testGaussianElimination();
testInterpolation();
testIntegration();
testRoot();
testTimeStep();
testRSignificanceTest();
testChiSquared();
testMultipleRegression();
test_golden_fit();
testSimplex();
testSimplex2();
testSimplex3();
test_simplex_fit();
return (EXIT_SUCCESS);
}
