#include <iostream>
using std::cout;
using std::cerr;
using std::cin;
using std::endl;
using std::hex;

#include <valarray>
using std::valarray;

#include <fstream>
using std::fstream;

#include <string>
using std::string;

#include <list>
using std::list;

#include <cmath>

#include "TObject.h"
#include "TVector3.h"
#include "TRandom.h"
#include "Math/Vector3D.h"
using ROOT::Math::XYZVector;

#include "Math/Point3D.h"
using ROOT::Math::XYZPoint;

#include "Math/Transform3D.h"
using ROOT::Math::Transform3D;

#include "Math/RotationX.h"
using ROOT::Math::RotationX;
#include "Math/RotationY.h"
using ROOT::Math::RotationY;
#include "Math/RotationZ.h"
using ROOT::Math::RotationZ;
#include "Math/Rotation3D.h"
using ROOT::Math::Rotation3D;

#include "DrcPhoton.h"
#include "DrcOptReflAbs.h"
#include "DrcSurfAbs.h"
#include "DrcOptMatAbs.h"
#include "DrcOptDev.h"


//----------------------------------------------------------------------
DrcOptDev::DrcOptDev()
{
  m_copyNumber  = 0;
  m_name        = "unnamed DrcOptDev";
  m_photonTrace = false;
  m_verbosity   = 0;
}
//----------------------------------------------------------------------
void DrcOptDev::copy(const DrcOptDev& d)
{
  m_copyNumber        = d.m_copyNumber;
  m_name              = d.m_name;
  m_verbosity         = d.m_verbosity;
  m_photonTrace       = d.m_photonTrace;
  m_photonTraceStream = d.m_photonTraceStream;

  m_listSurf.clear();
  list<DrcSurfAbs*>::const_iterator isurf;
  for(isurf=d.m_listSurf.begin(); isurf != d.m_listSurf.end(); ++isurf) 
    {
      const DrcSurfAbs* tmp = (*isurf);
      DrcSurfAbs* tmp1 = tmp->clone();
      m_listSurf.push_back(tmp1);
    }
}
//----------------------------------------------------------------------
DrcOptDev::DrcOptDev(const DrcOptDev& d)
{
  if (d.m_verbosity>=1) cout<<"  DrcOptDev::DrcOptDev"
			    <<"(const DrcOptDev&) name,copy: "
			    <<d.m_name<<" "<<d.m_copyNumber<<endl;
  copy(d);
}
//----------------------------------------------------------------------
DrcOptDev& DrcOptDev::operator=(const DrcOptDev& d)
{
  if (d.m_verbosity>=1) cout<<"  DrcOptDev::operator="
			    <<"(const DrcOptDev&) name,copy: "
			    <<d.m_name<<" "<<d.m_copyNumber<<endl;
  if (&d != this) copy(d);
  return *this;

}
//----------------------------------------------------------------------
DrcOptDev::~DrcOptDev()
{
  list<DrcSurfAbs*>::const_iterator isurf;

  for(isurf=m_listSurf.begin(); isurf != m_listSurf.end(); ++isurf) 
    {
      delete (*isurf);
   }
}
//----------------------------------------------------------------------
void DrcOptDev::setCopyNumber(int i)
{
  m_copyNumber = i;

  list<DrcSurfAbs*>::const_iterator isurf;

  for(isurf=m_listSurf.begin(); isurf != m_listSurf.end(); ++isurf) 
    {
      (*isurf)->setCopyNumber(i);
   }


}
//----------------------------------------------------------------------
void DrcOptDev::addSurface(const DrcSurfAbs& surf)
{
  DrcSurfAbs* tmp = surf.clone();

  m_listSurf.push_back(tmp);
}
//----------------------------------------------------------------------
void DrcOptDev::print(fstream& stream) const
{
  list<DrcSurfAbs*>::const_iterator isurf;

  for(isurf=m_listSurf.begin(); isurf != m_listSurf.end(); ++isurf) 
    {
      (*isurf)->print(stream);
   }

  // prepare tracing of photons
  m_photonTrace = true;
  m_photonTraceStream = &stream;
}
//----------------------------------------------------------------------
DrcOptMatAbs& DrcOptDev::optMaterial() const
{
  static DrcOptMatAbs* dummy = 0;
  //  //cerr<<" *** DrcOptDev::optMaterial: This function must be overloaded when\n";
  //cerr<<"                             using volumes with refraction. Abort.\n";
  //exit(EXIT_FAILURE);
  return *dummy;
}
//----------------------------------------------------------------------
void DrcOptDev::addTransform(const Transform3D& trans)
{
  list<DrcSurfAbs*>::const_iterator i;
	      
  for (i=m_listSurf.begin(); 
       i != m_listSurf.end(); 
       ++i)
    {
      const_cast<DrcSurfAbs*>(*i)->addTransform(trans);
    }

}

//----------------------------------------------------------------------
void DrcOptDev::print()
{
  list<DrcSurfAbs*>::const_iterator i;
	    
  for (i=m_listSurf.begin(); 
       i != m_listSurf.end(); 
       ++i)
    {
      cout<<"\t\t\t"<<hex<<int(*i)<<" "<<(*i)->name()<<endl;
    }
}
//----------------------------------------------------------------------
void DrcOptDev::setPrintColor(int col)
{
  list<DrcSurfAbs*>::const_iterator i;
	    
  for (i=m_listSurf.begin(); 
       i != m_listSurf.end(); 
       ++i)
    {
      (*i)->setPrintColor(col);
    }

}