56 G4Exception(
"DicompFileMgr.:BuildMaterials",
59 (
"Compression factor = " + std::to_string(fCompress)
60 +
" has to be a divisor of Number of voxels X = " + std::to_string(
fNoVoxelsX)
61 +
" and Y " + std::to_string(
fNoVoxelsY)).c_str());
66 for(
int ir = 0; ir <
fNoVoxelsY; ir += fCompress ) {
67 for(
int ic = 0; ic <
fNoVoxelsX; ic += fCompress ) {
69 int isumrMax = std::min(ir+fCompress,
fNoVoxelsY);
70 int isumcMax = std::min(ic+fCompress,
fNoVoxelsX);
71 for(
int isumr = ir; isumr < isumrMax; isumr ++ ) {
72 for(
int isumc = ic; isumc < isumcMax; isumc ++ ) {
76 meanHV /= (isumrMax-ir)*(isumcMax-ic);