#include <XtoFFbarFilter.h>

Inheritance diagram for XtoFFbarFilter:

Public Member Functions
virtual bool	filter (edm::Event &iEvent, const edm::EventSetup &iSetup)

	XtoFFbarFilter (const edm::ParameterSet &)

	~XtoFFbarFilter ()

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	endJob ()

bool	found (const std::vector< int > &v, int j)

bool	foundXtoFFbar (const reco::GenParticleRef &moth, const std::vector< int > &idMotherX, const std::vector< int > &idDaughterF)

Private Attributes
edm::Handle < reco::GenParticleCollection >	genParticles_

std::vector< int >	idDaughterF_

std::vector< int >	idDaughterG_

std::vector< int >	idMotherX_

std::vector< int >	idMotherY_

unsigned int	rejectedEvents_

bool	requireY_

edm::InputTag	src_

unsigned int	totalEvents_

double	xSumCtau_

double	xSumPt_

double	xSumR_

unsigned int	xTotal_

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

XtoFFbarFilter A GenParticle-based filter that searches for X -> f fbar where X and f are any particle you choose. Can optionally also require a second decay Y -> g g-bar to be present in same event.

Author: Ian Tomalin, RAL

Definition at line 30 of file XtoFFbarFilter.h.

Constructor & Destructor Documentation

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

Definition at line 6 of file XtoFFbarFilter.cc.

References idDaughterG_, idMotherY_, and requireY_.

                                                              :
   src_(iConfig.getParameter<edm::InputTag>("src")),
   idMotherX_(iConfig.getParameter<vector<int> >("idMotherX")),
   idDaughterF_(iConfig.getParameter<vector<int> >("idDaughterF")),
   idMotherY_(iConfig.getParameter<vector<int> >("idMotherY")),
   idDaughterG_(iConfig.getParameter<vector<int> >("idDaughterG")),
   xTotal_(0), xSumPt_(0.), xSumR_(0.), xSumCtau_(0.), totalEvents_(0),rejectedEvents_(0)
 { 
   // Note if if not searching for Y --> g-gbar.
   requireY_ = (idMotherY_.size() > 0 && idDaughterG_.size() > 0);
 }

XtoFFbarFilter::~XtoFFbarFilter ( )

inline

Definition at line 33 of file XtoFFbarFilter.h.

33 {}

Member Function Documentation

void XtoFFbarFilter::endJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 103 of file XtoFFbarFilter.cc.

References gather_cfg::cout, rejectedEvents_, totalEvents_, xSumCtau_, xSumPt_, xSumR_, and xTotal_.

                             {
   cout<<endl;
   cout<<"=== XtoFFbarFilter statistics of selected X->ffbar or Y->ggbar"<<endl;
   if (xTotal_ > 0) {
     cout<<"===   mean X & Y Pt = "<<xSumPt_/xTotal_<<" GeV and transverse decay length = "<<xSumR_/xTotal_<<" cm"<<endl;
     cout<<"===   mean c*tau = "<<xSumCtau_/xTotal_<<" cm"<<endl;
   } else {
     cout<<"===   WARNING: NONE FOUND !"<<endl;
   }
   cout<<"===   events rejected = "<<rejectedEvents_<<"/"<<totalEvents_<<endl;                     
 }

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

virtual

Implements edm::EDFilter.

Definition at line 19 of file XtoFFbarFilter.cc.

References foundXtoFFbar(), genParticles_, edm::Event::getByLabel(), idDaughterF_, idDaughterG_, idMotherX_, idMotherY_, j, rejectedEvents_, requireY_, src_, and totalEvents_.

                                                                           {
 
   iEvent.getByLabel(src_, genParticles_);
 
   totalEvents_++;
 
   unsigned int numX = 0;
   unsigned int numY = 0;
   unsigned int numXorY = 0;
 
   for (unsigned int j = 0; j < genParticles_->size(); j++) {
     GenParticleRef moth(genParticles_, j);
 
     // Is it X -> f fbar ?
     bool isXtoFFbar = this->foundXtoFFbar(moth, idMotherX_, idDaughterF_);
     if (isXtoFFbar) numX++; 
 
     if (!requireY_) {
 
       // Has X been found already ?
       if (numX >= 1) return true;
 
     } else {
 
       // Is it Y -> g gbar ?
       bool isYtoGGbar = this->foundXtoFFbar(moth, idMotherY_, idDaughterG_);
       if (isYtoGGbar) numY++;
       if (isXtoFFbar || isYtoGGbar) numXorY++;
 
       // Have X and Y been found already ?
       if (numX >= 1 && numY >= 1 && numXorY >= 2) return true;
     }
   }
 
   rejectedEvents_++;
   //  cout<<"REJECTED "<<totalEvents_<<endl;
   return false;
 }

bool XtoFFbarFilter::found	(	const std::vector< int > &	v,
		int	j
	)

inlineprivate

Definition at line 40 of file XtoFFbarFilter.h.

References spr::find(), and j.

Referenced by pf_01_customizeAll.SeqVisitor::__init__(), pf_01_customizeSimulation.SeqVisitor::__init__(), pf_01_customizeAll.SeqVisitor::enter(), pf_01_customizeSimulation.SeqVisitor::enter(), foundXtoFFbar(), pf_01_customizeAll.SeqVisitor::prepareSearch(), and pf_01_customizeSimulation.SeqVisitor::prepareSearch().

                                              {
     return std::find(v.begin(), v.end(), j) != v.end();
   }

bool XtoFFbarFilter::foundXtoFFbar	(	const reco::GenParticleRef &	moth,
		const std::vector< int > &	idMotherX,
		const std::vector< int > &	idDaughterF
	)

private

Definition at line 59 of file XtoFFbarFilter.cc.

References found(), i, j, rho, xSumCtau_, xSumPt_, xSumR_, and xTotal_.

Referenced by filter().

                                                   {
   // Check if given particle "moth" is X-->f fbar
   bool isXtoFFbar = false;
   int pdgIdMoth = moth->pdgId();
   double rho = -9.9e9;
 
   if (this->found(idMother, pdgIdMoth)) {
     bool foundF    = false;
     bool foundFbar = false;
     unsigned int nDau = moth->numberOfDaughters();
 
     for (unsigned int i = 0; i < nDau; i++) {
       GenParticleRef dau = moth->daughterRef(i);
       int pdgIdDau = dau->pdgId();
       if (this->found(idDaughter, -pdgIdDau))  foundFbar = true; 
       if (this->found(idDaughter,  pdgIdDau)) {foundF    = true; 
 
         // Just for statistics, get transverse decay length.
         // (To be really accurate, should do it w.r.t. P.V., but couldn't be bothered ...)
     // This is the normal case
         rho = dau->vertex().Rho();
         // Unfortunately, duplicate particles can appear in the event record. Handle this as follows:
         for (unsigned int j = 0; j < dau->numberOfDaughters(); j++) {
           GenParticleRef granddau = dau->daughterRef(j);        
           if (granddau->pdgId() == pdgIdDau) rho = granddau->vertex().Rho();
         }
       }
     }
     if (foundF && foundFbar) isXtoFFbar = true;
   }
 
   if (isXtoFFbar) {
     // Get statistics 
     xTotal_++;
     xSumPt_ += moth->pt();
     xSumR_  += rho;
     xSumCtau_ += rho*(moth->mass()/(moth->pt()+0.01)); // protection against unlikely case Pt = 0.
   }
 
   return isXtoFFbar;
 }

Member Data Documentation

edm::Handle<reco::GenParticleCollection> XtoFFbarFilter::genParticles_

private

Definition at line 60 of file XtoFFbarFilter.h.

Referenced by filter().

std::vector<int> XtoFFbarFilter::idDaughterF_

private

Definition at line 55 of file XtoFFbarFilter.h.

Referenced by filter().

std::vector<int> XtoFFbarFilter::idDaughterG_

private

Definition at line 57 of file XtoFFbarFilter.h.

Referenced by filter(), and XtoFFbarFilter().

std::vector<int> XtoFFbarFilter::idMotherX_

private

Definition at line 54 of file XtoFFbarFilter.h.

Referenced by filter().

std::vector<int> XtoFFbarFilter::idMotherY_

private

Definition at line 56 of file XtoFFbarFilter.h.

Referenced by filter(), and XtoFFbarFilter().

unsigned int XtoFFbarFilter::rejectedEvents_

private

Definition at line 68 of file XtoFFbarFilter.h.

Referenced by endJob(), and filter().

bool XtoFFbarFilter::requireY_

private

Definition at line 58 of file XtoFFbarFilter.h.

Referenced by filter(), and XtoFFbarFilter().

edm::InputTag XtoFFbarFilter::src_

private

Definition at line 53 of file XtoFFbarFilter.h.

Referenced by filter().

unsigned int XtoFFbarFilter::totalEvents_

private

Definition at line 67 of file XtoFFbarFilter.h.

Referenced by endJob(), and filter().

double XtoFFbarFilter::xSumCtau_

private

Definition at line 66 of file XtoFFbarFilter.h.

Referenced by endJob(), and foundXtoFFbar().

double XtoFFbarFilter::xSumPt_

private

Definition at line 64 of file XtoFFbarFilter.h.

Referenced by endJob(), and foundXtoFFbar().

double XtoFFbarFilter::xSumR_

private

Definition at line 65 of file XtoFFbarFilter.h.

Referenced by endJob(), and foundXtoFFbar().

unsigned int XtoFFbarFilter::xTotal_

private

Definition at line 63 of file XtoFFbarFilter.h.

Referenced by endJob(), and foundXtoFFbar().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation