Personal tools
You are here: Home Virgo Data Analysis GRB GRB070729
Document Actions

GRB070729

by Elena Cuoco last modified 2009-03-25 11:12
**GRB070729**


GRB trigger time:  GPS 869703967 (Jul 29 2007 00:25:53 UTC)

--------------

Posizione:

RA(J2000)  =  03h 45m 17s
Dec(J2000) = -39d 18' 41"

--------------

Detector  F_+  F_x  F_+^2 + F_x^2

H1  -0.83363  -0.17174  0.72444
H2  -0.83363  -0.17174  0.72444
L1  0.91529  0.31438  0.93658
V1  -0.046994  -0.59314  0.35403

---------------

Regione on-source: [869703847,869704027] (i.e. [GPS - 120 s,GPS + 60 s]).

Regione off source: [869698567,869709367] (i.e. [GPS - 1.5 ore, GPS + 1.5 ore]).

--------------

Segnali simulati: i dettagli si trovano in
https://ldas-jobs.ligo.caltech.edu/~eharstad/ExtTrigMDCSetS5GRB.html, mentre i
frame veri e propri sono salvati sul cluster LIGO del CIT, come descritto alla
fine di questa pagina web.


--------------


Situazione MDC (Alessandra):


- frame: si trovano in /archive/frames/S5/ExtTrig-MDC/GRB070729, e direi che
gia' che ci siamo convinene scaricare tutti i segnali sine gaussian, sia quelli
polarizzati circolarmente che linearmente, a tutte le frequenze (che vanno da
100 Hz a 1000 Hz). Poi gli americani ci chiedevano di iniziare l'analisi dalla
sine gaussian a 1000 Hz, quindi magari partiamo da quella. Per Alberto: le
cartelle che vanno scaricate sono quelle chiamate SG*Q8d9_OFF, che contengono
le iniezioni per la regione off-source che ha noi ci interessa.

- i file log che spiegano cosa c'e' dentro i frame file si trovano in
/data/node5/eharstad/ExtTrigMDC/GRB070729/. Per Alberto: anche qui occorre
scaricare tutte le cartelle del tipo SG*Q8d9_OFF.


- ogni frame porta il nome del rivelatore, es. H1, H2, L1, V1, V2. Virgo pur
essendo uno solo ha sia V1 che V2 perche' i V1 sono campionati a 20 KHz, V2 a
16 kHz. Penso che sia meglio utilizzare il V2 per uniformita' con i frame LIGO.

- in ogni frame ci sono 12 canali. Ogni canale contiene una iniezione che viene
iniettata a tempi random all'interno del frame

