GenericDauHepMCFilter Class Reference

#include <GenericDauHepMCFilter.h>

Inheritance diagram for GenericDauHepMCFilter:

Public Member Functions
virtual bool	filter (const HepMC::GenEvent *evt)
	GenericDauHepMCFilter (const edm::ParameterSet &)
	~GenericDauHepMCFilter ()
Public Member Functions inherited from BaseHepMCFilter
	BaseHepMCFilter ()
virtual	~BaseHepMCFilter ()

Private Attributes
bool	chargeconju
std::vector< int >	dauIDs
double	maxetacut
double	maxptcut
double	minetacut
double	minptcut
int	ndaughters
int	particleID

Detailed Description

Description: Filter events using MotherId and ChildrenIds infos

Implementation: <Notes on="" implementation>="">

Definition at line 35 of file GenericDauHepMCFilter.h.

Constructor & Destructor Documentation

GenericDauHepMCFilter::GenericDauHepMCFilter

(

const edm::ParameterSet &

iConfig

)

Definition at line 12 of file GenericDauHepMCFilter.cc.

                                                                           :
        particleID(iConfig.getParameter<int>("ParticleID")),
        chargeconju(iConfig.getParameter<bool>("ChargeConjugation")),
        ndaughters(iConfig.getParameter<int>("NumberDaughters")),
        dauIDs(iConfig.getParameter<std::vector<int> >("DaughterIDs")),
        minptcut(iConfig.getParameter<double>("MinPt")),
        maxptcut(iConfig.getParameter<double>("MaxPt")),
        minetacut(iConfig.getParameter<double>("MinEta")),
        maxetacut(iConfig.getParameter<double>("MaxEta"))
{

}

GenericDauHepMCFilter::~GenericDauHepMCFilter

(

)

Definition at line 26 of file GenericDauHepMCFilter.cc.

{

}

Member Function Documentation

bool GenericDauHepMCFilter::filter ( const HepMC::GenEvent *  evt )

virtual

Implements BaseHepMCFilter.

Definition at line 37 of file GenericDauHepMCFilter.cc.

References chargeconju, dauIDs, i, maxetacut, maxptcut, minetacut, minptcut, ndaughters, AlCaHLTBitMon_ParallelJobs::p, and particleID.

{
  
  for ( HepMC::GenEvent::particle_const_iterator p = evt->particles_begin();
                  p != evt->particles_end(); ++p ) {

          if( (*p)->pdg_id() != particleID ) continue ;

          int ndauac = 0;
          int ndau = 0;     
          if ( (*p)->end_vertex() ) {     
                  for ( HepMC::GenVertex::particle_iterator 
                                  des=(*p)->end_vertex()->particles_begin(HepMC::children);
                                  des != (*p)->end_vertex()->particles_end(HepMC::children);
                                  ++des ) {
                          ++ndau;       
                          //       cout << "   -> Daughter = " << (*des)->pdg_id() << endl;
                          for( unsigned int i=0; i<dauIDs.size(); ++i) {
                                  if( (*des)->pdg_id() != dauIDs[i] ) continue ;
                                  if(   (*des)->momentum().perp() >  minptcut  &&
                                                  (*des)->momentum().perp() <  maxptcut  &&
                                                  (*des)->momentum().eta()  >  minetacut && 
                                                  (*des)->momentum().eta()  <  maxetacut ) {
                                          ++ndauac;
                                          //           cout << "     -> pT = " << (*des)->momentum().perp() << endl;
                                          break;
                                  } 
                          }                            
                  }
          }  
          if( ndau ==  ndaughters && ndauac == ndaughters ) {
                  return true;
          }    

  }


  if (chargeconju) {

          for ( HepMC::GenEvent::particle_const_iterator p = evt->particles_begin();
                          p != evt->particles_end(); ++p ) {

                  if( (*p)->pdg_id() != -particleID ) continue ;
                  int ndauac = 0;
                  int ndau = 0;     
                  if ( (*p)->end_vertex() ) {
                          for ( HepMC::GenVertex::particle_iterator 
                                          des=(*p)->end_vertex()->particles_begin(HepMC::children);
                                          des != (*p)->end_vertex()->particles_end(HepMC::children);
                                          ++des ) {
                                  ++ndau;
                                  for( unsigned int i=0; i<dauIDs.size(); ++i) {
                                          bool has_antipart = !(dauIDs[i]==22 || dauIDs[i]==23);
                                          int IDanti = has_antipart ? -dauIDs[i] : dauIDs[i];
                                          if( (*des)->pdg_id() != IDanti ) continue ;
                                          if(   (*des)->momentum().perp() >  minptcut  &&
                                                          (*des)->momentum().perp() <  maxptcut  &&
                                                          (*des)->momentum().eta()  >  minetacut && 
                                                          (*des)->momentum().eta()  <  maxetacut ) {
                                                  ++ndauac;
                                                  break;
                                          } 
                                  }                            
                          }
                  }
                  if( ndau ==  ndaughters && ndauac == ndaughters ) {
                          return true;
                  }    
          }

  }    

  return false;

}

Member Data Documentation

bool GenericDauHepMCFilter::chargeconju

private

Definition at line 48 of file GenericDauHepMCFilter.h.

Referenced by filter().

std::vector<int> GenericDauHepMCFilter::dauIDs

private

Definition at line 50 of file GenericDauHepMCFilter.h.

Referenced by filter().

double GenericDauHepMCFilter::maxetacut

private

Definition at line 54 of file GenericDauHepMCFilter.h.

Referenced by filter().

double GenericDauHepMCFilter::maxptcut

private

Definition at line 52 of file GenericDauHepMCFilter.h.

Referenced by filter().

double GenericDauHepMCFilter::minetacut

private

Definition at line 53 of file GenericDauHepMCFilter.h.

Referenced by filter().

double GenericDauHepMCFilter::minptcut

private

Definition at line 51 of file GenericDauHepMCFilter.h.

Referenced by filter().

int GenericDauHepMCFilter::ndaughters

private

Definition at line 49 of file GenericDauHepMCFilter.h.

Referenced by filter().

int GenericDauHepMCFilter::particleID

private

Definition at line 47 of file GenericDauHepMCFilter.h.

Referenced by filter().