Page Comparison

...

• Noting that, by eye, the binned residual data for some of the surveys with a relatively wide redshift range had apparent oscillations in redshift (probably most apparent in PS1MD), I tried fitting a shifted sine function to all of the surveys for which we have >2 SN.  The results are shown below.  Basics of the fit include:
  • The fitting function was of the form: Res(z) = C + A * sin(k * z + \phi) where C, A, k, and \phi are the fitted parameters
  • The domain of the fitting parameters was largely unconstrained
  • The fit was performed on the unbinned data for each survey 
  • The fits for each survey were entirely independent
• In the below plots:
  • The top panel shows the unbinned data and the bottom panel shows the binned data for each survey
    • There were 20 bins of equal redshift space over the redshift range covered by the particular survey in question
  • There is no particular order to the grid set up (it would also not be hard to rearrange the grid, if there is a reason a new layout would be helpful)
  • The fit lines are shown in red
  • The best fit parameters (rounded to 5 decimal points) are shown in the red text underneath each line
    • In order, the parameters in the text are (see equation definition above):
      • C: the constant, overall residual shift
      • A: the amplitude 
      • k: the wavenumber 
      • \phi: the phase shift
• I see no consistent signal in the below plots.  However, I am not yet convinced that there is nothing here.  I think we just need to:
  • (a) be smarter about our approach to fitting, perhaps starting with a better function
  • (b) recognize that most of the surveys cover too small a range of redshifts to give much information about larger variations
  • (c) Try to focus on SN in matching regions of sky to be sure we aren't entangling extinction artifacts with redshift signals (though I do note that we saw now obvious extinction dependence in our previous analysis)
    

  Fitting to a single mode gaussian:

  • Noting that, by eye, many of the fields for the PS SN show an apparent dip in the residual around z~0.3, I tried fitting that data to a single mode Gaussian binned to the SN from each survey in each field individually 
  • Basics of the fit include:
    • The fitting function was of the form: Res(z) = C + A * e ^ (-(z - \mu)^2 / (2 w^2)) where C, A, \mu, and w are the fitted parameters
    • The domain of the fitting parameters was largely unconstrained
    • The fit was performed on the unbinned data for each survey 
    • The fits for each survey were entirely independent
    • All fits were primed with the following values to try to guide them to the object of interest: (C, A, \mu, w) = (-0.05, -0.1, 0.3, 0.05)
  • In the below plots:
    • In each plot, the upper panel shows the unbinned data and the lower panel shows the binned data 
      • binning scheme is 20 bins of equal redshift space, determined for each survey individually 
        • so in a field with multiple surveys, the bins for each survey are different and likely do not cover the whole displayed z-range; only the range covered by the survey in question 
    • The best fit parameters (rounded to 5 decimal points) are shown in the text underneath.  
      • In order, the parameters in the text are (see equation definition above):
        • C: the constant, overall residual shift
        • A: the constant scale
        • \mu the center
        • w the width
    • The colors of the fit curves and text match the colors of the points
  • Some of the fits did pick out the dip I've been looking at 
    • For the PS1MD data, Fields 0, 1, 3(?), 5 and 6 identify the dip 
    • For the SNLS data, Field 0 identified the dip 
    • For the SDSS data, field 9 identified the dip
  • I think some of the other fields could identify the dip as well if we apply some constraints on the range of the fit, or make it a bit more flexible (maybe allow for a polynomial addition instead of a constant shift?)
    • However, we should think about how much a fit could be made to identify a dip anywhere, with sufficient prompting... 

• To compare the PS1MD fields (those above) to the SDSS fields, I generated a plot for just the SDSS SN and fit a single mode Gaussian to them 
  • The fitting information is like that above, with similar constraints
  • The fit did select a dip at a similar position to that noted above (though somewhat narrower).  This data mays also show evidence of an earlier dip (around z ~0.18) that I do not note in the PS1MD fields above (perhaps because of insufficient sampling) 
    Image Added