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


std::ostream & operator<<(std::ostream & os, const DDMaterial & mat)
{ 
  static int level=0;
  ++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 
                  << " : " << f.first;
    }            
    //--level;
  } 
  else
    os << "* material not declared * ";  
  --level;   
  return os;
}