| Geant4 Cross Reference |
1 { 1 {
2 // Read reference data in granero.txt 2 // Read reference data in granero.txt
3 FILE *fg1=fopen("granero.txt", "r"); 3 FILE *fg1=fopen("granero.txt", "r");
4 Int_t n_points_granero =13; 4 Int_t n_points_granero =13;
5 Float_t x1[n_points_granero], y1[n_points_gran 5 Float_t x1[n_points_granero], y1[n_points_granero];
6 Float_t x, y; 6 Float_t x, y;
7 Int_t ncols_granero; 7 Int_t ncols_granero;
8 Int_t nlines1 =0; 8 Int_t nlines1 =0;
9 9
10 while(1) 10 while(1)
11 { 11 {
12 ncols_granero = fscanf(fg1,"%f %f",&x, &y); 12 ncols_granero = fscanf(fg1,"%f %f",&x, &y);
13 if (ncols_granero<0) break; 13 if (ncols_granero<0) break;
14 // std::cout << "x " << x << std::endl; 14 // std::cout << "x " << x << std::endl;
15 x1[nlines1]=x; 15 x1[nlines1]=x;
16 y1[nlines1]=y; 16 y1[nlines1]=y;
17 nlines1++; 17 nlines1++;
18 } 18 }
19 19
20 fclose(fg1); 20 fclose(fg1);
21 21
22 // Read the results of the brachytherapy advan 22 // Read the results of the brachytherapy advanced example
23 // FlexiSorceMacro.mac with 280 M events 23 // FlexiSorceMacro.mac with 280 M events
24 FILE *fg2=fopen("geant4.txt", "r"); 24 FILE *fg2=fopen("geant4.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 // std::cout << "x " << x << std::endl; 34 // std::cout << "x " << x << std::endl;
35 x2[nlines2]=x; 35 x2[nlines2]=x;
36 y2[nlines2]=y; 36 y2[nlines2]=y;
37 nlines2++; 37 nlines2++;
38 } 38 }
39 39
40 fclose(fg2); 40 fclose(fg2);
41 41
42 TGraph *gr1 = new TGraph (nlines1, x1, y1); 42 TGraph *gr1 = new TGraph (nlines1, x1, y1);
43 TGraph *gr2 = new TGraph (nlines2, x2, y2); 43 TGraph *gr2 = new TGraph (nlines2, x2, y2);
44 44
45 TCanvas *c1 = new TCanvas("c1","Graph Draw Opt 45 TCanvas *c1 = new TCanvas("c1","Graph Draw Options",
46 200,10,600,400); 46 200,10,600,400);
47 47
48 gPad->SetLogy(); 48 gPad->SetLogy();
49 49
50 // draw the graph with axis, continuous line, 50 // draw the graph with axis, continuous line, and put
51 // a * at each point 51 // a * at each point
52 gr1->SetTitle("Dose rate distribution"); 52 gr1->SetTitle("Dose rate distribution");
53 gr1-> GetXaxis()->SetTitle("Distance from the 53 gr1-> GetXaxis()->SetTitle("Distance from the centre (cm)");
54 gr1->GetYaxis()->SetTitle("Normalised dose rat 54 gr1->GetYaxis()->SetTitle("Normalised dose rate distribution");
55 gr1->SetLineWidth(1); 55 gr1->SetLineWidth(1);
56 gr1->SetMarkerColor(1); 56 gr1->SetMarkerColor(1);
57 gr1->SetMarkerStyle(20); 57 gr1->SetMarkerStyle(20);
58 gr1->Draw("AP"); 58 gr1->Draw("AP");
59 59
60 gr2->SetLineWidth(1); 60 gr2->SetLineWidth(1);
61 gr2->SetMarkerColor(2); 61 gr2->SetMarkerColor(2);
62 gr2->SetMarkerStyle(21); 62 gr2->SetMarkerStyle(21);
63 gr2->SetMarkerSize(0.5); 63 gr2->SetMarkerSize(0.5);
64 gr2->SetLineColor(2); 64 gr2->SetLineColor(2);
65 gr2->Draw("CP"); 65 gr2->Draw("CP");
66 66
67 TLegend *leg = new TLegend(0.3, 0.5, 0.6, 0.8) 67 TLegend *leg = new TLegend(0.3, 0.5, 0.6, 0.8);
68 leg->SetFillColor(0); 68 leg->SetFillColor(0);
69 leg->AddEntry(gr1, "Reference data", "lp"); 69 leg->AddEntry(gr1, "Reference data", "lp");
70 leg->AddEntry(gr2, "Geant4 - 280 M events", "l 70 leg->AddEntry(gr2, "Geant4 - 280 M events", "lp");
71 leg->Draw(); 71 leg->Draw();
72 72
73 73
74 } 74 }
75 75