Classes
struct	EtMiss

Typedefs
typedef ap_int< kInternalPtWidth+kEtExtra >	Et_t

typedef ap_int< kInternalEtaWidth >	eta_t

typedef ap_uint< kMHTSize >	MHT_t

typedef ap_uint< kMHTPhiSize >	MHTphi_t

typedef ap_uint< 5 >	ntracks_t

typedef ap_int< kInternalPhiWidth >	phi_t

typedef ap_uint< kInternalPtWidth >	pt_t

Functions
EtMiss	cordicSqrt (Et_t x, Et_t y, int cordicSteps, std::vector< l1tmhtemu::MHTphi_t > atanLUT, std::vector< Et_t > magNormalisationLUT)

template<typename T >
T	digitizeSignedValue (double value, unsigned int nBits, double lsb)

std::vector< MHTphi_t >	generateaTanLUT (int cordicSteps)

std::vector< phi_t >	generateCosLUT (unsigned int size)

std::vector< Et_t >	generatemagNormalisationLUT (int cordicSteps)

Variables
const unsigned int	kEtExtra {10}

const unsigned int	kInternalEtaWidth {l1t::TkJetWord::TkJetBitWidths::kGlbEtaSize}

const unsigned int	kInternalPhiWidth {l1t::TkJetWord::TkJetBitWidths::kGlbPhiSize}

const unsigned int	kInternalPtWidth {l1t::TkJetWord::TkJetBitWidths::kPtSize}

const float	kMaxCosLUTPhi {M_PI}

const float	kMaxMHT {4096}

const float	kMaxMHTPhi {2 * M_PI}

const unsigned int	kMHTBins = 1 << kMHTSize

const unsigned int	kMHTPhiBins = 1 << kMHTPhiSize

const unsigned int	kMHTPhiSize {14}

const unsigned int	kMHTSize {15}

const unsigned int	kPhiBins = 1 << kInternalPhiWidth

const double	kStepEta {M_PI / (720)}

const double	kStepMHT = (l1tmhtemu::kMaxMHT / l1tmhtemu::kMHTBins)

const double	kStepMHTPhi = (l1tmhtemu::kMaxMHTPhi / l1tmhtemu::kMHTPhiBins)

const double	kStepPhi {M_PI / (720)}

const double	kStepPt {0.25}

Typedef Documentation

◆ Et_t

typedef ap_int<kInternalPtWidth + kEtExtra> l1tmhtemu::Et_t

Definition at line 40 of file L1TkHTMissEmulatorProducer.h.

◆ eta_t

typedef ap_int<kInternalEtaWidth> l1tmhtemu::eta_t

Definition at line 37 of file L1TkHTMissEmulatorProducer.h.

◆ MHT_t

typedef ap_uint<kMHTSize> l1tmhtemu::MHT_t

Definition at line 41 of file L1TkHTMissEmulatorProducer.h.

◆ MHTphi_t

typedef ap_uint<kMHTPhiSize> l1tmhtemu::MHTphi_t

Definition at line 42 of file L1TkHTMissEmulatorProducer.h.

◆ ntracks_t

typedef ap_uint<5> l1tmhtemu::ntracks_t

Definition at line 35 of file L1TkHTMissEmulatorProducer.h.

◆ phi_t

typedef ap_int<kInternalPhiWidth> l1tmhtemu::phi_t

Definition at line 38 of file L1TkHTMissEmulatorProducer.h.

◆ pt_t

typedef ap_uint<kInternalPtWidth> l1tmhtemu::pt_t

Definition at line 36 of file L1TkHTMissEmulatorProducer.h.

Function Documentation

◆ cordicSqrt()

EtMiss l1tmhtemu::cordicSqrt	(	Et_t	x,
		Et_t	y,
		int	cordicSteps,
		std::vector< l1tmhtemu::MHTphi_t >	atanLUT,
		std::vector< Et_t >	magNormalisationLUT
	)

inline

Definition at line 104 of file L1TkHTMissEmulatorProducer.h.

References l1tmhtemu::EtMiss::Et, dqmMemoryStats::float, kMHTBins, kMHTPhiBins, kStepMHT, kStepPhi, kStepPt, l1tmhtemu::EtMiss::Phi, Validation_hcalonly_cfi::sign, and x.

