Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
extended/medical/dna/range/plot.C
이 파일의 문서화 페이지로 가기
1 //
2 // *********************************************************************
3 // To execute this macro under ROOT after your simulation ended,
4 // 1 - launch ROOT (usually type 'root' at your machine's prompt)
5 // 2 - type '.X plot.C' at the ROOT session prompt
6 // 3 - OR directly type 'root plot.C'
7 // *********************************************************************
8 
9 {
10 
11 gROOT->Reset();
12 gStyle->SetPalette(1);
13 gROOT->SetStyle("Plain");
14 gStyle->SetOptStat(00000);
15 
16 c1 = new TCanvas ("c1","Range",60,60,700,700);
17 c1->Divide(1,1);
18 c1->cd(1);
19 gPad->SetLogx();
20 gPad->SetLogy();
21 
22 TH2F * h2 = new TH2F("h2","",2,9.99,1e+4,2,1e-1,1e+4);
23 h2->Draw();
24 h2->GetXaxis()->SetLabelSize(0.025);
25 h2->GetYaxis()->SetLabelSize(0.025);
26 h2->GetXaxis()->SetTitleSize(0.035);
27 h2->GetYaxis()->SetTitleSize(0.035);
28 h2->GetXaxis()->SetTitleOffset(1.4);
29 h2->GetYaxis()->SetTitleOffset(1.4);
30 h2->GetXaxis()->SetTitle("E (eV)");
31 h2->GetYaxis()->SetTitle("Distance (nm)");
32 
33 FILE * fp = fopen("range.txt","r");
34 
38 
39 TNtuple *ntuple = new TNtuple("ntuple","range","e:track:strack:proj:sproj:pene:spene");
40 
41 while (1)
42 {
43  ncols = fscanf(fp,"%f %f %f %f %f %f %f",&e,&track,&strack,&proj,&sproj,&pene,&spene);
44  if (ncols < 0) break;
45  ntuple->Fill(e,track,strack,proj,sproj,pene,spene);
46  nlines++;
47 }
48 
49 fclose(fp);
50 
51 ntuple->SetLineWidth(3);
52 ntuple->SetLineColor(2);
53 ntuple->Draw("track:e","","Lsame");
54 ntuple->SetLineColor(3);
55 ntuple->Draw("pene:e","","Lsame");
56 ntuple->SetLineColor(4);
57 ntuple->Draw("proj:e","","Lsame");
58 }
TTree * ntuple
nlines
G4float Float_t
fclose(fg1)
G4int Int_t
FILE * fp
TH1F * h2
Int_t ncols