domenica, gennaio 21, 2007

macro di test di ROOT

Io ci metto questa macro, aspetto le altre!


  • Chi posta la lettura di un istogramma da file di testo?
  • Chi posta il fit gaussiano con funzione definita dall'utente, compresa l'impostazione dei parametri iniziali?
  • Chi posta la lettura del TGraph da file di testo?
  • Chi posta il fit del TGraph e il modo di recuperare dall'oggetto TF1, via TGraph, i valori dei parametri del fit?
  • Chi posta la lettura di un TH1F da un file di ROOT?
  • Chi posta la lettura di una TNtuple da un file di root ed il suo disegno?
  • Oppure la produzione e il rimepimento di un TGraphErrors, con ciascun punto determinato da una funzione nota, un po' spostato sia in x che in y con un numero casuale e con un errore determinato altrettanto casualmente, sia in x che in y?
  • ... ed il fit di questo TGraph?


testRootHisto.C



// macro di test di un istogramma ROOT
void testRootHisto ()
{
// creo un nuovo istogramma
TH1F * histoProva = new TH1F ("histoProva","istogramma di test",10,0,10) ;
// nomeOggettto nomeXroot titolo bins,min,max
histoProva->Fill (1) ;
histoProva->Fill (3) ;
histoProva->Fill (4) ;
histoProva->Fill (6.3) ;
histoProva->Fill (6.8) ;
histoProva->Fill (3) ;

// per visualizzare nel box di statistica tutte le statistiche
// relative all'istogramma (media, rms, overfflow, underflow...)
gStyle->SetOptStat (1111111) ;
// per visualizzare nel box di statistica tutte le statistiche
// relative ad eventuali fit fatti all'istogramma
gStyle->SetOptFit (1111) ;

// preparo la cornice dove disegnare l'istogramma

TCanvas * c1 = new TCanvas ;
// imposto opzioni grafiche dell'istogramma: il colore di riempimento ...

histoProva->SetFillColor (8) ;
// ... ed il titolo per l'asse x

histoProva->GetXaxis ()->SetTitle ("asse x") ;
// disegno l'istogramma

histoProva->Draw () ;
histoProva->Fit ("gaus") ;


return ;
}

Nessun commento: