Skip to content

Commit

Permalink
Add warning if TGeo materialScan hits TGeoTessellated volumes. See gi…
Browse files Browse the repository at this point in the history
…thub issue AIDASoft#1322
  • Loading branch information
MarkusFrankATcernch committed Sep 17, 2024
1 parent 82db94b commit a331842
Showing 1 changed file with 20 additions and 9 deletions.
29 changes: 20 additions & 9 deletions DDRec/src/MaterialScan.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -22,6 +22,7 @@
#include <DD4hep/DD4hepUnits.h>
#include <DD4hep/Detector.h>
#include <DD4hep/Printout.h>
#include <DD4hep/Shapes.h>

/// C/C++ include files
#include <cstdio>
Expand Down Expand Up @@ -193,6 +194,7 @@ void MaterialScan::print(const Vector3D& p0, const Vector3D& p1, double epsilon)
::printf(" | Layer \\ %-16s [g/mole] [g/cm3] [cm] [cm] [cm] [cm] [cm] [cm] ( cm, cm, cm)\n","");
::printf("%s",line);
MaterialVec materials;
std::set<Solid> tessellated_solids;
for( unsigned i=0, n=placements.size(); i<n; ++i){
TGeoNode* pv = placements[i].first.ptr();
double length = placements[i].second;
Expand Down Expand Up @@ -223,6 +225,11 @@ void MaterialScan::print(const Vector3D& p0, const Vector3D& p1, double epsilon)
mname += " -> ";
mname += next_mat->GetName();
}
Volume vol(placements[i].first->GetVolume());
Solid shape(vol.solid());
if ( shape->IsA() == TGeoTessellated::Class() ) {
tessellated_solids.insert(shape);
}
if ( 0 == i ) {
::printf(fmt, "(start)" , curr_mat->GetName(), curr_mat->GetZ(), curr_mat->GetA(),
curr_mat->GetDensity(), curr_mat->GetRadLen()/dd4hep::cm, curr_mat->GetIntLen()/dd4hep::cm,
Expand All @@ -237,18 +244,13 @@ void MaterialScan::print(const Vector3D& p0, const Vector3D& p1, double epsilon)
p0[0]/dd4hep::cm, p0[1]/dd4hep::cm, p0[2]/dd4hep::cm);
}
else {
double next_dens = next_mat ? next_mat->GetDensity() : 0e0;
double next_rad = next_mat ? next_mat->GetRadLen() : 0e0;
double next_int = next_mat ? next_mat->GetIntLen() : 0e0;
double next_Z = next_mat ? next_mat->GetZ() : 0e0;
double next_A = next_mat ? next_mat->GetA() : 0e0;
::printf(fmt, std::to_string(i+1).c_str(), mname.c_str(), next_Z, next_A,
next_dens, next_rad/dd4hep::cm, next_int/dd4hep::cm,
::printf(fmt, std::to_string(i+1).c_str(), mname.c_str(), next_mat->GetZ(), next_mat->GetA(),
next_mat->GetDensity(), next_mat->GetRadLen()/dd4hep::cm, next_mat->GetIntLen()/dd4hep::cm,
length/dd4hep::cm, path_length/dd4hep::cm, sum_x0, sum_lambda,
end[0]/dd4hep::cm, end[1]/dd4hep::cm, end[2]/dd4hep::cm);
}
}
printf("%s",line);
::printf("%s",line);
const MaterialData& avg = matMgr.createAveragedMaterial(materials);
const char* fmt = avg.radiationLength() >= 1e5 ? fmt2 : fmt1;
::printf(fmt,"","Average Material",avg.Z(),avg.A(),avg.density(),
Expand All @@ -257,7 +259,16 @@ void MaterialScan::print(const Vector3D& p0, const Vector3D& p1, double epsilon)
path_length/avg.radiationLength(),
path_length/avg.interactionLength(),
end[0]/dd4hep::cm, end[1]/dd4hep::cm, end[2]/dd4hep::cm);
printf("%s",line);
::printf("%s",line);

if ( !tessellated_solids.empty() ) {
::printf(" | WARNING: %ld tessellated shape were encountered during the volume traversal:\n",
tessellated_solids.size());
for(auto shape : tessellated_solids )
::printf(" | \t\t %s\n", shape.name());
::printf(" | WARNING: The results of this material scan are unreliable!\n");
::printf("%s",line);
}
}

/// Scan along a line and print the materials traversed
Expand Down

0 comments on commit a331842

Please sign in to comment.