#include <L1Trigger/L1CaloTrigger/plugin/Phase1L1TJetProducer.cc>

Inheritance diagram for Phase1L1TJetProducer:

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

	~Phase1L1TJetProducer () override

Public Member Functions inherited from edm::one::EDProducer<>
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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 &	operator= (EDConsumerBase &&)=default

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

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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
std::vector< reco::CaloJet >	_buildJetsFromSeeds (const std::vector< std::tuple< int, int >> &seeds) const

std::vector< reco::CaloJet >	_buildJetsFromSeedsWithPUSubtraction (const std::vector< std::tuple< int, int >> &seeds, bool killZeroPt) const

reco::CaloJet	buildJetFromSeed (const std::tuple< int, int > &seed) const

template<class Container >
void	fillCaloGrid (TH2F &caloGrid, const Container &triggerPrimitives)

std::vector< std::tuple< int, int > >	findSeeds (float seedThreshold) const
	Finds the seeds in the caloGrid, seeds are saved in a vector that contain the index in the TH2F of each seed. More...

float	getTowerEnergy (int iEta, int iPhi) const
	Get the energy of a certain tower while correctly handling phi periodicity in case of overflow. More...

void	produce (edm::Event &, const edm::EventSetup &) override

void	subtract9x9Pileup (reco::CaloJet &jet) const

Private Attributes
std::unique_ptr< TH2F >	caloGrid_

std::vector< double >	etaBinning_

edm::EDGetTokenT< edm::View< reco::Candidate > >	inputCollectionTag_

unsigned int	jetIEtaSize_

unsigned int	jetIPhiSize_

size_t	nBinsEta_

unsigned int	nBinsPhi_

std::string	outputCollectionName_

double	phiLow_

double	phiUp_

bool	puSubtraction_

double	seedPtThreshold_

bool	vetoZeroPt_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	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

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Detailed Description

Description: Produces jets with a phase-1 like sliding window algorithm using a collection of reco::Candidates in input.

INPUT PARAMETERS ***
etaBinning: vdouble with eta binning (allows non-homogeneous binning in eta)
nBinsPhi: uint32, number of bins in phi
phiLow: double, min phi (typically -pi)
phiUp: double, max phi (typically +pi)
jetIEtaSize: uint32, jet cluster size in ieta
jetIPhiSize: uint32, jet cluster size in iphi
seedPtThreshold: double, threshold of the seed tower
puSubtraction: bool, runs chunky doughnut pile-up subtraction, 9x9 jet only
outputCollectionName: string, tag for the output collection
vetoZeroPt: bool, controls whether jets with 0 pt should be save. 
It matters if PU is ON, as you can get negative or zero pt jets after it.
inputCollectionTag: inputtag, collection of reco::candidates used as input to the algo

Definition at line 56 of file Phase1L1TJetProducer.cc.

Constructor & Destructor Documentation

◆ Phase1L1TJetProducer()

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

explicit

Definition at line 102 of file Phase1L1TJetProducer.cc.

     :  // getting configuration settings
       inputCollectionTag_{
           consumes<edm::View<reco::Candidate>>(iConfig.getParameter<edm::InputTag>("inputCollectionTag"))},
       etaBinning_(iConfig.getParameter<std::vector<double>>("etaBinning")),
       nBinsEta_(etaBinning_.size() - 1),
       nBinsPhi_(iConfig.getParameter<unsigned int>("nBinsPhi")),
       phiLow_(iConfig.getParameter<double>("phiLow")),
       phiUp_(iConfig.getParameter<double>("phiUp")),
       jetIEtaSize_(iConfig.getParameter<unsigned int>("jetIEtaSize")),
       jetIPhiSize_(iConfig.getParameter<unsigned int>("jetIPhiSize")),
       seedPtThreshold_(iConfig.getParameter<double>("seedPtThreshold")),
       puSubtraction_(iConfig.getParameter<bool>("puSubtraction")),
       vetoZeroPt_(iConfig.getParameter<bool>("vetoZeroPt")),
       outputCollectionName_(iConfig.getParameter<std::string>("outputCollectionName")) {
   caloGrid_ =
       std::make_unique<TH2F>("caloGrid", "Calorimeter grid", nBinsEta_, etaBinning_.data(), nBinsPhi_, phiLow_, phiUp_);
   caloGrid_->GetXaxis()->SetTitle("#eta");
   caloGrid_->GetYaxis()->SetTitle("#phi");
   produces<std::vector<reco::CaloJet>>(outputCollectionName_).setBranchAlias(outputCollectionName_);
 }

