69 sprintf(a,
"%d", seed);
224 std::cout<<
"n. of events collected in the whole ROG phantom for all geometries: "<<
nTotalEvents<<
G4endl;
225 std::cout<<
"n. of events collected in the whole ROG phantom for the current geometry: "<<
nSingleTotalEvents<<
G4endl;
234 out.open(Filename, std::ios::out);
235 out <<
"Sensitive Detector-Voxels. Total number of events, [mm]->centreX centreY centreZ HalfSizeX HalfSizeY HalfSizeZ minX maxX, minY maxY, minZ maxZ, Dx, Dy, Dz, nX, nY, nZ: \n";
251 out <<
"Phys Volume x [mm], y [mm], z [mm], ix, iy, iz, Dose [Gy], Dose2 [Gy^2], nEvents" <<
G4endl;
259 if (voxels[ix][iy][iz].nEvents>0)
261 out << voxels[ix][iy][iz].
volumeId <<
'\t';
262 out << voxels[ix][iy][iz].
pos.
getX()/
mm <<
'\t';
263 out << voxels[ix][iy][iz].
pos.
getY()/
mm <<
'\t';
264 out << voxels[ix][iy][iz].
pos.
getZ()/
mm <<
'\t';
284 static unsigned int indGeom=0;
286 sprintf(cT,
"%d",indGeom);
void set(double x, double y, double z)
G4bool ProcessHits(G4Step *aStep, G4TouchableHistory *ROHist)
G4int getIdFromVolumeName(G4String name)
std::vector< ExP01TrackerHit * > a
G4LogicalVolume * GetLogicalVolume() const
G4int GetReplicaNumber(G4int depth=0) const
G4StepPoint * GetPreStepPoint() const
G4String fullOutFileDataSingle
static constexpr double mm
G4VPhysicalVolume * GetPhysicalVolume() const
void setFullOutFileDataSingle(G4String val)
static constexpr double kg
G4Material * GetMaterial() const
G4double halfXVoxelDimensionY
const G4String & GetName() const
G4int NumberOfVoxelsAlongY
G4double halfXVoxelDimensionX
G4int saving_in_ROG_Voxels_every_events
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
G4double GetTotalEnergyDeposit() const
G4double halfXVoxelDimensionZ
G4double depEnergyNormError
void saveData(G4String Filename, Svoxel ***voxels)
std::vector< SvolumeNameId > volumeNameIdLink
CML2SDWithVoxels(G4String name, G4int saving_in_ROG_Voxels_every_events, G4int seed, G4String ROGOutFile, G4bool bSaveROG, G4ThreeVector centre, G4ThreeVector halfSize, G4int NumberOfVoxelsAlongX, G4int NumberOfVoxelsAlongY, G4int NumberOfVoxelsAlongZ)
static constexpr double joule
G4int NumberOfVoxelsAlongZ
G4int NumberOfVoxelsAlongX
G4double GetDensity() const