Obliczenia do PHD

0. Sciaganie danych SWIFT-a do DB :

fill_swift_grb! - loads grb/alert and grb_lc tables
get_swift_grb!

get_gcn_all.sh

1. Obliczenia czasu obserwacji znormalizowanego do 4 PI

- akceptance : select calc_obs_time_total(20061027,20061027,1000);
   FOV=33x33 - 2004/2005 = 0.121839 RAD -> 0.121839 / 4PI = 0.009695987
   FOV=20x20 - 2006
pusicic missing astro - skatalogowac ponownie 2004_2005 - zeby miec pelna akceptancje
na razie moge to przyblizyc srednia na dzien i zobaczyc
   piastrometry k2a_040819_034.mag k2a_040819_034.ast a.cfg -ord 4 -pixscale 59.5 -cat $DATADIR/cat/act -err 0.3 -errf 0.6 -verb -fi 180
klatka posumowana po 20 zajmuje 271 sec - 2*20 sec=40sec czas martwy -> 231 sec czysty obrazek
   select (select min(ttime_ut) from frame f2 where f2.idaynight=f1.idaynight and
   f2.ttime_ut>f1.ttime_ut)-ttime_ut from frame f1;
fix_old_astro! - dorabia brakujaca astro na danych z 2004-2005
   [pi@heplx47 ASAS_PIPELINE]$ /disk02/results/ASAS_PIPELINE/fix_old_astro.out

W tej chwili do obliczania czasu obserwacji uzylem baz pidb@heplx40 i 2004_2005@heplx48 , maja one pole stars i mozna po nim odroznic czy byly chmury czy nie
Wywolania mozna wykonac w ten sposob :

    -- PARAMETERS: start_night end_night min_star_count exptime
    select calc_obs_time_total(20061027,20061027,1000,10);
    OUTPUT : in days

   zeby policzyc znormalizowane do 4PI :
    -- PARAMETERS :
   -- 1 start_night - starting night
-- 2 end_night   - end night
   -- 3 min_star_count - minimal number of stars on image to count frame
   -- 4 exposition time - real exposition time ( without dead time - 240sec for sum 20 )
-- 5 fov [in degrees]
CREATE OR REPLACE FUNCTION calcobstime_norm( INTEGER, INTEGER, INTEGER, INTEGER, DOUBLE PRECISION ) RETURNS double p

select calcobstime_norm(20061027,20061027,1000,10,21);

Z tych obliczen wyszlo mi ze nasz czas obserwacji 4PI nieba wynoski okolo 2.4 dnia , co daje 4-6 GRB w polu widzenia

2. Stronka z grb :

3. Badania i flagowanie wartosci GRB

CRITERIA :
   SURE               - points with mag <= 13 mag
   OPTIMISTIC   - points with mag <= 15 mag
   QUASI-OPTIMISTIC - extrapolation <= 15 mag
poprawienie ladowania grblc
wyczyscilem grb/alert na pidb@heplx40 z swift,hete i integrala
oznaczenia wedlug :
   grb_swift_verylate - glc_time>1 hour
   grb_swift_early - glc_time<60 sec
   grb_swift_medium - 60 < glc_time < 15min
   grb_swift_late - 15min < glc_time < 1hour
ladowanie krzywych blasku, uwagi :
1. zrobilem zmienianie :
  - '--' -> '-'
  r' -> r_prim
  i' -> i_prim
  k' -> k_prim
proba automatyzacji wyceniania blyskow i krzywych wedlug mozliwosci obserwacyjnych pi :

   OPTIMISTIC : update grb set ot_flag=0 where ot_flag is null and ot_time<120 and min_mag<15;
   PESIMISTIC : update grb set ot_flag=2 where (ot_flag is null or ot_flag=-1) and min_mag>19;
     PESIMISTIC : update grb set ot_flag=2 where (ot_flag is null or ot_flag=-1) and mag_first_limit>=18 AND has_ot='f' AND first_limit_time<=(120.00/(24.00*3600.0));
DODAC : dodac auto-fixing pola grb.name i wyliczanie ot_flag - tak jak mam powyzej te update

4. Wyniki szacunkowe ze SWIFT-a

Zroione sql-em :
select date_part('year',date),count(*) from alert a1 where seq_num=
=(select max(seq_num) from alert a2 where a2.trg_num=a1.trg_num)
and alert_status=2 and source_id=5 group by date_part('year',date) order by date_part('year',date);

YEAR	GRB#
2005	3
2005	125
2006	119
2007	14 ( on 20070310 )

                   SURE               - points with mag <= 13 mag
               OPTIMISTIC   - points with mag <= 15 mag
               QUASI-OPTIMISTIC - extrapolation <= 15 mag

Przyjmujac kryteria jak poprzednio : SURE , OPT , QUASI-OPT dla blyskow ktore maja wczesny odpowiednik optyczny mamy :

QUALITY	GRB# / year
SURE	3
OPTIMISTIC	20
QUASI-OPTIMISTIC	15

Z tego wynika ze z blyskow obserwowanych przez SWIFTa gdbysmy patrzyli przez caly czas w jego FOV to mamy rocznie :

QUALITY	GRB OC# / year
SURE	1-2
OPTIMISTIC	10
QUASI-OPTIMISTIC	7.5

Czemu QUASI-OPT jest mniej niz OPT ?
Co oznacza ze OT z potwierdzeniem ze SWIFTA mamy ~ 1/4 ( musze to jeszcze sprawdzic ) czyli 0.25-1 / 2.5 / 1.875 , natomiast takich blyskow ze powinnismy je widziec bedzie polowa
tzn 0.5-1 / 5 / 3.75 , moze moga byc czasowo potwierdzone przez jakiegos inneg satelite ( konus, integral etc ... )

4. Badania i ulepszenia algorytmow do flar

FlareStudies

Pokazywanie plotow dla gwiazdy - wraz ze znalezionym zakresem wiekszosci pomiarow wybranym przez flare-findera

jest sporo plotow tutaj /pi20/msok/badania/flares/200701/PLOT/
plot_star.sh

5. Cone distance calculation

- dla przypadkow policzyc cone distance :
   /home/8/msok/doc/PHD/Results/online_flashes
   satcone 21h44m40s "-41o49'" 20050807_235916 -ra_hms -dec_page
skrypcik w perlu ( parse_pi_ot.pl )
    parse_pi_ot.pl single_tab2004_2005.txt

6. Wyznaczanie efektywnosci algorytmow on-line

1. Confirmation on next ( ccdsingle )

wycinki sampli :

cd /disk03/results/msok/simul/ccdsingle/Sample/
get_samples.sh /disk03/results/fits/20060527/

