#include <MuonRawDigiTranslator.h>

Static Public Member Functions
static int	calcHwEta (const uint32_t &, const unsigned, const unsigned)

static void	fillMuon (Muon &, uint32_t, uint32_t, int, unsigned int)

static void	fillMuon (Muon &, uint64_t, int, unsigned int)

static uint64_t	generate64bitDataWord (const Muon &)

static void	generatePackedDataWords (const Muon &, uint32_t &, uint32_t &)

Static Public Attributes
static const unsigned	absEtaAtVtxShift_ = 23

static const unsigned	absEtaMask_ = 0xFF

static const unsigned	absEtaShift_ = 21

static const unsigned	chargeShift_ = 2

static const unsigned	chargeValidShift_ = 3

static const unsigned	etaAtVtxSignShift_ = 31

static const unsigned	etaSignShift_ = 29

static const unsigned	isoMask_ = 0x3

static const unsigned	isoShift_ = 0

static const unsigned	phiAtVtxShift_ = 0

static const unsigned	phiMask_ = 0x3FF

static const unsigned	phiShift_ = 11

static const unsigned	ptMask_ = 0x1FF

static const unsigned	ptShift_ = 10

static const unsigned	qualMask_ = 0xF

static const unsigned	qualShift_ = 19

static const unsigned	tfMuonIndexMask_ = 0x7F

static const unsigned	tfMuonIndexShift_ = 4

Detailed Description

Definition at line 7 of file MuonRawDigiTranslator.h.

Member Function Documentation

int l1t::MuonRawDigiTranslator::calcHwEta	(	const uint32_t &	raw,
		const unsigned	absEtaShift,
		const unsigned	etaSignShift
	)

static

Definition at line 110 of file MuonRawDigiTranslator.cc.

References absEtaMask_.

Referenced by fillMuon().

 {
   // eta is coded as two's complement
   int abs_eta = (raw >> absEtaShift) & absEtaMask_;
   if ((raw >> etaSignShift) & 0x1) {
      return abs_eta - (1 << (etaSignShift - absEtaShift));
   } else {
      return abs_eta;
   }
 }

void l1t::MuonRawDigiTranslator::fillMuon	(	Muon &	mu,
		uint32_t	raw_data_00_31,
		uint32_t	raw_data_32_63,
		int	fed,
		unsigned int	fw
	)

static

Definition at line 5 of file MuonRawDigiTranslator.cc.

Referenced by fillMuon(), l1t::stage2::IntermediateMuonUnpacker::unpackBx(), and l1t::stage2::MuonUnpacker::unpackBx().

 {
   int hwPt = (raw_data_00_31 >> ptShift_) & ptMask_;
   if (hwPt > 0) {
     mu.setHwPt(hwPt);
     mu.setHwQual((raw_data_00_31 >> qualShift_) & qualMask_);
     mu.setHwIso((raw_data_32_63 >> isoShift_) & isoMask_);
     // charge is coded as -1^chargeBit
     mu.setHwCharge((raw_data_32_63 >> chargeShift_) & 0x1);
     mu.setHwChargeValid((raw_data_32_63 >> chargeValidShift_) & 0x1);
     mu.setTfMuonIndex((raw_data_32_63 >> tfMuonIndexShift_) & tfMuonIndexMask_);
 
     // The position of the eta and phi coordinates in the RAW data changed between the 2016 run and the 2017 run.
     // Eta and phi at the muon system are replaced by eta and phi at the vertex
     // Eta and phi at the muon system are moved to spare bits
     // The uGMT (FED 1402) or uGT (FED 1404) FW versions are used to determine the era.
     if ((fed == 1402 && fw < 0x4010000) || (fed == 1404 && fw < 0x10A6)) {
       // coordinates at the muon system are in 2016 where in 2017 eta and phi at the vertex are
       mu.setHwEta(calcHwEta(raw_data_00_31, absEtaAtVtxShift_, etaAtVtxSignShift_));
       mu.setHwPhi((raw_data_00_31 >> phiAtVtxShift_) & phiMask_);
 
       // set the coordiantes at vertex to be the same as the coordinates at the muon system
       mu.setHwEtaAtVtx(mu.hwEta());
       mu.setHwPhiAtVtx(mu.hwPhi());
       // deltas are 0
       mu.setHwDEtaExtra(0);
       mu.setHwDPhiExtra(0);
     } else {
       // coordinates at the muon system
       mu.setHwEta(calcHwEta(raw_data_32_63, absEtaShift_, etaSignShift_));
       mu.setHwPhi((raw_data_32_63 >> phiShift_) & phiMask_);
 
       // coordinates at the vertex
       mu.setHwEtaAtVtx(calcHwEta(raw_data_00_31, absEtaAtVtxShift_, etaAtVtxSignShift_));
       mu.setHwPhiAtVtx((raw_data_00_31 >> phiAtVtxShift_) & phiMask_);
       // deltas
       mu.setHwDEtaExtra(mu.hwEtaAtVtx() - mu.hwEta());
       int dPhi = mu.hwPhiAtVtx() - mu.hwPhi();
       if (mu.hwCharge() == 1 && dPhi > 0) {
         dPhi -= 576;
       } else if (mu.hwCharge() == 0 && dPhi < 0) {
         dPhi += 576;
       }
       mu.setHwDPhiExtra(dPhi);
     }
 
     math::PtEtaPhiMLorentzVector vec{(mu.hwPt()-1)*0.5, mu.hwEta()*0.010875, mu.hwPhi()*0.010908, 0.0};
     mu.setP4(vec);
     // generate a muon at the vertex to extract the physical eta and phi coordinates
     math::PtEtaPhiMLorentzVector vecAtVtx{(mu.hwPt()-1)*0.5, mu.hwEtaAtVtx()*0.010875, mu.hwPhiAtVtx()*0.010908, 0.0};
     Muon muAtVtx;
     muAtVtx.setP4(vecAtVtx);
     mu.setEtaAtVtx(muAtVtx.eta());
     mu.setPhiAtVtx(muAtVtx.phi());
     if (mu.hwChargeValid()) {
       mu.setCharge(1 - 2 * mu.hwCharge());
     } else {
       mu.setCharge(0);
     }
   }
 }

