funcnts - Online in der Cloud

PROGRAMM:

NAME/FUNKTION

funcnts – Photonen in bestimmten Regionen zählen, mit BKGD-Subtraktion

ZUSAMMENFASSUNG

Funktionen [Schalter] [source_region] [bkgd_file] [bkgd_region⎪bkgd_value]

OPTIONAL

-e "source_exposure[;bkgd_exposure]"
# Quelle (bkgd) Passt das Belichtungsbild unter Verwendung passender Dateien an
-w „source_exposure[;bkgd_exposure]“
# Quelle (bkgd) Passt das Belichtungsbild mithilfe der WCS-Transformation an
-t „source_timecorr[;bkgd_timecorr]“
# Quelle (bkgd) Zeitkorrekturwert oder Header-Parametername
-g # Ausgabe im schönen G-Format
-G # Ausgabe im %.14g-Format (maximale Genauigkeit)
-i "[column;]int1;int2..." # spaltenbasierte Intervalle
-m # stimmt mit einzelnen Quell- und Hintergrundregionen überein
-p # Ausgabe in Pixel, auch wenn wcs vorhanden ist
-r # gibt innere/äußere Radien (und Winkel) für Ringe (und Pandas) aus
-s # summierte Werte ausgeben
-v "scol[;bcol]" # src- und bkgd-Wertspalten für Tabellen
-T # Ausgabe im Starbase/RDB-Format
-z # Ausgaberegionen mit Nullfläche

BESCHREIBUNG

Funktionen zählt Photonen in den angegebenen Quellregionen und meldet die Ergebnisse für jede
Region. Regionen werden mithilfe des Spatial Region Filtering-Mechanismus angegeben. Photonen sind
Wird auch in den angegebenen bkgd-Bereichen gezählt, die auf dieselbe oder eine andere Datendatei angewendet werden
Datendatei. (Alternativ kann ein konstanter Hintergrundwert in Anzahl/Pixel**2 angegeben werden
angegeben.) Die bkgd-Regionen werden entweder eins zu eins mit Quellregionen gepaart oder gepoolt
und nach Fläche normalisiert und dann von den Quellzahlen in jeder Region subtrahiert.
Zu den angezeigten Ergebnissen gehören die bkgd-subtrahierten Zählungen in jeder Region sowie der Fehler
von den Zählungen, der Fläche in jeder Region und der Oberflächenhelligkeit (Zähler/Fläche**2)
für jede Region berechnet.

Das erste Argument des Programms gibt das FITS-Eingabebild, das Array oder die Rohereignisdatei an
herstellen. Wenn „stdin“ angegeben ist, werden Daten aus der Standardeingabe gelesen. Verwenden Sie Funtools
Klammernotation zur Angabe von FITS-Erweiterungen, Bildabschnitten und Filtern.

Das optionale zweite Argument ist der Deskriptor der Quellregion. Wenn keine Region angegeben ist,
das gesamte Feld wird genutzt.

Die Hintergrundargumente können eine von zwei Formen annehmen, je nachdem, ob sie separat sind
Hintergrunddatei angegeben ist. Soll die Quelldatei auch als Hintergrund verwendet werden, ist die
Das dritte Argument kann entweder der Hintergrundbereich oder ein konstanter Wert sein
Hintergrundanzahl/Pixel. Alternativ kann das dritte Argument eine Hintergrunddatendatei sein,
In diesem Fall ist das vierte Argument die Hintergrundregion. Wenn kein drittes Argument vorhanden ist
angegeben, wird ein konstanter Wert von 0 verwendet (dh kein Hintergrund).

Zusammenfassend sind die folgenden Befehlsargumente gültig:

[sh] funcnts sfile # zählt in der Quelldatei
[sh] funcnts sfile sregion # zählt in der Quellregion
[sh] funcnts sfile sregion bregion # bkgd reg. stammt aus der Quelldatei
[sh] funcnts sfile sregion bvalue # bkgd reg. ist konstant
[sh] funcnts sfile sregion bfile bregion # bkgd reg. stammt aus einer separaten Datei

