Simulazioni Montecarlo per un telescopio di neutrini sottomarino

Per il telescopio di neutrini ANTARES, sul quale è basata questa tesi, il principale problema nella ricostruzione delle tracce dei muoni provenienti dal basso, i quali possono essere generati solo da neutrini provenienti dall’emisfero sud che interagiscono nei pressi del rivelatore, è rappresentato dai muoni atmosferici (in particolare dagli sciami muonici). La luce Cherenkov da essi prodotta può essere utilizzata per ricostruire erroneamente tracce dirette verso l’alto (upgoing). Lo scopo di questa tesi è la valutazione, tramite un metodo Montecarlo, della probabilità che i muoni atmosferici vengano ricostruiti come muoni diretti verso l’alto. L’esposizione è organizzata in quattro capitoli: nei primi due vengono descritti la fisica dell’astronomia a neutrini e i principi di funzionamento di un telescopio di neutrini, con particolare attenzione al telescopio ANTARES.Nel terzo capitolo è riportata la descrizione del software utilizzato per simulare lo sciame di muoni atmosferici e la conseguente risposta del rivelatore. Infine, nel quarto capitolo, sono presentati i grafici eseguiti con ROOT che riassumono ed evidenziano i risultati della simulazione.

Neutrini dal piano galattico con ANTARES

