36 #include "Analysis.hh"
51 fincident_map =
new TH2D(
"incident map",
"Incident Distributuon",
58 fincident_x_hist =
new TH1D(
"incident x",
"Incident X", 100, -5., 5.);
62 fdose_map =
new TH2D(
"dose map",
"Dose Distribution",
69 fdose_hist =
new TH1D(
"dose",
"Dose Distribution", 500, 0., 50.);
80 delete fincident_x_hist;
94 fincident_map->
Reset();
95 fincident_x_hist->
Reset();
105 TFile*
file =
new TFile(fname.c_str(),
"RECREATE",
"Geant4 ROOT analysis");
107 fincident_map-> Write();
108 fincident_x_hist-> Write();
110 fdose_hist-> Write();
122 if ( ! fincidentFlag ) {
124 fincident_x_hist->
Fill(p.
x()/
cm);
126 fincidentFlag =
true;
136 if ( std::abs(p.
y()) < dxy ) {
139 if ( std::abs(p.
x()) < dxy ) {
h1 GetXaxis() -> SetLabelSize(0.025)
hist SetTitle("Cytoplasm intensity on transverse section")
ntupleExperimental Fill(depthExp, EdepExp)
void FillDose(const G4ThreeVector &p, G4double dedx)
h1 GetYaxis() -> SetLabelSize(0.025)
static Analysis * GetAnalysis()
void Save(const G4String &fname)
ntupleExperimental Reset()
void FillIncident(const G4ThreeVector &p)