//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Implementation of class McIdCollection
//      see McIdCollection.hh for details
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Sebastian Neubert    TUM            (original author)
//
//
//-----------------------------------------------------------

// This Class' Header ------------------
#include "McIdCollection.h"

// C/C++ Headers ----------------------
#include <algorithm>

// Collaborating Class Headers --------


// Class Member definitions -----------

ClassImp(McIdCollection)

McIdCollection::McIdCollection(){;}

McIdCollection::~McIdCollection()
{
  _MCIDs.clear();
}

McIdCollection::McIdCollection(const McIdCollection& coll)
 : _MCIDs(coll._MCIDs)
{}


bool
McIdCollection::AddID(const McId& ID){
  vector<McId>::iterator it=find(_MCIDs.begin(),_MCIDs.end(),ID);
  if(it==_MCIDs.end()){
    _MCIDs.push_back(ID);  
    return true;
  }
  else {
    (*it).IncWeight(ID.weight());
    return false;
  }
};


bool
McIdCollection::AddIDCollection(const McIdCollection& coll){
  unsigned int n=coll.nIDs();
  bool ret=false;
  for(unsigned int i=0; i<n; ++i)ret|=AddID(coll.ID(i));
  return ret;
}


unsigned int
McIdCollection::TotalWeight() const {
  unsigned int n=nIDs();
  unsigned int result=0;
  for(unsigned int i=0; i<n;++i){
    result+=_MCIDs[i].weight();
  }
  return result;
}

double
McIdCollection::MaxRelWeight() const {
  double result=0;
  double total=0;
  unsigned int n=nIDs();
  for(unsigned int i=0; i<n;++i){
    double w=(double)_MCIDs[i].weight();
    result=std::max(result,w);
    total+=w;
  }
  return result/total;
}

McId 
McIdCollection::DominantID() const {
  unsigned int n=nIDs();
  unsigned int index=0;
  unsigned int maxw=0;
  for(unsigned int i=0; i<n;++i){
    unsigned int w=_MCIDs[i].weight();
    if(maxw<w){
      maxw=w;
      index=i;
    }
  }
  return _MCIDs[index];
}

std::ostream& operator<< (std::ostream& s, const McIdCollection& me){
  s << "MCID-Collection:\n";
  for(unsigned int i=0; i<me._MCIDs.size(); ++i)s<<me._MCIDs[i]<<std::endl;
  return s;
}