OutputMagneticFieldDDToDDL Class Reference

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 ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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 34 of file OutputMagneticFieldDDToDDL.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 72 of file OutputMagneticFieldDDToDDL.cc.

References gather_cfg::cout, alignBH_cfg::fixed, 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 93 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 243 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 262 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 281 of file OutputMagneticFieldDDToDDL.cc.

References DDLogicalPart::attachedSpecifics(), and DEFINE_FWK_MODULE.

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 42 of file OutputMagneticFieldDDToDDL.cc.

{}

void OutputMagneticFieldDDToDDL::beginJob ( void  )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 40 of file OutputMagneticFieldDDToDDL.cc.

References iEvent.

{}

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

override

Definition at line 101 of file OutputMagneticFieldDDToDDL.cc.

References addToMatStore(), addToSolStore(), addToSpecStore(), alignBH_cfg::fixed, edm::EventSetup::get(), DDCompactView::graph(), mps_fire::i, DDCoreToDDXMLOutput::logicalPart(), m_fname, m_rotNumSeed, m_xos, DDCoreToDDXMLOutput::material(), DDLogicalPart::material(), DDCoreToDDXMLOutput::ns_, MillePedeFileConverter_cfg::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";
}

void OutputMagneticFieldDDToDDL::endJob ( void  )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 44 of file OutputMagneticFieldDDToDDL.cc.

References genParticles_cff::map.

{}

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

inlineoverride

Definition at line 43 of file OutputMagneticFieldDDToDDL.cc.

{}

Member Data Documentation

std::string OutputMagneticFieldDDToDDL::m_fname

private

Definition at line 52 of file OutputMagneticFieldDDToDDL.cc.

Referenced by beginRun(), and OutputMagneticFieldDDToDDL().

int OutputMagneticFieldDDToDDL::m_rotNumSeed

private

Definition at line 51 of file OutputMagneticFieldDDToDDL.cc.

Referenced by beginRun(), and OutputMagneticFieldDDToDDL().

std::ostream* OutputMagneticFieldDDToDDL::m_xos

private

Definition at line 53 of file OutputMagneticFieldDDToDDL.cc.

Referenced by beginRun(), OutputMagneticFieldDDToDDL(), and ~OutputMagneticFieldDDToDDL().