#include <MyEDFilter/STFilter/src/STFilter.cc>

Inheritance diagram for STFilter:

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

	~STFilter () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override

void	endJob () override

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

Private Attributes
unsigned int	accepted_events

edm::ParameterSet	conf_

int	DEBUGLVL

std::string	fOutputFileName

TH1D *	hbEta

TH1D *	hbEtaFiltered

TH1D *	hbPt

TH1D *	hbPtFiltered

edm::InputTag	hepMCProductTag_

TFile *	hOutputFile

unsigned int	input_events

bool	m_produceHistos

double	pTMax_

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

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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 &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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:

used for single top t-channel events generated with MadEvent and matched the "Karlsruhe way" filter on 2->2 process events, where the crucial candidate (the 2nd b quark) is not available until parton showering is done filter criterion: transverse momentum of 2nd b quark "pT < pTMax" -> event accepted! How-To: include STFilter.cfg in your .cfg, replace pTMax by the desired value, include module "STFilter" in outpath Implementation: <Notes on="" implementation>="">

Definition at line 13 of file STFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 of file STFilter.cc.

References accepted_events, conf_, DEBUGLVL, fOutputFileName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), input_events, m_produceHistos, pTMax_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                  :
   hepMCProductTag_(iConfig.getUntrackedParameter<edm::InputTag>("hepMCProductTag",edm::InputTag("generator","unsmeared"))) {
 
   pTMax_ = iConfig.getParameter<double>("pTMax");
   edm::LogInfo("SingleTopMatchingFilter")<<"+++ maximum pt of associated-b  pTMax = "<<pTMax_;
   DEBUGLVL = iConfig.getUntrackedParameter<int>("debuglvl", 0); // get debug level
   input_events=0; accepted_events=0; // counters
   m_produceHistos = iConfig.getParameter<bool>("produceHistos"); // produce histograms?
   fOutputFileName = iConfig.getUntrackedParameter<std::string>("histOutFile"); // get name of output file with histograms
   conf_ = iConfig;
 }

STFilter::~STFilter ( )

override

Definition at line 43 of file STFilter.cc.

43 {}

Member Function Documentation

void STFilter::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 115 of file STFilter.cc.

References conf_, fOutputFileName, edm::ParameterSet::getParameter(), hbEta, hbEtaFiltered, hbPt, hbPtFiltered, hOutputFile, and m_produceHistos.

                         {
   if(m_produceHistos){ // initialize histogram output file
     edm::ParameterSet Parameters;
     edm::LogInfo("SingleTopMatchingFilter)")<<"beginJob : creating histogram file: "<<fOutputFileName.c_str()<<std::endl;
     hOutputFile   = new TFile( fOutputFileName.c_str(), "RECREATE" ) ; hOutputFile->cd();
     // book histograms
     Parameters =  conf_.getParameter<edm::ParameterSet>("TH1bPt");
     hbPt = new TH1D( "bPt", "Pt of 2nd b quark", Parameters.getParameter<int32_t>("Nbinx"), Parameters.getParameter<double>("xmin"), Parameters.getParameter<double>("xmax"));
     Parameters =  conf_.getParameter<edm::ParameterSet>("TH1bEta");
     hbEta = new TH1D( "bEta", "Eta of 2nd b quark", Parameters.getParameter<int32_t>("Nbinx"), Parameters.getParameter<double>("xmin"), Parameters.getParameter<double>("xmax"));
     Parameters =  conf_.getParameter<edm::ParameterSet>("TH1bPtFiltered");
     hbPtFiltered = new TH1D( "bPtFiltered", "Pt of 2nd b quark filtered",  Parameters.getParameter<int32_t>("Nbinx"), Parameters.getParameter<double>("xmin"), Parameters.getParameter<double>("xmax"));
     Parameters =  conf_.getParameter<edm::ParameterSet>("TH1bEtaFiltered");
     hbEtaFiltered = new TH1D( "bEtaFiltered", "Eta of 2nd b quark filtered",  Parameters.getParameter<int32_t>("Nbinx"), Parameters.getParameter<double>("xmin"), Parameters.getParameter<double>("xmax"));
   }
 }

void STFilter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 132 of file STFilter.cc.

References accepted_events, gather_cfg::cout, relativeConstraints::error, dedxEstimators_cff::fraction, hbEta, hbEtaFiltered, hbPt, hbPtFiltered, hOutputFile, input_events, m_produceHistos, and mathSSE::sqrt().

                       {
  if(m_produceHistos){
    hOutputFile->cd();
    hbPt->Write(); hbEta->Write(); hbPtFiltered->Write(); hbEtaFiltered->Write();
    hOutputFile->Write() ; hOutputFile->Close() ;} // Write out histograms to file, then close it
    double fraction = (double) accepted_events / (double) input_events;
    double percent =  100. * fraction; double error =  100. * sqrt( fraction*(1-fraction) / (double) input_events );
    std::cout<<"STFilter ++ accepted_events/input_events = "<<accepted_events<<"/"<<input_events<<" = "<<fraction<<std::endl;
    std::cout<<"STFilter ++ efficiency  = "<<percent<<" % +/- "<<error<<" %"<<std::endl;
 }

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

