37 #include "Analysis.hh"
58 fincident_x_hist(0), fincident_map(0), fdose_hist(0), fdose_map(0),
59 fdose_prof(0), fdose_map_prof(0), fdose_map3d(0)
72 5 *
cm, 50, -5 *
cm, 5 *
cm,
"cm",
"cm");
74 = mgr->
CreateH1(
"dose",
"Dose distribution", 500, 0, 50 *
cm,
"cm");
76 200, -10 *
cm, 10 *
cm,
"cm",
"cm");
78 50 *
cm, 20, -10 *
cm, 10 *
cm, 20, -10 *
cm, 10 *
cm,
"cm",
"cm",
"cm");
80 0, 100 *
MeV,
"cm",
"MeV");
82 30 *
cm, 80, -4 *
cm, 4 *
cm, 0, 100 *
MeV,
"cm",
"cm",
"MeV");
142 if (std::abs(p.
y()) < dxy)
150 if (std::abs(p.
x()) < dxy)
G4int CreateH3(const G4String &name, const G4String &title, G4int nxbins, G4double xmin, G4double xmax, G4int nybins, G4double ymin, G4double ymax, G4int nzbins, G4double zmin, G4double zmax, const G4String &xunitName="none", const G4String &yunitName="none", const G4String &zunitName="none", const G4String &xfcnName="none", const G4String &yfcnName="none", const G4String &zfcnName="none", const G4String &xbinSchemeName="linear", const G4String &ybinSchemeName="linear", const G4String &zbinSchemeName="linear")
G4bool FillH2(G4int id, G4double xvalue, G4double yvalue, G4double weight=1.0)
static constexpr double MeV
static constexpr double mm
G4int CreateH2(const G4String &name, const G4String &title, G4int nxbins, G4double xmin, G4double xmax, G4int nybins, G4double ymin, G4double ymax, const G4String &xunitName="none", const G4String &yunitName="none", const G4String &xfcnName="none", const G4String &yfcnName="none", const G4String &xbinSchemeName="linear", const G4String &ybinSchemeName="linear")
G4bool OpenFile(const G4String &fileName="")
G4bool FillP1(G4int id, G4double xvalue, G4double yvalue, G4double weight=1.0)
G4bool FillH3(G4int id, G4double xvalue, G4double yvalue, G4double zvalue, G4double weight=1.0)
G4int CreateH1(const G4String &name, const G4String &title, G4int nbins, G4double xmin, G4double xmax, const G4String &unitName="none", const G4String &fcnName="none", const G4String &binSchemeName="linear")
void FillIncident(const G4ThreeVector &p)
void Save(const G4String &fname)
G4ThreadLocal Analysis * the_analysis
static Analysis * GetAnalysis()
void FillDose(const G4ThreeVector &p, G4double dedx)
G4int CreateP2(const G4String &name, const G4String &title, G4int nxbins, G4double xmin, G4double xmax, G4int nybins, G4double ymin, G4double ymax, G4double zmin=0, G4double zmax=0, const G4String &xunitName="none", const G4String &yunitName="none", const G4String &zunitName="none", const G4String &xfcnName="none", const G4String &yfcnName="none", const G4String &zfcnName="none", const G4String &xbinSchemeName="linear", const G4String &ybinSchemeName="linear")
static constexpr double cm
static G4ThreadLocal G4int fincidentFlag
G4bool FillP2(G4int id, G4double xvalue, G4double yvalue, G4double zvalue, G4double weight=1.0)
G4bool FillH1(G4int id, G4double value, G4double weight=1.0)
static constexpr double GeV
G4int CreateP1(const G4String &name, const G4String &title, G4int nbins, G4double xmin, G4double xmax, G4double ymin=0, G4double ymax=0, const G4String &xunitName="none", const G4String &yunitName="none", const G4String &xfcnName="none", const G4String &yfcnName="none", const G4String &xbinSchemeName="linear")