-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathRecirculatingLoop.m
262 lines (217 loc) · 10.4 KB
/
RecirculatingLoop.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
classdef RecirculatingLoop < handle
%RECIRCULATINGLOOP Recirculating Optical Fiber Loop
% This class defines the RecirculatingLoop object, which
% characterizes a recirculating optical fiber loop.
% Copyright (c) 2018 Dario Pilori, Politecnico di Torino <[email protected]>
% SPDX-License-Identifier: MIT
%% Read-only properties
properties (SetAccess = private, GetAccess = public)
tloop % propagation delay of each loop (seconds)
max_loops % maximum number of loops to measure (integer)
cur_loop % current loop
aom_ppg % GPIB object of PPG that controls AOMs
trigger_ppg % GPIB object of PPG that controls trigger
initialized = false % true if loop has been initialized successfully
scrambler_delay = 0 % polarization scrambler delay (seconds)
end
%% Public properties
properties
start_delay = 0.25 % trigger start delay (fraction of tloop)
stop_advance = 0.15 % trigger stop advance (fraction of tloop)
end
%% Private properties
properties (Access = private)
tfill % time to fill the loop (seconds)
f_trig % trigger frequency (Hz)
end
%% Methods
methods
%% Constructor
function obj = RecirculatingLoop(addr1,addr2,tloop,max_loops)
%RECIRCULATINGLOOP Initializes the loop
% This function initializes a new recirculation loop
% Parameters:
% addr1 := GPIB address of first DG535 (loop+fill switch)
% addr2 := GPIB address of second DG535 (trigger)
% tloop := Propagation delay of loop (c/neff*L, seconds )
% n_loops := Number of loops to propagate
validateattributes(addr1,{'numeric'},{'scalar','nonnegative','integer'})
validateattributes(addr2,{'numeric'},{'scalar','nonnegative','integer'})
validateattributes(tloop,{'numeric'},{'scalar','positive'})
validateattributes(max_loops,{'numeric'},{'scalar','positive','integer'})
obj.initialized = false; % at beginning loop is not initialized
obj.tloop = tloop;
obj.max_loops = max_loops;
obj.tfill = 1.5*obj.tloop; % time to fill the loop
obj.f_trig = 1/(obj.tfill+max_loops*obj.tloop); % trigger time
obj.cur_loop = 1; % select first loop by default
obj.aom_ppg = gpib('ni',0,addr1); % set first PPG GPIB
obj.trigger_ppg = gpib('ni',0,addr2); % set second PPG GPIB
end
%% Initialize loop
function obj = Initialize(obj)
%INITIALIZE Initialize loop
% This function initializes the loop, re-setting the two PPGs
% from scratch. Must be run the first time.
if (obj.start_delay+obj.stop_advance)>=1
error('start_delay and stop_advance are too large');
end
%% AOMs PPG
% Open connection
fopen(obj.aom_ppg);
fprintf(obj.aom_ppg,'CL'); % clear instrument
% Set to high-impedance load
fprintf(obj.aom_ppg,'TZ 1,1');
fprintf(obj.aom_ppg,'TZ 4,1');
fprintf(obj.aom_ppg,'TZ 7,1');
% Set TTL output
fprintf(obj.aom_ppg,'OM 1,0');
fprintf(obj.aom_ppg,'OM 4,0');
fprintf(obj.aom_ppg,'OM 7,0');
% Set trigger
fprintf(obj.aom_ppg,'TM 0'); % internal trigger
fprintf(obj.aom_ppg,['TR 0,',num2str(obj.f_trig,'%f')]); % set trigger frequency
% Set delays
fprintf(obj.aom_ppg,'DT 2,1,0');
fprintf(obj.aom_ppg,['DT 3,2,',num2str(obj.tfill,'%E')]);
fprintf(obj.aom_ppg,'DT 5,2,0');
fprintf(obj.aom_ppg,'DT 6,3,0');
% Close connection
fclose(obj.aom_ppg);
%% Trigger PPG
fopen(obj.trigger_ppg);
fprintf(obj.trigger_ppg,'CL'); % clear instrument
% Set to high-impedance load
fprintf(obj.trigger_ppg,'TZ 1,1');
fprintf(obj.trigger_ppg,'TZ 4,1');
% Set TTL output
fprintf(obj.trigger_ppg,'OM 1,0');
fprintf(obj.trigger_ppg,'OM 4,0');
% Set trigger
fprintf(obj.trigger_ppg,'TM 1'); % Set external trigger
fprintf(obj.trigger_ppg,'TZ 0,1'); % set trigger impedance
% Set delays
fprintf(obj.trigger_ppg,['DT 2,1,',num2str(obj.tfill+...
(obj.cur_loop-1+obj.start_delay)*obj.tloop,'%E')]);
fprintf(obj.trigger_ppg,['DT 3,2,',...
num2str(obj.tloop*(1-obj.start_delay-obj.stop_advance),'%E')]);
fprintf(obj.trigger_ppg,['DT 5,1,',num2str(obj.tfill-obj.scrambler_delay,'%E')]);
fprintf(obj.trigger_ppg,['DT 6,5,',num2str(100e-6,'%E')]);
% Close connection
fclose(obj.trigger_ppg);
%% Set to initialized
obj.initialized = true;
%% Suppress output
if nargout == 0
clear obj
end
end
%% Reconfigure loop
function obj = ReconfigureLoop(obj,tloop,max_loops)
%RECONFIGURELOOP Change loop parameters
% Use this function to change the loop params (tloop and
% max_loops)
validateattributes(tloop,{'numeric'},{'scalar','positive'})
validateattributes(max_loops,{'numeric'},{'scalar','positive','integer'})
%% Check if loop is OK
if ~obj.initialized
error('Loop not initialized; please initialize');
end
if (obj.start_delay+obj.stop_advance)>=1
error('start_delay and stop_advance are too large');
end
%% Calculate parameters
obj.initialized = false;
obj.tloop = tloop;
obj.max_loops = max_loops;
obj.tfill = 1.5*obj.tloop; % time to fill the loop
obj.f_trig = 1/(obj.tfill+max_loops*obj.tloop); % trigger time
obj.cur_loop = 1; % select first loop by default
%% Set up first PPG
fopen(obj.aom_ppg);
fprintf(obj.aom_ppg,['TR 0,',num2str(obj.f_trig,'%f')]); % set trigger frequency
fprintf(obj.aom_ppg,'DT 2,1,0');
fprintf(obj.aom_ppg,['DT 3,2,',num2str(obj.tfill,'%E')]);
fprintf(obj.aom_ppg,'DT 5,2,0');
fprintf(obj.aom_ppg,'DT 6,3,0');
fclose(obj.aom_ppg);
%% Set up second PPG
fopen(obj.trigger_ppg);
fprintf(obj.trigger_ppg,['DT 2,1,',num2str(obj.tfill+...
(obj.cur_loop-1+obj.start_delay)*obj.tloop,'%E')]);
fprintf(obj.trigger_ppg,['DT 3,2,',...
num2str(obj.tloop*(1-obj.start_delay-obj.stop_advance),'%E')]);
fprintf(obj.trigger_ppg,['DT 5,1,',num2str(obj.tfill-obj.scrambler_delay,'%E')]);
fprintf(obj.trigger_ppg,['DT 6,5,',num2str(100e-6,'%E')]);
fclose(obj.trigger_ppg);
obj.initialized = true;
%% Suppress output
if nargout == 0
clear obj
end
end
%% Select loop
function obj = SelectLoop(obj,cur_loop)
%SELECTLOOP Select loop
% Use this function to select a loop between 0 (B2B) to
% max_loops
validateattributes(cur_loop,{'numeric'},{'scalar','nonnegative',...
'<=',obj.max_loops})
%% Check
if ~obj.initialized
error('Loop not initialized; please initialize');
end
if (obj.start_delay+obj.stop_advance)>=1
error('start_delay and stop_advance are too large');
end
%% Calculate parameters
obj.cur_loop = cur_loop; % select first loop by default
%% Set up second PPG
fopen(obj.trigger_ppg);
fprintf(obj.trigger_ppg,['DT 2,1,',num2str(obj.tfill+...
(obj.cur_loop-1+obj.start_delay)*obj.tloop,'%E')]);
fprintf(obj.trigger_ppg,['DT 3,2,',...
num2str(obj.tloop*(1-obj.start_delay-obj.stop_advance),'%E')]);
fclose(obj.trigger_ppg);
%% Suppress output
if nargout == 0
clear obj
end
end
%% Change scrambler delay
function obj = SetScramblerDelay(obj,value)
%SCRAMBLER_DELAY Set scrambler delay
% Use this function to change the time that the polarization
% scrambler needs to react to a trigger
validateattributes(value,{'numeric'},{'scalar','nonnegative'});
%% Set delay on PPG
fopen(obj.trigger_ppg);
fprintf(obj.trigger_ppg,['DT 5,1,',num2str(obj.tfill-value,'%E')]);
fprintf(obj.trigger_ppg,['DT 6,5,',num2str(100e-6,'%E')]);
fclose(obj.trigger_ppg);
%% Set property
obj.scrambler_delay = value;
%% Suppress output
if nargout == 0
clear obj
end
end
%% Set property methods
% Check start_delay
function set.start_delay(obj, value)
validateattributes(value,{'numeric'},{'scalar','>=',0,'<=',1});
obj.start_delay = value;
end
% Check stop advance
function set.stop_advance(obj, value)
validateattributes(value,{'numeric'},{'scalar','>=',0,'<=',1});
obj.stop_advance = value;
end
%% Destructor
function delete(obj)
% Delete handles to PPG
delete(obj.trigger_ppg);
delete(obj.aom_ppg);
end
end
end