Summary of chop-nod analysis telecon of March 10, 2010

John, Giles, Lingzhen

_____________________________________________________________

Lingzhen's work on galaxies:

Decided to do NGC 253 first.

use the file called newlist for the smoothed tau correction and the pointing correction:

novak@zamin:/home/johnv/SHARP/2006dec/process/ngc253$ ls
combine.fits  medbg.fits	 ngc253_5_cmb.fits  ngc253_tmp.ps
docombine     newlist		 ngc253_6_cmb.fits  pointcorr_output.txt
linbgnd.pro   ngc253_1_cmb.fits  ngc253_7_cmb.fits  README
list1	      ngc253_2_cmb.fits  ngc253_dec06.eps   signal2noise
list.pt       ngc253_3_cmb.fits  ngc253_dec06.pdf
list.tau      ngc253_4_cmb.fits  ngc253_dec06.ps

use the rgm file called rgm_new.dat:

novak@zamin:/home/johnv/SHARP/2006dec/calibs/sharp_rgm$ ls
34401.fits.bz2	doit2				     SHARC2PIXEL.34401
34401.txt	doit3				     SHARC2PIXEL.34402
34402.fits.bz2	listall				     SHARC2PIXEL.34403
34402.txt	list.cmp			     SHARC2PIXEL.34404
34403.fits.bz2	rgm_350_34401_2006dec05_ORIIRC2.dat  SHARC2PIXEL.dat
34403.txt	rgm_350_34402_2006dec05_ORIIRC2.dat  solve34401.fits.bz2
34404.fits.bz2	rgm_350_34403_2006dec05_ORIIRC2.dat  solve34402.fits.bz2
34404.txt	rgm_350_34404_2006dec05_ORIIRC2.dat  solve34403.fits.bz2
doit0		rgm_new.dat			     solve34404.fits.bz2
doit1		rgm_old.dat			     temp.list

For i.p. correction, see Giles' January 2008 memo on the analysis logbook

NOTE: Also see Darren's latest e-mail sent just after the telecon (to sharp-software list) for more 
info on smoothed tau.

_____________________________________________________________

Giles' work on IC 348, L 1527, and B 335:

Giles asked John some questions about polsharp5 :

Q: How to make the table?  

A: see polsharp.out (in "IDL bin" on Giles' mac):

   x      y      P       dP    Phi   dPhi   P/dP
  3.6   21.2    9.1     3.4   75.0    8.8   2.72
 23.3   36.5    6.6     2.4  -81.5    9.2   2.72
 23.5   24.3    1.2     0.4   63.4   10.5   2.72
 21.6   32.3    4.8     1.6   82.6    8.8   3.02
 28.8   15.7    7.6     2.8  -11.8   10.1   2.67
 25.8   32.9    4.7     1.9   66.9   11.1   2.52
 24.4   37.0    7.5     2.7   74.2    9.4   2.78
 32.0   14.9   13.8     5.5   -0.0    9.5   2.50
 32.6   18.2   12.4     3.8   -5.3    7.5   3.29
 32.4   18.9    5.8     2.3   -8.2   10.8   2.53
 29.0   34.2    7.5     2.2   53.8    7.5   3.42
 33.4   27.4    8.6     2.2   36.7    6.4   3.92  **
 34.5   31.0    6.6     2.5   26.3   10.2   2.65
 38.1   28.1    9.1     3.3   25.1    9.6   2.75
 36.9   36.5   18.1     5.4   19.9    6.8   3.34

NOTE: The first column has a different zero and different units than one might want for a table in a 
paper, but related linearly.

Q: How to use a different intensity map?

A: Look at the .pro file itself for clues:

;         imagefile = 'filename' -- use this FITS image file for flux map
;                      rather than that in fitsfile
;         hzfile = 'filename' --- Hertz polarization file (ASCII bigfit format)
;         sharp2 = 'filename' -- another SHARP FITS file output
;                   from sharp_combine with vectors to be overplotted.
;                   N.B.: color=n will have to be set in order to distinguish
;                         the vectors in the 2 different SHARP files

Note: if you use a sharp_combine file it will use the *first* fits image in the file

Q: how to do bias correction and efficiency correction?

A: Again, see the .pro file:

;          eff = [e1, e2] --- SHARP polarization efficiencies (fractional) for
;                primary input file (e1) and secondary input file
;                (e2; using sharp2=filename option)
;          /debias --- apply typical debias formula in eqs. (11)-(12) of
;                      Vaillancourt (2006, PASP, 118, 1340) to all data 

_____________________________________________________________