Remeis-Wiki - User contributions [en]

Skivergnuegen 2023

2023-01-05T17:15:18Z

Falkner:

'''Remeis Skifahren 2023''' The same procedure as every year ;-)
Hier werden die wichtigsten organisatorischen Dinge festgehalten. Jeder darf hier gerne editieren; wer keinen Account für das Wiki hat, muss sich bitte an jemanden wenden, der Zugang hat.

[[File:154854 silvretta-montafon-panorama.jpeg|right|700px]]

== Allgemeine Infos ==

*Ferienwohnung(en):
**Wir kommen wieder unter im [https://www.grandau.at/montafon/zimmer/ferienhaus-enzian Haus Enzian] des [https://www.grandau.at/ Sporthotels Grandau].
**Adresse des Sporthotel: Montafonerstraße 274a, 6791 St. Gallenkirch, Österreich. Unsere Unterkünfte liegen im Türkeiweg.

*Skigebiet: Montafon (Vorarlberg)
**Skigebiet Karte: https://winter.intermaps.com/montafon?lang=de
**Skipass Preise: https://www.silvretta-montafon.at/de/onlineshop/ticket-uebersicht
**Wer nicht alle 7 Tage Skifahren möchte: Es gibt Angebote wie z.B. 5 aus 6, mit einem solchen Skipass kann man innerhalb von 6 Tagen an 5 beliebigen Tagen skifahren.

*Skiverleih: [http://www.sportharry.at/ Sport Harry], direkt an der Talstation.

== Unterkunft ==
=== Aufteilung ===

Dieses Jahr gibt es nur die Hütte mit folgenden Zimmern (wir sind insgesamt 14, also bleiben keine Plätze frei).

{| class="wikitable"
! scope="col" | Zimmer
! scope="col" | Insassen
|-
! scope="row" | 2-er
| Amy, Aafia

|-
! scope="row" | 3-er
| Jakob, Max, Katya
|-
! scope="row" | 4-er
| Eva, Steffen, Flo, Christian
|-
! scope="row" | 5-er
| Basti, Philipp
|}

== Fahrer ==
Fahrer können sich hier eintragen und Eckdaten angeben. Diejenige, die mitfahren möchten, sprechen sich mit den Fahrern ab und tragen sich ebenfalls ein. Klärt bitte auch die Gepäcklage, ggf. kann ein anderer Fahrer z.B. Skiausrüstung mitnehmen.

'''Auf der Suche nach einer Mitfahrgelegenheit: '''

{| class="wikitable"
! scope="col" | Fahrer
! scope="col" | # Plätze
! scope="col" | Anreise
! scope="col" | Mitfahrer Anreise
! scope="col" | Abreise
! scope="col" | Mitfahrer Abreise
! scope="col" | Skimitnahme
! scope="col" | Kommentar
|-
|Eva
|0(3)
|Samstag
|Steffen+Christian
|Samstag?
|Steffen+Christian
|eigene ja
|Auto ist voll
|-
|Max
|0(3)
|Samstag (from Arlberg)
|Jakob+Katya
|Samstag
|Jakob+Katya
|eigene ja
|Auto ist voll
|-
|Basti
|2(2)
|Sunday (early morning)
|
|Friday (lunch time)
|
|Ja
|I will be at the airport NUE (alternative meeting point would be Nuremberg trainstation)
|-

|-

|-

|-
|}

== Schwarzes Brett ==

=== Essen ===
Hier gibt's den Essensplan, soweit übernommen von den letzten Jahren. Da wir sehr viele Leute sind, sollten wir wieder vorab grob planen, was wir kochen wollen. Wir können einige Dinge in Deutschland besorgen und mitnehmen, es gib aber auch einen nahe gelegen Supermarkt. Generelle Dinge, wie z.B. Gewürze oder Aufstriche für Frühstücksbrötchen könnte man auch mitbringen.

Wer Kochvorschläge oder andere Ideen hat, gerne unten eintragen und kommentieren. Bedenkt, dass die Zubereitung (relativ) einfach sein sollte. Da sich die Liste in den letzten Jahren gut bewaehrt hat, behalten wir sie bei, aber die Reihenfolge/Tage kann man sicher noch hin und her bewegen. Einige (wenige) Sachen müssen wir dann manchmal trotzdem noch beim Spar kaufen, da ja auch der Kühlschrank nur begrenzt Platz hat.

'''Essenplan:'''
*Samstag: Spaghetti Bolognese + veg. Bolognese
*Sonntag: Käse Spätzle
*Montag: Burritos
*Dienstag: Risotto
*Mittwoch: Gulasch
*Donnerstag: Kartoffeln mit Kräuter-Quark
*Freitag: Dal mit Reis (od. Curry mit Suesskartoffeln )

Einkaufsliste von 2022 (letztes Jahr): https://docs.google.com/document/d/166b_mIXxmoAgcv1upIGremd_36E1fmwwxsQADP2uEv8/

@Christian: Vegane Option entwender selbst mitbringen oder hier mit reinschreiben (falls man die einfach bekommt)

=== Alkohol/Party ===

=== Spiele ===

Immer lustig sind gemeinsame Spieleabende. Also, wer im Besitz von Gesellschaftsspielen ist, bringt diese gerne mit! Um eine Übersicht zu bekommen auch bitte hier eintragen:

''* Spiele hier eintragen''

=== Sonstiges ===

=== Alternativprogramm ===
* Rodeln (http://www.montafon.at/de/urlaubswelten/echte_naturliebhaber/rodeln)
* Schneeschuhwandern (http://www.montafon.at/schneeschuhwanderungen)
* Therme (http://www.montafon.at/schwimmen), z.B. http://www.aqua-dome.at/de (ca. 130km entfernt!)

[[Category:Internal]]

PhD hat for Maria

2019-05-21T14:32:20Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle ('''**Max**''')
** a beeper with button to start alarm (Jakob takes a look at this)
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike (with trailer) + figure (maybe fierfighter figure by playmobile?) with colored hair and shoes of different color ('''**Basti**''')
* whiskey (ardberg supernova): klopfer bottle ('''**Simon/Basti**''')

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS?] + accretion disk also sheep but colored wool?)
** BH as half sphere (styrofoam size of a tennis ball) ('''**Max**''')
** accretion disk made of colored wool ('''**Katya**''')
** clumpy wind: wool sheep on wire [similar like [https://i.pinimg.com/originals/3a/1e/e8/3a1ee866215e17f927e184652dc9a895.jpg Shaun das /Schaf]] ('''**Andrea**''')
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) ('''**Christian**''')
* Web admin: extensive fights with WordPress -> small print for the side of the hat ('''**Matthias**''')
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator ('''**Max**''')

[[Category:Internal]]

PhD hat for Maria

2019-05-11T18:10:29Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle ('''**Max**''')
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike (with trailer) + figure (maybe fierfighter figure by playmobile?) with colored hair and shoes of different color ('''**Basti**''')
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS?] + accretion disk also sheep but colored wool?)
** BH as half sphere (styrofoam size of a tennis ball)
** accretion disk made of colored wool
** clumpy wind: wool sheep on wire [similar like [https://i.pinimg.com/originals/3a/1e/e8/3a1ee866215e17f927e184652dc9a895.jpg Shaun das /Schaf]] ('''**Andrea**''')
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) ('''**Christian**''')
* Web admin: extensive fights with WordPress -> small print for the side of the hat ('''**Matthias**''')
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-11T18:10:06Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle ('''**Max**''')
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike (with trailer) + figure (maybe fierfighter figure by playmobile?) with colored hair and shoes of different color ('''**Basti**''')
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS?] + accretion disk also sheep but colored wool?)
** BH as half sphere (styrofoam size of a tennis ball)
** accretion disk made of colored wool
** clumpy wind: wool sheep on wire [similar like [https://i.pinimg.com/originals/3a/1e/e8/3a1ee866215e17f927e184652dc9a895.jpg Shaun das /Schaf]] ('''**Andrea**''')
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) ('''**Christian**''')
* Web admin: extensive fights with WordPress -> small print for the side of the hat ('''**Matthias**'''s)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-11T18:08:47Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike (with trailer) + figure (maybe fierfighter figure by playmobile?) with colored hair and shoes of different color (**Basti**)
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS?] + accretion disk also sheep but colored wool?)
** BH as half sphere (styrofoam size of a tennis ball)
** accretion disk made of colored wool
** clumpy wind: wool sheep on wire [similar like [https://i.pinimg.com/originals/3a/1e/e8/3a1ee866215e17f927e184652dc9a895.jpg Shaun das /Schaf]] (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-11T17:58:32Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike (with trailer) + figure (maybe fierfighter figure by playmobile?) with colored hair and shoes of different color (**Basti**)
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep [similar like [https://i.pinimg.com/originals/3a/1e/e8/3a1ee866215e17f927e184652dc9a895.jpg Shaun das /Schaf]] (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-09T07:09:05Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike with trailer + figure with colored hair and shoes of different color
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-09T07:08:01Z

Falkner:

== PhD Hat for Maria ==

====== list of ideas ======

* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

* line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike with trailer + figure with colored hair and shoes of different color
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-09T07:06:57Z

Falkner:

== PhD Hat for Maria ==
====== list of ideas ======

==== Private ====

* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

** line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

== building plan ==

====== Private ======

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike with trailer + figure with colored hair and shoes of different color
* whiskey (ardberg supernova)

====== Work ======

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-09T07:06:09Z

Falkner:

====== PhD Hat for Maria ======

====== list of ideas ======

==== Private ====

* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

** line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

====== building plan ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* Cycling from Erlangen to Remeis
** ESTW smokestack + canal leading to observatory + bike with trailer + figure with colored hair and shoes of different color
* whiskey (ardberg supernova)

==== Work ====

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor: best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat; Images: /home/falkner/work/PhDhat)
* first-aider @ remeis: first-aid kit & defibrillator

[[Category:Internal]]

PhD hat for Maria

2019-05-07T07:19:09Z

Falkner:

====== PhD Hat for Maria ======

====== list of ideas ======

==== Private ====

* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* sheep!
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

** line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations

====== building plan ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* mit dem Fahrrad nach Remeis (von Erlangen): ESTW smokestack + canal leading to observatory + bike with trailer
* whiskey (ardberg supernova)

==== Work ====

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor
** best of CCD protocols printout [ccd protocol layout]
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat, still looking for the images)
** https://pbs.twimg.com/media/DfVoxMjWsAEsz12.jpg:small
** https://visit.cern/sites/visits.web.cern.ch/files/images/image/shop-21.jpg

[[Category:Internal]]

PhD hat for Maria

2019-05-06T13:46:47Z

Falkner:

====== PhD Hat for Maria ======

====== list of ideas ======

==== Private ====

* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* sheep!
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

** line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat, still looking for the images)

====== building plan ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* mit dem Fahrrad nach Remeis (von Erlangen): ESTW smokestack + canal leading to observatory + bike with trailer

==== Work ====

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor
** best of CCD protocols printout [ccd protocol layout]

[[Category:Internal]]

PhD hat for Maria

2019-04-30T11:53:07Z

Falkner:

====== PhD Hat for Maria ======

====== list of ideas ======

==== Private ====

* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* sheep!
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

** line measurements (remeis intern poll)
* timing analysis (power spectra)
* constellations
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat, still looking for the images)

