Home > example_scripts > pipelinedDatav10.m

pipelinedDatav10

PURPOSE ^

function used for recording cbass data for Descartes map making

SYNOPSIS ^

function d = pipelindData(d)

DESCRIPTION ^

function used for recording cbass data for Descartes map making

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

Generated on Sun 14-Jun-2015 17:12:45 by m2html © 2005