[go: up one dir, main page]
Menu ▾ ▴

[73c35b]: / init / COR_box.m  Maximize  Restore  History

Download this file

64 lines (52 with data), 2.3 kB 

 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
% COR_box: defines lag profiles for GEN-remote background and calibration 

% GUISDAP v.1.60 96-05-27 Copyright Asko Huuskonen and Markku Lehtinen

%

% input parameters

% ra      : First result memory address to use

% ri      : result memory address increment

% vc      : virtual channel number

% type    : background ('b'), calibration ('c') or signal ('s')

% gating : number of crossed products added together

% overlap: number of products skipped between gates (positive skipping)

%         : if negative, then same crossed products are used for successive gates

% N_gates : number of gates calculated

% lags    : Lags for which output needed in us, e.g. 0:10:250;

% code    : user specified code number

%

% See also: CORR_box, COR_check, COR_caltemp, COR_status, pulse_times

%function COR_lp(ra,ri,vc,type,gating,overlap,N_gates,lags,code)

 function COR_lp(ra,ri,vc,type,gating,overlap,N_gates,lags,code)



global vc_adcint vc_sampling vc_samplingend vc_ba bm_samples bm_next

global lp_t1 lp_t2 lp_h lp_ra lp_nfir lp_fir lp_dec lp_T lp_dt p_dtau

global lp_nt lp_vc lp_ri lp_bcs lp_code lp_ind lp_indold ra_next ra_prev



if (type~='b' & type~='c'),

  error(' Unknown data type in COR_box')

end  



ra_prev=ra;



adcint=vc_adcint(vc); lags=lags/p_dtau; 

COR_check(lags,adcint)



N_lags=length(lags);

index=(lp_ind+1):(lp_ind+N_lags);



lp_T(index)=COR_caltemp(type)*ones(1,N_lags);



% find first when pulses are transmitted and samples are taken

pulsetimes=pulse_times(vc);

t1=0;% The remote timing is still not OK 

%t1=sample_times(vc,type2ind(type));

lp_t1(index)=t1(1)*ones(1,N_lags);

lp_h(index)=lp_t1(index)-pulsetimes(1);

  

% finally parameters common to background, calibration and signal

lp_t2(index)=lp_t1(index)+lags;

lp_dt(index)=adcint*ones(1,N_lags);

lp_nfir(index)=gating*ones(1,N_lags);

lp_fir(1:gating,index)=ones(gating,length(index));

lp_dec(index)=(gating+overlap)*ones(1,N_lags);

lp_nt(index)=N_gates*ones(1,N_lags);

lp_vc(index)=vc*ones(1,N_lags);

lp_ra(index)=ra+(0:N_lags-1);

lp_ri(index)=ri*ones(1,N_lags);

lp_bcs(index)=type*ones(1,N_lags);

lp_code(index)=code*ones(1,N_lags);

lp_ind=lp_ind+N_lags;





% first result memory location not used

ra_next=max(lp_ra(index)+(lp_nt(index)-1).*lp_ri(index))+1;



COR_status('  COR_box',vc,type,index)