//
// $Id: SiCenter.cpp,v 1.8 2001/02/03 22:53:35 hobbs Exp $
//
// File: SiCenter.cpp
// Purpose: 
// Created: 17-JAN-1998   Emanuela Barberis
//
// $Revision: 1.8 $
//
//
#include <iostream>
#include "SiCenter.hpp"
#include <geometry_system/components/CartesianCoordinate.hpp>

// Global definitions
using namespace std;
using namespace dgs;

bool operator<(const SiCenter& lhs, const SiCenter &rhs) {
  return lhs.get_position()[2] < rhs.get_position()[2];
}
// Constructors/Destructors
SiCenter::SiCenter()
{
}

SiCenter::SiCenter(const vector<SiCentral>& centrals,
		   const vector<float>& z_positions): _centrals(centrals)
{

  for( int ic=0 ; ic<_centrals.size() ; ic++ ) {
    CartesianCoordinate myz(0.0,0.0,z_positions[ic]);
    position_child_at(myz,_centrals[ic]);
    int nb = _centrals[ic].get_barrel_count();
    int nd = _centrals[ic].get_disk_count();
    int k;
    
    for(k=1; k<=nb; k++) {
       _barrel_map_centralID.push_back(ic+1);
       _barrel_map_barrelID.push_back(k);
    }

    for( k=1; k<=nd; k++) {
       _disk_map_centralID.push_back(ic+1);
       _disk_map_diskID.push_back(k);
    }

  }

}
const list<GeometryElement*> SiCenter::get_children()
{
  list<GeometryElement*> kids;
  for( int i=0 ; i<_centrals.size() ; i++ ) kids.push_back( &_centrals[i] );
  return kids;
}

std::ostream& operator <<(std::ostream& os, const SiCenter& me)
{  int i;
  os << "SiCenter at " << me.get_position() << " has " << me._centrals.size() 
     << " centrals" << endl;
  for(i=0 ; i<me._centrals.size() ; i++)  os << "  " << me._centrals[i];

  return os;
}


bool operator ==(const SiCenter& lhs, const SiCenter &rhs)
{
  if( lhs.get_position() != rhs.get_position() ) return false;
  if( lhs._centrals.size() != rhs._centrals.size() ) return false;
  for( int i=0 ; i<lhs._centrals.size() ; i++ ) 
    if( lhs._centrals[i] != rhs._centrals[i] ) return false;
  return true;
}

bool operator !=(const SiCenter& lhs, const SiCenter &rhs) {
  return !(lhs==rhs);
}