#include <RecoJets/FFTJetProducers/plugins/FFTJetDijetFilter.cc>

Inheritance diagram for FFTJetDijetFilter:

Public Types
typedef fftjet::ProximityClusteringTree< fftjet::Peak, long >	ClusteringTree

typedef fftjet::SparseClusteringTree< fftjet::Peak, long >	SparseTree

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
	FFTJetDijetFilter (const edm::ParameterSet &)

	~FFTJetDijetFilter () 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 Types
typedef reco::PattRecoTree< float, reco::PattRecoPeak< float > >	StoredTree

Private Member Functions
void	beginJob () override

template<class Ptr >
void	checkConfig (const Ptr &ptr, const char *message) const

void	endJob () override

	FFTJetDijetFilter ()=delete

	FFTJetDijetFilter (const FFTJetDijetFilter &)=delete

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

FFTJetDijetFilter &	operator= (const FFTJetDijetFilter &)=delete

double	peakPt (const fftjet::Peak &peak) const

Private Attributes
ClusteringTree *	clusteringTree

double	completeEventScale

std::unique_ptr< fftjet::AbsDistanceCalculator< fftjet::Peak > >	distanceCalc

double	fixedScale

std::unique_ptr< std::vector< double > >	iniScales

bool	insertCompleteEvent

double	maxPeakEta

double	maxThirdJetFraction

double	min1to0PtRatio

double	minDeltaPhi

double	minPt0

double	minPt1

std::vector< unsigned >	nodes

unsigned long	nPassed

unsigned long	nRejected

std::vector< fftjet::Peak >	peaks

double	ptConversionFactor

SparseTree *	sparseTree

edm::InputTag	treeLabel

edm::EDGetTokenT< StoredTree >	treeToken

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: selects good dijet events looking only at the clustering tree

Implementation: [Notes on implementation]

Definition at line 51 of file FFTJetDijetFilter.cc.

Member Typedef Documentation

typedef fftjet::ProximityClusteringTree<fftjet::Peak,long> FFTJetDijetFilter::ClusteringTree

Definition at line 54 of file FFTJetDijetFilter.cc.

typedef fftjet::SparseClusteringTree<fftjet::Peak,long> FFTJetDijetFilter::SparseTree

Definition at line 55 of file FFTJetDijetFilter.cc.

typedef reco::PattRecoTree<float,reco::PattRecoPeak<float> > FFTJetDijetFilter::StoredTree

private

Definition at line 61 of file FFTJetDijetFilter.cc.

Constructor & Destructor Documentation

FFTJetDijetFilter::FFTJetDijetFilter ( const edm::ParameterSet & ps )

explicit

Definition at line 124 of file FFTJetDijetFilter.cc.

References checkConfig(), clusteringTree, distanceCalc, fftjetcms::fftjet_DistanceCalculator_parser(), fftjetcms::fftjet_ScaleSet_parser(), edm::ParameterSet::getParameter(), iniScales, jetUpdater_cfi::sort, sparseTree, fftjetdijetfilter_cfi::TreeDistanceCalculator, treeLabel, and treeToken.

     : init_param(edm::InputTag, treeLabel),
       init_param(double, ptConversionFactor),
       init_param(double, fixedScale),
       init_param(double, completeEventScale),
       init_param(double, min1to0PtRatio),
       init_param(double, minDeltaPhi),
       init_param(double, maxThirdJetFraction),
       init_param(double, minPt0),
       init_param(double, minPt1),
       init_param(double, maxPeakEta),
       init_param(bool, insertCompleteEvent),
       clusteringTree(nullptr),
       sparseTree(nullptr)
 {
     // Parse the set of scales
     iniScales = fftjet_ScaleSet_parser(
         ps.getParameter<edm::ParameterSet>("InitialScales"));
     checkConfig(iniScales, "invalid set of scales");
     std::sort(iniScales->begin(), iniScales->end(), std::greater<double>());
 
     // Parse the distance calculator
     const edm::ParameterSet& TreeDistanceCalculator(
         ps.getParameter<edm::ParameterSet>("TreeDistanceCalculator"));
     distanceCalc = fftjet_DistanceCalculator_parser(TreeDistanceCalculator);
     checkConfig(distanceCalc, "invalid tree distance calculator");
 
     // Create the clustering tree
     clusteringTree = new ClusteringTree(distanceCalc.get());
     sparseTree = new SparseTree();
 
     treeToken = consumes<StoredTree>(treeLabel);
 }