====== building plan ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter
** firetruck with hose as tassel : Paper model + cord + nozzle (**Max**)
** a beeper with button to start alarm
* mit dem Fahrrad nach Remeis (von Erlangen): ESTW smokestack + canal leading to observatory + bike with trailer

==== Work ====

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk also sheep but colored wool?)
** wool sheep (**Andrea**)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* CCD lab instructor
** best of CCD protocols printout [ccd protocol layout]

[[Category:Internal]]

PhD hat for Maria

2019-04-30T11:24:51Z

Falkner:

====== PhD Hat for Maria ======

====== list of ideas ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* sheep!
* mit dem Fahrrad nach Remeis (von Erlangen)
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk?)
** line measurements (remeis intern poll)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* timing analysis (power spectra)
* CCD lab instructor
* constellations
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat, still looking for the images)

====== building plan ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat => accretion disk (*Basti*)
* firefighter (beeper goes of at inconvenient times)
** this actually includes paramedic training (with ride-alongs in ambulances)
** first-aid @ Remeis
* colored hair
* shoes and socks of different color
* sheep!
* mit dem Fahrrad nach Remeis (von Erlangen)
* kids as little thesis helpers
* card/board games
* wiskey

==== Work ====

* Cyg X-1 (Draft Idea: build BH + Clumpy wind {clumps are sheep shaped} [wool & LEDS] + accretion disk?)
** line measurements (remeis intern poll)
* extended work on paper: LightCurve with words/day or git commits/day (get git repo) (**Christian**)
* timing analysis (power spectra)
* CCD lab instructor
* constellations
* Web admin: extensive fights with WordPress -> small print for the side of the hat (**Matthias**)
* Teilchenweltvermittler (-> little "alien" particle stickers for the side of the hat, still looking for the images)

[[Category:Internal]]

PhD hat for Maria

2019-04-23T07:32:24Z

Falkner: /* Work */

====== PhD Hat for Maria ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat (*Basti*)
* firefighter (beeper goes of at inconvenient times)
* colored hair
* shoes and socks of different color

==== Work ====

* Cyg X-1
* extended work on paper
* line measurements (remeis intern poll)
* timing analysis (power spectra)

[[Category:Internal]]

PhD hat for Maria

2019-04-23T07:32:05Z

Falkner: /* Work */

====== PhD Hat for Maria ======

==== Private ====

* colored hat: glue rainbow-paper on the hat base (and top) for a colorful hat (*Basti*)
* firefighter (beeper goes of at inconvenient times)
* colored hair
* shoes and socks of different color

==== Work ====

* Cyg X-1
* work on paper
* line measurements (remeis intern poll)
* timing analysis (power spectra)

[[Category:Internal]]

PhD hat for Maria

2019-04-23T07:29:02Z

Falkner: /* Work */

PhD hat for Maria

2019-04-23T07:28:46Z

Falkner: /* Private */

PhD hat for Maria

2019-04-23T07:23:47Z

Falkner: Created page with " ====== PhD Hat for Maria ====== ==== Private ==== ==== Work ==== Category:Internal"

====== PhD Hat for Maria ======

==== Private ====

==== Work ====

[[Category:Internal]]

ERosita Coding standards

2018-05-02T11:31:39Z

Falkner: Created page with " Category:Internal ===== NRTA library coding standards ===== * **variables**: written in lowercase <code c> int name_of_variable;</code> * **constants**: marked w..."

[[Category:Internal]]

===== NRTA library coding standards =====

