/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/GeneratorInterface/HiGenCommon/src/MultiCandGenEvtSelector.cc

Go to the documentation of this file.

#include <iostream>
#include "GeneratorInterface/HiGenCommon/interface/MultiCandGenEvtSelector.h"

MultiCandGenEvtSelector::MultiCandGenEvtSelector(const edm::ParameterSet& iConfig)
   : BaseHiGenEvtSelector(iConfig)
{
   ptMin_ = iConfig.getParameter<double>("ptMin");
   etaMax_ = iConfig.getParameter<double>("etaMax");
   pdg_ = iConfig.getParameter<int>("pdg");
   st_ = iConfig.getParameter<int>("status");
   nTrig_ = iConfig.getParameter<int>("minimumCandidates");
}

bool MultiCandGenEvtSelector::filter(HepMC::GenEvent * evt){
   std::cout<<"Di Muon Fired"<<std::endl;
   
   int found = 0;
   HepMC::GenEvent::particle_const_iterator begin = evt->particles_begin();
   HepMC::GenEvent::particle_const_iterator end = evt->particles_end();
   for(HepMC::GenEvent::particle_const_iterator it = begin; it != end; ++it){
      if(selectParticle(*it, st_, pdg_, ptMin_, etaMax_)) found++;
      if(found == nTrig_) return true;      
   }
   
   return false;
}