DEFINE INSTRUMENT Guide_P(sa_pos=0.5,fp=0.5)
   
DECLARE
%{
 

 double R=0.86;
 double Qc=0.0217;
 double R0=0.995;
 double W=1.0/350.0;
 double M=2.5;
 double alpha;
 
 
%}

INITIALIZE
%{

  alpha=(R0-R)/Qc/(M-1);

%}

TRACE

COMPONENT a1 = Arm()
  AT (0,0,0) ABSOLUTE

COMPONENT csource_PSI = Source_gen(
   h = 0.135, w = 0.08, dist=0.5, xw = 0.05, yh = 0.15,
   Lmin = 1.0, Lmax = 10.0, T1= 296.16, I1=8.5E11,
   T2=40.68, I2=5.2E11)
AT (0,0,0) RELATIVE a1  ROTATED (0,0,0) RELATIVE a1


COMPONENT Lmonsource1 = L_monitor(
    nchan = 200, filename = "L_source.dat", xwidth=0.05,
    yheight = 0.15, Lmin = 2.0, Lmax = 10.0)
  AT (0, 0, 0.5) RELATIVE a1

COMPONENT PSD_source = PSD_monitor(
    nx = 100, ny = 100, filename = "PSD_source.dat",
    xwidth = 0.05, yheight = 0.2)
  AT (0, 0, 0.5) RELATIVE a1

COMPONENT ell_guide1 = Guide_tapering(
     w1 = 0.05, h1=0.15,
     linw = 0.0, loutw = fp, l=36.0, linh=0.0,
     louth = fp, option="parabolical", R0 = R0, 
     Qcx = 0.0217, Qcy = 0.0217,alphax = alpha, 
     alphay = alpha, W = W,
     mx = M, my = M, segno=200)
AT (0,0,0.5) RELATIVE a1 ROTATED (0,0,0) RELATIVE a1


COMPONENT GuideEnd = Slit( 
        xmin = -0.1 , xmax = 0.1,
        ymin = -0.1, ymax = 0.1)
  AT (0,0,36.5) RELATIVE a1 ROTATED (0,0,0) RELATIVE a1

COMPONENT sa_arm = Arm()
   AT (0, 0, sa_pos) RELATIVE GuideEnd 

COMPONENT PSD_sa1 = PSD_monitor(
        xmin = -0.01, xmax = 0.01, ymin = - 0.01, ymax = 0.01,
        nx=50, ny=50, filename="PSD_sa1.dat")
  AT (0, 0, 0.0001) RELATIVE sa_arm

COMPONENT PSD_sa2 = PSD_monitor(
        xmin = -0.005, xmax = 0.005, ymin = - 0.005, ymax = 0.005,
        nx=50, ny=50, filename="PSD_sa2.dat")
  AT (0, 0, 0.0001) RELATIVE PREVIOUS

COMPONENT PSD_full = PSD_monitor(
        xmin = -0.1, xmax = 0.1, ymin = - 0.1, ymax = 0.1,
        nx=50, ny=50, filename="PSD_full.dat")
  AT (0, 0, 0.0001) RELATIVE PREVIOUS

COMPONENT div_sa = Divergence_monitor(
    nh = 100, nv = 100, filename = "Div_sa.dat", xmin = - 0.04,
    xmax = 0.04, ymin = -0.04, ymax = 0.04, h_maxdiv = 3,
    v_maxdiv = 3)
  AT (0, 0, 0.0001) RELATIVE PREVIOUS

COMPONENT Lambda_guide = L_monitor(
    nchan = 100, filename = "lambda_sa.dat", xwidth = 0.05,
    yheight = 0.15, Lmin = 1.5, Lmax = 11.0)
  AT (0, 0, 0.0001) RELATIVE PREVIOUS


END