EEmcTTMatch.cxx

// $Id: EEmcTTMatch.cxx,v 1.7 2004/05/14 01:47:22 zolnie Exp $
#include <iostream>
#include <ostream>

#include "TList.h"

#include "StMuDSTMaker/COMMON/StMuEvent.h"
#include "StMuDSTMaker/COMMON/StMuTrack.h"

#include "EEmcTower.h"
#include "EEmcTTMatch.h"

#if !defined(ST_NO_NAMESPACES)
using std::ostream;
#endif


ClassImp(EEmcTTMatch);

//_____________________________________________________________________________
EEmcTTMatch::EEmcTTMatch()
{
  mTower  = NULL;
  mTracks = new TList();
}

//_____________________________________________________________________________
EEmcTTMatch::~EEmcTTMatch()
{ 
  if(mTracks) {
    mTracks->Clear();
    delete mTracks; 
  }
}

//_____________________________________________________________________________
void
EEmcTTMatch::Clear(Option_t *opt)
{
  mTower = NULL;
  mTracks->Clear();
}

//_____________________________________________________________________________
void 
EEmcTTMatch::Add(EEmcTower *t) 
{ 
  mTower = t; 
}

//_____________________________________________________________________________
void 
EEmcTTMatch::Add(StMuTrack *t) 
{
  mTracks->Add(t);
}


//_____________________________________________________________________________
Int_t
EEmcTTMatch::Matches()
{
  return mTracks->GetSize();
}


//_____________________________________________________________________________
Bool_t
EEmcTTMatch::ExtrapolateToZ(const StMuTrack *track, const double   z, TVector3 &r)
{
  const double kMinDipAngle   = 1.0e-13;

  StPhysicalHelixD   helix  = track->helix();
  double             dipAng = helix.dipAngle();
  double             z0     = helix.origin().z();
  if(dipAng<kMinDipAngle) return kFALSE; 
  double s  = ( z - z0 ) / sin(dipAng)  ;
  StThreeVectorD hit = helix.at(s);
  r.SetXYZ(hit.x(),hit.y(),hit.z());
  return   kTRUE;
}


//_____________________________________________________________________________
ostream& 
EEmcTTMatch::Out(ostream &out ) const
{
  out << *mTower ;
  StMuTrack *track;
  TIter nextTrack(mTracks);
  while((track=(StMuTrack *)nextTrack()))  out << *track;
  return out;
}




// ===========================================================================================
ostream&  operator<<(ostream &out, const StMuTrack    &t  )  {
  out << "<StMuTrack";
  out << " ORIGIN=\""   << t.helix().origin() << "\"";
  out << " MOMENTUM=\"" << t.momentum()   << "\"";
  out << "/>\n"; 
  return out;
}
// ===========================================================================================
ostream&  operator<<(ostream &out, const EEmcTTMatch  &m  )  {
  return m.Out(out);
}