References edm::ParameterSet::getParameter().

◆ ~Phase1L1TJetProducer()

Phase1L1TJetProducer::~Phase1L1TJetProducer ( )

override

Definition at line 124 of file Phase1L1TJetProducer.cc.

124 {}

Member Function Documentation

◆ _buildJetsFromSeeds()

std::vector< reco::CaloJet > Phase1L1TJetProducer::_buildJetsFromSeeds ( const std::vector< std::tuple< int, int >> & seeds ) const

private

Definition at line 306 of file Phase1L1TJetProducer.cc.

                                                       {
   // For each seed take a grid centered on the seed of the size specified by the user
   // Sum the pf in the grid, that will be the pt of the l1t jet. Eta and phi of the jet is taken from the seed.
   std::vector<reco::CaloJet> jets;
   for (const auto& seed : seeds) {
     reco::CaloJet jet = buildJetFromSeed(seed);
     jets.push_back(jet);
   }
   return jets;
 }

References buildJetFromSeed(), metsig::jet, singleTopDQM_cfi::jets, fileCollector::seed, and InitialStep_cff::seeds.

Referenced by produce().

◆ _buildJetsFromSeedsWithPUSubtraction()

std::vector< reco::CaloJet > Phase1L1TJetProducer::_buildJetsFromSeedsWithPUSubtraction	(	const std::vector< std::tuple< int, int >> &	seeds,
		bool	killZeroPt
	)		const

private

Definition at line 290 of file Phase1L1TJetProducer.cc.

                                                                        {
   // For each seed take a grid centered on the seed of the size specified by the user
   // Sum the pf in the grid, that will be the pt of the l1t jet. Eta and phi of the jet is taken from the seed.
   std::vector<reco::CaloJet> jets;
   for (const auto& seed : seeds) {
     reco::CaloJet jet = buildJetFromSeed(seed);
     subtract9x9Pileup(jet);
     //killing jets with 0 pt
     if ((vetoZeroPt_) && (jet.pt() <= 0))
       continue;
     jets.push_back(jet);
   }
   return jets;
 }

References buildJetFromSeed(), metsig::jet, singleTopDQM_cfi::jets, fileCollector::seed, InitialStep_cff::seeds, subtract9x9Pileup(), and vetoZeroPt_.

Referenced by produce().

◆ buildJetFromSeed()

reco::CaloJet Phase1L1TJetProducer::buildJetFromSeed ( const std::tuple< int, int > & seed ) const

private

Definition at line 265 of file Phase1L1TJetProducer.cc.

                                                                                        {
   int iEta = std::get<0>(seed);
   int iPhi = std::get<1>(seed);
  
   int etaHalfSize = (int)jetIEtaSize_ / 2;
   int phiHalfSize = (int)jetIPhiSize_ / 2;
  
   float ptSum = 0;
   // Scanning through the grid centered on the seed
   for (int etaIndex = -etaHalfSize; etaIndex <= etaHalfSize; etaIndex++) {
     for (int phiIndex = -phiHalfSize; phiIndex <= phiHalfSize; phiIndex++) {
       ptSum += getTowerEnergy(iEta + etaIndex, iPhi + phiIndex);
     }
   }
  
   // Creating a jet with eta phi centered on the seed and momentum equal to the sum of the pt of the components
   reco::Candidate::PolarLorentzVector ptVector;
   ptVector.SetPt(ptSum);
   ptVector.SetEta(caloGrid_->GetXaxis()->GetBinCenter(iEta));
   ptVector.SetPhi(caloGrid_->GetYaxis()->GetBinCenter(iPhi));
   reco::CaloJet jet;
   jet.setP4(ptVector);
   return jet;
 }

References caloGrid_, getTowerEnergy(), L1TowerCalibrationProducer_cfi::iEta, createfilelist::int, metsig::jet, jetIEtaSize_, jetIPhiSize_, and fileCollector::seed.

Referenced by _buildJetsFromSeeds(), and _buildJetsFromSeedsWithPUSubtraction().

◆ fillCaloGrid()

template<class Container >

void Phase1L1TJetProducer::fillCaloGrid	(	TH2F &	caloGrid,
		const Container &	triggerPrimitives
	)

private

Definition at line 319 of file Phase1L1TJetProducer.cc.

                                                                                           {
   //Filling the calo grid with the primitives
   for (const auto& primitiveIterator : triggerPrimitives) {
     caloGrid.Fill(static_cast<float>(primitiveIterator.eta()),
                   static_cast<float>(primitiveIterator.phi()),
                   static_cast<float>(primitiveIterator.pt()));
   }
 }

◆ fillDescriptions()

void Phase1L1TJetProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 328 of file Phase1L1TJetProducer.cc.

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("inputCollectionTag", edm::InputTag("l1pfCandidates", "Puppi"));
   desc.add<std::vector<double>>("etaBinning");
   desc.add<unsigned int>("nBinsPhi", 72);
   desc.add<double>("phiLow", -M_PI);
   desc.add<double>("phiUp", M_PI);
   desc.add<unsigned int>("jetIEtaSize", 7);
   desc.add<unsigned int>("jetIPhiSize", 7);
   desc.add<double>("seedPtThreshold", 5);
   desc.add<bool>("puSubtraction", false);
   desc.add<string>("outputCollectionName", "UncalibratedPhase1L1TJetFromPfCandidates");
   desc.add<bool>("vetoZeroPt", true);
   descriptions.add("Phase1L1TJetProducer", desc);
 }

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and M_PI.

◆ findSeeds()

std::vector< std::tuple< int, int > > Phase1L1TJetProducer::findSeeds ( float seedThreshold ) const

private

Finds the seeds in the caloGrid, seeds are saved in a vector that contain the index in the TH2F of each seed.

Definition at line 215 of file Phase1L1TJetProducer.cc.

                                                                                        {
   int nBinsX = caloGrid_->GetNbinsX();
   int nBinsY = caloGrid_->GetNbinsY();
  
   std::vector<std::tuple<int, int>> seeds;
  
   int etaHalfSize = (int)jetIEtaSize_ / 2;
   int phiHalfSize = (int)jetIPhiSize_ / 2;
  
   // for each point of the grid check if it is a local maximum
   // to do so I take a point, and look if is greater than the points around it (in the 9x9 neighborhood)
   // to prevent mutual exclusion, I check greater or equal for points above and right to the one I am considering (including the top-left point)
   // to prevent mutual exclusion, I check greater for points below and left to the one I am considering (including the bottom-right point)
  
   for (int iPhi = 1; iPhi <= nBinsY; iPhi++) {
     for (int iEta = 1; iEta <= nBinsX; iEta++) {
       float centralPt = caloGrid_->GetBinContent(iEta, iPhi);
       if (centralPt < seedThreshold)
         continue;
       bool isLocalMaximum = true;
  
       // Scanning through the grid centered on the seed
       for (int etaIndex = -etaHalfSize; etaIndex <= etaHalfSize; etaIndex++) {
         for (int phiIndex = -phiHalfSize; phiIndex <= phiHalfSize; phiIndex++) {
           if ((etaIndex == 0) && (phiIndex == 0))
             continue;
           if (phiIndex > 0) {
             if (phiIndex > -etaIndex) {
               isLocalMaximum = ((isLocalMaximum) && (centralPt > getTowerEnergy(iEta + etaIndex, iPhi + phiIndex)));
             } else {
               isLocalMaximum = ((isLocalMaximum) && (centralPt >= getTowerEnergy(iEta + etaIndex, iPhi + phiIndex)));
             }
           } else {
             if (phiIndex >= -etaIndex) {
               isLocalMaximum = ((isLocalMaximum) && (centralPt > getTowerEnergy(iEta + etaIndex, iPhi + phiIndex)));
             } else {
               isLocalMaximum = ((isLocalMaximum) && (centralPt >= getTowerEnergy(iEta + etaIndex, iPhi + phiIndex)));
             }
           }
         }
       }
       if (isLocalMaximum) {
         seeds.emplace_back(iEta, iPhi);
       }
     }
   }
  
   return seeds;
 }

References caloGrid_, getTowerEnergy(), L1TowerCalibrationProducer_cfi::iEta, createfilelist::int, jetIEtaSize_, jetIPhiSize_, InitialStep_cff::seeds, and UEAnalysisJets_cfi::seedThreshold.

Referenced by produce().

◆ getTowerEnergy()

float Phase1L1TJetProducer::getTowerEnergy	(	int	iEta,
		int	iPhi
	)		const

private

Get the energy of a certain tower while correctly handling phi periodicity in case of overflow.

