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)

 

  1. Creare la curva di luce del GRB
      1. In una working directory copiare il workdata file contenente il GRB;
      2. Lanciare IBAS_select selezionando il picco del GRB (gli step 2 e 3 servono per costruire la PIF);
      3. Lanciare IBAS_deconvolve (ATTENZIONE: non ci deve essere nessuna PIF.fits nella dir di lavoro);
      4. Collegarsi a Ginevra per lanciare converttime:
        1. ssh -X -l ibas login.isdc.unige.ch
          ssh -X isdcsf4
          source .ibaslogin 
        2. converttime  da UTC (poco prima di burst) a OBT e poi mettere i valori come input a make_lc_plot, e.g.
          Log_1  : Input Time(UTC): 2007-03-09T10:00:00.000 Output Time(OBT): 145348921502228
          per controllare i tempi dei workdata, in IDL: print, 145348921502228.d0/1024./1024,format='(f20.10)'
      5. in IDL (esempio per GRB070309)  :
                      burst='070309'
                      infile='/IL_MIO_PATH/workdata__EVTS_0138615346_0138615646.fits'
                      utstring='10:00:00'
                      obt0 =145348921502228
                      obt0 = double(obt0)/2.d0^20.
                      .run  /zeus/sandro/ibas/analysis/make_lc_plot
                      make_lc_plot,burst,infile,'PIF.fits',obt0,utstring,1    ;(mettere "0" al posto di "1" se non si vuole usare PIF.fits:   make_lc_plot,burst,infile,'PIF.fits',obt0,utstring,0)                
      6. qdp LCPLOT070309.qdp (sistemarlo e salvarlo come .ps:  hardcopy lc.ps/PS)
      7. xv  lc.ps e salvare jpg con qualita' massima ---> YYMMDDlc.jpg
  2. Copiare il file jpg in /ariel/www/ibas/products/
  3. In /ariel/www/ibas
    creare ed editare il file YYMMDD.html e cambiare il link a YYMMDDlc.jpg
  4. cp  /ariel/www/ibas/IBAS_Results.html /ariel/www/ibas/IBAS_Results.html.old ed editare  /ariel/www/ibas/IBAS_Results.html
    <tr> <td> YYMMDD </td>
    <td> 02:54:21 </td> RA
    <td> + 69:07:34</td> DEC
    <td> 2</td> UNC
    <td> 0.3 </td> Peak FLUX
    <td> 00:05:56</td> UT
    <td> 40</td> T90
    <td> X/O?</td> AG
    <td> -</td>
    <td> <a href="http://gcn.gsfc.nasa.gov/gcn/other/YYMMDD.gcn3">here</a></td>
    <td>-</td>
    <td><a href="http://ibas.iasf-milano.inaf.it/YYMMDD.html">Here</a></td>

USEFUL LINKS

OPERATIONS
ALERTS SU OPS.
ISDC Planning Page
INTEGRAL short term plan
SIMBAD
GCN Circulars