Sirri, Gabriele

Almae Matris Studiorum Campus

Analisi Statistica dei dati nella Fisica Nucl. e Subnucl. [Laboratorio ]Gabriele SirriIstituto Nazionale di Fisica Nucleare2015.05.18 CalendarioLaboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 2• Lunedì 23 febbraio 2015 14-16 M. SioliMARZO• Lunedì 2 marzo 2015 14-16 M. SioliGiovedì 5 marzo 2015 11-13 T. Chiarusi• Lunedì 9 marzo 2015 14-16 M. SioliGiovedì 12 marzo 2015 11-13 M. Sioli• Lunedì 16 marzo 2015 14-16 M. SioliGiovedì 19 marzo 2014 11-13 T. Chiarusi• Lunedì 23 marzo 2015 14-16 M. SioliGiovedì 26 marzo 2015 11-13 M. Sioli Giovedì 26 marzo 2015 16-18 G. Sirri• Lunedì 30 marzo 2015 14-16 M. SioliAPRILE• Mercol. 8 aprile 2015 10-13 M. Sioli/T.Chiar.• Giovedì 9 aprile 2015 11-13 G. Sirri• Lunedì 13 aprile 2015 14-16 M. Sioli Giovedì 16 aprile 2015 11-13 T. ChiarusiGiovedì 16 aprile 2015 16-18 G. Sirri • Lunedì 20 aprile 2015 14-16 M. Sioli Giovedì 23 aprile 2015 11-12 G. SirriGiovedì 23 aprile 2015 12-13 T. Chiarusi • Giovedì 30 aprile 2015 11-13 G. Sirri Giovedì 30 aprile 2015 16-18 G. SirriMAGGIO• Lunedì 4 maggio 2015 14-16 M. Sioli Giovedì 7 maggio 2015 11-13 T. Chiarusi • Lunedì 11 maggio 2015 14-16 M. Sioli Giovedì 14 maggio 2015 11-13 M. SioliGiovedì 14 maggio 2015 16-18 T. Chiarusi• Lunedì 18 maggio 2015 14-16 G. SirriGiovedì 21 maggio 2015 11-13 G. Sirri Giovedì 21 maggio 2015 14-16 M. Sioli Tutte le lezioni in Aula C, via Irnerio2015.05.182015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 3 Esercitazioni (soluzioni)Le soluzioni di Es. 3 e 4 NON SARANNO PUBBLICATE su AMSCampus, finchè non mi saranno inviate la soluzioni da tutti.Chi volesse può chiederle privatamente via mail.accesso riservato agli iscritti a gabriele.sirri2.ASD-2015con password (richiedetela via mail)Homework2015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 42015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 5Preparazione Es.6 per giovedìSeguite il Tutorial :https://twiki.cern.ch/twiki/bin/view/RooStats/RooStatsTutorialsAugust2012Costruite una macro per la creazione di un modello per un conteggio poissoniano di segnale piu’ fondo assumendo che il fondo sia distribuito secondo gauss.Non c’è bisogno di spedire la soluzione via mail. Completeremo l’esercizio giovedi’.Exercise time (Esercizio 5) 2015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 6ROOSTATS : https://twiki.cern.ch/twiki/bin/view/RooStatsshort tutorial: https://twiki.cern.ch/twiki/bin/view/RooStats/RooStatsTutorialsAugust2012Esercizio 5 - roostats_ex1[1] roostats_ex1.C Riprendiamo l’esercizio della lezione precedente(gaussian signal + exponential background, extended p.d.f.) .Trovate un template già fatto qui :  http://campus.unibo.it/186979/In makemodel : c’è la costruzione del modello , la generazione di un set di  1000 dati simulato. Il workspace è salvato in «model.root».Alcune modifiche rispetto alla lezione precedente :Range e Valori iniziali dei parametri sono stati modificati.«mean» e «sigma» sono fissati come costanti.I dati sono importati nel workspace e salvati su file.In usemodel : legge il modello dal file, esegue il fit del modello e disegna dati e risultato del fitin un plot.2015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 7ROOSTATS : https://twiki.cern.ch/twiki/bin/view/RooStatsshort tutorial: https://twiki.cern.ch/twiki/bin/view/RooStats/RooStatsTutorialsAugust2012Si modifichi il modello:- Specificare i componenti del modello per i tool statistici di roostat:  osservabile e parametro di interesse. - Utilizzare il numero di eventi di segnale come unico parametro di interesse. - Fissare costanti tutti gli altri parametri del modello.- Importare la configurazione nel workspace e salvare su file.Si modifichi l’uso del modello:- leggere il modelConfig dal workspace esempio: ModelConfig* mc = (ModelConfig*) w.obj("ModelConfig");- calcolare un Confidence Interval utilizzando il ProfileLikelihoodCalculator- Disegnare il profilo della likelihood e sovrapporre l’intervallo- calcolare la discovery significance utilizzando il profilelikelihoodcalculatorcome test di ipotesi- scrivere sulla console i limiti dell’intervallo e la significativitàEsercizio 5 - roostats_ex1ROOSTATS : https://twiki.cern.ch/twiki/bin/view/RooStatsshort tutorial: https://twiki.cern.ch/twiki/bin/view/RooStats/RooStatsTutorialsAugust20122015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 8[2] roostats_ex2.CAggiungere l’intervallo calcolato con Feldman-Cousinsuggerimento : modificare solo usemodel()guardare il codice in :  $ROOTSYS/tuturials/roostats/IntervalExamples.C[3] roostats_ex3.CDefinire tau e Nb come nuisance parametersripetere i test dell’ex. 1 e confrontarlisuggerimento modificare makemodel()  (vedere il tutorial) Esercizio 5 - roostats_ex2 , roostats_ex3ROOSTATS : https://twiki.cern.ch/twiki/bin/view/RooStatsshort tutorial: https://twiki.cern.ch/twiki/bin/view/RooStats/RooStatsTutorialsAugust20122015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 92015.05.18 Laboratorio Analisi Statistica dei Dati per HEP - G.Sirri - 2015 10Suggerimento:usate questocostruttoreAnziché questohttps://root.cern.ch/root/html/RooStats__ProfileLikelihoodCalculator.htmlGuardate la Class ReferenceOvvero :passategli come argomento il ModelConfigenon singolarmenteil Modello e il POI

[7] RooStats 2

http://campus.unibo.it/191676/1/ASD-lab-2015-lesson-7.pdf

[7] RooStats 2

Abstract

Similar works

Full text

Available Versions

Almae Matris Studiorum Campus