169 const double tol(0.01);
170 for (
unsigned int i = 0;
i <
layers_.size();
i++) {
172 double routF =
rMax(zi);
176 for (
int ly = laymin; ly < laymax; ++ly) {
185 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalModule: Layer " << ly <<
":" <<
ii <<
" Front " << zi <<
", " << routF
186 <<
" Back " << zo <<
", " << rinB <<
" superlayer thickness " <<
layerThick_[
i];
193 double rmax = routF *
cos(
alpha) - tol;
194 std::vector<double> pgonZ, pgonRin, pgonRout;
195 pgonZ.emplace_back(-0.5 *
thick_[
ii]);
196 pgonZ.emplace_back(0.5 *
thick_[
ii]);
197 pgonRin.emplace_back(rinB);
198 pgonRin.emplace_back(rinB);
199 pgonRout.emplace_back(rmax);
200 pgonRout.emplace_back(rmax);
208 for (
unsigned int k = 0;
k < pgonZ.size(); ++
k)
209 edm::LogVerbatim(
"HGCalGeom") <<
"[" <<
k <<
"] z " << pgonZ[
k] <<
" R " << pgonRin[
k] <<
":" << pgonRout[
k];
215 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalModule: " << solid.
name() <<
" Tubs made of " << matName
216 <<
" of dimensions " << rinB <<
", " << routF <<
", " << 0.5 *
thick_[
ii]
228 <<
module.name() <<
" at " <<
r1 <<
" with " <<
rot;
237 << thickTot <<
" of all its components **** ERROR ****\n";
240 << thickTot <<
" of the components\n";
Log< level::Info, true > LogVerbatim
std::vector< std::string > names_
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string ©no, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=nullptr)
std::vector< std::string > materials_
DDMaterial is used to define and access material information.
constexpr NumType convertRadToDeg(NumType radians)
DDName is used to identify DDD entities uniquely.
std::string to_string(const V &value)
Log< level::Error, false > LogError
std::vector< int > copyNumber_
A DDSolid represents the shape of a part.
Represents a uniquely identifyable rotation matrix.
U second(std::pair< T, U > const &p)
std::vector< int > layerSense_
std::vector< int > layerType_
std::vector< double > slopeB_
Cos< T >::type cos(const T &t)
std::vector< double > thick_
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
static constexpr double tol_
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
std::vector< int > layers_
void positionSensitive(DDLogicalPart &glog, double rin, double rout, DDCompactView &cpv)
Log< level::Warning, false > LogWarning
std::vector< double > layerThick_
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = ':'
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
static DDSolid polyhedra(const DDName &name, int sides, double startPhi, double deltaPhi, const std::vector< double > &z, const std::vector< double > &rmin, const std::vector< double > &rmax)
Creates a polyhedra (refere to Geant3 or Geant4 documentation)