Inheritance diagram for L1MetPfProducer:

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

	~L1MetPfProducer () override

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

bool	wantsStreamLuminosityBlocks () const noexcept final

bool	wantsStreamRuns () const noexcept final

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

ModuleDescription const &	moduleDescription () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > 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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Types
typedef ap_fixed< 32, 16 >	input_t

typedef ap_ufixed< pt_t::width, 0 >	inv_t

typedef ap_int< 12 >	phi_t

typedef ap_fixed< 2 pt_t::width, 2 pt_t::iwidth, AP_RND, AP_SAT >	pt2_t

typedef ap_ufixed< 14, 12, AP_RND, AP_WRAP >	pt_t

typedef ap_fixed< pt_t::width+1, pt_t::iwidth+1, AP_RND, AP_SAT >	pxy_t

typedef ap_fixed< 32, 16 >	result_t

Private Member Functions
void	CalcMetHLS (const std::vector< float > &pt, const std::vector< float > &phi, reco::Candidate::PolarLorentzVector &metVector) const

void	CalcMlMet (const std::vector< float > &pt, const std::vector< float > &eta, const std::vector< float > &phi, const std::vector< float > &puppiWeight, const std::vector< int > &pdgId, const std::vector< int > &charge, reco::Candidate::PolarLorentzVector &metVector) const

int	EncodePdgId (int pdgId) const

void	PhiFromXY (pxy_t px, pxy_t py, phi_t &phi, bool debug=false) const

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

void	Project (pt_t pt, phi_t phi, pxy_t &pxy, bool isX, bool debug=false) const

Private Attributes
edm::EDGetTokenT< vector< l1t::PFCandidate > >	_l1PFToken

const phi_t	hwPi_ = round(M_PI / phiLSB_)

const phi_t	hwPiOverTwo_ = round(M_PI / (2 * phiLSB_))

int	maxCands_ = 128

std::shared_ptr< hls4mlEmulator::Model >	model

std::string	modelVersion_

bool	useMlModel_

Static Private Attributes
static constexpr int	dropBits_ = 2

static constexpr int	dropFactor_ = (1 << dropBits_)

static constexpr int	invTableBits_ = 10

static constexpr int	invTableSize_ = (1 << invTableBits_)

static constexpr float	maxPt_ = ((1 << pt_t::width) - 1) * ptLSB_

static constexpr int	numCatInputs_ = 2

static constexpr int	numContInputs_ = 4

static constexpr int	numInputs_ = numContInputs_ + numPxPyInputs_ + numCatInputs_

static constexpr int	numPxPyInputs_ = 2

static constexpr float	phiLSB_ = M_PI / 720

static constexpr float	ptLSB_ = 0.25

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
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)

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

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

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

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

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

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

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

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

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 21 of file L1MetPfProducer.cc.

Member Typedef Documentation

◆ input_t

typedef ap_fixed<32, 16> L1MetPfProducer::input_t

private

Definition at line 59 of file L1MetPfProducer.cc.

◆ inv_t

typedef ap_ufixed<pt_t::width, 0> L1MetPfProducer::inv_t

private

Definition at line 47 of file L1MetPfProducer.cc.

◆ phi_t

typedef ap_int<12> L1MetPfProducer::phi_t

private

Definition at line 35 of file L1MetPfProducer.cc.

◆ pt2_t

typedef ap_fixed<2 * pt_t::width, 2 * pt_t::iwidth, AP_RND, AP_SAT> L1MetPfProducer::pt2_t

private

Definition at line 41 of file L1MetPfProducer.cc.

◆ pt_t

typedef ap_ufixed<14, 12, AP_RND, AP_WRAP> L1MetPfProducer::pt_t

private

Definition at line 34 of file L1MetPfProducer.cc.

◆ pxy_t

typedef ap_fixed<pt_t::width + 1, pt_t::iwidth + 1, AP_RND, AP_SAT> L1MetPfProducer::pxy_t

