@@ -76,7 +76,8 @@ void ZdcGeometry::newCell(const GlobalPoint& f1,
7676 const DetId& detId) {
7777 const CaloGenericDetId cgid (detId);
7878#ifdef EDM_ML_DEBUG
79- edm::LogVerbatim (" ZDCGeom" ) << " ZDCGeometry " << HcalZDCDetId (detId) << " Generic ID " << std::hex << cgid.rawId () << std::dec << " ZDC? " << cgid.isZDC ();
79+ edm::LogVerbatim (" ZDCGeom" ) << " ZDCGeometry " << HcalZDCDetId (detId) << " Generic ID " << std::hex << cgid.rawId ()
80+ << std::dec << " ZDC? " << cgid.isZDC ();
8081#endif
8182 assert (cgid.isZDC ());
8283
@@ -95,9 +96,9 @@ const CaloCellGeometry* ZdcGeometry::getGeometryRawPtr(uint32_t index) const {
9596}
9697
9798void ZdcGeometry::getSummary (CaloSubdetectorGeometry::TrVec& tVec,
98- CaloSubdetectorGeometry::IVec& iVec,
99- CaloSubdetectorGeometry::DimVec& dVec,
100- CaloSubdetectorGeometry::IVec& /* dins*/ ) const {
99+ CaloSubdetectorGeometry::IVec& iVec,
100+ CaloSubdetectorGeometry::DimVec& dVec,
101+ CaloSubdetectorGeometry::IVec& /* dins*/ ) const {
101102 tVec.reserve (m_validIds.size () * numberOfTransformParms ());
102103 iVec.reserve (numberOfShapes () == 1 ? 1 : m_validIds.size ());
103104 dVec.reserve (numberOfShapes () * numberOfParametersPerShape ());
@@ -112,51 +113,53 @@ void ZdcGeometry::getSummary(CaloSubdetectorGeometry::TrVec& tVec,
112113 Tr3D tr;
113114 std::shared_ptr<const CaloCellGeometry> ptr (cellGeomPtr (i));
114115#ifdef EDM_ML_DEBUG
115- edm::LogVerbatim (" ZDCGeom" ) << " ZDCGeometry:Summary " << i << " :" << HcalZDCDetId::kSizeForDenseIndexingRun1 << " Pointer " << ptr << " :" << (nullptr != ptr);
116+ edm::LogVerbatim (" ZDCGeom" ) << " ZDCGeometry:Summary " << i << " :" << HcalZDCDetId::kSizeForDenseIndexingRun1
117+ << " Pointer " << ptr << " :" << (nullptr != ptr);
116118#endif
117119 if (i < static_cast <int32_t >(HcalZDCDetId::kSizeForDenseIndexingRun1 ))
118120 assert (nullptr != ptr);
119121 if (ptr != nullptr ) {
120122 ptr->getTransform (tr, (Pt3DVec*)nullptr );
121123
122124 if (Tr3D () == tr) { // for preshower there is no rotation
123- const GlobalPoint& gp (ptr->getPosition ());
124- tr = HepGeom::Translate3D (gp.x (), gp.y (), gp.z ());
125+ const GlobalPoint& gp (ptr->getPosition ());
126+ tr = HepGeom::Translate3D (gp.x (), gp.y (), gp.z ());
125127 }
126128
127129 const CLHEP::Hep3Vector tt (tr.getTranslation ());
128130 tVec.emplace_back (tt.x ());
129131 tVec.emplace_back (tt.y ());
130132 tVec.emplace_back (tt.z ());
131133 if (6 == numberOfTransformParms ()) {
132- const CLHEP::HepRotation rr (tr.getRotation ());
133- const ROOT::Math::Transform3D rtr (rr.xx (), rr.xy (), rr.xz (), tt.x (), rr.yx (), rr.yy (), rr.yz (), tt.y (), rr.zx (), rr.zy (), rr.zz (), tt.z ());
134- ROOT::Math::EulerAngles ea;
135- rtr.GetRotation (ea);
136- tVec.emplace_back (ea.Phi ());
137- tVec.emplace_back (ea.Theta ());
138- tVec.emplace_back (ea.Psi ());
134+ const CLHEP::HepRotation rr (tr.getRotation ());
135+ const ROOT::Math::Transform3D rtr (
136+ rr.xx (), rr.xy (), rr.xz (), tt.x (), rr.yx (), rr.yy (), rr.yz (), tt.y (), rr.zx (), rr.zy (), rr.zz (), tt.z ());
137+ ROOT::Math::EulerAngles ea;
138+ rtr.GetRotation (ea);
139+ tVec.emplace_back (ea.Phi ());
140+ tVec.emplace_back (ea.Theta ());
141+ tVec.emplace_back (ea.Psi ());
139142 }
140143
141144 const CCGFloat* par (ptr->param ());
142145
143146 unsigned int ishape (9999 );
144147 for (unsigned int ivv (0 ); ivv != parVecVec ().size (); ++ivv) {
145- bool ok (true );
146- const CCGFloat* pv (&(*parVecVec ()[ivv].begin ()));
147- for (unsigned int k (0 ); k != numberOfParametersPerShape (); ++k) {
148- ok = ok && (fabs (par[k] - pv[k]) < 1 .e -6 );
149- }
150- if (ok) {
151- ishape = ivv;
152- break ;
153- }
148+ bool ok (true );
149+ const CCGFloat* pv (&(*parVecVec ()[ivv].begin ()));
150+ for (unsigned int k (0 ); k != numberOfParametersPerShape (); ++k) {
151+ ok = ok && (fabs (par[k] - pv[k]) < 1 .e -6 );
152+ }
153+ if (ok) {
154+ ishape = ivv;
155+ break ;
156+ }
154157 }
155158 assert (9999 != ishape);
156159
157160 const unsigned int nn ((numberOfShapes () == 1 ) ? (unsigned int )1 : m_validIds.size ());
158161 if (iVec.size () < nn)
159- iVec.emplace_back (ishape);
162+ iVec.emplace_back (ishape);
160163 }
161164 }
162165}
0 commit comments