ReWiA 2.0 html manual


HISTO          tworzenie histogramu

Zapis   HISTO destsht type { (val1 val2 step) | (srcsht width) }

Parametry

   destsht
      typu Sheet,  określa numer  widma, które będzie zawierać histogram widma
      bieżącego.  Jeżeli  widmo  o  numerze  sht  poprzednio  istniało  będzie
      skasowane. Dla type=0 widmo będzie miało jedno pole (YPTR=0), dla type=1
      dwa pola  (YPTR=0, EPTR=1).  Oba pola  typu REAL*4 (6). Dla type=2 widmo
      będzie miało 4 pola : dwa pierwsze typu jak srcsht, dwa ostatnie typu 6.
      Pola będą  zawierać  :  długość  fali,  współczynnik  korekcyjny,  ilość
      punktów w  widmie bieżącym zawarta w danym obszarze uśredniania energii,
      tj. w  zakresie od (długość fali - width / 2) do (długość fali + width /
      2)  oraz   średnia  energia.   Współczynnik  korekcyjny   liczony   jest
      następująco :

          absolutna energia (X)
wsp_kor(X) = -------------------------------------------
             śr. energia widma od X-width/2 do X+width/2

      Średnia energia to po prostu suma Y podzielona przez ich ilość.

   type
      typu Bool, określa rodzaj histogramu. Dla type=0 tworzony jest histogram
      współrzędnych Y,  dla  type=1  histogram  współrzędnych  X,  dla  type=2
      histogram współrzędnych  Y oraz  współczynniki korekcyjne  potrzebne  do
      absolutnej kalibracji  energii. W  przypadku widm dwuwymiarowych komenda
      HISTO oblicza  histogram ze współrzędnych Z. Skutkiem tego dla tych widm
      sensowny jest  jedynie histogram  typu 0.  Dla type=0 lub type=1 robiony
      jest histogram ze stałym krokiem i komenda pyta się o zakres i krok. Dla
      type=2 komenda  zapyta się o numer widma zawierającego absolutny rozkład
      energii oraz szerokość obszaru uśredniania.

   val1, val2   
      typu Real, określają zakres z którego będzie tworzony histogram
   step
      typu Real,  określa rozdzielczość  tworzonego histogramu.  Powinien  być
      większy od zera.

   srcsht
      typu Sheet, wymagany tylko dla type=2. Określa numer widma zawierającego
      absolutny rozkład  energii gwiazdy standardowej. Wymaga się aby widmo to
      miało dwa  pola typu  Real*4 lub Real*8, oraz XPTR=0 i  YPTR=1. Pierwsze
      pole powinno  zawierać długości fali, a drugie średnią energię gwiazdy w
      erg/cm^2/s/ w otoczeniu tej długości. Szerokość obszaru z którego brana
      jest średnia odpowiada parametrowi width.

   width
      typu Real, szerokość obszaru uśredniania (patrz opis srcsht)

   Punkty widma  wynikowego (histogramu)  zawierają informację  ile razy  dana
współrzędna miała  wartość mieszczącą  się w określonym przedziale. Dla type=0
analizowany jest  rozkład Y,  dla type=1  rozkład X.  Dodatkowo dla  type=1  w
oddzielnym  polu   zapisywana  jest   suma  współrzędnych  Y  punktów  których
współrzędne X są w określonym przedziale. Granice przedziału dla i-tego punktu
histogramu obliczane są następująco (tylko dla type=0 lub type=1) :

        DolnaGranicai = val1 + step * i
        GórnaGranicai = val1 + step * (i + 1)

Przedział jest domknięty od dołu. Widmo wynikowe będzie miało długość :

        int((val2 - val1) / step) + 1

Współczynniki skalujące będą miały wartości :

        x0 = val1 + step / 2
        x1 = step

Przykłady

   - obliczenie średniego poziomu energii widma w przedziałach co 100 A :

        LOAD 0 "R23.SHT"                        
        HISTO 1 1 x(0,0) x(0,last(0)) 100       
        SETA 1:0 1:0:i/1:1:i                    
        USE 0 1