to powoduje zrobionie sampli , trzeba je obejrzec skryptem viewlist! czy maja jakis sens. Jest tam robione le_samp_filter ktore zapisuje do pliku list tylko
dobre sample bez zadnych smieci czy innych problemow
Nastepnie trzeba przygotowac pliki konfiguracyjne do zapuszczenia symulacji. Parametry symulacji znajduja sie w pliku mc.cfg i tam powinny byc zmieniane
Nalezy plik mc.cfg zaincludowac na poczatku ccd.cfg , piszac :
    %LOAD% mc.cfg
Nastepnie przeedytowac plik mc.cfg ustawic :
   # ustawic ze uzywa plikow ( to juz jest w offline.cfg , ale dla pewnosci ) :
CCD_MONTE_CARLO=1

# on how many images sample shuold be pasted ( for ccdsingle put 1, for ccddouble put 0 )
CCD_PUT_OBJ_ON_N_FRAMES=1
             # if samples coordinates on second camera should be calculated from RADEC ( only for ccddouble )
             CCD_PUT_SAMPLE_BY_RA_DEC=1

            # path to directory where samples are stored :
        CCD_IMAGES_SAMPLE_DIR=Samples/


           # radius of matching sample with found event, sample is considered as
           # identified if there exists found final event closer then this value :
           CCD_MC_GEN_IDENT_REDIAL=10

            # special option to check only generated events :
          CCD_CHECK_GEN_ONLY=0

    # sample putting options :
      # in case >0 specifies |r_put-r_sample|<D - distance in which sample can be used
          # - this is to use sample from same part of image :
          # -1 disables this check and sample from any position can be put
          CCD_PUT_SAMPLE_CUT_OUT_IN_WIDTH=-1

          # this option is for ccddouble, in order to put image of same star on second
          # camera
          CCD_PUT_SECOND_SAMPLE_BY_NAME=0

          # if logining of sample position and frame is performed
          CCD_LOG_SAMPLE_PUT=0

            # in order to use samples from frames close to current frame on which sample
          # is pasted enable this option :
          CCD_PUT_TAKEN_FROM_FRAME_IN_RANGE=0

           # number of samples to be put on single image :
       CCD_N_SAMPLES_TO_PUT_ON_FRAME=1
przygotowac plik mcinput.txt o nastepujacej postaci :

# Magnitudes to test
9.00,10.00

# paramters to test :
CCD_NEW_LAPLACE_TRESHOLD_IN_SIGMA 4.00 8.00 1.00 DOUBLE
inne dane :
20060519 :
get_samples.sh /disk03/results/fits/20060519/

tabelka z wynikami :

Night	Test location	Results ( mag=9 )	INFO
20060519	heplx48:/disk03/results/msok/simul/ccdsingle/20060519	-	NO PIXEL SHIFT ( taken on old pi2 ) ! Verfy small number of images Camera focusing, clouds + tests - cannot be used for eff/bkg tests
20060520	heplx48:/disk03/results/msok/simul/ccdsingle/20060520	eff_bkg_9.00.txt laplace=4 laplace=12	NO PIXEL SHIFT ( taken on old pi2 ) ! Very cloudy data only frames upto 00350 are acceptable for testing.
20060521	heplx48:/disk03/results/msok/simul/ccdsingle/20060521	???	NO PIXEL SHIFT ( taken on old pi2 ) ! very poor data
20060522		???	Data with pixel shift ( taken on DRAGON ) Almost no images - NO TESTS
20060523		???	Data with pixel shift ( taken on DRAGON ) Very cloudy data - NO TESTS
20060524		???	Data with pixel shift ( taken on DRAGON ) Very cloudy data - NO TESTS
20060525		???	Data with pixel shift ( taken on DRAGON ) quite cloudy - only ~500 images
20060526	heplx48:/disk03/results/msok/simul/ccdsingle/20060526/RESULTS_bigtest_CCD_CHECK_PREV_OF_MAX_PIXEL=0 heplx43:/pi20/msok/praca/MY/badania/eff_vs_bkg/20060526/no_tv_on_max_pixel_in_cluster	eff_bkg_9.00.txt laplace=4 laplace=12	Data with pixel shift ( taken on DRAGON ) quite good data CCD_CHECK_PREV_OF_MAX_PIXEL=0 Very high background ( probably due to pixel shift )
20060526	heplx48:/disk03/results/msok/simul/ccdsingle/20060526/RESULTS_bigtest_CCD_CHECK_PREV_OF_MAX_PIXEL=1 heplx43:/pi20/msok/praca/MY/badania/eff_vs_bkg/20060526/with_tv_on_max_pixel_in_cluster	eff_bkg_9.00.txt laplace=4 laplace=12	same , but CCD_CHECK_PREV_OF_MAX_PIXEL=1 Which is Tv cut applied to MaxPixel in cluster Widac ze jest kilku procentowy spadek efektywnosci Natomiast spadek liczby przypadkow tla nie jest za duzy wiec chyba mozna tego ciecia wogole nie uzywac ! Natomiast dla laplace=12 to spadku efektywnosci praktycznie nie ma !
20060527	heplx48:/disk03/results/msok/simul/ccdsingle/20060527/RESULTS_bigtest	eff_bkg_9.00.txt	Data with pixel shift ( taken on DRAGON )

20070409		???
20070425	heplx48:/disk03/results/msok/simul/ccdsingle/20070425/RESULTS_no_tv_on_maxpixel	eff_bkg_9.00.txt laplace=4 laplace=12 (plots done by plot_eff_bkg! )	Done on data 20070425 - images taken with shutter in normal mode This test was without checking of Tv on max pixel - there is no so much background so it can be used ( CCD_CHECK_PREV_OF_MAX_PIXEL=0 )
20070425 NOTE : THIS PLOTS WERE TAKEN TO PHD	heplx48:/disk03/results/msok/simul/ccdsingle/20070425/RESULTS_with_tv_on_maxpixel	MAG=9 : eff_bkg_9.00.txt laplace=4 laplace=12 MAG=10 : eff_bkg_10.00.txt laplace=4 laplace=12 MAG=11 eff_bkg_11.00.txt laplace=4 laplace=12 AFTER getDist change : Mag=9 eff_bkg_9.00.txt laplace=4	Same data - 20070425 - images taken with shutter in normal mode. But Tv was checked on maxpixel ( CCD_CHECK_PREV_OF_MAX_PIXEL=1 )
20070425	heplx48:/disk03/results/msok/simul/ccdsingle/20070425_once/RESULTS/tn6_tv2_bigstar5/	Gen# = 622 (622) GenIdent# = 440 Eff = 0.71	test with normal settings from ccd cfg where CCD_BIG_STAR_MAX_MAGNITUDO=5
20070425	heplx48:/disk03/results/msok/simul/ccdsingle/20070425_once/RESULTS/tn6_tv2_bigstar6	Gen# = 622 (622) GenIdent# = 437 Eff = 0.70	CCD_BIG_STAR_MAX_MAGNITUDO=6 Only 3 generated events were lost due to more strict cut

Script for making Bkg vs Eff plot is : eff_vs_bkg_multi.C , example of usage :
.x eff_vs_bkg_multi.C("eff_bkg_9.00_lap4.txt",10,20000,1)
also script : plot_eff_bkg!
See : /pi20/msok/praca/MY/badania/eff_vs_bkg

2. Coincidence ( ccddouble )

Coincidence of 2 cameras on single parallactic mount was tested in LCO, testy w heplx48:/disk03/results/msok/simul/ccddouble/20060526/FINAL_RESULTS/ , results are in table below :

Night	Test location	Results ( mag=9 )	INFO
20060526	heplx48:/disk03/results/msok/simul/ccddouble/20060526/FINAL_RESULTS/	MAG=9 eff_bkg_9.00.txt laplace=4 laplace=12	Standard set of paramters only Tn and Tv tested
20070526 (details )	heplx48:/disk03/results/msok/simul/ccddouble/20060526/FINAL_RESULTS_DETAIL	MAG=9 eff_bkg_9.00.txt laplace=4 laplace=12	Detail test of smaller paramter space

UWAGA :
Aktualnie skrypt synch_eff_bkg.sh sciaga wyniki testow efektywnosci z pi2@pidb na heplx40@pidb oraz updatuje tabelke NightStat o pola z testami Eff. Niestety nie jest to miarodajne , bo
wklejanie sampla na chmury tak naprawde nie jest madre i zwykle jest on wykrywany !!!!
W tej sytuacji trzeba inaczej szacowac efektywnosc, poprzez sprawdzanie efektywnosci wykrywania gwiazd na klatce i nastepnie wyznaczenie efektywnosci kolejnych ciec. Lub oszacowanie
ile efektywnosci jest tracone na kolejne ciecia

7. Wyznaczenie efektywnosci redukcji i katalogowania

To zostalo zrobione przy uzyciu katalogu TYCHO , zostaly do tego uzyte gwiazdy z katalogu TYCHO-2 i sprawdzilem ile z tych gwiazd my widzimy , rozne wariacje :

zbiorcza efektywnosc na przdzialach jasnosci 0-12 i 0-14
eff vs mag ( dla 1 , 10 , 20 , 30 i 48 klatek )
eff vs num images

Po pierwsze trzeba zaladowac dane do bazy zasianej przy pomocy katalogu TYCHO :

zalozyc baze ( np. photo_eff_test@heplx48 )
init_pidb_tycho!
run_cataloging_local! 20070512 - - photo_eff_test

Do celu tych testow napisalem odpowiednie skrypty i ich uzylem sa to :

    init_pidb_tycho!
    run_cataloging_local!
    check_frame_eff!
    check_eff_vs_image_num!
    check_photo_eff_vs_mag!
    check_photo_eff_vs_mag_all!
    check_eff_vs_xy!
    plot_eff_vs_mag!
    plot_eff_vs_xy!
    plot_trash_vs_mag!

SQL Procedures :
CalcPhotoEffVsMag
CalcPhotoEffOnFrame
CalcPhotoEffVsXY
ad2xy_x / ad2xy_y ( C / pgsql procedure )
ad2xy_frame_x / ad2xy_frame_y ( C / pgsql procedure )

Uzywaja one funkcji w pg/sql : CalcPhotoEffVsMag oraz CalcPhotoEffOnFrame ( do tego celu rowniesz powstala biblioteka libastsql.so oraz funkcje : ad2xy_x / ad2xy_frame_x )
Przy okazji przy katalogowaniu zmienilem warunek w selekcie gwiazd z bazy na coalesce(no_measurements,0)>=0 , oczywiscie mozna to zmienic paramterem , ale default jest >=0
bo jak mam gwiazdy z TYCHO to one na razie maja no_measurements=NULL.
Uzycie funkcji SQL :
             select   CalcPhotoEffOnFrame( 9, 56, 0.0 , 14.00 );
             select CalcPhotoEffVsMag( 9, 56, 0.0 , 14.00, 0.5 );

BADANIA glownie w heplx48:/disk03/results/msok/piphoto_efficiency na razie baza photo_eff_test ( 20070512 - aver20 )
Noce :

20070512 ( otwarta migawka ) , baza photo_eff_test , zakres klatek z h>=15 deg id_frm>=9 AND id_frm<=39 , field=0800+20 ( star# ~22000-25000) :
Po pierwsze zbiorcza efektywnosc na jednej klatce id_frm=12 ( k2a_070512_011.ast ) dla :
         eff( mag=0 - mag=11 ) = 0.875815085158151
         eff( mag=0 - mag=12 ) = 0.795994065281899
Praktycznie juz bez zmian ( bo malo tam jest gwiazd wiec nie maja wplywu ) :
         eff( mag=0 - mag=13 ) = 0.742578869362693
         eff( mag=0 - mag=14 ) = 0.740797028031071
         eff( mag=0 - mag=15 ) = 0.740765756258179

Nastepnie robie ploty "Eff vs mag" dla roznej liczby klatek ( skrypty "nohup check_photo_eff_vs_mag! 9 56 > out 2>&1 &" , plot_eff_vs_mag! )

      Obecnie ploty sa tutaj heplx48:/disk03/results/msok/piphoto_efficiency/ ( oraz heplx43:/pi20/msok/praca/MY/badania/photo_eff/eff_vs_mag/eff_vs_mag/ )

         All plots

         Efficiency plots :

             Efficiency vs magnitudo
           Efficiency vs number of images
             Efficiency vs XY ( 200x200 pixels )

         Objects not having counterpart in TYCHO :

             Bkg vs magnitudo
             Bkg vs number of images (~3500 / per klatka )
   Bkg vs xy ( single image )
             Number of new objects vs image no

     TODO :

Sprawdzic co to sa te smieci ktore sie lapia - byc moze wiekszosc to : galaktyki - galaxy catalog needed, opened shutter ( check night with normal shutter ) etc ...
Tresh vs image no - zrobic fit od pewnego momentu i zobaczyc ile jest Trash(image_no=0) - to powinny bys obiekty spoza TYCHO na tym polu

20070425 ( normal shutter ) , zakres klatek z h>=15 deg id_frm>=9 AND id_frm<=47

Instrukcja ladowania :
1. cd /disk03/results/msok/piphoto_efficiency/20070425_normal_shutter
2. nohup init_pidb_tycho! - photo_eff_test_20070425 > init_pidb_tycho.log 2>&1 &
3. Wybor klatek :
  photo_eff_test_20070425=> select id_frm,sobject,spathtofile,shutter_mode from frame order by id_frm;
  id_frm | sobject |    spathtofile     | shutter_mode
  --------+-----------+--------------------+--------------
        1 | S0720+20 | k2a_070425_000.ast |            1
        2 | S0720+20 | k2a_070425_001.ast |            1
        3 | S0720+20 | k2a_070425_002.ast |            1
        4 | S0720+20 | k2a_070425_003.ast |            1
        5 | S0720+20 | k2a_070425_004.ast |            1
        6 | S0720+20 | k2a_070425_005.ast |            1
        7 | S0720+20 | k2a_070425_006.ast |            1
        8 | S0800+20 | k2a_070425_007.ast |            1
        9 | S0800+20 | k2a_070425_008.ast |            1
       10 | S0800+20 | k2a_070425_009.ast |            2
       11 | S0800+20 | k2a_070425_010.ast |            2
       12 | S0800+20 | k2a_070425_011.ast |            2
       13 | S0800+20 | k2a_070425_012.ast |            2
       14 | S0800+20 | k2a_070425_013.ast |            2 | 1177545070
       15 | S0800+20 | k2a_070425_014.ast |            2 | 1177545341 <== tez schodek
       16 | S0800+20 | k2a_070425_015.ast |            2
       17 | S0800+20 | k2a_070425_016.ast |            2
       18 | S0800+20 | k2a_070425_017.ast |            2
       19 | S0800+20 | k2a_070425_018.ast |            2
       20 | S0800+20 | k2a_070425_019.ast |            2
       21 | S0800+20 | k2a_070425_020.ast |            2
       22 | S0800+20 | k2a_070425_021.ast |            2
       23 | S0800+20 | k2a_070425_022.ast |            2
       24 | S0800+20 | k2a_070425_023.ast |            2
       25 | S0800+20 | k2a_070425_024.ast |            2
       26 | S0800+20 | k2a_070425_025.ast |            2
       27 | S0800+20 | k2a_070425_026.ast |            2
       28 | S0800+20 | k2a_070425_027.ast |            2
       29 | S0800+20 | k2a_070425_028.ast |            2
       30 | S0800+20 | k2a_070425_029.ast |            2
       31 | S0800+20 | k2a_070425_030.ast |            2
       32 | S0800+20 | k2a_070425_031.ast |            2
       33 | S0800+20 | k2a_070425_032.ast |            2 | 1177550243
       34 | S0800+20 | k2a_070425_033.ast |            2 | 1177550533 <==== tez schodek w trash plocie, PLANE OR SATELLITE !!!
       35 | S0800+20 | k2a_070425_034.ast |            2
       36 | S0800+20 | k2a_070425_035.ast |            2
       37 | S0800+20 | k2a_070425_036.ast |            2
       38 | S0800+20 | k2a_070425_037.ast |            2
       39 | S0800+20 | k2a_070425_038.ast |            2   | 1177551894
       40 | S0800+20 | k2a_070425_039.ast |            2   | 1177552166   <==== schodek na liczbie "nowych" , PLANE or SATELLITE again
       41 | S0800+20 | k2a_070425_040.ast |            2
       42 | S0800+20 | k2a_070425_041.ast |            2
       43 | S0800+20 | k2a_070425_042.ast |            2
       44 | S0800+20 | k2a_070425_043.ast |            2
       45 | S0800+20 | k2a_070425_044.ast |            2
       46 | S0800+20 | k2a_070425_045.ast |            2
       47 | S0800+20 | k2a_070425_046.ast |            2
       48 | S0800+20 | k2a_070425_047.ast |            2
       49 | S0800+20 | k2a_070425_048.ast |            2
       50 | S0800+20 | k2a_070425_049.ast |            2
       51 | S0735-70 | k2a_070425_050.ast |            1
       52 | S0735-70 | k2a_070425_051.ast |            1
       53 | S0735-70 | k2a_070425_052.ast |            1
       54 | S0735-70 | k2a_070425_053.ast |            1
       55 | S0735-70 | k2a_070425_054.ast |            1
       56 | S1242+30 | k2a_070425_055.ast |            1
       57 | S1242+30 | k2a_070425_056.ast |            1
       58 | S1242+30 | k2a_070425_057.ast |            1
       59 | S1242+30 | k2a_070425_058.ast |            1
       60 | S1242+30 | k2a_070425_059.ast |            1
       61 | S1242+30 | k2a_070425_060.ast |            1
       62 | S1242+30 | k2a_070425_061.ast |            1
       63 | S1739-30 | k2a_070425_062.ast |            1
       64 | S1739-30 | k2a_070425_063.ast |            1
       65 | S1520-20 | k2a_070425_064.ast |            1
       66 | S1520-20 | k2a_070425_065.ast |            1
       67 | S1520-20 | k2a_070425_066.ast |            1
       68 | S1520-20 | k2a_070425_067.ast |            1
       69 | S1520-20 | k2a_070425_068.ast |            1
       70 | S1520-20 | k2a_070425_069.ast |            1
       71 | S1520-20 | k2a_070425_070.ast |            1
       72 | S1840-20 | k2a_070425_071.ast |            1
       73 | S1840-20 | k2a_070425_072.ast |            1
       74 | S1840-20 | k2a_070425_073.ast |            1
       75 | S1242+30 | k2a_070425_080.ast |            1
       76 | S1242+30 | k2a_070425_081.ast |            1
       77 | S1242+30 | k2a_070425_082.ast |            1
       78 | S1242+30 | k2a_070425_083.ast |            1
       79 | S1242+30 | k2a_070425_084.ast |            1
       80 | S1242+30 | k2a_070425_085.ast |            1
       81 | S1242+30 | k2a_070425_086.ast |            1
       82 | S1715-40 | k2a_070425_087.ast |            1
       83 | S1715-40 | k2a_070425_088.ast |            1
       84 | GRB070426 | k2a_070425_089.ast |            2
       85 | GRB070426 | k2a_070425_090.ast |            2
       86 | GRB070426 | k2a_070425_091.ast |            2
       87 | GRB070426 | k2a_070425_092.ast |            2
       88 | GRB070426 | k2a_070425_093.ast |            2
       89 | GRB070426 | k2a_070425_094.ast |            2
       90 | S1715-40 | k2a_070425_095.ast |            2
  (90 rows)
  
  photo_eff_test_20070425=>
4. Czyli mamy ciag klatek z pola S0800+20 klateki od id_frm=10 do 50 ( czyli 41 klatek )
5. Puszczam :
  
  cd /disk03/results/msok/piphoto_efficiency/results/20070425/
  nohup ./photo_eff_20070425! > out 2>&1 &
6. I zrobione :
Total efficiency on ( wartosci podobne jak na 20070512 )
             Single frame ( id_frm=10 )
             Single frame ( id_frm=20 )
             Single frame ( id_frm=30)
             All 51 frames : 51 frames in 0-12 mag range , 51 frames in 0-15 mag range

All Plots

          Efficiency plots :

             Efficiency vs magnitudo ( all )
           Efficiency vs number of images
             Efficiency vs XY ( 200x200 pixels )

         Objects not having counterpart in TYCHO :

             Bkg vs magnitudo ( all frames , single image )
             Bkg vs number of images (~3500 / per klatka )
   Bkg vs xy ( single image )
             Bkg vs number of images , new vs number of images    - steps ( see list of images ) are due to SATELLITES or PLANES - I checked on fits - it is nice, in other parts the steep of line
                                                                                                       is constant
20070604 - shutter opened, but different field (0851-70) , star# ~30000-35000
During this night data was collected with permanently opened shutter, but on different field and it can be seen on plot showing number of background events vs imageno ! - This is due to fact that there is much more stars and also more stars which can produce the "open shutter" strip.
All Plots

          Efficiency plots :

             Efficiency vs magnitudo ( all )
           Efficiency vs number of images
             Efficiency vs XY ( 200x200 pixels )

         Objects not having counterpart in TYCHO :

             Bkg vs magnitudo ( all frames , single image )
             Bkg vs number of images (~3500 / per klatka )
   Bkg vs xy ( single image )
             Bkg vs number of images , new vs number of images    - steps ( see list of images ) are due to SATELLITES or PLANES - I checked on fits - it is nice, in other parts the steep of line
                                                                                                       is constant

Jesli chodzi o nachylenie plotu, ktore jest wyraznie wieksze niz poprzednio to okazuje sie ze jest to wynik tego ze na tym polu ( 0851-70 ) jest znacznie wiecej gwiazd (~33000) i przez to tez wiecej paskow od otwartej migawki co powoduje, ze jest tez znacznie wiecej przypadkow "nowych" od tych paskow. Stad tez ~60 nowych/klatke na tych zdjeciach.
20070604 - after shutter opened correction

Dane zebrane z otwarta migawka poprawione przy pomocy skryptu fix_shutter! . Widac ze liczbna nowy obiektow / klatke zmniejszyla sie ~ x 2 .
All Plots

       Efficiency plots :

             Efficiency vs magnitudo ( all )
           Efficiency vs number of images
             Efficiency vs XY ( 200x200 pixels )

         Objects not having counterpart in TYCHO :

             Bkg vs magnitudo ( all frames , single image )
             Bkg vs number of images (~3500 / per klatka )
   Bkg vs xy ( single image )
             Bkg vs number of images , new vs number of images    - steps ( see list of images ) are due to SATELLITES or PLANES - I checked on fits - it is nice, in other parts the steep of line
                                                                                                       is constant
eff vs image no - 20070425,20070512,20070604,20070604_sc on single plot
Wiele nocy dla zadanego pola , select :

select sobject,count(*) from frame group by sobject having count(*)>100 order by count(*) desc;
    Dla aver20 :
        S0800+20 |   839
J0234-50 |   451
S0914-80 |   258
    S1000+10 |   205
O0840+20 |   181
S1122-70 |   175
R0720-20 |   174
S1200-10 |   162
S1200+30 |   162
Jednak posluzy mi do tego baza aver20_tycho , ktora wlasnie sie laduje tutaj : heplx48:/disk03/results/PI_PIPELINE/all/
All plots
Wynki na heplx48:/disk03/results/msok/piphoto_efficiency/results/aver20_tycho/
baza aver20_tycho@heplx48 - wielki test bkg_vs_image i eff_vs_image
ewenutlanie mozna jeszcze potestowac na tych nocach :
- Do testow : 20070509-43 , 20070511-59 , 20070522 -48, 20070529-77 , 20070604-57->OK !!
select sobject,count(*) from frame where idaynight=20070509 and falt>15 group by sobject;
POPRAWKI w funkcji : EffAndBkgForField

28799

16766

Eff vs IMAGE_STARS
Script plot_eff_vs_imagestars!
TEST heplx48:/disk03/results/msok/piphoto_efficiency/results/aver20_tycho/0800+20/eff_vs_stars
nohup ./plot_eff_vs_imagestars\! - - 9 56 > out2 2>&1 &
RESULTS :
Small test for 20070512 ( photo_eff_test database ) - NOT CONCLUSIVE

SPRAWDZENIE DLA POLA 0800+20 :
- Eff vs istarcount ?
DZIALA na aver20_tycho@heplx48 :
heplx48:/disk03/results/msok/piphoto_efficiency/results/aver20_tycho/0800+20/eff_vs_stars/NEW/ALL
Zrobic ploty i do pracy !!!

DLA SINGLE IMAGE :
DZIALA : /disk03/results/msok/piphoto_efficiency/results/single_tycho
Number of NEW/BKG vs image number :
- cd $SRCDIR/ccd/scripts/shell/
   cd /disk03/results/msok/piphoto_efficiency/results/aver20_tycho/0800+20a
   check_bkg_vs_image_num! 20050101 20100101 aver20_tycho "0800+20" "Number of new objects not
      present in TYCHO star catalog" "Number of objects not present in TYCHO star catalog" 0
   STARASZNIE WOLNO DZIALA !!! :
      /disk03/results/msok/piphoto_efficiency/results/aver20_tycho/0800+20
      Counting PI/TYCHO stars in this range
      Zajmuje najwiecej ...
   Nadal wolno dziala ...

   heplx48:/disk03/results/msok/piphoto_efficiency/results/aver20_tycho/0800+20/NEW/
   heplx43:/pi20/msok/praca/MY/badania/nova_eff/new_vs_image/NEW
   check_new_vs_image_many! check_new_vs_image!
   tutaj po prosut new_star vs image

ZAPISKI :

      nohup ccdsingle -putsample=9.00 >putsample.out 2>&1
          pieffbkg

           TESTY zrobic na nocach : 20060519 , 20060520, 20060521
        Bo na nich nie ma pixel shifta !!!
        Dziala : /disk03/results/msok/simul/ccdsingle/20060520/
        Moze dociagnac : 20070318/
        Moze wziasc jedna noc z otwarta migawka !

8. Krzywa sigma_mag VS mag

Krzywa dla probki gwiazd obserwowanych 20070606 jest tutaj. Zostalo to wykonane w heplx43:/pi20/msok/praca/MY/badania/sigmamag_vs_mag/20070606 przy pomocy skryptu : plot_sigmamag_vs_mag! . Przyjemna jest parametryzacja , ktora pozwoli generawac gwiazdy do wyznaczenia efektywnosci algorytmow do flar i nowych.

9. Results of subsequent cuts

Od 20070529 wlacznie mozna juz robic ladne i sensowne ploty

nova algo :
cd heplx48:/disk03/results/msok/plots/nova_cuts/
plot_nova_cuts!
A oto wyniki dla 5 nocy : 20070603,20070604,20070605,20070606,20070607
All plots

a tu zrobie wszystkie : heplx48:/disk03/results/msok/plots/nova_cuts/all/
flare algo :

All plots

10. On-line efficiency vs istarcount

Do tego celu zrobilem tabele EventGen w ktorej jest kazdy event wygenerowany z wartoscia czy zostal wykryty czy nie. Ale to nie dziala dla chmur, bo jak naklejamy blysk na chmury to nie jest to dobry test, bo on i tak zostanie wykryty !!!!
Mam pomysl : zeby to zrobic tak ze zrobie test eff na bezchmurnej nocy , potem zrobie fast photo i spradze jej efektywnosc per klatka ( czy po calosci ) i porownam te 2 liczby
generalnie to eff_algo < eff_fast photo o inne straty na innych cieciach, ale stosunek powinien byc staly , czyli potem z eff_fast photo mozna wyliczyc eff_online_algo ...

- tabela EventGen - dodac zaladowac genevents_0.log do niej i miec tam
eg_night i eg_frame i juz mozna query zrobic do eff vs istarcount !!!
zautomatyzowac i bedzie co trzeba ...
sprawdzic : /opt/pi/dev/pisys/daq/ndir/data/20070512/offline/
zrobic zebi pieffbkg generowal odpowiedniego sql-a
   heplx40:/lhome/piwww/www/pi0/daq/events/200706/20070618/offline
policzyc eff vs stars tam - 2 histogramy , potem podzielic biny jeden
przez drugi : divide_histo.C
calkiem ladny plot, choc efekt calkiem odwrotny, ale to chyba jasne ...
     http://grb.fuw.edu.pl/pi0/daq/events/200706/20070618/offline/eff_vs_stars.gif
wklejanie nie uwzglednia chmur, wkleja i jest OK zawsze, bo na chmury nakleja a nie pod ...
!!! Hmm no wlasnie nie jest wiec miarodajna !!!!!!
To jednyna szansa , ze test z bazy danych pokaze jak to zalzy od liczby
gwiazd tamto jest dosc miarodajne i da sie przeniesc na on-line
Mozna skatalogowac po Tn i bedzie wynik :)
   Skrypt : plot_eff_vs_stars!
