# Difference: 20160111MSciPhotometryAnalysisOb1part5 (1 vs. 11)

#### Revision 1114 Jan 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"
Deleted:
<
<

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
• To generate the positions of stars we uniformly generated x and y pixel positions.
• To create a background we generated Poisson distributed random variables around the mean of 160, which corresponded to a typical value for images of Location 1.
• For each position we wanted to create a Gaussian distribution that represents a star centered on the star's x,y position.
• As the standard deviations should be the same for each star we decided to generate only one value for for all stars. We also assumed that the standard deviations are the same in both directions so we generated only one for and .
• To experiment we uniformly generated values for and for each individual star.
• We then for each pixel in the image calculated the Gaussian value at that particular pixel for each star.
• Resulting image for 50 stars is shown below:

### Comparing stars between different filters

• In order to compare the quality of our data, we looked at the data with the largest signal to noise ratio, in this case data taken with a green filter.
• For this, we found the largest and smallest stars, as determined by the amount of pixels which constituted them.
• The coordinates of the stars were recorded, so they can later be compared to the same stars in other filters.
Changed:
<
<
• Below fits can be seen for the maximum and minimum size stars for the green filter:
>
>
• For the fitting the backgrounds were removed, therefore the error on the no background data is the following (assuming error is Poissonian):
• Below fits can be seen for the maximum and minimum size stars for all filters:

Changed:
<
<
>
>
• The plots are projections (they are summed in rows and columns), so to plot the errors the errors were summed in the following way:
• where the is the error on the data point (pixel value).

### Introductory estimation

Changed:
<
<
• Then set about finding for the above stars for no filters and the green filter.
>
>
• Then set about finding for the above stars for all filters.

Changed:
<
<
 Visual Brightest Star Visual Dimmest Star Green Filter Brightest Star Green Filter Dimmest Star Degrees of Freedom 53 53 53 53 3640000 61100 1930000 67700 reduced 68600 1150 36500 1277
>
>
 Visual Brightest Star Visual Dimmest Star Green Filter Brightest Star Green Filter Dimmest Star Red Filter Brightest Star Red Filter Dimmest Star} Blue Filter Brightest Star Blue Filter Dimmest Star

• In general, a good should be around equal to the number of degrees of freedom,

Changed:
<
<
• In the above case, they all seem to give very poor agreement between the fit and the data, as the values are all very large in comparison to the number of degrees of freedom.

>
>
• From the plots it is obvious that the Poisson errors are not correct as they are too small for the brightest star.

• Also the brightest stars are poorly fitted and therefore have high in comparison to dimmest stars.

### Testing curve_fit for better error estimation

• In order to better obtain errors from curve_fit (as it tends to give outlandish errors), we set out to create a Gaussian of known parameters, and then fit a general Gaussian curve to this.
 Fractional Error When Guessing Incorectly Fractional Error When Guessing Exactly
Line: 79 to 81

%ENDLATEX%

