Root Stefano Lacaprara slides rubate a Tommaso Boccali

Root Stefano Lacaprara (slides rubate a Tommaso Boccali)

Cos'e' l'analisi dati L'analisi dati consiste nell'estrarre da una (decentemente confusa) serie di dati alcune (poche) quantita’ di interesse fisico Visto che tutte le quantita’ misurabili sono affette da rumore (strumentale, stocastico, ecc), e’ necessario adottare procedure per eliminare tale componente

Inoltre…. non sempre le quantita’ misurate sono pure, cioe’ non affette da misurazioni di quantita’ che poco c’entrano con quello che vogliamo misurare per esempio: voglio misurare la temperature di una stanza ma se c’e’ gente dentro la temperatura cambia, per cui dovrei avere modo di capire se la stanza e’ occupata oppure no non sempre posso farlo, devo cercare di capire quale sia la temperatura media della stanza senza saperlo

Misure le 3 componenti di una misura possono essere divise in segnale fondo misure che vengono da fenomeni diversi da quello che vogliamo misurare rumore quello che vorremmo misurare componente aggiuntiva dovuta tipicamente agli strumenti di misura (o alla meccanica quantistica…) Misura = segnale + fondo + rumore

Come estrarre il segnale? serve: una modellazione del comportamento del segnale una modellazione del comportamento del fondo una modellazione del rumore un FIT!

esempio …

posso dire a occhio ma la simuolazione lo deve dconfermare che il segnale e’ gaussiano il fondo e’ circa piatto il rumore c’e’… da quel plot voglio estrarre il valore medio della gaussiana … come fare?

FIT Ci sono molti modi, ma il soldoni sapete gia’cosa si prendono dei dati si prende una funzione matematica con dei parametri liberi si cambiano I parametri fino a quando I dati e la funzione sono il piu’simili possibile e qui sta il casino …

ma se voglio cose piu’complesse? Se non mi basta fare fit a gaussiane ecc? Uso ROOT (o un altro pacchetto specifico di analisi dati) Per esempio: root permette di usare una funzione in C++ (con parametri liberi) per effettuare il fit

Come fare? definire una funzione C++ del tipo double funzione(double * x, double *par){ return par[0]*x[0]+par[1]; } puntatore, consideratelo un vettore delle coordinate se e’una funzione ad una sola variabile, x[0] e’x puntatore, consideratelo un vettore dei parametri liberi

ho quindi definito una funzione ad una variabile (x[0]) che ha due parametri liberi (par[0] e par[1]) in pratica funzione = P 0*x+P 1 sto facendo il fit ad una retta!

a questo punto … devo insegnare a root questa funzione TF 1 funz(“funzione. Mia”, funzione, 0, 100, 2); intervallo di validita’ nome con cui root la conoscera’ nome della funzione che ho creato in C++ numero di parametri liberi

A questo punto isto. Fit(“funzione. Mia”); esegue il fit… facciamo un tentativo completo!

Fisica delle alte energie ho eseguito una misura in fisica delle alte energie, devo estrarre il valore della massa di una particella, e la sua vita media posso modellizzare il segnale dalla particella come una gaussiana, il cui valore medio e’la massa e la cui sigma e’ legata alla vita media la misura ha del fondo dovuto a fenomeni che al momento non ci interessano, ma il fondo posso supporlo constante (a meno del rumore)

ho generato I dati per voi… come sovrapposizione di una gaussiana, il segnale e di una costante wget http: //www. pi. infn. it/~boccali/Milano/Massa/dati. dat lo ho creato con wget http: //www. pi. infn. it/~boccali/Milano/Massa/file. c

Root … Fit di funzioni semplici gaussiane rette parabole ecc Stupido, dato un istogramma, aprire il Fit Panel oppure isto. Fit(“gaus”)

ora… facciamo un programma che legga questi dati, come se venissero dallo strumento di misura li metta in un istogramma defininiamo una funzione gaussiana+costate, quindi con 3+1 parametri liberi facciamo il fit…

Una volta fatto il fit… possiamo avere in C++ il risultato, con funz. Get. Parameter(2); // legge par[2] funz. Get. Par. Error(2); // legge l’errore su par[2]

prego … provare

ha funzionato? probabilmente no … il fatto e’ che la procedura di fit non puo’essere piu’di tanto furba, per cui devo dare di valori di partenza “non inverosimili” Prima di fare il fit, posso fare funz. Set. Parameter(2, 66); inizializza il par[2] a 66 prima di fare il fit

se la sigma viene negativa? possiamo dare dei limiti ai parametri liberi per esempio funz. Set. Par. Limits(3, 0, 100); impongo ch par[3] possa andare da 0 a 100 solamente