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);
146 (*m_xos) <<
"<PosPartSection label=\"" << ns_ <<
"\">\n";
148 for (; git != gend; ++git) {
150 if (lpStore.find(ddLP) != lpStore.end()) {
153 lpStore.insert(ddLP);
159 auto cit = git->begin();
160 auto cend = git->end();
161 for (; cit != cend; ++cit) {
163 if (lpStore.find(ddcurLP) != lpStore.end()) {
166 lpStore.insert(ddcurLP);
169 rotStore.insert(gra.edgeData(cit->second)->ddrot());
175 (*m_xos) <<
"</PosPartSection>\n";
177 (*m_xos) << std::scientific << std::setprecision(18);
178 std::set<DDMaterial>::const_iterator it(matStore.begin()), ed(matStore.end());
179 (*m_xos) <<
"<MaterialSection label=\"" << ns_ <<
"\">\n";
180 for (; it != ed; ++it) {
181 if (!it->isDefined().second)
185 (*m_xos) <<
"</MaterialSection>\n";
187 (*m_xos) <<
"<RotationSection label=\"" << ns_ <<
"\">\n";
188 (*m_xos) <<
std::fixed << std::setprecision(18);
189 std::set<DDRotation>::iterator rit(rotStore.begin()),
red(rotStore.end());
190 for (; rit !=
red; ++rit) {
191 if (!rit->isDefined().second)
193 if (rit->toString() !=
":") {
198 (*m_xos) <<
"</RotationSection>\n";
200 (*m_xos) <<
std::fixed << std::setprecision(18);
201 std::set<DDSolid>::const_iterator sit(solStore.begin()), sed(solStore.end());
202 (*m_xos) <<
"<SolidSection label=\"" << ns_ <<
"\">\n";
203 for (; sit != sed; ++sit) {
204 if (!sit->isDefined().second)
208 (*m_xos) <<
"</SolidSection>\n";
210 std::set<DDLogicalPart>::iterator lpit(lpStore.begin()), lped(lpStore.end());
211 (*m_xos) <<
"<LogicalPartSection label=\"" << ns_ <<
"\">\n";
212 for (; lpit != lped; ++lpit) {
213 if (!lpit->isDefined().first)
218 (*m_xos) <<
"</LogicalPartSection>\n";
220 (*m_xos) <<
std::fixed << std::setprecision(18);
221 std::map<DDsvalues_type, std::set<const DDPartSelection*> >::const_iterator mit(specStore.begin()),
222 mend(specStore.end());
223 (*m_xos) <<
"<SpecParSection label=\"" << ns_ <<
"\">\n";
224 for (; mit != mend; ++mit) {
227 (*m_xos) <<
"</SpecParSection>\n";
math::Graph< DDLogicalPart, DDPosData * > Graph
edm::ESGetToken< DDCompactView, IdealMagneticFieldRecord > ddToken_
std::vector< double >::size_type index_type
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.