//
// This is a rip-off of MultiTrace.  Code common between MultiTimeDistribution
// and MultiTrace needs to be factored out.
//
// $Id: MultiTimeDistribution.h 21160 2012-05-15 13:43:26Z javierg $

#ifndef _utl_MultiTimeDistribution_h_
#define _utl_MultiTimeDistribution_h_

#include <utl/TimeDistribution.h>
#include <utl/MultiObject.h>


namespace utl {


  template<typename T>
  class LabeledTimeDistribution : public LabeledObject<TimeDistribution<T> > {

  public:
    LabeledTimeDistribution(const int label, TimeDistribution<T>& td) :
      LabeledObject<TimeDistribution<T> >(label, td) { }

    LabeledTimeDistribution(const std::pair<int, TimeDistribution<T>*> pair) :
      LabeledObject<TimeDistribution<T> >(pair) { }

    TimeDistribution<T>& GetTimeDistribution()
    { return LabeledObject<TimeDistribution<T> >::GetObject(); }
    const TimeDistribution<T>& GetTimeDistribution() const
    { return LabeledObject<TimeDistribution<T> >::GetObject(); }

  };


  /** \brief A collection of TimeDistribution<T>,which provides methods to access
      different sources.

      \author Stefano Argiro', Darko Veberic
      \date Jan 2003
      \date 14 Nov 2005 revision
      \version $Id: MultiTimeDistribution.h 21160 2012-05-15 13:43:26Z javierg $

      <ul>
        <li> MultiTimeDistribution::GetTimeDistribution(sevt::StationConstants::eElectron) gives the electron component,
        <li> MultiTimeDistribution::GetTimeDistribution(sevt::StationConstants::eMuon) gives the muon component
        <li> MultiTimeDistribution::GetTimeDistribution() or
        <li> MultiTimeDistribution::GetTimeDistribution(sevt::StationConstants::eTotal) gives the total
      </ul>
      For the FD case examples are:
      <ul>
        <li> MultiTimeDistribution::GetTimeDistribution(fevt::FdConstants::eFluorDirect)
        <li> MultiTimeDistribution::GetTimeDistribution(fevt::FdConstants::eCherRayleighScattered)
      </ul>

      \param T is the type stored in the TimeDistribution (double, int etc)
      \warning: this class is for developer convenience, not
      meant for the average user
      \ingroup stl time
  */

  template<typename T>
  class MultiTimeDistribution :
    public MultiObject<TimeDistribution<T>, LabeledTimeDistribution<T> > {

  public:
    /// Returns the trace for /par source
    TimeDistribution<T>& GetTimeDistribution(const int label = 0)
    { return MultiObject<TimeDistribution<T>, LabeledTimeDistribution<T> >::GetObject(label); }
    const TimeDistribution<T>& GetTimeDistribution(const int label = 0) const
    { return MultiObject<TimeDistribution<T>, LabeledTimeDistribution<T> >::GetObject(label); }

    /** Creates a new TimeDistribution of size /par size and binning
     * /par binsize and assign source /par source */
    void AddTimeDistribution(const double binSize, const int label)
    { this->AddObject(new TimeDistribution<T>(binSize), label); }

    /** Copies trace in the internal list
     *  and assign source /par source -The original can be deleted */
    void AddTimeDistribution(const TimeDistribution<T>& td, const int label)
    { this->AddObject(td, label); }

  };

  // To keep the forward declarations consistent
  // modify TimeDistribution-fwd when chaniging/adding below

  typedef utl::MultiTimeDistribution<double> MultiTimeDistributionD;
  typedef utl::MultiTimeDistribution<int> MultiTimeDistributionI;

}


#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// compile-command: "cd ..; make -k"
// End: