(* ::Package:: *)

Needs["ErrorBarPlots`"]

(* ::Section:: *)

(*Plot successrate over time*)

(* 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: edges

3: HH tri seq

4: {ccm1 #failed-attempts, ccm1 seq}

5: {ccm2 #failed-attempts, ccm2 seq}

*)

nlabels=Map["n = "<>ToString[#]&,gdata[[1,All,1,1,1]]];

taulabels=Map["\[Tau] = "<>ToString[#]&,gdata[[All,1,1,1,2]]];

(* ::Subsection:: *)

(*Canonical dataset (old code for multiple ccm runs per run)*)

(* For export *)

selectedRun=gdata[[1,1,1]];

measureSkip=1;