---

DDLCompositeMaterial.cc

Go to the documentation of this file.

/***************************************************************************
                          DDLCompositeMaterial.cc  -  description
                             -------------------
    begin                : Wed Oct 31 2001
    email                : case@ucdhep.ucdavis.edu
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *           DDDParser sub-component of DDD                                *
 *                                                                         *
 ***************************************************************************/



// -------------------------------------------------------------------------
// Includes
// -------------------------------------------------------------------------
#include "DDLCompositeMaterial.h"
#include "DDLElementRegistry.h"
#include "DDXMLElement.h"

// DDCore dependencies
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Base/interface/DDdebug.h"
#include "DetectorDescription/Core/interface/DDMaterial.h"
#include "DetectorDescription/Base/interface/DDException.h"

#include "DetectorDescription/ExprAlgo/interface/ExprEvalSingleton.h"

//#include <strstream>
#include <string>

// Default constructor.
DDLCompositeMaterial::DDLCompositeMaterial()
{
}

// Default destructor.
DDLCompositeMaterial::~DDLCompositeMaterial()
{
}

// to initialize the CompositeMaterial, clear all rMaterials in case some other 
// rMaterial was used for some other element.
void DDLCompositeMaterial::preProcessElement (const std::string& type, const std::string& nmspace)
{
  // fyi: no need to clear MaterialFraction because it is cleared at the end of each
  // CompositeMaterial
  DDLElementRegistry::getElement("rMaterial")->clear();
}

void DDLCompositeMaterial::processElement (const std::string& type, const std::string& nmspace)
{
  DCOUT_V('P', "DDLCompositeMaterial::processElement started");


  ExprEvalInterface & ev = ExprEvalSingleton::instance();
  DDXMLAttribute atts = getAttributeSet();

  DDName ddn = getDDName(nmspace);
  DDMaterial mat;

  mat = DDMaterial(ddn, ev.eval(nmspace, atts.find("density")->second));
  
  // Get references to relevant DDL elements that are needed.
  DDXMLElement* myMF = DDLElementRegistry::getElement("MaterialFraction");
  DDXMLElement* myrMaterial = DDLElementRegistry::getElement("rMaterial");

  // Get the names from those elements and also the namespace for the reference element.
  // The parent element CompositeMaterial MUST be in the same namespace as this fraction.
  // additionally, because it is NOT a reference, we do not try to dis-entangle the namespace.
  // That is, we do not use the getName() which searches the name for a colon, but instead use
  // the "raw" name attribute.

  // TO DO:  sfractions assumes that the values are "mixture by weight" (I think)
  // we need to retrieve the fraction attributes and then check the method
  // attribute of the CompositeMaterial to determine if the numbers go straight
  // in to the DDCore or if the numbers need to be manipulated so that the 
  // sfractions are really mixtures by weight going into the DDCore.  Otherwise,
  // the DDCore has to know which are which and right now it does not.

  if (myMF->size() != myrMaterial->size())
    {
      std::string msg = "/nDDLCompositeMaterial::processElement found that the ";
      msg += "number of MaterialFractions does not match the number ";
      msg += "of rMaterial names for ";
      msg += ddn.ns() + ":" + ddn.name() + ".";
      throwError(msg);
    }
  for (size_t i = 0; i < myrMaterial->size(); ++i)
    {
      atts = myMF->getAttributeSet(i);
      mat.addMaterial(myrMaterial->getDDName(nmspace, "name", i)
                      , ev.eval(nmspace, atts.find("fraction")->second));
    }

  // clears and sets new reference to THIS material.
  DDLMaterial::setReference(nmspace);
  myMF->clear();
  clear();
  // print it.
  DCOUT_V('P', "DDLCompositeMaterial::processElement completed " << mat);
}

---