PartonShowerCsHepMCFilter.cc

Go to the documentation of this file.

#include "GeneratorInterface/Core/interface/PartonShowerCsHepMCFilter.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include <iostream>
#include "HepPDT/ParticleID.hh"


using namespace edm;
using namespace std;


//constructor
PartonShowerCsHepMCFilter::PartonShowerCsHepMCFilter(const edm::ParameterSet& iConfig) 
{

}


//destructor
PartonShowerCsHepMCFilter::~PartonShowerCsHepMCFilter()
{

}

//
// member functions
//

// ------------ method called to produce the data  ------------
bool PartonShowerCsHepMCFilter::filter(const HepMC::GenEvent* evt)
{
  
  // loop over gen particles
  for ( HepMC::GenEvent::particle_const_iterator p = evt->particles_begin(); p != evt->particles_end(); ++p ){
    
    // check only status 2 particles
    if( (*p)->status()==2 ){
      // if one of the status 2 particles is a C-hadron, accept the event
      HepPDT::ParticleID pid((*p)->pdg_id());
      if( pid.hasCharm() ){
        return true; // accept event
      }
    }
    
  }

  return false; // skip event

}