FFTJetDijetFilter::~FFTJetDijetFilter ( )

override

Definition at line 159 of file FFTJetDijetFilter.cc.

References clusteringTree, and sparseTree.

 {
     delete sparseTree;
     delete clusteringTree;
 }

FFTJetDijetFilter::FFTJetDijetFilter ( )

privatedelete

FFTJetDijetFilter::FFTJetDijetFilter ( const FFTJetDijetFilter & )

privatedelete

Member Function Documentation

void FFTJetDijetFilter::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 166 of file FFTJetDijetFilter.cc.

References nPassed, and nRejected.

 {
     nPassed = 0;
     nRejected = 0;
 }

template<class Ptr >

void FFTJetDijetFilter::checkConfig	(	const Ptr &	ptr,
		const char *	message
	)		const

inlineprivate

Definition at line 72 of file FFTJetDijetFilter.cc.

Referenced by FFTJetDijetFilter().

     {
         if (ptr.get() == nullptr)
             throw cms::Exception("FFTJetBadConfig") << message << std::endl;
     }

void FFTJetDijetFilter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 173 of file FFTJetDijetFilter.cc.

 {
 //     std::cout << "In FTJetDijetFilter::endJob: nPassed = " << nPassed
 //               << ", nRejected = " << nRejected << std::endl;
 }

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

overrideprivate

Definition at line 181 of file FFTJetDijetFilter.cc.

References funct::abs(), clusteringTree, completeEventScale, DEFINE_FWK_MODULE, reco::deltaPhi(), fftjetcms::densePeakTreeFromStorable(), PVValHelper::eta, fixedScale, edm::Event::getByToken(), mps_fire::i, iniScales, input, insertCompleteEvent, reco::PattRecoTree< ScaleType, Cluster >::isSparse(), maxPeakEta, maxThirdJetFraction, min1to0PtRatio, minDeltaPhi, minPt0, minPt1, nodes, nPassed, nRejected, peakPt(), peaks, phi, hiDetachedQuadStep_cff::pt1, jetUpdater_cfi::sort, fftjetcms::sparsePeakTreeFromStorable(), sparseTree, and treeToken.

 {
     edm::Handle<StoredTree> input;
     iEvent.getByToken(treeToken, input);
 
     // Convert the stored tree into a normal FFTJet clustering tree
     // and extract the set of peaks at the requested scale
     const double eventScale = insertCompleteEvent ? completeEventScale : 0.0;
     if (input->isSparse())
     {
         sparsePeakTreeFromStorable(*input, iniScales.get(),
                                    eventScale, sparseTree);
         sparseTree->sortNodes();
         fftjet::updateSplitMergeTimes(*sparseTree, sparseTree->minScale(),
                                       sparseTree->maxScale());
         const unsigned usedLevel = sparseTree->getLevel(fixedScale);
         sparseTree->getLevelNodes(usedLevel, &nodes);
         const unsigned numNodes = nodes.size();
         peaks.clear();
         peaks.reserve(numNodes);
         for (unsigned i=0; i<numNodes; ++i)
             peaks.push_back(sparseTree->uncheckedNode(nodes[i]).getCluster());
     }
     else
     {
         densePeakTreeFromStorable(*input, iniScales.get(),
                                   eventScale, clusteringTree);
         const unsigned usedLevel = clusteringTree->getLevel(fixedScale);
         double actualScale = 0.0;
         long dummyInfo;
         clusteringTree->getLevelData(usedLevel,&actualScale,&peaks,&dummyInfo);
     }
 
     // Get out if we don't have two clusters
     const unsigned nClusters = peaks.size();
     if (nClusters < 2)
     {
         ++nRejected;
         return false;
     }
     std::sort(peaks.begin(), peaks.end(), std::greater<fftjet::Peak>());
 
     // Calculate all quantities needed to make the pass/fail decision
     const double pt0 = peakPt(peaks[0]);
     const double pt1 = peakPt(peaks[1]);
     const double dphi = reco::deltaPhi(peaks[0].phi(), peaks[1].phi());
     const double ptratio = pt1/pt0;
     double thirdJetFraction = 0.0;
     if (nClusters > 2)
         thirdJetFraction = peakPt(peaks[2])/(pt0 + pt1);
 
     // Calculate the cut
     const bool pass = pt0 > minPt0 &&
                       pt1 > minPt1 &&
                       ptratio > min1to0PtRatio &&
                       std::abs(dphi) > minDeltaPhi &&
                       thirdJetFraction < maxThirdJetFraction &&
                       std::abs(peaks[0].eta()) < maxPeakEta &&
                       std::abs(peaks[1].eta()) < maxPeakEta;
     if (pass)
         ++nPassed;
     else
         ++nRejected;
     return pass;
 }

FFTJetDijetFilter& FFTJetDijetFilter::operator= ( const FFTJetDijetFilter & )

privatedelete

double FFTJetDijetFilter::peakPt ( const fftjet::Peak & peak ) const

inlineprivate

Definition at line 78 of file FFTJetDijetFilter.cc.

References fftjetdijetfilter_cfi::ptConversionFactor, and alignCSCRings::s.

Referenced by filter().

     {
         const double s = peak.scale();
         return ptConversionFactor*s*s*peak.magnitude();
     }

Member Data Documentation

ClusteringTree* FFTJetDijetFilter::clusteringTree

private

Definition at line 104 of file FFTJetDijetFilter.cc.

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

double FFTJetDijetFilter::completeEventScale

private

Definition at line 88 of file FFTJetDijetFilter.cc.

Referenced by filter().

std::unique_ptr<fftjet::AbsDistanceCalculator<fftjet::Peak> > FFTJetDijetFilter::distanceCalc

private

Definition at line 98 of file FFTJetDijetFilter.cc.

Referenced by FFTJetDijetFilter().

double FFTJetDijetFilter::fixedScale

private

Definition at line 87 of file FFTJetDijetFilter.cc.

Referenced by filter().

std::unique_ptr<std::vector<double> > FFTJetDijetFilter::iniScales

private

Definition at line 101 of file FFTJetDijetFilter.cc.

Referenced by FFTJetDijetFilter(), and filter().

bool FFTJetDijetFilter::insertCompleteEvent

private

Definition at line 95 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::maxPeakEta

private

Definition at line 94 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::maxThirdJetFraction

private

Definition at line 91 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::min1to0PtRatio

private

Definition at line 89 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::minDeltaPhi

private

Definition at line 90 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::minPt0

private

Definition at line 92 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::minPt1

private

Definition at line 93 of file FFTJetDijetFilter.cc.

Referenced by filter().

std::vector<unsigned> FFTJetDijetFilter::nodes

private

Definition at line 111 of file FFTJetDijetFilter.cc.

Referenced by filter().

unsigned long FFTJetDijetFilter::nPassed

private

Definition at line 114 of file FFTJetDijetFilter.cc.

Referenced by beginJob(), and filter().

unsigned long FFTJetDijetFilter::nRejected

private

Definition at line 115 of file FFTJetDijetFilter.cc.

Referenced by beginJob(), and filter().

std::vector<fftjet::Peak> FFTJetDijetFilter::peaks

private

Definition at line 108 of file FFTJetDijetFilter.cc.

Referenced by filter().

double FFTJetDijetFilter::ptConversionFactor

private

Definition at line 86 of file FFTJetDijetFilter.cc.

SparseTree* FFTJetDijetFilter::sparseTree

private

Definition at line 105 of file FFTJetDijetFilter.cc.

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

edm::InputTag FFTJetDijetFilter::treeLabel

private

Definition at line 85 of file FFTJetDijetFilter.cc.

Referenced by FFTJetDijetFilter().

edm::EDGetTokenT<StoredTree> FFTJetDijetFilter::treeToken

private

Definition at line 117 of file FFTJetDijetFilter.cc.

Referenced by FFTJetDijetFilter(), and filter().

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation