Geant4 Cross Reference |
1 { 2 // For Flexi source 3 // eventually convert from ASCII output files to ROOT files 4 /* 5 gROOT->ProcessLine(".L convert.C"); 6 ReadASCII("EnergyDeposition_Flexi_livermore.out", "brachytherapy_Flexi_livermore.root"); 7 ReadASCII("EnergyDeposition_Flexi_penelope.out", "brachytherapy_Flexi_penelope.root"); 8 ReadASCII("EnergyDeposition_Flexi_opt0.out", "brachytherapy_Flexi_opt0.root"); 9 ReadASCII("EnergyDeposition_Flexi_opt3.out", "brachytherapy_Flexi_opt3.root"); 10 ReadASCII("EnergyDeposition_Flexi_opt4.out", "brachytherapy_Flexi_opt4.root"); 11 */ 12 13 // g(r) is calculated from the root files 14 gROOT->ProcessLine(".L TG43_relative_dose.C"); 15 Read("Flexi", "livermore"); 16 Read("Flexi", "penelope"); 17 Read("Flexi", "opt0"); 18 Read("Flexi", "opt3"); 19 Read("Flexi", "opt4"); 20 //gROOT->ProcessLine(".x compare.C"); 21 // or... 22 // plot all results with the alternative EM physics constructors 23 gROOT->ProcessLine(".x compare_all.C"); 24 // use compare.C when only one physics approach is used 25 26 // For Oncura source 27 /*gROOT->ProcessLine(".L TG43_relative_dose.C"); 28 Read("Oncura", "livermore"); 29 Read("Oncura", "penelope"); 30 Read("Oncura", "opt0"); 31 Read("Oncura", "opt3"); 32 Read("Oncura", "opt4"); 33 gROOT->ProcessLine(".x compare_6711.C"); 34 //or ... 35 // gROOT->ProcessLine(".x compare_6711_all.C"); 36 */ 37 } 38 39