* **variables**: written in lowercase <code c>
int name_of_variable;</code>
* **constants**: marked with 'c' at the beginning <code c>
const int cname_of_variable;</code>
* **statics**: marked with 's' at the beginning <code c>
static int sname_of_variable;</code>
* **globals**: marked with 'g' at the beginning <code c>
int gname_of_variable;</code>
* **pointers (and arrays)**: marked with '*p' at the beginning <code c>
int *pname_of_pointer;
long pnaxes[2] = {xsize, ysize};</code>
* **loop variables**: use double letters, starting with 'ii' and 'jj' <code c>
int ii, jj;
for(ii = 0; ii < 50; ii++) {
for(jj = 0; jj < 10; jj++) {
printf("do something...\n");
}
}</code>
* **functions**: all functions start with 'ero_' and have '_' between separate words <code c>
void ero_name_of_function(void);</code>
* **structures**: start with capital letter <code c>
struct Structure {
int blah;
char some_string[80];
(void)(*pero_some_function)(int, char*);
};</code>
* **enums**: to be favored over preprocessor directives (#define ...); name of enum prefixes each entry <code c>
enum List {
LIST_ENTRY_A = 1,
LIST_ENTRY_B = 2,
LIST_ENTRY_C = 3
};</code>
* **typedefs**: objects to end with '_t' <code c>
typedef unsigned int size_t;</code>

Cinema

2018-05-02T11:28:46Z

Falkner:

[[Category:Internal]]

===== Remeis Movie Club =====

The schedule for upcoming movies at the Remeis theater. Movies usually start around 7-8 pm. Dinner plans are made on short notice (default is pizza).\\
Please do not just edit any dates, they're usually coordinated between some people already when they are entered into the table.

{|
| Week of June 3
|-
| Tue, June 4
| Waterworld
| Cloud Atlas
|-
| Week of June 10
|-
| Wed, June 12
| Game of Thrones S03E7-10
|-
| Week of June 17
|-
| Wed, June 19
| End of Watch
|-
| Week of July 8
|-
| TBD
| Terminator Series
|}

Also planned (feel free to add wishes):

True Grit\\
Thor\\
Gattaca\\
Highlander\\
Back to the Future I-III\\
\\
Iron Man 3 [October]\\
Marvel's Agents of S.H.I.E.L.D. [Series coming in fall 2013]\\
\\
Terminator\\
A clockwork orange\\
Dracula\\
Some like it hot\\

Cinema

2018-05-02T11:27:59Z

Falkner:

[[Category:Internal]]

===== Remeis Movie Club =====

The schedule for upcoming movies at the Remeis theater. Movies usually start around 7-8 pm. Dinner plans are made on short notice (default is pizza).\\
Please do not just edit any dates, they're usually coordinated between some people already when they are entered into the table.

{|
| Week of June 3
|-
| Tue, June 4
| Waterworld
| Cloud Atlas
|-
| Week of June 10
|-
| Wed, June 12
| Game of Thrones S03E7-10
|-
| Week of June 17
|-
| Wed, June 19
|-
| End of Watch
|-
| Week of July 8
|-
| TBD
| Terminator Series
|}

Also planned (feel free to add wishes):

True Grit\\
Thor\\
Gattaca\\
Highlander\\
Back to the Future I-III\\
\\
Iron Man 3 [October]\\
Marvel's Agents of S.H.I.E.L.D. [Series coming in fall 2013]\\
\\
Terminator\\
A clockwork orange\\
Dracula\\
Some like it hot\\

Cinema

2018-05-02T11:23:04Z

Falkner: Created page with " Category:Internal ===== Remeis Movie Club ===== The schedule for upcoming movies at the Remeis theater. Movies usually start around 7-8 pm. Dinner plans are made on s..."

[[Category:Internal]]

===== Remeis Movie Club =====

The schedule for upcoming movies at the Remeis theater. Movies usually start around 7-8 pm. Dinner plans are made on short notice (default is pizza).\\
Please do not just edit any dates, they're usually coordinated between some people already when they are entered into the table.

^ Week of June 3 ^^^
| Tue, June 4 | Waterworld | Cloud Atlas |
^ Week of June 10 ^^^
| Wed, June 12 | Game of Thrones S03E7-10 ||
^ Week of June 17 ^^^
| Wed, June 19 | End of Watch ||
^ Week of July 8 ^^^
| TBD | Terminator Series ||

\\

Also planned (feel free to add wishes):

True Grit\\
Thor\\
Gattaca\\
Highlander\\
Back to the Future I-III\\
\\
Iron Man 3 [October]\\
Marvel's Agents of S.H.I.E.L.D. [Series coming in fall 2013]\\
\\
Terminator\\
A clockwork orange\\
Dracula\\
Some like it hot\\

Intern:tanami:sources

2018-05-02T11:21:13Z

Falkner: Created page with "====== TANAMI sources ======"

====== TANAMI sources ======

TANAMI

2018-05-02T11:20:26Z

Falkner: Created page with " Category:Internal ====== TANAMI ====== TANAMI sources ===== Radio ===== Radio Analysis Spectral Index Maps ===== Optical/UV ===== //Swif..."

[[Category:Internal]]

====== TANAMI ======
[[intern:tanami:sources | TANAMI sources]]

===== Radio =====

Radio Analysis
Spectral Index Maps

===== Optical/UV =====
//Swift///UVOT
//XMM///OM

===== X-ray =====

===== Gamma-ray =====
Fermi

Proposals: XMM

2018-05-02T11:17:22Z

Falkner: Created page with "====== XMM Proposals ====== == Simultaneous XMM-Newton and Multiwavelength Observations of the Nucleus of CentaurusA == PI: Cornelia Müller, J. Wilms, M. Kadler, R. Ojha, A..."

====== XMM Proposals ======

== Simultaneous XMM-Newton and Multiwavelength Observations of the Nucleus of CentaurusA ==
PI: Cornelia Müller, J. Wilms, M. Kadler, R. Ojha, A. Markowitz, M. Böck, T. Beuchert, T. Dauser, F. Krauß, E. Litzinger, K. Mannheim, R. Rothschild, F. Spanier, M. Weidinger and the TANAMI Collaboration (AO12)

== Constraining the nature of \pmn: Is it a $\gamma$-loud young radio galaxy? ==
PI: Cornelia Müller, W. McConville, J. Wilms, M. Kadler, R. Ojha,F. Krauß, T. Dauser, A. Markowitz, M. Böck, C.C. Cheung, L. Stawarz, M.\,Dutka (AO12)

Proposals: Suzaku

2018-05-02T11:16:45Z

Falkner: Created page with " ====== Suzaku Proposals ====== ===== Accreting Pulsar Sample ===== == Cyclotron Lines in Transient Pulsars I: Probing the B-Field == approved 2012, Suzaku AO-7, NASA\\ PI..."

====== Suzaku Proposals ======

===== Accreting Pulsar Sample =====

== Cyclotron Lines in Transient Pulsars I: Probing the B-Field ==
approved 2012, Suzaku AO-7, NASA\\
PI: K. Pottschmidt\\
CoIs: D.M. Marcu, R.E. Rothschild, S. Suchy, E. Rivers,J. Wilms, I. Kreykenbohm, F. Fuerst, F. W. Schwarm, L. Barragan, M. Kuehnel, S. Mueller,
I. Caballero, P. Kretschmar, V.A. McBride, G. Schoenherr, M. DeCesar, R. Staubert, A. Santangelo, D. Klochkov, V. Doroshenko,
M. Finger, C. Wilson-Hodge, A. Camero, P.A. Jenke, C. Ferrigno, K. Makishima, T. Mihara, M. Nakajima, T. Enoto, Y. Terada, W., Iwakiri

== Cyclotron Lines in Transient Pulsars II: New Lines ==
approved 2012, Suzaku AO-7, NASA\\
PI: K. Pottschmidt\\
CoIs: D.M. Marcu, R.E. Rothschild, S. Suchy, E. Rivers,J. Wilms, I. Kreykenbohm, F. Fuerst, F. W. Schwarm, L. Barragan, M. Kuehnel, S. Mueller,
I. Caballero, P. Kretschmar, V.A. McBride, G. Schoenherr, M. DeCesar, R. Staubert, A. Santangelo, D. Klochkov, V. Doroshenko,
M. Finger, C. Wilson-Hodge, A. Camero, P.A. Jenke, C. Ferrigno, K. Makishima, T. Mihara, M. Nakajima, T. Enoto, Y. Terada, W., Iwakiri

===== Cyg X-1 =====

== Cygnus X-1 with Suzaku and NuSTAR ==
approved 2012, Suzaku AO-7, NASA\\
PI: J.A. Tomsick\\
CoIs: M. Nowak, J. Wilms, K. Pottschmidt, D. Stern, F. Harrison

===== Her X-1 =====

== Study of the cyclotron line feature in Hercules X-1: the profile and the short- and long-term variability ==
approved 2012, Suzaku AO-7, NASA\\
PI: B.W. Grefenstette\\
CoIs: J. Tomsick, D. Klochkov, D. Stern, J. Wilms, R. Staubert

===== 4U 1538-522 =====

== The First Suzaku Observation of 4U1538-522 ==
approved 2012, Suzaku AO-7, NASA\\
PI: R.E Rothschild\\
CoIs: K. Pottschmidt, J. Wilms, P. Hemphill, I. Kreykenbohm, J.J. Rodes, J.M. Torrejon, S. Martinez-Nunez, D. Marcu

===== GX 1+4 =====

== TOO OBSERVATIONS OF THE NEUTRON STAR / RED GIANT BINARY SYSTEM GX 1+4 / V2116 OPHIUCHI ==
approved 2012, Suzaku AO-7, NASA\\
PI: A.Camero-Arranz\\
CoIs: M. Finger, C. Wilson-Hogde, P.Jenke1, K. Pottschmidt, R.E. Rothschild, P. Kretschmar, I. Caballero, F. Fuerst, I. Kreykenbohm, A. González-Galán

===== EXO 2030+375 =====

== Unveil accretion onset during a normal outburst of EXO2030+375 ==
approved 2012, Suzaku AO-7, ESA\\
PI: C. Ferrigno\\
CoIs: D. Klochkov, V. Doroshenko, A. Santangelo, C. Wilson-Hodge, A. Camero-Arranz, J. Wilms, I. Kreykenbohm, K. Pottschmidt, E. Bozzo

===== X Per =====

== Cracking the X Per ==
approved 2012, Suzaku AO-7, ESA\\
PI: V. Doroshenko\\
CoIs: A. Santangelo, I. Kreykenbohm, D. Klochkov, J. Wilms, K. Pottschmidt, R.E. Rothschild

===== 4U 1909+07 =====
== Studying accretion in the uncelebrated HMXB 4U 1909+07 ==
approved 2010, Suzaku AO-5, ESA\\
PI: F. Fuerst\\
CoIs: I. Kreykenbohm, J. Wilms, L. Barragan, K. Pottschmidt, R. E. Rothschild, S. Suchy,
M. Hanke, P. Kretschmar, F. Schwarm, Y. Terada, and the MAGNET collaboration
//observation performed 2010-11-03 for ~20ksec, spectrum shown in Felix' PhD thesis//

Proposals: Chandra

2018-05-02T11:15:21Z

Falkner: Created page with "====== Chandra Proposals ====== ===== Cyg X-1 ===== == Filling the gap in understanding the wind structure of HDE 226868 / Cyg X-1: == PI: Katja Pottschmidt, CoIs: I. Miško..."

====== Chandra Proposals ======

===== Cyg X-1 =====

== Filling the gap in understanding the wind structure of HDE 226868 / Cyg X-1: ==
PI: Katja Pottschmidt, CoIs: I. Miškovičová, D. Marcu, N. Hell, J. Wilms, M.A. Nowak, L. Oskinova, N. Schulz, F. Fürst, J. Torrejón, M. Hanke, G.V. Brown, D. Liedahl, Chandra Guest Observer Program, AO15\\
//The observation was proposed and won as a ToO but has not been triggered yet as of the time of submission. The present proposal is a resubmission of this and the AO12--AO14 programs.//

== Filling the gap in understanding the wind structure of HDE 226868 / Cyg X-1: ==
PI: Katja Pottschmidt, CoIs: I. Miškovičová, D. Marcu, N. Hell, J. Wilms, M.A. Nowak, L. Oskinova, N. Schulz, F. Fürst, J. Torrejón, M. Hanke, G.V. Brown, D. Liedahl, Chandra Guest Observer Program, AO14\\
//The observation was proposed and won as a ToO but has not been triggered yet as of the time of submission.//

== Filling the gap in understanding the wind structure of HDE 226868 / Cyg X-1: ==
PI: Katja Pottschmidt, CoIs: I. Miškovičová, D. Marcu, M. Hanke, J. Wilms, M.A. Nowak, L. Oskinova, N. Schulz, F. Fürst, J. Torrejón, Chandra Guest Observer Program, AO13\\
//The observation was scheduled in 2012 February but we canceled it before execution due to Cyg X-1 returning
into the soft state and it was not rescheduled.//

== Filling the gap in understanding the wind structure of HDE 226868 / Cyg X-1: ==
PI: Manfred Hanke, CoIs: J. Wilms, M.A. Nowak, L. Oskinova, K. Pottschmidt, N. Schulz,
F. Fürst, I. Miškovičová, J. Torrejón, Chandra Guest Observer Program, AO12\\
//Observation performed in 2011 January and analyzed. No lines were seen due to strong photoionization caused by the soft state of the source. Data given to the group of J. McClintock in order to help with the determination of the spin of the black hole from continuum fitting (paper submitted).//

== Cygnus X-1 Viewed Outside of its Secondary's Wind ==
PI: Michael A. Nowak, Co-Is: J. Wilms, M. Hanke, N. Schulz, C. Canizares, K. Pottschmidt, Chandra Guest Observer Program, AO11 \\
//Observation performed, analyzed, results presented at several conferences and in the PhD thesis of M. Hanke, refereed publication submitted (jointly with the AO9 observation).//

== Joint XMM-Newton/Chandra/RXTE Observations of Dips in Cyg X-1 ==
PI: Jörn Wilms, CoIs: M.A. Nowak, K. Pottschmidt, S. Fritz, N.S. Schulz, C. Canizares, Chandra Guest Observer Program, AO9\\
//Observation performed, analyzed, results presented at several conferences and in the PhD thesis of M. Hanke, refereed publication submitted (jointly with the AO11 observation).//

== Cygnus X-1: X-ray Dips, Distance, and the X-ray Halo ==
PI: Katja Pottschmidt, CoIs: J. Wilms, M.A. Nowak, J. Lee, R. McCray, A. Allen, P. Predehl, Chandra Guest Observer Program, 2002--2003.\\
//Observation performed, two papers published, multiple conference posters.//

===== 3A 1954+319 =====
== Exposing the symbiosis of 3A 1954+319 ==
PI: K. Pottschmidt\\
CoIs: D. M. Marcu, F. Fuerst, R. H. D. Corbet, C. B. Markwardt, J. Wilms, S. Suchy, R. E. Rothschild, and the MAGNET collaboration
//observation performed in a very low luminosity state//

Proposals

2018-05-02T11:14:30Z

Falkner: Created page with " Category:Internal ===== Proposals ===== When submitting a proposal for satellite observations usually a list of earlier __successful__ proposals is added. To make it e..."

[[Category:Internal]]

===== Proposals =====

When submitting a proposal for satellite observations usually a list of earlier __successful__ proposals is added. To make it easier to find those relevant for the individual observatory and the scientific aim of the current proposal, here is a list of proposals our group was involved in. Please include title, PI/CoIs, cycle and publications based on the observation.

* [[intern:Proposals:Chandra | Chandra]]
* [[intern:Proposals:Suzaku | Suzaku]]
* [[intern:Proposals:XMM | XMM]]

Orbit plot of GX 301-2 (xfig example)

2018-04-13T16:12:19Z

Falkner: Created page with " Category:SLxfig ====== Orbit plot of GX 301-2 ====== <pre> require("isisscripts") ; % require("/scratch2/fuerst/isisscripts/share/isisscripts") ; % ()=evalfile("/scratc..."

[[Category:SLxfig]]

====== Orbit plot of GX 301-2 ======
<pre>
require("isisscripts") ;
% require("/scratch2/fuerst/isisscripts/share/isisscripts") ;
% ()=evalfile("/scratch2/fuerst/isisscripts/src/plot/xfig_draw_orbit.sl") ;
()=evalfile("xfig_draw_orbit.sl") ;

% orbit parameters of GX 301-2 according to Koh et al (1997).
variable asini = 368.3 ; % lt-sec (semi major axis) 368.3
variable i = 66; % degree to rad 66
variable ecc = 0.46; %
variable omega = 310.4 ; % degrees
variable rstar = 62 ; % solar radii in lt-sec
variable c= 2.998e8 ; %speed of light in cm/s
variable rsun=1.3914e9/2. ; %solar radius in c

%calculating the parameters of the ellipse
%(necessary to get the ASM lightcurve and its coordinate system)
% orbit parameters
i = i /360.*2*PI ; % degree to rad
variable a = asini/sin(i) ;
variable b = sqrt( a^2 * (1-ecc^2) ); % semi-minor axis, lt-sec
omega = omega / 360. * 2 * PI ;

%mean and eccentric anomaly for a full circle / ellipse
variable M = [0:2*PI:#1024] ; %mean anomaly
variable eccanom = KeplerEquation (M, ecc) ;

%widht and height of th plot
variable W = 12 ;
variable H = W ;

%get the orbital phase of the XMM obs and the two Suzaku obs
variable orbphase1 = orbitalphase (55024.103 , "GX301-2_Koh1997") ;
variable orbphase2 = orbitalphase (55024.603 , "GX301-2_Koh1997") ;
variable orbphase_s1b = orbitalphase (54703.552280 , "GX301-2_Koh1997") ;
variable orbphase_s1f= orbitalphase (54704.0036689 , "GX301-2_Koh1997") ;
variable orbphase_s2b = orbitalphase (54836.439629 , "GX301-2_Koh1997") ;
variable orbphase_s2f = orbitalphase (54838.0419328 , "GX301-2_Koh1997") ;

%load the ASM lightcurve and fold into ont the orbit
variable asmlc = RXTE_ASM_lightcurve ("gx301-2") ;
variable asmpp = pfold(asmlc.time, asmlc.rate, 41.498 ;nbins = 80, t0 = 48802.79) ;

%min and max values of the ASM lightcurve
variable max4cord = 5.0 ;
variable min4cord = 0.0 ;
%normalize ASM lightcurve to these values above
variable pplc = (asmpp.value -min4cord)/(max4cord-min4cord);

%%%calculating the ellipse which is modulated by the ASM lightcurve
variable gxx2, gxy2, gxx2b, gxy2b ;
variable gxxall = Double_Type[0];
variable gxyall = Double_Type[0];

%how far away from the orbit should the ASM lightcurve appear?
variable minscal = 0.04 ;
variable maxscal = 0.4 ;
variable j ;
_for j (0, length(pplc)-1, 1) {
(gxx2,gxy2) = ellipse(a+asini*(minscal+maxscal*pplc[j]), b+asini*(minscal+maxscal*pplc[j]),
omega, eccanom[where(asmpp.bin_lo[j] < M/2/PI <= asmpp.bin_hi[j])]) ;
gxxall = [gxxall,gxx2] ;
gxyall = [gxyall,gxy2] ;
}
gxxall = [gxxall,gxxall[0]];
gxyall = [gxyall,gxyall[0]];

%defining eccentric and true anomaly to be used in the coordinate grid
%for the ASM lightcurve (phase indicator lines)
variable eccanom_cgrid = KeplerEquation([0.0:1:0.1]*2*PI, ecc) ;
variable theta_cgrid = 2.*atan( sqrt( (1.+ecc)/(1.-ecc) )*tan(eccanom_cgrid/2.) ); % true anomaly from eccentric anomaly

%all ellipses are shifted so that the star sits at [0,0].
variable shiftx = a*ecc*cos(omega) ;
variable shifty = a*ecc*sin(omega) ;

%define some colors:
variable darkgray = rgb2hex(0.4,0.4,0.4 ; str) ; % used for the phase labels
variable darkestgray = rgb2hex(0.15,0.15,0.15 ; str) ; % used for the ASM labels
variable asmcol = rgb2hex(0.1,0.55,0.1 ; str) ; %

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
%opening the plot and defining the world coordinate system
variable xfgx = xfig_plot_new(W,H) ;
xfgx.world(-720,500,-500, 720) ;
variable axf, bxf ;
%plotting the orbit and highlighting the three observations
(axf, bxf) = xfig_draw_orbit(xfgx, asini, i/2/PI*360., ecc, omega/2/PI*360., rstar ; phases={[orbphase1,orbphase2]},
phase_color=["red", "blue", "blue"], phase_width=4, scale=1.0, trueanom,zerolines=0, width=3, depth=45,
forward_arrow, arrow_width=8,arrow_height=12,arrow_type=2) ;

%Earth indicating arrow
xfgx.plot([0,0],[-145,-425] ; line=1, forward_arrow, depth = 10, arrow_thickness =2) ;
xfgx.xylabel(0,-445, "To Earth" ; size ="\scriptsize"R);
%render now without ASM stuff or phase indicators.
% xfgx.render("plots/gx301_skizze2.eps") ;

%plotting the coodinate grid around the ellipse for the ASM lightcurve
variable xlab, ylab ;
variable labscal = 0.88 ;

% _for j (0,1, 1){
% %phase indicator lines
% xfgx.plot([0,cos(theta_cgrid[j]+omega)]*1500,[0,sin(theta_cgrid[j]+omega)]*1500 ; color="gray", line=3) ;
% %their labels
% (xlab, ylab) = ellipse(a, b, omega, eccanom_cgrid[j]+0.07) ;
% xfgx.xylabel((xlab-shiftx)*labscal, (ylab-shifty)*labscal,sprintf("%.1f", j/10. ) ;
% rotate=(theta_cgrid[j]+omega)/2/PI*360., color=darkgray, size ="\scriptsize"R ) ;
% }
% _for j (2,length(theta_cgrid)-1, 1){
% %phase indicator lines
% xfgx.plot([0,cos(theta_cgrid[j]+omega)]*1500,[0,sin(theta_cgrid[j]+omega)]*1500 ; color="gray", line=3, depth=81) ;
% %their labels
% (xlab, ylab) = ellipse(a, b, omega, eccanom_cgrid[j]-0.07) ;
% xfgx.xylabel((xlab-shiftx)*labscal, (ylab-shifty)*labscal,sprintf("%.1f", j/10. ) ;
% rotate=(theta_cgrid[j]+omega)/2/PI*360.+180, color=darkgray, size ="\scriptsize"R,depth=20 ) ;
% }

%calculating and plotting the grid for the ASM lightcurve
variable xgrid, ygrid,g ;
%how many lines to we want to have?
variable gridpart = 1./4 ;
variable xp1, yp1, xp2, yp2 ;
for (g=0 ; g<=1 ; g+=gridpart)
{
(xgrid, ygrid) = ellipse(a+(minscal+g*maxscal)*asini, b+(minscal+g*maxscal)*asini, omega, [0:2*PI:#256]) ;
(xp1, yp1) = ellipse(a+(minscal+g*maxscal)*asini, b+(minscal+g*maxscal)*asini, omega, 0.3*2*PI) ;
xfgx.plot(xgrid-shiftx, ygrid-shifty; color="gray", line=3,depth=80) ;
%labeling the ASM coordinate system
xfgx.xylabel(xp1-shiftx,yp1-shifty,sprintf("%.2f",g*max4cord+min4cord) ;size ="\scriptsize"R, rotate=-34, depth=7, color=darkestgray ) ;
}
%writing out the units
(xp1, yp1) = ellipse(a+(minscal+g*maxscal)*asini, b+(minscal+g*maxscal)*asini, omega, 0.3*2*PI) ;
xfgx.xylabel(xp1-a*ecc*cos(omega),yp1-a*ecc*sin(omega), "[cts\,sec$^{-1}$]"R ;size ="\scriptsize"R, rotate=-34, depth=7, color=darkestgray ) ;

%plotting the ASM lightcurve
xfgx.plot(gxxall-shiftx,gxyall-shifty; color=asmcol, depth = 10, width=3) ;

%calculating and plotting a bezier curve pointing to the XMM observation
variable bez = bezier({[-600,-350],[-550,-120],[-200,-590],[-142,-250]}, 100) ;
xfgx.plot(bez.x,bez.y ; forward_arrow, color="red", depth = 5) ;
xfgx.xylabel(-600,-375,"XMM" ; size = "\scriptsize"R, color="red") ;
xfgx.xylabel(-600,-400,"47\\,ksec Obs." ; size = "\scriptsize"R, color="red") ;

xfgx.xylabel(0,0,"\\centering Wray\\\\ 977" ; depth=0);

%calculating and plotting a bezier curve pointing to the Suzaku observation
% variable bezs1 = bezier({[-100,280],[-100,460],[-200,100],[-250,505]}, 100) ;
% xfgx.plot(bezs1.x,bezs1.y ; forward_arrow, color="blue", depth=5) ;
% variable bezs1 = bezier({[-80,280],[-30,460],[-50,200],[105,400]}, 100) ;
% xfgx.plot(bezs1.x,bezs1.y ; forward_arrow, color="blue", depth=5) ;
% xfgx.xylabel(-100,260,"Suzaku" ; size = "\scriptsize"R, color="blue") ;

%rendering it all
xfgx.render("plots/gx301_skizze_nosuz2.png") ;

</pre>

Chi-squared parabola (xfig example)

2018-04-13T16:10:40Z

Falkner: Created page with " Category:SLxfig ====== Chi-squared parabola ====== <pre> require("isisscripts"); variable t = [-3:3:.01]; define f(x){ return -x^2+9; } variable L = f(t); variab..."

[[Category:SLxfig]]

====== Chi-squared parabola ======
<pre>
require("isisscripts");

variable t = [-3:3:.01];

define f(x){
return -x^2+9;
}

variable L = f(t);

variable X = xfig_plot_new;
X.plot(t,L);
%X.axis(;off);
X.xlabel(`$\theta$`);
X.ylabel(`$ln L$`);

variable sig1_lo,sig1_hi,sig2_lo,sig2_hi,sig3_lo,sig3_hi;
sig1_lo = min(t[where(L>max(L)-0.5)]);
sig1_hi = max(t[where(L>max(L)-0.5)]);
sig2_lo = min(t[where(L>max(L)-2)]);
sig2_hi = max(t[where(L>max(L)-2)]);
sig3_lo = min(t[where(L>max(L)-4.5)]);
sig3_hi = max(t[where(L>max(L)-4.5)]);

X.yaxis(; major=[max(L),max(L)-0.5,max(L)-2,max(L)-4.5],
ticlabels=[`$L_{\text{max}}$`,
`$L_{\text{max}}$-0.5`,
`$L_{\text{max}}$-2`,
`$L_{\text{max}}$-4.5`],
ticlabels2=0);
X.x1axis(; major=[sig1_lo,sig2_lo,sig3_lo,0,sig1_hi,sig2_hi,sig3_hi],
ticlabels=[`$-1\sigma$`,`$-2\sigma$`,
`$-3\sigma$`,`$\hat{\theta}=0$`,`$+1\sigma$`,
`$+2\sigma$`,`$+3\sigma$`],
);

%X.world1(min_max(t;pad=0.1),min_max(L;pad=0.1));
X.plot([min(t[where(L>max(L)-0.5)]),max(t[where(L>max(L)-0.5)])],[max(L)-0.5,max(L)-0.5];color="black",width=3,
backward_arrow,forward_arrow,arrow_type=2,arrow_style=3,arrow_width=6);
X.plot([min(t[where(L>max(L)-2)]),max(t[where(L>max(L)-2)])],[max(L)-2,max(L)-2];color="black",width=3,
backward_arrow,forward_arrow,arrow_type=1,arrow_style=3,arrow_width=6);
X.plot([min(t[where(L>max(L)-4.5)]),max(t[where(L>max(L)-4.5)])],[max(L)-4.5,max(L)-4.5];color="black",width=3,
backward_arrow,forward_arrow,arrow_type=3,arrow_style=1,arrow_width=6);
% lines 1 sigma
X.plot([min(t[where(L>max(L)-0.5)]),min(t[where(L>max(L)-0.5)])],[0,max(L)-0.5]);
X.plot([max(t[where(L>max(L)-0.5)]),max(t[where(L>max(L)-0.5)])],[0,max(L)-0.5]);
% lines 2 sigma
X.plot([min(t[where(L>max(L)-2)]),min(t[where(L>max(L)-2)])],[0,max(L)-2]);
X.plot([max(t[where(L>max(L)-2)]),max(t[where(L>max(L)-2)])],[0,max(L)-2]);
% lines 3 sigma
X.plot([min(t[where(L>max(L)-4.5)]),min(t[where(L>max(L)-4.5)])],[0,max(L)-4.5]);
X.plot([max(t[where(L>max(L)-4.5)]),max(t[where(L>max(L)-4.5)])],[0,max(L)-4.5]);
% line max, max-0.5, max-2, max-4.5
X.plot([-4,0,],[max(L),max(L)];line=1);
X.plot([-4,0,],[max(L)-0.5,max(L)-0.5];line=1);
X.plot([-4,0,],[max(L)-2,max(L)-2];line=1);
X.plot([-4,0,],[max(L)-4.5,max(L)-4.5];line=1);

variable sig_1 = xfig_new_text(`$68.3\,\%$`);
variable sig_2 = xfig_new_text(`$90\,\%$`);
variable sig_3 = xfig_new_text(`$99\,\%$`);

X.add_object(sig_1,0,max(L)-1);
X.add_object(sig_2,0,max(L)-2.5);
X.add_object(sig_3,0,max(L)-5);

X.render("chisqr_parabola.eps");
</pre>

Chi-squared contours (xfig example)

2018-04-13T16:09:09Z

Falkner: Created page with " Category:SLxfig ====== Chi-squared contours ====== <pre> require("isisscripts"); variable x = [0:2.2*PI:.01]; variable alpha,x0,y0,range,col; variable a,b; variable..."

[[Category:SLxfig]]

====== Chi-squared contours ======

<pre>
require("isisscripts");

variable x = [0:2.2*PI:.01];

variable alpha,x0,y0,range,col;
variable a,b;
variable i;

define ellipse(x,i,alpha,x0,y0,a,b){
return {x0+cos(alpha)*(a-i*1.5)*cos(x)-sin(alpha)*(b-i*1.5)*sin(x),
y0+sin(alpha)*(a-i*1.5)*cos(x)+cos(alpha)*(b-i*1.5)*sin(x)};
}

x0=0.;y0=0.;
alpha=PI/4.;
range=10.;
a=4.,b=10.;

col=["red","blue","green4"];

vmessage("sin(alpha)=%f, cos(alpha=%f)",sin(alpha),cos(alpha));

variable X=xfig_plot_new(15,15);
X.world([-range:range],[-range:range]);

variable x_par_0 = ellipse(x,0,alpha,x0,y0,a,b)[0];
variable y_par_0 = ellipse(x,0,alpha,x0,y0,a,b)[1];
variable x_par_1 = ellipse(x,1,alpha,x0,y0,a,b)[0];
variable y_par_1 = ellipse(x,1,alpha,x0,y0,a,b)[1];
variable x_par_2 = ellipse(x,2,alpha,x0,y0,a,b)[0];
variable y_par_2 = ellipse(x,2,alpha,x0,y0,a,b)[1];

X.yaxis(; major=[min(y_par_0),min(y_par_1),min(y_par_2),0,max(y_par_0),max(y_par_1),max(y_par_2)],
ticlabels=[`$3\sigma$`,
`$-2\sigma$`,
`$-1\sigma$`,
`$\hat{b}=0$`,
`$+3\sigma$`,
`$+2\sigma$`,
`$+1\sigma$`],
ticlabels2=0);
X.y2axis(;off);
X.x2axis(;off);
X.x1axis(; major=[min(x_par_0),min(x_par_1),min(x_par_2),0,max(x_par_0),max(x_par_1),max(x_par_2)],
ticlabels=[`$-3\sigma$`,
`$-2\sigma$`,
`$-1\sigma$`,
`$\hat{a}=0$`,
`$+3\sigma$`,
`$+2\sigma$`,
`$+1\sigma$`],
ticlabels2=0);

X.plot([-10,0],[0,0];line=1);
X.plot([0,0],[-10,0];line=1);

for (i=0;i<3;i++){
variable x_par = ellipse(x,i,alpha,x0,y0,a,b)[0];
variable y_par = ellipse(x,i,alpha,x0,y0,a,b)[1];

X.plot(x_par,y_par;color=col[i],width=3);
X.plot([min(x_par),min(x_par)],[-10,y_par[where(x_par==min(x_par))]];line=1,color=col[i]);
X.plot([max(x_par),max(x_par)],[-10,y_par[where(x_par==max(x_par))]];line=1,color=col[i]);
X.plot([-10,x_par[where(y_par==min(y_par))]],[min(y_par),min(y_par)];line=1,color=col[i]);
X.plot([-10,x_par[where(y_par==max(y_par))]],[max(y_par),max(y_par)];line=1,color=col[i]);
X.plot([0,0],[0,0];sym="x");
}

X.render("chisqr_ellipse.eps");
</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:06:56Z

Falkner: /* Peridodic Table of the Elements */

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

[[File:periodictable_nh.png|1000px]]

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:06:35Z

Falkner: /* Peridodic Table of the Elements */

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

[[File:periodictable_nh.png|500x]]

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:06:24Z

Falkner: /* Peridodic Table of the Elements */

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

[[File:periodictable_nh.png|1000x]]

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:06:17Z

Falkner: /* Peridodic Table of the Elements */

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

[[File:periodictable_nh.png|1000]]

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:06:03Z

Falkner: /* Peridodic Table of the Elements */

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

[[File:periodictable_nh.png]]

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:05:47Z

Falkner: /* Peridodic Table of the Elements */

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

[[File:periodictable_nh.png?nolink&1000]]

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Periodic Table of the Elements (xfig example)

2018-04-13T16:05:11Z

Falkner: Created page with " Category:SLxfig ==== Peridodic Table of the Elements ==== {{:isis:slxfig:periodictable_nh.png?nolink&1000|}} This script loads a datafile (ptdata.sl) with the symbols..."

[[Category:SLxfig]]

==== Peridodic Table of the Elements ====

{{:isis:slxfig:periodictable_nh.png?nolink&1000|}}

This script loads a datafile (ptdata.sl) with the symbols, names, and other properties in form of an array of lists. The array is built such that the index of an entry equals the proton number //Z//. Some of the properties (e.g., groups like metals / nonmetals, state, ...) are used as names for color definitions.

The data file format is as:
<pre>
variable elements = [{"Symbol","Element Name","\centering \small Atomic\\Weight"R,"black","white"},
{"H", "Hydrogen", "1.01" ,"gas", "alkali"},
{"He", "Helium", "4.00" ,"gas", "noblegas"},
...
];
</pre>

Additionally, the code is written in a way that allows to easily switch between color and black/white versions of the table via a simple integer variable ''incolor'', which triggers different color definitions.

<pre>
require("isisscripts");

% periodic table entries:
require("ptdata.sl");
%-> variable elements

variable incolor = 1;
if(incolor){
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("gas" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0.4 , 0.4 , 0.4 ) );
xfig_new_color("alkali" , rgb2hex(1. , 0. , 0. ) );
xfig_new_color("alkaline" , rgb2hex(0. , 1. , 0. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 0. ) );
xfig_new_color("metal" , rgb2hex(1. , 0.5 , 0. ) );
xfig_new_color("nonmetal" , rgb2hex(0. , .4 , 0.8 ) );
xfig_new_color("noblegas" , rgb2hex(0. , 0. , 1. ) );
xfig_new_color("unknown" , rgb2hex(0.75 , 0.75 , .75 ) );
}else{
xfig_new_color("solid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("liquid" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("gas" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("synthetic" , rgb2hex(0. , 0. , 0. ) );
xfig_new_color("alkali" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("alkaline" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("transmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("metal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("nonmetal" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("noblegas" , rgb2hex(1. , 1. , 1. ) );
xfig_new_color("unknown" , rgb2hex(1. , 1. , 1. ) );
}

% a map for the positions of the individual elements (Z) in the periodic table
variable rows = {
[1,0*ones(17),2],
[3,4,0*ones(11),[5:10]],
[11,12,0*ones(11),[13:18]],
[19,20,0,[21:36]],
[37,38,0,[39:54]],
[55,56,0,[71:86]],
[87,88,0,[103:118]],
[57:70],
[89:102]
};

% the dimensions for the cells of the table
variable wcell = 2.15;
variable hcell = 2.5;

define cell(Z){
variable xc = xfig_plot_new(wcell,hcell);
xc.world(0,wcell,0,hcell);
xc.axes(;off);
xc.plot([0,0,wcell,wcell,0],[0,hcell,hcell,0,0];depth=15);
variable bg = xfig_new_polyline([0,0,wcell,wcell],[0,hcell,hcell,0];
closed,areafill=30,fillcolor=elements[Z][4],color="white");
xc.add_object(bg,0,0,-0.5,-0.5;depth=50);
xc.xylabel(0.5*wcell,0.1*hcell, sprintf("\usefont{T1}{arial}{b}{n}\scriptsize %s"R,elements[Z][1]),0,0;depth=10);
xc.xylabel(0.5*wcell,0.58*hcell, elements[Z][0] ,0,0; size="\Large"R, color=elements[Z][3],depth=10);
if(Z==0){
xc.xylabel(0.05*wcell,0.87*hcell, "Z" ,-0.5,0 ;depth=10);
}else{
xc.xylabel(0.05*wcell,0.87*hcell, "$Z"$ ,-0.5,0 ;depth=10);
}
xc.xylabel(0.5*wcell,0.31*hcell, elements[Z][2] ,0,0 ;depth=10);
return xc;
}

variable pl = xfig_plot_new(19*wcell,10*hcell);
pl.world(1,19*wcell,10*hcell,1);
pl.axes(;off);

variable ii, rr,cc;
_for ii(1,7,1){
cc=1;
foreach rr (rows[ii-1]){
if(rr==0){ cc++; continue; }
pl.add_object(cell(rr), cc*wcell, ii*hcell, -0.5,-0.5);
cc++;
}
}

_for ii(8,9,1){
cc=4;
foreach rr (rows[ii-1]){
pl.add_object(cell(rr), cc*wcell, (ii+1)*hcell, -0.5,-0.5);
cc++;
}
}

pl.xylabel(3.5*wcell,5.5*hcell,"57--70"R,0,0);
pl.xylabel(3.5*wcell,5.7*hcell,"\large *"R,0,0);
pl.xylabel(3.5*wcell,6.5*hcell,"89--102"R,0,0);
pl.xylabel(3.5*wcell,6.7*hcell,"\large **"R,0,0);

pl.xylabel(3.8*wcell,8.5*hcell,"{\large *}\,Lanthanoids"R,0.5,0);
pl.xylabel(3.8*wcell,9.5*hcell,"{\large **}\,Actinoids"R,0.5,0);

% cell description:
variable sc = 1.2;
wcell *=sc; hcell *=sc;
pl.add_object(cell(0),5.4*wcell/sc,2.1*hcell/sc,-0.5,-0.5);
wcell /=sc; hcell /=sc;

if(incolor){
pl.xylabel((5.35+sc)*wcell,1.4*hcell, "\Huge\{"R,-0.5,0);
pl.xylabel((5.65+sc)*wcell,1.125*hcell, "solid", -0.5,0;color="solid");
pl.xylabel((5.65+sc)*wcell,1.325*hcell, "liquid", -0.5,0;color="liquid");
pl.xylabel((5.65+sc)*wcell,1.525*hcell, "gas", -0.5,0;color="gas");
pl.xylabel((5.65+sc)*wcell,1.65*hcell, "synthetic", -0.5,0;color="synthetic");

variable desc = ["alkali metals", "alkaline earth metals", "transitional metals",
"other metals", "nonmetals", "noble gases", "unknown group"];
variable cl = ["alkali","alkaline","transmetal","metal","nonmetal",
"noblegas","unknown"];
variable bx;

_for ii(0,6,1){
bx = xfig_new_polyline([0,0,wcell/2.,wcell/2.],[0,hcell/4.,hcell/4.,0];
closed,areafill=30,fillcolor=cl[ii],color="white");
pl.add_object(bx,9.75*wcell,(0.75+ii*0.26)*hcell,0.5,0;depth=50);
pl.xylabel(9.85*wcell,(0.75+ii*0.26)*hcell,desc[ii],-0.5,0);
}

}

pl.render("pt.pdf");

</pre>

Plotting Model Components (xfig example)

2018-04-13T16:01:28Z

Falkner:

[[Category:SLxfig]]

====== Plotting Model Components ======

Even in ISIS, plotting individual model components is somewhat tricky. Once upon a time, there was a useful description to be found [[http://www.black-hole.eu/media/summerschool2/X-ray_Spectra_Part_II.pdf|here]] on pp. 39-43. Unfortunately, since then, the fancy plot commands changed but not the presentation. I will try here to give you a short overview of how to do this at the moment.

===== The plot in ISIS =====

=== Loading the data ===
<pre>

load_data("ridge_epoch_4_2.pha"); % load your spectra
flux_corr(1); % the usual setup things
group(1; min_sn=5, min_chan=1, bounds=3, unit="keV");
ignore_en(1, ,3);
ignore_en(1,20, );

</pre>

=== Starting the plot ===
<pre>
open_plot; % open a plot window...
multiplot([5,2]); % ...to get the two panels for spectrum and residuals
mpane(1); % start plotting in the upper panel
xlog; ylog; % if that's the usual way you want to have your spectra

load_par("2011-03-09_15:00:21_best.par"); % load the best-fit parameters for the whole model
() = eval_counts; % evaluate the model

plot_data(1; dsym=-4, dcol=1, decol=1, mcol=15, xrange={2.9,20.5}, yrange={0.005,1.5}, no_reset=1);
% the command to plot your spectrum. note that the
% "no_reset" command now has to be used as a qualifier,
% different from Mike's presentation. also, the commands
% for setting the ranges have changed (from "xrng" to "xrange").

write_plot("write_data_allin"); % dump the spectrum and model into an ASCII file. for a data
% plot it will end in .dat

</pre>

and now a single component of the model:

<pre>
ff = get_fit_fun(); % save the complete model

fit_fun("gauss(1)"); % now only the component you are interested in
% **ATTENTION**: next time, reload the best-fit parameters before
% adding another model component: "fit_fun(ff)"!!!
() = eval_counts; % evaluate the model (= the single component)
plot_data(1; dsym=-4, dcol=1, decol=1, mcol=4, oplt=1, no_reset=1); % plot the component in a different color
write_plot("write_data_gauss"); % again, dump the plot into another ASCII file

</pre>

Of course the colors are unimportant for the moment, as you will have to do them again when doing the plot in xfig, but now the example is complete if you don't have the time for xfig and just want a quick way to plot.

=== The second panel ===

<pre>
mpane(2); % move to the lower panel
ylin; % for residuals, a linear y scale

load_par("2011-03-09_15:00:21_best.par"); % again, reload the best-fit parameters
() = eval_counts; % evaluate the model

plot_residuals(1; res=1, rsym=0, rcol=15, recol=1, xrange={2.9,20.5}, yrange={-3,20}, no_reset=1);
% plot only the residuals. qualifier are the
% same as for "plot_data"
write_plot("write_res_allin"); % dump the plot into an ASCII file. for residual plots
% it will end in .res
</pre>

Now the residuals you get when the component is missing in the model (in this case: an iron line complex consisting of three gaussian lines):

<pre>
set_par([5,8,11],0); % set the norm of your component to zero
() = eval_counts; % evaluate the model
plot_residuals(1; res=1, rsym=0, rcol=4, recol=1, oplt=1, no_reset=1);
% plot the residuals
write_plot("write_res_gauss.dat"); % dump the plot into an ASCII file
close_plot; % close the plot window when finished
</pre>

Now you have all the files you need to make a nice xfig plot.

===== Reading and plotting the data with xfig =====

=== Initializing the plot ===
<pre>
require("xfig"); % you need this to do the nice plots

variable p = Struct_Type[2]; % this way, you can do a multiplot
p[0] = xfig_plot_new(12.6,6.4); % the size of the upper panel
p[1] = xfig_plot_new(12.6,2.6); % size of the lower panel. proportions are the same as in ISIS plots.

p[0].world (2.9,20.5,0.005,1.5; xlog=1, ylog=1); % ranges and scales of the upper panel
p[0].x1axis(;ticlabels=0); % no ticlabels between the two panels
p[0].y2axis(;ticlabels=0); % no ticlabels at the right side of the plot
p[0].ylabel("Counts/s/keV"); % ylabel of the upper panel

p[1].world (2.9,20.5,-5,41; xlog=1); % ranges and scales of the lower panel
p[1].ylabel(`$\chi$`); % ylabel of the lower panel
p[1].xlabel("Energy [keV]"); % xlabel of the whole plot
</pre>

=== Read and plot the ASCII files: upper panel ===
<pre>
variable lo, hi, y, e, m; % five variables for the five columns of the .dat file
(lo, hi, y, e, m) = readcol ("write_data_allin_0.dat",1,2,3,4,5); % read all five columns: bin_lo, bin_hi, value, error, model
variable h = [lo, hi[-1]]; % add the upper limit of the last bin to the lower limits of the bins
variable x = (lo+hi)/2; % an array with the central points of each bin

p[0].hplot (h,m; color="blue4"); % make a histogram plot of the model
% see "p[0].hplot(;help)" for help:
% if length(h)=length(m)+1, the last element is taken as the upper
% limit of the last bin. if both arguments have the same length,
% the last bin will have the same size as the one before.
% **ATTENTION** at the moment, hplot can only plot a histogram with
% errorbars when all arrays have the same length!!!

p[0].plot (x,y,e; sym="point"); % a normal plot to get data points and errorbars
</pre>

now the model component:

<pre>
variable lo, hi, y, e, m; % clear the variables the previous variables
(lo, hi, y, e, m) = readcol ("write_data_gauss_0.dat",1,2,3,4,5); % read the ASCII file with the component
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting

p[0].hplot (h,m; color="green4"); % plot a histogram of the model component
</pre>

=== Read and plot the ASCII files: lower panel ===
<pre>
variable lo, hi, y, el, eu; % five variables for the columns of the .res file
(lo, hi, y, el, eu) = readcol ("write_res_allin_0.res",1,2,3,4,5); % read all five columns: bin_lo, bin_hi, value, err_lo, err_hi
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting
variable x = (lo+hi)/2; % the center of each bin

p[1].hplot (h,y; color="blue4"); % plot the residuals histogram of the whole model

p[1].plot (x,y, ,{y-el,eu-y}; sym="point", color="blue4", size=-1, eb_factor=0);
% optional: add the errorbars, if you did the plot with an
% individually sized last bin (i.e., added the last bin_hi, see above)
</pre>

now the residuals without the component:

<pre>
variable lo, hi, y, el, eu; % clear the variables
(lo, hi, y, el, eu) = readcol ("write_res_gauss_0.res",1,2,3,4,5); % read the .res file
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting
variable x = (lo+hi)/2; % center of each bin

p[1].hplot (h,y; color="green4"); % plot the residuals histogram without the component\

p[1].plot (x,y, ,{y-el,eu-y}; sym="point", color="green4", size=-1, eb_factor=0);
% optional: add the errorbars
</pre>

=== Render the plot ===
<pre>
xfig_multiplot (p).render ("plot_components.png");
</pre>

===== The Result =====
[[File:plot_individual_components.png]]

To remove the annoying vertical lines at the first and last bin of the histogram, you can use the (undocumented) qualifiers ''y_last'' and ''y_first'' of ''hplot''. If you set them to the value of the first or the last bin respectively, the vertical lines vanish:
p[0].hplot (h,m; color="blue4", y_first=m[0], y_last=m[-1]);

Plotting Model Components (xfig example)

2018-04-13T15:58:43Z

Falkner: /* The Result */

[[Category:SLxfig]]

====== Plotting Model Components ======

Even in ISIS, plotting individual model components is somewhat tricky. Once upon a time, there was a useful description to be found [[http://www.black-hole.eu/media/summerschool2/X-ray_Spectra_Part_II.pdf|here]] on pp. 39-43. Unfortunately, since then, the fancy plot commands changed but not the presentation. I will try here to give you a short overview of how to do this at the moment.

===== The plot in ISIS =====

=== Loading the data ===
<code>

load_data("ridge_epoch_4_2.pha"); % load your spectra
flux_corr(1); % the usual setup things
group(1; min_sn=5, min_chan=1, bounds=3, unit="keV");
ignore_en(1, ,3);
ignore_en(1,20, );

</code>

=== Starting the plot ===
<code>
open_plot; % open a plot window...
multiplot([5,2]); % ...to get the two panels for spectrum and residuals
mpane(1); % start plotting in the upper panel
xlog; ylog; % if that's the usual way you want to have your spectra

load_par("2011-03-09_15:00:21_best.par"); % load the best-fit parameters for the whole model
() = eval_counts; % evaluate the model

plot_data(1; dsym=-4, dcol=1, decol=1, mcol=15, xrange={2.9,20.5}, yrange={0.005,1.5}, no_reset=1);
% the command to plot your spectrum. note that the
% "no_reset" command now has to be used as a qualifier,
% different from Mike's presentation. also, the commands
% for setting the ranges have changed (from "xrng" to "xrange").

write_plot("write_data_allin"); % dump the spectrum and model into an ASCII file. for a data
% plot it will end in .dat

</code>

and now a single component of the model:

<code>
ff = get_fit_fun(); % save the complete model

fit_fun("gauss(1)"); % now only the component you are interested in
% **ATTENTION**: next time, reload the best-fit parameters before
% adding another model component: "fit_fun(ff)"!!!
() = eval_counts; % evaluate the model (= the single component)
plot_data(1; dsym=-4, dcol=1, decol=1, mcol=4, oplt=1, no_reset=1); % plot the component in a different color
write_plot("write_data_gauss"); % again, dump the plot into another ASCII file

</code>

Of course the colors are unimportant for the moment, as you will have to do them again when doing the plot in xfig, but now the example is complete if you don't have the time for xfig and just want a quick way to plot.

=== The second panel ===

<code>
mpane(2); % move to the lower panel
ylin; % for residuals, a linear y scale

load_par("2011-03-09_15:00:21_best.par"); % again, reload the best-fit parameters
() = eval_counts; % evaluate the model

plot_residuals(1; res=1, rsym=0, rcol=15, recol=1, xrange={2.9,20.5}, yrange={-3,20}, no_reset=1);
% plot only the residuals. qualifier are the
% same as for "plot_data"
write_plot("write_res_allin"); % dump the plot into an ASCII file. for residual plots
% it will end in .res
</code>

Now the residuals you get when the component is missing in the model (in this case: an iron line complex consisting of three gaussian lines):

<code>
set_par([5,8,11],0); % set the norm of your component to zero
() = eval_counts; % evaluate the model
plot_residuals(1; res=1, rsym=0, rcol=4, recol=1, oplt=1, no_reset=1);
% plot the residuals
write_plot("write_res_gauss.dat"); % dump the plot into an ASCII file
close_plot; % close the plot window when finished
</code>

Now you have all the files you need to make a nice xfig plot.

===== Reading and plotting the data with xfig =====

=== Initializing the plot ===
<code>
require("xfig"); % you need this to do the nice plots

variable p = Struct_Type[2]; % this way, you can do a multiplot
p[0] = xfig_plot_new(12.6,6.4); % the size of the upper panel
p[1] = xfig_plot_new(12.6,2.6); % size of the lower panel. proportions are the same as in ISIS plots.

p[0].world (2.9,20.5,0.005,1.5; xlog=1, ylog=1); % ranges and scales of the upper panel
p[0].x1axis(;ticlabels=0); % no ticlabels between the two panels
p[0].y2axis(;ticlabels=0); % no ticlabels at the right side of the plot
p[0].ylabel("Counts/s/keV"); % ylabel of the upper panel

p[1].world (2.9,20.5,-5,41; xlog=1); % ranges and scales of the lower panel
p[1].ylabel(`$\chi$`); % ylabel of the lower panel
p[1].xlabel("Energy [keV]"); % xlabel of the whole plot
</code>

=== Read and plot the ASCII files: upper panel ===
<code>
variable lo, hi, y, e, m; % five variables for the five columns of the .dat file
(lo, hi, y, e, m) = readcol ("write_data_allin_0.dat",1,2,3,4,5); % read all five columns: bin_lo, bin_hi, value, error, model
variable h = [lo, hi[-1]]; % add the upper limit of the last bin to the lower limits of the bins
variable x = (lo+hi)/2; % an array with the central points of each bin

p[0].hplot (h,m; color="blue4"); % make a histogram plot of the model
% see "p[0].hplot(;help)" for help:
% if length(h)=length(m)+1, the last element is taken as the upper
% limit of the last bin. if both arguments have the same length,
% the last bin will have the same size as the one before.
% **ATTENTION** at the moment, hplot can only plot a histogram with
% errorbars when all arrays have the same length!!!

p[0].plot (x,y,e; sym="point"); % a normal plot to get data points and errorbars
</code>

now the model component:

<code>
variable lo, hi, y, e, m; % clear the variables the previous variables
(lo, hi, y, e, m) = readcol ("write_data_gauss_0.dat",1,2,3,4,5); % read the ASCII file with the component
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting

p[0].hplot (h,m; color="green4"); % plot a histogram of the model component
</code>

=== Read and plot the ASCII files: lower panel ===
<code>
variable lo, hi, y, el, eu; % five variables for the columns of the .res file
(lo, hi, y, el, eu) = readcol ("write_res_allin_0.res",1,2,3,4,5); % read all five columns: bin_lo, bin_hi, value, err_lo, err_hi
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting
variable x = (lo+hi)/2; % the center of each bin

p[1].hplot (h,y; color="blue4"); % plot the residuals histogram of the whole model

p[1].plot (x,y, ,{y-el,eu-y}; sym="point", color="blue4", size=-1, eb_factor=0);
% optional: add the errorbars, if you did the plot with an
% individually sized last bin (i.e., added the last bin_hi, see above)
</code>

now the residuals without the component:

<code>
variable lo, hi, y, el, eu; % clear the variables
(lo, hi, y, el, eu) = readcol ("write_res_gauss_0.res",1,2,3,4,5); % read the .res file
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting
variable x = (lo+hi)/2; % center of each bin

p[1].hplot (h,y; color="green4"); % plot the residuals histogram without the component\

p[1].plot (x,y, ,{y-el,eu-y}; sym="point", color="green4", size=-1, eb_factor=0);
% optional: add the errorbars
</code>

=== Render the plot ===
<code>
xfig_multiplot (p).render ("plot_components.png");
</code>

===== The Result =====
[[File:plot_individual_components.png]]

To remove the annoying vertical lines at the first and last bin of the histogram, you can use the (undocumented) qualifiers ''y_last'' and ''y_first'' of ''hplot''. If you set them to the value of the first or the last bin respectively, the vertical lines vanish:
p[0].hplot (h,m; color="blue4", y_first=m[0], y_last=m[-1]);

Plotting Model Components (xfig example)

2018-04-13T15:57:38Z

Falkner: Created page with " Category:SLxfig ====== Plotting Model Components ====== Even in ISIS, plotting individual model components is somewhat tricky. Once upon a time, there was a useful des..."

[[Category:SLxfig]]

====== Plotting Model Components ======

Even in ISIS, plotting individual model components is somewhat tricky. Once upon a time, there was a useful description to be found [[http://www.black-hole.eu/media/summerschool2/X-ray_Spectra_Part_II.pdf|here]] on pp. 39-43. Unfortunately, since then, the fancy plot commands changed but not the presentation. I will try here to give you a short overview of how to do this at the moment.

===== The plot in ISIS =====

=== Loading the data ===
<code>

load_data("ridge_epoch_4_2.pha"); % load your spectra
flux_corr(1); % the usual setup things
group(1; min_sn=5, min_chan=1, bounds=3, unit="keV");
ignore_en(1, ,3);
ignore_en(1,20, );

</code>

=== Starting the plot ===
<code>
open_plot; % open a plot window...
multiplot([5,2]); % ...to get the two panels for spectrum and residuals
mpane(1); % start plotting in the upper panel
xlog; ylog; % if that's the usual way you want to have your spectra

load_par("2011-03-09_15:00:21_best.par"); % load the best-fit parameters for the whole model
() = eval_counts; % evaluate the model

plot_data(1; dsym=-4, dcol=1, decol=1, mcol=15, xrange={2.9,20.5}, yrange={0.005,1.5}, no_reset=1);
% the command to plot your spectrum. note that the
% "no_reset" command now has to be used as a qualifier,
% different from Mike's presentation. also, the commands
% for setting the ranges have changed (from "xrng" to "xrange").

write_plot("write_data_allin"); % dump the spectrum and model into an ASCII file. for a data
% plot it will end in .dat

</code>

and now a single component of the model:

<code>
ff = get_fit_fun(); % save the complete model

fit_fun("gauss(1)"); % now only the component you are interested in
% **ATTENTION**: next time, reload the best-fit parameters before
% adding another model component: "fit_fun(ff)"!!!
() = eval_counts; % evaluate the model (= the single component)
plot_data(1; dsym=-4, dcol=1, decol=1, mcol=4, oplt=1, no_reset=1); % plot the component in a different color
write_plot("write_data_gauss"); % again, dump the plot into another ASCII file

</code>

Of course the colors are unimportant for the moment, as you will have to do them again when doing the plot in xfig, but now the example is complete if you don't have the time for xfig and just want a quick way to plot.

=== The second panel ===

<code>
mpane(2); % move to the lower panel
ylin; % for residuals, a linear y scale

load_par("2011-03-09_15:00:21_best.par"); % again, reload the best-fit parameters
() = eval_counts; % evaluate the model

plot_residuals(1; res=1, rsym=0, rcol=15, recol=1, xrange={2.9,20.5}, yrange={-3,20}, no_reset=1);
% plot only the residuals. qualifier are the
% same as for "plot_data"
write_plot("write_res_allin"); % dump the plot into an ASCII file. for residual plots
% it will end in .res
</code>

Now the residuals you get when the component is missing in the model (in this case: an iron line complex consisting of three gaussian lines):

<code>
set_par([5,8,11],0); % set the norm of your component to zero
() = eval_counts; % evaluate the model
plot_residuals(1; res=1, rsym=0, rcol=4, recol=1, oplt=1, no_reset=1);
% plot the residuals
write_plot("write_res_gauss.dat"); % dump the plot into an ASCII file
close_plot; % close the plot window when finished
</code>

Now you have all the files you need to make a nice xfig plot.

===== Reading and plotting the data with xfig =====

=== Initializing the plot ===
<code>
require("xfig"); % you need this to do the nice plots

variable p = Struct_Type[2]; % this way, you can do a multiplot
p[0] = xfig_plot_new(12.6,6.4); % the size of the upper panel
p[1] = xfig_plot_new(12.6,2.6); % size of the lower panel. proportions are the same as in ISIS plots.

p[0].world (2.9,20.5,0.005,1.5; xlog=1, ylog=1); % ranges and scales of the upper panel
p[0].x1axis(;ticlabels=0); % no ticlabels between the two panels
p[0].y2axis(;ticlabels=0); % no ticlabels at the right side of the plot
p[0].ylabel("Counts/s/keV"); % ylabel of the upper panel

p[1].world (2.9,20.5,-5,41; xlog=1); % ranges and scales of the lower panel
p[1].ylabel(`$\chi$`); % ylabel of the lower panel
p[1].xlabel("Energy [keV]"); % xlabel of the whole plot
</code>

=== Read and plot the ASCII files: upper panel ===
<code>
variable lo, hi, y, e, m; % five variables for the five columns of the .dat file
(lo, hi, y, e, m) = readcol ("write_data_allin_0.dat",1,2,3,4,5); % read all five columns: bin_lo, bin_hi, value, error, model
variable h = [lo, hi[-1]]; % add the upper limit of the last bin to the lower limits of the bins
variable x = (lo+hi)/2; % an array with the central points of each bin

p[0].hplot (h,m; color="blue4"); % make a histogram plot of the model
% see "p[0].hplot(;help)" for help:
% if length(h)=length(m)+1, the last element is taken as the upper
% limit of the last bin. if both arguments have the same length,
% the last bin will have the same size as the one before.
% **ATTENTION** at the moment, hplot can only plot a histogram with
% errorbars when all arrays have the same length!!!

p[0].plot (x,y,e; sym="point"); % a normal plot to get data points and errorbars
</code>

now the model component:

<code>
variable lo, hi, y, e, m; % clear the variables the previous variables
(lo, hi, y, e, m) = readcol ("write_data_gauss_0.dat",1,2,3,4,5); % read the ASCII file with the component
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting

p[0].hplot (h,m; color="green4"); % plot a histogram of the model component
</code>

=== Read and plot the ASCII files: lower panel ===
<code>
variable lo, hi, y, el, eu; % five variables for the columns of the .res file
(lo, hi, y, el, eu) = readcol ("write_res_allin_0.res",1,2,3,4,5); % read all five columns: bin_lo, bin_hi, value, err_lo, err_hi
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting
variable x = (lo+hi)/2; % the center of each bin

p[1].hplot (h,y; color="blue4"); % plot the residuals histogram of the whole model

p[1].plot (x,y, ,{y-el,eu-y}; sym="point", color="blue4", size=-1, eb_factor=0);
% optional: add the errorbars, if you did the plot with an
% individually sized last bin (i.e., added the last bin_hi, see above)
</code>

now the residuals without the component:

<code>
variable lo, hi, y, el, eu; % clear the variables
(lo, hi, y, el, eu) = readcol ("write_res_gauss_0.res",1,2,3,4,5); % read the .res file
variable h = [lo, hi[-1]]; % add the last bin_hi for histogram plotting
variable x = (lo+hi)/2; % center of each bin

p[1].hplot (h,y; color="green4"); % plot the residuals histogram without the component\

p[1].plot (x,y, ,{y-el,eu-y}; sym="point", color="green4", size=-1, eb_factor=0);
% optional: add the errorbars
</code>

=== Render the plot ===
<code>
xfig_multiplot (p).render ("plot_components.png");
</code>

===== The Result =====
{{:isis:slxfig:plot_individual_components.png|}}

To remove the annoying vertical lines at the first and last bin of the histogram, you can use the (undocumented) qualifiers ''y_last'' and ''y_first'' of ''hplot''. If you set them to the value of the first or the last bin respectively, the vertical lines vanish:
p[0].hplot (h,m; color="blue4", y_first=m[0], y_last=m[-1]);

Minimalistic SimFit example

2018-04-12T14:38:04Z

Falkner: Created page with "=== Minimalistic SimFit example === The aim of this example is to demonstrate the functionality of the SimFit routines with a minimalistic Example. <code> require("isisscri..."

=== Minimalistic SimFit example ===

The aim of this example is to demonstrate the functionality of the SimFit routines with a minimalistic Example.

<code>
require("isisscripts");

% Define Data
Minimum_Stat_Err = 1e-20;
seed_random(1);

variable N = 30;
variable bin_lo = [1:2:#N+1][[:-2]];
variable bin_hi = [1:2:#N+1][[1:]];

variable data1 = struct{
bin_lo = bin_lo,
bin_hi = bin_hi,
value = (30+5*sin(2*PI*bin_lo) + grand(N)-0.5) * (bin_hi-bin_lo),
err = (9e-1*[1.:1.:#N]) * (bin_hi-bin_lo)
};
variable data2 = struct{
bin_lo = bin_lo,
bin_hi = bin_hi,
value = (35+10*sin(2*PI*bin_lo) + grand(N)-0.5) * (bin_hi-bin_lo),
err = (10e-1*[1.:1.:#N]) * (bin_hi-bin_lo)
};
variable data3 = struct{
bin_lo = bin_lo,
bin_hi = bin_hi,
value = (35+10*sin(2*PI*bin_lo) + grand(N)-0.5) * (bin_hi-bin_lo),
err = (10e-1*[1.:1.:#N]) * (bin_hi-bin_lo)
};
fancy_plot_unit("A"); % Nur zum plotten, falls Daten kein Spektrum sind!
%
% Initialize sim fit
variable sfit = simultaneous_fit();
% Add data
() = sfit.add_data(data1);
() = sfit.add_data([data2,data3]);

% Set fit function
sfit.fit_fun("sinwave2(%)");

% Define global Parameter: Tie parameter sinwave2(%).c of all groups to
% that of group 1
sfit.set_global("sinwave2(%).b");
sfit.set_global("sinwave2(%).c");
sfit.set_par("sinwave2(%).d",1,0,0,1e6);

% load Parameter:
variable file_org = "simpi";
variable file_new = "simpi_new";

variable file;
if( wherefirst(__argv == "org") != NULL ){
file = file_org;
}
else{
file = file_new;
}
vmessage("\n Loading %s.par & %s.model.",file,file);
load_par(file+".par");
sfit.load_model(file+".model");

define save_newpar(){
vmessage("\n SAVING %s.par & %s.model.\n",file_new,file_new);
save_par(file_new+".par");
%sfit.save_model(file_new+".model");
};

if( wherefirst(__argv == "run") != NULL ){
variable runjob = sfit.fit_pars_run_job("simpi_0000.job", ["simpi_0001.job", "simpi_0002.job"],
"simpi_fitpars/*.fits", &save_newpar;
agree,maxiter=3,wait=10);
}
</code>

<code>
#!/bin/bash
#
#PBS -S /bin/bash -V
#PBS -l nodes=2
#PBS -l arch=x86_64
#PBS -l walltime=00:01:00
#PBS -N simpi_fitpars_0000
#PBS -o /home/falkner/Public/simpi_fitpars/simpi_fitpars/simpi_fitpars_0000.out
#PBS -e /home/falkner/Public/simpi_fitpars/simpi_fitpars/simpi_fitpars_0000.err

cd /home/falkner/Public/simpi_fitpars
COMMAND="mpiexec isis-script simpi_fitpars.sl 0"
/usr/bin/nice -n +15 ${COMMAND}
</code>

<code>
#!/bin/bash
#
#PBS -S /bin/bash -V
#PBS -l nodes=2
#PBS -l arch=x86_64
#PBS -l walltime=00:01:00
#PBS -N simpi_fitpars_0001
#PBS -o /home/falkner/Public/simpi_fitpars/simpi_fitpars/simpi_fitpars_0001.out
#PBS -e /home/falkner/Public/simpi_fitpars/simpi_fitpars/simpi_fitpars_0001.err

cd /home/falkner/Public/simpi_fitpars
COMMAND="mpiexec isis-script simpi_fitpars.sl 1"
/usr/bin/nice -n +15 ${COMMAND}
</code>

<code>
#!/bin/bash
#
#PBS -S /bin/bash -V
#PBS -l nodes=2
#PBS -l arch=x86_64
#PBS -l walltime=00:01:00
#PBS -N simpi_fitpars_0002
#PBS -o /home/falkner/Public/simpi_fitpars/simpi_fitpars/simpi_fitpars_0002.out
#PBS -e /home/falkner/Public/simpi_fitpars/simpi_fitpars/simpi_fitpars_0002.err

cd /home/falkner/Public/simpi_fitpars
COMMAND="mpiexec isis-script simpi_fitpars.sl 2"
/usr/bin/nice -n +15 ${COMMAND}
</code>

<code>
require("isisscripts");
_traceback=1;

% TIME LOG
tic;vmessage("%s : Calling: '%s'",time,strjoin(__argv," "));

%%% initialize the simultaneous fit, adding the %%%
%%% data and defining/loading the model here %%%

()=evalfile ("simpi.sl");

% get command line arguments (last argument has to be the group index!)
variable grp = atoi(__argv[-1]);

%%% set the directory for the output fits-files here %%%
%%% make sure that this directory exists! %%%
variable parpath = "/home/falkner/Public/simpi_fitpars/simpi_fitpars";
% check on existence of parpath
if (access(parpath, F_OK) != 0) {
vmessage("error: directory %s does not exist, aborting...", parpath);
exit;
}

% calculate uncertainties
%%% replace 'simfit' by the name of your simultaneous-fit-structure %%%
variable fpars = sfit.mpi_fit_pars( grp ;
dir=parpath,
basefilename = sprintf("simpi_fitpars_%04d", grp)
);

% TIME LOG
variable elaptime = toc;
vmessage("%s : Finished '%s' after %d s < %d %s",
time,strjoin(__argv," "),nint(elaptime),time_array(elaptime)[-1]+1,
["sec","min","h","d","y"][length(time_array(elaptime))-1] );
</code>

Simfit

2018-04-12T14:35:12Z

Falkner: Created page with " Category:Isis / Slang ===== SimFit ===== Example: Minimalistic example"

[[Category:Isis / Slang]]

===== SimFit =====
[[isis:simfit:example|Example: Minimalistic example]]