(* ::Package:: *)

gsraw=Import[NotebookDirectory[]<>"cpp/graphdata_spectrum.m"];
(* gsraw=SortBy[gsraw,{#[[1,1]]&,#[[1,2]]&}]; (* Sort by n and then by tau. The {} forces a *stable* sort because otherwise Mathematica sorts also on triangle count and other things. *) *)


gdata=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
(* Data format: *)
(* gdata[[ tau index, n index, run index , datatype index ]] *)
(* datatype index:
1: {n,tau}
2: etmt
*)
tauvalues=gdata[[All,1,1,1,2]];
nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]];
taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]];


ListPlot[gdata[[1,1,1,{4,6}]]]
ListPlot[gdata[[1,1,1,{5,7}]]]


Histogram[gdata[[1,1,1,{4,6}]],50]
Histogram[gdata[[1,1,1,{5,7}]],50]



