plot.C
Go to the documentation of this file.
1 // *********************************************************************
2 // To execute this macro under ROOT after your simulation ended,
3 // 1 - launch ROOT (usually type 'root' at your machine's prompt)
4 // 2 - type '.X plot.C' at the ROOT session prompt
5 // 3 - OR type directly 'root plot.C'
6 // *********************************************************************
7 
8 {
9 gROOT->Reset();
10 
11 gStyle->SetPalette(1);
12 gROOT->SetStyle("Plain");
13 gStyle->SetOptStat(000000);
14 
15 c1 = new TCanvas ("c1","",60,60,800,800);
16 c1->Divide(1,1);
17 
18 FILE * fp = fopen("wvalue.txt","r");
19 
20 Float_t E,nbioni,snbioni,w,sw;
21 Int_t ncols = 0;
22 Int_t nlines = 0;
23 
24 TNtuple *ntuple = new TNtuple("ntuple","w","E:nbioni:snbioni:w:sw");
25 while (1)
26 {
27  ncols = fscanf(fp,"%f %f %f %f %f",&E,&nbioni,&snbioni,&w,&sw);
28  if (ncols < 0) break;
29  ntuple->Fill(E,nbioni,snbioni,w,sw);
30  nlines++;
31 }
32 fclose(fp);
33 
34 c1->cd(1);
35 gPad->SetLogx();
36 gPad->SetLogy();
37 
38 TH2F * h2 = new TH2F ("h2","",2,9.99,1e3,2,9.99,1e3);
39 h2->Draw();
40 ntuple->SetMarkerStyle(20);
41 ntuple->SetMarkerSize(1.);
42 ntuple->Draw("w:E","","LPsame");
43 
44 h2->GetXaxis()->SetLabelSize(0.025);
45 h2->GetYaxis()->SetLabelSize(0.025);
46 h2->GetXaxis()->SetTitleSize(0.035);
47 h2->GetYaxis()->SetTitleSize(0.035);
48 h2->GetXaxis()->SetTitleOffset(1.4);
49 h2->GetYaxis()->SetTitleOffset(1.4);
50 h2->GetXaxis()->SetTitle("E (eV)");
51 h2->GetYaxis()->SetTitle("W (eV)");
52 
53 }
Float_t nbioni
Definition: plot.C:20
TNtuple * ntuple
Definition: plot.C:20
Float_t snbioni
Definition: plot.C:20
Float_t E
Definition: plot.C:20
FILE * fp
Definition: plot.C:35
TH1F * h2
Definition: plot.C:45
fclose(fp)
Float_t sw
Definition: plot.C:20
c1
Definition: plot.C:28
Int_t nlines
Definition: plot.C:52
Float_t w
Definition: plot.C:20
Int_t ncols
Definition: plot.C:53