ScGMTRawToDigi Class Reference

#include <ScGMTRawToDigi.h>

Inheritance diagram for ScGMTRawToDigi:

Public Member Functions
	ScGMTRawToDigi (const edm::ParameterSet &)
	~ScGMTRawToDigi () override
Public Member Functions inherited from edm::stream::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
const EDProducer &	operator= (const EDProducer &)=delete

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

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override
void	unpackOrbit (const unsigned char *buf, size_t len)

Private Attributes
bool	debug_ = false
int	nMuonsOrbit_
std::vector< std::vector< l1ScoutingRun3::Muon > >	orbitBuffer_
edm::EDGetToken	rawToken
bool	skipInterm_ = true
edm::InputTag	srcInputTag

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 24 of file ScGMTRawToDigi.h.

Constructor & Destructor Documentation

ScGMTRawToDigi()

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

explicit

Definition at line 3 of file ScGMTRawToDigi.cc.

References debug_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), nMuonsOrbit_, orbitBuffer_, rawToken, skipInterm_, and srcInputTag.

                                                             {
  using namespace edm;
  srcInputTag = iConfig.getParameter<InputTag>("srcInputTag");
  skipInterm_ = iConfig.getParameter<bool>("skipInterm");
  debug_ = iConfig.getUntrackedParameter<bool>("debug", false);

  // initialize orbit buffer for BX 1->3564;
  orbitBuffer_ = std::vector<std::vector<l1ScoutingRun3::Muon>>(3565);
  for (auto& bxVec : orbitBuffer_) {
    bxVec.reserve(8);
  }
  nMuonsOrbit_ = 0;

  produces<l1ScoutingRun3::MuonOrbitCollection>("Muon").setBranchAlias("MuonOrbitCollection");
  rawToken = consumes<SDSRawDataCollection>(srcInputTag);
}

~ScGMTRawToDigi()

ScGMTRawToDigi::~ScGMTRawToDigi ( )

override

Definition at line 20 of file ScGMTRawToDigi.cc.

{};

Member Function Documentation

fillDescriptions()

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

static

Definition at line 159 of file ScGMTRawToDigi.cc.

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

                                                                                {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("srcInputTag", edm::InputTag("rawDataCollector"));
  desc.add<bool>("skipInterm", true);
  desc.addUntracked<bool>("debug", false);
  descriptions.add("ScGMTRawToDigi", desc);
}

produce()

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

overrideprivate

Definition at line 22 of file ScGMTRawToDigi.cc.

References gather_cfg::cout, FEDRawData::data(), debug_, SDSRawDataCollection::FEDData(), SDSNumbering::GmtSDSID, iEvent, eostools::move(), nMuonsOrbit_, orbitBuffer_, rawToken, FEDRawData::size(), and unpackOrbit().

                                                                          {
  using namespace edm;

  Handle<SDSRawDataCollection> ScoutingRawDataCollection;
  iEvent.getByToken(rawToken, ScoutingRawDataCollection);

  const FEDRawData& sourceRawData = ScoutingRawDataCollection->FEDData(SDSNumbering::GmtSDSID);
  size_t orbitSize = sourceRawData.size();

  std::unique_ptr<l1ScoutingRun3::MuonOrbitCollection> unpackedMuons(new l1ScoutingRun3::MuonOrbitCollection);

  if ((sourceRawData.size() == 0) && debug_) {
    std::cout << "No raw data for GMT FED\n";
  }

  // unpack current orbit and store data into the orbitBufferr
  unpackOrbit(sourceRawData.data(), orbitSize);

  // fill orbit collection and clear the Bx buffer vector
  unpackedMuons->fillAndClear(orbitBuffer_, nMuonsOrbit_);

  // store collection in the event
  iEvent.put(std::move(unpackedMuons), "Muon");
}

unpackOrbit()

void ScGMTRawToDigi::unpackOrbit ( const unsigned char *  buf, size_t  len  )

private

Definition at line 47 of file ScGMTRawToDigi.cc.