private

Definition at line 40 of file L1MetPfProducer.cc.

◆ result_t

typedef ap_fixed<32, 16> L1MetPfProducer::result_t

private

Definition at line 60 of file L1MetPfProducer.cc.

Constructor & Destructor Documentation

◆ L1MetPfProducer()

L1MetPfProducer::L1MetPfProducer ( const edm::ParameterSet & cfg )

explicit

Definition at line 83 of file L1MetPfProducer.cc.

References web.browse_db::loader, model, modelVersion_, and useMlModel_.

     : _l1PFToken(consumes<std::vector<l1t::PFCandidate>>(cfg.getParameter<edm::InputTag>("L1PFObjects"))),
       maxCands_(cfg.getParameter<int>("maxCands")),
       modelVersion_(cfg.getParameter<std::string>("modelVersion")) {
   produces<std::vector<l1t::EtSum>>();
   useMlModel_ = (modelVersion_.length() > 0);
   if (useMlModel_) {
     hls4mlEmulator::ModelLoader loader(modelVersion_);
     model = loader.load_model();
   }
 }

◆ ~L1MetPfProducer()

L1MetPfProducer::~L1MetPfProducer ( )

override

Definition at line 340 of file L1MetPfProducer.cc.

340 {}

Member Function Documentation

◆ CalcMetHLS()

void L1MetPfProducer::CalcMetHLS	(	const std::vector< float > &	pt,
		const std::vector< float > &	phi,
		reco::Candidate::PolarLorentzVector &	metVector
	)		const

private

Definition at line 199 of file L1MetPfProducer.cc.

References ALCARECOEcalPhiSym_cff::float, mps_fire::i, maxPt_, SiStripPI::min, phi, lst_math::Phi_mpi_pi(), PhiFromXY(), phiLSB_, Project(), DiDispStaMuonMonitor_cfi::pt, mathSSE::sqrt(), and parallelization::uint.

Referenced by produce().

                                                                                    {
   pxy_t hw_px = 0;
   pxy_t hw_py = 0;
   pxy_t hw_sumx = 0;
   pxy_t hw_sumy = 0;
 
   for (uint i = 0; i < pt.size(); i++) {
     pt_t hw_pt = min(pt[i], maxPt_);
     phi_t hw_phi = float(TVector2::Phi_mpi_pi(phi[i]) / phiLSB_);
 
     Project(hw_pt, hw_phi, hw_px, true);
     Project(hw_pt, hw_phi, hw_py, false);
 
     hw_sumx = hw_sumx - hw_px;
     hw_sumy = hw_sumy - hw_py;
   }
 
   pt2_t hw_met = pt2_t(hw_sumx) * pt2_t(hw_sumx) + pt2_t(hw_sumy) * pt2_t(hw_sumy);
   hw_met = sqrt(int(hw_met));  // stand-in for HLS::sqrt
 
   phi_t hw_met_phi = 0;
   PhiFromXY(hw_sumx, hw_sumy, hw_met_phi);
 
   metVector.SetPt(hw_met.to_double());
   metVector.SetPhi(hw_met_phi.to_double() * phiLSB_);
   metVector.SetEta(0);
 }

◆ CalcMlMet()

void L1MetPfProducer::CalcMlMet	(	const std::vector< float > &	pt,
		const std::vector< float > &	eta,
		const std::vector< float > &	phi,
		const std::vector< float > &	puppiWeight,
		const std::vector< int > &	pdgId,
		const std::vector< int > &	charge,
		reco::Candidate::PolarLorentzVector &	metVector
	)		const

private

Definition at line 156 of file L1MetPfProducer.cc.