Albo szybka fotometrie z takim samym progiem jak w on-line Tn=5 i
sprawdzic jak tam zalezy efektywnosc od istarcount

Oto moje zapiski, tabela EventGen sciaga sie sama na hoza .

11. Wyznaczenie efektywnosci algorytmu do nowych

dodalem funkcje xy2ad_a i xy2ad_d w pidb_ast.sql i fitslib/astsql/libastsql.so
zrobilem baze nova_simul i sciagnalem pare nocy z aver20 @ heplx48
przetestowalem ja na bazie nova_simul @ heplx43

W tej chwili robi sie to tak ze trzeba miec 2 bazy :

glowna baza aver20 ( w tej chwili uzywam aver20_tycho@heplx40 )
mala pusta baza do ktorej beda wrzucone wygenerowane sztuczne gwiazdy ( w tej chwili nova_simul@heplx43 )

Nastepnie jest gotowy skrypt, ktory mozna wykonac na jednej nocy :

gen_nova_and_anal.sh 20070512

lub inny skrypt ktory bierze liste nocy i wykonuje ten pierwszy :

gen_nova_and_anal_many.sh night_list

Wyniki mam na razie tutaj : heplx43:/pi20/msok/praca/MY/badania/nova_eff/ , skrypt doit_new!

I na razie puszczalem analize z wymaganiem co najmniej 5 obserwacji pola inaczej byla za mala efektywnosc bo taka jest ta baza aver20_tycho

