# Difference: 20160214MSciPhotometryAnalysisObs3part2 (1 vs. 2)

#### Revision 218 Feb 2016 - ElenaCukanovaite

Line: 1 to 1

 META TOPICPARENT name="StewartBoogertPhotometry2015"

### Testing the fitting

• First we wanted to test that the function we are fitting is actually normalised.
Line: 72 to 72

 1.71e+06 3222 Degrees of freedom 3591
>
>

### Matching stars between images

• Used Pythagoras to estimate the closest star image to image. Only used stars whose were between 0.5 < < 1.5
• Went through images and found the image with highest number of stars.
• This was the base image that all other images were compared to.
• It worked fine for images in focus and blue filter since they were taken close enough in time for individual stars to not have moved much.
• In red filter however, there was double matching, when more than one star was matched to the same base star as shown below:

• Improved the algortihm by:
• first calculate the minimum distances, some of them will be wronged as stars will be wrongly matched.
• Use those matches to work out the median x and y shifts.
• Then check between matches if the x and y shifts match the median x and y shifts to 1 pixel, if they do its a match.
• Below is the same red image and matches:

### HR diagram for Cr39

• Using the matching and ML fits we got the following HR diagrams:

• The stars in V-R HR that seem to have the same V-R are indexed 26 and 27. The star with V close to -8 is indexed 2, it appeared 3 times. Star with V close to -7.5 appeared 3 times.
• The HR diagrams are terrible.
• From SIMBAD they should look like this:

-- ElenaCukanovaite - 14 Feb 2016
Line: 124 to 163

 META FILEATTACHMENT attachment="latexe95444d662b6089e30433993bcc95746.png" attr="h" comment="" date="1455479439" name="latexe95444d662b6089e30433993bcc95746.png" user="zxap014" version="1" attachment="latexb80262d0b0ab437f7d128be6dc25f7c0.png" attr="h" comment="" date="1455479440" name="latexb80262d0b0ab437f7d128be6dc25f7c0.png" user="zxap014" version="1" attachment="latex92101250b7e949896dc7365a6d63e4e4.png" attr="h" comment="" date="1455479569" name="latex92101250b7e949896dc7365a6d63e4e4.png" user="zxap014" version="1"
>
>
 META FILEATTACHMENT attachment="red_first.png" attr="" comment="" date="1455825044" name="red_first.png" path="red_first.png" size="525683" user="zxap014" version="1" attachment="red_first_closeup.png" attr="h" comment="" date="1455825265" name="red_first_closeup.png" path="red_first_closeup.png" size="351034" user="zxap014" version="2" attachment="red_first_closeup_two.png" attr="h" comment="" date="1455825198" name="red_first_closeup_two.png" path="red_first_closeup_two.png" size="445456" user="zxap014" version="2" attachment="focus_first.png" attr="h" comment="" date="1455825059" name="focus_first.png" path="focus_first.png" size="530273" user="zxap014" version="1" attachment="red_second.png" attr="" comment="" date="1455825970" name="red_second.png" path="red_second.png" size="98693" user="zxap014" version="2" attachment="red_second_closeup.png" attr="" comment="" date="1455825975" name="red_second_closeup.png" path="red_second_closeup.png" size="40623" user="zxap014" version="2" attachment="red_second_closeup_two.png" attr="h" comment="" date="1455825980" name="red_second_closeup_two.png" path="red_second_closeup_two.png" size="44952" user="zxap014" version="2" attachment="hr_b_v.png" attr="h" comment="" date="1455826422" name="hr_b_v.png" path="hr_b_v.png" size="31012" user="zxap014" version="1" attachment="hr_v_r.png" attr="h" comment="" date="1455826427" name="hr_v_r.png" path="hr_v_r.png" size="31135" user="zxap014" version="1" attachment="B_V_simbad.png" attr="" comment="" date="1455827326" name="B_V_simbad.png" path="B_V_simbad.png" size="50510" user="zxap014" version="1" attachment="V_R_simbad.png" attr="" comment="" date="1455827330" name="V_R_simbad.png" path="V_R_simbad.png" size="34826" user="zxap014" version="1" attachment="latex9bf7a2667af9a16245d4de60644f5686.png" attr="h" comment="" date="1455827418" name="latex9bf7a2667af9a16245d4de60644f5686.png" user="zxap014" version="1" attachment="latex555fe965e3be97eb5df382bf33d5eb4e.png" attr="h" comment="" date="1455827418" name="latex555fe965e3be97eb5df382bf33d5eb4e.png" user="zxap014" version="1"

#### Revision 114 Feb 2016 - ElenaCukanovaite

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

### Testing the fitting

• First we wanted to test that the function we are fitting is actually normalised.
• The skewed Gaussian function is:
• Below are some skewed Gaussians for different skews:
• A normalized function would have values between 0 and 1, but here the values are between 0 and 14. This is because the function is summed over all columns and rows.
• If one looks at an individual column, it can be seen that the function is indeed normalised.
• Middle column:
• At the far columns one can see:

• Also, the formula was extracted from this reference ( doi: 10.1093/biomet/83.4.715 ) which derives a normalized skewed bivariate Gaussian.
• We wanted to look at a good star and see how the fit looks like in as much detail as possible, as the summation along rows and columns seems to distort the details of fitting.
• The fit is done by minimizing the in the following way:
• So we fit for 9 parameters.
• We do not integrate, instead we evaluate the function at the pixel and multiply by pixel area.

• We fitted this star:
• Column by column the fit looks like this:
• The plot parameters are shown below:
 0.0509 33.5 33.8 5.021 5.373 -1.184 -1.384 0.0014 1.71e+06 3222 Degrees of freedom 3591

-- ElenaCukanovaite - 14 Feb 2016

 META FILEATTACHMENT attachment="summed.png" attr="" comment="" date="1455477180" name="summed.png" path="summed.png" size="142981" user="zxap014" version="1" attachment="003_150.png" attr="" comment="" date="1455477197" name="003_150.png" path="003_150.png" size="105656" user="zxap014" version="1" attachment="005_250.png" attr="" comment="" date="1455477202" name="005_250.png" path="005_250.png" size="84196" user="zxap014" version="1" attachment="001_50.png" attr="" comment="" date="1455477207" name="001_50.png" path="001_50.png" size="83770" user="zxap014" version="1" attachment="crop.png" attr="" comment="" date="1455477684" name="crop.png" path="crop.png" size="59764" user="zxap014" version="1" attachment="complete_fit.png" attr="" comment="" date="1455477691" name="complete_fit.png" path="complete_fit.png" size="114384" user="zxap014" version="1" attachment="column_15.png" attr="h" comment="" date="1455478119" name="column_15.png" path="column_15.png" size="38923" user="zxap014" version="2" attachment="column_16.png" attr="h" comment="" date="1455478126" name="column_16.png" path="column_16.png" size="41241" user="zxap014" version="2" attachment="column_17.png" attr="h" comment="" date="1455478154" name="column_17.png" path="column_17.png" size="39862" user="zxap014" version="2" attachment="column_18.png" attr="h" comment="" date="1455478158" name="column_18.png" path="column_18.png" size="43716" user="zxap014" version="2" attachment="column_19.png" attr="h" comment="" date="1455478163" name="column_19.png" path="column_19.png" size="42753" user="zxap014" version="2" attachment="column_20.png" attr="h" comment="" date="1455478167" name="column_20.png" path="column_20.png" size="43058" user="zxap014" version="2" attachment="column_21.png" attr="h" comment="" date="1455478171" name="column_21.png" path="column_21.png" size="45588" user="zxap014" version="2" attachment="column_22.png" attr="h" comment="" date="1455478175" name="column_22.png" path="column_22.png" size="45782" user="zxap014" version="2" attachment="column_23.png" attr="h" comment="" date="1455478178" name="column_23.png" path="column_23.png" size="50482" user="zxap014" version="2" attachment="column_24.png" attr="h" comment="" date="1455478183" name="column_24.png" path="column_24.png" size="44552" user="zxap014" version="2" attachment="column_25.png" attr="h" comment="" date="1455478187" name="column_25.png" path="column_25.png" size="47111" user="zxap014" version="2" attachment="column_26.png" attr="h" comment="" date="1455478192" name="column_26.png" path="column_26.png" size="48647" user="zxap014" version="2" attachment="column_27.png" attr="h" comment="" date="1455478196" name="column_27.png" path="column_27.png" size="49285" user="zxap014" version="2" attachment="column_28.png" attr="h" comment="" date="1455478200" name="column_28.png" path="column_28.png" size="52047" user="zxap014" version="2" attachment="column_29.png" attr="h" comment="" date="1455478204" name="column_29.png" path="column_29.png" size="52552" user="zxap014" version="2" attachment="column_30.png" attr="h" comment="" date="1455478208" name="column_30.png" path="column_30.png" size="53859" user="zxap014" version="2" attachment="column_31.png" attr="h" comment="" date="1455478212" name="column_31.png" path="column_31.png" size="53637" user="zxap014" version="2" attachment="column_32.png" attr="h" comment="" date="1455478215" name="column_32.png" path="column_32.png" size="53943" user="zxap014" version="2" attachment="column_33.png" attr="h" comment="" date="1455478220" name="column_33.png" path="column_33.png" size="53177" user="zxap014" version="2" attachment="column_34.png" attr="h" comment="" date="1455478223" name="column_34.png" path="column_34.png" size="51721" user="zxap014" version="2" attachment="column_35.png" attr="h" comment="" date="1455478227" name="column_35.png" path="column_35.png" size="51468" user="zxap014" version="2" attachment="column_36.png" attr="h" comment="" date="1455478230" name="column_36.png" path="column_36.png" size="56885" user="zxap014" version="2" attachment="column_37.png" attr="h" comment="" date="1455478234" name="column_37.png" path="column_37.png" size="45663" user="zxap014" version="2" attachment="column_38.png" attr="h" comment="" date="1455478238" name="column_38.png" path="column_38.png" size="51061" user="zxap014" version="2" attachment="column_39.png" attr="h" comment="" date="1455478242" name="column_39.png" path="column_39.png" size="48392" user="zxap014" version="2" attachment="column_40.png" attr="h" comment="" date="1455478246" name="column_40.png" path="column_40.png" size="44720" user="zxap014" version="2" attachment="column_41.png" attr="h" comment="" date="1455478249" name="column_41.png" path="column_41.png" size="43718" user="zxap014" version="2" attachment="column_42.png" attr="h" comment="" date="1455478254" name="column_42.png" path="column_42.png" size="42595" user="zxap014" version="2" attachment="column_43.png" attr="h" comment="" date="1455478257" name="column_43.png" path="column_43.png" size="40670" user="zxap014" version="2" attachment="column_44.png" attr="h" comment="" date="1455478262" name="column_44.png" path="column_44.png" size="41082" user="zxap014" version="2" attachment="column_45.png" attr="h" comment="" date="1455478266" name="column_45.png" path="column_45.png" size="42211" user="zxap014" version="2" attachment="latexe69822a46019b6a4e4cd8d38f0184b29.png" attr="h" comment="" date="1455478614" name="latexe69822a46019b6a4e4cd8d38f0184b29.png" user="zxap014" version="1" attachment="latex2f796cc9e36e7a28f035aa71bffa8204.png" attr="h" comment="" date="1455479436" name="latex2f796cc9e36e7a28f035aa71bffa8204.png" user="zxap014" version="1" attachment="latexc1ee5e4b7bbc4b3de0498005a0d5d5dd.png" attr="h" comment="" date="1455479436" name="latexc1ee5e4b7bbc4b3de0498005a0d5d5dd.png" user="zxap014" version="1" attachment="latex2c1db3ab700e3511de9f696230f6dca6.png" attr="h" comment="" date="1455479437" name="latex2c1db3ab700e3511de9f696230f6dca6.png" user="zxap014" version="1" attachment="latex609b1f9aaf260c1506b685af6cec7404.png" attr="h" comment="" date="1455479437" name="latex609b1f9aaf260c1506b685af6cec7404.png" user="zxap014" version="1" attachment="latexe7b9e3aacf9c90156004f9381aba5cc4.png" attr="h" comment="" date="1455479437" name="latexe7b9e3aacf9c90156004f9381aba5cc4.png" user="zxap014" version="1" attachment="latex2c74254b08260f2622cce7e27b09e566.png" attr="h" comment="" date="1455479438" name="latex2c74254b08260f2622cce7e27b09e566.png" user="zxap014" version="1" attachment="latex44fd95b4604601049034012851eedcf0.png" attr="h" comment="" date="1455479438" name="latex44fd95b4604601049034012851eedcf0.png" user="zxap014" version="1" attachment="latexa023df2e6bfceb5d5e63e9f68f847525.png" attr="h" comment="" date="1455479439" name="latexa023df2e6bfceb5d5e63e9f68f847525.png" user="zxap014" version="1" attachment="latexe95444d662b6089e30433993bcc95746.png" attr="h" comment="" date="1455479439" name="latexe95444d662b6089e30433993bcc95746.png" user="zxap014" version="1" attachment="latexb80262d0b0ab437f7d128be6dc25f7c0.png" attr="h" comment="" date="1455479440" name="latexb80262d0b0ab437f7d128be6dc25f7c0.png" user="zxap014" version="1" attachment="latex92101250b7e949896dc7365a6d63e4e4.png" attr="h" comment="" date="1455479569" name="latex92101250b7e949896dc7365a6d63e4e4.png" user="zxap014" version="1"

Copyright © 2008-2022 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