In questa tesi si cerca un flusso diffuso di neutrini nel piano Galattico. ANTARES, un telescopio di neutrini che si trova nel Mar Mediterraneo, ha un’area di visibilit`a tale da poter osservare questa parte di cielo, quindi isolare l’emissione di neutrini per trovarne le sorgenti. Sono confrontate diverse regioni del cielo, selezionate in modo tale che abbiano le stesse dimensioni. In particolare si comparano la regione corrispondente al Centro Galattico, detta on-zone, e altre regioni, dove secondo i modelli analizzati ci si aspetta un minor segnale, dette off-zone.

Confronto tra la direzione di arrivo dei neutrini di ANTARES e IceCube

[Abstract]. Questo lavoro è finalizzato al confronto delle rivelazioni di canditati eventi di neutrino cosmico da parte di due esperimenti ANTARES ed ICECUBE. I due apparati, entrambi basati sulla rivelazione di radiazione Cherenkov generata da prodotti di interazione CC dei neutrini, si differenziano per la geometria del detector e per il mezzo in cui viene prodotta la radiazione, e hanno quindi sensibilità e risoluzione diverse. L'analisi dei dati è basata sulla rappresentazione delle direzioni di arrivo degli eventi su mappa celeste nei due sistemi di coordinate galattico ed equatoriale. La sovrapposizione della distribuzioni degli eventi ha permesso di evidenziare le zone di visibilità sulla mappa celeste caratteristiche di ciascun esperimento. Una successiva analisi statistica ha analizzato la frequenza delle occorenze di eventi del dataset di ANTARES in intorni di raggio crescente di eventi di ICECUBE, per valutare eventuali eventi identificabili nei due dataset. Questa analisi ha evidenziato la necessità di una selezione preventiva degli eventi del dataset di ANTARES, in base a criteri non solo topologici, in modo da permettere una eventuale identificazione di eventi coincidenti.

Correlazioni tra UHECRs di Auger e neutrini di ANTARES

Questo lavoro ha lo scopo di individuare una strada nella ricerca di un confronto tra canditati eventi di neutrino cosmico da parte dell'esperimento ANTARES ,nel Mar Mediterraneo a sud est della costa di Tolone, in Francia e gli Ultra High Energy Cosmic Rays dell'osservatorio Pierre Auger in Argentina. Si utilizzano dati raccolti dal telescopio di neutrini ANTARES tra il 2007 e il 2013, con energie dell'ordine del TeV, e quelli di Auger ottenuti tra il 2004 e il 2014, con energie dei raggi cosmici dell'ordine dell'EeV. Entrambi gli apparati sono basati sulla rivelazione di radiazione Cherenkov, ma l'osservatorio di Auger è dotato di un secondo rivelatore di luminescenza che lavora in sincrono per ridurre gli errori sistematici associati ai singoli rivelatori. Le distribuzioni degli eventi in base alle direzioni di arrivo e dell' energia osservata, sono i parametri necessari per la ricerca di una correlazione. Dalla sovrapposizione delle distribuzioni possiamo evidenziare le zone di visibilità caratteristiche di ciascun esperimento, e condurre un eventuale successiva analisi finalizzata ad una identificazione univoca di eventi coincidenti.

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completely absent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involved parts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method is introduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that the theoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approach has reasonable properties from a compositional point of view. In particular, it is “natural” in the sense that it recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in the same paper a substitution method for missing values on compositional data sets is introduced

Valoración de activos financieros utilizando Montecarlo: cuando el mundo no es tan normal.

Valorar activos financieros cuando el mundo no es tan normal como asumen muchos modelos financieros exige un método flexible para operar con diversas distribuciones, el cual, adicionalmente, pueda incorporar discontinuidades como las que se dan en procesos estocásticos de salto. El método Monte Carlo cumple con esos requisitos, además de generar una buena aproximación y ser eficiente, lo cual lo convierte en el más adecuado para aquellos casos en los cuales no se cumple el supuesto de normalidad. Este artículo explica como se aplica este método para la valoración de activos financieros, particularmente opciones financieras, cuando el activo subyacente sigue un proceso de volatilidad estocástica o de salto-difusión.

Reliability assessment of energy limited systems using sequential montecarlo simulation

This paper explains why the reliability assessment of energy limited systems requires more detailed models for primary generating resources availability, internal and external generating dispatch and customer demand than the ones commonly used for large power systems and presents a methodology based on the full sequential Montecarlo simulation technique with AC power flow for their long term reliability assessment which can properly include these detailed models. By means of a real example, it is shown how the simplified modeling traditionally used for large power systems leads to pessimistic predictions if it is applied to an energy limited system and also that it cannot predict all the load point adequacy problems. © 2006 IEEE.

Caracterización computacional de un detector de germanio hiperpuro de rango extendido (HPGe- XtRa) con simulación Montecarlo y optimización con un algoritmo evolutivo

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Studies of the atmospheric muon flux with the ANTARES detector

The thesis main topic is the determination of the vertical component of the atmospheric muon flux as a function of the sea depth at the ANTARES site. ANTARES is a Cherenkov neutrino telescope placed at 2500m depth in the Mediterranean Sea at 40 km from the southern cost of France. In order to retrieve back the physical flux from the experimental data a deconvolution algorithm has been perform which takes into consideration the trigger inefficiensies and the reconstruction errors on the zenith angle. The obtained results are in good agreement with other ANTARES indipendent analysis.

Tecniche MonteCarlo per l'analisi di sistemi a rete

The assessment of the RAMS (Reliability, Availability, Maintainability and Safety) performances of system generally includes the evaluations of the “Importance” of its components and/or of the basic parameters of the model through the use of the Importance Measures. The analytical equations proposed in this study allow the estimation of the first order Differential Importance Measure on the basis of the Birnbaum measures of components, under the hypothesis of uniform percentage changes of parameters. The aging phenomena are introduced into the model by assuming exponential-linear or Weibull distributions for the failure probabilities. An algorithm based on a combination of MonteCarlo simulation and Cellular Automata is applied in order to evaluate the performance of a networked system, made up of source nodes, user nodes and directed edges subjected to failure and repair. Importance Sampling techniques are used for the estimation of the first and total order Differential Importance Measures through only one simulation of the system “operational life”. All the output variables are computed contemporaneously on the basis of the same sequence of the involved components, event types (failure or repair) and transition times. The failure/repair probabilities are forced to be the same for all components; the transition times are sampled from the unbiased probability distributions or it can be also forced, for instance, by assuring the occurrence of at least a failure within the system operational life. The algorithm allows considering different types of maintenance actions: corrective maintenance that can be performed either immediately upon the component failure or upon finding that the component has failed for hidden failures that are not detected until an inspection; and preventive maintenance, that can be performed upon a fixed interval. It is possible to use a restoration factor to determine the age of the component after a repair or any other maintenance action.

Realisierung der IceCube MonteCarlo-Produktion im WLCG und Untersuchung von Auswirkungen meteorologischer Parameter auf die Erzeugung atmosphärischer Myonen

Das am Südpol gelegene Neutrinoteleskop IceCube detektiert hochenergetische Neutrinos über die schwache Wechselwirkung geladener und neutraler Ströme. Die Analyse basiert auf einem Vergleich mit Monte-Carlo-Simulationen, deren Produktion global koordiniert wird. In Mainz ist es erstmalig gelungen, Simulationen innerhalb der Architektur des Worldwide LHC Computing Grid (WLCG) zu realisieren, was die Möglichkeit eröffnet, Monte-Carlo-Berechnungen auch auf andere deutsche Rechnerfarmen (CEs) mit IceCube-Berechtigung zu verteilen. Atmosphärische Myonen werden mit einer Rate von über 1000 Ereignissen pro Sekunde aufgezeichnet. Eine korrekte Interpretation dieses dominanten Signals, welches um einen Faktor von 10^6 reduziert werden muss um das eigentliche Neutrinosignal zu extrahieren, ist deswegen von großer Bedeutung. Eigene Simulationen mit der Software-Umgebung CORSIKA wurden durchgeführt um die von Energie und Einfallswinkel abhängige Entstehungshöhe atmosphärischer Myonen zu bestimmen. IceCube Myonraten wurden mit Wetterdaten des European Centre for Medium-Range Weather Forcasts (ECMWF) verglichen und Korrelationen zwischen jahreszeitlichen sowie kurzzeitigen Schwankungen der Atmosphärentemperatur und Myonraten konnten nachgewiesen werden. Zudem wurde eine Suche nach periodischen Effekten in der Atmosphäre, verursacht durch z.B. meteorologische Schwerewellen, mit Hilfe einer Fourieranalyse anhand der IceCube-Daten durchgeführt. Bislang konnte kein signifikanter Nachweis zur Existenz von Schwerewellen am Südpol erbracht werden.

Valutazione della fluenza fotonica di un tubo a raggi X tramite metodi semiempirici basati su simulazioni Montecarlo

Il lavoro svolto in questa tesi si inserisce in una ricerca iniziata pochi anni fa dal Dott. Bontempi. Riguarda la stima della fluenza fotonica tramite due metodi diversi da quelli adottati nella letteratura. I metodi proposti dalla letteratura valutano operativamente la fluenza fotonica ad alta intensità. Sono basati sulla correlazione tra fluenza fotonica e quantità dosimetriche, come l’esposizione o l’Air Kerma. Altre correlano la valutazione dell’HVL con la fluenza dei fotoni. Sebbene queste metodologie siano semplici da implementare, la valutazione della fluenza dei fotoni è ottenuta tramite il concetto di equivalenza del fascio monocromatico. Questo causa discrepanze nella valutazione della reale quantità di fotoni emessa dalla sorgente di raggi X. Il Dott. Bontempi ha sviluppato due diverse metodologie per il calcolo della fluenza: il primo metodo è basato sul calcolo del Kerma unito ad un modello semi empirico per il computo dello spettro dei raggi X. Il secondo metodo è invece basato sulla valutazione della funzione nNPS registrata in un’immagine a raggi X e quella ottenuta simulando il sistema sorgente-rivelatore.