12. Wyznaczenie efektywnosci algorytmu do flar

TESTY :

heplx43:/pi20/msok/praca/MY/badania/flare_eff

aver20_tycho@heplx48 :
select sobject,count(*) from frame where idaynight=20070616 group by sobject;
   sobject | count
   ----------+-------
I0240-10 |    19

aver20_tycho=> select id_frm from frame where idaynight=20070616 and sobject like '%1845-40';

aver20_tycho=> select star from measurements,stars where id_frm=3029 and star=id and no_measurements>=50;
star
---------
1626796
1626815
1626907
1626931

SCRIPTS :
   gen_flare.sh
   gen_flare_events.pl -night=20070616 -max_flare_count=1 -verb
   gen_flare.sh
   gen_flare_events.pl
   gen_flare_and_anal.sh

13. Straty efektywnosci na kolejnych cieciach :

Do tego zrobilem skrypty :
      show_cuts_effloss_coinc!
      show_cuts_effloss_confirm!
      plot_effloss!

    calc_eff_loss! - sporo polecialo na IfMore ??? Jak to czemu - sprawdzic i
      calc_eff_loss_slt_confirm!
      calc_eff_loss_slt_coinc!

Wyprobowalem jest w heplx43:/pi20/msok/praca/MY/badania/online_eff/EffLoss/ jest tez skrypt plot_effloss! ktory sciaga dane z pi2 i robi ploty

14. Efektywnosc wykrywania gwiazd VS istarcount oraz VS ratio ( image_stars/cat_stars)

Wymyslilem ze mozna sciagnac co 25 klatke ( wczesniej mialem co 100 ) i wyznaczyc srednia efektywnosc na pojedynczej klatce dla wielu nocy ( biorac do sredniej tylko klatki co maja >=1000 klatek ). Skrypt do sciagania jest zrobiony i dziala na heplx48 :
   nohup get_few_night_images_many.sh - - 25 > get_few_night_images_many_25.out 2>&1 &
Potem mozna skatalogwac i juz . Dane sa w heplx48:/disk03/results/fits , skrypt do katalogowania do bazy single_tycho_tresh5_more to : cat_single_tycho!
Moje notatki :
   - sciagnac co 100 klatke z Chile i sprawdzic jak sie zmienia efektywnosc na
        pojedynczej klatce w zal od liczby gwiazd na klatce.
   Ocenic eff kolejnych ciec i sprawdzic czy jest w miare stala.
   Potem wyznaczyc srednia efektywnosc po dlugim kawalku danych
   Puscic na heplx48 : get_few_night_images.sh

DZIALA : /disk03/results/fits/get_few_night_images_many.out
TESTOWAC : czy z crona dziala : /opt/pi/dev/pisys/log/get_few_night_images

WHERE DONE :
heplx48: /disk03/results/msok/piphoto_efficiency/results/single_tycho_tresh5_more/NEW
heplx43:/pi20/msok/praca/MY/badania/nova_eff/eff_vs_stars/single_tycho_tresh5_more