References funct::abs(), Matriplex::atan2(), ALCARECOTkAlJpsiMuMu_cff::charge, funct::cos(), EncodePdgId(), PVValHelper::eta, Matriplex::hypot(), mps_fire::i, input, maxCands_, model, numContInputs_, numInputs_, numPxPyInputs_, EgammaValidation_cff::pdgId, phi, DiDispStaMuonMonitor_cfi::pt, mps_fire::result, funct::sin(), and parallelization::uint.

Referenced by produce().

                                                                                   {
   const int inputSize = maxCands_ * numInputs_;
 
   input_t input[800];
   result_t result[2];
 
   // initialize with zeros (for padding)
   for (int i = 0; i < inputSize; i++) {
     input[i] = 0;
   }
 
   for (uint i = 0; i < pt.size(); i++) {
     // input_cont
     input[i * numContInputs_] = pt[i];
     input[i * numContInputs_ + 1] = eta[i];
     input[i * numContInputs_ + 2] = phi[i];
     input[i * numContInputs_ + 3] = puppiWeight[i];
     // input_pxpy
     input[(maxCands_ * numContInputs_) + (i * numPxPyInputs_)] = pt[i] * cos(phi[i]);
     input[(maxCands_ * numContInputs_) + (i * numPxPyInputs_) + 1] = pt[i] * sin(phi[i]);
     // input_cat0
     input[maxCands_ * (numContInputs_ + numPxPyInputs_) + i] = EncodePdgId(pdgId[i]);
     // input_cat1
     input[maxCands_ * (numContInputs_ + numPxPyInputs_ + 1) + i] = (abs(charge[i]) <= 1) ? (charge[i] + 2) : 0;
   }
 
   model->prepare_input(input);
   model->predict();
   model->read_result(result);
 
   double met_px = -result[0].to_double();
   double met_py = -result[1].to_double();
   metVector.SetPt(hypot(met_px, met_py));
   metVector.SetPhi(atan2(met_py, met_px));
   metVector.SetEta(0);
 }

◆ EncodePdgId()

int L1MetPfProducer::EncodePdgId ( int pdgId ) const

private

Definition at line 139 of file L1MetPfProducer.cc.

References funct::abs(), and EgammaValidation_cff::pdgId.

Referenced by CalcMlMet().

                                                 {
   switch (abs(pdgId)) {
     case 211:  // charged hadron (pion)
       return 1;
     case 130:  // neutral hadron (kaon)
       return 2;
     case 22:  // photon
       return 3;
     case 13:  // muon
       return 4;
     case 11:  // electron
       return 5;
     default:
       return 0;
   }
 }

◆ fillDescriptions()

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

static

Definition at line 95 of file L1MetPfProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                  {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("L1PFObjects", edm::InputTag("L1PFProducer", "l1pfCandidates"));
   desc.add<int>("maxCands", 128);
   desc.add<std::string>("modelVersion", "");
   descriptions.add("L1MetPfProducer", desc);
 }

◆ PhiFromXY()

void L1MetPfProducer::PhiFromXY	(	pxy_t	px,
		pxy_t	py,
		phi_t &	phi,
		bool	debug = `false`
	)		const

private

Definition at line 266 of file L1MetPfProducer.cc.

References a, Matriplex::atan2(), b, debug, dropFactor_, ALCARECOEcalPhiSym_cff::float, hwPi_, hwPiOverTwo_, invTableSize_, M_PI, maxPt_, phi, phiLSB_, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, x, and y.

Referenced by CalcMetHLS().

                                                                                 {
   if (px == 0 && py == 0) {
     phi = 0;
     return;
   }
   if (px == 0) {
     phi = py > 0 ? hwPiOverTwo_ : phi_t(-hwPiOverTwo_);
     return;
   }
   if (py == 0) {
     phi = px > 0 ? phi_t(0) : phi_t(-hwPi_);
     return;
   }
 
   // get q1 coordinates
   pt_t x = px > 0 ? pt_t(px) : pt_t(-px);  //px>=0 ? px : -px;
   pt_t y = py > 0 ? pt_t(py) : pt_t(-py);  //px>=0 ? px : -px;
   // transform so a<b
   pt_t a = x < y ? x : y;
   pt_t b = x < y ? y : x;
 
   if (b.to_double() > maxPt_ / dropFactor_)
     b = maxPt_ / dropFactor_;
   // map [0,max/4) to inv table size
   int index = round((b.to_double() / (maxPt_ / dropFactor_)) * invTableSize_);
   float bcheck = (float(index) / invTableSize_) * (maxPt_ / dropFactor_);
   inv_t inv_b = 1. / ((float(index) / invTableSize_) * (maxPt_ / dropFactor_));
 
   inv_t a_over_b = a * inv_b;
 
   if (debug) {
     printf("  a, b = %f, %f;   index, inv = %d, %f; ratio = %f \n",
            a.to_double(),
            b.to_double(),
            index,
            inv_b.to_double(),
            a_over_b.to_double());
     printf("    bcheck, 1/bc = %f, %f  -- %d  %f  %d  \n", bcheck, 1. / bcheck, invTableSize_, maxPt_, dropFactor_);
   }
 
   int atanTableBits_ = 7;
   int atanTableSize_ = (1 << atanTableBits_);
   index = round(a_over_b.to_double() * atanTableSize_);
   phi = atan(float(index) / atanTableSize_) / phiLSB_;
 
   if (debug) {
     printf("    atan index, phi = %d, %f (%f rad)  real atan(a/b)= %f  \n",
            index,
            phi.to_double(),
            phi.to_double() * (M_PI / hwPi_.to_double()),
            atan(a.to_double() / b.to_double()));
   }
 
   // rotate from (0,pi/4) to full quad1
   if (y > x)
     phi = hwPiOverTwo_ - phi;  //phi = pi/2 - phi
   // other quadrants
   if (px < 0 && py > 0)
     phi = hwPi_ - phi;  // Q2 phi = pi - phi
   if (px > 0 && py < 0)
     phi = -phi;  // Q4 phi = -phi
   if (px < 0 && py < 0)
     phi = -(hwPi_ - phi);  // Q3 composition of both
 
   if (debug) {
     printf("    phi hw, float, real = %f, %f    (%f rad from x,y = %f, %f) \n",
            phi.to_double(),
            phi.to_double() * (M_PI / hwPi_.to_double()),
            atan2(py.to_double(), px.to_double()),
            px.to_double(),
            py.to_double());
   }
 }

◆ produce()

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

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 103 of file L1MetPfProducer.cc.

References _l1PFToken, CalcMetHLS(), CalcMlMet(), ALCARECOTkAlJpsiMuMu_cff::charge, PVValHelper::eta, mps_fire::i, iEvent, createfilelist::int, L1Analysis::kTotalHt, maxCands_, susyDQM_cfi::metCollection, eostools::move(), EgammaValidation_cff::pdgId, phi, DiDispStaMuonMonitor_cfi::pt, and useMlModel_.

                                                                                               {
   edm::Handle<l1t::PFCandidateCollection> l1PFCandidates;
   iEvent.getByToken(_l1PFToken, l1PFCandidates);
 
   std::vector<float> pt;
   std::vector<float> eta;
   std::vector<float> phi;
   std::vector<float> puppiWeight;
   std::vector<int> pdgId;
   std::vector<int> charge;
 
   for (int i = 0; i < int(l1PFCandidates->size()) && (i < maxCands_ || maxCands_ < 0); i++) {
     const auto& l1PFCand = l1PFCandidates->at(i);
     pt.push_back(l1PFCand.pt());
     eta.push_back(l1PFCand.eta());
     phi.push_back(l1PFCand.phi());
     puppiWeight.push_back(l1PFCand.puppiWeight());
     pdgId.push_back(l1PFCand.pdgId());
     charge.push_back(l1PFCand.charge());
   }
 
   reco::Candidate::PolarLorentzVector metVector;
 
   if (useMlModel_) {
     CalcMlMet(pt, eta, phi, puppiWeight, pdgId, charge, metVector);
   } else {
     CalcMetHLS(pt, phi, metVector);
   }
 
   l1t::EtSum theMET(metVector, l1t::EtSum::EtSumType::kTotalHt, 0, 0, 0, 0);
 
   auto metCollection = std::make_unique<std::vector<l1t::EtSum>>(0);
   metCollection->push_back(theMET);
   iEvent.put(std::move(metCollection));
 }

