source("abcd.m")
########################################## 
lambda= 1.064E-6 
Ltot= 0.55
r0= 1e100
w0= 0.00675
f1=0.1526;
lns1.abcd=abcd_lens( f1) ;
lns1.x= 0.2;
f2=0.019;
lns2.abcd=abcd_lens( f2) ;
lns2.x= lns1.x+ f1 + f2;
lns3.abcd=abcd_lens( 0.0629) ;
lns3.x= 0.500;
# befrore we calculate final waste we need to take in account propagation through the thick mirror
# since we try to match to the waste right in front of front coupling mirror
r_at_mirror=.01;
w_at_mirror=83.729e-6;
q_at_mirror=waste_r2q(w_at_mirror, r_at_mirror, lambda);
n_fused_cilica=1.44968;
R_of_mirror=.01;
mirror_thickness=6.35e-3;
f_of_mirrr_formed_lense=1/((n_fused_cilica-1)*(1/(-R_of_mirror)));
abcd_of_mirrr_formed_lense=abcd_lens(f_of_mirrr_formed_lense);
abcd_of_propagation_through_mirror=abcd_free_space(mirror_thickness/n_fused_cilica);
abcd_of_thick_mirror=abcd_of_propagation_through_mirror*abcd_of_mirrr_formed_lense;
# q right in front of front coupling mirror of the cavity 
q_infront_of_mirror=q_after_element(q_at_mirror, abcd_of_thick_mirror);
wf=q2waste(q_infront_of_mirror,lambda)
# we need minus for rf since we will solve injection to the cavity problem
rf= - q2radius(q_infront_of_mirror,lambda)
##########################################