void l1t::MuonRawDigiTranslator::fillMuon	(	Muon &	mu,
		uint64_t	dataword,
		int	fed,
		unsigned int	fw
	)

static

Definition at line 68 of file MuonRawDigiTranslator.cc.

References fillMuon().

 {
   fillMuon(mu, (uint32_t)(dataword & 0xFFFFFFFF), (uint32_t)((dataword >> 32) & 0xFFFFFFFF), fed, fw);
 }

uint64_t l1t::MuonRawDigiTranslator::generate64bitDataWord ( const Muon & mu )

static

Definition at line 100 of file MuonRawDigiTranslator.cc.

References generatePackedDataWords().

 {
   uint32_t lsw;
   uint32_t msw;
 
   generatePackedDataWords(mu, lsw, msw);
   return (((uint64_t)msw) << 32) + lsw;
 }

void l1t::MuonRawDigiTranslator::generatePackedDataWords	(	const Muon &	mu,
		uint32_t &	raw_data_00_31,
		uint32_t &	raw_data_32_63
	)

static

Definition at line 74 of file MuonRawDigiTranslator.cc.

References absEtaAtVtxShift_, absEtaMask_, absEtaShift_, chargeShift_, chargeValidShift_, etaAtVtxSignShift_, etaSignShift_, l1t::Muon::hwCharge(), l1t::Muon::hwChargeValid(), l1t::L1Candidate::hwEta(), l1t::Muon::hwEtaAtVtx(), l1t::L1Candidate::hwIso(), l1t::L1Candidate::hwPhi(), l1t::Muon::hwPhiAtVtx(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), isoMask_, isoShift_, phiAtVtxShift_, phiMask_, phiShift_, ptMask_, ptShift_, qualMask_, qualShift_, l1t::Muon::tfMuonIndex(), tfMuonIndexMask_, and tfMuonIndexShift_.

Referenced by generate64bitDataWord(), and l1t::stage2::MuonPacker::pack().

 {
   int abs_eta = mu.hwEta();
   if (abs_eta < 0) {
     abs_eta += (1 << (etaSignShift_ - absEtaShift_));
   }
   int abs_eta_at_vtx = mu.hwEtaAtVtx();
   if (abs_eta_at_vtx < 0) {
     abs_eta_at_vtx += (1 << (etaAtVtxSignShift_ - absEtaAtVtxShift_));
   }
   raw_data_00_31 = (mu.hwPt() & ptMask_) << ptShift_
                  | (mu.hwQual() & qualMask_) << qualShift_
                  | (abs_eta_at_vtx & absEtaMask_) << absEtaAtVtxShift_
                  | (mu.hwEtaAtVtx() < 0) << etaAtVtxSignShift_
                  | (mu.hwPhiAtVtx() & phiMask_) << phiAtVtxShift_;
 
   raw_data_32_63 = mu.hwCharge() << chargeShift_
                  | mu.hwChargeValid() << chargeValidShift_
                  | (mu.tfMuonIndex() & tfMuonIndexMask_) << tfMuonIndexShift_
                  | (mu.hwIso() & isoMask_) << isoShift_
                  | (abs_eta & absEtaMask_) << absEtaShift_
                  | (mu.hwEta() < 0) << etaSignShift_
                  | (mu.hwPhi() & phiMask_) << phiShift_;
 }