notebook.community

Edit and run 





In [1]:



    


import numpy as np
from astropy.table import Table
from astropy.io.ascii import InconsistentTableError
import matplotlib.pyplot as plt
from scipy import stats
import os
from glob import glob
from scipy.optimize import least_squares



    











In [15]:



    


def model2(theta, mag, y, zp=22.5, sigdet=5, k=1):
    a, b, maglim = theta 
    teff = np.exp(a + b * (maglim - 21.)) 

    # Compute flux/limit. 
    F = 10**(-0.4 * (mag - zp)) 
    Flim = 10**(-0.4 * (maglim - zp)) 

    # Compute noise. 
    Fnoise = (Flim / sigdet)**2 * k * teff - Flim
    magerr = 2.5 / np.log(10.) * np.sqrt((1. + Fnoise / F) / (F * k * teff))

    return np.sum(np.abs(y - magerr))
    #return magerr - y

def model3(theta, mag, y, a=4.56, b=1, zp=22.5, sigdet=5, k=1):
    maglim = theta 
    teff = np.exp(a + b * (maglim - 21.)) 

    # Compute flux/limit. 
    F = 10**(-0.4 * (mag - zp)) 
    Flim = 10**(-0.4 * (maglim - zp)) 

    # Compute noise. 
    Fnoise = (Flim / sigdet)**2 * k * teff - Flim
    magerr = 2.5 / np.log(10.) * np.sqrt((1. + Fnoise / F) / (F * k * teff))

    return magerr - y

def model4(theta, mag, y, zp=22.5, sigdet=5, k=1):
    teff, maglim = theta 
    #teff = np.exp(a + b * (maglim - 21.)) 

    # Compute flux/limit. 
    F = 10**(-0.4 * (mag - zp)) 
    Flim = 10**(-0.4 * (maglim - zp)) 

    # Compute noise. 
    Fnoise = (Flim / sigdet)**2 * k * teff - Flim
    magerr = 2.5 / np.log(10.) * np.sqrt((1. + Fnoise / F) / (F * k * teff))

    #return -np.sum(np.abs(y - magerr))
    return magerr - y



    











In [3]:



    


datapath = f'{os.environ["HOME"]}/Projects/planckClusters/data/extern/SDSS'
files = glob(f'{datapath}/PSZ**/*.csv')



    











In [13]:



    


files



    


















    
Out[13]:








['/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G009.04+31.09/PSZ2_G009.04+31.09_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G023.05+20.52/PSZ2_G023.05+20.52_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G029.87-17.81/PSZ2_G029.87-17.81_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G030.70+09.44/PSZ2_G030.70+09.44_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G031.41-19.16/PSZ2_G031.41-19.16_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G032.31+66.07/PSZ2_G032.31+66.07_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G036.36+16.01/PSZ2_G036.36+16.01_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G044.21+52.13/PSZ2_G044.21+52.13_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G045.96-26.94/PSZ2_G045.96-26.94_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G048.47+34.86/PSZ2_G048.47+34.86_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G051.48-30.87/PSZ2_G051.48-30.87_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G052.08+46.13/PSZ2_G052.08+46.13_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G053.80+36.49/PSZ2_G053.80+36.49_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G058.31+41.96/PSZ2_G058.31+41.96_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G058.42-33.50/PSZ2_G058.42-33.50_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G066.59-58.51/PSZ2_G066.59-58.51_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G067.21-20.49/PSZ2_G067.21-20.49_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G069.47-29.06/PSZ2_G069.47-29.06_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G071.67-42.76/PSZ2_G071.67-42.76_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G076.81-32.57/PSZ2_G076.81-32.57_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G077.67+30.59/PSZ2_G077.67+30.59_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G084.69-58.60/PSZ2_G084.69-58.60_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G086.07-41.99/PSZ2_G086.07-41.99_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G091.27-38.62/PSZ2_G091.27-38.62_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G092.11-33.73/PSZ2_G092.11-33.73_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G092.46-35.22/PSZ2_G092.46-35.22_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G092.69+59.92/PSZ2_G092.69+59.92_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G093.04-32.38/PSZ2_G093.04-32.38_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G093.41-16.26/PSZ2_G093.41-16.26_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G093.71-30.90/PSZ2_G093.71-30.90_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G093.94+13.75/PSZ2_G093.94+13.75_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G094.31-11.31/PSZ2_G094.31-11.31_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G095.00-37.14/PSZ2_G095.00-37.14_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G098.38+77.22/PSZ2_G098.38+77.22_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G098.39+57.68/PSZ2_G098.39+57.68_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G098.62+51.76/PSZ2_G098.62+51.76_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G100.22+33.81/PSZ2_G100.22+33.81_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G104.15-38.85/PSZ2_G104.15-38.85_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G107.83-45.45/PSZ2_G107.83-45.45_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G110.69-46.25/PSZ2_G110.69-46.25_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G112.07-39.86/PSZ2_G112.07-39.86_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G112.54+59.53/PSZ2_G112.54+59.53_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G112.69+33.37/PSZ2_G112.69+33.37_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G115.34-54.89/PSZ2_G115.34-54.89_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G117.38-52.47/PSZ2_G117.38-52.47_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G121.87-45.97/PSZ2_G121.87-45.97_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G125.84-18.72/PSZ2_G125.84-18.72_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G126.36-19.11/PSZ2_G126.36-19.11_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G126.62-53.42/PSZ2_G126.62-53.42_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G126.72-21.03/PSZ2_G126.72-21.03_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G129.99-22.42/PSZ2_G129.99-22.42_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G130.64+37.16/PSZ2_G130.64+37.16_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G131.15-14.72/PSZ2_G131.15-14.72_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G030.98+22.43/PSZ1_G030.98+22.43_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G031.91+67.94/PSZ1_G031.91+67.94_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G040.33-16.55/PSZ1_G040.33-16.55_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G131.19+14.48/PSZ2_G131.19+14.48_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G131.27-25.82/PSZ2_G131.27-25.82_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G134.26-44.28/PSZ2_G134.26-44.28_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G135.94-68.22/PSZ2_G135.94-68.22_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G136.02-47.15/PSZ2_G136.02-47.15_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G136.31+54.67/PSZ2_G136.31+54.67_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G136.33-44.53/PSZ2_G136.33-44.53_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G137.24+53.93/PSZ2_G137.24+53.93_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G137.58+53.88/PSZ2_G137.58+53.88_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G144.84-35.16/PSZ2_G144.84-35.16_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G145.25+50.84/PSZ2_G145.25+50.84_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G146.10-55.55/PSZ2_G146.10-55.55_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G146.13+40.97/PSZ2_G146.13+40.97_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G146.16-48.92/PSZ2_G146.16-48.92_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G147.17+42.67/PSZ2_G147.17+42.67_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G148.60-48.61/PSZ2_G148.60-48.61_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G149.73+24.49/PSZ2_G149.73+24.49_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G150.64-14.21/PSZ2_G150.64-14.21_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G152.40+75.00/PSZ2_G152.40+75.00_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G152.47+42.11/PSZ2_G152.47+42.11_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G153.56+36.82/PSZ2_G153.56+36.82_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G153.68+36.96/PSZ2_G153.68+36.96_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G157.07-33.63/PSZ2_G157.07-33.63_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G160.94+44.85/PSZ2_G160.94+44.85_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G165.39+09.22/PSZ2_G165.39+09.22_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G167.43-53.67/PSZ2_G167.43-53.67_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G171.48+16.17/PSZ2_G171.48+16.17_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G176.07-26.95/PSZ2_G176.07-26.95_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G177.03+32.64/PSZ2_G177.03+32.64_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G179.33-22.22/PSZ2_G179.33-22.22_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G179.45-43.92/PSZ2_G179.45-43.92_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G185.72-32.23/PSZ2_G185.72-32.23_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G186.50-13.45/PSZ2_G186.50-13.45_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G189.79-37.25/PSZ2_G189.79-37.25_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G191.57+58.88/PSZ2_G191.57+58.88_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G192.40-67.89/PSZ2_G192.40-67.89_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G202.61-26.26/PSZ2_G202.61-26.26_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G202.66+66.98/PSZ2_G202.66+66.98_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G203.32+08.91/PSZ2_G203.32+08.91_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G203.71+50.82/PSZ2_G203.71+50.82_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G210.71+63.08/PSZ2_G210.71+63.08_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G218.58+08.71/PSZ2_G218.58+08.71_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G231.41+77.48/PSZ2_G231.41+77.48_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G235.96+38.21/PSZ2_G235.96+38.21_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G237.68+57.83/PSZ2_G237.68+57.83_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G252.45+73.44/PSZ2_G252.45+73.44_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G254.52+62.52/PSZ2_G254.52+62.52_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G269.02+46.30/PSZ2_G269.02+46.30_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G270.88+37.23/PSZ2_G270.88+37.23_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G271.53+36.41/PSZ2_G271.53+36.41_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G310.81+83.91/PSZ2_G310.81+83.91_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G318.46+83.79/PSZ2_G318.46+83.79_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G320.94+83.69/PSZ2_G320.94+83.69_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G321.94+75.57/PSZ2_G321.94+75.57_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G328.96+71.97/PSZ2_G328.96+71.97_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G347.96+80.46/PSZ2_G347.96+80.46_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G041.70+21.65/PSZ1_G041.70+21.65_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G001.00+25.71/PSZ1_G001.00+25.71_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G015.42+58.42/PSZ1_G015.42+58.42_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G018.09+26.07/PSZ1_G018.09+26.07_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G031.41+28.75/PSZ1_G031.41+28.75_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G032.76+42.33/PSZ1_G032.76+42.33_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G044.82-31.66/PSZ1_G044.82-31.66_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G046.98+66.62/PSZ1_G046.98+66.62_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G047.44+37.39/PSZ1_G047.44+37.39_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G050.01-16.88/PSZ1_G050.01-16.88_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G051.42-26.16/PSZ1_G051.42-26.16_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G053.50+09.56/PSZ1_G053.50+09.56_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G054.59-18.18/PSZ1_G054.59-18.18_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G055.83-41.64/PSZ1_G055.83-41.64_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G058.77-26.14/PSZ1_G058.77-26.14_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G063.92-16.75/PSZ1_G063.92-16.75_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G064.83+35.79/PSZ1_G064.83+35.79_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G065.13+57.53/PSZ1_G065.13+57.53_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G066.01-23.30/PSZ1_G066.01-23.30_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G078.39+46.13/PSZ1_G078.39+46.13_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G081.56+31.03/PSZ1_G081.56+31.03_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G083.35+76.41/PSZ1_G083.35+76.41_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G084.04+58.75/PSZ1_G084.04+58.75_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G087.25-41.86/PSZ1_G087.25-41.86_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G090.14-49.71/PSZ1_G090.14-49.71_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G090.48+46.89/PSZ1_G090.48+46.89_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G091.73-30.23/PSZ1_G091.73-30.23_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G091.93+35.48/PSZ1_G091.93+35.48_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G092.41-37.39/PSZ1_G092.41-37.39_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G094.95-36.72/PSZ1_G094.95-36.72_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G100.46-61.45/PSZ1_G100.46-61.45_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G102.86-31.07/PSZ1_G102.86-31.07_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G103.50+31.36/PSZ1_G103.50+31.36_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G103.56-39.35/PSZ1_G103.56-39.35_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G105.91-38.39/PSZ1_G105.91-38.39_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G106.81-36.44/PSZ1_G106.81-36.44_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G108.26+48.66/PSZ1_G108.26+48.66_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G108.90-52.04/PSZ1_G108.90-52.04_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G109.09-52.45/PSZ1_G109.09-52.45_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G112.60-39.30/PSZ1_G112.60-39.30_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G115.59-44.47/PSZ1_G115.59-44.47_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G125.54-56.25/PSZ1_G125.54-56.25_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G126.44-70.36/PSZ1_G126.44-70.36_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G128.75-17.97/PSZ1_G128.75-17.97_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G133.50-46.77/PSZ1_G133.50-46.77_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G134.08-44.61/PSZ1_G134.08-44.61_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G134.75-56.53/PSZ1_G134.75-56.53_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G135.92+76.21/PSZ1_G135.92+76.21_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G138.60-10.85/PSZ1_G138.60-10.85_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G140.10+50.09/PSZ1_G140.10+50.09_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G144.99+54.39/PSZ1_G144.99+54.39_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G146.00-49.42/PSZ1_G146.00-49.42_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G148.20+23.49/PSZ1_G148.20+23.49_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G150.33-20.04/PSZ1_G150.33-20.04_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G151.44-64.86/PSZ1_G151.44-64.86_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G153.41+36.58/PSZ1_G153.41+36.58_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G153.87+41.05/PSZ1_G153.87+41.05_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G158.58-52.47/PSZ1_G158.58-52.47_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G159.29-24.66/PSZ1_G159.29-24.66_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G165.76+31.15/PSZ1_G165.76+31.15_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G165.94+50.48/PSZ1_G165.94+50.48_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G171.01+15.93/PSZ1_G171.01+15.93_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G178.10+18.58/PSZ1_G178.10+18.58_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G182.49-57.09/PSZ1_G182.49-57.09_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G185.93-31.21/PSZ1_G185.93-31.21_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G186.54+85.60/PSZ1_G186.54+85.60_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G197.20-59.04/PSZ1_G197.20-59.04_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G198.06-12.52/PSZ1_G198.06-12.52_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G199.70+37.01/PSZ1_G199.70+37.01_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G203.42-04.03/PSZ1_G203.42-04.03_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G206.52-26.37/PSZ1_G206.52-26.37_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G206.64-21.17/PSZ1_G206.64-21.17_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G212.80+46.65/PSZ1_G212.80+46.65_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G213.27+78.35/PSZ1_G213.27+78.35_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G218.54+13.26/PSZ1_G218.54+13.26_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G223.80+58.50/PSZ1_G223.80+58.50_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G236.51+48.39/PSZ1_G236.51+48.39_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G240.42+77.58/PSZ1_G240.42+77.58_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G249.01+73.75/PSZ1_G249.01+73.75_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G264.01+59.77/PSZ1_G264.01+59.77_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G268.86+55.92/PSZ1_G268.86+55.92_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G286.25+62.68/PSZ1_G286.25+62.68_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G332.30+72.17/PSZ1_G332.30+72.17_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G346.57+42.71/PSZ1_G346.57+42.71_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G355.14+55.96/PSZ1_G355.14+55.96_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G022.03+17.75/PSZ2_G022.03+17.75_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G029.80-17.40/PSZ2_G029.80-17.40_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G032.77+19.69/PSZ2_G032.77+19.69_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G039.86+18.70/PSZ2_G039.86+18.70_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G046.39+11.71/PSZ2_G046.39+11.71_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G048.39-16.78/PSZ2_G048.39-16.78_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G067.72-24.13/PSZ2_G067.72-24.13_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G092.34+14.22/PSZ2_G092.34+14.22_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G092.64+20.78/PSZ2_G092.64+20.78_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G095.75-11.16/PSZ2_G095.75-11.16_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G185.45-32.01/PSZ2_G185.45-32.01_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G185.50-33.94/PSZ2_G185.50-33.94_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G185.68+09.82/PSZ2_G185.68+09.82_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G216.25+10.10/PSZ2_G216.25+10.10_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G227.30+09.00/PSZ2_G227.30+09.00_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ2_G349.18+38.66/PSZ2_G349.18+38.66_SDSS_catalog.csv',
 '/home/boada/Projects/planckClusters/data/extern/SDSS/PSZ1_G153.56+36.32/PSZ1_G153.56+36.32_SDSS_catalog.csv']

















In [19]:



    


band = 'i'
maglim = []
for f in files:
    try:
        cat = Table.read(f)
    except InconsistentTableError:
        cat = Table.read(f, header_start=1, data_start=2)
    # filter out the stars
    cat = cat[cat['type'] == 3]
    if len(cat) < 100:
        continue
    result = least_squares(model2, (4, 1, 20), loss='cauchy', f_scale=np.std(cat[f'fiberMagErr_{band}']),
                       args=(cat[f'fiberMag_{band}'], cat[f'fiberMagErr_{band}']), 
                       bounds=([3,0.5, 18], [5,2, 25]), kwargs={'zp':22.5})
    maglim.append(result['x'][2])



    











In [20]:



    


np.median(maglim)



    


















    
Out[20]:








22.13957126364067

















In [27]:



    


plt.hist(np.array(maglim) - 22.2, bins=np.linspace(-1, 1, 50))
#plt.xlim(-0.5, 0.5)
plt.xlabel('maglim - sdsslim')
ax = plt.gca()
plt.text(0.8, 0.8, f'{band}', transform=ax.transAxes, fontsize=30)
plt.grid()
plt.savefig('zlim.png', bbox='tight')



    


















    
























In [22]:



    


np.std(maglim)



    


















    
Out[22]:








0.13765118825514713

















In [23]:



    


len(files)



    


















    
Out[23]:








214

















In [24]:



    


len(maglim)



    


















    
Out[24]:








212

















In [18]:



    


len(cat)



    


















    
Out[18]:








1

















In [25]:



    


maglim



    


















    
Out[25]:








[22.13989924068213,
 22.03319256312883,
 22.105963250491456,
 22.25526379254132,
 22.088863590996187,
 22.201904303744683,
 22.243263159895474,
 22.201931345167402,
 22.017902011343466,
 22.09824226183278,
 22.1538954447528,
 22.390817524841477,
 22.288800414309957,
 22.15410523200863,
 22.078550243566035,
 22.078202189505227,
 22.06355278151773,
 22.05328919686633,
 22.16955901747614,
 22.008635940579747,
 22.138669357323028,
 21.99644098067351,
 21.996540256400653,
 22.001573729756014,
 22.140893479607335,
 22.18746436039054,
 22.118382241927154,
 22.169503644771538,
 22.230933410081278,
 22.272224931706386,
 22.302156111991348,
 22.192801919122846,
 22.217116539931208,
 22.34816144341298,
 22.143976718407682,
 22.314935047140203,
 22.11094259918202,
 21.995646135901847,
 21.98593710907667,
 22.03198498924528,
 22.05334865221346,
 22.1901625320884,
 21.947836677472395,
 22.14178617815363,
 22.133994396238993,
 22.02238845993229,
 22.2382110909303,
 22.286932148917714,
 22.005368573000148,
 22.280413685761108,
 22.194307382788622,
 22.029979919258697,
 22.29995001507309,
 22.139243286599203,
 22.20862789886526,
 22.11569736290571,
 22.087460835890514,
 22.31359010725374,
 22.034850948281527,
 22.099163549305306,
 22.104468005807934,
 22.28744505220428,
 22.002309864179836,
 22.120615753307096,
 22.126573662387862,
 22.07743556218359,
 22.123830354683818,
 22.147497574274084,
 22.194692052515613,
 21.97832660717114,
 21.91611634480617,
 21.959000788990448,
 22.16987668302287,
 22.264119396353543,
 22.23800050377299,
 22.01431962225958,
 22.199845002987335,
 22.2544042703315,
 22.045406811522486,
 22.324701577214416,
 21.771146262307774,
 22.13354545367421,
 21.73693757369795,
 22.105545440961247,
 22.195660104486517,
 22.286997591668978,
 22.17203630128716,
 22.04011602847732,
 22.005607768230426,
 22.253524038194257,
 22.22045978085639,
 22.107405558274184,
 22.31183338061261,
 22.186375312789195,
 22.08646760341893,
 22.107089092265454,
 22.009706891328843,
 21.90408430975012,
 22.25842899899832,
 22.13222150635493,
 22.12996574253558,
 22.13783998377934,
 22.23948894629359,
 22.219804716370106,
 22.157858713597506,
 22.15730583298902,
 22.275494304515515,
 22.272811776400914,
 22.27229885792802,
 22.200162288181158,
 22.165245701140122,
 22.114443626298407,
 22.338315735855616,
 21.854269495336755,
 22.06493851726256,
 22.085561479204184,
 22.17597933458561,
 22.204047497606066,
 22.197521904298817,
 22.167765923389677,
 22.197367524711765,
 22.258205063541347,
 22.22074254636119,
 22.20375585745528,
 21.95317385188259,
 22.067614064172698,
 22.080018163219627,
 21.893470444049363,
 22.255722873981256,
 22.244486969279418,
 22.216023373376544,
 22.24606310122821,
 22.141393735623485,
 22.26151652338132,
 22.21807099854541,
 22.01473622956524,
 22.183720944517265,
 22.093664701668388,
 21.828095175720357,
 22.131234809324738,
 22.162758729301967,
 22.222666298032,
 22.142689854931053,
 22.197431578739927,
 22.033160090578143,
 21.98204878660101,
 22.00313591405844,
 22.307601407028567,
 22.024678712490438,
 22.08172654529274,
 22.071291197915592,
 22.04330655022171,
 22.069464377330693,
 22.192118376463377,
 22.00237429318566,
 22.25237073084867,
 22.092368843922937,
 22.012405607212816,
 22.133391120495425,
 22.230948613160113,
 22.157961421092306,
 22.21725645482593,
 22.327312494231176,
 21.980263221040914,
 22.110211690373703,
 22.11546708726473,
 22.156906563392187,
 22.193576815650047,
 21.969232610302758,
 22.21198148829942,
 21.885472640730132,
 22.254930732263475,
 22.150517672117545,
 21.731733227169805,
 22.321817039447858,
 22.031601867564806,
 22.018068130515164,
 22.11611001485142,
 22.156018844310015,
 21.745066091717568,
 22.247554050923128,
 22.33525690026594,
 21.72859359684438,
 21.729185157397918,
 22.252786484846578,
 22.073195020517435,
 21.80718903384758,
 22.116716251661312,
 22.340988090018804,
 22.272700214297743,
 22.050697437726054,
 21.92201063107949,
 22.338033357621,
 22.222347437031978,
 22.121744527474586,
 22.119810178886933,
 22.28384295499266,
 21.52172397310003,
 22.12107049784167,
 22.25470556699994,
 22.121158053928017,
 22.157606716759233,
 22.208210655532756,
 22.31158227901672,
 22.365383643903925,
 22.049457604306713,
 21.905731698556313,
 22.014029609496795,
 21.868337014511564,
 22.255468490507216,
 22.0833038934787,
 22.1990251841942]

















In [ ]:

Content source: boada/planckClusters

Similar notebooks:

notebook.community | gallery | about