#include <PhysicsTools/FlavorHistoryFilter/src/FlavorHistoryFilter.cc>

Inheritance diagram for FlavorHistoryFilter:

Public Types
typedef reco::FlavorHistory::FLAVOR_T	flavor_type

typedef std::vector< int >	flavor_vector

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

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

	~FlavorHistoryFilter () 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	endJob () override

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

Additional Inherited Members
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:

This now filters events hierarchically. Previously this was done at the python configuration level, which was cumbersome for users to use.

Now, the hierarchy is:

Create prioritized paths to separate HF composition samples.

These are exclusive priorities, so sample "i" will not overlap with "i+1". Note that the "dr" values below correspond to the dr between the matched genjet, and the sister genjet.

1) W+bb with >= 2 jets from the ME (dr > 0.5) 2) W+b or W+bb with 1 jet from the ME 3) W+cc from the ME (dr > 0.5) 4) W+c or W+cc with 1 jet from the ME 5) W+bb with 1 jet from the parton shower (dr == 0.0) 6) W+cc with 1 jet from the parton shower (dr == 0.0)

These are the "trash bin" samples that we're throwing away:

7) W+bb with >= 2 partons but 1 jet from the ME (dr == 0.0) 8) W+cc with >= 2 partons but 1 jet from the ME (dr == 0.0) 9) W+bb with >= 2 partons but 2 jets from the PS (dr > 0.5) 10)W+cc with >= 2 partons but 2 jets from the PS (dr > 0.5)

And here is the true "light flavor" sample:

11) Veto of all the previous (W+ light jets)

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

Definition at line 76 of file FlavorHistoryFilter.h.

Member Typedef Documentation

typedef reco::FlavorHistory::FLAVOR_T FlavorHistoryFilter::flavor_type

Definition at line 78 of file FlavorHistoryFilter.h.

typedef std::vector<int> FlavorHistoryFilter::flavor_vector

Definition at line 79 of file FlavorHistoryFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 31 of file FlavorHistoryFilter.cc.

References b_me_, b_ps_, b_ps_comp_, bb_me_, bb_me_comp_, c_me_, c_ps_, c_ps_comp_, cc_me_, cc_me_comp_, dr_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), pathToSelect_, and verbose.

                                                                        :
   bsrcToken_           ( consumes<FlavorHistoryEvent>( iConfig.getParameter<edm::InputTag>("bsrc" ) ) ),
   csrcToken_           ( consumes<FlavorHistoryEvent>( iConfig.getParameter<edm::InputTag>("csrc" ) ) )
 {
   if ( iConfig.exists("pathToSelect") )
     pathToSelect_ =        iConfig.getParameter<int>    ("pathToSelect");
   else
     pathToSelect_ = -1;
 
   // This is the "interface" delta R with which to decide
   // where to take the event from
   dr_                =  iConfig.getParameter<double> ("dr" ) ;
   bool verbose        ( iConfig.getParameter<bool>   ("verbose") );
 
   // Set up the boundaries.
   // dr0 = 0.0
   // dr1 = set by user
   // dr2 = infinity
   double dr0 = 0.0;
   double dr1 = dr_;
   double dr2 = 99999.0;
 
 
   // These are the processes that can come from the matrix element calculation
   std::vector<int> me_ids;
   me_ids.push_back(2);   // flavor excitation
   me_ids.push_back(3);   // flavor creation
 
   // These are the processes that can come from the parton shower calculation
   std::vector<int> ps_ids;
   ps_ids.push_back(1);   // gluon splitting
 
 
   // To select bb->2 events from matrix element... Path 1
   bb_me_ = new FlavorHistorySelectorUtil( 5,
                       2,
                       me_ids,
                       dr1,
                       dr2,
                       verbose );
 
   // To select  b->1 events from matrix element... Path 2
   b_me_  = new FlavorHistorySelectorUtil( 5,
                       1,
                       me_ids,
                       dr0,
                       dr0,
                       verbose );
 
 
   // To select cc->2 events from matrix element... Path 3
   cc_me_ = new FlavorHistorySelectorUtil( 4,
                       2,
                       me_ids,
                       dr1,
                       dr2,
                       verbose );
 
   // To select  c->1 events from matrix element... Path 4
   c_me_  = new FlavorHistorySelectorUtil( 4,
                       1,
                       me_ids,
                       dr0,
                       dr0,
                       verbose );
 
   // To select bb->2 events from parton shower ... Path 5
   b_ps_  = new FlavorHistorySelectorUtil( 5,
                       1,
                       ps_ids,
                       dr0,
                       dr1,
                       verbose );
 
 
   // To select cc->2 events from parton shower ... Path 6
   c_ps_  = new FlavorHistorySelectorUtil( 4,
                       1,
                       ps_ids,
                       dr0,
                       dr1,
                       verbose );
 
   // To select bb->1 events from matrix element... Path 7
   bb_me_comp_ = new FlavorHistorySelectorUtil( 5,
                            2,
                            me_ids,
                            dr0,
                            dr1,
                            verbose );
 
   // To select cc->1 events from matrix element... Path 8
   cc_me_comp_ = new FlavorHistorySelectorUtil( 4,
                            2,
                            me_ids,
                            dr0,
                            dr1,
                            verbose );
 
   // To select bb->2 events from parton shower ... Path 9
   b_ps_comp_  = new FlavorHistorySelectorUtil( 5,
                            2,
                            ps_ids,
                            dr1,
                            dr2,
                            verbose );
 
   // To select cc->1 events from parton shower ... Path 10
   c_ps_comp_  = new FlavorHistorySelectorUtil( 4,
                            2,
                            ps_ids,
                            dr1,
                            dr2,
                            verbose );
 
   // The veto of all of these is               ... Path 11
 
 
   // This will write 1-11 (the path number), or 0 if error.
   produces<unsigned int>();
 }

FlavorHistoryFilter::~FlavorHistoryFilter ( )

override

Definition at line 154 of file FlavorHistoryFilter.cc.

References b_me_, b_ps_, b_ps_comp_, bb_me_, bb_me_comp_, c_me_, c_ps_, c_ps_comp_, cc_me_, and cc_me_comp_.

 {
 
   if ( bb_me_ ) delete bb_me_;
   if (  b_me_ ) delete b_me_;
   if ( cc_me_ ) delete cc_me_;
   if (  c_me_ ) delete c_me_;
   if (  b_ps_ ) delete b_ps_;
   if (  c_ps_ ) delete c_ps_;
 
   if ( bb_me_comp_ ) delete bb_me_comp_;
   if ( cc_me_comp_ ) delete cc_me_comp_;
   if (  b_ps_comp_ ) delete b_ps_comp_;
   if (  c_ps_comp_ ) delete c_ps_comp_;
 
 
 
 
 }

Member Function Documentation

void FlavorHistoryFilter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 249 of file FlavorHistoryFilter.cc.

References DEFINE_FWK_MODULE.

249 {

250 }

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

overrideprivate

Definition at line 177 of file FlavorHistoryFilter.cc.

References b_me_, b_ps_, b_ps_comp_, bb_me_, bb_me_comp_, bsrcToken_, c_me_, c_ps_, c_ps_comp_, cc_me_, cc_me_comp_, csrcToken_, reco::FlavorHistoryEvent::deltaR(), runTauDisplay::dr, reco::FlavorHistoryEvent::flavorSource(), edm::Event::getByToken(), reco::FlavorHistoryEvent::highestFlavor(), eostools::move(), reco::FlavorHistoryEvent::nb(), reco::FlavorHistoryEvent::nc(), pathToSelect_, edm::Event::put(), TauGenJetsDecayModeSelectorAllHadrons_cfi::select, reco::FlavorHistorySelectorUtil::select(), and corrVsCorr::selection.

 {
 
   // Get the flavor history
   Handle<FlavorHistoryEvent > bFlavorHistoryEvent;
   iEvent.getByToken(bsrcToken_,bFlavorHistoryEvent);
 
   Handle<FlavorHistoryEvent > cFlavorHistoryEvent;
   iEvent.getByToken(csrcToken_,cFlavorHistoryEvent);
 
   auto selection = std::make_unique<unsigned int>(0);
 
   // Get the number of matched b-jets in the event
   unsigned int nb = bFlavorHistoryEvent->nb();
   // Get the number of matched c-jets in the event
   unsigned int nc = cFlavorHistoryEvent->nc();
   // Get the two flavor sources. The highest takes precedence
   // over the rest.
   FlavorHistory::FLAVOR_T bFlavorSource = bFlavorHistoryEvent->flavorSource();
   FlavorHistory::FLAVOR_T cFlavorSource = cFlavorHistoryEvent->flavorSource();
   FlavorHistory::FLAVOR_T flavorSource = FlavorHistory::FLAVOR_NULL;
   // Get the highest flavor in the event
   unsigned int highestFlavor = 0;
   // Get the delta r between the two heavy flavor matched jets.
   double dr = -1;
 
   // Preference is in increasing priority:
   //  1: gluon splitting
   //  2: flavor excitation
   //  3: flavor creation (matrix element)
   //  4: flavor decay
   if ( bFlavorSource >= cFlavorSource ) {
     flavorSource = bFlavorHistoryEvent->flavorSource();
     highestFlavor = bFlavorHistoryEvent->highestFlavor();
     dr = bFlavorHistoryEvent->deltaR();
   }
   else {
     flavorSource = cFlavorHistoryEvent->flavorSource();
     highestFlavor = cFlavorHistoryEvent->highestFlavor();
     dr = cFlavorHistoryEvent->deltaR();
   }
 
 
   *selection = 0;
   // Now make hierarchical determination
   if      ( bb_me_     ->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 1;
   else if (  b_me_     ->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 2;
   else if ( cc_me_     ->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 3;
   else if (  c_me_     ->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 4;
   else if (  b_ps_     ->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 5;
   else if (  c_ps_     ->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 6;
   else if ( bb_me_comp_->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 7;
   else if ( cc_me_comp_->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 8;
   else if (  b_ps_comp_->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 9;
   else if (  c_ps_comp_->select( nb, nc, highestFlavor, flavorSource, dr ) ) *selection = 10;
   else *selection = 11;
 
   bool pass = false;
   if ( pathToSelect_ > 0 ) {
     pass = (*selection > 0 && *selection == static_cast<unsigned int>(pathToSelect_ ) );
   } else {
     pass = true;
   }
 
   iEvent.put(std::move(selection) );
 
   return pass;
 
 }

Member Data Documentation

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::b_me_

private

Definition at line 95 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::b_ps_

private

Definition at line 98 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::b_ps_comp_

private

Definition at line 102 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::bb_me_

private

Definition at line 94 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::bb_me_comp_

private

Definition at line 100 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

edm::EDGetTokenT<reco::FlavorHistoryEvent> FlavorHistoryFilter::bsrcToken_

private

Definition at line 89 of file FlavorHistoryFilter.h.

Referenced by filter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::c_me_

private

Definition at line 97 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::c_ps_

private

Definition at line 99 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::c_ps_comp_

private

Definition at line 103 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::cc_me_

private

Definition at line 96 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

reco::FlavorHistorySelectorUtil* FlavorHistoryFilter::cc_me_comp_

private

Definition at line 101 of file FlavorHistoryFilter.h.

Referenced by filter(), FlavorHistoryFilter(), and ~FlavorHistoryFilter().

edm::EDGetTokenT<reco::FlavorHistoryEvent> FlavorHistoryFilter::csrcToken_

private

Definition at line 90 of file FlavorHistoryFilter.h.

Referenced by filter().

double FlavorHistoryFilter::dr_

private

Definition at line 92 of file FlavorHistoryFilter.h.

Referenced by FlavorHistoryFilter().

int FlavorHistoryFilter::pathToSelect_

private

Definition at line 91 of file FlavorHistoryFilter.h.

Referenced by filter(), and FlavorHistoryFilter().

Private Attributes
reco::FlavorHistorySelectorUtil *	b_me_

reco::FlavorHistorySelectorUtil *	b_ps_

reco::FlavorHistorySelectorUtil *	b_ps_comp_

reco::FlavorHistorySelectorUtil *	bb_me_

reco::FlavorHistorySelectorUtil *	bb_me_comp_

edm::EDGetTokenT< reco::FlavorHistoryEvent >	bsrcToken_

reco::FlavorHistorySelectorUtil *	c_me_

reco::FlavorHistorySelectorUtil *	c_ps_

reco::FlavorHistorySelectorUtil *	c_ps_comp_

reco::FlavorHistorySelectorUtil *	cc_me_

reco::FlavorHistorySelectorUtil *	cc_me_comp_

edm::EDGetTokenT< reco::FlavorHistoryEvent >	csrcToken_

double	dr_

int	pathToSelect_

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation