PhysicsNtuple/Event.cxx

Go to the documentation of this file.

// $Id: Event.cxx,v 1.24 2009/10/15 18:30:13 jyuko Exp $

// C++
#include <cassert>
#include <climits>
#include <cmath>
#include <iostream>

// Local
#include "PhysicsNtuple/Shower.h"
#include "PhysicsNtuple/Track.h"

// Local
#include "PhysicsNtuple/Event.h"

using namespace std;

//-------------------------------------------------------------------
Anp::Event::Event()
   :fBasic(),
    fNu(),
    fVtx(),
    index(-1000),
    gev(-1.0e6),
    pid(-1.0e6),
    fData(),
    fShower(),
    fTrack()
{
}

//-------------------------------------------------------------------
Anp::Event::~Event()
{
}

//----------------------------------------------------------------------
void Anp::Event::Clear()
{
   fBasic.Clear();
   fNu.Clear();
   fVtx.Clear();

   index = -10000;
   gev   = -1.0e6;
   pid   = -1.0e6;

   fData.clear();
   fShower.clear();
   fTrack.clear();
}

//----------------------------------------------------------------------
void Anp::Event::ClearData()
{ 
   fData.clear();
}

//---------------------------------------------------------------------- 
bool Anp::Event::Add(const int key, const float data)
{ 
   //
   // Add data at key, if key does not already exits 
   //

   assert(std::abs(key) < SHRT_MAX && "key is out of range");
   
   if(KeyExists(key))
   {
      cerr << "Event::Add(" << key << ") - key already exists" << endl;
      return false;
   }
   
   fData.push_back(Anp::Data(static_cast<short>(key), data));
   return true;
}

//----------------------------------------------------------------------
bool Anp::Event::FirstShower(const short index) const
{
   if(fShower.empty())
   {
      return false;
   }
   
   return (fShower.front() == index);
}

//----------------------------------------------------------------------
float Anp::Event::DataAt(const short key) const
{
   //
   // Return data at key, if key exits 
   //

   const DataIter it = find(fData.begin(), fData.end(), key);
   if(it == fData.end())
   {
      return -999.0;
   }
   
   return it -> Data();
}

//----------------------------------------------------------------------
void Anp::Event::Print(ostream &os) const
{
   os << "Event::Print" << endl
      << "   index = " << index << endl
      << "   gev = " << gev << endl
      << "   nshower = " << GetNShowers() << endl
      << "   ntrack = " << GetNTracks() << endl;
   
   for(DataVec::const_iterator dit = fData.begin(); dit != fData.end(); ++dit)
   {
      os  << "   (key, data) = (" << dit -> Key() << ", " << dit -> Data() << ")" << endl;
   }
   
   fBasic.Print(os);
   fNu.Print(os);
   fVtx.Print(os);
}

//----------------------------------------------------------------------
bool Anp::operator==(const Event  &event, const Shower &shower)
{
   return event.MatchShower(shower.ShowerIndex());
}

//----------------------------------------------------------------------
bool Anp::operator==(const Shower &shower, const Event &event)
{
   return event.MatchShower(shower.ShowerIndex());
}

//----------------------------------------------------------------------
bool Anp::operator==(const Track &track, const Event &event)
{
   return event.MatchTrack(track.TrackIndex());
}

//----------------------------------------------------------------------
bool Anp::operator==(const Event &event, const Track &track)
{
   return event.MatchTrack(track.TrackIndex());
}

---