Hinweis: Im Gegensatz zu anderen Funtools-Programmen werden Quell- und Hintergrundregionen als angegeben
separate Argumente in der Befehlszeile, anstatt sie als Teil in Klammern zu setzen
der Quell- und Hintergrunddateinamen. Dies liegt daran, dass Regionen in Funktionen nicht einfach sind
werden als Datenfilter, aber auch zur Berechnung von Flächen, Belichtung usw. verwendet
Geben Sie den Quellbereich innerhalb der Klammern an (verwenden Sie ihn also einfach als Filter), anstatt ihn anzugeben
Wenn Sie es als Argument zwei verwenden, zählt das Programm weiterhin nur Photonen, die den Regionsfilter passieren.
Die Flächenberechnung wird jedoch auf dem gesamten Feld durchgeführt, da Feld() lernen muss die
Standardquellregion. Dies ist selten das gewünschte Verhalten. Andererseits mit FITS
Bei Binärtabellen ist es oft sinnvoll, einen Spaltenfilter in die Klammern des Dateinamens einzufügen
Innerhalb der Region werden nur Ereignisse gezählt, die mit dem Spaltenfilter übereinstimmen.

Um beispielsweise die Zählungen innerhalb eines Radius von 22 Pixeln von der Mitte des zu extrahieren
Passt in die Binärtabelle snr.ev und subtrahiert den ermittelten Hintergrund aus demselben Bild darin
ein Ring mit Radien von 50-100 Pixeln:

[sh] funcnts snr.ev „circle(502,512,22)“ „annulus(502,512,50,100)“
# Quelle
# Datendatei: snr.ev
# Grad/Pixel: 0.00222222
# Hintergrund
# Datendatei: snr.ev
# Spalteneinheiten
# Bereich: Bogensekunden**2
# surf_bri: cnts/arcsec**2
# surf_err: cnts/arcsec**2

# vom Hintergrund subtrahierte Ergebnisse
reg net_counts Fehler Hintergrundfehler Bereich surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 3826.403 66.465 555.597 5.972 96831.98 0.040 0.001

# Folgende Quell- und Hintergrundkomponenten wurden verwendet:
Quellregion(en)
----------------
Kreis(502,512,22)

reg zählt Pixel
---- ------------ ---------
1 4382.000 1513

Hintergrundregion(en)
--------------------
Ring (502,512,50,100)

reg zählt Pixel
---- ------------ ---------
alle 8656.000 23572

Die Flächeneinheiten für die Ausgabespalten mit den Bezeichnungen „area“, „surf_bri“ (Oberflächenhelligkeit) und
„surf_err“ wird entweder in Bogensekunden angegeben (sofern entsprechende WCS-Informationen enthalten sind).
Datendatei-Header) oder in Pixel. Wenn die Datendatei WCS-Informationen enthält, Sie aber keine Arc-
zweite Einheiten, verwenden Sie die -p Wechseln Sie, um die Ausgabe in Pixeln zu erzwingen. Auch Regionen mit einer Fläche von Null
sind normalerweise nicht in der primären (vom Hintergrund subtrahierten) Tabelle enthalten, sind aber enthalten
in den sekundären Quell- und Bkgd-Tabellen. Wenn Sie möchten, dass diese Regionen in die einbezogen werden
Primärtabelle verwenden Sie die -z Schalter.

Beachten Sie, dass ein einfacher sed-Befehl die vom Hintergrund subtrahierten Ergebnisse zur weiteren Verwendung extrahiert
Analyse:

[sh] cat funcnts.sed
1,/---- .*/d
/^$/,$d

[sh] sed -f funcnts.sed funcnts.out
1 3826.403 66.465 555.597 5.972 96831.98 0.040 0.001

Wenn separate Quell- und Hintergrunddateien angegeben werden, Funktionen Ich werde versuchen, mich zu normalisieren
Passen Sie den Hintergrundbereich so an, dass die Hintergrundpixelgröße mit der Quellpixelgröße übereinstimmt
Größe. Diese Normalisierung kann nur erfolgen, wenn die entsprechenden WCS-Informationen vorhanden sind
in beiden Dateien enthalten sind (z. B. Grad-/Pixelwerte in CDELT). Wenn eine der beiden Dateien dies nicht tut
die erforderlichen Größeninformationen enthalten, wird die Normalisierung nicht durchgeführt. In diesem Fall,
Es liegt in der Verantwortung des Benutzers, sicherzustellen, dass die Pixelgrößen bei beiden gleich sind
Dateien.

