#-----------------------------------------------------------------------
# An Advanced Virgo Plus input file for Finesse 3

# Source: current finesse-virgo's model
# Modified with: https://git.ligo.org/virgo/isc/parameters/-/blob/377982cc0fd35a561a2f0d9356ca847072d39c57/mechanics.m
# Modified with: https://git.ligo.org/virgo/isc/parameters/-/blob/c7e2cf038a096fbb5a20cc58ff01650cec5a96fa/advanced_virgo.m
# Modified with: https://git.ligo.org/virgo/isc/parameters/-/blob/7de842d4f3ebe3744a608e2f19d6d0a15a437cdf/advanced_virgo_current.m
# Modified with: https://git.ligo.org/virgo/isc/parameters/-/blob/a89dfb3b8556c9796b64531da93343f56e773e6b/detection_current.m
# Modified with: https://git.ligo.org/virgo/isc/optickle/-/blob/f577e6f073c0dda2ef9db3d2fbab27021f406182/models/DRITFscatter.m
# Modified with: https://git.ligo.org/virgo/virgoapp/simulinkNoiseBudget/-/blob/6e7e1b95b7817cdb6cc67d5748281ed1d7ffc805/DRITF/scatterRadiationPressureTune.m
# Modified with: https://git.ligo.org/virgo/virgoapp/simulinkNoiseBudget/-/blob/6e7e1b95b7817cdb6cc67d5748281ed1d7ffc805/DRITF/scatterRadiationPressure.m

# Defines general variables needed for the model creation here

# Constant
variable	speed_of_light	299792458

# Material
variable	silica	1.44963

# Injection
variable	wavelength	1.064e-6	units='m'
variable	frequency	speed_of_light/wavelength	units='Hz'

# Schnupp asymmetry
variable	schnupp_asymetry	0.23	units='m'

# Cavity length
variable	length_PRCL	11.952	units='m'
variable	length_SRCL	11.952	units='m'

# Laser
laser	laser	P=25	f=frequency

# Power Recycling
mirror	PR	T=0.04835	L=37.5e-6	Rc=-1430

space	space_laser_PR	laser.p1	PR.p1	L=11.67

# Beamsplitter

beamsplitter	BS	T=0.5012	L=37.5e-6	alpha=-45.0

space	space_PR_BS	PR.p2	BS.p1	L=5.925

# North arm

mirror	NI	T=0.01377	L=37.5e-6	Rc=-1420
mirror	NE	T=4.4e-6	L=37.5e-6	Rc=1683

space	space_NI_NE	NI.p2	NE.p1	L=2999.8

space	space_BS_NI	BS.p3	NI.p1	L=length_PRCL-space_PR_BS.L+schnupp_asymetry/2

# West arm

mirror	WI	T=0.01377	L=37.5e-6	Rc=-1420
mirror	WE	T=4.3e-6	L=37.5e-6	Rc=1683

space	space_WI_WE	WI.p2	WE.p1	L=2999.8

space	space_BS_WI	BS.p2	WI.p1	L=length_PRCL-space_PR_BS.L-schnupp_asymetry/2

# SR

mirror	SR	T=0.40	L=37.5e-6	Rc=1430

space	space_BS_SR	BS.p4	SR.p1	L=length_SRCL-(space_BS_WI.L+space_BS_NI.L)/2

# power detectors

power_detector_dc	NE_p1	NE.p1.o
power_detector_dc	NI_p1	NI.p1.o
power_detector_dc	NE_p2	NE.p2.o
power_detector_dc	NI_p2	NI.p2.o

power_detector_dc	WE_p1	WE.p1.o
power_detector_dc	WI_p1	WI.p1.o
power_detector_dc	WE_p2	WE.p2.o
power_detector_dc	WI_p2	WI.p2.o

power_detector_dc	PR_p1	PR.p1.o
power_detector_dc	PR_p2	PR.p2.o

power_detector_dc	SR_p1	SR.p1.o
power_detector_dc	SR_p2	SR.p2.o

# degree of freedom

degree_of_freedom	NORTH_ARM	NE.dofs.z	+1
degree_of_freedom	WEST_ARM	WE.dofs.z	+1
degree_of_freedom	PRCL	PR.dofs.z	+1
degree_of_freedom	SRCL	SR.dofs.z	+1

degree_of_freedom MICH	NI.dofs.z	-1	WI.dofs.z	+1	NE.dofs.z	-1	WE.dofs.z	+1
degree_of_freedom DARM	NE.dofs.z	-1	WE.dofs.z	+1
degree_of_freedom CARM	NE.dofs.z	+1	WE.dofs.z	+1
degree_of_freedom	DARM_Fz	NE.dofs.F_z	-1	WE.dofs.F_z	+1

# cavities

cavity	PRCL_north	PR.p2.o	via=WE.p1.i	priority=2
cavity	PRCL_west	PR.p2.o	via=NE.p1.i	priority=2

cavity	SRCL_north	SR.p1.o	via=WE.p1.i	priority=1
cavity	SRCL_west	SR.p1.o	via=WE.p1.i	priority=1

readout_rf	B1	SR.p2.o	output_detectors=True