a=Select[Rest[Import["~/doc/research/pub_wolff-field/data/correlation-times/ising-3d.dat","Data"]],#≠{}&];

δ3=4.78984;β3=0.326419;ν3=0.629971;z3=0.28;

p1=Show[ListLogLogPlot[{(#3 / #2)((#1)^(δ3 β3/ν3)),Around[#4,#5]Around[#6,#7]((#1)^(-z3 -3))}&@@@#&/@GatherBy[a,First],Frame->True,FrameLabel->{"h (L^(-β δ / ν))","τ (L^(-z))"},Joined->True,LabelStyle->{FontSize->12}],LogLogPlot[0.55(x^(-z3 ν3/(β3 δ3))),{x,(10^(-6)),(10^6)},PlotStyle->Red],Epilog->{Inset[Style["3D Ising in uniform field,",12],{-7,-1.3}],Inset["T = (Subscript[T, c])",{-12.2,-1.45}]}]

β2Vector3D=0.3470;α2Vector3D=−0.011;δ2Vector3D=(2-α2Vector3D)/β2Vector3D-1;ν2Vector3D=0.6703;γ2Vector3D=1.3169;z2=0.17;

b=Select[Rest[Import["~/doc/research/pub_wolff-field/data/correlation-times/planar-pert.dat","Data"]],#≠{}&&#[[6]]>0&];

p2=Show[ListLogLogPlot[{(#3 / #2)((#1)^(δ2Vector3D β2Vector3D/ν2Vector3D)),Around[#4,#5]Around[#6,#7]((#1)^(-z2 -3))}&@@@#&/@GatherBy[b,First],Frame->True,FrameLabel->{"h (L^(-β δ / ν))","τ (L^(-z))"},Joined->True,LabelStyle->{FontSize->12}],LogLogPlot[1.5(x^(-z2 ν2Vector3D/(β2Vector3D δ2Vector3D))),{x,8 (10^(-7)),2 (10^6)},PlotStyle->Red],Epilog->{Inset[Style["3D XY in uniform vector field,",12],{-3.5,-0.07}],Inset["T = (Subscript[T, c])",{-10.9,-0.15}]}]

Export["~/doc/research/pres/statphys27/figs/times_3dising_new_alert.pdf",p1];

Export["~/doc/research/pres/statphys27/figs/times_3dxy_new_alert.pdf",p2];