IBAS ON OFFICE NETWORK
COMPLETE SETUP AND START ON OFFICE NETWORK
STOP AND RESTART AT PERIGEE
EXTRA WORKING SESSION
IBAS FILES
ALERT ANALYSIS
FILE TRANSFER ISDC-->IASF
ATTITUDE AND TIME CONVERSIONS
GRB ANALYSIS
IBAS WEB PAGE UPDATE (LIGHTCURVE)
VARIA
PROCESSI
IDL TRICKS
RICEVERE ALERTS A MILANO
CONTROLLO DEI TRIGGERS E DEGLI ALERTS
ISDC CATALOG
BLACKLIST
USO DI THOR
IBAS WEB PAGES UPDATE (LIGHTCURVES)
CURRENT THRESHOLDS
USEFUL LINKS
PHONE NUMBERS
COMPLETE SETUP AND START ON OFFICE NETWORK
Causa passaggio alle macchine virtuali occorre prima fare login su isdc-nx00 (o isdc-nx01, con gli alias isdc o isdc2) e poi da qui ci si puo' loggare su isdc-ibas00 ( o isdc-ibas01, con gli alias ibas00 o ibas01);
Attenzione ai backup settimanali che bloccano le macchine
isdc-ibas00 TUTTI I LUNEDI alle 9.30
isdc-ibas01 TUTTI I GIOVEDI alle 9.30
quindi usare la macchina adeguata che sia accesa per tutta la rivoluzione
ssh -X -l ibas isdc-nx00.isdc.unige.ch
ssh -X -l ibas isdc-ibas00
source new_ibas_installation/.tcshrc
(oppure usare gli alias isdc su zeus e poi ibas00 o ibas01)
Far partire l'infrastruttura di IBAS (data receipt e alert daemon) con
~/bin/ibas_infrastructure_restart.csh
(dovrebbe fermare tutti i processi attivi e rilanciarli ma nel dubbio si puo' anche fermarli via
top. se solo mocsimul sembra fermo, si puo' anche solo rilanciare mocsimul e vedere se i monitors
processano i nuovi inputs. in caso negativo e' meglio fermare tutto con top (monitors etc) e rilanciare da capo con ibas_infrastructure_restart.csh)
oppure, se in remoto da casa con
~/bin/ibas_infrastructure_restart_remote.csh
Il primo apre 2 finestre che si possono killare con Ctrl-C.
Da una nuova finestra (sulla stessa macchina isdc-ibasXX), andare nella directory con i dati ( /isdc/pvphase/nrt/pck_1/RRRR ) e far partire il mocsimulator (dopo aver killato con top eventuali altre istanze di mocsimul) :
mocsimul psdio_file_name=tm*.fits.gz tx_time_step=0.05 &
Il file di telemetria e' basato sul tempo UTC ed e' del tipo tmYYYYMMDDobt.fits.gz, e.g.
tm2006112012104200.fits.gz. Una volta visto partire mocsimul si puo' chiudere questa finestra.
Tornare sulla prima finestra e far partire i monitors con :
S96imonitors1.new start
S96imonitors2.new start
ed eventualmente il display della curva di luce con :
~/bin/lc2.csh
~/bin/lc.csh
Controllare che vengano creati i workdata in
/unsaved_data/ibas/run/work/ibas/i_monitor_save
si può fare con il comando
~/bin/controlla
Controllare che ci sia abbastanza spazio su /unsaved_data ed eventualmente cancellare i vecchi log-files
Monitorare i triggers con
trigger_tool
p
NOTE
STOP AND RESTART AT PERIGEE
Fermare i monitors
S96imonitors1.new stop
S96imonitors2.new stop
Controllare che non ci sia nulla ancora attivo (vecchi processi)
S96imonitors1.new status
S96imonitors2.new status
Uccidere eventuali jobs rimasti con
kill #job
a) controllare che non ci siano triggers interessanti
trigger_tool
T
100
Q
b) pulire i files della rivoluzione precedente (i workdata e i prodotti salvati in caso di trigger)
pulibas
c) rinominare e gzippare i files della rivoluzione precedente appena conclusasi (RRRR)
~/bin/rinomina.csh RRRR
Monitorare i triggers con
trigger_tool
p
EXTRA WORKING SESSION
Tutta questa parte va fatta sulla macchina ibas00 o ibas01 su cui non sta girando la versione standard
cd work/temp
Rimuovere i workdata, *log, vecchi ma non cancellare le directories
i_monitor_rate_*/ e ontrigger/
Far partire l'infrastruttura IBAS
~/bin/ibas_infrastructure_restart.csh
Da una nuova finestra, andare nella directory con i dati (orbita piu' recente su cui si vuole lanciare
la sessione IBAS, /isdc/pvphase/nrt/pck_1/RRRR) e far partire il mocsimulator con velocita' maggiore
dello standard:
mocsimul psdio_file_name=tm*.fits.gz tx_time_step=0.02 &
Il file di telemetria e' basato sul tempo UTC ed e' del tipo tmYYYYMMDDobt.fits.gz, e.g. tm2006112012104200.fits.gz. Una volta visto partire mocsimul si puo' chiudere questa finestra.
Far partire un monitor_img, un monitor_rate e monitor_save
~/bin/imonitors-linux.csh
(l'analogo di S96imonitors*.new start)
Se e' necessario cambiare i parametri sono in
~/pfiles/ibas_imonitor_img.par
~/pfiles/ibas_imonitor_rate.par
Far partire il display della curva di luce
~/bin/lc2.csh
~/bin/lc.csh
I workdata saranno creati nel direttorio corrente.
Per fermare il tutto dare Ctrl-C sulle finestre dei monitors, ma poi usare top per gli altri processi
IBAS FILES
Le parameter files dei monitors sono in
/unsaved_data/ibas/run/work/ibas/……/*.par
La parameter file di rttmd sta in
$ISDC_ENV/pfiles
La parameter file di ibasalertd sta in
/home/isdc/ibas/pfiles
I log files di alertd sono in
/unsaved_data/ibas/run/log/ibas/ibas_alertd
OFFICE NETWORK: /unsaved_data/ibas/run/work/ibas/
i_monitor_save/ workdata__EVTS_OBTstart_OBTend.fits (~300 sec per file)
i_monitor_img_1/ 5s 15-30 keV monitor=300
i_monitor_img_2/ 10s 15-30 keV monitor=301
i_monitor_img_3/ 40s 15-200 keV monitor=302
i_monitor_img_4/ 10s 15-40 keV monitor=303
i_monitor_img_5/ 100s 15-200 keV monitor=304
i_monitor_img_6/ 10s 15-200 keV monitor=305 (dalla rev.1000)
i_monitor_img_7/ 20s 15-200 keV monitor=306 (dalla rev.1000)
i_monitor_rate_1/ 2,8,32 ms 15-200 keV monitor=200
i_monitor_rate_2/ 0.08, 0.32, 1.28, 5.12 s 15-200 keV monitor=201
i_monitor_rate_3/ 0.08, 0.32, 1.28, 5.12 s 15-30 keV monitor=202
i_monitor_rate_7/ 20,80,320 ms 40-200 keV monitor=206
i_monitor_rate_8/ 10,40,160,640 ms 40-200 keV monitor=207 (dalla rev.1000)
NOTA: gli operatori stoppano e riavviano copy_nrt che da` subito accesso sull'office
network agli ultimi files di attitudine e telemetria
OPERATIONS NETWORK
i_monitor_img_1/ 10s 15-200 keV monitor=300
i_monitor_img_2/ 20s 15-200 keV monitor=301
i_monitor_img_3/ 40s 15-200 keV monitor=302
i_monitor_img_4/ 10s 15-40 keV monitor=303
i_monitor_img_5/ 100s 15-200 keV monitor=304 (with strong threshold = 20)
i_monitor_rate_1/ 2,8,32 ms 15-200 keV monitor=200
i_monitor_rate_2/ 0.08, 0.32, 1.28, 5.12 s 15-200 keV monitor=201
i_monitor_rate_3/ 0.01, 0.04, 0.16, 0.64 s 40 -200 keV monitor=202
CURRENT THRESHOLDS:
possible/strong 5.5 / 8. Office
6.5 / 8. Operat. (tranne il mon. 304 che ha la strong thr a 20, vedi sopra)
FILE TRANSFER ISDC-->IASF
Dall'ISDC fare:
sftp sandro@zeus.lambrate.inaf.it
mput *.gz
ATTITUDE AND TIME CONVERSIONS
Ricavare l'attitudine del satellite al momento del GRB con
GetAttitude
ssh -X -l ibas crab.unige.ch
ssh -X isdcsf4
source .ibaslogin
NOTA: GetAttitude richiede il tempo UTC del GRB e la "Attitude type".
Quando non ci sono ancora i files di correlazione dare la
stessa ora del giorno prima e aggiungere 86400.
GetAttitude
-----> In output vengono date "S/C_ATT" e "INSTR_ATT". Usare "INSTR_ATT" in IBAS_deconvolve.
converttime
Per usare converttime bisogna essere collegati a Ginevra (e.g. isdcsf4)
ssh -X -l ibas login.isdc.unige.ch
ssh -X isdcsf4
source .ibaslogin
converttime
Per passare da OBT a UTC (e viceversa) bisogna passare per IJD (e.g. prima da UTC a IJD poi da IJD a OBT).
NOTA: OBT nei triggers di IBAS e' in secondi e va trasformato per usarlo in converttime come segue:
idl
obt_trigger= 129363347
obt_converttime=double(obt_trigger)*1024.D0*1024.D0
print, obt_converttime, format='(d30.0)'
exit
[Per il contrario e.g. print, 221995795643133.d0/1024./1024,format='(f20.10)']
Usare obt_converttime come input a converttime e una volta ottenuto il tempo del GRB in UTC tornare a GetAttitude.
GRB ANALYSIS
L'analisi di un GRB consiste di 3 parti: [1] la creazione del file di eventi in un dato
intervallo di tempo (IBAS_select, che crea anche una curva di luce), [2] la creazione
dell'immagine nello stesso intervallo (IBAS_deconvolve) e [3] l'estrazione dello
spettro (IBAS_spectrum). Sono disponibili alcuni tools in IDL.
source /red/Software/isdc/dev6.3/.ibaslogin
cp /red/Software/isdc/dev6.3/pfiles/IBAS_* ~/pfiles/
IBAS_select [ multienergy=yes ] [ usepif=yes ]
Prima di IBAS_deconvolve ricordarsi di lanciare IBAS_select sulla durata del GRB!!!!
IBAS_deconvolve
IBAS_spectrum nchan=64 binning=2
[binfilename=/zeus/sandro/ibas/isgr_rmf_64_osa3_t3.txt]
(matrice a 64 ch binnata log a partire da 13 keV)
Lo spettro risultante e' IBAS_RAW_GRB_....
xspec11
data IBAS_RAW_GRB_....
resp /zeus/sandro/ibas/isgr_rmf_64_osa3.fits.gz
arf /zeus/sandro/ibas/isgr_arf_rsp_0005.fits
ignore **-20.
cpd /xs
setplo en
plo lda
mo po
renorm
fit 100
plo lda del
flux 20. 200.
Ricordarsi che i flussi derivati vanno moltiplicati per 2!
Analisi IDL
- per vedere la curva di luce e selezionare eventi e detector
siena,det
salvare gli eventi corrispondenti al GRB
- per trovare il ghost corretto
findburstposition,det
- per leggere i files di eventi
evt_read,'workdata_EVTS_0028943713_0028944013.fits',y,z,t,pha,rise
nrg_corr2,y,z,pha,rise,e,pi
PROCESSI
per vedere i processi usare top
u user
n XX XX=n proc. da mostrare
k -9 PID per killare
q per uscire
Ctrl-C killa
Ctrl-Z sospende per riprendere si usano:
bg (mette in BKG) oppure
fg (riporta in foreground)
IDL TRICKS
conversioni
print,sixty(23.23)
23.0000 13.0000 47.9984
print,ten(23,13,47.9984)
23.230000
distanze angolari
gcirc
Calling sequence: GCIRC,U,RA1,DC1,RA2,DC2,DIS
U = 0 ==> Everything in radians
U = 1 ==> RAx decimal hours, DCx decimal degrees, DIS arc sec
crabra = ten ( 5.d0 , 34. , 32 )
crabdec = ten ( 22.d0 , 0. , 52.1)
gcirc,1,crabra,crabdec,grb_ra/15.d0,grb_dec,d & print,d/60.
off-axis angle
con copy/paste dall'Alert message copiarsi
grb_ra grb_dec point_ra point_dec
gcirc,1,point_ra/15.,point_dec,grb_ra/15.d0,grb_dec,d & print,d/3600.,' deg'
oppure dal trigger message copy/paste grb_y e grb_z
print,'Y=', atan(grb_y)*180./!pi,' Z=',atan(grb_z)*180./!pi,' deg'
distanze da strong Sources
;##### Cyg_X-1
ra = ten( 19,58,21.7) & dec = ten( +35,12,05.8)
gcirc,1,ra,dec,grb_ra/15.d0,grb_dec,d & print,d/60.,' arcmin'
;##### Crab
ra = ten( 05,34,32.0) & dec = ten( +22,00,52.1)
gcirc,1,ra,dec,grb_ra/15.d0,grb_dec,d & print,d/60.,' arcmin'
;##### Vela_X-1
ra = ten( 09,02,06.9) & dec = ten( -40,33,16.9)
gcirc,1,ra,dec,grb_ra/15.d0,grb_dec,d & print,d/60.,' arcmin'
;##### Sco_X-1
ra = ten( 16,19,55.1) & dec = ten( -15,38,25.0)
gcirc,1,ra,dec,grb_ra/15.d0,grb_dec,d & print,d/60.,' arcmin'
;##### Her_X-1
ra = ten( 16,57,49.8) & dec = ten( +35,20,32.6)
gcirc,1,ra,dec,grb_ra/15.d0,grb_dec,d & print,d/60.,' arcmin'
;#### SGR 1806-20
ra = ten(18, 08, 39.3) & dec = ten(-20,24,39.5)
gcirc,1,ra,dec,grb_ra/15.d0,grb_dec,d & print,d/60.,' arcmin'
to read images
fxread,nomefile,array ; se nella prima estensione
fits_read,nomefile,array,exten_no=1 ; altrimenti
RICEVERE ALERTS A MILANO
***questa parte è obsoleta ****
per controllare connessione da una Linux (zeus) alla macchina IBAS su operation network :
ibas_ping 129.194.168.222 1966 1944 (IBAS clients list in /isdc/WWW/Soft/ibas/config_files a ginevra)
fermarlo con Ctrl-C (attenzione: non funziona se alert reciever sta girando)
per ricevere alert messages:
ibas_alert_receiver 129.194.168.222 1966 1944
per ricevere alert messages e mandarli su un log file:
ibas_alert_receiver 129.194.168.222 1966 1944 dump_alert.csh &
ti manda gli alerts sia su schermo sia in un logfile che si chiama alerts.log e viene creato nel direttorio dove lanci il programma (viene creato quando viene ricevuto il primo alert)
CONTROLLO DEI TRIGGERS E DEGLI ALERTS
ssh -X -l ibas crab.unige.ch
source .ibaslogin
trigger_tool [tail trg_tail=40]
alert_tool
ISDC CATALOG
A Milano sono in : /zeus/sandro/ibas/
con numerazione progressiva dal piu` vecchio al piu` recente.
La region file corripondente e' sempre nella stessa dir.
e.g. /zeus/sandro/ibas/gnrl_refr_cat_0010.fits
/zeus/sandro/ibas/region_31.reg
BLACKLIST
OFFICE NETWORK
${ISDC_SITE}/run/work/ibas/ibas_alertd/etc/black_list.txt
e dovrebbe essercene anche una copia in ~diego/blacklist
OPERATIONAL NETWORK
dovrebbe essercene una copia in /isdc/WWW/Soft/ibas/config_files
USO DI THOR
**** questa parte obsoleta ***
internet ---KPPP
dopo aver dato da root il comando
>tcsh
>ifconfig eth0 down
#basi02% (#=Shift3)
per copiare i files da bacchus fare
scp user@macchina...:file .
IBAS WEB PAGE UPDATE (LIGHTCURVES)
USEFUL LINKS
OPERATIONS
ALERTS SU OPS.
ISDC Planning Page
INTEGRAL short term plan
SIMBAD
GCN Circulars