36 if (ebg->
parMgr() ==
nullptr)
44 for (doSubDets = fv.firstChild(); doSubDets; doSubDets = fv.nextSibling()) {
47 if (answer !=
"EcalTBH4BeamHits")
51 const DDSolid& solid = fv.logicalPart().solid();
56 <<
"\n CaloGeometryEcalTBHodoscope::upDate(...): currently only box fiber shapes supported ";
57 edm::LogWarning(
"EcalTBHodoscopeGeometry") <<
"Wrong shape for sensitive volume!" << solid;
66 fv.rotation().GetComponents(
x,
y,
z);
67 CLHEP::Hep3Vector hx(
x.X(),
x.Y(),
x.Z());
68 CLHEP::Hep3Vector hy(
y.X(),
y.Y(),
y.Z());
69 CLHEP::Hep3Vector hz(
z.X(),
z.Y(),
z.Z());
70 CLHEP::HepRotation hrot(hx, hy, hz);
71 CLHEP::Hep3Vector htran(fv.translation().X(), fv.translation().Y(), fv.translation().Z());
73 const HepGeom::Transform3D ht3d(hrot,
76 const HepGeom::Point3D<float>
ctr(ht3d * HepGeom::Point3D<float>(0, 0, 0));
80 std::vector<CCGFloat> vv;
81 vv.reserve(
pv.size() + 1);
82 for (
unsigned int i(0);
i !=
pv.size(); ++
i) {
92 ebg->
newCell(refPoint, refPoint, refPoint, pP, detId);
virtual void newCell(const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)=0
const std::vector< double > & parameters(void) const
Give the parameters of the solid.
void allocatePar(ParVec::size_type n, unsigned int m)
A DDSolid represents the shape of a part.
std::string getDDDString(std::string s, DDFilteredView *fv)
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DD3Vector
CaloCellGeometry::CCGFloat CCGFloat
static const CCGFloat * getParmPtr(const std::vector< CCGFloat > &vd, ParMgr *mgr, ParVecVec &pvv)
static const CCGFloat k_ScaleFromDDDtoGeant
DDSolidShape shape(void) const
The type of the solid.
CaloCellGeometry::CornersMgr * cornersMgr()
unsigned int getDetIdForDDDNode(const DDFilteredView &fv)
void allocateCorners(CaloCellGeometry::CornersVec::size_type n)
Log< level::Warning, false > LogWarning