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 [ta,A,G]=calculateAlpha(d); 0011 catch 0012 disp('Failed calculateAlpha') 0013 end 0014 try 0015 d=applyAlpha(d,ta,A,G) 0016 catch 0017 disp('Failed the applyAlpha') 0018 end 0019 try 0020 [r,A]=rFactorCorrection(d,1) 0021 catch 0022 disp('Failed the rFactorCorrection') 0023 end 0024 try 0025 d=calculateStokes(d,r); 0026 catch 0027 disp('Failed the calculateStokes') 0028 end 0029 try 0030 d = opacity(d); 0031 catch 0032 disp('Failed opacity') 0033 end 0034 disp('Pixellating the data') 0035 % d=pixdat(d,nside); 0036 % d=aveHealpix(d,nside); 0037 0038 0039 % I = d.antenna0.healpixDat(:,2); 0040 % Q1 = d.antenna0.healpixDat(:,3); 0041 % Q2 = d.antenna0.healpixDat(:,5); 0042 % U1 = d.antenna0.healpixDat(:,4); 0043 % U2 = d.antenna0.healpixDat(:,6); 0044 % N_OBS = d.antenna0.healpixDat(:,14); 0045 % sig_I = d.antenna0.healpixDat(:,8); 0046 % sig_Q1 = d.antenna0.healpixDat(:,9); 0047 % sig_Q2 = d.antenna0.healpixDat(:,11); 0048 % sig_U1 = d.antenna0.healpixDat(:,10); 0049 % sig_U2 = d.antenna0.healpixDat(:,12); 0050 % d.antenna0.writeHealpixDat= [I Q1 Q2 U1 U2 N_OBS sig_I sig_Q1 sig_Q1 sig_U1 sig_U2]; 0051 %unix('rm filename.fits'); 0052 %cbass_write_image('filename.fits',[I Q1 Q2 U1 U2 N_OBS],nside,npix,'start','stop'); 0053 0054 return