DDMaterial.cc

Go to the documentation of this file.

#include <ostream>
#include <string>
#include <vector>

#include "DetectorDescription/Core/interface/DDBase.h"
#include "DetectorDescription/Core/interface/DDMaterial.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/src/Material.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "FWCore/Utilities/interface/Exception.h"

using DDI::Material;
using namespace geant_units::operators;

DDMaterial::DDMaterial() : DDBase< DDName, std::unique_ptr<Material>>() { }

DDMaterial::DDMaterial( const DDName & name )
  : DDBase< DDName, std::unique_ptr<Material>>()
{ 
  create( name );
}

DDMaterial::DDMaterial( const DDName & name, double z, double a, double d )
  : DDBase< DDName, std::unique_ptr<Material>>()
{ 
  create( name, std::make_unique<Material>( z, a, d ));
}

DDMaterial::DDMaterial( const DDName & name, double density )
  : DDBase< DDName, std::unique_ptr<Material>>()
{ 
  create( name, std::make_unique<Material>( 0, 0, density ));
}

int
DDMaterial::addMaterial( const DDMaterial & m, double fm )
{  
  if( m.ddname() == ddname()) {
    throw cms::Exception("DDException") << "DDMaterial::addMaterial(..): name-clash\n        trying to add material " << m << " to itself! ";
  }  
  rep().addMaterial( m, fm );
  return rep().noOfConstituents();
}

int
DDMaterial::noOfConstituents() const
{
   return rep().noOfConstituents();
}

DDMaterial::FractionV::value_type DDMaterial::constituent( int i ) const 
{ 
  return rep().constituent( i );
}

double
DDMaterial::a() const
{
  return rep().a(); 
}

double
DDMaterial::z() const
{
  return rep().z(); 
}

double DDMaterial::density() const
{
  return rep().density(); 
}

namespace {
  std::ostream &doStream(std::ostream & os, const DDMaterial & mat, int level)
  {
    ++level; 
    if (mat) {
      os << '[' << mat.name() <<']' << " z=" << mat.z()
     << " a=" << convertUnitsTo(1._g_per_mole, mat.a()) << "*g/mole"
     << " d=" << convertUnitsTo(1._g_per_cm3 , mat.density()) << "*g/cm3";
      std::string s(2*level,' ');
      for (int i=0; i<mat.noOfConstituents(); ++i) {
         DDMaterial::FractionV::value_type f = mat.constituent(i);
         os << std::endl << s << i+1 << " : fm=" << f.second
                    << " : ";
         doStream(os, f.first, level);
      }
    }
    else
      os << "* material not declared * ";
    --level;
    return os;
  }
}

std::ostream & operator<<(std::ostream & os, const DDMaterial & mat)
{ 
  return doStream( os, mat, 0 );
}