Geant4 Cross Reference |
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 auto 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 } 54