Definition at line 126 of file Phase1L1TJetProducer.cc.

                                                                    {
   // We return the pt of a certain bin in the calo grid, taking account of the phi periodicity when overflowing (e.g. phi > phiSize), and returning 0 for the eta out of bounds
  
   int nBinsEta = caloGrid_->GetNbinsX();
   int nBinsPhi = caloGrid_->GetNbinsY();
   while (iPhi < 1) {
     iPhi += nBinsPhi;
   }
   while (iPhi > nBinsPhi) {
     iPhi -= nBinsPhi;
   }
   if (iEta < 1) {
     return 0;
   }
   if (iEta > nBinsEta) {
     return 0;
   }
   return caloGrid_->GetBinContent(iEta, iPhi);
 }

References caloGrid_, L1TowerCalibrationProducer_cfi::iEta, hgcalTowerMapProducer_cfi::nBinsEta, and Phase1L1TJetProducer_cfi::nBinsPhi.

Referenced by buildJetFromSeed(), findSeeds(), and subtract9x9Pileup().

◆ produce()

void Phase1L1TJetProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDProducerBase.

Definition at line 146 of file Phase1L1TJetProducer.cc.

                                                                                 {
   edm::Handle<edm::View<reco::Candidate>> inputCollectionHandle;
   iEvent.getByToken(inputCollectionTag_, inputCollectionHandle);
   // dumping the data
   caloGrid_->Reset();
   fillCaloGrid<>(*(caloGrid_), *inputCollectionHandle);
  
   const auto& seedsVector = findSeeds(seedPtThreshold_);  // seedPtThreshold = 6
   auto l1jetVector = puSubtraction_ ? _buildJetsFromSeedsWithPUSubtraction(seedsVector, vetoZeroPt_)
                                     : _buildJetsFromSeeds(seedsVector);
  
   // sort by pt
   std::sort(l1jetVector.begin(), l1jetVector.end(), [](const reco::CaloJet& jet1, const reco::CaloJet& jet2) {
     return jet1.pt() > jet2.pt();
   });
  
   auto l1jetVectorPtr = std::make_unique<std::vector<reco::CaloJet>>(l1jetVector);
   iEvent.put(std::move(l1jetVectorPtr), outputCollectionName_);
  
   return;
 }

References _buildJetsFromSeeds(), _buildJetsFromSeedsWithPUSubtraction(), caloGrid_, findSeeds(), iEvent, inputCollectionTag_, eostools::move(), outputCollectionName_, reco::LeafCandidate::pt(), puSubtraction_, seedPtThreshold_, and vetoZeroPt_.

◆ subtract9x9Pileup()

void Phase1L1TJetProducer::subtract9x9Pileup ( reco::CaloJet & jet ) const

private

Definition at line 168 of file Phase1L1TJetProducer.cc.

                                                                    {
   // these variables host the total pt in each sideband and the total pileup contribution
   float topBandPt = 0;
   float leftBandPt = 0;
   float rightBandPt = 0;
   float bottomBandPt = 0;
   float pileUpEnergy;
  
   // hold the jet's x-y (and z, as I have to use it, even if 2D) location in the histo
   int xCenter, yCenter, zCenter;
   // Retrieving histo-coords for seed
   caloGrid_->GetBinXYZ(caloGrid_->FindFixBin(jet.eta(), jet.phi()), xCenter, yCenter, zCenter);
  
   // Computing pileup
   for (int x = x_scroll_min; x <= x_scroll_max; x++) {
     for (int y = y_scroll_min; y < y_scroll_max; y++) {
       // top band, I go up 5 squares to reach the bottom of the top band
       // +x scrolls from left to right, +y scrolls up
       topBandPt += getTowerEnergy(xCenter + x, yCenter + (5 + y));
       // left band, I go left 5 squares (-5) to reach the bottom of the top band
       // +x scrolls from bottom to top, +y scrolls left
       leftBandPt += getTowerEnergy(xCenter - (5 + y), yCenter + x);
       // right band, I go right 5 squares (+5) to reach the bottom of the top band
       // +x scrolls from bottom to top, +y scrolls right
       rightBandPt += getTowerEnergy(xCenter + (5 + y), yCenter + x);
       // right band, I go right 5 squares (+5) to reach the bottom of the top band
       // +x scrolls from bottom to top, +y scrolls right
       bottomBandPt += getTowerEnergy(xCenter + x, yCenter - (5 + y));
     }
   }
   // adding bands and removing the maximum band (equivalent to adding the three minimum bands)
   pileUpEnergy = topBandPt + leftBandPt + rightBandPt + bottomBandPt -
                  std::max(topBandPt, std::max(leftBandPt, std::max(rightBandPt, bottomBandPt)));
  
   //preparing the new 4-momentum vector
   reco::Candidate::PolarLorentzVector ptVector;
   // removing pu contribution
   float ptAfterPUSubtraction = jet.pt() - pileUpEnergy;
   ptVector.SetPt((ptAfterPUSubtraction > 0) ? ptAfterPUSubtraction : 0);
   ptVector.SetEta(jet.eta());
   ptVector.SetPhi(jet.phi());
   //updating the jet
   jet.setP4(ptVector);
   jet.setPileup(pileUpEnergy);
   return;
 }

