| Geant4 Cross Reference |
1 // ******************************************* 1
2 // To execute this macro under ROOT after your
3 // 1 - launch ROOT (usually type 'root' at y
4 // 2 - type '.X plot.C' at the ROOT session
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 auto c1 = new TCanvas ("c1","",60,60,800,800);
16 c1->Divide(1,1);
17
18 FILE * fp = fopen("mfp.txt","r");
19
20 Float_t E,mfp,smfp;
21 Int_t ncols = 0;
22 Int_t nlines = 0;
23
24 TNtuple *ntuple = new TNtuple("ntuple","mfp","
25 while (1)
26 {
27 ncols = fscanf(fp,"%f %f %f",&E,&mfp,&smfp);
28 if (ncols < 0) break;
29 ntuple->Fill(E,mfp,smfp);
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,1e4,2,0.0
39 h2->Draw();
40 ntuple->SetMarkerStyle(20);
41 ntuple->SetMarkerSize(1.);
42 ntuple->Draw("mfp: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("MFP (nm)");
52 }
53