diff --git a/triangle_gcm_initial_analysis.m b/triangle_gcm_initial_analysis.m
index 0fc898477bbfe8bdda80c1f049b143089224cba6..93f7e75a5fe39524925a904fed6725e3988cef08 100644
--- a/triangle_gcm_initial_analysis.m
+++ b/triangle_gcm_initial_analysis.m
@@ -142,13 +142,18 @@ Export[NotebookDirectory[]<>"plots/ccm_construction_successrate9.pdf",plot3]
 (*CCMu rates only*)
 
 
-successrates=Map[#[[3,2]]/100&,gdata,{3}];
-successrates[[1,1]]=gdata2[[1,1,All,2]]; (* New datafile test *)
-successrates=Map[Clip[#,{0,0.99999}]&,successrates,{3}];
+successrates=Map[#[[3,2]]/100&,gdata,{3}]; (* Old datafile *)
 legends=Map["\[Tau] = "<>ToString[#[[1,1,2]]]<>" ; avg = "<>ToString[NumberForm[N[Mean[#[[All,3,2]]/100]],3]]&,gdata,{2}];
 
-datasets={successrates[[1,1]], successrates[[5,1]], successrates[[9,1]]};
-selectedLegends={legends[[1,1]],legends[[5,1]],legends[[9,1]]};
+
+successrates=gdata2[[All,All,All,2]]; (* New datafile *)
+legends=Map["\[Tau] = "<>ToString[#[[1,1,2]]]<>" ; avg = "<>ToString[NumberForm[N[Mean[#[[All,2]]]],{4,4}]]&,gdata2,{2}];
+
+
+successrates=Map[Clip[#,{0,0.99999}]&,successrates,{3}];
+
+datasets={successrates[[1,1]], successrates[[2,1]], successrates[[3,1]]};
+selectedLegends={legends[[1,1]],legends[[2,1]],legends[[3,1]]};
 
 plot1=Histogram[datasets,{-0.025,1.025,0.05},"Probability",
 PlotRange->{{0,1},{0,1}},