57 G4double maxEnergy = 100;
60 auto analysisManager = G4AnalysisManager::Instance();
63 analysisManager->CreateH1(
"energyParticle",
64 "Primary energy;E_{MC} (GeV);Entries", 256, 0,
66 analysisManager->CreateH1(
"energyDeposited",
67 "Deposited energy;E_{MC} (GeV);Entries", 256, 0,
69 analysisManager->CreateH1(
71 "Ratio of energy deposited to primary;E_{dep} / E_{MC};Entries", 1024, 0,
73 analysisManager->CreateH1(
"time",
"Simulation time; time (s);Entries", 2048,
75 analysisManager->CreateH1(
76 "longProfile",
"Longitudinal profile;t (mm);#LTE#GT (MeV)", cellNumZ,
77 -0.5 * cellSizeZ, (cellNumZ - 0.5) * cellSizeZ);
78 analysisManager->CreateH1(
79 "transProfile",
"Transverse profile;r (mm);#LTE#GT (MeV)", cellNumRho,
80 -0.5 * cellSizeRho, (cellNumRho - 0.5) * cellSizeRho);
81 analysisManager->CreateH1(
83 "First moment of longitudinal distribution;#LT#lambda#GT (mm);Entries",
84 1024, -0.5 * cellSizeZ,
85 cellNumZ * cellSizeZ / 2);
86 analysisManager->CreateH1(
"transFirstMoment",
87 "First moment of transverse distribution;#LTr#GT "
89 1024, -0.5 * cellSizeRho,
90 cellNumRho * cellSizeRho /
92 analysisManager->CreateH1(
94 "Second moment of longitudinal distribution;#LT#lambda^{2}#GT "
97 std::pow(cellNumZ * cellSizeZ, 2) /
99 analysisManager->CreateH1(
101 "Second moment of transverse distribution;#LTr^{2}#GT (mm^{2});Entries",
103 std::pow(cellNumRho * cellSizeRho, 2) /
105 analysisManager->CreateH1(
106 "hitType",
"hit type;type (0=full, 1= fast);Entries", 2, -0.5, 1.5);
109 analysisManager->OpenFile();