Ouvrir un fichier ROOT et afficher ses histogrammes.

Il existe deux méthodes pour ouvrir un fichier root et afficher son contenu: la ligne de commande ou l'interface graphique.

A la ligne de commande:

Pour ouvrir le fichier "hsimple.root", il faut taper:

 TFile *f=new TFile("hsimple.root")

Pour visualiser sont contenu, il faut taper:

f->ls()

On voit la liste des objets contenus dans le fichier comme dans cet exemple:

TFile**         hsimple.root    Demo ROOT file with histograms
 TFile*         hsimple.root    Demo ROOT file with histograms
  KEY: TH1F     hpx;1   This is the px distribution
  KEY: TH2F     hpxpy;1 py vs px
  KEY: TProfile hprof;1 Profile of pz versus px
  KEY: TNtuple  ntuple;1        Demo ntuple

Pour afficher l'histogramme 1D (TH1F) nommé hpx, il suffit de taper:

f->Get("hpx")->Draw()

ou plus simplement

hpx->Draw()

La deuxième méthode est plus rapide, mais moins sure que la première.

Pour visualiser les noms des variables contenues dans le Ntuple nommé "ntuple", il faut taper:

f->Get("ntuple")->Print()

ou plus simplement

ntuple->Print()

La deuxième méthode est plus rapide, mais moins sure que la première. Le résultat doit ressembler à ce qui suit:

******************************************************************************
*Tree    :ntuple    : Demo ntuple                                            *
*Entries :    25000 : Total =          504510 bytes  File  Size =     399823 *
*        :          : Tree compression factor =   1.26                       *
******************************************************************************
*Br    0 :px        :                                                        *
*Entries :    25000 : Total  Size=     100835 bytes  File Size  =      89056 *
*Baskets :        3 : Basket Size=      32000 bytes  Compression=   1.08     *
*............................................................................*
*Br    1 :py        :                                                        *
*Entries :    25000 : Total  Size=     100835 bytes  File Size  =      89061 *
*Baskets :        3 : Basket Size=      32000 bytes  Compression=   1.08     *
*............................................................................*
*Br    2 :pz        :                                                        *
*Entries :    25000 : Total  Size=     100835 bytes  File Size  =      87525 *
*Baskets :        3 : Basket Size=      32000 bytes  Compression=   1.10     *
*............................................................................*
*Br    3 :random    :                                                        *
*Entries :    25000 : Total  Size=     100871 bytes  File Size  =      85475 *
*Baskets :        3 : Basket Size=      32000 bytes  Compression=   1.12     *
*............................................................................*
*Br    4 :i         :                                                        *
*Entries :    25000 : Total  Size=     100826 bytes  File Size  =      30093 *
*Baskets :        3 : Basket Size=      32000 bytes  Compression=   3.19     *
*............................................................................*

Pour afficher un histogramme 1D de la variable "px" avec le ntuple, il faut taper:

f->Get("ntuple")->Draw("px")

ou plus simplement

ntuple->Draw("px")

Pour afficher un histogramme 2D des variables "px" et "py" pour les valeurs de "pz" inferieures à 2, il faut taper:

f->Get("ntuple")->Draw("py:px","pz < 2","box")

ou plus simplement

ntuple->Draw("py:px","pz < 2","box")

Pour plus de renseignements sur les NTuples, ou plus exactement les objets TNtuple et TTree, on pourra consulter les pages suivantes:

A l'interface graphique:

Il faut d'abord lancer un navigateur en tapant:

TBrowser b

Vous devrier voir une fenetre ressemblant a ceci:

Allez ensuite dans le menu "Fichier" puis activer l'item "Open...". Double cliquer ensuite sur le fichier voulu.

Cliquer ensuite dans le dossier "ROOT files" dans la colonne de gauche du navigateur, puis sélectionnez le fichier désiré en cliquant dessus. Normalement, le contenu du fichier s'affiche dans la fenetre de droite.

Pour afficher un histogramme, il suffit de double-cliquer dessus. Pour le Ntuple, pointez-le et cliquez sur le bouton droit. Sélectionnez l'item "StartViewer".

La suite des opérations est décrite dans comment afficher le contenu d'un TTree (TNtuple) avec le TreeViewer.