overrideprivate

Definition at line 45 of file STFilter.cc.

References funct::abs(), accepted_events, edm::Event::getByLabel(), edm::HepMCProduct::GetEvent(), hbEta, hbEtaFiltered, hbPt, hbPtFiltered, hepMCProductTag_, mps_fire::i, input_events, funct::m, m_produceHistos, AlCaHLTBitMon_ParallelJobs::p, parents, and pTMax_.

                                                                    {
   using namespace edm;
   
   bool accEvt = false;
   bool lo = false;
 
   int secBcount = 0;
   double pTSecB = 100;
   double etaSecB = 100;
 
   ++input_events;
   
   Handle<HepMCProduct> evt;
   iEvent.getByLabel(hepMCProductTag_, evt);
   const HepMC::GenEvent * myEvt = evt->GetEvent();   // GET EVENT FROM HANDLE 
   
   bool bQuarksOpposite = false;
   for(HepMC::GenEvent::particle_const_iterator i = myEvt->particles_begin();
       (*i)->status() == 3; ++i) { // abort after status 3 particles
     // logic:
     // - in 2->2 matrix elements, the incoming (top production)
     //   b quark and the outgoing (top decay) b quark have same sign,
     //   so we flip the bQuarksOpposite flag twice -> false
     //   (opposite-sign b quark comes from the shower and has status 2)
     // - in 2->3 matrix elements, we have two outgoing b quarks with status
     //   3 and opposite signs -> true
     if ((*i)->pdg_id() == -5)
       bQuarksOpposite = !bQuarksOpposite;
   }
 
   // ---- 22 or 23? ----
   if (!bQuarksOpposite) // 22 
     lo = true;
   else
     accEvt = true;      // 23
   
   // ---- filter only 22 events ----
   if (lo){
     for (HepMC::GenEvent::particle_const_iterator p = myEvt->particles_begin(); p!=myEvt->particles_end(); ++p){
       // ---- look in shower for 2nd b quark ----
       if ((*p)->status() == 2 && abs((*p)->pdg_id()) == 5){  
     // ---- if b quark is found, loop over its parents ----
     for (HepMC::GenVertex::particle_iterator m = (*p)->production_vertex()->particles_begin(HepMC::parents); 
          m != (*p)->production_vertex()->particles_end(HepMC::parents); ++m){
       // ---- found 2ndb-candidate in shower ---- // ---- check mother of this candidate ----
       if(abs((*m)->barcode()) < 5){        
         if(secBcount == 1) break;
         secBcount ++;
         pTSecB = (*p)->operator HepMC::FourVector().perp();
         etaSecB = (*p)->operator HepMC::FourVector().eta();
       }
     }
       }
     }
     if(pTSecB < pTMax_) accEvt = true;
     // fill histos if requested
     if(m_produceHistos){
       hbPt->Fill(pTSecB);
       hbEta->Fill(etaSecB);
       if(accEvt){
         hbPtFiltered->Fill(pTSecB);
         hbEtaFiltered->Fill(etaSecB);
       }
     }
   }
   if(accEvt) ++accepted_events;
   return accEvt;
 }

Member Data Documentation

unsigned int STFilter::accepted_events

private

Definition at line 27 of file STFilter.h.

Referenced by endJob(), filter(), and STFilter().

edm::ParameterSet STFilter::conf_

private

Definition at line 36 of file STFilter.h.

Referenced by beginJob(), and STFilter().

int STFilter::DEBUGLVL

private

Definition at line 24 of file STFilter.h.

Referenced by STFilter().

std::string STFilter::fOutputFileName

private

Definition at line 33 of file STFilter.h.

Referenced by beginJob(), and STFilter().

TH1D* STFilter::hbEta

private

Definition at line 31 of file STFilter.h.

Referenced by beginJob(), endJob(), and filter().

TH1D* STFilter::hbEtaFiltered

private

Definition at line 31 of file STFilter.h.

Referenced by beginJob(), endJob(), and filter().

TH1D* STFilter::hbPt

private

Definition at line 30 of file STFilter.h.

Referenced by beginJob(), endJob(), and filter().

TH1D* STFilter::hbPtFiltered

private

Definition at line 30 of file STFilter.h.

Referenced by beginJob(), endJob(), and filter().

edm::InputTag STFilter::hepMCProductTag_

private

Definition at line 37 of file STFilter.h.

Referenced by filter().

TFile* STFilter::hOutputFile

private

Definition at line 34 of file STFilter.h.

Referenced by beginJob(), and endJob().

unsigned int STFilter::input_events

private

Definition at line 26 of file STFilter.h.

Referenced by endJob(), filter(), and STFilter().

bool STFilter::m_produceHistos

private

Definition at line 29 of file STFilter.h.

Referenced by beginJob(), endJob(), filter(), and STFilter().

double STFilter::pTMax_

private

Definition at line 22 of file STFilter.h.

Referenced by filter(), and STFilter().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation