Nova ( on sum 20 ) algorithm:

1. Ciecia :

2. Odrzucanie gdy w blasku jasnej gwiazdy

•  wyznaczenie promienia odrzucania gdy bliska gwiazda, trzeba to zrobic tak :
  
  uzywajac szybkiej fotometrii , wyznaczyc max_radius vs normalized magnitudo , hmm trzeba wiec zrovic matchowanie
  i normalizacje z katalogu ....
  Ok tymczasowo uzyje star.shape zeby wyliczyc max radius i wpisze to opcjonalnie do maga ( tzn wyliczanie bedzie opcjonalne ! )
  dodalem do piphoto opcje -calc_max_radius, dzieki ktorej mozna zapisac do maga w polu shape maxymalny promien gwiazdy
  ( liczony od srodka ciezkosci gwiazdy )
  tymaczsowo zmienilem program piaddast2 ( zapisalem wersje w pliku main.cpp.radius ) zeby pole sharp wpisywal do pola
  Measurments.error i z tego moge potem zrobic wykres zaleznosci magnitude vs error !
  Test tego katalogowania zrobilem do bazy radius_test , wykonalem to w /pi20/msok/badania/nova_flares/20060911/radius_test/
  Plot wykonalem skryptem : radius_vs_mag!
  Dostalem dosc ciekawy plot .
  Dziwna jest jego dolna galaz , ale ona najwyrazniej ( objerzalem fit-a ) bierze sie z duzych gwiazd ktore sa wpisane kilka
  razy do bazy !!! To jest blad szybkiej fotometri , ktorej uzywalem do wykonania tego plotu , pytanie dlaczego tak sie zrobilo !!!
  Pokazal to select :
       select star,ccdx,ccdy,ra,dec,magnitude,error from measurements where error*36<200 and magnitude<6;
    star  |  ccdx   |  ccdy   |     ra      |     dec     | magnitude | error
  -------+---------+---------+-------------+-------------+-----------+--------
      91 | 1036.91 | 47.3238 | 23.31604839 | -9.60967927 |   5.31844 | 3.1087
     467 | 1063.23 | 91.0504 | 23.29846282 | -9.18170552 |   4.77721 | 2.7739
     548 | 1113.33 | 100.588 | 23.26495287 | -9.08769136 |   4.48721 | 3.3597
    1634 | 1089.67 |  238.97 | 23.28085791 | -7.72738939 |   4.76314 | 2.9883
    1671 | 1381.26 | 241.175 | 23.08607928 | -7.69390064 |   5.65933 | 3.2605
    2660 | 1447.05 | 354.033 | 23.04234095 | -6.57554059 |   5.96585 | 3.0546
    3085 | 1152.95 | 408.602 |  23.2386629 | -6.04819064 |   4.62456 | 3.6024
    3984 | 1026.31 | 501.701 | 23.32336675 | -5.12412935 |   5.64929 | 3.7647
    5413 | 1122.12 | 664.935 | 23.25943525 |  -3.4960268 |    5.7647 | 4.9368
    6111 | 311.175 | 736.648 | 23.79908727 | -2.76214843 |   5.51003 |  3.357
    8884 | 137.922 | 1025.55 | 23.91304905 |   0.1090013 |   5.17253 | 2.5497
    9720 | 1525.09 | 1110.32 |  22.9909497 |  0.96273125 |   5.44885 | 3.9665
   10012 | 839.092 | 1139.58 | 23.44883384 |   1.2553205 |   5.47054 | 2.5823
   10548 | 460.616 | 1192.85 |  23.7007306 |  1.77965815 |   4.88338 | 3.6189
   10884 | 602.732 | 1224.68 | 23.60643263 |  2.10186238 |   5.87114 | 3.7453
   10900 | 1294.25 |  1226.5 | 23.14476955 |  2.12728019 |   5.53561 | 3.5055
   11660 | 210.183 | 1309.87 | 23.86606158 |  2.93068714 |   5.32216 | 2.8763
   11733 | 1493.56 | 1315.36 | 23.01195119 |  3.01142592 |   5.75485 | 2.5873
   11984 | 1082.82 | 1341.78 | 23.28609584 |   3.2818336 |   4.62246 | 5.2837
  
• Jesli wziasc pod uwage tylko ta gorna galaz plotu to mamy zaleznosc, nie rozumiem tylko dlaczego nasza szybka fotometria dziala
  na tyle zle ze pojawia sie ta dolna galaz ! Jak patrzylem na ta gwiazde (1036,47) to dla niej radius powinien wyjsc jakies
  10 a jest 3 i dlaczego tak jest ???
  Jakis bez sensu maly cluster wychodzi dla tej gwiazdy see : /pi20/msok/badania/nova_flares/20060911/.gdbinit !!!
  ma tylko 23 pixle !!!
  Juz rozumiem, cluster jest szukany po laplacjanie .... Hmm to chyba nie jest najlepsze rozwiazanie, pojawiaja sie <0 wartosci !!!
• jest juz niezle tu mam ladny wynik :
      heplx43:/pi20/msok/badania/nova_flares/20060911/radius_test2
   trzeba go wrzucic na strone i fitnac prosta i uzyc zaleznosci !!!
  plot zaleznosci Star_radius vs Magnitude ( more ) dla clustra znajdywanego po
  raw data z Tresh=500 ADU ( na badanej klatce na sztywno )
      
  
  
• Okazuje sie ze alogorytm fast photometry z clustrem liczonym po raw data
  daje za duzo pixli w clustrze przez to pozycja jest slabo okresona
  i astrometria slabo zbiega, jak dalem recznie wysoki prog na cluster = 500 ADU
  to zbiegala znacznie lepiej, moze trzeba dac >5 sigma zeby bylo dobrze
  
• I fitted straight line and obtained result :
  star_radius = - 107.26 * MAGNITUDE + 967.51
  fit result is here. in order to have upper border of star_radius vs magnitudo , I shifted
  fitted line 50 arcsec up and obtained the following result
  Where formula is :
  star_radius = - 107.26 * MAGNITUDE + 1050
  Takiej formuly uzylem w skrypcie do_flareevents.pl zeby odrzucac przypadki gdy w poblizu jest duza gwiazda, promien 'bliskosci' jest wyliaczny z tego wzor
  
  
• sprawdzenie po jakiej liczbie obserwacji pola mozna juz wykrywac nowe obiekty :
  Zrobilem rozklad liczby nowych w funkcji liczby obserwacji pola z tego wynika ze min_obs_field=50 jest ok, choc wlasciwie powinno byc kolo 80
  
  
• ze wzgledu na uzycie tycho do odrzucania gwiazd staly zrobilem rozklad jasnosci gwiazd w tycho . Ustalilem ze tylko gwiazdy jasniejsze niz 11mag moga byc
  odrzucone jesli sa w tycho, zeby nie bylo ze nie wykryjemy jakies flary ponizej naszego zasiegu !