DDMaterial.cc

Go to the documentation of this file.

#include "DetectorDescription/Base/interface/DDdebug.h"
#include "CLHEP/Units/GlobalSystemOfUnits.h"
#include "DetectorDescription/Core/interface/DDMaterial.h"
#include "DetectorDescription/Core/src/Material.h"

using DDI::Material;


//DDBase<DDName,Material*>::StoreT::pointer_type 
//  DDBase<DDName,Material*>::StoreT::instance_ = 0;


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);
  DCOUT('c', std::string(" -> ADDED MATERIAL=") + m.toString() );
  return rep().noOfConstituents();
}


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


DDMaterial::FractionV::value_type DDMaterial::constituent(int i) const 
{ 
  //const DDMaterialImpl::Fraction & f(rep().constituent(i));
  //return std::make_pair(DDMaterial(f.first,true), f.second);  
  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(); 
}

// void DDMaterial::clear()
// {
//   StoreT::instance().clear();
// }
// private ctor for internal use
//DDMaterial::DDMaterial(DDRedirect<DDMaterialImpl>* p, bool dummy)
// : DDRegistered<DDMaterialImpl>(p,true)
//{ 
//  DCOUT_V('C',"DC: (redir) init=" << rep_ ); 
//}   

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);
}