Alpha and beta fixed bugs

code/RDownScriptalphabeta.m

@@ -30,10 +30,10 @@ SetDirectory[rootpath<>"code"];
 SetOptions[InputNotebook[], AutoGeneratedPackage -> Automatic]
 
 
-mysxscase="/Users/RayneLiu/Desktop/Research/AEI/Simulations/SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.440";
-tshift=0;
-npoints=100;
-ntones=7;
+mysxscase=rootpath<>"/SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.440";
+tshift=19;
+npoints=200;
+ntones=1;
 allspectrum = {x0,x1,x2,x3,x4,x5,x6,x7};
 spectrum = allspectrum[[;;ntones+1]];
 \[Omega]fact=0.1;
@@ -41,7 +41,7 @@ spectrum = allspectrum[[;;ntones+1]];
 params = {tshift, ntones, \[Omega]fact, \[Tau]fact}
 
 
-combinations = DeleteCases[Tuples[Range[0, ntones],2],Alternatives@@Subsets[Reverse[Range[0, ntones]],{2}]]
+combinations = DeleteCases[Tuples[Range[1, ntones],2],Alternatives@@Subsets[Reverse[Range[0, ntones]],{2}]]
 
 
 (* Change them accordingly *)
@@ -138,7 +138,9 @@ ansatz=OvertoneModelV2[ntones,{\[Eta],\[Chi]1,\[Chi]2},t0,"Fit\[Alpha]"->combina
 cfit=NonlinearModelFit[data,ansatz,spectrum,t];
 cfitd=Transpose[{data[[All,1]],Normal[cfit]/.t->data[[All,1]]}];
 comp\[Alpha]s = spectrum/.cfit["BestFitParameters"];
-{{randomvar1, randomvar2}, combinations[[k]], (*ansatz, cfit["BestFitParameters"], comp\[Alpha]s,*) Re[Sqrt[comp\[Alpha]s * Conjugate[comp\[Alpha]s]]], Mod[Arg[comp\[Alpha]s], 2 Pi], {tshift, 1-EasyMatchT[data,cfitd,t0,t0+90]}},{j,npoints}], {k, Length[combinations]}];
+mismatchlist = {{randomvar1, randomvar2}, combinations[[k]], {tshift, 1-EasyMatchT[data,cfitd,t0,t0+90]}};
+ampandphs = Transpose[{Re[Sqrt[comp\[Alpha]s * Conjugate[comp\[Alpha]s]]], Mod[Arg[comp\[Alpha]s], 2 Pi]}];
+Join[mismatchlist, ampandphs], {j,npoints}], {k, Length[combinations]}];
 
 
 tab\[Alpha]\[Beta]s=Flatten[Table[Table[randomvar1= RandomReal[{-\[Omega]fact, \[Omega]fact}];randomvar2 = RandomReal[{-\[Tau]fact, \[Tau]fact}];
@@ -146,7 +148,9 @@ ansatz=OvertoneModelV2[ntones,{\[Eta],\[Chi]1,\[Chi]2},t0,"Fit\[Alpha]"->{k},"Fi
 cfit=NonlinearModelFit[data,ansatz,spectrum,t];
 cfitd=Transpose[{data[[All,1]],Normal[cfit]/.t->data[[All,1]]}];
 comp\[Alpha]\[Beta]s = spectrum/.cfit["BestFitParameters"];
-{{randomvar1, randomvar2}, {k, l}, (*ansatz, cfit["BestFitParameters"], comp\[Alpha]\[Beta]s,*) Re[Sqrt[comp\[Alpha]\[Beta]s * Conjugate[comp\[Alpha]\[Beta]s]]], Mod[Arg[comp\[Alpha]\[Beta]s], 2 Pi], {tshift, 1-EasyMatchT[data,cfitd,t0,t0+90]}}, {j,npoints}], {k, Range[0, ntones]}, {l, Range[0, ntones]}], 1];
+mismatchlist = {{randomvar1, randomvar2}, {k, l}, {tshift, 1-EasyMatchT[data,cfitd,t0,t0+90]}};
+ampandphs = Transpose[{Re[Sqrt[comp\[Alpha]\[Beta]s * Conjugate[comp\[Alpha]\[Beta]s]]], Mod[Arg[comp\[Alpha]\[Beta]s], 2 Pi]}];
+Join[mismatchlist, ampandphs], {j,npoints}], {k, Range[1, ntones]}, {l, Range[1, ntones]}], 1];
 
 
 tab\[Beta]s=Table[Table[randomvar1= RandomReal[{-\[Tau]fact, \[Tau]fact}];randomvar2 = RandomReal[{-\[Tau]fact, \[Tau]fact}];
@@ -154,13 +158,15 @@ ansatz=OvertoneModelV2[ntones,{\[Eta],\[Chi]1,\[Chi]2},t0,"Fit\[Alpha]"->{},"Fit
 cfit=NonlinearModelFit[data,ansatz,spectrum,t];
 cfitd=Transpose[{data[[All,1]],Normal[cfit]/.t->data[[All,1]]}];
 comp\[Beta]s = spectrum/.cfit["BestFitParameters"];
-{{randomvar1, randomvar2}, combinations[[k]], (*ansatz, cfit["BestFitParameters"], comp\[Beta]s,*) Re[Sqrt[comp\[Beta]s * Conjugate[comp\[Beta]s]]], Mod[Arg[comp\[Beta]s], 2 Pi], {tshift, 1-EasyMatchT[data,cfitd,t0,t0+90]}},{j,npoints}], {k, Length[combinations]}];
+mismatchlist = {{randomvar1, randomvar2}, combinations[[k]], {tshift, 1-EasyMatchT[data,cfitd,t0,t0+90]}};
+ampandphs = Transpose[{Re[Sqrt[comp\[Beta]s * Conjugate[comp\[Beta]s]]], Mod[Arg[comp\[Beta]s], 2 Pi]}];
+Join[mismatchlist, ampandphs], {j,npoints}], {k, Length[combinations]}];
 
 
-Export[rootpath<>"plots/params.fits", params]
-Export[rootpath<>"plots/freq_bfitpars_mmatch_as.fits", tab\[Alpha]s]
-Export[rootpath<>"plots/freq_bfitpars_mmatch_abs.fits", tab\[Alpha]\[Beta]s]
-Export[rootpath<>"plots/freq_bfitpars_mmatch_bs.fits", tab\[Beta]s]
+Export[rootpath<>"plots/params_n="<>ToString[ntones]<>".fits", params]
+Export[rootpath<>"plots/freq_bfitpars_mmatch_as_n="<>ToString[ntones]<>".fits", tab\[Alpha]s]
+Export[rootpath<>"plots/freq_bfitpars_mmatch_abs_n="<>ToString[ntones]<>".fits", tab\[Alpha]\[Beta]s]
+Export[rootpath<>"plots/freq_bfitpars_mmatch_bs_n="<>ToString[ntones]<>".fits", tab\[Beta]s]