| Geant4 Cross Reference |
1 { 1 {
2 // Read reference data in dolan.txt 2 // Read reference data in dolan.txt
3 FILE * fg1=fopen("dolan.txt", "r"); 3 FILE * fg1=fopen("dolan.txt", "r");
4 Int_t n_points_dolan = 8; 4 Int_t n_points_dolan = 8;
5 Float_t x1[n_points_dolan], y1[n_points_dolan] 5 Float_t x1[n_points_dolan], y1[n_points_dolan];
6 Float_t x, y; 6 Float_t x, y;
7 Int_t ncols_dolan; 7 Int_t ncols_dolan;
8 Int_t nlines1 = 0; 8 Int_t nlines1 = 0;
9 9
10 while(1) 10 while(1)
11 { 11 {
12 ncols_dolan = fscanf(fg1, "%f %f", &x, &y); 12 ncols_dolan = fscanf(fg1, "%f %f", &x, &y);
13 if (ncols_dolan<0) break; 13 if (ncols_dolan<0) break;
14 x1[nlines1]= x; 14 x1[nlines1]= x;
15 y1[nlines1] = y; 15 y1[nlines1] = y;
16 nlines1++; 16 nlines1++;
17 } 17 }
18 18
19 fclose(fg1); 19 fclose(fg1);
20 20
21 // Read the results of the brachyadvanced exam 21 // Read the results of the brachyadvanced example
22 // OncuraSourceMacro.mac with 20 B events 22 // OncuraSourceMacro.mac with 20 B events
23 23
24 FILE *fg2=fopen("geant4_6711_dose.txt", "r"); 24 FILE *fg2=fopen("geant4_6711_dose.txt", "r");
25 Int_t n_points_geant4 = 398; 25 Int_t n_points_geant4 = 398;
26 Float_t x2[n_points_geant4], y2[n_points_geant 26 Float_t x2[n_points_geant4], y2[n_points_geant4];
27 Int_t ncols_geant4; 27 Int_t ncols_geant4;
28 Int_t nlines2 = 0; 28 Int_t nlines2 = 0;
29 29
30 while(1) 30 while(1)
31 { 31 {
32 ncols_geant4 = fscanf(fg2, "%f %f", &x, &y); 32 ncols_geant4 = fscanf(fg2, "%f %f", &x, &y);
33 if (ncols_geant4<0) break; 33 if (ncols_geant4<0) break;
34 x2[nlines2] = x; 34 x2[nlines2] = x;
35 y2[nlines2] = y; 35 y2[nlines2] = y;
36 nlines2++; 36 nlines2++;
37 } 37 }
38 38
39 fclose(fg2); 39 fclose(fg2);
40 40
41 TGraph *gr1 = new TGraph (nlines1, x1, y1); 41 TGraph *gr1 = new TGraph (nlines1, x1, y1);
42 TGraph *gr2 = new TGraph (nlines2, x2, y2); 42 TGraph *gr2 = new TGraph (nlines2, x2, y2);
43 43
44 TCanvas *c1 = new TCanvas("c1", "Graph Draw Op 44 TCanvas *c1 = new TCanvas("c1", "Graph Draw Options", 200, 10, 600, 400);
45 45
46 gPad->SetLogy(); 46 gPad->SetLogy();
47 47
48 // Draw the graph with axis, continuous line, 48 // Draw the graph with axis, continuous line, and put a '*' at each point
49 49
50 gr1->SetTitle("Dose rate distribution"); 50 gr1->SetTitle("Dose rate distribution");
51 gr1->GetXaxis()->SetTitle("Distance from the c 51 gr1->GetXaxis()->SetTitle("Distance from the centre (cm)");
52 gr1->GetYaxis()->SetTitle("Normalised dose rat 52 gr1->GetYaxis()->SetTitle("Normalised dose rate distribution");
53 gr1->SetLineWidth(1); 53 gr1->SetLineWidth(1);
54 gr1->SetMarkerColor(1); 54 gr1->SetMarkerColor(1);
55 gr1->SetMarkerStyle(20); 55 gr1->SetMarkerStyle(20);
56 gr1->Draw("AP"); 56 gr1->Draw("AP");
57 57
58 gr2->SetLineWidth(1); 58 gr2->SetLineWidth(1);
59 gr2->SetMarkerColor(2); 59 gr2->SetMarkerColor(2);
60 gr2->SetMarkerStyle(21); 60 gr2->SetMarkerStyle(21);
61 gr2->SetMarkerSize(0.5); 61 gr2->SetMarkerSize(0.5);
62 gr2->SetLineColor(2); 62 gr2->SetLineColor(2);
63 gr2->Draw("CP"); 63 gr2->Draw("CP");
64 64
65 TLegend *leg = new TLegend(0.3, 0.5, 0.6, 0.8) 65 TLegend *leg = new TLegend(0.3, 0.5, 0.6, 0.8);
66 leg->SetFillColor(0); 66 leg->SetFillColor(0);
67 leg->AddEntry(gr1, "Reference data", "lp"); 67 leg->AddEntry(gr1, "Reference data", "lp");
68 leg->AddEntry(gr2, "Geant4 - 20B Events", "lp" 68 leg->AddEntry(gr2, "Geant4 - 20B Events", "lp");
69 leg->Draw(); 69 leg->Draw();
70 70
71 } 71 }
72 72
73 73
74 74