DDLCompositeMaterial.cc

Go to the documentation of this file.

#include "DetectorDescription/Parser/src/DDLCompositeMaterial.h"

#include <stddef.h>
#include <map>
#include <utility>
#include "DetectorDescription/Core/interface/DDMaterial.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/interface/ClhepEvaluator.h"
#include "DetectorDescription/Parser/interface/DDLElementRegistry.h"
#include "DetectorDescription/Parser/src/DDLMaterial.h"
#include "DetectorDescription/Parser/src/DDXMLElement.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

class DDCompactView;

DDLCompositeMaterial::DDLCompositeMaterial( DDLElementRegistry* myreg )
  : DDLMaterial( myreg )
{}

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

void
DDLCompositeMaterial::processElement( const std::string& name, const std::string& nmspace, DDCompactView& cpv )
{
  ClhepEvaluator & ev = myRegistry_->evaluator();
  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.
  auto myMF = myRegistry_->getElement( "MaterialFraction" );
  auto myrMaterial = myRegistry_->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, cpv );
  myMF->clear();
  clear();
}