src/figures/cantilever-sim/v-dep.asy

   1 import wtk_graph;
   2
   3 size(15cm,10cm,IgnoreAspect);
   4
   5 scale(Log, Linear);
   6 real vscale = 1e9;
   7 real fscale = 1e12;
   8
   9 real vmin = 0.2e-6;
  10 real vmax = 5e-6;
  11
  12 /* f(x) = A + log10(x) + B */
  13 real fn_logxliny(real x, real[] params) {
  14   return params[0] * log10(x) + params[1];
  15 }
  16
  17 graphFile("v-dep.d/v_dep_127_8", vscale, fscale, p=phard, mpath=m8,
  18           t=math(units("127","pN/nm"))+", 8 domains", dots=true);
  19 graphFile("v-dep.d/v_dep_27_8", vscale, fscale, p=pmed, mpath=m8,
  20           t=math(units("27","pN/nm"))+", 8 domains", dots=true);
  21 graphFile("v-dep.d/v_dep_127_30", vscale, fscale, p=phard, mpath=m30,
  22           t=math(units("127","pN/nm"))+", 30 domains", dots=true);
  23 graphFile("v-dep.d/v_dep_27_30", vscale, fscale, p=pmed, mpath=m30,
  24           t=math(units("27","pN/nm"))+", 30 domains", dots=true);
  25 graphFile("v-dep.d/v_dep_0.1_1", vscale, fscale, p=psoft, mpath=m1,
  26           t=math(units("0.1","pN/nm"))+", 1 domain", dots=true);
  27 graphFile("v-dep.d/v_dep_0.1_30", vscale, fscale, p=psoft, mpath=m30,
  28           t=math(units("0.1","pN/nm"))+", 30 domains", dots=true);
  29 fitFile("v-dep.d/v_dep_127_8.fit.dat", f=fn_logxliny,
  30         xmin=1e-9, xmax=1e-5, xscale=vscale, yscale=fscale, p=phard);
  31 fitFile("v-dep.d/v_dep_27_8.fit.dat", f=fn_logxliny,
  32         xmin=1e-9, xmax=1e-5, xscale=vscale, yscale=fscale, p=pmed);
  33 fitFile("v-dep.d/v_dep_127_30.fit.dat", f=fn_logxliny,
  34         xmin=1e-9, xmax=1e-5, xscale=vscale, yscale=fscale, p=phard);
  35 fitFile("v-dep.d/v_dep_27_30.fit.dat", f=fn_logxliny,
  36         xmin=1e-9, xmax=1e-5, xscale=vscale, yscale=fscale, p=pmed);
  37 fitFile("v-dep.d/v_dep_0.1_1.fit.dat", f=fn_logxliny,
  38         xmin=1e-9, xmax=1e-5, xscale=vscale, yscale=fscale, p=psoft);
  39 fitFile("v-dep.d/v_dep_0.1_30.fit.dat", f=fn_logxliny,
  40         xmin=1e-9, xmax=1e-5, xscale=vscale, yscale=fscale, p=psoft);
  41
  42 ylimits(150, 340, crop=true);
  43
  44 label(sLabel("Pulling speed dependence"), point(N), N);
  45 xaxis(sLabel("Pulling speed (nm/s)"), BottomTop, LeftTicks);
  46 yaxis(sLabel("Unfolding force (pN)"), LeftRight, RightTicks);
  47
  48 add(legend(), point(E), 20E, UnFill);