Wenn mehr als eine Hintergrundregion angegeben ist, Funktionen werde sie alle kombinieren
in eine einzelne Region und verwenden Sie diese Hintergrundregion, um den hintergrundsubtrahierten zu erzeugen
Ergebnisse für jede Quellregion. Der -m (passt zu mehreren Hintergründen) Switch Tells Funktionen
um eine Eins-zu-Eins-Entsprechung zwischen Hintergrund- und Quellregionen herzustellen, anstatt
Verwendung einer einzigen kombinierten Hintergrundregion. Der Standardfall ist beispielsweise die Kombination von 2
Hintergrundregionen in eine einzelne Region umwandeln und diese Region dann auf jede Quelle anwenden
Regionen:

[sh] funcnts snr.ev "annulus(502,512,0,22,n=2)" "annulus(502,512,50,100,n=2)"
# Quelle
# Datendatei: snr.ev
# Grad/Pixel: 0.00222222
# Hintergrund
# Datendatei: snr.ev
# Spalteneinheiten
# Bereich: Bogensekunden**2
# surf_bri: cnts/arcsec**2
# surf_err: cnts/arcsec**2

# vom Hintergrund subtrahierte Ergebnisse
reg net_counts Fehler Hintergrundfehler Bereich surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 3101.029 56.922 136.971 1.472 23872.00 0.130 0.002
2 725.375 34.121 418.625 4.500 72959.99 0.010 0.000

# Folgende Quell- und Hintergrundkomponenten wurden verwendet:
Quellregion(en)
----------------
Ring(502,512,0,22,n=2)

reg zählt Pixel
---- ------------ ---------
1 3238.000 373
2 1144.000 1140

Hintergrundregion(en)
--------------------
Ring(502,512,50,100,n=2)

reg zählt Pixel
---- ------------ ---------
alle 8656.000 23572

Beachten Sie, dass die grundlegende Regionsfilterregel lautet: „Jedes Photon wird einmal gezählt, kein Photon jedoch.“
mehr als einmal gezählt" gilt weiterhin bei Verwendung von The -m um Hintergrundregionen anzupassen. Das
Das heißt, wenn sich zwei Hintergrundbereiche überlappen, werden die überlappenden Pixel nur in einem gezählt
von ihnen. Im schlimmsten Fall, wenn zwei Hintergrundregionen dieselbe Region sind, wird die
Der erste erhält alle Zählungen und Flächen, der zweite keine.

Verwendung der -m Schalter verursacht Funktionen um jede der beiden Hintergrundregionen unabhängig zu verwenden
mit jeder der beiden Quellregionen:

[sh] funcnts -m snr.ev "annulus(502,512,0,22,n=2)" "ann(502,512,50,100,n=2)"
# Quelle
# Datendatei: snr.ev
# Grad/Pixel: 0.00222222
# Hintergrund
# Datendatei: snr.ev
# Spalteneinheiten
# Bereich: Bogensekunden**2
# surf_bri: cnts/arcsec**2
# surf_err: cnts/arcsec**2

# vom Hintergrund subtrahierte Ergebnisse
reg net_counts Fehler Hintergrundfehler Bereich surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 3087.015 56.954 150.985 2.395 23872.00 0.129 0.002
2 755.959 34.295 388.041 5.672 72959.99 0.010 0.000

# Folgende Quell- und Hintergrundkomponenten wurden verwendet:
Quellregion(en)
----------------
Ring(502,512,0,22,n=2)

reg zählt Pixel
---- ------------ ---------
1 3238.000 373
2 1144.000 1140

Hintergrundregion(en)
--------------------
ann(502,512,50,100,n=2)

reg zählt Pixel
---- ------------ ---------
1 3975.000 9820
2 4681.000 13752

Beachten Sie, dass die meisten Gleitkommagrößen im „f“-Format angezeigt werden. Du kannst ändern
Dies in das „g“-Format mit dem -g schalten. Dies kann nützlich sein, wenn jedes Pixel gezählt wird
ist sehr klein oder sehr groß. Wenn Sie maximale Präzision wünschen und sich nicht darum kümmern
Säulen gut aneinanderreihen, verwenden -G, das alle Gleitkommawerte als %.14g ausgibt.

Bei der Zählung von Photonen mithilfe der Annulus- und Panda-Form (Kuchen und Annuli) ist dies häufig der Fall
nützlich, um Zugriff auf die Radien (und Panda-Winkel) für jede einzelne Region zu haben. Der -r
Der Schalter fügt der Ausgabetabelle Radien- und Winkelspalten hinzu:

