DNeutralParticle_factory.cc

Go to the documentation of this file.

// $Id$
//
//    File: DNeutralParticle_factory.cc
// Created: Tue Aug  9 14:29:24 EST 2011
// Creator: pmatt (on Linux ifarml6 2.6.18-128.el5 x86_64)
//


#include <iostream>
#include <iomanip>
using namespace std;

#include "DNeutralParticle_factory.h"
using namespace jana;

//------------------
// init
//------------------
jerror_t DNeutralParticle_factory::init(void)
{
   return NOERROR;
}

//------------------
// brun
//------------------
jerror_t DNeutralParticle_factory::brun(jana::JEventLoop *locEventLoop, int32_t runnumber)
{
   return NOERROR;
}

//------------------
// evnt
//------------------
jerror_t DNeutralParticle_factory::evnt(jana::JEventLoop *locEventLoop, uint64_t eventnumber)
{
   vector<const DNeutralParticleHypothesis*> locNeutralParticleHypotheses;
   locEventLoop->Get(locNeutralParticleHypotheses);

   map<const DNeutralShower*, vector<const DNeutralParticleHypothesis*> > locHypothesesByShower;
   for(size_t loc_i = 0; loc_i < locNeutralParticleHypotheses.size(); loc_i++)
      locHypothesesByShower[locNeutralParticleHypotheses[loc_i]->Get_NeutralShower()].push_back(locNeutralParticleHypotheses[loc_i]);

   for(auto& locPair : locHypothesesByShower)
   {
      DNeutralParticle* locNeutralParticle = new DNeutralParticle();
      locNeutralParticle->dNeutralShower = locPair.first;
      locNeutralParticle->dNeutralParticleHypotheses = locPair.second;
      _data.push_back(locNeutralParticle);
   }

   return NOERROR;
}

//------------------
// erun
//------------------
jerror_t DNeutralParticle_factory::erun(void)
{
   return NOERROR;
}

//------------------
// fini
//------------------
jerror_t DNeutralParticle_factory::fini(void)
{
   return NOERROR;
}