Wyznaczylem srednia efektywnosc wykrywania gwiazd na pojedynczej klatce ( z wielu nocy i co 25 klatka z nocy ) , eff_aver=0.489 ( poprzednio jak mialem co 100 klatke i
zle opcje katalogowania - tzn tylko pola dlugo obsewowane mi sie skatalogowaly to wyszlo ~0.54 ). Oto zapiski :

        - srednia eff z wielu klatek wyznaczyc, metoda single_tycho_tresh5 i juz !
       Zajac sie tym :
            /disk03/results/msok/piphoto_efficiency/results/single_tycho_tresh5/FIELD
        LADNY PLOT - juz mozna je do pracy wrzucic - to jest niezla paramteryzacja
        efektywnosci wykrywania gwiazd
           /pi20/msok/praca/MY/badania/nova_eff/eff_vs_stars/single_tycho_tresh5/eff_vs_fract.txt.eps
        Pytanie jak go rozszerzyc, tzn moze scany sciagnac ?
        I podkatalog FIELD tez piekne ploty !

            /disk03/results/msok/piphoto_efficiency/results/single_tycho_tresh5
             awk -v eff_avg=0 -v n=0 '{eff_avg=eff_avg+$2;n++;}END{print "Avarage
               Eff = "eff_avg/n;}' eff_vs_stars.txt
               Avarage Eff = 0.541194

        Sciaga sie jeszcze troche klatek beda w /disk03/results/fits
        mozna je bedzie skatalogowac do single_tycho_tresh5_more


        TESTY :

        Juz calkiem ladnie widac linie : eff vs fract :
            heplx43:/pi20/msok/praca/MY/badania/nova_eff/eff_vs_stars/single_tycho

       inne pole tu zbadam, oraz wiele roznych klatek z dowolnego pola :
       /disk03/results/msok/piphoto_efficiency/results/single_tycho_tresh5
       Mozna tez single ze scanow z wielu klatek wziasc i zrobic z nich Eff

       co 25 klatka :
            /disk03/results/msok/piphoto_efficiency/results/single_tycho_tresh5_more/

JESZCZE TROCHE NOTATEK :

- DZIALA : single image : eff vs istarcount :
   heplx48:/disk03/results/msok/piphoto_efficiency/results/single_tycho
   heplx48:/disk03/results/msok/piphoto_efficiency/results/single_tycho_tresh5

   pidb! :
      heplx43:/pi20/msok/praca/MY/badania/nova_eff/eff_vs_stars/other_fields
      select substr(sobject,2) as field,count(*) from frame where falt>40.00
      and idaynight>=20070604 group by field order by count(*) desc;
      1656-30 17
      1330-40 8
      select distinct idaynight from frame where sobject like '%1656-30' and
      falt>40.00 and idaynight>=20070604 order by idaynight;


      skrypt : prep_scp!


   Trzeba by chyba inna probke wziasc ? - na pewno tresh=5.00 trzeba puscic

15. Efektywnosc novaeff dla roznych MinObsField

Zostala zrobiona na heplx48 w /disk03/results/msok/nova/eff/nova_eff/ ( min_obs_field10/ , min_obs_field5/ min_obs_field20/ , min_obs_field30/ )
( heplx43: /pi20/msok/praca/MY/badania/nova_eff/new_vs_image/NEW/ALL/ )
Zapiski :

czesc 1 :

- generator novych :
   ./gen_nova_events.pl -night=20070614 -verb > out 2>&1
i juz mozna sprawdzic efektywnosc !!!
   /pi20/msok/praca/MY/badania/nova_eff
   NEW : heplx48:/disk03/results/msok/nova/eff/nova_eff gen_nova_and_anal.sh
         /disk03/results/msok/nova/eff/nova_eff/min_obs_field20
sprawdzic ze tutaj to nie zalezy od jasnosci
czesc 2 :

       - PHD , zrobic te /pi20/msok/praca/MY/badania/nova_eff/new_vs_image/NEW/ALL/
          na 1 plocie !!!
        Zrobic testy nova eff dla min_field_obs=20 i 30 , sprawdzic roznice
            jeszcze moze 10 przetestowac

16. Efektywnosc wykrywania flar

Aktualne testy na heplx48:/disk03/results/msok/simul/flares/TEST/20070616/I0240-10/min_points_above_3/min_points_above_1/
Glowny skrypt : gen_flare_and_anal.sh

Moje notatki :
       - KONTYNUOWAC :
        generator flar :
       /pi20/msok/praca/MY/badania/flare_eff/NEW/


        sprawdzic paramtery tej wykrytej flary i na poczatek z takimi sprawdzac
       sa tutaj heplx43:/pi20/msok/praca/MY/badania/flare_eff/FLARE_FIT/
       script : fit_flare.C
            20061127fast@heplx49 star=21878
            20061127ap0 - poprawna baza z ladna krzywa
        Dofitowac jakis model :
            /pi20/msok/praca/MY/badania/flare_eff/FLARE_FIT

        skrypt : gen_flare_and_anal.sh
        przetestowac najpierw : gen_flare.sh - sprawdzic paramtery
       ZROBIC TEST : /pi20/msok/praca/MY/badania/flare_eff/NEW/MANY/ANAL/

17. Przeglad przypadkow ( wyznaczanie jasnosci )

Lista watpliwych przypadkow - niektore na pewno trzeba usunac z listy ostatecznej przypadkow :

   2005.05.16 & 06:26:04 & 13h58m52s & -26\deg08' & 25026&63 & - & Konus 6:41:01
   2005.05.16 & 00:26:50 & 13h03m11s & -42\deg17' & 7888&47 & - & \\
   2005.05.02 & 08:36:34 & 13h53m22s & -11\deg23' & 59060&43 & - & \\
   2005.04.29 & 08:27:37 & 14h37m26s & -08\deg35' & 35516&11 & - & LEDA 1001633
   2005.02.15 & 06:51:10 & 11h06m27s & +19\deg57' & 15588&09 & - & \\
   2005.01.09 & 01:45:03 & 06h15m02s & +18\deg45' & 7730&47 & - & \\
   2004.12.13 & 05:23:02 & 05h26m40s & -01\deg16' & 9123&02 & - & \\

Wyznaczanie jasnosci :

       bad astro : CCD_BAD_ASTRO_DX CCD_BAD_ASTRO_DY
calc_event_mag.sh

przyklad :

       calc_event_mag.sh 20040912 0.24584 2.83025 00993 2 2 200 0

Stara metoda ( bardziej ogolna ) :

    - wyznaczanie mag_max starych przypadkow :
   heplx43:/opt/pi/dev/pisys/daq/ndir/data/FINAL_EVENTS/2005/Frame00606/Cam1
   do_parts_ast! list0 17.6792 14.3705 0 20050731_606_0_3_0_0_event

   do_event_lc! list0 17.6792 14.3705 0 20050731 606 0 2

18. Wyznaczenie limitu na OT dla blyskow z algorytmu wykrywania nowych na sumach po 20 ( skala czasowa ~ 12*20 = 240sec )