[sh] funcnts -r snr.ev "annulus(502,512,0,22,n=2)" "ann(502,512,50,100,n=2)"
# Quelle
# Datendatei: snr.ev
# Grad/Pixel: 0.00222222
# Hintergrund
# Datendatei: snr.ev
# Spalteneinheiten
# Bereich: Bogensekunden**2
# surf_bri: cnts/arcsec**2
# surf_err: cnts/arcsec**2
# Radien: Bogensekunden
# Winkel: Grad

# vom Hintergrund subtrahierte Ergebnisse
reg net_counts Fehlerhintergrund Fehlerbereich surf_bri surf_err Radius1 Radius2 Winkel1 Winkel2
---- ------------ --------- ------------ --------- ---- ----- --------- --------- --------- --------- --------- ---------
1 3101.029 56.922 136.971 1.472 23872.00 0.130 0.002 0.00 88.00 NA NA
2 725.375 34.121 418.625 4.500 72959.99 0.010 0.000 88.00 176.00 NA NA

# Folgende Quell- und Hintergrundkomponenten wurden verwendet:
Quellregion(en)
----------------
Ring(502,512,0,22,n=2)

reg zählt Pixel
---- ------------ ---------
1 3238.000 373
2 1144.000 1140

Hintergrundregion(en)
--------------------
ann(502,512,50,100,n=2)

reg zählt Pixel
---- ------------ ---------
alle 8656.000 23572

Radien werden in Pixel- oder Bogensekundeneinheiten angegeben (abhängig vom Vorhandensein von WCS-Informationen).
während die Winkelwerte (falls vorhanden) in Grad angegeben sind. Diese Spalten können zum Plotten verwendet werden
Radialprofile. Zum Beispiel das Drehbuch funcnts.plot in der Funtools-Distribution) wird
Zeichnen Sie ein Radialprofil mit gnuplot (Version 3.7 oder höher). Eine vereinfachte Version davon
Das Skript ist unten dargestellt:

#!/ Bin / sh

if [ x"$1" = xgnuplot ]; Dann
if [ x`which gnuplot 2>/dev/null` = x ]; Dann
echo „FEHLER: Gnuplot nicht verfügbar“
Beenden Sie 1
fi
awk '
BEGIN{HEADER=1; DATEN=0; DATEIEN=""; XLABEL="unbekannt"; YLABEL="unbekannt"}
HEADER==1{
if( $1 == "#" && $2 == "data" && $3 == "file:" ){
if( FILES != "" ) FILES = FILES ","
DATEIEN = DATEIEN $4
}
sonst if( $1 == "#" && $2 == "radii:" ){
XLABEL = 3 $
}
sonst if( $1 == "#" && $2 == "surf_bri:" ){
YLABEL = 3 $
}
sonst if( $1 == "----" ){
printf "set nokey; set title \"funcnts(%s)\"\n", FILES
printf "set xlabel \" radius(%s)\"\n", XLABEL
printf "set ylabel \"surf_bri(%s)\"\n", YLABEL
print „plot \“-\“ using 3:4:6:7:8 with boxerrorbars“
HEADER = 0
DATEN = 1
weiter
}
}
DATEN==1{
if( NF == 12 ){
Drucken Sie $9, $10, ($9+$10)/2, $7, $8, $7-$8, $7+$8, $10-$9
}
sonst {
wunsch
}
}
' ⎪ gnuplot -persist - 1>/dev/null 2>&1

elif [ x"$1" = xds9 ]; Dann
awk '
BEGIN{HEADER=1; DATEN=0; XLABEL="unbekannt"; YLABEL="unbekannt"}
HEADER==1{
if( $1 == "#" && $2 == "data" && $3 == "file:" ){
if( FILES != "" ) FILES = FILES ","
DATEIEN = DATEIEN $4
}
sonst if( $1 == "#" && $2 == "radii:" ){
XLABEL = 3 $
}
sonst if( $1 == "#" && $2 == "surf_bri:" ){
YLABEL = 3 $
}
sonst if( $1 == "----" ){
printf „funcnts(%s) radius(%s) surf_bri(%s) 3\n“, FILES, XLABEL, YLABEL
HEADER = 0
DATEN = 1
weiter
}
}
DATEN==1{
if( NF == 12 ){
Drucken Sie 9 $, 7 $, 8 $
}
sonst {
wunsch
}
}
'
sonst
echo "funcnts -r ... ⎪ funcnts.plot [ds9⎪gnuplot]"
Beenden Sie 1
fi

Also laufen Funktionen und zeichnen Sie die Ergebnisse mit gnuplot (Version 3.7 oder höher) auf, verwenden Sie:

funcnts -r snr.ev "annulus(502,512,0,50,n=5)" ... ⎪ funcnts.plot gnuplot

Das -s (Summe) Schalterursachen Funktionen um eine zusätzliche Tabelle mit summierten (integrierten)
im Hintergrund subtrahierte Werte, zusammen mit der Standardtabelle der einzelnen Werte:

[sh] funcnts -s snr.ev "annulus(502,512,0,50,n=5)" "annulus(502,512,50,100)"
# Quelle
# Datendatei: snr.ev
# Grad/Pixel: 0.00222222
# Hintergrund
# Datendatei: snr.ev
# Spalteneinheiten
# Bereich: Bogensekunden**2
# surf_bri: cnts/arcsec**2
# surf_err: cnts/arcsec**2

# summierte, vom Hintergrund subtrahierte Ergebnisse
bis zu net_counts Fehler Hintergrundfehler Bereich surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 2880.999 54.722 112.001 1.204 19520.00 0.148 0.003
2 3776.817 65.254 457.183 4.914 79679.98 0.047 0.001
3 4025.492 71.972 1031.508 11.087 179775.96 0.022 0.000
4 4185.149 80.109 1840.851 19.786 320831.94 0.013 0.000
5 4415.540 90.790 2873.460 30.885 500799.90 0.009 0.000

# vom Hintergrund subtrahierte Ergebnisse
reg zählt Fehler Hintergrundfehler Bereich surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 2880.999 54.722 112.001 1.204 19520.00 0.148 0.003
2 895.818 35.423 345.182 3.710 60159.99 0.015 0.001
3 248.675 29.345 574.325 6.173 100095.98 0.002 0.000
4 159.657 32.321 809.343 8.699 141055.97 0.001 0.000
5 230.390 37.231 1032.610 11.099 179967.96 0.001 0.000

# Folgende Quell- und Hintergrundkomponenten wurden verwendet:
Quellregion(en)
----------------
Ring(502,512,0,50,n=5)

reg zählt Pixel sumcnts sumpix
---- ------------ --------- ------------ ---------
+1 2993.000 305 2993.000
+2 1241.000 940 4234.000
+3 823.000 1564 5057.000
+4 969.000 2204 6026.000
+5 1263.000 2812 7289.000

Hintergrundregion(en)
--------------------
Ring (502,512,50,100)

reg zählt Pixel
---- ------------ ---------
alle 8656.000 23572

Das -t und -e Mit Schaltern können Zeit- und Belichtungskorrekturen vorgenommen werden.
zu den Daten. Bitte beachten Sie, dass diese Korrekturen qualitativ verwendet werden sollen, da
Die Anwendung genauerer Korrekturfaktoren ist ein komplexer und auftragsabhängiger Aufwand.
Der Algorithmus zur Anwendung dieser einfachen Korrekturen lautet wie folgt:

C = Rohzahlen in der Quellregion
Ac= Fläche der Quellregion
Tc=Belichtungszeit für Quelldaten
Ec= Durchschnittliche Exposition in der Quellregion, aus der Expositionskarte

B = Rohzählungen im Hintergrundbereich
Ab= Fläche der Hintergrundregion
Tb= (Belichtungs-)Zeit für Hintergrunddaten
Eb = Durchschnittliche Belichtung im Hintergrundbereich, aus der Belichtungskarte

Dann beträgt die Nettoanzahl in der Quellregion

Netto= C - B * (Ac*Tc*Ec)/(Ab*Tb*Eb)

mit der Standardausbreitung von Fehlern für den Fehler im Netz. Der Nettopreis wäre dann

Nettorate = Netto/(Ac*Tc*Ec)

Die durchschnittliche Belichtung in jeder Region wird durch Summieren der Pixelwerte in der berechnet
Belichtungskarte für die angegebene Region erstellen und dann durch die Anzahl der darin enthaltenen Pixel dividieren
Region. Belichtungskarten werden häufig mit einem Blockfaktor > 1 generiert (z. B. Block 4 bedeutet das).
jedes Belichtungspixel enthält 4x4 Pixel bei voller Auflösung) und Funktionen werde mich damit befassen
automatisch sperren. Verwendung der -e Mit dem Schalter können Sie sowohl Quelle als auch Hintergrund bereitstellen
Belichtungsdateien (getrennt durch „;“), wenn Sie über separate Quell- und Hintergrunddatendateien verfügen.
Wenn Sie keine Hintergrundbelichtungsdatei bereitstellen, gehen Sie zu separaten Hintergrunddaten
Datei, Funktionen geht davon aus, dass die Belichtung bereits auf die Hintergrunddatendatei angewendet wurde.
Darüber hinaus wird davon ausgegangen, dass der Pixelfehler in der Hintergrunddatendatei Null ist.

NB: Der -e Der Schalter geht davon aus, dass die Belichtungskarte die Bilddatei überlagert genau, außer
für den Blockfaktor. Jedes Pixel im Bild wird um den Blockfaktor skaliert, um darauf zugreifen zu können
entsprechenden Pixel in der Belichtungskarte. Wenn Ihre Belichtungskarte nicht genau übereinstimmt
mit dem Bild, do nicht - -e Belichtungskorrektur. In diesem Fall ist es immer noch möglich
um eine Belichtungskorrektur durchzuführen if Sowohl das Bild als auch die Belichtungskarte verfügen über ein gültiges WCS
Informationen: Verwenden Sie die -w Schalten Sie so um, dass die Transformation vom Bildpixel zur Belichtung erfolgt
Pixel verwendet die WCS-Informationen. Das heißt, jedes Pixel im Bildbereich wird sein
zuerst von Bildkoordinaten in Himmelskoordinaten umgewandelt, dann von Himmelskoordinaten in
Belichtungskoordinaten. Bitte beachten Sie, dass mit -w kann die Bearbeitungszeit verlängern
die Belichtungskorrektur erheblich.

Mithilfe von kann eine Zeitkorrektur sowohl auf Quell- als auch auf Hintergrunddaten angewendet werden -t Schalter.
Der Wert für die Korrektur kann entweder eine numerische Konstante oder der Name eines Headers sein
Parameter in der Quelldatei (oder Hintergrunddatei):

[sh] funcnts -t 23.4 ... # Nummer für Quelle
[sh] funcnts -t "LIVETIME;23.4" ... # Parameter für Quelle, numerisch für Hintergrund

Wenn eine Zeitkorrektur angegeben ist, wird diese auch auf die Nettozählungen angewendet (siehe
Algorithmus oben), so dass die Einheiten der Oberflächenhelligkeit cnts/Fläche**2/Sek. werden.

Das -i Zum Ausführen wird ein (Intervall-)Schalter verwendet Funktionen auf mehreren spaltenbasierten Intervallen mit
nur ein einziger Durchgang durch die Daten. Es ist gleichbedeutend mit Laufen Funktionen mehrmals
wobei jedes Mal ein anderer Spaltenfilter zu den Quell- und Hintergrunddaten hinzugefügt wird. Für jede
Intervall, das volle Funktionen Die Ausgabe wird generiert, wobei ein Zeilenvorschubzeichen (^L) eingefügt wird
zwischen jedem Lauf. Darüber hinaus enthält die Ausgabe für jedes Intervall das Intervall
Spezifikation im Header. Intervalle sind für die Erzeugung der Röntgenhärte sehr nützlich
Verhältnisse effizient zu gestalten. Natürlich werden sie nur unterstützt, wenn die Eingabedaten enthalten sind
in einer Tabelle.

Für die Intervallangabe werden zwei Formate unterstützt. Das allgemeinste Format ist semi-
Durch Doppelpunkte getrennte Liste von Filtern, die als Intervalle verwendet werden sollen:

funcnts -i "pha=1:5;pha=6:10;pha=11:15" snr.ev "circle(502,512,22)" ...

Vom Konzept her entspricht dies dem Laufen Funktionen drei Mal:

funcnts snr.ev'[pha=1:5]' "circle(502,512,22)"
funcnts snr.ev'[pha=6:10]' "circle(502,512,22)"
funcnts snr.ev'[pha=11:15]' "circle(502,512,22)"

Verwenden Sie jedoch die -i Der Switch erfordert nur einen Datendurchlauf.

Beachten Sie, dass komplexe Filter zur Angabe von Intervallen verwendet werden können:

funcnts -i "pha=1:5&&pi=4;pha=6:10&&pi=5;pha=11:15&&pi=6" snr.ev ...

Das Programm lässt die Daten einfach nacheinander durch jeden Filter laufen und generiert drei Funktionen
Ausgaben, getrennt durch das Zeilenumbruchzeichen.

Tatsächlich besteht die Absicht darin, Intervalle für Härteverhältnisse zu unterstützen, die spezifiziert sind
Filter müssen überhaupt keine Intervalle sein. Es muss auch kein „Intervall“-Filter vorhanden sein
mit einem anderen verwandt. Zum Beispiel:

funcnts -i "pha=1:5;pi=6:10;energy=11:15" snr.ev "circle(502,512,22)" ...

ist gleichbedeutend mit Laufen Funktionen dreimal mit unabhängigen Filterspezifikationen.

Für den einfachen Fall, dass eine einzelne Spalte verwendet wird, wird ein zweites Intervallformat unterstützt
um mehrere homogene Intervalle für diese Spalte anzugeben. In diesem Format ein Spaltenname
wird zuerst angegeben, gefolgt von Intervallen:

funcnts -i "pha;1:5;6:10;11:15" snr.ev "circle(502,512,22)" ...

Dies entspricht dem ersten Beispiel, erfordert jedoch weniger Tipparbeit. Der Funktionen Programm
wird einfach „pha=" vor jedem der angegebenen Intervalle voranstellen. (Beachten Sie, dass dieses Format
enthält nicht das Zeichen „=“ im Spaltenargument.)

Normalerweise, wenn Funktionen wird auf einer FITS-Binärtabelle (oder einer Rohereignistabelle) ausgeführt, eine
Die Integralzahl wird für jede Zeile (Ereignis) akkumuliert, die in einem bestimmten Bereich enthalten ist. Der -v
„scol[;bcol]“ (Wertspalte) Der Schalter akkumuliert die Zählungen unter Verwendung des Werts aus der
angegebene Spalte für das angegebene Ereignis. Wenn nur eine einzelne Spalte angegeben ist, wird diese verwendet
sowohl die Quell- als auch die Hintergrundregion. Zwei separate Spalten, getrennt durch ein Semikolon,
kann für Quelle und Hintergrund angegeben werden. Dazu kann der spezielle Token „$none“ verwendet werden
Geben Sie an, dass eine Wertespalte für das eine, aber nicht für das andere verwendet werden soll. Zum Beispiel,
'pha;$none' verwendet die pha-Spalte für die Quelle, verwendet jedoch ganzzahlige Zählungen für die
Hintergrund, während '$none;pha' das Gegenteil bewirkt. Wenn die Wertespalte vom Typ ist
logisch, dann ist der verwendete Wert 1 für T und 0 für F. Es werden Wertespalten verwendet, z
Zum Beispiel, um Wahrscheinlichkeiten anstelle von Integralzahlen zu integrieren.

Besitzt das -T (rdb table) Schalter verwendet wird, entspricht die Ausgabe der Starbase/rdb-Datenbank
Format: Zwischen den Spalten werden Tabulatoren anstelle von Leerzeichen und ein Zeilenvorschub eingefügt
zwischen Tabellen eingefügt.

Beachten Sie schließlich, dass Funktionen ist ein Image-Programm, obwohl es direkt auf FITS ausgeführt werden kann
Binärtabellen. Dies bedeutet, dass zur Gewährleistung eine Bildfilterung auf die Zeilen angewendet wird
dass die gleichen Ergebnisse erzielt werden, unabhängig davon, ob eine Tabelle oder eine entsprechende Klasse gruppiert wird
Bild verwendet wird. Aus diesem Grund ist die Anzahl der gefundenen Zählungen jedoch gering Funktionen kann
unterscheiden sich von der Anzahl der Ereignisse, die mithilfe von Zeilenfilterprogrammen gefunden wurden, z Fundisp or
lustig Weitere Informationen zu diesen Unterschieden finden Sie in der Diskussion über Region
Grenzen.

Nutzen Sie Funktionen online über die Dienste von onworks.net