In This Package:

GenPruneTool.cc

Go to the documentation of this file.

#include "GenPruneTool.h"
#include "Context/TimeStamp.h"
#include "Event/GenHeader.h"

using namespace DayaBay;

GenPruneTool::GenPruneTool(const std::string& type,
                                         const std::string& name, 
                                         const IInterface* parent)
    : GaudiTool(type,name,parent)
{
  declareInterface< IGenPruneTool >(this) ;

  declareProperty("MaxVertices", m_maxVertices = -1,
                  "Negative: no touch. Otherwise, the number of vertices must be <= the MaxVertices");
}

GenPruneTool::~GenPruneTool()
{}

StatusCode GenPruneTool::initialize()
{
  return StatusCode::SUCCESS;
}

StatusCode GenPruneTool::finalize()
{
  return StatusCode::SUCCESS;
}


StatusCode GenPruneTool::prune(DayaBay::GenHeader* genHeader)
{
  if (m_maxVertices < 0) return StatusCode::SUCCESS;

  HepMC::GenEvent* event = genHeader->event();
  if (!event) {
    debug() << "GenHeader has null event, my work here is already done" << endreq;
    return StatusCode::SUCCESS;
  }

  int nverts = genHeader->event()->vertices_size();
  if (nverts <= m_maxVertices) {
    debug() << "GenHeader has " << nverts << " <= " << m_maxVertices
            << ", my work here is already done."
            << endreq;
    return StatusCode::SUCCESS;
  }

  genHeader->setEvent(0);
  debug() << "nullified GenHeader's event pointer" << endreq;
  return StatusCode::SUCCESS;
}

---