16 {
17
18 const McEventCollection* mcCollptr;
19 StatusCode stat = this->
getDH( mcCollptr );
20 if ( !stat ) { return StatusCode::FAILURE; }
21
22
23 McEventCollection::const_iterator itr;
24 for ( itr = mcCollptr->begin(); itr != mcCollptr->end(); ++itr )
25 {
26
27 const HepMC::GenEvent* genEvt = ( *itr );
28 if ( genEvt == 0 ) return StatusCode::FAILURE;
29
30 HepMC::GenEvent::particle_const_iterator it = genEvt->particles_begin();
31 HepMC::GenEvent::particle_const_iterator en = genEvt->particles_end();
32 for ( ; it != en; ++it )
33 {
34 if ( ( *selector )( *it ) ) { mcParticles.push_back( *it ); }
35 }
36 }
37
38 return StatusCode::SUCCESS;
39}
StatusCode getDH(const McEventCollection *&) const