References cms::cuda::assert(), visDQMUpload::buf, l1ScoutingRun3::ugmt::block::bx, nano_mu_digi_cff::bx, gather_cfg::cout, debug_, TauDecayModes::dec, PVValHelper::dxy, PVValHelper::eta, HLT_2024v10_cff::eta1, HLT_2024v10_cff::eta2, l1ScoutingRun3::ugmt::muon::extra, l1ScoutingRun3::ugmt::muon::f, RecoTauValidation_cfi::header, mps_fire::i, hcalRecHitTable_cff::ieta, hcaltpdigi_cfi::ieta1, hcaltpdigi_cfi::ieta2, hcalRecHitTable_cff::iphi, l1ScoutingRun3::ugmt::block::mu, HLT_2024v10_cff::muon, nMuonsOrbit_, l1ScoutingRun3::ugmt::block::orbit, orbitBuffer_, PVValHelper::phi, l1ScoutingRun3::printMuon(), DiDispStaMuonMonitor_cfi::pt, l1tGTMenu_lepSeeds_cff::qual, l1ScoutingRun3::ugmt::muon::s, skipInterm_, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by produce().

                                                                     {
  using namespace l1ScoutingRun3;

  // reset counters
  nMuonsOrbit_ = 0;

  size_t pos = 0;

  while (pos < len) {
    assert(pos + 4 <= len);

    // get BX header
    uint32_t header = *((uint32_t*)(buf + pos));
    pos += 4;
    // count mA and mB
    uint32_t mAcount = (header & header_masks::mAcount) >> header_shifts::mAcount;
    uint32_t mBcount = (header & header_masks::mBcount) >> header_shifts::mBcount;

    // declare block to read
    ugmt::block* bl = (ugmt::block*)(buf + pos);
    pos += 4 + 4 + (mAcount + mBcount) * 12;
    assert(pos <= len);

    uint32_t orbit = bl->orbit & 0x7FFFFFFF;
    uint32_t bx = bl->bx;

    if (debug_) {
      std::cout << "GMT Orbit " << orbit << ", BX -> " << bx << ", nMuons -> " << mAcount + mBcount << std::endl;
    }

    // Unpack muons for this BX
    orbitBuffer_[bx].reserve(mAcount + mBcount);

    for (unsigned int i = 0; i < mAcount + mBcount; i++) {
      uint32_t interm = (bl->mu[i].extra >> ugmt::shiftsMuon::interm) & ugmt::masksMuon::interm;
      if ((interm == 1) && (skipInterm_)) {
        if (debug_) {
          std::cout << " -> Excluding intermediate muon\n";
        }
        continue;
      }

      uint32_t index = (bl->mu[i].s >> ugmt::shiftsMuon::index) & ugmt::masksMuon::index;
      uint32_t ietaextu = (bl->mu[i].f >> ugmt::shiftsMuon::etaext) & ugmt::masksMuon::etaextv;
      int32_t ietaext;
      if (((bl->mu[i].f >> ugmt::shiftsMuon::etaext) & ugmt::masksMuon::etaexts) != 0) {
        ietaext = ietaextu -= 256;
      } else {
        ietaext = ietaextu;
      }

      // extract pt and quality and apply cut if required
      int32_t iptuncon = (bl->mu[i].s >> ugmt::shiftsMuon::ptuncon) & ugmt::masksMuon::ptuncon;
      int32_t ipt = (bl->mu[i].f >> ugmt::shiftsMuon::pt) & ugmt::masksMuon::pt;
      if ((ipt - 1) < 0) {
        continue;
      }
      uint32_t qual = (bl->mu[i].f >> ugmt::shiftsMuon::qual) & ugmt::masksMuon::qual;

      // extract integer value for extrapolated phi
      int32_t iphiext = ((bl->mu[i].f >> ugmt::shiftsMuon::phiext) & ugmt::masksMuon::phiext);

      // extract integer value for extrapolated phi
      int32_t idxy = ((bl->mu[i].s >> ugmt::shiftsMuon::dxy) & ugmt::masksMuon::dxy);

      // extract iso bits and charge
      uint32_t iso = (bl->mu[i].s >> ugmt::shiftsMuon::iso) & ugmt::masksMuon::iso;
      int32_t chrg = 0;
      if (((bl->mu[i].s >> ugmt::shiftsMuon::chrgv) & ugmt::masksMuon::chrgv) == 1)
        chrg = ((bl->mu[i].s >> ugmt::shiftsMuon::chrg) & ugmt::masksMuon::chrg) == 1 ? -1 : 1;

      // extract eta and phi at muon station
      int32_t iphi = (bl->mu[i].s >> ugmt::shiftsMuon::phi) & ugmt::masksMuon::phi;
      uint32_t ieta1 = (bl->mu[i].extra >> ugmt::shiftsMuon::eta1) & ugmt::masksMuon::eta;
      uint32_t ieta2 = (bl->mu[i].extra >> ugmt::shiftsMuon::eta2) & ugmt::masksMuon::eta;

      uint32_t ieta_u;
      int32_t ieta;
      // checking if raw eta should be taken from muon 1 or muon 2
      if ((bl->mu[i].extra & 0x1) == 0) {
        ieta_u = ieta1;
      } else {
        ieta_u = ieta2;
      }

      // two's complement
      if (ieta_u > 256) {
        ieta = ieta_u - 512;
      } else {
        ieta = ieta_u;
      }

      // increment muon counter
      nMuonsOrbit_++;

      l1ScoutingRun3::Muon muon(ipt, ieta, iphi, qual, chrg, chrg != 0, iso, index, ietaext, iphiext, iptuncon, idxy);

      orbitBuffer_[bx].push_back(muon);

      if (debug_) {
        std::cout << "--- Muon " << i << " ---\n";
        std::cout << "  Raw f:     0x" << std::hex << bl->mu[i].f << std::dec << "\n";
        std::cout << "  Raw s:     0x" << std::hex << bl->mu[i].s << std::dec << "\n";
        std::cout << "  Raw extra: 0x" << std::hex << bl->mu[i].extra << std::dec << "\n";
        printMuon(muon);
      }

    }  // end of bx

  }  // end orbit while loop
}

Member Data Documentation

debug_

bool ScGMTRawToDigi::debug_ = false

private

Definition at line 41 of file ScGMTRawToDigi.h.

Referenced by produce(), ScGMTRawToDigi(), and unpackOrbit().

nMuonsOrbit_

int ScGMTRawToDigi::nMuonsOrbit_

private

Definition at line 39 of file ScGMTRawToDigi.h.

Referenced by produce(), ScGMTRawToDigi(), and unpackOrbit().

orbitBuffer_

std::vector<std::vector<l1ScoutingRun3::Muon> > ScGMTRawToDigi::orbitBuffer_

private

Definition at line 38 of file ScGMTRawToDigi.h.

Referenced by produce(), ScGMTRawToDigi(), and unpackOrbit().

rawToken

edm::EDGetToken ScGMTRawToDigi::rawToken

private

Definition at line 44 of file ScGMTRawToDigi.h.

Referenced by produce(), and ScGMTRawToDigi().

skipInterm_

bool ScGMTRawToDigi::skipInterm_ = true

private

Definition at line 42 of file ScGMTRawToDigi.h.

Referenced by ScGMTRawToDigi(), and unpackOrbit().

srcInputTag

edm::InputTag ScGMTRawToDigi::srcInputTag

private

Definition at line 43 of file ScGMTRawToDigi.h.

Referenced by ScGMTRawToDigi().