PythiaFilterTTBar Class Reference

#include <GeneratorInterface/GenFilter/src/PythiaFilterTTBar.cc>

Inheritance diagram for PythiaFilterTTBar:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	PythiaFilterTTBar (const edm::ParameterSet &)
	~PythiaFilterTTBar ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Attributes
unsigned int	decayType_
std::string	label_
unsigned int	leptonFlavour_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
typedef WorkerT< EDFilter >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Description: edmFilter to select a TTBar decay channel

Implementation: decayType: 1 + leptonFlavour: 0 -> Semi-leptonic leptonFlavour: 1 -> Semi-e leptonFlavour: 2 -> Semi-mu leptonFlavour: 3 -> Semi-tau decayType: 2 -> di-leptonic (no seperate channels implemented yet)

decayType: 3 -> fully-hadronic

Definition at line 56 of file PythiaFilterTTBar.h.

Constructor & Destructor Documentation

PythiaFilterTTBar::PythiaFilterTTBar ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 16 of file PythiaFilterTTBar.cc.

                                                                   : 
  label_(iConfig.getUntrackedParameter("moduleLabel",std::string("generator"))),
  decayType_(iConfig.getUntrackedParameter("decayType",1)),
  leptonFlavour_(iConfig.getUntrackedParameter("leptonFlavour",0))       
{
}

PythiaFilterTTBar::~PythiaFilterTTBar

(

)

Definition at line 24 of file PythiaFilterTTBar.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

bool PythiaFilterTTBar::filter ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDFilter.

Definition at line 39 of file PythiaFilterTTBar.cc.

References accept(), decayType_, edm::Event::getByLabel(), label_, leptonFlavour_, AlCaHLTBitMon_ParallelJobs::p, and ntuplemaker::status.

{

  bool accept=false;

  edm::Handle<edm::HepMCProduct> evt;
  iEvent.getByLabel(label_, evt);

  const HepMC::GenEvent * myGenEvent = evt->GetEvent();

  unsigned int iE=0, iMu=0, iTau=0;

  unsigned int iNuE=0, iNuMu=0, iNuTau=0;

  unsigned int iLep=0, iNu=0;


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

    int pdgID = (*p)->pdg_id();
    
    int status = (*p)->status();

    if ( status == 3 ) {

      // count the final state leptons

      if (fabs(pdgID) == 11)
    iE++;

      if(fabs(pdgID) == 13)
    iMu++;
      
      if(fabs(pdgID) == 15)
    iTau++;

      // count the final state neutrinos
      
      if(fabs(pdgID) == 12)
    iNuE++;

      if(fabs(pdgID) == 14)
    iNuMu++;
      
      if(fabs(pdgID) == 16)
    iNuTau++;

    }

  }

  iLep = (iE+iMu+iTau);
  iNu = (iNuE+iNuMu+iNuTau);

  if (decayType_ == 1) { // semi-leptonic decay
    
    // l = e,mu,tau

    if (leptonFlavour_ == 0 && iLep == 1 && iNu == 1)
      accept=true;
    
    // l = e

    else if (leptonFlavour_ == 1 && iE == 1 && iNuE == 1 && iLep == 1 && iNu == 1)
      accept=true;

    // l = mu

    else if (leptonFlavour_ == 2 && iMu == 1 && iNuMu == 1 && iLep == 1 && iNu == 1)
      accept=true;
    
    // l = tau 

    else if (leptonFlavour_ == 3 && iTau == 1 && iNuTau == 1 && iLep == 1 && iNu == 1)
      accept=true;

  }

  else if (decayType_ == 2) { // di-leptonic decay (inclusive)

    if (iLep == 2 && iNu == 2) 
      accept=true;
    
  }
  
  else if (decayType_ == 3) { // fully-hadronic decay

    if (iLep == 0 && iNu == 0)
      accept=true;
  }

  else
    accept=false;
  
  
  return accept;
}

Member Data Documentation

unsigned int PythiaFilterTTBar::decayType_

private

Definition at line 67 of file PythiaFilterTTBar.h.

Referenced by filter().

std::string PythiaFilterTTBar::label_

private

Definition at line 65 of file PythiaFilterTTBar.h.

Referenced by Mixins._Labelable::_findDependencies(), Mixins._TypedParameterizable::dumpPython(), and filter().

unsigned int PythiaFilterTTBar::leptonFlavour_

private

Definition at line 69 of file PythiaFilterTTBar.h.

Referenced by filter().