DDMaterial.cc

Go to the documentation of this file.

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

#include "CLHEP/Units/GlobalSystemOfUnits.h"
#include "CLHEP/Units/SystemOfUnits.h"
#include "DetectorDescription/Core/interface/Store.h"
#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 "FWCore/Utilities/interface/Exception.h"

using DDI::Material;

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


DDMaterial::DDMaterial(const DDName & name) : DDBase<DDName,Material*>()
{ 
  prep_ = StoreT::instance().create(name);
}
 

DDMaterial::DDMaterial(const DDName & name, double z, double a, double d)
 : DDBase<DDName,Material*>()
{ 
  prep_ = StoreT::instance().create(name, new Material(z, a, d));
}
              

DDMaterial::DDMaterial(const DDName & name, double density)
 : DDBase<DDName,Material*>()
{ 
  prep_ = StoreT::instance().create(name, new 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=" << mat.a()/g*mole << "*g/mole"
                       << " d=" << mat.density()/g*cm3 << "*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);
}