DDLPgonGenerator.cc

Go to the documentation of this file.

#include "DetectorDescription/Parser/src/DDLPgonGenerator.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "DetectorDescription/Core/interface/ClhepEvaluator.h"
#include "DetectorDescription/Parser/interface/DDLElementRegistry.h"
#include "DetectorDescription/Parser/src/DDLSolid.h"
#include "DetectorDescription/Parser/src/DDXMLElement.h"

#include <cstddef>
#include <map>
#include <utility>
#include <vector>

class DDCompactView;

DDLPgonGenerator::DDLPgonGenerator( DDLElementRegistry* myreg )
  : DDLSolid( myreg )
{}

void
DDLPgonGenerator::preProcessElement( const std::string& name, const std::string& nmspace, DDCompactView& cpv )
{
  myRegistry_->getElement("XYPoint")->clear();
  myRegistry_->getElement("ZXYSection")->clear();
}

// Upon encountering an end Extruded Polygone tag, process the XYPoints
// element and extract the x and y std::vectors to feed into DDCore. Then, clear
// the XYPoints and clear this element.
void
DDLPgonGenerator::processElement( const std::string& name,
                  const std::string& nmspace,
                  DDCompactView& cpv )
{
  auto myXYPoints = myRegistry_->getElement( "XYPoint" );
  auto myZXYSection = myRegistry_->getElement( "ZXYSection" );

  ClhepEvaluator & ev = myRegistry_->evaluator();
  DDXMLAttribute atts;

  // get x and y
  std::vector<double> x, y;
  for( size_t i = 0; i < myXYPoints->size(); ++i )
  {
    atts = myXYPoints->getAttributeSet( i );
    auto xit = atts.find( "x" );
    if( xit != atts.end())
      x.emplace_back( ev.eval( nmspace, xit->second ));
    auto yit = atts.find( "y" );
    if( yit != atts.end())
      y.emplace_back( ev.eval( nmspace, yit->second ));
  }
  assert( x.size() == y.size());

  // get zSection information
  std::vector<double> z, zx, zy, zscale;

  for( size_t i = 0; i < myZXYSection->size(); ++i )
  {
    atts = myZXYSection->getAttributeSet( i );
    auto zit = atts.find( "z" );
    if( zit != atts.end())
      z.emplace_back( ev.eval( nmspace, zit->second ));
    auto xit = atts.find( "x" );
    if( xit != atts.end())
      zx.emplace_back( ev.eval( nmspace, xit->second ));
    auto yit = atts.find( "y" );
    if( yit != atts.end())
      zy.emplace_back( ev.eval( nmspace, yit->second ));
    auto sit = atts.find( "scale" );
    if( sit != atts.end())
      zscale.emplace_back( std::stod( sit->second ));
  }
  assert( z.size() == zx.size());
  assert( z.size() == zy.size());
  assert( z.size() == zscale.size());

  atts = getAttributeSet();
  if( name == "ExtrudedPolygon" )
  {
    DDSolid extrudedpolygon = 
      DDSolidFactory::extrudedpolygon( getDDName( nmspace ),
                       x, y,
                       z, zx, zy, zscale );
  } else
  {
    std::string msg = "\nDDLPgonGenerator::processElement was called with incorrect name of solid: " + name;
    throwError( msg );
  }
  DDLSolid::setReference( nmspace, cpv );

  // clear out XYPoint element accumulator and ZXYSections
  myXYPoints->clear();
  myZXYSection->clear();
  clear();
}