Date post: | 03-May-2015 |
Category: |
Documents |
Upload: | marietta-mariotti |
View: | 214 times |
Download: | 0 times |
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
1
Stato e Prospettive del Calcolo di MAGIC
Major
Atmospheric
Gamma
Imaging
Cherenkovtelescope
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
2
Outline
• Cosa e’ MAGIC • Stato attuale di MAGIC • Quantita’ di dati attesi in un
anno • Analisi dati • Produzione MC necessaria • Produzione MC al CNAF • Conclusioni
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
3
MAGIC
• O.R.M. 2200m s.l.m.• Struttura leggera in
fibra di carbonio• Specchio parabolico
diametro 17m f1 • Camera con 577
PMT• Trasmissione segnali
via fibra ottica• FADC @ 300MHz
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
4
Stato attuale
Inaugurato 10/10/2003
Regolare presa dati da Febbraio 2004
Rate ≈ 200Hz -> 1KHz => Eth ≈ 80GeV ->
30GeVTotalizzate 112 ore di
osservazione Visto segnale Crab e
Markarian 421
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
5
Flusso dati reali
DAQ @ LaPalmaFileserver Tapes x 2
Wuerzburg Data Center
Regional Data Center 1
Regional Data Center 2
Regional Data Center 3
Tapes ?MERPP
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
6
Eventi attesi
Rate segnale ~ 1Hz Rate fondo = 500Hz DutyCycle = 900h/1anno = 10%
Eventi = Rate * 1 anno * DutyCycle
Eventi segnale = 3.2 MeventiEventi fondo = 1.6 Geventi
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
7
Necessita’ di Storage
Neventi = 1.6 Gevs (@ 500Hz)
RawDATA = 1.6Gevs * 18.3KB/evs = 27TB
=> Circa 140 LTO2 tapes x 2
Rootified = 27TB * 0.53 = 14 TB
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
8
Data preprocessing
MERPP(rootification)Raw Data
Sottrazione dei piedistalli+
Calibrazione
Calcolo dei parametri di Hillas
?
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
9
Velocita’ del data preprocessing
MERPPing 410evs/sCalibrazione 225evs/s Calcolo Parametri Hillas 980evs/s
Ttot = 7.9 ms/evento
Velocita’ = 125 eventi/s
Su XEON @ 3.0 GHz
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
10
CPU per il data processing
Rate acquisizione 500->1000 Hz
Rate preprocessing 125 Hz
4 -> 8 CPU tipo XEON @ 3.0 GHz
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
11
Analisi
L’analisi non termina con calcolo Hillas:separazione gamma/adroni, stima energia, ecc.
Metodo false sorgenti. T x 30 Metodi analisi alternativi T x 100 (?)
Velocita’ dell’analisi < 125evs/s
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
12
Il MC di MAGIC 1 (responsabilita’
italiana)Simulazione sciami
atmosferici CORSIKA (f77)
Riflessione + Assorbimento atmosferico
Reflector (c)
Fotoconversione +Simulazione elettronica
Camera (c++)
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
13
Produzione Standard
• Gammas 10 GeV < E < 30 TeV R < 400 m Spectral Index 2.6
<E> = 26GeV
• Hadrons 30 GeV < E < 30 TeV R < 300 m Spectral Index 2.75 p 2.62 He
<E> = 70GeV
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
14
Dati MC necessari
No test beam => Dipendenti da MCAlmeno tanti sciami di segnale simulati
quanti quelli attesi (10 x ?). Almeno 1/100 (1/10 ?) degli sciami di
fondo.
3.2 - 32M eventi di segnale 160 - 16 M eventi di fondo
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
15
Dati MC necessari
Efficienza di trigger per gamma 7% Efficienza di trigger per adroni 0.15%
3.2 Meventi => 45 Meventi prima del trigger
16Meventi => 11 Geventi prima del trigger
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
16
Velocita’ MC
Corsika 900evs/hReflector 57000evs/h Camera ~ 10^6evs/h
Ttot = 4 s/evento
Velocita’ = 900 eventi/h
Su XEON @ 3.0 GHz
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
17
CPU per MC
11Geventi @ 900evs/h => 1.2*10^7h
1.2*10^7h/1anno = 1200cpu
Possiamo usare lo stesso evento di fondo piu’ volte con diverso parametro di impatto
≈50 CPU
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
18
Storage per MC
• Corsika output 28KB/evs 10.8KB/evs• Reflector output 7.6KB/evs 1.3KB/evs
Conserviamo l’output del Reflector
Gamma Adroni
45Meventi gamma => 320 GB
11Geventi adronici => 13 TB
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
19
CONDOR
Circa 300 CPUs disponibili, 30 Alpha + 270 PC
0
50
100
150
200
250
Time a.u. (3 days)
CP
Us
In media 105 CPUs simulano sciami per MAGIC
Utilizzo intensivo di CONDOR negli ultimi 3 anni, usando script csh => 90Meventi prodotti
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
20
Produzione al CNAF
• Test iniziati il 5 Marzo con 1 nodo• Passaggio semplice da CONDOR a CNAF• Script PBS + generatore di steering card• Adesso 5 nodi dedicati + 5 nel pool
comune• 1TB di spazio disco
(Si ringrazia lo staff del CNAF per il prezioso aiuto)
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
21
La produzione al CNAF oggi
• Circa 2 mesi di esperienza– Nessun problema di rilievo
• Stiamo usando ~ 60 CPUs• Generati finora circa 150Kevs gamma
HE equivalenti a circa 30 Mevs gamma STD – Circa raddoppiata la potenza disponibile su
CONDOR– Il CNAF produce il 100% del MC di MAGIC 1
27/05/2004 C.Bigongiari & M.Pieraccini INFN Workshop Castiadas (CA)
22
Conclusioni
CPU x produzione MC ~50
Storage x MC 13TB
MC di MAGIC 1 interamente prodotto al CNAF
CPU x Analisi Dati ~10
Storage x Dati 27TB (@ production centre)