CEvent/StationSimData.cc

Go to the documentation of this file.

#include <cevt/StationSimData.h>
// #include <cevt/StationTriggerData.h>
// #include <cevt/StationGPSData.h>
#include <utl/Particle.h>
#include <utl/ShadowPtr.h>
#include <sstream>
#include <utl/ErrorLogger.h>

using namespace utl;
using namespace cevt;
using namespace std;


StationSimData::StationSimData() :
  fNumberOfMuons(0), 
  fNumberOfElectrons(0), 
  fNumberOfElectronsFromMuons(0),
  fNumberOfElectronsFromPhotons(0)
{
}


StationSimData::~StationSimData()
{
}


void
StationSimData::AddParticle(const utl::Particle& particle)
{
  fParticles.push_back(ParticlePtr(new Particle(particle)));
  if (particle.GetType() == utl::Particle::eMuon || particle.GetType() == utl::Particle::eAntiMuon )
    ++fNumberOfMuons;

  if (particle.GetType() == utl::Particle::eElectron || particle.GetType() == utl::Particle::ePositron ) {
    ++fNumberOfElectrons;

    if (particle.HasParent()) {
      if (particle.GetParent().GetType() == utl::Particle::eMuon || 
          particle.GetParent().GetType() == utl::Particle::eAntiMuon)
        ++fNumberOfElectronsFromMuons;

      if (particle.GetParent().GetType() == utl::Particle::ePhoton) 
        ++fNumberOfElectronsFromPhotons;
    }
  }
}


void
StationSimData::ClearParticleList()
{
  fParticles.clear();
}


void
StationSimData::AddInjectedParticle(const utl::Particle& particle)
{
  fInjectedParticles.push_back(ParticlePtr(new Particle(particle)));
}


void
StationSimData::ClearInjectedParticleList()
{
  fInjectedParticles.clear();
}


// CBT
// void
// StationSimData::MakeTriggerAndGPSData(const TimeStamp& lsTime)
// {
//   if (HasTriggerData(lsTime)) {
//     WARNING("Trigger and GPS data already present!");
//     return;
//   }

//   TriggerGPSPair tg;
//   tg.first = new StationTriggerData;
//   tg.second = new StationGPSData;
//   fTriggerAndGPSData.insert(make_pair(lsTime, tg));
//   fTriggerTimes.push_back(lsTime);
// }


// StationGPSData&
// StationSimData::GetGPSData(const TimeStamp& lsTime)
// {
//   if (!HasTriggerData(lsTime)) {
//     ostringstream err;
//     err << "Tried to get unavailable GPS data at time " << lsTime;
//     ERROR(err);
//     throw NonExistentComponentException(err.str());
//   }

//   TriggerGPSPair& tg = fTriggerAndGPSData.find(lsTime)->second;
//   return *tg.second;
// }


// const StationGPSData&
// StationSimData::GetGPSData(const TimeStamp& lsTime)
//   const
// {
//   if (!HasTriggerData(lsTime)) {
//     ostringstream err;
//     err << "Tried to get unavailable GPS data at time " << lsTime;
//     ERROR(err);
//     throw NonExistentComponentException(err.str());
//   }

//   const TriggerGPSPair& tg = fTriggerAndGPSData.find(lsTime)->second;
//   return *tg.second;
// }


// StationTriggerData&
// StationSimData::GetTriggerData(const TimeStamp& lsTime)
// {
//   if (!HasTriggerData(lsTime)) {
//     ostringstream err;
//     err << "Tried to get unavailable trigger data at time " << lsTime;
//     ERROR(err);
//     throw NonExistentComponentException(err.str());
//   }

//   TriggerGPSPair& tg = fTriggerAndGPSData.find(lsTime)->second;
//   return *tg.first;
// }


// const StationTriggerData&
// StationSimData::GetTriggerData(const TimeStamp& lsTime)
//   const
// {
//   if (!HasTriggerData(lsTime)) {
//     ostringstream err;
//     err << "Tried to get unavailable trigger data at time " << lsTime;
//     ERROR(err);
//     throw NonExistentComponentException(err.str());
//   }

//   const TriggerGPSPair& tg = fTriggerAndGPSData.find(lsTime)->second;
//   return *tg.first;
// }


// bool
// StationSimData::HasTriggerData(const TimeStamp& lsTime)
//   const
// {
//   const ConstTriggerTimeIterator end = fTriggerTimes.end();
//   return find(fTriggerTimes.begin(), end, lsTime) != end;
// }


// bool
// StationSimData::HasGPSData(const TimeStamp& lsTime)
//   const
// {
//   const ConstTriggerTimeIterator end = fTriggerTimes.end();
//   return find(fTriggerTimes.begin(), end, lsTime) != end;
// }


// End CBT

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