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