36 #include "Analysis.hh"
56 fincident_map =
new TH2D(
"incident map",
"Incident Distributuon",
67 fdose_map =
new TH2D(
"dose map",
"Dose Distribution",
74 fdose_hist =
new TH1D(
"dose",
"Dose Distribution", 500, 0., 50.);
111 TFile*
file =
new TFile(fname.c_str(),
"RECREATE",
"Geant4 ROOT analysis");
142 if ( std::abs(p.
y()) < dxy ) {
145 if ( std::abs(p.
x()) < dxy ) {
ntupleExperimental Fill(depthExp, EdepExp)
static constexpr double mm
gr1 GetXaxis() -> SetTitle("Distance from the centre (cm)")
gr1 SetTitle("Dose rate distribution")
void FillIncident(const G4ThreeVector &p)
#define G4MUTEX_INITIALIZER
void Save(const G4String &fname)
ntupleExperimental Reset()
G4ThreadLocal Analysis * the_analysis
gr1 GetYaxis() -> SetTitle("Normalised dose rate distribution")
static Analysis * GetAnalysis()
void FillDose(const G4ThreeVector &p, G4double dedx)
static constexpr double cm
static G4ThreadLocal G4int fincidentFlag
static constexpr double GeV