◆ Project()

void L1MetPfProducer::Project	(	pt_t	pt,
		phi_t	phi,
		pxy_t &	pxy,
		bool	isX,
		bool	debug = `false`
	)		const

private

Definition at line 229 of file L1MetPfProducer.cc.

References funct::cos(), l1tPhase1JetProducer_cfi::cosPhi, hwPi_, hwPiOverTwo_, M_PI, phi, DiDispStaMuonMonitor_cfi::pt, objects.METAnalyzer::pxy, funct::sin(), and l1tPhase1JetProducer_cfi::sinPhi.

Referenced by CalcMetHLS().

                                                                                         {
   /*
       Convert pt and phi to px (py)
       1) Map phi to the first quadrant to reduce LUT size
       2) Lookup sin(phiQ1), where the result is in [0,maxPt]
       which is used to encode [0,1].
       3) Multiply pt by sin(phiQ1) to get px. Result will be px*maxPt, but
       wrapping multiplication is 'mod maxPt' so the correct value is returned.
       4) Check px=-|px|.
     */
 
   // set phi to first quadrant
   phi_t phiQ1 = (phi > 0) ? phi : phi_t(-phi);  // Q1/Q4
   if (phiQ1 >= hwPiOverTwo_)
     phiQ1 = hwPi_ - phiQ1;
 
   if (phiQ1 > hwPiOverTwo_) {
     edm::LogWarning("L1MetPfProducer") << "unexpected phi (high)";
     phiQ1 = hwPiOverTwo_;
   } else if (phiQ1 < 0) {
     edm::LogWarning("L1MetPfProducer") << "unexpected phi (low)";
     phiQ1 = 0;
   }
   if (isX) {
     typedef ap_ufixed<14, 12, AP_RND, AP_WRAP> pt_t;  // LSB is 0.25 and max is 4 TeV
     ap_ufixed<pt_t::width, 0> cosPhi = cos(phiQ1.to_double() / hwPiOverTwo_.to_double() * M_PI / 2);
     pxy = pt * cosPhi;
     if (phi > hwPiOverTwo_ || phi < -hwPiOverTwo_)
       pxy = -pxy;
   } else {
     ap_ufixed<pt_t::width, 0> sinPhi = sin(phiQ1.to_double() / hwPiOverTwo_.to_double() * M_PI / 2);
     pxy = pt * sinPhi;
     if (phi < 0)
       pxy = -pxy;
   }
 }

Member Data Documentation

◆ _l1PFToken

edm::EDGetTokenT<vector<l1t::PFCandidate> > L1MetPfProducer::_l1PFToken

private

Definition at line 29 of file L1MetPfProducer.cc.

Referenced by produce().

◆ dropBits_

constexpr int L1MetPfProducer::dropBits_ = 2

staticprivate

Definition at line 50 of file L1MetPfProducer.cc.

◆ dropFactor_

constexpr int L1MetPfProducer::dropFactor_ = (1 << dropBits_)

staticprivate

Definition at line 51 of file L1MetPfProducer.cc.

Referenced by PhiFromXY().

◆ hwPi_

const phi_t L1MetPfProducer::hwPi_ = round(M_PI / phiLSB_)

private

Definition at line 44 of file L1MetPfProducer.cc.

Referenced by PhiFromXY(), and Project().

◆ hwPiOverTwo_

const phi_t L1MetPfProducer::hwPiOverTwo_ = round(M_PI / (2 * phiLSB_))

private

Definition at line 45 of file L1MetPfProducer.cc.

Referenced by PhiFromXY(), and Project().

◆ invTableBits_

constexpr int L1MetPfProducer::invTableBits_ = 10

staticprivate

Definition at line 52 of file L1MetPfProducer.cc.

◆ invTableSize_

constexpr int L1MetPfProducer::invTableSize_ = (1 << invTableBits_)

staticprivate

Definition at line 53 of file L1MetPfProducer.cc.

Referenced by PhiFromXY().

◆ maxCands_

int L1MetPfProducer::maxCands_ = 128

private

Definition at line 31 of file L1MetPfProducer.cc.

Referenced by CalcMlMet(), and produce().

◆ maxPt_

constexpr float L1MetPfProducer::maxPt_ = ((1 << pt_t::width) - 1) * ptLSB_

staticprivate

Definition at line 43 of file L1MetPfProducer.cc.

Referenced by CalcMetHLS(), and PhiFromXY().

◆ model

std::shared_ptr<hls4mlEmulator::Model> L1MetPfProducer::model

private

Definition at line 57 of file L1MetPfProducer.cc.

Referenced by CalcMlMet(), L1MetPfProducer(), and querying.connection::search_everything().

◆ modelVersion_

std::string L1MetPfProducer::modelVersion_

private

Definition at line 58 of file L1MetPfProducer.cc.

Referenced by L1MetPfProducer().

◆ numCatInputs_

constexpr int L1MetPfProducer::numCatInputs_ = 2

staticprivate

Definition at line 63 of file L1MetPfProducer.cc.

◆ numContInputs_

constexpr int L1MetPfProducer::numContInputs_ = 4

staticprivate

Definition at line 61 of file L1MetPfProducer.cc.

Referenced by CalcMlMet().

◆ numInputs_

constexpr int L1MetPfProducer::numInputs_ = numContInputs_ + numPxPyInputs_ + numCatInputs_

staticprivate

Definition at line 64 of file L1MetPfProducer.cc.

Referenced by CalcMlMet().

◆ numPxPyInputs_

constexpr int L1MetPfProducer::numPxPyInputs_ = 2

staticprivate

Definition at line 62 of file L1MetPfProducer.cc.

Referenced by CalcMlMet().

◆ phiLSB_

constexpr float L1MetPfProducer::phiLSB_ = M_PI / 720

staticprivate

Definition at line 37 of file L1MetPfProducer.cc.

Referenced by CalcMetHLS(), and PhiFromXY().

◆ ptLSB_

constexpr float L1MetPfProducer::ptLSB_ = 0.25

staticprivate

Definition at line 36 of file L1MetPfProducer.cc.

◆ useMlModel_

bool L1MetPfProducer::useMlModel_

private

Definition at line 56 of file L1MetPfProducer.cc.

Referenced by L1MetPfProducer(), and produce().

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ input_t

◆ inv_t

◆ phi_t

◆ pt2_t

◆ pt_t

◆ pxy_t

◆ result_t

Constructor & Destructor Documentation

◆ L1MetPfProducer()

◆ ~L1MetPfProducer()

Member Function Documentation

◆ CalcMetHLS()

◆ CalcMlMet()

◆ EncodePdgId()

◆ fillDescriptions()

◆ PhiFromXY()

◆ produce()

◆ Project()

Member Data Documentation

◆ _l1PFToken

◆ dropBits_

◆ dropFactor_

◆ hwPi_

◆ hwPiOverTwo_

◆ invTableBits_

◆ invTableSize_

◆ maxCands_

◆ maxPt_

◆ model

◆ modelVersion_

◆ numCatInputs_

◆ numContInputs_

◆ numInputs_

◆ numPxPyInputs_

◆ phiLSB_

◆ ptLSB_

◆ useMlModel_