def mk_lightcurve_fig( ): from hstsnpipe.tools import snana from hstsnpipe.tools.figs import plotsetup plotsetup.fullpaperfig( figsize=[8,4] ) pl.clf() simIa = snana.SimTable( 'snIa_zgrid' ) FLTMATRIX = simIa.FLT.reshape( simIa.LCMATRIX.shape ) bandlist = FLTMATRIX[0,0,0,0,:,0] def iband(band) : return( np.where( bandlist==band )[0][0] ) iz18 = np.argmin( np.abs( simIa.z-1.8 ) ) ilp0 = np.argmin( np.abs( simIa.x1 ) ) iF160W = iband('H') iF140W = iband('N') iF110W = iband('M') muIa = 5. ax1 = pl.subplot(121) ax2 = pl.subplot(122, sharex=ax1 ) # snIaMAGH = simIa.LCMATRIX[ ilp0, 0, 0, iz18, iF160W, : ] snIaMAGN = vega2ab( simIa.LCMATRIX[ ilp0, 0, 0, iz18, iF140W, : ], 'F140W' ) snIaMAGM = vega2ab( simIa.LCMATRIX[ ilp0, 0, 0, iz18, iF110W, : ], 'F140W' ) snIaMJD = simIa.TOBS[ iz18, : ] # ax1.plot( snIaMJD, snIaMAGH - 2.5*np.log10( muIa ), 'k-', lw=2, marker=' ', ) ax1.plot( snIaMJD, snIaMAGN - 2.5*np.log10( muIa ), 'r-', lw=2, marker=' ', ) ax2.plot( snIaMJD, snIaMAGM - 2.5*np.log10( muIa ), 'r-', lw=2, marker=' ', ) ax1.invert_yaxis() simII = snana.SimTable( 'snII_zgrid' ) FLTMATRIX = simII.FLT.reshape( simII.LCMATRIX.shape ) bandlist = FLTMATRIX[0,0,0,0,:,0] def iband(band) : return( np.where( bandlist==band )[0][0] ) iz20 = np.argmin( np.abs( simII.z-2.0 ) ) iF140W = iband('N') iF110W = iband('M') muII = 12. ilpII=13 # select one representative II-P model snIIMAGN = vega2ab( simII.LCMATRIX[ ilpII, 0, 0, iz20, iF140W, : ], 'F140W' ) snIIMAGM = vega2ab( simII.LCMATRIX[ ilpII, 0, 0, iz20, iF110W, : ], 'F110W' ) snIIMJD = simII.TOBS[ iz20, : ] ax1.plot( snIIMJD, snIIMAGN - 2.5*np.log10( muII ), 'b--', lw=2, marker=' ', ) ax2.plot( snIIMJD, snIIMAGM - 2.5*np.log10( muII ), 'b--', lw=2, marker=' ', ) ax1.axhline( vega2ab( 24.6, 'f140w'), color='k',ls=':', lw=2, ) ax2.axhline( vega2ab( 25.3, 'f110w'), color='k',ls=':', lw=2, ) # ax2.text(0.95,0.95,'F160W', fontsize='large', ha='right',va='top', transform=ax2.transAxes ) ax1.text(0.95,0.95,'F140W', fontsize='large', ha='right',va='top', transform=ax1.transAxes ) ax2.text(0.95,0.95,'F110W', fontsize='large', ha='right',va='top', transform=ax2.transAxes ) ax1.text( 115, 25.65, '30 min SNAP\n detection limit', ha='right',va='bottom') # ax2.text( 115, 26.0, '30 min SNAP\n detection limit', ha='right',va='bottom') ax1.text( 38, 25.3, 'Type Ia SN\n z=1.8, $\mu$=5\n t$_{vis}$=75 days', ha='left',va='bottom', color='r', fontsize='large') # ax2.text( 25, 25.5, 'Type II-P SN\n z=2.0, $\mu$=12\n t$_{vis}$=50 days', ha='left',va='center', color='b', fontsize='large') # ax1.text( 67, 25.6, 'Type II-P SN\n z=2.0, $\mu$=12\n t$_{vis}$=50 days', ha='left',va='top', color='b', fontsize='large', backgroundcolor='w' ) ax2.text( 67.5, 26.15, 'Type II-P SN\n z=2.0, $\mu$=12\n t$_{vis}$=50 days', ha='left',va='top', color='b', fontsize='large', backgroundcolor='w' ) ax1.set_xlim(-35,119) ax1.set_ylim(27, 27-2.51) ax2.set_ylim(27.4, 27.4-2.51) ax2.yaxis.set_ticks_position('right') ax2.yaxis.set_ticks_position('both') ax2.yaxis.set_label_position('right') fig = pl.gcf() fig.subplots_adjust( left=0.13, right=0.87, bottom=0.12, top=0.95, wspace=0.05 ) ax2.set_xlabel('Observer-frame Time [days]') ax1.set_xlabel('Observer-frame Time [days]') ax1.set_ylabel('Observed F140W magnitue [AB]') ax2.set_ylabel('Observed F110W magnitue [AB]', labelpad=20, rotation=-90)
import djs_angle_match import numpy as np import matplotlib.pyplot as plt import matplotlib import plotsetup from matplotlib import gridspec plotsetup.fullpaperfig() co=0 xlabel=['g mag','MJD/100','Airmass','Ang. Dist. (Deg)'] fil2=['g','r','i','z'] fields=['053','054','055','057','058','059'] fields=['0','053','054','056','057','058','059'] fbin=[0,0.1,0.2,0.3,0.5,0.7,1.0,1.3,2.0] if (1>0): co=0 l1=0 jj=0 gs1 = gridspec.GridSpec(1, 4) gs1.update(bottom=0.17, top=0.95, wspace=0.0) ax1= plt.subplot(gs1[0]) ax2= plt.subplot(gs1[1]) ax3= plt.subplot(gs1[2]) ax4= plt.subplot(gs1[3]) ax=[ax1,ax2,ax3,ax4] for k1 in range(0,2): for k2 in range(0,2): if (co==0):fil,mag,color,off1,off1e,num = np.loadtxt('/project/rkessler/dscolnic/DougComp/ps1star/ncomp_'+fil2[l1]+'_'+fields[jj]+'_1_mag.txt', usecols=(0,1,2,3,4,5), unpack=True, dtype='string')
def mkTvisHistFigs(): import plotsetup plotsetup.fullpaperfig(1, [8,3]) datfile1 = 'lensed_galaxies.txt' dat1 = ascii.read(datfile1, header_start=-1) mu1 = dat1['mu1'] z1 = dat1['z'] dt1 = dat1['dt1'] datfile2 = 'lensed_galaxies2.txt' dat2 = ascii.read(datfile2) mu2 = dat2['col6'] z2 = dat2['col5'] dt2 = dat2['col1'] dterr2p = dat2['col2']/11.2 dterr2m = dat2['col3']/11.2 dterr2 = dat2['col4']/11.2 ax1 = pl.subplot(131) z = np.append( z1, z2 ) dt = np.append( dt1, dt2 ) i5 = np.where( (dt>0) & (dt<5) )[0] print( len(i5) / float(len(dt) ) ) print( len(i5) ) randoff = (np.random.randn( len(z) ) -0.5 ) * 0.1 z = z+randoff pl.hist( z , normed=False, bins=np.arange(0.1,4.51,0.2), color='darkcyan')#, alpha=0.3 ) # pl.hist( z[i5] , bins=np.arange(0.1,4.51,0.2), color='darkorange')#, alpha=0.3 ) ax1.set_xlabel('Redshift') ax1.set_ylabel(r'Number of Lensed Galaxy Images') ax2 = pl.subplot(132) pl.hist( dt[i5] , bins=np.arange(0.0,5.01,0.5), color='darkcyan')#, alpha=0.3 ) ax2.set_xlabel(r'Time Delay $\Delta$t') ax3 = pl.subplot(133) ipos = np.where( (dt2>0) ) # err = dterr2[i5] / dt2[i5] err = dterr2[ipos] / dt2[ipos] pl.hist( err ,normed=False, bins=np.arange(0.0,3.0,0.05), color='darkcyan' )#, alpha=0.3 ) ax3.set_xlabel(r'Time Delay \%Error') # pl.hist( z[i5] , bins=np.arange(0.1,4.51,0.2), color='darkorange', alpha=0.3 ) #ax3 = pl.axes( [0.5,0.5,0.4,0.4], transform=ax2.transAxes ) #ax3.hist( dt, bins=np.arange(0.0,0.1,0.02), color='darkorange', alpha=0.3 ) fig = pl.gcf() fig.subplots_adjust(left=0.08, bottom=0.18,right=0.95, top=0.90, wspace=0.05,) ax1.set_xlim(0.3,4.6) ax2.set_xlim([0.01, 4.99]) ax3.set_xlim([0.01, 0.99]) # ax1.set_xticklabels([1,2,3,4]) ax2.set_yticklabels([]) ax3.set_yticklabels([])