Go to the source code of this file.
|
| while (1) |
|
| fclose (fg1) |
|
| fclose (fg2) |
|
| fclose (fg3) |
|
| fclose (fg4) |
|
| fclose (fg5) |
|
| fclose (fg6) |
|
std::cout<< "Livermore"<< std::endl;for(Int_t j=0;j< nlines1;j++){std::cout<< x1[j]<< ", "<< ratio_liv[j]<< std::endl;}std::cout<< "penelope"<< std::endl;for(Int_t j=0;j< nlines1;j++){std::cout<< x1[j]<< ", "<< ratio_pen[j]<< std::endl;}std::cout<< "opt0"<< std::endl;for(Int_t j=0;j< nlines1;j++){std::cout<< x1[j]<< ", "<< ratio_opt0[j]<< std::endl;}std::cout<< "opt3"<< std::endl;for(Int_t j=0;j< nlines1;j++){std::cout<< x1[j]<< ", "<< ratio_opt3[j]<< std::endl;}std::cout<< "opt4"<< std::endl;for(Int_t j=0;j< nlines1;j++){std::cout<< x1[j]<< ", "<< ratio_opt4[j]<< std::endl;}TCanvas *c1=new TCanvas("c1","Graph Draw Options", 200, 10, 600, 400);gPad-> | SetLogy () |
|
gr1 | SetTitle ("Dose rate distribution") |
|
gr1 | GetXaxis () -> SetTitle("Distance from the centre (cm)") |
|
gr1 | GetYaxis () -> SetTitle("Normalised dose rate distribution") |
|
gr1 | SetLineWidth (1) |
|
gr1 | SetMarkerColor (1) |
|
gr1 | SetMarkerStyle (20) |
|
gr1 | Draw ("AP") |
|
gr2 | SetMarkerColor (2) |
|
gr2 | SetMarkerStyle (21) |
|
gr2 | SetMarkerSize (0.5) |
|
gr2 | SetLineColor (2) |
|
gr2 | Draw ("CP") |
|
gr3 | SetLineWidth (0.3) |
|
gr3 | SetMarkerColor (3) |
|
gr3 | SetLineColor (3) |
|
gr4 | SetMarkerColor (4) |
|
gr4 | SetLineColor (4) |
|
gr5 | SetMarkerColor (6) |
|
gr5 | SetLineColor (6) |
|
gr6 | SetMarkerColor (8) |
|
gr6 | SetLineColor (8) |
|
leg | SetFillColor (0) |
|
leg | AddEntry (gr1,"Reference data","lp") |
|
leg | AddEntry (gr2,"Geant4 - Oncura - Livermore","lp") |
|
leg | AddEntry (gr3,"Geant4 - Penelope","lp") |
|
leg | AddEntry (gr4,"Geant4 - Standard opt0","lp") |
|
leg | AddEntry (gr5,"Geant4 - Standard opt3","lp") |
|
leg | AddEntry (gr6,"Geant4 - Standard opt4","lp") |
|
leg | Draw () |
|
|
Int_t | n_points_dolan = 8 |
|
Float_t | x1 [n_points_dolan] |
|
Float_t | y1 [n_points_dolan] |
|
Float_t | x |
|
Float_t | y |
|
Int_t | ncols_dolan |
|
Int_t | nlines1 = 0 |
|
FILE * | fg2 =fopen("geant4_dose_Oncura_livermore.txt", "r") |
|
Int_t | n_points_geant4 = 398 |
|
Float_t | x2 [n_points_geant4] |
|
Float_t | y2 [n_points_geant4] |
|
Float_t | ratio_liv [n_points_dolan] |
|
Int_t | ncols_geant4 |
|
Int_t | nlines2 = 0 |
|
FILE * | fg3 =fopen("geant4_dose_Oncura_penelope.txt", "r") |
|
Float_t | x3 [n_points_geant4] |
|
Float_t | y3 [n_points_geant4] |
|
Float_t | ratio_pen [n_points_dolan] |
|
Int_t | ncols_geant4_penelope |
|
Int_t | nlines3 =0 |
|
FILE * | fg4 =fopen("geant4_dose_Oncura_opt0.txt", "r") |
|
Float_t | x4 [n_points_geant4] |
|
Float_t | y4 [n_points_geant4] |
|
Float_t | ratio_opt0 [n_points_dolan] |
|
Int_t | ncols_geant4_opt0 |
|
Int_t | nlines4 =0 |
|
FILE * | fg5 =fopen("geant4_dose_Oncura_opt3.txt", "r") |
|
Float_t | x5 [n_points_geant4] |
|
Float_t | y5 [n_points_geant4] |
|
Float_t | ratio_opt3 [n_points_dolan] |
|
Int_t | ncols_geant4_opt3 |
|
Int_t | nlines5 =0 |
|
FILE * | fg6 =fopen("geant4_dose_Oncura_opt4.txt", "r") |
|
Float_t | x6 [n_points_geant4] |
|
Float_t | y6 [n_points_geant4] |
|
Float_t | ratio_opt4 [n_points_dolan] |
|
Int_t | ncols_geant4_opt4 |
|
Int_t | nlines6 =0 |
|
TGraph * | gr1 = new TGraph (nlines1, x1, y1) |
|
TGraph * | gr2 = new TGraph (nlines2, x2, y2) |
|
TGraph * | gr3 = new TGraph (nlines3, x3, y3) |
|
TGraph * | gr4 = new TGraph (nlines4, x4, y4) |
|
TGraph * | gr5 = new TGraph (nlines5, x5, y5) |
|
TGraph * | gr6 = new TGraph (nlines6, x6, y6) |
|
TLegend * | leg = new TLegend(0.3, 0.5, 0.6, 0.8) |
|
leg AddEntry |
( |
gr1 |
, |
|
|
"Reference data" |
, |
|
|
"lp" |
|
|
) |
| |
leg AddEntry |
( |
gr2 |
, |
|
|
"Geant4 - Oncura - Livermore" |
, |
|
|
"lp" |
|
|
) |
| |
leg AddEntry |
( |
gr3 |
, |
|
|
"Geant4 - Penelope" |
, |
|
|
"lp" |
|
|
) |
| |
leg AddEntry |
( |
gr4 |
, |
|
|
"Geant4 - Standard opt0" |
, |
|
|
"lp" |
|
|
) |
| |
leg AddEntry |
( |
gr5 |
, |
|
|
"Geant4 - Standard opt3" |
, |
|
|
"lp" |
|
|
) |
| |
leg AddEntry |
( |
gr6 |
, |
|
|
"Geant4 - Standard opt4" |
, |
|
|
"lp" |
|
|
) |
| |
gr1 GetXaxis |
( |
| ) |
-> SetTitle("Distance from the centre (cm)") |
gr5_ratio SetLineColor |
( |
6 |
| ) |
|
gr5_ratio SetLineWidth |
( |
0. |
3 | ) |
|
std::cout<< "Livermore" << std::endl;for (Int_t j=0; j < nlines1; j++){ std::cout << x1[j] << ", " << ratio_liv[j] << std::endl;} std::cout<< "penelope" << std::endl;for (Int_t j=0; j < nlines1; j++) { std::cout << x1[j] << ", " << ratio_pen[j] << std::endl;} std::cout<< "opt0" << std::endl;for (Int_t j=0; j < nlines1; j++) { std::cout << x1[j] << ", " << ratio_opt0[j] << std::endl;} std::cout<< "opt3" << std::endl;for (Int_t j=0; j < nlines1; j++) { std::cout << x1[j] << ", " << ratio_opt3[j] << std::endl;} std::cout<< "opt4" << std::endl;for (Int_t j=0; j < nlines1; j++) { std::cout << x1[j] << ", " << ratio_opt4[j] << std::endl;} TCanvas *c1 = new TCanvas("c1", "Graph Draw Options", 200, 10, 600, 400);gPad-> SetLogy |
( |
| ) |
|
gr3_ratio SetMarkerColor |
( |
3 |
| ) |
|
gr5_ratio SetMarkerColor |
( |
6 |
| ) |
|
gr2 SetMarkerSize |
( |
0. |
5 | ) |
|
gr1 SetMarkerStyle |
( |
20 |
| ) |
|
gr2 SetMarkerStyle |
( |
21 |
| ) |
|
gr1 SetTitle |
( |
"Dose rate distribution" |
| ) |
|
FILE* fg2 =fopen("geant4_dose_Oncura_livermore.txt", "r") |
FILE* fg3 =fopen("geant4_dose_Oncura_penelope.txt", "r") |
FILE* fg4 =fopen("geant4_dose_Oncura_opt0.txt", "r") |
FILE* fg5 =fopen("geant4_dose_Oncura_opt3.txt", "r") |
FILE* fg6 =fopen("geant4_dose_Oncura_opt4.txt", "r") |
TLegend* leg = new TLegend(0.3, 0.5, 0.6, 0.8) |
Int_t n_points_geant4 = 398 |
Int_t ncols_geant4_penelope |