################################################################################
                            Project TBASEDSB030110
################################################################################
 
General Comments:

IF tests on optical links

Weather: light snow; Temp 25 F; high winds (shutdown most of time)

IF noise --> ODM 2 --> CM 1,3   --> SF 1,3 --> AC sampler 0,4
IF noise --> ODM 4 --> CM 5,7   --> SF 5,7 --> AC sampler 1,5
IF noise --> ODM 6 --> CM 10,12 --> SF 2,4 --> AC sampler 2,6
IF noise --> ODM 8 --> CM 14,16 --> SF 6,8 --> AC sampler 3,7

ODM 2: original external modulator
ODM 4: out of service
ODM 6: new external modulator
ODM 8: direct modulator

Bank A: 1W2-001-800  (ch1 -> odm 2, ch2 -> odm 4)
Bank B: 1W2-123-800  (ch1 -> odm 6, ch2 -> odm 8)
Bank C: 1W2-245-800  (ch1 -> odm 2, ch2 -> odm 4)
Bank D: 1W2-367-800  (ch1 -> odm 6, ch2 -> odm 8)

Note: Initial problems getting the spectrometer to download software.
This was caused by the computer switch being set to a debugging mode
(switch flipped to 167 instead of sparc in the front panel).

Note: The telescope was staitionary during scan 2 except for the
following periods:

22:35:53 - 22:42:05  moved from survival to snow dump
18:05:15 - 18:10:05  moved from snow dump to survival


Ant:    1 if_noise             RA 2.92   Dec 54.99  Az 44.15  El 60.10 
AC :    1 2003_01_10_22:24:22A.fits      SigRef  0 0     Cal  0 1    
          Mode 1W2-0XY-800     Samplers 0 1
          BW 800 Ch 2048   Rx 2 Phase 2  Rec 4  
Comments: if_noise; test
--------------------------------------------------------------------------------
AC :    2 2003_01_10_22:39:47A.fits      SigRef  0 0     Cal  0 1    
          Mode 1W2-0XY-800     Samplers 0 1
          BW 800 Ch 2048   Rx 2 Phase 2  Rec 402
Comments: if_noise; duration 201 min
--------------------------------------------------------------------------------

================================================================================
                            Data Reduction [dsb]
================================================================================

# General scales for Fourier transform plots
s.pg.setyscale(-0.01, 0.2);
s.pg.setxscale(0.3, 0.55);
s.pg.setxscale(0.405, 0.426);


- Figure acf.ps
y := s.getAcf(proj, 2, bank=1, rcvr=1, rec_num=[1:100])
s.plot(ind(y), y, title='if_noise  Scan 2  Rec 1:100', ytitle='ACF', xtitle='Lags')
s.plotAcf(proj, 2, bank=4, rcvr=1, rec_num=[1:100], newplot=F, linecolor=3)         
s.plotAcf(proj, 2, bank=4, rcvr=2, rec_num=[1:100], newplot=F, linecolor=4)   


o Bank A, Rcvr 1

- Figure odm2_scans2-50_ft.ps, odm2_scans2-50_ft2.ps
y2 := s.getSR1(proj, 2, 1, 2, bank=1, rcvr=1);
y10 := s.getSR1(proj, 2, 1, 10, bank=1, rcvr=1);
y20 := s.getSR1(proj, 2, 1, 20, bank=1, rcvr=1);
y30 := s.getSR1(proj, 2, 1, 30, bank=1, rcvr=1);
y40 := s.getSR1(proj, 2, 1, 40, bank=1, rcvr=1);
y50 := s.getSR1(proj, 2, 1, 50, bank=1, rcvr=1);
f2 := s.ft(y2.sr[200:1848], 0.391);
f10 := s.ft(y10.sr[200:1848], 0.391);
f20 := s.ft(y20.sr[200:1848], 0.391);
f30 := s.ft(y30.sr[200:1848], 0.391);
f40 := s.ft(y40.sr[200:1848], 0.391);
f50 := s.ft(y50.sr[200:1848], 0.391);
s.plot(f2.x, f2.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 1  rcvr 1  rec 2');
s.plot(f10.x, f10.y, newplot=F, linecolor=3);
s.plot(f20.x, f20.y, newplot=F, linecolor=4);
s.plot(f30.x, f30.y, newplot=F, linecolor=5); 
s.plot(f40.x, f40.y, newplot=F, linecolor=6); 
s.plot(f50.x, f50.y, newplot=F, linecolor=7);

- Figure odm2_scans2-400_ft.ps, odm2_scans2-400_ft2.ps
y2 := s.getSR1(proj, 2, 1, 2, bank=1, rcvr=1);
y50 := s.getSR1(proj, 2, 1, 50, bank=1, rcvr=1);
y100 := s.getSR1(proj, 2, 1, 100, bank=1, rcvr=1);
y200 := s.getSR1(proj, 2, 1, 200, bank=1, rcvr=1);
y300 := s.getSR1(proj, 2, 1, 300, bank=1, rcvr=1);
y400 := s.getSR1(proj, 2, 1, 400, bank=1, rcvr=1);
f2 := s.ft(y2.sr[200:1848], 0.391);
f50 := s.ft(y50.sr[200:1848], 0.391);
f100 := s.ft(y100.sr[200:1848], 0.391);
f200 := s.ft(y200.sr[200:1848], 0.391);
f300 := s.ft(y300.sr[200:1848], 0.391);
f400 := s.ft(y400.sr[200:1848], 0.391);
s.plot(f2.x, f2.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 1  rcvr 1  rec 2');
s.plot(f50.x, f50.y, newplot=F, linecolor=3);
s.plot(f100.x, f100.y, newplot=F, linecolor=4);
s.plot(f200.x, f200.y, newplot=F, linecolor=5);
s.plot(f300.x, f300.y, newplot=F, linecolor=6);
s.plot(f400.x, f400.y, newplot=F, linecolor=7);

- Figure odm2_avg1-10_ft.ps, odm2_avg1-10_ft2.ps 
y1 := s.getSR1(proj, 2, 1, 2, bank=1, rcvr=1) 
f1 := s.ft(y1.sr[200:1848], 0.391); 
s.plot(f1.x, f1.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 1  rcvr 1')
y4 := s.getSR1(proj, 2, [1:4], [5:8], bank=1, rcvr=1)
f4 := s.ft(y4.sr[200:1848], 0.391);                   
s.plot(f4.x, f4.y, newplot=F, linecolor=3)
y8 := s.getSR1(proj, 2, [1:8], [9:16], bank=1, rcvr=1)
f8 := s.ft(y8.sr[200:1848], 0.391);                   
s.plot(f8.x, f8.y, newplot=F, linecolor=4)
y10 := s.getSR1(proj, 2, [1:10], [11:20], bank=1, rcvr=1)
f10 := s.ft(y10.sr[200:1848], 0.391);                   
s.plot(f10.x, f10.y, newplot=F, linecolor=5)

- Figure odm2_avg10-50_ft.ps, odm2_avg10-50_ft2.ps
y10 := s.getSR1(proj, 2, [1:10], [11:20], bank=1, rcvr=1)
f10 := s.ft(y10.sr[200:1848], 0.391);
s.plot(f10.x, f10.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 1  rcvr 1')
y20 := s.getSR1(proj, 2, [1:20], [21:40], bank=1, rcvr=1)
f20 := s.ft(y20.sr[200:1848], 0.391);                   
s.plot(f20.x, f20.y, newplot=F, linecolor=3)
y30 := s.getSR1(proj, 2, [1:30], [31:60], bank=1, rcvr=1)
f30 := s.ft(y30.sr[200:1848], 0.391);                   
s.plot(f30.x, f30.y, newplot=F, linecolor=4)
y40 := s.getSR1(proj, 2, [1:40], [41:80], bank=1, rcvr=1)
f40 := s.ft(y40.sr[200:1848], 0.391);                   
s.plot(f40.x, f40.y, newplot=F, linecolor=5)
y50 := s.getSR1(proj, 2, [1:50], [51:100], bank=1, rcvr=1)
f50 := s.ft(y50.sr[200:1848], 0.391);                   
s.plot(f50.x, f50.y, newplot=F, linecolor=6)

- Figure odm2_5grp_spec.ps, odm2_5grp_ft.ps
r1 := [11:20, 31:40, 51:60, 71:80, 91:100];
r2 := [111:120, 131:140, 151:160, 171:180, 191:200];
r3 := [211:220, 231:240, 251:260, 271:280, 291:300];
r4 := [311:320, 331:340, 351:360, 371:380, 391:400];
s1 := [1:10, 21:30, 41:50, 61:70, 81:90];
s2 := [101:110, 121:130, 141:150, 161:170, 181:190];
s3 := [201:210, 221:230, 241:250, 261:270, 281:290];
s4 := [301:310, 321:330, 341:350, 361:370, 381:390];
y1 := s.getSR1(proj, 2, r1, s1, bank=1, rcvr=1);
y2 := s.getSR1(proj, 2, r2, s2, bank=1, rcvr=1);
y3 := s.getSR1(proj, 2, r3, s3, bank=1, rcvr=1);
y4 := s.getSR1(proj, 2, r4, s4, bank=1, rcvr=1);
f1 := s.ft(y1.sr[200:1848], 0.391);
f2 := s.ft(y2.sr[200:1848], 0.391);
f3 := s.ft(y3.sr[200:1848], 0.391);
f4 := s.ft(y4.sr[200:1848], 0.391);
s.plot(y1.freq, y1.sr, xtitle='Frequency [MHz]', ytitle='(Sig-Ref)/Ref', title='IF Noise  Scan 2  Bank 1  rcvr 1', drop=200);
s.plot(y2.freq, y2.sr, newplot=F, linecolor=3);
s.plot(y3.freq, y3.sr, newplot=F, linecolor=4);
s.plot(y4.freq, y4.sr, newplot=F, linecolor=5);
s.pg.setyscale(-0.0015, 0.0015)
s.plot(f1.x, f1.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 1  rcvr 1');
s.plot(f2.x, f2.y, newplot=F, linecolor=3);
s.plot(f3.x, f3.y, newplot=F, linecolor=4);
s.plot(f4.x, f4.y, newplot=F, linecolor=5);
s.pg.setyscale(-0.001, 0.015);
s.pg.setxscale(0.3, 0.55);

- Figure odm2_20grp_spec.ps, odm2_20grp_ft.ps
y := y1;
y.sr := (y1.sr + y2.sr + y3.sr + y4.sr)/4.0;
f := s.ft(y.sr[200:1848], 0.391);
s.plot(y.freq, y.sr, xtitle='Frequency [MHz]', ytitle='(Sig-Ref)/Ref', title='IF Noise  Scan 2  Bank 1  rcvr 1', drop=200);
s.pg.setyscale(-0.0002, 0.0002);
s.plot(f.x, f.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 1  rcvr 1');
s.pg.setyscale(-0.001, 0.015);
s.pg.setxscale(0.3, 0.55);


o Bank D, Rcvr 1

- Figure odm6_scans2-50_ft.ps
y2 := s.getSR1(proj, 2, 1, 2, bank=4, rcvr=1);
y10 := s.getSR1(proj, 2, 1, 10, bank=4, rcvr=1);
y20 := s.getSR1(proj, 2, 1, 20, bank=4, rcvr=1);
y30 := s.getSR1(proj, 2, 1, 30, bank=4, rcvr=1);
y40 := s.getSR1(proj, 2, 1, 40, bank=4, rcvr=1);
y50 := s.getSR1(proj, 2, 1, 50, bank=4, rcvr=1);
f2 := s.ft(y2.sr[200:1848], 0.391);
f10 := s.ft(y10.sr[200:1848], 0.391);
f20 := s.ft(y20.sr[200:1848], 0.391);
f30 := s.ft(y30.sr[200:1848], 0.391);
f40 := s.ft(y40.sr[200:1848], 0.391);
f50 := s.ft(y50.sr[200:1848], 0.391);
s.plot(f2.x, f2.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 1  rec 2');
s.plot(f10.x, f10.y, newplot=F, linecolor=3);
s.plot(f20.x, f20.y, newplot=F, linecolor=4);
s.plot(f30.x, f30.y, newplot=F, linecolor=5); 
s.plot(f40.x, f40.y, newplot=F, linecolor=6); 
s.plot(f50.x, f50.y, newplot=F, linecolor=7);

- Figure odm6_scans2-400_ft.ps
y2 := s.getSR1(proj, 2, 1, 2, bank=4, rcvr=1);
y50 := s.getSR1(proj, 2, 1, 50, bank=4, rcvr=1);
y100 := s.getSR1(proj, 2, 1, 100, bank=4, rcvr=1);
y200 := s.getSR1(proj, 2, 1, 200, bank=4, rcvr=1);
y300 := s.getSR1(proj, 2, 1, 300, bank=4, rcvr=1);
y400 := s.getSR1(proj, 2, 1, 400, bank=4, rcvr=1);
f2 := s.ft(y2.sr[200:1848], 0.391);
f50 := s.ft(y50.sr[200:1848], 0.391);
f100 := s.ft(y100.sr[200:1848], 0.391);
f200 := s.ft(y200.sr[200:1848], 0.391);
f300 := s.ft(y300.sr[200:1848], 0.391);
f400 := s.ft(y400.sr[200:1848], 0.391);
s.plot(f2.x, f2.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 1  rec 2');
s.plot(f50.x, f50.y, newplot=F, linecolor=3);
s.plot(f100.x, f100.y, newplot=F, linecolor=4);
s.plot(f200.x, f200.y, newplot=F, linecolor=5);
s.plot(f300.x, f300.y, newplot=F, linecolor=6);
s.plot(f400.x, f400.y, newplot=F, linecolor=7);

- Figure odm6_avg10-50_ft.ps, odm6_avg10-50_ft2.ps
y10 := s.getSR1(proj, 2, [1:10], [11:20], bank=4, rcvr=1)
f10 := s.ft(y10.sr[200:1848], 0.391);
s.plot(f10.x, f10.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 1')
y20 := s.getSR1(proj, 2, [1:20], [21:40], bank=4, rcvr=1)
f20 := s.ft(y20.sr[200:1848], 0.391);                   
s.plot(f20.x, f20.y, newplot=F, linecolor=3)
y30 := s.getSR1(proj, 2, [1:30], [31:60], bank=4, rcvr=1)
f30 := s.ft(y30.sr[200:1848], 0.391);                   
s.plot(f30.x, f30.y, newplot=F, linecolor=4)
y40 := s.getSR1(proj, 2, [1:40], [41:80], bank=4, rcvr=1)
f40 := s.ft(y40.sr[200:1848], 0.391);                   
s.plot(f40.x, f40.y, newplot=F, linecolor=5)
y50 := s.getSR1(proj, 2, [1:50], [51:100], bank=4, rcvr=1)
f50 := s.ft(y50.sr[200:1848], 0.391);                   
s.plot(f50.x, f50.y, newplot=F, linecolor=6)

- Figure odm6_5grp_spec.ps, odm6_5grp_ft.ps
r1 := [11:20, 31:40, 51:60, 71:80, 91:100];
r2 := [111:120, 131:140, 151:160, 171:180, 191:200];
r3 := [211:220, 231:240, 251:260, 271:280, 291:300];
r4 := [311:320, 331:340, 351:360, 371:380, 391:400];
s1 := [1:10, 21:30, 41:50, 61:70, 81:90];
s2 := [101:110, 121:130, 141:150, 161:170, 181:190];
s3 := [201:210, 221:230, 241:250, 261:270, 281:290];
s4 := [301:310, 321:330, 341:350, 361:370, 381:390];
y1 := s.getSR1(proj, 2, r1, s1, bank=4, rcvr=1);
y2 := s.getSR1(proj, 2, r2, s2, bank=4, rcvr=1);
y3 := s.getSR1(proj, 2, r3, s3, bank=4, rcvr=1);
y4 := s.getSR1(proj, 2, r4, s4, bank=4, rcvr=1);
f1 := s.ft(y1.sr[200:1848], 0.391);
f2 := s.ft(y2.sr[200:1848], 0.391);
f3 := s.ft(y3.sr[200:1848], 0.391);
f4 := s.ft(y4.sr[200:1848], 0.391);
s.plot(y1.freq, y1.sr, xtitle='Frequency [MHz]', ytitle='(Sig-Ref)/Ref', title='IF Noise  Scan 2  Bank 4  rcvr 1', drop=200);
s.plot(y2.freq, y2.sr, newplot=F, linecolor=3);
s.plot(y3.freq, y3.sr, newplot=F, linecolor=4);
s.plot(y4.freq, y4.sr, newplot=F, linecolor=5);
s.pg.setyscale(-0.0015, 0.0015)
s.plot(f1.x, f1.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 1');
s.plot(f2.x, f2.y, newplot=F, linecolor=3);
s.plot(f3.x, f3.y, newplot=F, linecolor=4);
s.plot(f4.x, f4.y, newplot=F, linecolor=5);
s.pg.setyscale(-0.001, 0.015);
s.pg.setxscale(0.3, 0.55);

- Figure odm6_20grp_spec.ps, odm6_20grp_ft.ps
y := y1;
y.sr := (y1.sr + y2.sr + y3.sr + y4.sr)/4.0;
f := s.ft(y.sr[200:1848], 0.391);
s.plot(y.freq, y.sr, xtitle='Frequency [MHz]', ytitle='(Sig-Ref)/Ref', title='IF Noise  Scan 2  Bank 4  rcvr 1', drop=200);
s.pg.setyscale(-0.0005, 0.0005);
s.plot(f.x, f.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 1');
s.pg.setyscale(-0.001, 0.015);
s.pg.setxscale(0.3, 0.55);


o Bank D, Rcvr 2

- Figure odm8_scans2-50_ft.ps
y2 := s.getSR1(proj, 2, 1, 2, bank=4, rcvr=2);
y10 := s.getSR1(proj, 2, 1, 10, bank=4, rcvr=2);
y20 := s.getSR1(proj, 2, 1, 20, bank=4, rcvr=2);
y30 := s.getSR1(proj, 2, 1, 30, bank=4, rcvr=2);
y40 := s.getSR1(proj, 2, 1, 40, bank=4, rcvr=2);
y50 := s.getSR1(proj, 2, 1, 50, bank=4, rcvr=2);
f2 := s.ft(y2.sr[200:1848], 0.391);
f10 := s.ft(y10.sr[200:1848], 0.391);
f20 := s.ft(y20.sr[200:1848], 0.391);
f30 := s.ft(y30.sr[200:1848], 0.391);
f40 := s.ft(y40.sr[200:1848], 0.391);
f50 := s.ft(y50.sr[200:1848], 0.391);
s.plot(f2.x, f2.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 2  rec 2');
s.plot(f10.x, f10.y, newplot=F, linecolor=3);
s.plot(f20.x, f20.y, newplot=F, linecolor=4);
s.plot(f30.x, f30.y, newplot=F, linecolor=5); 
s.plot(f40.x, f40.y, newplot=F, linecolor=6); 
s.plot(f50.x, f50.y, newplot=F, linecolor=7);

- Figure odm8_scans2-400_ft.ps
y2 := s.getSR1(proj, 2, 1, 2, bank=4, rcvr=2);
y50 := s.getSR1(proj, 2, 1, 50, bank=4, rcvr=2);
y100 := s.getSR1(proj, 2, 1, 100, bank=4, rcvr=2);
y200 := s.getSR1(proj, 2, 1, 200, bank=4, rcvr=2);
y300 := s.getSR1(proj, 2, 1, 300, bank=4, rcvr=2);
y400 := s.getSR1(proj, 2, 1, 400, bank=4, rcvr=2);
f2 := s.ft(y2.sr[200:1848], 0.391);
f50 := s.ft(y50.sr[200:1848], 0.391);
f100 := s.ft(y100.sr[200:1848], 0.391);
f200 := s.ft(y200.sr[200:1848], 0.391);
f300 := s.ft(y300.sr[200:1848], 0.391);
f400 := s.ft(y400.sr[200:1848], 0.391);
s.plot(f2.x, f2.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 2  rec 2');
s.plot(f50.x, f50.y, newplot=F, linecolor=3);
s.plot(f100.x, f100.y, newplot=F, linecolor=4);
s.plot(f200.x, f200.y, newplot=F, linecolor=5);
s.plot(f300.x, f300.y, newplot=F, linecolor=6);
s.plot(f400.x, f400.y, newplot=F, linecolor=7);

- Figure odm8_avg10-50_ft.ps, odm8_avg10-50_ft2.ps
y10 := s.getSR1(proj, 2, [1:10], [11:20], bank=4, rcvr=2)
f10 := s.ft(y10.sr[200:1848], 0.391);
s.plot(f10.x, f10.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 2')
y20 := s.getSR1(proj, 2, [1:20], [21:40], bank=4, rcvr=2)
f20 := s.ft(y20.sr[200:1848], 0.391);                   
s.plot(f20.x, f20.y, newplot=F, linecolor=3)
y30 := s.getSR1(proj, 2, [1:30], [31:60], bank=4, rcvr=2)
f30 := s.ft(y30.sr[200:1848], 0.391);                   
s.plot(f30.x, f30.y, newplot=F, linecolor=4)
y40 := s.getSR1(proj, 2, [1:40], [41:80], bank=4, rcvr=2)
f40 := s.ft(y40.sr[200:1848], 0.391);                   
s.plot(f40.x, f40.y, newplot=F, linecolor=5)
y50 := s.getSR1(proj, 2, [1:50], [51:100], bank=4, rcvr=2)
f50 := s.ft(y50.sr[200:1848], 0.391);                   
s.plot(f50.x, f50.y, newplot=F, linecolor=6)


- Figure odm8_5grp_spec.ps, odm8_5grp_ft.ps
r1 := [11:20, 31:40, 51:60, 71:80, 91:100];
r2 := [111:120, 131:140, 151:160, 171:180, 191:200];
r3 := [211:220, 231:240, 251:260, 271:280, 291:300];
r4 := [311:320, 331:340, 351:360, 371:380, 391:400];
s1 := [1:10, 21:30, 41:50, 61:70, 81:90];
s2 := [101:110, 121:130, 141:150, 161:170, 181:190];
s3 := [201:210, 221:230, 241:250, 261:270, 281:290];
s4 := [301:310, 321:330, 341:350, 361:370, 381:390];
y1 := s.getSR1(proj, 2, r1, s1, bank=4, rcvr=2);
y2 := s.getSR1(proj, 2, r2, s2, bank=4, rcvr=2);
y3 := s.getSR1(proj, 2, r3, s3, bank=4, rcvr=2);
y4 := s.getSR1(proj, 2, r4, s4, bank=4, rcvr=2);
f1 := s.ft(y1.sr[200:1848], 0.391);
f2 := s.ft(y2.sr[200:1848], 0.391);
f3 := s.ft(y3.sr[200:1848], 0.391);
f4 := s.ft(y4.sr[200:1848], 0.391);
s.plot(y1.freq, y1.sr, xtitle='Frequency [MHz]', ytitle='(Sig-Ref)/Ref', title='IF Noise  Scan 2  Bank 4  rcvr 2', drop=200);
s.plot(y2.freq, y2.sr, newplot=F, linecolor=3);
s.plot(y3.freq, y3.sr, newplot=F, linecolor=4);
s.plot(y4.freq, y4.sr, newplot=F, linecolor=5);
s.pg.setyscale(-0.0015, 0.0015)
s.plot(f1.x, f1.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 2');
s.plot(f2.x, f2.y, newplot=F, linecolor=3);
s.plot(f3.x, f3.y, newplot=F, linecolor=4);
s.plot(f4.x, f4.y, newplot=F, linecolor=5);
s.pg.setyscale(-0.001, 0.015);
s.pg.setxscale(0.3, 0.55);

- Figure odm8_20grp_spec.ps, odm8_20grp_ft.ps
y := y1;
y.sr := (y1.sr + y2.sr + y3.sr + y4.sr)/4.0;
f := s.ft(y.sr[200:1848], 0.391);
s.plot(y.freq, y.sr, xtitle='Frequency [MHz]', ytitle='(Sig-Ref)/Ref', title='IF Noise  Scan 2  Bank 4  rcvr 2', drop=200);
s.pg.setyscale(-0.0006, 0.0006);
s.plot(f.x, f.y, xtitle='1/f [MHz\\u-1\\d]', ytitle='Power', title='Fourier Transform  IF Noise  Scan 2  Bank 4  rcvr 2');
s.pg.setyscale(-0.001, 0.015);
s.pg.setxscale(0.3, 0.55);


# Simulate sine wave

y50 := s.getSR1(proj, 2, 1, 50, bank=1, rcvr=1);
s.plot(y50.freq, y50.sr)  
f50 := s.ft(y50.sr[200:1848], 0.391);
s.plot(f50.x, f50.y)

# add sine wave with peak-to-peak of 0.01 (rms = 0.01/2.8 = 0.0036)
y50m := y50                                                         
s.plot(y50m.freq, y50m.sr)
for(i in 1:2049){y50m.sr[i] +:= 0.005*sin(as_float(i)*0.391*pi/10.0)} 
s.plot(y50m.freq, y50m.sr)                                          
f50m := s.ft(y50m.sr[200:1848], 0.391); 
s.plot(f50m.x, f50m.y)