#include <OutputMagneticFieldDDToDDL.h>

Inheritance diagram for OutputMagneticFieldDDToDDL:

Public Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &) override

void	beginJob () override

void	beginRun (edm::Run const &iEvent, edm::EventSetup const &) override

void	endJob () override

void	endRun (edm::Run const &iEvent, edm::EventSetup const &) override

	OutputMagneticFieldDDToDDL (const edm::ParameterSet &iConfig)

	~OutputMagneticFieldDDToDDL (void)

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	addToMatStore (const DDMaterial &mat, std::set< DDMaterial > &matStore)

void	addToSolStore (const DDSolid &sol, std::set< DDSolid > &solStore, std::set< DDRotation > &rotStore)

void	addToSpecStore (const DDLogicalPart &lp, std::map< const DDsvalues_type, std::set< const DDPartSelection * >, ddsvaluesCmp > &specStore)

Private Attributes
std::string	m_fname

int	m_rotNumSeed

std::ostream *	m_xos

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 23 of file OutputMagneticFieldDDToDDL.h.

Constructor & Destructor Documentation

OutputMagneticFieldDDToDDL::OutputMagneticFieldDDToDDL ( const edm::ParameterSet & iConfig )

explicit

Definition at line 33 of file OutputMagneticFieldDDToDDL.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_fname, m_rotNumSeed, m_xos, and AlCaHLTBitMon_QueryRunRegistry::string.

   : m_fname()
 {
   m_rotNumSeed = iConfig.getParameter<int>( "rotNumSeed" );
   m_fname = iConfig.getUntrackedParameter<std::string>( "fileName" );
   if( m_fname == "" )
   {
     m_xos = &std::cout;
   }
   else
   {
     m_xos = new std::ofstream( m_fname.c_str());
   }
   
   ( *m_xos ) << "<?xml version=\"1.0\"?>\n";
   ( *m_xos ) << "<DDDefinition xmlns=\"http://www.cern.ch/cms/DDL\"\n";
   ( *m_xos ) << " xmlns:xsi=\"http://www.w3.org/2001/XMLSchema-instance\"\n";
   ( *m_xos ) << "xsi:schemaLocation=\"http://www.cern.ch/cms/DDL ../../../DetectorDescription/Schema/DDLSchema.xsd\">\n";
   ( *m_xos ) << std::fixed << std::setprecision( 18 );
 }

OutputMagneticFieldDDToDDL::~OutputMagneticFieldDDToDDL ( void )

Definition at line 54 of file OutputMagneticFieldDDToDDL.cc.

References m_xos.

 {
   ( *m_xos ) << "</DDDefinition>\n";
   ( *m_xos ) << std::endl;
   m_xos->flush();
 }

Member Function Documentation

void OutputMagneticFieldDDToDDL::addToMatStore	(	const DDMaterial &	mat,
		std::set< DDMaterial > &	matStore
	)

private

Definition at line 204 of file OutputMagneticFieldDDToDDL.cc.

References DDMaterial::constituent(), cropTnPTrees::frac, and DDMaterial::noOfConstituents().

Referenced by beginRun().

 {
   matStore.insert( mat );
   if( mat.noOfConstituents() != 0 )
   {
     DDMaterial::FractionV::value_type frac;
     int findex( 0 );
     while( findex < mat.noOfConstituents())
     {
       if( matStore.find( mat.constituent( findex ).first ) == matStore.end())
       {
     addToMatStore( mat.constituent( findex ).first, matStore );
       }
       ++findex;
     }
   }
 }

void OutputMagneticFieldDDToDDL::addToSolStore	(	const DDSolid &	sol,
		std::set< DDSolid > &	solStore,
		std::set< DDRotation > &	rotStore
	)

private

Definition at line 223 of file OutputMagneticFieldDDToDDL.cc.

References ddintersection, ddsubtraction, ddunion, DDBooleanSolid::rotation(), DDSolid::shape(), DDBooleanSolid::solidA(), and DDBooleanSolid::solidB().