PHD - limit wynikajacy z algorytmu do nowych - z ciaglych obserwacji nieba
   wyliczyc ile dni obserwacji z dzialajacym algo do nowych mamy - i jaki ze tego wynika limit ?

   Tez trzeba by to wlaczyc !
   Moge to wyznaczyc sprawdzajac kiedy count(flareevents) where fl_type=1 > 0 i
   liczac wlasnie te noce :

     select idaynight,(select count(*) from flareevents where
   fl_night=idaynight and fl_type=1) as cnt from frame where (select count(*)
         from flareevents where fl_night=idaynight and fl_type=1)>0 group by
     idaynight order by idaynight;

         select idaynight,calcobstime_norm(idaynight,idaynight,1000,240,20,20)
         from frame where (select count(*) from flareevents where fl_night=idaynight
         and fl_type=1)>0 group by idaynight order by idaynight;

   wyznaczyelem w jakim zakresie nocy to dzialalo :

   SELECT : select count(*),min(idaynight),max(idaynight) from frame where (select count(*) from flareevents where fl_night=idaynight and fl_type=1)>0;
         WYNIK :
            count |   min    |   max
            -------+----------+----------
            32016 | 20060721 | 20070906
            (1 row)

   mozna zrobic sum(calcobstime_norm)
   DZIALA : [piwww@heplx60 ~]$ cat nova_algo_days.txt
   sumowac awkiem
      bash-2.05b$ sumcol\! nova_algo_days2.txt 2
      0.828143
   see heplx43:/home/8/msok/doc/PHD/Results/ResultsFigs/NOVA_ALGO

A drugi sposob , to taki ze po prosu zrobilem sum(calcobstime_norm) :

podobny limit, tu nie mamy kandydata na blysk wogole
   sprobowac w 2 SELECTACH sum(days) from (select calc .... )
   sprawdzic min_frames=5/10/2/1

     SELECT : select sum(obs_time) from (select
         calcobstime_norm(idaynight,idaynight,1000,240,20,10) as obs_time from frame
         where (select count(*) from flareevents where fl_night=idaynight and
         fl_type=1)>0 group by idaynight order by idaynight) as sel;
   WYNIK      = 0.829408472004836

   select sum(obs_time) from (select
         calcobstime_norm(idaynight,idaynight,1000,200,20,10) as obs_time from frame
         where (select count(*) from flareevents where fl_night=idaynight and
         fl_type=1)>0 group by idaynight order by idaynight) as sel;

                   =      0.693889581088023

19. Rozklad blyskow ktore maja OT we wspolrzednych galaktycznych

Nie widac tam zadnych zgrupowan, raczej jest tez ich za malo ( ~80 dla >=20050101 ) zeby mozna bylo cos definitywnie powiedziec
Badania robilem w heplx43:/pi20/msok/praca/MY/badania/grb_luminocity/gal_distrib/ oto ten rozklad
Co robilem :
    show_ot_in_gal.sh
    sel2txt! grb_with_ot.txt > g.txt
    awk '{print $2" "$3;}' g.txt > radec.txt
    java -jar /opt/pi/ext/dload/VOPlot/voplot1_4/binaries/VOPlot/voplot.jar

Zapiski :

     rozklad we wspolrzednych gal tych co maja widzialne OT - czy to nie wplyw
galaktycznej ekstynkcji ?
   Ciekawie to wyglada ale raczej o tyle ze wyglada jakby bylo wiecej
   przypadkow z OT w plaszczyznie gal ! czy to znaczy ze sa zle skojarzone
   z GRB jakies gwiazdy ???

20. Wyznaczenie jasnosci optycznej blyskow w zrodle

Przeliczylem D_lumin ( Omega_M=0.3, Omega_L=0.7 ) na podstawie zmierzonych z dla wszystkich mozliwych blyskow i wyszly takie ploty :

pierwsze 7200 sec ( distribution)

pierwsze 86400 sec

W tych plotach zostala uwzgledniona zmiana czasu na *(1+z) i ekstynkcja galaktyczna wedle stronki (calculator ekstynkcji ) Policzylem tez jak w funkcji z wyglada jasnosc w magnitudo blyskow na Ziemi dla roznych L od 10^35 do 10^40 erg/s/Hz , wyglada to tak.
Oba ploty w pliku ps.
Jak to sie robi podstawowe komendy :

grblc_anal -mag=12 -z=0.553 -filter=V
grblc_anal -filter=V,U,B,White,white,none -use_v_on_error
plot_grblumin.sh "-filter=V,U,B,White,white,none -use_v_on_error" 86400
mag_vs_z.C
load_extinction.pl - Skrypcik do sciagania wartosci ekstynkcji galaktycznej dla kazdego z blyskow

Przykladowe uzycie :
Zeby policzyc jasnosc w zrodle przy limicie 12 mag :
   grblc_anal -mag=12 -z=0.553 -filter=V
      Luminosity limit = 4.728709e+38[erg/s/Hz]

a przy limicie 14mag :
      grblc_anal -mag=14 -z=0.553 -filter=V
      Luminosity limit = 7.494499e+37[erg/s/Hz]

Potem juz sie wlasciwie samo wszystko liczy , skrypt w ROOT to mag_vs_z.C

Nie sa to moze zbyt madre wnioski , ale to co mi sie nasuwa to :

zasieg 14-15 mag pozwolilby nam zobaczyc znacznie wiecej blyskow , moze lepiej miec taki i troche mniejsze pole widzenia
latwiej zobaczyc blizsze blyski ( to oczywiste ) , czyli takie z mniejszym z , ale to oznacza ze sa one twardsze, SWIFT ma 15-350 keV i chyba przez to widzi raczej
dalsze blyski, dlatego moze bardziej obiecujacy i priorytetowy powinien byc GLAST. Wtedy bedzie wieksza szansa zobaczyc blizszy blysk ! :
1. SWIFT - 15-350 keV
2. BATSE - 20 keV do 1 MeV
3. GLAST - 10 MeV - 100 GeV !
990123 nie byl blisko ( z ~= 1.61 ) , ale byl bardzo jasny stad mial tak jasny OT

21 . UWAGA - ROOT ( narowy root-a )

Hatching problem z EPS :
Jest cos dziwnie z hatchingiem ( np histogramy - jak jest kilka to mozna tak zrobic ) potem ggv tego nie pokazuje w eps, trzeba zrobic :
w razie jak nie widac to kwestia ggv , trzeba wylaczyc anti-aliasing opis problemu : http://couet.home.cern.ch/couet/root/ht09.html
robi sie to Edit/Postscript Viewer Preferences/ i trzeba odhaczyc Anty-aliasing
UWAGA : po przelaczeniu na stale i jak jest dobrze w ggv to potem jest zle w ROOT !!!!!!!!!!!!!!!!!!!!!