Aug 2, 2012

Palomar 5 Project Summary


Files for this project can be found in /home/ch558/palomar5/

Membership Selection:
=========================
KECK
cmd: dist le 0.1 to isochrone
(isochrone dist modulus = 16.808, see selection.pro for extinction, etc)
velocity cuts: -70 < * < -40 km/s
HYDRA
cmd: dist le 0.1 to isochrone
(isochone dist modulus = 16.808, see hydraselection.pro for extinction, etc)
velocity cuts: -70 < * < -35 km/s
magnesium: on a plot of (mg1 + mg2) vs (mgb), members are < 2 from (0.15,1.0)
ODEN
velo: -80 < * < -30
(done in allvelo.pro and allcenter.pro) 

Overlapping data:
=========================
overlap.pro gives overlaps between keck and hydra data
doubles.pro finds doubles in the hydra data
oden.pro looks for overlap between odenchirken's data and keck+hydra


Relevant code:
=========================
allvelo.pro
Reads in keck, hydra, and oden data and combines it (requires that the alldata_Pal5.fits.gz and allhydra_Pal5.fits.gz files have usable x.member columns (that = 1 for members). Plots l*cos(b) versus velocity, creates bins. Plots the bins and the dispersion bins, and then does a linear fit to the velocity bins. Prints slope of fit line (the velocity gradient per degree) and the average dispersion.


allcenter.pro
Same as allvelo.pro, but only looks at the center of pal5 (limits made and can be changed in the binning section of the code... change max_covered's first definition to control start point of binning and then change n_bins and the binning increments to adjust end point)


overlap.pro
Finds overlap between hydra and keck


doubles.pro
Finds doubles in the hydra data


oden.pro
Looks for overlap between odenkirchen's data and keck+hydra


selection.pro
Does selection for keck data, including lots of plots showing the cuts. Requires the besancon model (mw-model.resu) and m92.iso isochrone)


hydraselection.pro
Same as selection.pro, but for hydra data. 


Note: selection.pro is for keck data, hydraselection.pro for hydra data. Both have commented-out sections at the end which write the "member" column back to the fits file


Data files:
=========================
alldata_Pal5.fits.gz
allhydra_Pal5.fits.gz
(both files have member=1 for members selected by [hydra]selection.pro)

oden_center.dat
oden_tail.dat
(membership cuts only done in allvelo.pro, allcenter.pro, "just in time")


hydradoubles.dat
(two-column file listing indicies for pairs in hydra data after a quality cut)


Misc Notes
=========================
Hydra systemic error calculated to be +/- 9.8, this has already been 
added into the fits file errors


What's Next
=========================
• The doubled stars (doubles.pro, overlap.pro, oden.pro) should be folded
into each other, rather than being counted multiple times.


• In allvelo.pro, try fitting the velocity gradient line to all of the
member stars, instead of binning them (likely will just give far too much 
weight to the center of Pal5, but easily tried)


• Membership probabilites instead of hard binary cuts. Right now the members are likely, but not secure.

No comments: