#include <GeneratorInterface/GenFilters/src/ComphepSingletopFilter.cc>

Inheritance diagram for ComphepSingletopFilter:

Public Member Functions
	ComphepSingletopFilter (const edm::ParameterSet &)

	~ComphepSingletopFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

virtual bool	filter (edm::Event &, const edm::EventSetup &)

Private Attributes
edm::InputTag	hepMCProductTag_

unsigned int	pass22

unsigned int	pass23

double	ptsep

unsigned int	read22

unsigned int	read23

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: a filter to match LO/NLO in comphep-generated singletop

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

Definition at line 17 of file ComphepSingletopFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file ComphepSingletopFilter.cc.

References edm::ParameterSet::getParameter(), and ptsep.

                                                                              :
     hepMCProductTag_(iConfig.getParameter<edm::InputTag>("hepMCProductTag")) {
     ptsep = iConfig.getParameter<double>("pTSep");
 }

ComphepSingletopFilter::~ComphepSingletopFilter ( )

Definition at line 39 of file ComphepSingletopFilter.cc.

39 {}

Member Function Documentation

void ComphepSingletopFilter::beginJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 42 of file ComphepSingletopFilter.cc.

References pass22, pass23, read22, and read23.

                                       {
     read22 = read23 = 0;
     pass22 = pass23 = 0;
 }

void ComphepSingletopFilter::endJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 48 of file ComphepSingletopFilter.cc.

References gather_cfg::cout, cmsPerfStripChart::format, pass22, pass23, read22, and read23.

                                     {
     using namespace std;
     cout << "Proc:     2-->2     2-->3     Total" << endl;
     cout << boost::format("Read: %9d %9d %9d") % read22 % read23 % (read22+read23)
          << endl;
     cout << boost::format("Pass: %9d %9d %9d") % pass22 % pass23 % (pass22+pass23)
          << endl;
 }

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

privatevirtual

Implements edm::EDFilter.

Definition at line 58 of file ComphepSingletopFilter.cc.

References funct::abs(), ecal_dqm_sourceclient-live_cfg::cerr, GenParticle::GenParticle, edm::Event::getByLabel(), hepMCProductTag_, NULL, parents, pass22, pass23, EnergyCorrector::pt, ptsep, read22, and read23.

                                  {
 
     using namespace std;
     using namespace HepMC;
   
     edm::Handle<edm::HepMCProduct> evt;
     iEvent.getByLabel(hepMCProductTag_, evt);
     const HepMC::GenEvent * myEvt = evt->GetEvent();
 //  myEvt->print();  // to print the record
 
     const GenParticle * gp_clep = NULL;
 
     for (GenEvent::particle_const_iterator it = myEvt->particles_begin();
          it != myEvt->particles_end(); ++it) {
         if ((*it)->status() == 3) {
             int abs_id = abs((*it)->pdg_id());
             if (abs_id==11 || abs_id==13 || abs_id==15) {
                 gp_clep = *it;
                 break;
             }
         }
     }
 
     if (! gp_clep) {
         cerr << "ERROR: ComphepSingletopFilter: no charged lepton" << endl;
         return false;
     }
 
     int id_bdec = 0;
     vector<const GenParticle *> vgp_bsec;
 
     GenVertex * gv_hard = gp_clep->production_vertex();
 
     for (GenVertex::particle_iterator it = gv_hard->particles_begin(children);
          it != gv_hard->particles_end(children); ++it) {
         int pdg_id = (*it)->pdg_id();
         if (abs(pdg_id) == 5) {
             if (pdg_id * (gp_clep->pdg_id()) < 0) {
                 id_bdec = pdg_id;
             } else {
                 vgp_bsec.push_back(*it);
             }
         }
     }
 
     bool process22 = (vgp_bsec.size() == 0);
 
     GenVertex * gv = NULL;
     if (process22) {
         for (GenVertex::particle_iterator it = gv_hard->particles_begin(parents);
              it != gv_hard->particles_end(parents); ++it) {
             if ((*it)->pdg_id() == id_bdec) {
                 gv = (*it)->production_vertex();
                 break;
             }
         }
         if (! gv) {
             cerr << "ERROR: ComphepSingletopFilter: HepMC inconsistency" << endl;
             myEvt->print();
             return false;
         }
     } else {
         gv = vgp_bsec.back()->end_vertex();
     }
     const GenParticle * gp;
     while (gv) {
         gp = NULL;
         for (GenVertex::particle_iterator it = gv->particles_begin(children);
              it != gv->particles_end(children); ++it) {
             if ((*it)->pdg_id() == -id_bdec) {
                 if (!gp || (*it)->momentum().perp2() > gp->momentum().perp2()) {
                     gp = *it;
                 }
             }
         }
         if (gp) {
             gv = gp->end_vertex();
             vgp_bsec.push_back(gp);
         } else {
             gv = NULL;
         }
     }
 
     if (vgp_bsec.size() == 0) {
         cerr << "ERROR: ComphepSingletopFilter: HepMC inconsistency" << endl;
         return false;
     }
 
     double pt = vgp_bsec.back()->momentum().perp();
     bool pass;
     if (process22) {
         read22 += 1;
         pass = pt < ptsep;
         if (pass) pass22 += 1;
     } else {
         read23 += 1;
         pass = ptsep <= pt;
         if (pass) pass23 += 1;
     }
 
     return pass;
 }