diff options
author | Eugeniy Mikhailov <evgmik@gmail.com> | 2014-05-13 15:35:26 -0400 |
---|---|---|
committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2014-05-13 15:58:56 -0400 |
commit | d1b62466d81b164656588ac1c55acff15ee9ea43 (patch) | |
tree | cce78f90b0768361c4a268b946848842e1e47456 /examples/A_beam_visualization_script.m | |
download | optics_toolkit-d1b62466d81b164656588ac1c55acff15ee9ea43.tar.gz optics_toolkit-d1b62466d81b164656588ac1c55acff15ee9ea43.zip |
initial
The optics_toolkit code taken from
http://mercury.pr.erau.edu/~greta9a1/downloads/index.html
the older version is also available at mathwork web site
http://www.mathworks.com/matlabcentral/fileexchange/15459-basic-paraxial-optics-toolkit
Diffstat (limited to 'examples/A_beam_visualization_script.m')
-rw-r--r-- | examples/A_beam_visualization_script.m | 601 |
1 files changed, 601 insertions, 0 deletions
diff --git a/examples/A_beam_visualization_script.m b/examples/A_beam_visualization_script.m new file mode 100644 index 0000000..e58da4c --- /dev/null +++ b/examples/A_beam_visualization_script.m @@ -0,0 +1,601 @@ +% This is a fairly elaborate example of how to use the beam visualization
+% tools and beam propagation tools in the "axial" subdirectory. Basically,
+% the functions in the "axial" subdirectory implement an ABCD matrix
+% propagation environment. The functions beamplot.m and lensplot.m are
+% simple examples of visualization tools used to fair effect by this
+% script. This script propagates the main laser beam at the LIGO Livingston
+% observatory from the modecleaner to the various output ports using the
+% ABCD approximation. Then a figure is drawn with the various beam paths
+% illustrated.
+
+%----------------------------------------------------------------------------------------------
+% VARIABLE CONTROL
+
+var_mcwaist=0; % Allow (1) or disallow (0) beam width variation at MC waist
+var_rmwidth=0; % Allow (1) or disallow (0) beam width variation insided PRM (at RM)
+var_deltaMMT2pos=1; % Allow (1) or disallow (0) MMT2 position variation
+
+postmc_power=3; % Input power [Watts]
+%----------------------------------------------------------------------------------------------
+fontsize=0.01;
+figure(1);
+
+for s=1:1
+switch s
+case 1
+ wmc=0.01; % MC waist width [inches] (only used if var_mcwaist=1)
+ wrm=4.8/2.54; % PRM beam width specified at RM [inches] (only used if var_rmwidth=1)
+ deltaMMT2pos=3.18; % MMT2, distance from design position [inches] pos. value lengthens telescope (only used if var_deltaMMT2pos=1)
+ linestyle='-';
+ linecolor='b';
+case 2
+ wmc=0.03;
+ wrm=4.8/2.54;
+ deltaMMT2pos=1.5;
+ linestyle='-';
+ linecolor='r';
+case 3
+ wmc=0.06413;
+ wrm=4.8/2.54;
+ deltaMMT2pos=3.18;
+ linestyle='-';
+ linecolor='g';
+end
+
+beam_propagate_llo_asbuilt;
+
+%-----------------------------------------------------------------------------
+% INPUT BEAM
+%-----------------------------------------------------------------------------
+x=0; y=86;
+resolution=5000;
+
+propdist=52-1.48;
+[h,x,y]=beamplot(qMCwaist,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_smalloptic,0,x,y,0,1064e-9/0.0254,'MMT1');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=533.55;
+[h,x,y]=beamplot(qMMT1FF,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_smalloptic,0,x,y,0,1064e-9/0.0254,'MMT2');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=565.99;
+[h,x,y]=beamplot(qMMT2FF,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'MMT3');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=637.41;
+[h,x,y]=beamplot(qMMT3FF,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'RM');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=173.13;
+[h,x,y]=beamplot(qRMFF,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'BS');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=191.40;
+[h,x,y]=beamplot(qBSFF,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'ITM');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%propdist=157284.7/10;
+%[h,x,y]=beamplot(qITMXFF,[0:propdist/resolution:propdist],x,y); hold on;
+%set(h,'Color',linecolor);
+%set(h,'LineStyle',linestyle);
+%[h,th]=lensplot(qd_largeoptic,0,x,y,0,1064e-9/0.0254,'ETM');
+%set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+%-----------------------------------------------------------------------------
+% REFL
+%-----------------------------------------------------------------------------
+
+% REFL In vac
+x=0; y=73;
+resolution=5000;
+
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'RM');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=637.41;
+[h,x,y]=beamplot(qRMrear,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'MMT3');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=565.99;
+[h,x,y]=beamplot(qMMT3FFreturning,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_smalloptic,0,x,y,0,1064e-9/0.0254,'MMT2',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=533.55;
+[h,x,y]=beamplot(qMMT2FFreturning,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_smalloptic,0,x,y,0,1064e-9/0.0254,'MMT1',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=17.21-1.48;
+[h,x,y]=beamplot(qMMT1FFreturning,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_faraday,0,x,y,0,1064e-9/0.0254,'Faraday',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=60.8;
+[h,x,y]=beamplot(qFaradayreturning,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'viewport',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% REFL Table, common path
+propdist=43;
+[h,x,y]=beamplot(qREFLviewport,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+xcommon=x; ycommon=y;
+
+% REFL Table, 62 MHz LSC path
+x=xcommon; y=ycommon;
+
+propdist=40.25;
+[h,x,y]=beamplot(qREFLL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2/25.4,0,x,y,0,1064e-9/0.0254,'60 MHz PD',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% REFL Table, 25 MHz LSC path
+x=xcommon; y=ycommon-4;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=30.4;
+[h,x,y]=beamplot(qREFLL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'~180 mm',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=6.5;
+[h,x,y]=beamplot(qREFLL2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2/25.4,0,x,y,0,1064e-9/0.0254,'25 MHz PD',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% REFL Table, camera path
+x=xcommon; y=ycommon-8;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=27;
+[h,x,y]=beamplot(qREFLL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(5/25.4,0,x,y,0,1064e-9/0.0254,'Camera',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% REFL Table, WFS3 path
+x=xcommon; y=ycommon-12;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=72.15;
+[h,x,y]=beamplot(qREFLL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'401 mm',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=20.75;
+[h,x,y]=beamplot(qREFLL3,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'-50 mm',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=4.75;
+[h,x,y]=beamplot(qREFLL4,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(10/25.4,0,x,y,0,1064e-9/0.0254,'WFS3',2);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+% REFL Table, WFS4 path
+x=xcommon; y=ycommon-16;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=65.55;
+[h,x,y]=beamplot(qREFLL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'~1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=37.75;
+[h,x,y]=beamplot(qREFLL5,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(10/25.4,0,x,y,0,1064e-9/0.0254,'WFS4',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+%-----------------------------------------------------------------------------
+% AS
+%-----------------------------------------------------------------------------
+
+% AS in vac
+x=0; y=59;
+resolution=5000;
+
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'BS');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=172.85;
+[h,x,y]=beamplot(qBSRFreturning,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT1,0,x,y,0,1064e-9/0.0254,'ASOT1');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=52.2;
+[h,x,y]=beamplot(qASOT1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT2,0,x,y,0,1064e-9/0.0254,'ASOT2',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=76.97;
+[h,x,y]=beamplot(qASOT2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_faraday,0,x,y,0,1064e-9/0.0254,'Faraday',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=84.09;
+[h,x,y]=beamplot(qASFaraday,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'viewport',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% AS Table, LSC (detect) path
+xlsc=x; ylsc=y;
+
+propdist=86;
+[h,x,y]=beamplot(qASviewport,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+xlsc=x; ylsc=y;
+
+propdist=53.05;
+[h,x,y]=beamplot(qASL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'229 mm',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=20.75;
+[h,x,y]=beamplot(qASL2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2/25.4,0,x,y,0,1064e-9/0.0254,'Detection PD',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% AS Table, Acquire path
+x=xlsc; y=ylsc-4;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=40.59;
+[h,x,y]=beamplot(qASL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2/25.4,0,x,y,0,1064e-9/0.0254,'Acquisition PD',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% AS Table, Camera path
+
+x=xlsc; y=ylsc-8;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=53.8;
+[h,x,y]=beamplot(qASL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(5/25.4,0,x,y,0,1064e-9/0.0254,'Camera',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+% AS Table, WFS1 path
+
+x=xlsc; y=ylsc-12;
+
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=39.55;
+[h,x,y]=beamplot(qASL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'-76 mm',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=14.25;
+[h,x,y]=beamplot(qASL3,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(10/25.4,0,x,y,0,1064e-9/0.0254,'WFS1',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%-----------------------------------------------------------------------------
+% POB
+%-----------------------------------------------------------------------------
+
+% POB in vac
+x=0; y=38;
+resolution=1000;
+
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'BS');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=525.97;
+[h,x,y]=beamplot(qBSFFexiting,[0:propdist/resolution:propdist],x,y,0); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT1,0,x,y,0,1064e-9/0.0254,'POBOT1');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=52.2;
+[h,x,y]=beamplot(qPOBOT1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT2,0,x,y,0,1064e-9/0.0254,'POBOT2',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=133.85;
+[h,x,y]=beamplot(qPOBOT2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'viewport',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+% POB lsc path
+
+propdist=82;
+[h,x,y]=beamplot(qPOBviewport,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'0.69 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+xlsc=x; ylsc=y;
+
+propdist=26.75;
+[h,x,y]=beamplot(qPOBL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2/25.4,0,x,y,0,1064e-9/0.0254,'POB PD',2);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+% POB SPOB path
+
+x=xlsc; y=ylsc-4;
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'0.57 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=38.5;
+[h,x,y]=beamplot(qPOBL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'50 mm',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=2.7;
+[h,x,y]=beamplot(qPOBL2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(0.1/25.4,0,x,y,0,1064e-9/0.0254,'1811',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%-----------------------------------------------------------------------------
+% POX
+%-----------------------------------------------------------------------------
+
+x=0; y=24;
+resolution=1000;
+
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'ITMX');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=358.77;
+[h,x,y]=beamplot(qITMXRF,[0:propdist/resolution:propdist],x,y,0); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT1,0,x,y,0,1064e-9/0.0254,'POXOT1');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=52.2;
+[h,x,y]=beamplot(qPOXOT1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT2,0,x,y,0,1064e-9/0.0254,'POXOT2',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=150.35;
+[h,x,y]=beamplot(qPOXOT2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'viewport',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=88.5;
+[h,x,y]=beamplot(qPOXviewport,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'0.69 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=27;
+[h,x,y]=beamplot(qPOXL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2/25.4,0,x,y,0,1064e-9/0.0254,'POX PD',2);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%-----------------------------------------------------------------------------
+% POY
+%-----------------------------------------------------------------------------
+
+% POY in vac
+x=0; y=10;
+resolution=1000;
+
+[h,th]=lensplot(d_largeoptic,0,x,y,0,1064e-9/0.0254,'ITMY');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=326.22;
+[h,x,y]=beamplot(qITMYRF,[0:propdist/resolution:propdist],x,y,0); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT1,0,x,y,0,1064e-9/0.0254,'POYOT1');
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=52.2;
+[h,x,y]=beamplot(qPOYOT1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(d_OT2,0,x,y,0,1064e-9/0.0254,'POYOT2',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=95.97;
+[h,x,y]=beamplot(qPOYOT2,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'viewport',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% POY Table, common path
+
+propdist=41;
+[h,x,y]=beamplot(qPOYviewport,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+xcommon=x; ycommon=y;
+
+%POY Table, WFS2 path
+x=xcommon; y=ycommon;
+
+%propdist=40.25;
+%[h,x,y]=beamplot(qPOYL1,[0:propdist/resolution:propdist],x,y); hold on;
+%set(h,'Color',linecolor);
+%set(h,'LineStyle',linestyle);
+%[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'-76.3 mm',1);
+%set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%propdist=29.7;
+%[h,x,y]=beamplot(qPOYL2,[0:propdist/resolution:propdist],x,y); hold on;
+%set(h,'Color',linecolor);
+%set(h,'LineStyle',linestyle);
+%[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'~220 mm',0);
+%set(th,'FontUnits','normalized','FontSize',fontsize);
+
+propdist=72.95;
+[h,x,y]=beamplot(qPOYL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(10/25.4,0,x,y,0,1064e-9/0.0254,'WFS2',2);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% POY Table, video camera path
+
+x=xcommon; y=ycommon-4;
+
+%[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'-76.3 mm',0);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%propdist=28.75;
+propdist=57.75;
+%[h,x,y]=beamplot(qPOYL2,[0:propdist/resolution:propdist],x,y); hold on;
+[h,x,y]=beamplot(qPOYL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(5/25.4,0,x,y,0,1064e-9/0.0254,'Camera',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+% POY Table, phasecamera path
+x=xcommon; y=ycommon-8;
+
+%[h,th]=lensplot(1,0,x,y,0,1064e-9/0.0254,'-76.3 mm',0);
+[h,th]=lensplot(2,0,x,y,0,1064e-9/0.0254,'1.1 m',0);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+%propdist=28.75;
+propdist=69;
+%[h,x,y]=beamplot(qPOYL2,[0:propdist/resolution:propdist],x,y); hold on;
+[h,x,y]=beamplot(qPOYL1,[0:propdist/resolution:propdist],x,y); hold on;
+set(h,'Color',linecolor);
+set(h,'LineStyle',linestyle);
+[h,th]=lensplot(0.1/25.4,0,x,y,0,1064e-9/0.0254,'1811',1);
+set(th,'FontUnits','normalized','FontSize',fontsize);
+
+
+end
+
+[h,th]=scalemarker(1800,15,393.70,'10 m','h');
+set(th,'FontUnits','normalized','FontSize',0.03);
+[h,th]=scalemarker(2000,20,3.9370,'0.1 m','v');
+set(th,'FontUnits','normalized','FontSize',0.03);
+set(gcf,'Color','w');
+set(gcf,'NumberTitle','off');
+set(gcf,'Name','LLO beams');
+xlabel('Inches');
+ylabel('Inches');
+set(gca,'Visible','on');
+set(gca,'YLim',[0 95]);
+set(gca,'XLim',[-100,2400]);
+set(gca,'Box','on');
+set(gca,'Visible','on');
+%title('Input and output beams');
+hold off;
+
|