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


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

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

// Collaborating Class Headers --------
//#include "CLHEP/GenericFunctions/Gaussian.hh"
#include "TpcDriftedElectron.h"


ClassImp(TpcAvalanche)

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

TpcAvalanche::TpcAvalanche()
  : fx(0),fy(0),ft(0),famp(0),
    fmcTrackId(0),fmcSecId(0),fmcHitId(0),findex(0),fDriftedElectronId(0)
{;}

TpcAvalanche::TpcAvalanche(const double X, 
			   const double Y,
			   const double T,
			   const double Amp,
			   TpcDriftedElectron* tpcDElectronp)
  : fx(X),fy(Y),ft(T),famp(Amp),
    fmcTrackId(tpcDElectronp->mcTrackId()),
    fmcSecId(tpcDElectronp->mcSecId()),
    fmcHitId(tpcDElectronp->mcHitId()),
    findex(0),fDriftedElectronId(tpcDElectronp->index())
{;}



bool operator== (const TpcAvalanche& lhs,const TpcAvalanche& rhs)
{
  double dx=lhs.fx - rhs.fx;
  double dy=lhs.fy - rhs.fy;
  double dt=lhs.ft - rhs.ft;
  double d2=dx*dx+dy*dy;
  return (d2<1E-10) && (fabs(dt)<1E-4);
}

std::ostream& operator<< (std::ostream& s, const TpcAvalanche& me){
  s << "TpcAvalanche:\n"
    << "   t    ="<<me.ft<<"\n"
    << "   x    ="<<me.fx<<"\n"
    << "   y    ="<<me.fy<<"\n"
    << "   Amp  ="<<me.famp<<"\n";
  return s;
}