function used for recording cbass data for Descartes map making
0001 function d = pipelindData(d) 0002 %function used for recording cbass data for Descartes map making 0003 0004 %h = gov.fnal.eag.healpix.PixTools; 0005 0006 %npix= h.Nside2Npix(nside); 0007 0008 disp('Here') 0009 try 0010 [t,A,G,T,horiz,equa,offStartPos,onEndPos,offEndPos,onStartPos]=calculateAlpha(d); 0011 catch 0012 disp('Failed calculateAlpha') 0013 end 0014 try 0015 d=applyAlpha(d,t,A,G); 0016 catch 0017 disp('Failed the applyAlpha') 0018 end 0019 try 0020 thisAlpha = [t A G T horiz equa]; 0021 d.correction.alpha.indices = [offStartPos' onEndPos' offEndPos' onStartPos']; 0022 d.correction.alpha.values = thisAlpha; 0023 %[r,A]=rFactorCorrection(d,1) 0024 r=calculateRfactor(d); 0025 catch 0026 disp('Failed the rFactorCorrection') 0027 end 0028 try 0029 d=calculateStokes2(d,r,1); 0030 catch 0031 disp('Failed the calculateStokes') 0032 end 0033 try 0034 d = opacity(d); 0035 catch 0036 disp('Failed opacity') 0037 end 0038 % disp('Pixellating the data') 0039 % d=pixdat(d,nside); 0040 % d=aveHealpix(d,nside); 0041 0042 0043 % I = d.antenna0.healpixDat(:,2); 0044 % Q1 = d.antenna0.healpixDat(:,3); 0045 % Q2 = d.antenna0.healpixDat(:,5); 0046 % U1 = d.antenna0.healpixDat(:,4); 0047 % U2 = d.antenna0.healpixDat(:,6); 0048 % N_OBS = d.antenna0.healpixDat(:,14); 0049 % sig_I = d.antenna0.healpixDat(:,8); 0050 % sig_Q1 = d.antenna0.healpixDat(:,9); 0051 % sig_Q2 = d.antenna0.healpixDat(:,11); 0052 % sig_U1 = d.antenna0.healpixDat(:,10); 0053 % sig_U2 = d.antenna0.healpixDat(:,12); 0054 % d.antenna0.writeHealpixDat= [I Q1 Q2 U1 U2 N_OBS sig_I sig_Q1 sig_Q1 sig_U1 sig_U2]; 0055 %unix('rm filename.fits'); 0056 %cbass_write_image('filename.fits',[I Q1 Q2 U1 U2 N_OBS],nside,npix,'start','stop'); 0057 0058 return