- le forme d'onda simulate sono normalizzate in modo da avere un'ampiezza hrss
alla sorgente di 1e-21 (cioe' prima degli F+ ed Fx del rivelatore), ed i frame
file ***gia' contengono il segnale moltiplicato per F+ ed Fx****. QUESTA COSA
E' IMPORTANTE: A DIFFERENZA DI QUANTO FACEVAMO CON LE INJECTION DI VIRGO SIAMO
UN PASSO AVANTI, CIOE' NON OCCORRE MOLTIPLICARE H+ PER F+, HX PER FX E POI
SOMMARE, MA L'OPERAZIONE LA HANNO GIA' FATTA GLI AMERICANI! (Penso che Elena
soprattutto sara' molto contenta di cio'...). Quindi l'unica operazione che
rimane da fare e' prendere le injection, moltiplicarle per i fattore di scala
di normalizzazione e sommarle ai dati.

---------------------------------


Posizione dei dati hrec di Virgo (Gianluca)


la versione che e' attualmente usata nell'analsi e' la v2, da utilizzare
con i periodi
vetati che puoi trovare utilizzando il Virgo Data base. Esiste anche una
versione
3, nella quale un errore che ha coinvolto pochi frames nella
ricostruzione e' corretto.
La versione 2 la puoi trovare a Bologna:
/storage/gpfs_virgo4/virgo/data/VSR1/hrec-v2

dove trovi anche i dati di LIGO:
/storage/gpfs_virgo4/virgo/data/LIGO/S5
La versione 3 la trovi a Cascina:

/virgoData/ffl/VSR1/HrecV3.ffl

Trovi tutto poi a Lione: ti puoi riferire alla pagina web:
https://workarea.ego-gw.it/ego2/virgo/data-analysis/data/data-in-lyon/data-in-lyon/

(La pagina web corrispondente per Bologna non e' aggiornata)
Per esempio, la versione V2 sta in
/hpss/in2p3.fr/group/virgo/Run/VSR1/proc/V2
I dati di H1 di LIGO in:
/hpss/in2p3.fr/group/virgo/Run/S5/H-H1_RDS_C03_L2/pack/

---------------------------------


Accesso ai dati su Grid (Alberto)


 SIMULAZIONI


Il dataset e' registrato a Roma (virgo-se01.roma1.infn.it) e replicato a Pisa (gridsrm.pi.infn.it).


Il base folder sul file catalog e' /grid/virgo/GRB/GRB070729.

i frame sono contenuti nelle seguenti cartelle:

$ lfc-ls /grid/virgo/GRB/GRB070729

SGC1000Q8d9_OFF

SGC100Q8d9_OFF

SGC150Q8d9_OFF

SGC1850Q8d9_OFF

SGC250Q8d9_OFF

SGC554Q8d9_OFF

SGL1000Q8d9_OFF

SGL100Q8d9_OFF

SGL150Q8d9_OFF

SGL1850Q8d9_OFF

SGL250Q8d9_OFF

SGL554Q8d9_OFF

WNB100ms175to325Hz_OFF

WNB1200ms300to700Hz_OFF


I frame hanno lunghezza di 256 s. Ogni cartella contiene i frame di ciascuno dei 6 rivelatori (G1, H1,
H2, L1, V1, V2). Il periodo coperto e' di 3h circa (43 file per rivelatore). Di seguito un esempio di
nome di un frame file:

V-V2_SGC1000Q8d9_OFF_GRB070729-869704607-256.gwf


Oltre ai singoli frame ho pensato di registrare tutti i frame corrispondenti a ciascun rivelatore,
compressi in un file .tar.gz. Questo dovrebbe semplificare lo scaricamento dei dati sui worker nodes,
assumendo che ogni job analizzera' una singola frequenza ed un singolo rivelatore: invece di scaricare
i 43 file corrispondenti bastera' scaricare il singolo tar.gz ed estrarre i file. Il file compresso e'
nella stessa cartella dei frame, ad es.:

/grid/virgo/GRB/GRB070729/SGC1000Q8d9_OFF/H-H1-SGC1000Q8d9_OFF_GRB070729.tgz


In /grid/virgo/GRB/GRB070729/ trovate poi il file md5sum.txt che contiene la checksum di tutti i file:

$ head md5sum.txt

4eb8f4c0a7843de4e1d4ef92877dcfaa  SGC1000Q8d9_OFF/G-G1_SGC1000Q8d9_OFF_GRB070729-869698463-256.gwf

4af6f81a51a36a28440f97698c251f7b  SGC1000Q8d9_OFF/G-G1_SGC1000Q8d9_OFF_GRB070729-869698719-256.gwf

7ae47fc360ab24f082203075e4648c65  SGC1000Q8d9_OFF/G-G1_SGC1000Q8d9_OFF_GRB070729-869698975-256.gwf

...


puo' essere utile per verificare la consistenza dei file appena scaricati.


In ciascuna directory ho inoltre messo il folder logs che contiene i log files. Su suggerimento di
Alessandra ho estratto dai log i ctrgps times (unendo i s con i ns) e li ho messi nei file:

$ lfc-ls /grid/virgo/GRB/GRB070729/SGC1000Q8d9_OFF/logs |grep ctrgps
SGC1000Q8d9_OFF_GRB070729-G1-ctrgps.txt
SGC1000Q8d9_OFF_GRB070729-H1-ctrgps.txt
SGC1000Q8d9_OFF_GRB070729-H2-ctrgps.txt


In ogni file trovate quindi tutti i tempi dei segnali per ciascun rivelatore:

$ cat GRB070729/SGC1000Q8d9_OFF/logs/SGC1000Q8d9_OFF_GRB070729-H2-ctrgps.txt
869698643.386934757
869698537.881748318
869698493.559396505
869698479.567933440
869698598.590461850
869698612.309402704
869698522.711085438
869698508.977373957
869698554.780500531
869698584.634704470
...


Sono 12 canali * 43 frames = 516 tempi


DATI

Come per i file delle simulazioni ho registrato sia i singoli frame file sia il tar gz di tutti i
frame di ciascun rivelatore (a seconda di come procederemo con l'analisi eliminero' il formato dei
dati che non ci servira').

  • VIRGO

La versione v3 dei dati del VSR1 si trova a Cascina (la lista completa e' in /virgoData/ffl/VSR1/HrecV3.ffl).

Il logical file name dei frame singoli e':

/grid/virgo/GRB/GRB070729/data/VIRGO/VSR1/HrecV3/V-HrecV3-869690000-10000.gwf
/grid/virgo/GRB/GRB070729/data/VIRGO/VSR1/HrecV3/V-HrecV3-869700000-10000.gwf

mentre il file .tgz che li contiene entrambi e':

/grid/virgo/GRB/GRB070729/data/VIRGO/VSR1/HrecV3/V-HrecV3.tgz

Ho poi registrato un file di testo con l'md5sum dei due file:

/grid/virgo/GRB/GRB070729/data/VIRGO/VSR1/HrecV3/V-HrecV3-md5sum.out

  • LIGO

I frame si trovano a Bologna (/storage/gpfs_virgo4/virgo/data/LIGO/S5/strain-L2/L*O/...). Ho trovato i
dati relativi a H1, H2 e L1. Ciascun frame ha lunghezza 128 s.

Sul catalogo ho mantenuto due directory distinte per LHO e LLO.

Il LFN dei frame singoli e':

/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H1_RDS_C03_L2-869698401-128.gwf
/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H1_RDS_C03_L2-869698529-128.gwf
...
/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H2_RDS_C03_L2-869698343-128.gwf
...
/grid/virgo/GRB/GRB070729/data/LIGO/LLO/L-L1_RDS_C03_L2-869698404-128.gwf
...


I file .tgz si chiamano:

/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H1_RDS_C03_L2.tgz
/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H2_RDS_C03_L2.tgz

/grid/virgo/GRB/GRB070729/data/LIGO/LLO/L-L1_RDS_C03_L2.tgz

Anche in questo caso ho messo a disposizione i file testo con le checksum dei singoli frame file:

/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H1_RDS_C03_L2-md5sum.out
/grid/virgo/GRB/GRB070729/data/LIGO/LHO/H-H2_RDS_C03_L2-md5sum.out
/grid/virgo/GRB/GRB070729/data/LIGO/LHO/L-L1_RDS_C03_L2-md5sum.out