>
>

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
• To generate the positions of stars we uniformly generated x and y pixel positions.
• To create a background we generated Poisson distributed random variables around the mean of 160, which corresponded to a typical value for images of Location 1.
• For each position we wanted to create a Gaussian distribution that represents a star centered on the star's x,y position.
• As the standard deviations should be the same for each star we decided to generate only one value for for all stars. We also assumed that the standard deviations are the same in both directions so we generated only one for and .
• To experiment we uniformly generated values for and for each individual star.
• We then for each pixel in the image calculated the Gaussian value at that particular pixel for each star.
• Resulting image for 50 stars is shown below:
-- ElenaCukanovaite - 11 Jan 2016
Line: 115 to 134

 META FILEATTACHMENT attachment="sigma_x_loc_2.png" attr="" comment="" date="1452788220" name="sigma_x_loc_2.png" path="sigma_x_loc_2.png" size="47938" user="zxap014" version="1" attachment="latex8906deced046f2a1528ba6815ff3285e.png" attr="h" comment="" date="1452788248" name="latex8906deced046f2a1528ba6815ff3285e.png" user="zxap014" version="1" attachment="latex970d0dad73601c12ab733877a694729d.png" attr="h" comment="" date="1452788520" name="latex970d0dad73601c12ab733877a694729d.png" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="elena_no_br.png" attr="" comment="" date="1452798724" name="elena_no_br.png" path="elena_no_br.png" size="101223" user="zxap014" version="1" attachment="elena_no_d.png" attr="" comment="" date="1452798729" name="elena_no_d.png" path="elena_no_d.png" size="102546" user="zxap014" version="1" attachment="elena_g_b.png" attr="" comment="" date="1452800430" name="elena_g_b.png" path="elena_g_b.png" size="106913" user="zxap014" version="1" attachment="elena_g_d.png" attr="" comment="" date="1452800436" name="elena_g_d.png" path="elena_g_d.png" size="100413" user="zxap014" version="1" attachment="elena_b_b.png" attr="" comment="" date="1452801085" name="elena_b_b.png" path="elena_b_b.png" size="106557" user="zxap014" version="1" attachment="elena_b_d.png" attr="" comment="" date="1452801089" name="elena_b_d.png" path="elena_b_d.png" size="106535" user="zxap014" version="1" attachment="elena_r_b.png" attr="" comment="" date="1452801476" name="elena_r_b.png" path="elena_r_b.png" size="105999" user="zxap014" version="1" attachment="elena_r_d.png" attr="" comment="" date="1452801480" name="elena_r_d.png" path="elena_r_d.png" size="107593" user="zxap014" version="1" attachment="latex9c5e0ac6869e420fd954e08d0b7f2665.png" attr="h" comment="" date="1452802864" name="latex9c5e0ac6869e420fd954e08d0b7f2665.png" user="zxap014" version="1" attachment="latex863c65acecd243ddd78cdb3dfd912be6.png" attr="h" comment="" date="1452802960" name="latex863c65acecd243ddd78cdb3dfd912be6.png" user="zxap014" version="1" attachment="latex160454a2d1b6d874920cd3f90f7dc56c.png" attr="h" comment="" date="1452802960" name="latex160454a2d1b6d874920cd3f90f7dc56c.png" user="zxap014" version="1" attachment="latex0a3a7b13ca4e35a8f0232242f6646a97.png" attr="h" comment="" date="1452803639" name="latex0a3a7b13ca4e35a8f0232242f6646a97.png" user="zxap014" version="1" attachment="latex77039efc6bd25fa1afbcea3bc8d4810e.png" attr="h" comment="" date="1452803639" name="latex77039efc6bd25fa1afbcea3bc8d4810e.png" user="zxap014" version="1" attachment="latex13af512f02812991f946f3d04c7543e3.png" attr="h" comment="" date="1452803640" name="latex13af512f02812991f946f3d04c7543e3.png" user="zxap014" version="1" attachment="latex3ac4935b323901c8834566634d3f941b.png" attr="h" comment="" date="1452803920" name="latex3ac4935b323901c8834566634d3f941b.png" user="zxap014" version="1" attachment="latex692bd06395f79b76b9565253eab2380e.png" attr="h" comment="" date="1452803921" name="latex692bd06395f79b76b9565253eab2380e.png" user="zxap014" version="1" attachment="latex5f494a7f51ab74988ac63a6ae09c93b5.png" attr="h" comment="" date="1452803921" name="latex5f494a7f51ab74988ac63a6ae09c93b5.png" user="zxap014" version="1" attachment="latex57e1a35350428aef54524306ab127632.png" attr="h" comment="" date="1452803921" name="latex57e1a35350428aef54524306ab127632.png" user="zxap014" version="1" attachment="latex65fa69e87f5b4352100c37e6f353016c.png" attr="h" comment="" date="1452803921" name="latex65fa69e87f5b4352100c37e6f353016c.png" user="zxap014" version="1" attachment="latex15c4cce6584d7462fbd78b9743ddaacf.png" attr="h" comment="" date="1452803921" name="latex15c4cce6584d7462fbd78b9743ddaacf.png" user="zxap014" version="1" attachment="latexa4de15709fa8a3f8e1d41cc33f5fc4e0.png" attr="h" comment="" date="1452803921" name="latexa4de15709fa8a3f8e1d41cc33f5fc4e0.png" user="zxap014" version="1" attachment="latexa8db5af7bf00df8a6cb2e5195a0334d2.png" attr="h" comment="" date="1452804244" name="latexa8db5af7bf00df8a6cb2e5195a0334d2.png" user="zxap014" version="1" attachment="latex748f0b68e770bdb0a987fc5d7afbce84.png" attr="h" comment="" date="1452804244" name="latex748f0b68e770bdb0a987fc5d7afbce84.png" user="zxap014" version="1" attachment="latex6f5d2de59ac4246801002d9f1c6f2bcf.png" attr="h" comment="" date="1452804244" name="latex6f5d2de59ac4246801002d9f1c6f2bcf.png" user="zxap014" version="1" attachment="latex0991a2c909c138cc7ac0837469363678.png" attr="h" comment="" date="1452804245" name="latex0991a2c909c138cc7ac0837469363678.png" user="zxap014" version="1" attachment="latexac47bbb4e0810287d874606c0a90686c.png" attr="h" comment="" date="1452804245" name="latexac47bbb4e0810287d874606c0a90686c.png" user="zxap014" version="1" attachment="latex443fc72ce248557c59acdbc7bcd46185.png" attr="h" comment="" date="1452804245" name="latex443fc72ce248557c59acdbc7bcd46185.png" user="zxap014" version="1" attachment="latexdf5a1b0fdd8b199e23d17bf7c8f71245.png" attr="h" comment="" date="1452804245" name="latexdf5a1b0fdd8b199e23d17bf7c8f71245.png" user="zxap014" version="1" attachment="latex0f02f753d4ff10388f0f33d7f857b5d9.png" attr="h" comment="" date="1452804245" name="latex0f02f753d4ff10388f0f33d7f857b5d9.png" user="zxap014" version="1" attachment="latex1e6feae7f2af4455d204e39199b43c5d.png" attr="h" comment="" date="1452804245" name="latex1e6feae7f2af4455d204e39199b43c5d.png" user="zxap014" version="1" attachment="latex2751e1e3f5f019842265aebf4f8381be.png" attr="h" comment="" date="1452804245" name="latex2751e1e3f5f019842265aebf4f8381be.png" user="zxap014" version="1" attachment="latexa96edb988d74698ee8a7b02b75c27728.png" attr="h" comment="" date="1452804245" name="latexa96edb988d74698ee8a7b02b75c27728.png" user="zxap014" version="1" attachment="latex1de82062afa31fe0c3f9716b81cf056b.png" attr="h" comment="" date="1452804245" name="latex1de82062afa31fe0c3f9716b81cf056b.png" user="zxap014" version="1" attachment="latexd6e67894e3d21074bd1ca2793fe1bbcf.png" attr="h" comment="" date="1452804246" name="latexd6e67894e3d21074bd1ca2793fe1bbcf.png" user="zxap014" version="1" attachment="latexacba13189d2c8a790b59469857f92864.png" attr="h" comment="" date="1452804246" name="latexacba13189d2c8a790b59469857f92864.png" user="zxap014" version="1" attachment="latex7603b29ae16ad93ca012826c77fcc1e8.png" attr="h" comment="" date="1452804246" name="latex7603b29ae16ad93ca012826c77fcc1e8.png" user="zxap014" version="1" attachment="latex6fb13499d69248af4e445596037315a9.png" attr="h" comment="" date="1452804246" name="latex6fb13499d69248af4e445596037315a9.png" user="zxap014" version="1" attachment="latex1460d1d9337a7ab95a5fb9f129c86429.png" attr="h" comment="" date="1452804246" name="latex1460d1d9337a7ab95a5fb9f129c86429.png" user="zxap014" version="1" attachment="latex48e07142b9b00e8fe57832a19447254b.png" attr="h" comment="" date="1452804246" name="latex48e07142b9b00e8fe57832a19447254b.png" user="zxap014" version="1" attachment="latex104f9b2b1c5ebe79b8850d3b84db5135.png" attr="h" comment="" date="1452804246" name="latex104f9b2b1c5ebe79b8850d3b84db5135.png" user="zxap014" version="1" attachment="latex1f4d96afa7698d9cc55622761b9d952f.png" attr="h" comment="" date="1452804246" name="latex1f4d96afa7698d9cc55622761b9d952f.png" user="zxap014" version="1" attachment="latex16727e49767daa49cb47c5f62ee5ba74.png" attr="h" comment="" date="1452804246" name="latex16727e49767daa49cb47c5f62ee5ba74.png" user="zxap014" version="1" attachment="latexd9d1a187e14ae4b59ea046b7e5f10fca.png" attr="h" comment="" date="1452804246" name="latexd9d1a187e14ae4b59ea046b7e5f10fca.png" user="zxap014" version="1" attachment="latex6ead7341b4d1a804b0996b16fba4d07f.png" attr="h" comment="" date="1452804247" name="latex6ead7341b4d1a804b0996b16fba4d07f.png" user="zxap014" version="1" attachment="latex85199bca4d5ba4b7dd6f0d76cd61e91c.png" attr="h" comment="" date="1452804247" name="latex85199bca4d5ba4b7dd6f0d76cd61e91c.png" user="zxap014" version="1" attachment="latex3f307b6209d8d03996efd0f24c55e2ae.png" attr="h" comment="" date="1452805398" name="latex3f307b6209d8d03996efd0f24c55e2ae.png" user="zxap014" version="1" attachment="latexfe80bf9a1d72da94143732bf9338fffa.png" attr="h" comment="" date="1452805399" name="latexfe80bf9a1d72da94143732bf9338fffa.png" user="zxap014" version="1" attachment="latexeb81fd2f5a1552f70d3db5e483e2c2ce.png" attr="h" comment="" date="1452805399" name="latexeb81fd2f5a1552f70d3db5e483e2c2ce.png" user="zxap014" version="1" attachment="latex62a47f352751572198f9cf7f780f6e57.png" attr="h" comment="" date="1452805399" name="latex62a47f352751572198f9cf7f780f6e57.png" user="zxap014" version="1" attachment="latex50098bdc4423bf73994617856b69b78c.png" attr="h" comment="" date="1452805399" name="latex50098bdc4423bf73994617856b69b78c.png" user="zxap014" version="1" attachment="latexebb5cf071e487f7d320df5c8eda1b490.png" attr="h" comment="" date="1452805399" name="latexebb5cf071e487f7d320df5c8eda1b490.png" user="zxap014" version="1" attachment="latex955a3dd52654073c5b895b06c4f51b37.png" attr="h" comment="" date="1452805399" name="latex955a3dd52654073c5b895b06c4f51b37.png" user="zxap014" version="1" attachment="latex834bd2150b3d1f1846d39e0c4a944804.png" attr="h" comment="" date="1452805400" name="latex834bd2150b3d1f1846d39e0c4a944804.png" user="zxap014" version="1" attachment="latexad8de4c1a1939f6c479b5ab708362358.png" attr="h" comment="" date="1452805400" name="latexad8de4c1a1939f6c479b5ab708362358.png" user="zxap014" version="1" attachment="latex6656f542dbf52b9b6bcee55362787ed0.png" attr="h" comment="" date="1452805400" name="latex6656f542dbf52b9b6bcee55362787ed0.png" user="zxap014" version="1" attachment="latex4a5829f0e469b80d910bd66dd520c0d7.png" attr="h" comment="" date="1452805400" name="latex4a5829f0e469b80d910bd66dd520c0d7.png" user="zxap014" version="1" attachment="latexe77413f1890f1cc20b6ef49ce7c66903.png" attr="h" comment="" date="1452805400" name="latexe77413f1890f1cc20b6ef49ce7c66903.png" user="zxap014" version="1" attachment="latexdf7daf8505b6fa5055d98bdd54768864.png" attr="h" comment="" date="1452805400" name="latexdf7daf8505b6fa5055d98bdd54768864.png" user="zxap014" version="1" attachment="latex3c6dadb4b4a95f614cce21d899463e52.png" attr="h" comment="" date="1452805400" name="latex3c6dadb4b4a95f614cce21d899463e52.png" user="zxap014" version="1" attachment="latexc979e42c368b644371e756fe49178471.png" attr="h" comment="" date="1452805400" name="latexc979e42c368b644371e756fe49178471.png" user="zxap014" version="1" attachment="latex3f74765ff1450382b2a672ca00db7c55.png" attr="h" comment="" date="1452805401" name="latex3f74765ff1450382b2a672ca00db7c55.png" user="zxap014" version="1" attachment="latex5ff4bb3c1c341a64b95069696dfe89be.png" attr="h" comment="" date="1452805401" name="latex5ff4bb3c1c341a64b95069696dfe89be.png" user="zxap014" version="1" attachment="latex547a477bcf0736fdea6395bd12bb4f99.png" attr="h" comment="" date="1452805401" name="latex547a477bcf0736fdea6395bd12bb4f99.png" user="zxap014" version="1" attachment="latex0166807616121cd36320d00278ca9661.png" attr="h" comment="" date="1452805401" name="latex0166807616121cd36320d00278ca9661.png" user="zxap014" version="1" attachment="latexbe3d0e7988be0ca52cc4234d3822466a.png" attr="h" comment="" date="1452805401" name="latexbe3d0e7988be0ca52cc4234d3822466a.png" user="zxap014" version="1" attachment="latex939550d6d5dd99405596abf5f5de9599.png" attr="h" comment="" date="1452805401" name="latex939550d6d5dd99405596abf5f5de9599.png" user="zxap014" version="1" attachment="latexed1131c899bd8e3d0caf07f9d43bade7.png" attr="h" comment="" date="1452805401" name="latexed1131c899bd8e3d0caf07f9d43bade7.png" user="zxap014" version="1" attachment="latex870130fad6df59cda86584366a5d78db.png" attr="h" comment="" date="1452805401" name="latex870130fad6df59cda86584366a5d78db.png" user="zxap014" version="1" attachment="latex903b1ae0128b79968cd7b64361b9463c.png" attr="h" comment="" date="1452805401" name="latex903b1ae0128b79968cd7b64361b9463c.png" user="zxap014" version="1" attachment="latexfd4802c8eaf7348c291d5b88d8d0e07f.png" attr="h" comment="" date="1452805402" name="latexfd4802c8eaf7348c291d5b88d8d0e07f.png" user="zxap014" version="1" attachment="latex108639b99a78a6988da07a90e10e46bc.png" attr="h" comment="" date="1452805402" name="latex108639b99a78a6988da07a90e10e46bc.png" user="zxap014" version="1" attachment="latex36b2de544e4ed31d7b16685fd46d4df9.png" attr="h" comment="" date="1452805402" name="latex36b2de544e4ed31d7b16685fd46d4df9.png" user="zxap014" version="1" attachment="latex91932de15bf1e05c9a15eb1522613672.png" attr="h" comment="" date="1452805402" name="latex91932de15bf1e05c9a15eb1522613672.png" user="zxap014" version="1" attachment="latexe2b2018c80f65d28543c9a3ee08589fd.png" attr="h" comment="" date="1452805402" name="latexe2b2018c80f65d28543c9a3ee08589fd.png" user="zxap014" version="1" attachment="latex8a2c964a70ca179fddbdb0b5a800818b.png" attr="h" comment="" date="1452805402" name="latex8a2c964a70ca179fddbdb0b5a800818b.png" user="zxap014" version="1" attachment="latex97f52cc96ce093f09de06291e437399d.png" attr="h" comment="" date="1452805402" name="latex97f52cc96ce093f09de06291e437399d.png" user="zxap014" version="1" attachment="latexb7243637df846e5922e9722ad0315ae0.png" attr="h" comment="" date="1452805402" name="latexb7243637df846e5922e9722ad0315ae0.png" user="zxap014" version="1" attachment="latex8fc115f0fc580e02fbc97b35f90c337e.png" attr="h" comment="" date="1452805403" name="latex8fc115f0fc580e02fbc97b35f90c337e.png" user="zxap014" version="1" attachment="latex9f776835c4ecc6e12c9a3c2115b962e5.png" attr="h" comment="" date="1452805403" name="latex9f776835c4ecc6e12c9a3c2115b962e5.png" user="zxap014" version="1" attachment="latex6232fd1e809d21de413d19be8b641acb.png" attr="h" comment="" date="1452805403" name="latex6232fd1e809d21de413d19be8b641acb.png" user="zxap014" version="1" attachment="latexa48f8765d051aab5a9702f6a5ab61b7a.png" attr="h" comment="" date="1452805403" name="latexa48f8765d051aab5a9702f6a5ab61b7a.png" user="zxap014" version="1" attachment="latex63dfac519476134233fa316a7f4c530b.png" attr="h" comment="" date="1452805403" name="latex63dfac519476134233fa316a7f4c530b.png" user="zxap014" version="1" attachment="latex26c60d32b1687ba6b7b91b8cbbcc9191.png" attr="h" comment="" date="1452805403" name="latex26c60d32b1687ba6b7b91b8cbbcc9191.png" user="zxap014" version="1" attachment="latexc90264e54bd971c2d1fa8fb6552d2551.png" attr="h" comment="" date="1452805403" name="latexc90264e54bd971c2d1fa8fb6552d2551.png" user="zxap014" version="1" attachment="latexca2b42bd0e11b404e11b86daa24f2c87.png" attr="h" comment="" date="1452805403" name="latexca2b42bd0e11b404e11b86daa24f2c87.png" user="zxap014" version="1" attachment="latexd2cb5ab31a5b6e7eae8eb43d9bce2edd.png" attr="h" comment="" date="1452805404" name="latexd2cb5ab31a5b6e7eae8eb43d9bce2edd.png" user="zxap014" version="1" attachment="latex657f37490f4abd05098e98380cb2f4ec.png" attr="h" comment="" date="1452805404" name="latex657f37490f4abd05098e98380cb2f4ec.png" user="zxap014" version="1" attachment="latexea59f1a115f575f220ec918ed7701aa2.png" attr="h" comment="" date="1452805495" name="latexea59f1a115f575f220ec918ed7701aa2.png" user="zxap014" version="1" attachment="latexbb0b8cef2877c9470839627b0a140a8a.png" attr="h" comment="" date="1452805496" name="latexbb0b8cef2877c9470839627b0a140a8a.png" user="zxap014" version="1"

#### Revision 1014 Jan 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 66 to 66
%ENDLATEX%
• We cropped an area of 20x20 pixels around each star and found the standard deviations within the area.
• Histograms below are of the ratio of for Location 1 images:
Changed:
<
<
>
>
• for Location 1 images:
• Ratio of for Location 2 images:
• for Location 2 images:
• Therefore for both locations we decided to set the requirement that an objects is only a star if its ratio:
-- ElenaCukanovaite - 11 Jan 2016
Line: 98 to 110

 META FILEATTACHMENT attachment="latexf6c3d0a704b97e5e951289a087d84905.png" attr="h" comment="" date="1452703930" name="latexf6c3d0a704b97e5e951289a087d84905.png" user="zxap012" version="1" attachment="latex31059ec47f2cfa24cf629a760cced8de.png" attr="h" comment="" date="1452703930" name="latex31059ec47f2cfa24cf629a760cced8de.png" user="zxap012" version="1" attachment="latex3ef82bbef12e2dfd934e2a1fbd5c0c4a.png" attr="h" comment="" date="1452703930" name="latex3ef82bbef12e2dfd934e2a1fbd5c0c4a.png" user="zxap012" version="1"
>
>
 META FILEATTACHMENT attachment="sigma_x_loc_1.png" attr="" comment="" date="1452788039" name="sigma_x_loc_1.png" path="sigma_x_loc_1.png" size="67878" user="zxap014" version="1" attachment="location_2.png" attr="" comment="" date="1452788166" name="location_2.png" path="location_2.png" size="51693" user="zxap014" version="1" attachment="sigma_x_loc_2.png" attr="" comment="" date="1452788220" name="sigma_x_loc_2.png" path="sigma_x_loc_2.png" size="47938" user="zxap014" version="1" attachment="latex8906deced046f2a1528ba6815ff3285e.png" attr="h" comment="" date="1452788248" name="latex8906deced046f2a1528ba6815ff3285e.png" user="zxap014" version="1" attachment="latex970d0dad73601c12ab733877a694729d.png" attr="h" comment="" date="1452788520" name="latex970d0dad73601c12ab733877a694729d.png" user="zxap014" version="1"

#### Revision 913 Jan 2016 - DavidHadden

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 36 to 36

 reduced 68600 1150 36500 1277
• In general, a good should be around equal to the number of degrees of freedom,

• In the above case, they all seem to give very poor agreement between the fit and the data, as the values are all very large in comparison to the number of degrees of freedom.

>
>

### Testing curve_fit for better error estimation

• In order to better obtain errors from curve_fit (as it tends to give outlandish errors), we set out to create a Gaussian of known parameters, and then fit a general Gaussian curve to this.
 Fractional Error When Guessing Incorectly Fractional Error When Guessing Exactly Correlation Coefficient 0 Amplitude 0 X Centre 0 Y Centre 0 0 0 Offset 0

### Setting the threshold

• We decided to use our mode image to set the threshold.
Line: 82 to 92

 META FILEATTACHMENT attachment="latex13817a18e09b89ec91da62d244cd50fd.png" attr="h" comment="" date="1452698022" name="latex13817a18e09b89ec91da62d244cd50fd.png" user="zxap012" version="1" attachment="stan_dev_n_t.png" attr="" comment="" date="1452701616" name="stan_dev_n_t.png" path="stan_dev_n_t.png" size="71821" user="zxap014" version="1" attachment="latex52f206289e28eb74c4ebe60962d3aeee.png" attr="h" comment="" date="1452701656" name="latex52f206289e28eb74c4ebe60962d3aeee.png" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="latex03bd5c628ea405f98c42cdf729ae7d6c.png" attr="h" comment="" date="1452703889" name="latex03bd5c628ea405f98c42cdf729ae7d6c.png" user="zxap012" version="1" attachment="latex48029a73862d18811ef42a352f3d20a9.png" attr="h" comment="" date="1452703889" name="latex48029a73862d18811ef42a352f3d20a9.png" user="zxap012" version="1" attachment="latex9ea3b761b6756b1abd8e33db2ab29e28.png" attr="h" comment="" date="1452703889" name="latex9ea3b761b6756b1abd8e33db2ab29e28.png" user="zxap012" version="1" attachment="latexf6c3d0a704b97e5e951289a087d84905.png" attr="h" comment="" date="1452703930" name="latexf6c3d0a704b97e5e951289a087d84905.png" user="zxap012" version="1" attachment="latex31059ec47f2cfa24cf629a760cced8de.png" attr="h" comment="" date="1452703930" name="latex31059ec47f2cfa24cf629a760cced8de.png" user="zxap012" version="1" attachment="latex3ef82bbef12e2dfd934e2a1fbd5c0c4a.png" attr="h" comment="" date="1452703930" name="latex3ef82bbef12e2dfd934e2a1fbd5c0c4a.png" user="zxap012" version="1"

#### Revision 813 Jan 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 55 to 55
%ENDLATEX%
• We cropped an area of 20x20 pixels around each star and found the standard deviations within the area.
Changed:
<
<
>
>
• Histograms below are of the ratio of for Location 1 images:
-- ElenaCukanovaite - 11 Jan 2016

 META FILEATTACHMENT attachment="latex814eaaf2163b7785af28f5b62b89870c.png" attr="h" comment="" date="1452516527" name="latex814eaaf2163b7785af28f5b62b89870c.png" user="zxap014" version="1"
Line: 78 to 80

 META FILEATTACHMENT attachment="None_1_min.png" attr="" comment="" date="1452690783" name="None_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_min.png" size="112431" user="zxap012" version="1" attachment="monte_carlo.png" attr="" comment="" date="1452692332" name="monte_carlo.png" path="monte_carlo.png" size="284643" user="zxap014" version="1" attachment="latex13817a18e09b89ec91da62d244cd50fd.png" attr="h" comment="" date="1452698022" name="latex13817a18e09b89ec91da62d244cd50fd.png" user="zxap012" version="1"
>
>
 META FILEATTACHMENT attachment="stan_dev_n_t.png" attr="" comment="" date="1452701616" name="stan_dev_n_t.png" path="stan_dev_n_t.png" size="71821" user="zxap014" version="1" attachment="latex52f206289e28eb74c4ebe60962d3aeee.png" attr="h" comment="" date="1452701656" name="latex52f206289e28eb74c4ebe60962d3aeee.png" user="zxap014" version="1"

#### Revision 713 Jan 2016 - DavidHadden

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 27 to 27

>
>

### Introductory estimation

• Then set about finding for the above stars for no filters and the green filter.

 Visual Brightest Star Visual Dimmest Star Green Filter Brightest Star Green Filter Dimmest Star Degrees of Freedom 53 53 53 53 3640000 61100 1930000 67700 reduced 68600 1150 36500 1277
• In general, a good should be around equal to the number of degrees of freedom,

• In the above case, they all seem to give very poor agreement between the fit and the data, as the values are all very large in comparison to the number of degrees of freedom.

### Setting the threshold

• We decided to use our mode image to set the threshold.
• The mode image is shown below:
Line: 67 to 77

 META FILEATTACHMENT attachment="None_1_max.png" attr="" comment="" date="1452690774" name="None_1_max.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_max.png" size="115849" user="zxap012" version="1" attachment="None_1_min.png" attr="" comment="" date="1452690783" name="None_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_min.png" size="112431" user="zxap012" version="1" attachment="monte_carlo.png" attr="" comment="" date="1452692332" name="monte_carlo.png" path="monte_carlo.png" size="284643" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="latex13817a18e09b89ec91da62d244cd50fd.png" attr="h" comment="" date="1452698022" name="latex13817a18e09b89ec91da62d244cd50fd.png" user="zxap012" version="1"

#### Revision 613 Jan 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 16 to 16

• To experiment we uniformly generated values for and for each individual star.
• We then for each pixel in the image calculated the Gaussian value at that particular pixel for each star.
• Resulting image for 50 stars is shown below:
>
>

### Comparing stars between different filters

• In order to compare the quality of our data, we looked at the data with the largest signal to noise ratio, in this case data taken with a green filter.
• For this, we found the largest and smallest stars, as determined by the amount of pixels which constituted them.
• The coordinates of the stars were recorded, so they can later be compared to the same stars in other filters.
• Below fits can be seen for the maximum and minimum size stars for the green filter:
Changed:
<
<
>
>

### Setting the threshold

• We decided to use our mode image to set the threshold.
• The mode image is shown below:
Line: 43 to 44
$$\sigma_y = \sqrt{\frac{\sum_{i = 1}^N (y_i - y_{\rm center})^2p_i}{\sum_{i=1}^N p_i}}$$

%ENDLATEX%

Changed:
<
<
• We found the following:
>
>
• We cropped an area of 20x20 pixels around each star and found the standard deviations within the area.
-- ElenaCukanovaite - 11 Jan 2016

 META FILEATTACHMENT attachment="latex814eaaf2163b7785af28f5b62b89870c.png" attr="h" comment="" date="1452516527" name="latex814eaaf2163b7785af28f5b62b89870c.png" user="zxap014" version="1"
Line: 64 to 66

 META FILEATTACHMENT attachment="red_1_min.png" attr="" comment="" date="1452690764" name="red_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\red_1_min.png" size="120231" user="zxap012" version="1" attachment="None_1_max.png" attr="" comment="" date="1452690774" name="None_1_max.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_max.png" size="115849" user="zxap012" version="1" attachment="None_1_min.png" attr="" comment="" date="1452690783" name="None_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_min.png" size="112431" user="zxap012" version="1"
>
>
 META FILEATTACHMENT attachment="monte_carlo.png" attr="" comment="" date="1452692332" name="monte_carlo.png" path="monte_carlo.png" size="284643" user="zxap014" version="1"

#### Revision 513 Jan 2016 - DavidHadden

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 21 to 21

• For this, we found the largest and smallest stars, as determined by the amount of pixels which constituted them.
• The coordinates of the stars were recorded, so they can later be compared to the same stars in other filters.
• Below fits can be seen for the maximum and minimum size stars for the green filter:
Deleted:
<
<

Changed:
<
<
>
>

### Setting the threshold

• We decided to use our mode image to set the threshold.
• The mode image is shown below:
Line: 50 to 52

 META FILEATTACHMENT attachment="latexca7d89b9f7485aa401282d7fb7de4764.png" attr="h" comment="" date="1452516941" name="latexca7d89b9f7485aa401282d7fb7de4764.png" user="zxap014" version="1" attachment="latex271f02032cbb9e3d532d351c7f56632c.png" attr="h" comment="" date="1452517091" name="latex271f02032cbb9e3d532d351c7f56632c.png" user="zxap014" version="1" attachment="latex6bcf65c3353b2591f79fa6700e7a1a41.png" attr="h" comment="" date="1452517091" name="latex6bcf65c3353b2591f79fa6700e7a1a41.png" user="zxap014" version="1"
Deleted:
<
<
 META FILEATTACHMENT attachment="green_1_min.png" attr="" comment="" date="1452528352" name="green_1_min.png" path="green_1_min.png" size="106945" user="zxap012" version="1" attachment="green_1_max.png" attr="" comment="" date="1452528363" name="green_1_max.png" path="green_1_max.png" size="121556" user="zxap012" version="1"

 META FILEATTACHMENT attachment="latex4428cf23170176410a86bab987365e82.png" attr="h" comment="" date="1452530379" name="latex4428cf23170176410a86bab987365e82.png" user="zxap014" version="1" attachment="30_by_30.png" attr="" comment="" date="1452684705" name="30_by_30.png" path="30_by_30.png" size="54145" user="zxap014" version="1" attachment="mode_sec.png" attr="" comment="" date="1452685331" name="mode_sec.png" path="mode_sec.png" size="47403" user="zxap014" version="1"
Deleted:
<
<
 META FILEATTACHMENT attachment="latex81d3f9ff4723bbb8aff4fbb5c4a81d23.png" attr="h" comment="" date="1452685532" name="latex81d3f9ff4723bbb8aff4fbb5c4a81d23.png" user="zxap014" version="1"

 META FILEATTACHMENT attachment="latex1b62d92470af61325478d128b7895042.png" attr="h" comment="" date="1452685555" name="latex1b62d92470af61325478d128b7895042.png" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="blue_1_max.png" attr="" comment="" date="1452690714" name="blue_1_max.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\blue_1_max.png" size="122947" user="zxap012" version="1" attachment="blue_1_min.png" attr="" comment="" date="1452690725" name="blue_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\blue_1_min.png" size="118640" user="zxap012" version="1" attachment="green_1_max.png" attr="" comment="" date="1452690738" name="green_1_max.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\green_1_max.png" size="124301" user="zxap012" version="1" attachment="green_1_min.png" attr="" comment="" date="1452690746" name="green_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\green_1_min.png" size="109560" user="zxap012" version="1" attachment="red_1_max.png" attr="" comment="" date="1452690756" name="red_1_max.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\red_1_max.png" size="122642" user="zxap012" version="1" attachment="red_1_min.png" attr="" comment="" date="1452690764" name="red_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\red_1_min.png" size="120231" user="zxap012" version="1" attachment="None_1_max.png" attr="" comment="" date="1452690774" name="None_1_max.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_max.png" size="115849" user="zxap012" version="1" attachment="None_1_min.png" attr="" comment="" date="1452690783" name="None_1_min.png" path="Y:\Major Project\20151113_DH_EC\Analysis\Loaction_1\max_min_tests\None_1_min.png" size="112431" user="zxap012" version="1"

#### Revision 413 Jan 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 24 to 24

>
>

### Setting the threshold

• We decided to use our mode image to set the threshold.
• The mode image is shown below:
• We took the Poissonian error (the square root) of each section.
• We used the error as our threshold by setting the threshold:

### Standard deviation of the stars

• We previously tried to show that all stars in the image have the same standard deviations.
• To find our standard deviations we did the following:
Line: 44 to 53

 META FILEATTACHMENT attachment="green_1_min.png" attr="" comment="" date="1452528352" name="green_1_min.png" path="green_1_min.png" size="106945" user="zxap012" version="1" attachment="green_1_max.png" attr="" comment="" date="1452528363" name="green_1_max.png" path="green_1_max.png" size="121556" user="zxap012" version="1" attachment="latex4428cf23170176410a86bab987365e82.png" attr="h" comment="" date="1452530379" name="latex4428cf23170176410a86bab987365e82.png" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="30_by_30.png" attr="" comment="" date="1452684705" name="30_by_30.png" path="30_by_30.png" size="54145" user="zxap014" version="1" attachment="mode_sec.png" attr="" comment="" date="1452685331" name="mode_sec.png" path="mode_sec.png" size="47403" user="zxap014" version="1" attachment="latex81d3f9ff4723bbb8aff4fbb5c4a81d23.png" attr="h" comment="" date="1452685532" name="latex81d3f9ff4723bbb8aff4fbb5c4a81d23.png" user="zxap014" version="1" attachment="latex1b62d92470af61325478d128b7895042.png" attr="h" comment="" date="1452685555" name="latex1b62d92470af61325478d128b7895042.png" user="zxap014" version="1"

#### Revision 311 Jan 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 21 to 21

• For this, we found the largest and smallest stars, as determined by the amount of pixels which constituted them.
• The coordinates of the stars were recorded, so they can later be compared to the same stars in other filters.
• Below fits can be seen for the maximum and minimum size stars for the green filter:
Changed:
<
<
>
>

Changed:
<
<
>
>

### Standard deviation of the stars

• We previously tried to show that all stars in the image have the same standard deviations.
• To find our standard deviations we did the following:
%BEGINLATEX% $$\sigma_x = \sqrt{\frac{\sum_{i = 1}^N (x_i - x_{\rm center})^2 p_i}{\sum_{i=1}^N p_i}}$$

