Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // To execute this macro under ROOT after your 4 // 1 - launch ROOT (usually type 'root' at y 5 // 2 - type '.X plot.C' at the ROOT session 6 // 3 - OR directly type 'root plot.C' 7 // ******************************************* 8 9 { 10 gROOT->Reset(); 11 gStyle->SetPalette(1); 12 gROOT->SetStyle("Plain"); 13 gStyle->SetOptStat(00000); 14 15 auto c1 = new TCanvas ("c1","Range",60,60,700, 16 c1->Divide(1,1); 17 18 c1->cd(1); 19 gPad->SetLogx(); 20 gPad->SetLogy(); 21 22 TH2F * h2 = new TH2F("h2","",2,9.99,1e+4,2,1e 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 35 Float_t e,track,strack,proj,sproj,pene,spene; 36 Int_t ncols = 0; 37 Int_t nlines = 0; 38 39 TNtuple *ntuple = new TNtuple("ntuple","range" 40 41 while (1) 42 { 43 ncols = fscanf(fp,"%f %f %f %f %f %f %f",&e, 44 if (ncols < 0) break; 45 ntuple->Fill(e,track,strack,proj,sproj,pene, 46 nlines++; 47 } 48 fclose(fp); 49 50 ntuple->SetLineWidth(3); 51 52 ntuple->SetLineColor(2); 53 ntuple->Draw("track:e","","Lsame"); 54 55 ntuple->SetLineColor(3); 56 ntuple->Draw("pene:e","","Lsame"); 57 58 ntuple->SetLineColor(4); 59 ntuple->Draw("proj:e","","Lsame"); 60 } 61