Referenced by L1TkHTMissEmulatorProducer::produce(), and L1TrackerEtMissEmulatorProducer::produce().

                                                                 {
     Et_t new_x = 0;
     Et_t new_y = 0;
 
     MHTphi_t phi = 0;
     MHTphi_t new_phi = 0;
     bool sign = false;
 
     EtMiss ret_etmiss;
 
     if (x >= 0 && y >= 0) {
       phi = 0;
       sign = true;
       //x = x;
       //y = y;
     } else if (x < 0 && y >= 0) {
       phi = kMHTPhiBins >> 1;
       sign = false;
       x = -x;
       //y = y;
     } else if (x < 0 && y < 0) {
       phi = kMHTPhiBins >> 1;
       sign = true;
       x = -x;
       y = -y;
     } else {
       phi = kMHTPhiBins;
       sign = false;
       //x = x;
       y = -y;
     }
 
     for (int step = 0; step < cordicSteps; step++) {
       if (y < 0) {
         new_x = x - (y >> step);
         new_y = y + (x >> step);
       } else {
         new_x = x + (y >> step);
         new_y = y - (x >> step);
       }
 
       if ((y < 0) == sign) {
         new_phi = phi - atanLUT[step];
       } else {
         new_phi = phi + atanLUT[step];
       }
 
       x = new_x;
       y = new_y;
       phi = new_phi;
     }
 
     float sqrtval = (float(x * magNormalisationLUT[cordicSteps - 1]) / float(kMHTBins)) * float(kStepPt * kStepPhi);
 
     ret_etmiss.Et = std::floor(sqrtval / l1tmhtemu::kStepMHT);
     ret_etmiss.Phi = phi;
 
     return ret_etmiss;
   }

◆ digitizeSignedValue()

template<typename T >

T l1tmhtemu::digitizeSignedValue	(	double	value,
		unsigned int	nBits,
		double	lsb
	)

Definition at line 59 of file L1TkHTMissEmulatorProducer.h.

References funct::abs(), and hgcalVFEProducer_cfi::lsb.

                                                                       {
     T digitized_value = std::floor(std::abs(value) / lsb);
     T digitized_maximum = (1 << (nBits - 1)) - 1;  // The remove 1 bit from nBits to account for the sign
     if (digitized_value > digitized_maximum)
       digitized_value = digitized_maximum;
     if (value < 0)
       digitized_value = (1 << nBits) - digitized_value;  // two's complement encoding
     return digitized_value;
   }

◆ generateaTanLUT()

std::vector<MHTphi_t> l1tmhtemu::generateaTanLUT ( int cordicSteps )

inline

Definition at line 80 of file L1TkHTMissEmulatorProducer.h.

References kMHTPhiBins, M_PI, and funct::pow().

Referenced by L1TkHTMissEmulatorProducer::L1TkHTMissEmulatorProducer().

                                                               {  // Fill atan LUT with integer values
     std::vector<MHTphi_t> atanLUT;
     atanLUT.reserve(cordicSteps);
     for (int cordicStep = 0; cordicStep < cordicSteps; cordicStep++) {
       atanLUT.push_back(MHTphi_t(round((kMHTPhiBins * atan(pow(2, -1 * cordicStep))) / (2 * M_PI))));
     }
     return atanLUT;
   }

◆ generateCosLUT()

std::vector<phi_t> l1tmhtemu::generateCosLUT ( unsigned int size )

inline

Definition at line 69 of file L1TkHTMissEmulatorProducer.h.

References funct::cos(), kInternalPhiWidth, kStepPhi, and findQualityFiles::size.

Referenced by L1TkHTMissEmulatorProducer::L1TkHTMissEmulatorProducer().

                                                             {  // Fill cosine LUT with integer values
     float phi = 0;
     std::vector<phi_t> cosLUT;
     for (unsigned int LUT_idx = 0; LUT_idx < size; LUT_idx++) {
       cosLUT.push_back(digitizeSignedValue<phi_t>(cos(phi), l1tmhtemu::kInternalPhiWidth, l1tmhtemu::kStepPhi));
       phi += l1tmhtemu::kStepPhi;
     }
     cosLUT.push_back((phi_t)(0));  //Prevent overflow in last bin
     return cosLUT;
   }

◆ generatemagNormalisationLUT()

std::vector<Et_t> l1tmhtemu::generatemagNormalisationLUT ( int cordicSteps )

inline

Definition at line 89 of file L1TkHTMissEmulatorProducer.h.

References kMHTBins, funct::pow(), and heppy_batch::val.

Referenced by L1TkHTMissEmulatorProducer::L1TkHTMissEmulatorProducer().

                                                                       {
     float val = 1.0;
     std::vector<Et_t> magNormalisationLUT;
     for (int cordicStep = 0; cordicStep < cordicSteps; cordicStep++) {
       val = val / (pow(1 + pow(4, -1 * cordicStep), 0.5));
       magNormalisationLUT.push_back(Et_t(round(kMHTBins * val)));
     }
     return magNormalisationLUT;
   }