>
>
$$\sigma_y = \sqrt{\frac{\sum_{i = 1}^N (y_i - y_{\rm center})^2p_i}{\sum_{i=1}^N p_i}}$$

%ENDLATEX%

• We found the following:
-- ElenaCukanovaite - 11 Jan 2016

 META FILEATTACHMENT attachment="latex814eaaf2163b7785af28f5b62b89870c.png" attr="h" comment="" date="1452516527" name="latex814eaaf2163b7785af28f5b62b89870c.png" user="zxap014" version="1"
Line: 35 to 43

 META FILEATTACHMENT attachment="latex6bcf65c3353b2591f79fa6700e7a1a41.png" attr="h" comment="" date="1452517091" name="latex6bcf65c3353b2591f79fa6700e7a1a41.png" user="zxap014" version="1" attachment="green_1_min.png" attr="" comment="" date="1452528352" name="green_1_min.png" path="green_1_min.png" size="106945" user="zxap012" version="1" attachment="green_1_max.png" attr="" comment="" date="1452528363" name="green_1_max.png" path="green_1_max.png" size="121556" user="zxap012" version="1"
>
>
 META FILEATTACHMENT attachment="latex4428cf23170176410a86bab987365e82.png" attr="h" comment="" date="1452530379" name="latex4428cf23170176410a86bab987365e82.png" user="zxap014" version="1"

#### Revision 211 Jan 2016 - DavidHadden

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
Line: 16 to 16

• To experiment we uniformly generated values for and for each individual star.
• We then for each pixel in the image calculated the Gaussian value at that particular pixel for each star.
• Resulting image for 50 stars is shown below:
>
>

### Comparing stars between different filters

• In order to compare the quality of our data, we looked at the data with the largest signal to noise ratio, in this case data taken with a green filter.
• For this, we found the largest and smallest stars, as determined by the amount of pixels which constituted them.
• The coordinates of the stars were recorded, so they can later be compared to the same stars in other filters.
• Below fits can be seen for the maximum and minimum size stars for the green filter:

-- ElenaCukanovaite - 11 Jan 2016

 META FILEATTACHMENT attachment="latex814eaaf2163b7785af28f5b62b89870c.png" attr="h" comment="" date="1452516527" name="latex814eaaf2163b7785af28f5b62b89870c.png" user="zxap014" version="1"
Line: 24 to 33

 META FILEATTACHMENT attachment="latexca7d89b9f7485aa401282d7fb7de4764.png" attr="h" comment="" date="1452516941" name="latexca7d89b9f7485aa401282d7fb7de4764.png" user="zxap014" version="1" attachment="latex271f02032cbb9e3d532d351c7f56632c.png" attr="h" comment="" date="1452517091" name="latex271f02032cbb9e3d532d351c7f56632c.png" user="zxap014" version="1" attachment="latex6bcf65c3353b2591f79fa6700e7a1a41.png" attr="h" comment="" date="1452517091" name="latex6bcf65c3353b2591f79fa6700e7a1a41.png" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="green_1_min.png" attr="" comment="" date="1452528352" name="green_1_min.png" path="green_1_min.png" size="106945" user="zxap012" version="1" attachment="green_1_max.png" attr="" comment="" date="1452528363" name="green_1_max.png" path="green_1_max.png" size="121556" user="zxap012" version="1"

#### Revision 111 Jan 2016 - ElenaCukanovaite

Line: 1 to 1
>
>
 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Monte Carlo simulation of a 2D star field

• We have done some prelimenary work on the simulation.
• To generate the positions of stars we uniformly generated x and y pixel positions.
• To create a background we generated Poisson distributed random variables around the mean of 160, which corresponded to a typical value for images of Location 1.
• For each position we wanted to create a Gaussian distribution that represents a star centered on the star's x,y position.
• As the standard deviations should be the same for each star we decided to generate only one value for for all stars. We also assumed that the standard deviations are the same in both directions so we generated only one for and .
• To experiment we uniformly generated values for and for each individual star.
• We then for each pixel in the image calculated the Gaussian value at that particular pixel for each star.
• Resulting image for 50 stars is shown below:
-- ElenaCukanovaite - 11 Jan 2016

 META FILEATTACHMENT attachment="latex814eaaf2163b7785af28f5b62b89870c.png" attr="h" comment="" date="1452516527" name="latex814eaaf2163b7785af28f5b62b89870c.png" user="zxap014" version="1" attachment="latexd881f63bae4b1f73cb24e393705ff10b.png" attr="h" comment="" date="1452516941" name="latexd881f63bae4b1f73cb24e393705ff10b.png" user="zxap014" version="1" attachment="latex9d0de473d0caf20f0189c0fa2811e1b0.png" attr="h" comment="" date="1452516941" name="latex9d0de473d0caf20f0189c0fa2811e1b0.png" user="zxap014" version="1" attachment="latexca7d89b9f7485aa401282d7fb7de4764.png" attr="h" comment="" date="1452516941" name="latexca7d89b9f7485aa401282d7fb7de4764.png" user="zxap014" version="1" attachment="latex271f02032cbb9e3d532d351c7f56632c.png" attr="h" comment="" date="1452517091" name="latex271f02032cbb9e3d532d351c7f56632c.png" user="zxap014" version="1" attachment="latex6bcf65c3353b2591f79fa6700e7a1a41.png" attr="h" comment="" date="1452517091" name="latex6bcf65c3353b2591f79fa6700e7a1a41.png" user="zxap014" version="1"

Copyright © 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding RHUL Physics Department TWiki? Send feedback