References caloGrid_, getTowerEnergy(), metsig::jet, SiStripPI::max, x, x_scroll_max, x_scroll_min, y, y_scroll_max, and y_scroll_min.

Referenced by _buildJetsFromSeedsWithPUSubtraction().

Member Data Documentation

◆ caloGrid_

std::unique_ptr<TH2F> Phase1L1TJetProducer::caloGrid_

private

Definition at line 86 of file Phase1L1TJetProducer.cc.

Referenced by buildJetFromSeed(), findSeeds(), getTowerEnergy(), produce(), and subtract9x9Pileup().

◆ etaBinning_

std::vector<double> Phase1L1TJetProducer::etaBinning_

private

Definition at line 88 of file Phase1L1TJetProducer.cc.

◆ inputCollectionTag_

edm::EDGetTokenT<edm::View<reco::Candidate> > Phase1L1TJetProducer::inputCollectionTag_

private

Definition at line 84 of file Phase1L1TJetProducer.cc.

Referenced by produce().

◆ jetIEtaSize_

unsigned int Phase1L1TJetProducer::jetIEtaSize_

private

Definition at line 93 of file Phase1L1TJetProducer.cc.

Referenced by buildJetFromSeed(), and findSeeds().

◆ jetIPhiSize_

unsigned int Phase1L1TJetProducer::jetIPhiSize_

private

Definition at line 94 of file Phase1L1TJetProducer.cc.

Referenced by buildJetFromSeed(), and findSeeds().

◆ nBinsEta_

size_t Phase1L1TJetProducer::nBinsEta_

private

Definition at line 89 of file Phase1L1TJetProducer.cc.

◆ nBinsPhi_

unsigned int Phase1L1TJetProducer::nBinsPhi_

private

Definition at line 90 of file Phase1L1TJetProducer.cc.

◆ outputCollectionName_

std::string Phase1L1TJetProducer::outputCollectionName_

private

Definition at line 99 of file Phase1L1TJetProducer.cc.

Referenced by produce().

◆ phiLow_

double Phase1L1TJetProducer::phiLow_

private

Definition at line 91 of file Phase1L1TJetProducer.cc.

◆ phiUp_

double Phase1L1TJetProducer::phiUp_

private

Definition at line 92 of file Phase1L1TJetProducer.cc.

◆ puSubtraction_

bool Phase1L1TJetProducer::puSubtraction_

private

Definition at line 96 of file Phase1L1TJetProducer.cc.

Referenced by produce().

◆ seedPtThreshold_

double Phase1L1TJetProducer::seedPtThreshold_

private

Definition at line 95 of file Phase1L1TJetProducer.cc.

Referenced by produce().

◆ vetoZeroPt_

bool Phase1L1TJetProducer::vetoZeroPt_

private

Definition at line 97 of file Phase1L1TJetProducer.cc.

Referenced by _buildJetsFromSeedsWithPUSubtraction(), and produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Phase1L1TJetProducer()

◆ ~Phase1L1TJetProducer()

Member Function Documentation

◆ _buildJetsFromSeeds()

◆ _buildJetsFromSeedsWithPUSubtraction()

◆ buildJetFromSeed()

◆ fillCaloGrid()

◆ fillDescriptions()

◆ findSeeds()

◆ getTowerEnergy()

◆ produce()

◆ subtract9x9Pileup()

Member Data Documentation

◆ caloGrid_

◆ etaBinning_

◆ inputCollectionTag_

◆ jetIEtaSize_

◆ jetIPhiSize_

◆ nBinsEta_

◆ nBinsPhi_

◆ outputCollectionName_

◆ phiLow_

◆ phiUp_

◆ puSubtraction_

◆ seedPtThreshold_

◆ vetoZeroPt_