135 const double tol(0.01);
136 for (
unsigned int i = 0;
i <
layers.size();
i++) {
138 double routF =
rMax(zi);
139 int laymax = laymin +
layers[
i];
142 for (
int ly = laymin; ly < laymax; ++ly) {
152 <<
" Front " << zi <<
", " << routF
153 <<
" Back " << zo <<
", " << rinB
154 <<
" superlayer thickness " << layerThick[
i];
162 double rmax = routF *
cos(alpha) - tol;
163 std::vector<double> pgonZ, pgonRin, pgonRout;
164 pgonZ.emplace_back(-0.5 *
thick[ii]);
165 pgonZ.emplace_back(0.5 *
thick[ii]);
166 pgonRin.emplace_back(rinB);
167 pgonRin.emplace_back(rinB);
168 pgonRout.emplace_back(rmax);
169 pgonRout.emplace_back(rmax);
173 pgonZ, pgonRin, pgonRout);
178 <<
" sectors covering " 181 <<
" with " << pgonZ.size() <<
" sections";
182 for (
unsigned int k = 0;
k < pgonZ.size(); ++
k)
184 <<
" R " << pgonRin[
k]
185 <<
":" << pgonRout[
k];
194 <<
" Tubs made of " << matName
195 <<
" of dimensions " << rinB <<
", " 196 << routF <<
", " << 0.5 *
thick[
ii]
199 << glog.
name() <<
" number " <<
copy;
209 <<
" number " << copy <<
" positioned in " 210 << module.
name() <<
" at " <<
r1 217 if (fabs(thickTot - layerThick[
i]) < 0.00001) {
218 }
else if (thickTot > layerThick[i]) {
220 <<
"Thickness of the partition " << layerThick[
i]
221 <<
" is smaller than thickness " << thickTot
222 <<
" of all its components **** ERROR ****\n";
223 }
else if (thickTot < layerThick[i]) {
225 <<
"Thickness of the partition " << layerThick[
i]
226 <<
" does not match with " << thickTot <<
" of the components\n";
DDMaterial is used to define and access material information.
std::vector< int > copyNumber
DDName is used to identify DDD entities uniquely.
constexpr NumType convertRadToDeg(NumType radians)
std::vector< int > layerSense
constexpr long double piRadians(M_PI)
A DDSolid represents the shape of a part.
std::vector< int > layers
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Represents a uniquely identifyable rotation matrix.
U second(std::pair< T, U > const &p)
std::vector< double > slopeB
std::vector< std::string > names
std::vector< std::string > materials
Cos< T >::type cos(const T &t)
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
std::vector< int > layerType
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string ©no, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=0)
std::vector< double > thick
std::vector< double > layerThick
void positionSensitive(DDLogicalPart &glog, double rin, double rout, DDCompactView &cpv)
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = ':'
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)