101 edm::LogInfo(
"OutputMagneticFieldDDToDDL") <<
"OutputMagneticFieldDDToDDL::beginRun";
105 const auto& gra = pDD->
graph();
108 std::set<DDLogicalPart> lpStore;
109 std::set<DDMaterial> matStore;
110 std::set<DDSolid> solStore;
112 std::map<const DDsvalues_type, std::set<const DDPartSelection*>, ddsvaluesCmp> specStore;
113 std::set<DDRotation> rotStore;
118 size_t foundLastDot = rn.find_last_of(
'.');
119 size_t foundLastSlash = rn.find_last_of(
'/');
121 if (foundLastSlash > foundLastDot && foundLastSlash != std::string::npos) {
123 <<
"What? last . before last / in path for filename... this should die...";
125 if (foundLastDot != std::string::npos && foundLastSlash != std::string::npos) {
126 out.ns_ = rn.substr(foundLastSlash, foundLastDot);
127 }
else if (foundLastDot != std::string::npos) {
128 out.ns_ = rn.substr(0, foundLastDot);
131 <<
"What? no file name? Attempt at namespace =\"" <<
out.ns_ <<
"\" filename was " <<
m_fname;
137 (*m_xos) <<
std::fixed << std::setprecision(18);
145 (*m_xos) <<
"<PosPartSection label=\"" << ns_ <<
"\">\n";
147 for (; git != gend; ++git) {
149 if (lpStore.find(ddLP) != lpStore.end()) {
152 lpStore.insert(ddLP);
157 auto cit = git->begin();
158 auto cend = git->end();
159 for (; cit != cend; ++cit) {
161 if (lpStore.find(ddcurLP) != lpStore.end()) {
164 lpStore.insert(ddcurLP);
167 rotStore.insert(gra.edgeData(cit->second)->ddrot());
173 (*m_xos) <<
"</PosPartSection>\n";
175 (*m_xos) << std::scientific << std::setprecision(18);
176 std::set<DDMaterial>::const_iterator
it(matStore.begin()), ed(matStore.end());
177 (*m_xos) <<
"<MaterialSection label=\"" << ns_ <<
"\">\n";
178 for (;
it != ed; ++
it) {
179 if (!
it->isDefined().second)
183 (*m_xos) <<
"</MaterialSection>\n";
185 (*m_xos) <<
"<RotationSection label=\"" << ns_ <<
"\">\n";
186 (*m_xos) <<
std::fixed << std::setprecision(18);
187 std::set<DDRotation>::iterator rit(rotStore.begin()),
red(rotStore.end());
188 for (; rit !=
red; ++rit) {
189 if (!rit->isDefined().second)
191 if (rit->toString() !=
":") {
196 (*m_xos) <<
"</RotationSection>\n";
198 (*m_xos) <<
std::fixed << std::setprecision(18);
199 std::set<DDSolid>::const_iterator sit(solStore.begin()), sed(solStore.end());
200 (*m_xos) <<
"<SolidSection label=\"" << ns_ <<
"\">\n";
201 for (; sit != sed; ++sit) {
202 if (!sit->isDefined().second)
206 (*m_xos) <<
"</SolidSection>\n";
208 std::set<DDLogicalPart>::iterator lpit(lpStore.begin()), lped(lpStore.end());
209 (*m_xos) <<
"<LogicalPartSection label=\"" << ns_ <<
"\">\n";
210 for (; lpit != lped; ++lpit) {
211 if (!lpit->isDefined().first)
216 (*m_xos) <<
"</LogicalPartSection>\n";
218 (*m_xos) <<
std::fixed << std::setprecision(18);
219 std::map<DDsvalues_type, std::set<const DDPartSelection*> >::const_iterator mit(specStore.begin()),
220 mend(specStore.end());
221 (*m_xos) <<
"<SpecParSection label=\"" << ns_ <<
"\">\n";
222 for (; mit != mend; ++mit) {
225 (*m_xos) <<
"</SpecParSection>\n";
math::Graph< DDLogicalPart, DDPosData * > Graph
edm::ESGetToken< DDCompactView, IdealMagneticFieldRecord > ddToken_
void addToMatStore(const DDMaterial &mat, std::set< DDMaterial > &matStore)
Log< level::Error, false > LogError
Represents a uniquely identifyable rotation matrix.
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
void addToSolStore(const DDSolid &sol, std::set< DDSolid > &solStore, std::set< DDRotation > &rotStore)
Graph::const_adj_iterator adjl_iterator
void addToSpecStore(const DDLogicalPart &lp, std::map< const DDsvalues_type, std::set< const DDPartSelection *>, ddsvaluesCmp > &specStore)
const DDMaterial & material(void) const
Returns a reference object of the material this LogicalPart is made of.
Log< level::Info, false > LogInfo
const DDSolid & solid(void) const
Returns a reference object of the solid being the shape of this LogicalPart.
adj_list::const_iterator const_adj_iterator
const Graph & graph() const
Provides read-only access to the data structure of the compact-view.