LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
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("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","E:mfp:smfp");
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.099,1e2);
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 }
Float_t smfp
Definition: plot.C:20
Float_t mfp
Definition: plot.C:20
TNtuple * ntuple
Definition: plot.C:20
Float_t E
Definition: plot.C:20
FILE * fp
Definition: plot.C:34
TH1F * h2
Definition: plot.C:44
fclose(fp)
c1
Definition: plot.C:27
Int_t nlines
Definition: plot.C:51
Int_t ncols
Definition: plot.C:52