Referenced by beginRun().

 {
   solStore.insert( sol );
   if( sol.shape() == ddunion || sol.shape() == ddsubtraction || sol.shape() == ddintersection )
   {
     const DDBooleanSolid& bs ( sol );
     if( solStore.find(bs.solidA()) == solStore.end())
     {
       addToSolStore( bs.solidA(), solStore, rotStore );
     }
     if( solStore.find( bs.solidB()) == solStore.end())
     {
       addToSolStore( bs.solidB(), solStore, rotStore );
     }
     rotStore.insert( bs.rotation());
   }
 }

void OutputMagneticFieldDDToDDL::addToSpecStore	(	const DDLogicalPart &	lp,
		std::map< const DDsvalues_type, std::set< const DDPartSelection * >, ddsvaluesCmp > &	specStore
	)

private

Definition at line 242 of file OutputMagneticFieldDDToDDL.cc.

References DDLogicalPart::attachedSpecifics().

Referenced by beginRun().

 {
   std::vector<std::pair<const DDPartSelection*, const DDsvalues_type*> >::const_iterator spit( lp.attachedSpecifics().begin()), spend( lp.attachedSpecifics().end());
   for( ; spit != spend; ++spit )
   {
     specStore[ *spit->second ].insert( spit->first );
   }
 }

void OutputMagneticFieldDDToDDL::analyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &
	)

inlineoverridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 31 of file OutputMagneticFieldDDToDDL.h.

31 {}

void OutputMagneticFieldDDToDDL::beginJob ( void )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 29 of file OutputMagneticFieldDDToDDL.h.

29 {}

void OutputMagneticFieldDDToDDL::beginRun	(	edm::Run const &	iEvent,
		edm::EventSetup const &	es
	)

override

Definition at line 62 of file OutputMagneticFieldDDToDDL.cc.

References addToMatStore(), addToSolStore(), addToSpecStore(), edm::EventSetup::get(), i, DDCoreToDDXMLOutput::logicalPart(), m_fname, m_rotNumSeed, m_xos, DDCoreToDDXMLOutput::material(), DDLogicalPart::material(), DDCoreToDDXMLOutput::ns_, dbtoconf::out, DDCoreToDDXMLOutput::position(), alignCSCRings::r, cmsRelvalreport::red(), DDCoreToDDXMLOutput::rotation(), DDCoreToDDXMLOutput::solid(), DDLogicalPart::solid(), DDCoreToDDXMLOutput::specpar(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   edm::LogInfo( "OutputMagneticFieldDDToDDL" ) << "OutputMagneticFieldDDToDDL::beginRun";
 
   edm::ESTransientHandle<DDCompactView> pDD;
   es.get<IdealMagneticFieldRecord>().get( pDD );
 
   DDCompactView::DDCompactView::graph_type gra = pDD->graph();
   
   // Temporary stores:
   std::set<DDLogicalPart> lpStore;
   std::set<DDMaterial> matStore;
   std::set<DDSolid> solStore;
   
   std::map<const DDsvalues_type, std::set<const DDPartSelection*>, ddsvaluesCmp > specStore;
   std::set<DDRotation> rotStore;
 
   DDCoreToDDXMLOutput out;
   
   std::string rn = m_fname;
   size_t foundLastDot= rn.find_last_of('.');
   size_t foundLastSlash= rn.find_last_of('/');
 
   if( foundLastSlash > foundLastDot && foundLastSlash != std::string::npos )
   {
     edm::LogError( "OutputMagneticFieldDDToDDL" ) << "What? last . before last / in path for filename... this should die...";
   }
   if( foundLastDot != std::string::npos && foundLastSlash != std::string::npos )
   {
     out.ns_ = rn.substr( foundLastSlash, foundLastDot );
   }
   else if ( foundLastDot != std::string::npos )
   {
     out.ns_ = rn.substr(0, foundLastDot);
   }
   else
   {
     edm::LogError( "OutputMagneticFieldDDToDDL" ) << "What? no file name? Attempt at namespace =\"" << out.ns_ << "\" filename was " << m_fname;
   }
   
   edm::LogInfo( "OutputMagneticFieldDDToDDL" ) << "m_fname=" << m_fname << " namespace = " << out.ns_;
   std::string ns_ = out.ns_;
 
   ( *m_xos ) << std::fixed << std::setprecision( 18 );
   
   typedef  DDCompactView::graph_type::const_adj_iterator adjl_iterator;
 
   adjl_iterator git = gra.begin();
   adjl_iterator gend = gra.end();    
     
   DDCompactView::graph_type::index_type i=0;
   ( *m_xos) << "<PosPartSection label=\"" << ns_ << "\">\n";
   git = gra.begin();
   for( ; git != gend; ++git ) 
   {
     const DDLogicalPart & ddLP = gra.nodeData( git );
     if( lpStore.find(ddLP) != lpStore.end())
     {
       addToSpecStore( ddLP, specStore );
     }
     lpStore.insert( ddLP );
     addToMatStore( ddLP.material(), matStore );
     addToSolStore( ddLP.solid(), solStore, rotStore );
     ++i;
     if( git->size()) 
     {
       // ask for children of ddLP  
       DDCompactView::graph_type::edge_list::const_iterator cit  = git->begin();
       DDCompactView::graph_type::edge_list::const_iterator cend = git->end();
       for( ; cit != cend; ++cit ) 
       {
     const DDLogicalPart & ddcurLP = gra.nodeData( cit->first );
     if( lpStore.find(ddcurLP) != lpStore.end())
     {
       addToSpecStore( ddcurLP, specStore );
     }
     lpStore.insert( ddcurLP );
     addToMatStore( ddcurLP.material(), matStore );
     addToSolStore( ddcurLP.solid(), solStore, rotStore );
     rotStore.insert( gra.edgeData( cit->second )->rot_ );
     out.position( ddLP, ddcurLP, gra.edgeData( cit->second ), m_rotNumSeed, *m_xos );
       } // iterate over children
     } // if (children)
   } // iterate over graph nodes  
 
   ( *m_xos ) << "</PosPartSection>\n";
   
   ( *m_xos ) << std::scientific << std::setprecision( 18 );
   std::set<DDMaterial>::const_iterator it( matStore.begin()), ed( matStore.end());
   ( *m_xos) << "<MaterialSection label=\"" << ns_ << "\">\n";
   for( ; it != ed; ++it )
   {
     if( ! it->isDefined().second ) continue;
     out.material( *it, *m_xos );
   }
   ( *m_xos ) << "</MaterialSection>\n";
 
   ( *m_xos ) << "<RotationSection label=\"" << ns_ << "\">\n";
   ( *m_xos ) << std::fixed << std::setprecision( 18 );
   std::set<DDRotation>::iterator rit( rotStore.begin()), red( rotStore.end());
   for( ; rit != red; ++rit )
   {
     if( ! rit->isDefined().second ) continue;
     if( rit->toString() != ":" )
     {
       DDRotation r( *rit );
       out.rotation( r, *m_xos );
     }
   }
   ( *m_xos ) << "</RotationSection>\n";
 
   ( *m_xos ) << std::fixed << std::setprecision( 18 );
   std::set<DDSolid>::const_iterator sit( solStore.begin()), sed( solStore.end());
   ( *m_xos ) << "<SolidSection label=\"" << ns_ << "\">\n";
   for( ; sit != sed; ++sit)
   {
     if( ! sit->isDefined().second ) continue;  
     out.solid( *sit, *m_xos );
   }
   ( *m_xos ) << "</SolidSection>\n";
 
   std::set<DDLogicalPart>::iterator lpit( lpStore.begin()), lped( lpStore.end());
   ( *m_xos ) << "<LogicalPartSection label=\"" << ns_ << "\">\n";
   for( ; lpit != lped; ++lpit )
   {
     if( ! lpit->isDefined().first ) continue;  
     const DDLogicalPart & lp = *lpit;
     out.logicalPart( lp, *m_xos );
   }
   ( *m_xos ) << "</LogicalPartSection>\n";
 
   ( *m_xos ) << std::fixed << std::setprecision( 18 );
   std::map<DDsvalues_type, std::set<const DDPartSelection*> >::const_iterator mit( specStore.begin()), mend( specStore.end());
   ( *m_xos ) << "<SpecParSection label=\"" << ns_ << "\">\n";
   for( ; mit != mend; ++mit )
   {
     out.specpar( *mit, *m_xos );
   }
   ( *m_xos ) << "</SpecParSection>\n";
 }