# Script file: tweezer3D.lsf

run;

F=getresult("optical_force_volumetric","F_total");
Fx=F.F_totalx;
Fy=F.F_totaly;
Fz=F.F_totalz;
f=F.f;


plot(f/1e12,Fx,Fy,Fz,"f (THz)","Force (N)","Volumetric technique");
legend("Fx","Fy","Fz");
plot(f,Fx^(1/4),Fy^(1/4),Fz^(1/4),"f (THz)","f^4 dependence","Volumetric technique");
legend("Fx","Fy","Fz");



# compare results with MST monitor.  Results are less accurate because the forces 
# generated in this simulation are at the numerical limit of this technique.
F=getresult("optical_force_MST","F_total");
Fx=F.F_totalx;
Fy=F.F_totaly;
Fz=F.F_totalz;
f=F.f;

plot(f,Fx^(1/4),Fy^(1/4),Fz^(1/4),"f (THz)","f^4 dependence","MST technique");
legend("Fx","Fy","Fz");