CSCTFUnpacker Class Reference

#include <CSCTFUnpacker.h>

Inheritance diagram for CSCTFUnpacker:

Public Member Functions
	CSCTFUnpacker (const edm::ParameterSet &pset)
void	produce (edm::Event &e, const edm::EventSetup &c) override
	~CSCTFUnpacker (void) 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

Private Attributes
int	m_maxBX
int	m_minBX
CSCTriggerMapping *	mapping
edm::InputTag	producer
edm::EDGetTokenT< FEDRawDataCollection >	Raw_token
std::vector< int >	slot2sector
bool	swapME1strips
CSCTFEvent	tfEvent

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 22 of file CSCTFUnpacker.h.

Constructor & Destructor Documentation

CSCTFUnpacker()

CSCTFUnpacker::CSCTFUnpacker

(

const edm::ParameterSet &

pset

)

Definition at line 34 of file CSCTFUnpacker.cc.

                                                        : edm::stream::EDProducer<>(), mapping(nullptr) {
  LogDebug("CSCTFUnpacker|ctor") << "Started ...";
 
  // Edges of the time window, which LCTs are put into (unlike tracks, which are always centred around 0):
  m_minBX = pset.getParameter<int>("MinBX");  //3
  m_maxBX = pset.getParameter<int>("MaxBX");  //9
 
  // Swap: if(swapME1strips && me1b && !zplus) strip = 65 - strip; // 1-64 -> 64-1 :
  swapME1strips = pset.getParameter<bool>("swapME1strips");
 
  // Initialize slot<->sector assignment
  slot2sector = pset.getParameter<std::vector<int>>("slot2sector");
  LogDebug("CSCTFUnpacker|ctor") << "Verifying slot<->sector map from 'vint32 slot2sector'";
  for (int slot = 0; slot < 22; slot++)
    if (slot2sector[slot] < 0 || slot2sector[slot] > 12)
      throw cms::Exception("Invalid configuration") << "CSCTFUnpacker: sector index is set out of range (slot2sector["
                                                    << slot << "]=" << slot2sector[slot] << ", should be [0-12])";
  // Just for safety (in case of bad data):
  slot2sector.resize(32);
 
  // As we use standard CSC digi containers, we have to initialize mapping:
  std::string mappingFile = pset.getParameter<std::string>("mappingFile");
  if (mappingFile.length()) {
    LogDebug("CSCTFUnpacker|ctor") << "Define ``mapping'' only if you want to screw up real geometry";
    mapping = new CSCTriggerMappingFromFile(mappingFile);
  } else {
    LogDebug("CSCTFUnpacker|ctor") << "Generating default hw<->geometry mapping";
    class M : public CSCTriggerSimpleMapping {
      void fill(void) override {}
    };
    mapping = new M();
    for (int endcap = 1; endcap <= 2; endcap++)
      for (int station = 1; station <= 4; station++)
        for (int sector = 1; sector <= 6; sector++)
          for (int csc = 1; csc <= 9; csc++) {
            if (station == 1) {
              mapping->addRecord(endcap, station, sector, 1, csc, endcap, station, sector, 1, csc);
              mapping->addRecord(endcap, station, sector, 2, csc, endcap, station, sector, 2, csc);
            } else
              mapping->addRecord(endcap, station, sector, 0, csc, endcap, station, sector, 0, csc);
          }
  }
 
  producer = pset.getParameter<edm::InputTag>("producer");
 
  produces<CSCCorrelatedLCTDigiCollection>();
  produces<L1CSCTrackCollection>();
  produces<L1CSCStatusDigiCollection>();
  produces<CSCTriggerContainer<csctf::TrackStub>>("DT");
 
  Raw_token = consumes<FEDRawDataCollection>(producer);
}

References CSCTriggerMapping::addRecord(), makeMuonMisalignmentScenario::endcap, Exception, ntuplemaker::fill, LogDebug, m_maxBX, m_minBX, mapping, csctfunpacker_cfi::mappingFile, producer, muonDTDigis_cfi::pset, Raw_token, slot2sector, relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, and swapME1strips.

~CSCTFUnpacker()

CSCTFUnpacker::~CSCTFUnpacker ( void  )

override

Definition at line 87 of file CSCTFUnpacker.cc.

                              {
  if (mapping)
    delete mapping;
}

References mapping.

Member Function Documentation

produce()

void CSCTFUnpacker::produce ( edm::Event &  e, const edm::EventSetup &  c  )

override

Definition at line 92 of file CSCTFUnpacker.cc.

                                                               {
  // Get a handle to the FED data collection
  edm::Handle<FEDRawDataCollection> rawdata;
  e.getByToken(Raw_token, rawdata);
 
  // create the collection of CSC wire and strip digis as well as of DT stubs, which we receive from DTTF
  auto LCTProduct = std::make_unique<CSCCorrelatedLCTDigiCollection>();
  auto trackProduct = std::make_unique<L1CSCTrackCollection>();
  auto statusProduct = std::make_unique<L1CSCStatusDigiCollection>();
  auto dtProduct = std::make_unique<CSCTriggerContainer<csctf::TrackStub>>();
 
  for (int fedid = FEDNumbering::MINCSCTFFEDID; fedid <= FEDNumbering::MAXCSCTFFEDID; fedid++) {
    const FEDRawData& fedData = rawdata->FEDData(fedid);
    if (fedData.size() == 0)
      continue;
    //LogDebug("CSCTFUnpacker|produce");
    //if( monitor ) monitor->process((unsigned short*)fedData.data());
    unsigned int unpacking_status = tfEvent.unpack((unsigned short*)fedData.data(), fedData.size() / 2);
    if (unpacking_status == 0) {
      // There may be several SPs in event
      std::vector<const CSCSPEvent*> SPs = tfEvent.SPs_fast();
      // Cycle over all of them
      for (std::vector<const CSCSPEvent*>::const_iterator spItr = SPs.begin(); spItr != SPs.end(); spItr++) {
        const CSCSPEvent* sp = *spItr;
 
        L1CSCSPStatusDigi status;  
        status.sp_slot = sp->header().slot();
        status.l1a_bxn = sp->header().BXN();
        status.fmm_status = sp->header().status();
        status.track_cnt = sp->counters().track_counter();
        status.orbit_cnt = sp->counters().orbit_counter();
 
        // Finds central LCT BX
        // assumes window is odd number of bins
        int central_lct_bx = (m_maxBX + m_minBX) / 2;
 
        // Find central SP BX
        // assumes window is odd number of bins
        int central_sp_bx = int(sp->header().nTBINs() / 2);
 
        for (unsigned int tbin = 0; tbin < sp->header().nTBINs(); tbin++) {
          status.se |= sp->record(tbin).SEs();
          status.sm |= sp->record(tbin).SMs();
          status.bx |= sp->record(tbin).BXs();
          status.af |= sp->record(tbin).AFs();
          status.vp |= sp->record(tbin).VPs();
 
          for (unsigned int FPGA = 0; FPGA < 5; FPGA++)
            for (unsigned int MPClink = 0; MPClink < 3; ++MPClink) {
              std::vector<CSCSP_MEblock> lct = sp->record(tbin).LCT(FPGA, MPClink);
              if (lct.empty())
                continue;
 
              status.link_status[lct[0].spInput()] |= (1 << lct[0].receiver_status_frame1()) |
                                                      (1 << lct[0].receiver_status_frame2()) |
                                                      ((lct[0].aligment_fifo() ? 1 : 0) << 4);
              status.mpc_link_id |= (lct[0].link() << 2) | lct[0].mpc();
 
              int station = (FPGA ? FPGA : 1);
              int endcap = 0, sector = 0;
              if (slot2sector[sp->header().slot()]) {
                endcap = slot2sector[sp->header().slot()] / 7 + 1;
                sector = slot2sector[sp->header().slot()];
                if (sector > 6)
                  sector -= 6;
              } else {
                endcap = (sp->header().endcap() ? 1 : 2);
                sector = sp->header().sector();
              }
              int subsector = (FPGA > 1 ? 0 : FPGA + 1);
              int cscid = lct[0].csc();
 
              try {
                CSCDetId id = mapping->detId(endcap, station, sector, subsector, cscid, 0);
                // corrlcts now have no layer associated with them
                LCTProduct->insertDigi(id,
                                       CSCCorrelatedLCTDigi(0,
                                                            lct[0].vp(),
                                                            lct[0].quality(),
                                                            lct[0].wireGroup(),
                                                            (swapME1strips && cscid <= 3 && station == 1 &&
                                                                     endcap == 2 && lct[0].strip() < 65
                                                                 ? 65 - lct[0].strip()
                                                                 : lct[0].strip()),
                                                            lct[0].pattern(),
                                                            lct[0].l_r(),
                                                            (lct[0].tbin() + (central_lct_bx - central_sp_bx)),
                                                            lct[0].link(),
                                                            lct[0].BXN(),
                                                            0,
                                                            cscid));
              } catch (cms::Exception& e) {
                edm::LogInfo("CSCTFUnpacker|produce")
                    << e.what() << "Not adding digi to collection in event " << sp->header().L1A()
                    << " (endcap=" << endcap << ",station=" << station << ",sector=" << sector
                    << ",subsector=" << subsector << ",cscid=" << cscid << ",spSlot=" << sp->header().slot() << ")";
              }
            }
 
          std::vector<CSCSP_MBblock> mbStubs = sp->record(tbin).mbStubs();
          for (std::vector<CSCSP_MBblock>::const_iterator iter = mbStubs.begin(); iter != mbStubs.end(); iter++) {
            int endcap, sector;
            if (slot2sector[sp->header().slot()]) {
              endcap = slot2sector[sp->header().slot()] / 7 + 1;
              sector = slot2sector[sp->header().slot()];
              if (sector > 6)
                sector -= 6;
            } else {
              endcap = (sp->header().endcap() ? 1 : 2);
              sector = sp->header().sector();
            }
            const unsigned int csc2dt[6][2] = {{2, 3}, {4, 5}, {6, 7}, {8, 9}, {10, 11}, {12, 1}};
            DTChamberId id((endcap == 1 ? 2 : -2), 1, csc2dt[sector - 1][iter->id() - 1]);
            CSCCorrelatedLCTDigi base(0,
                                      iter->vq(),
                                      iter->quality(),
                                      iter->cal(),
                                      iter->flag(),
                                      iter->bc0(),
                                      iter->phi_bend(),
                                      tbin + (central_lct_bx - central_sp_bx),
                                      iter->id(),
                                      iter->bxn(),
                                      iter->timingError(),
                                      iter->BXN());
            csctf::TrackStub dtStub(base, id, iter->phi(), 0);
            dtProduct->push_back(dtStub);
          }
 
          std::vector<CSCSP_SPblock> tracks = sp->record(tbin).tracks();
          unsigned int trkNumber = 0;
          for (std::vector<CSCSP_SPblock>::const_iterator iter = tracks.begin(); iter != tracks.end();
               iter++, trkNumber++) {
            L1CSCTrack track;
            if (slot2sector[sp->header().slot()]) {
              track.first.m_endcap = slot2sector[sp->header().slot()] / 7 + 1;
              track.first.m_sector = slot2sector[sp->header().slot()];
              if (track.first.m_sector > 6)
                track.first.m_sector -= 6;
            } else {
              track.first.m_endcap = (sp->header().endcap() ? 1 : 2);
              track.first.m_sector = sp->header().sector();
            }
 
            track.first.m_lphi = iter->phi();
            track.first.m_ptAddress = iter->ptLUTaddress();
            track.first.m_fr = iter->f_r();
            track.first.m_ptAddress |= (iter->f_r() << 21);
 
            track.first.setStationIds(iter->ME1_id(), iter->ME2_id(), iter->ME3_id(), iter->ME4_id(), iter->MB_id());
            track.first.setTbins(
                iter->ME1_tbin(), iter->ME2_tbin(), iter->ME3_tbin(), iter->ME4_tbin(), iter->MB_tbin());
            track.first.setBx(iter->tbin() - central_sp_bx);
            track.first.setBits(iter->syncErr(), iter->bx0(), iter->bc0());
 
            track.first.setLocalPhi(iter->phi());
            track.first.setEtaPacked(iter->eta());
            track.first.setChargePacked(iter->charge());
 
            track.first.m_output_link = iter->id();
            if (track.first.m_output_link) {
              track.first.m_rank =
                  (iter->f_r() ? sp->record(tbin).ptSpy() & 0x7F : (sp->record(tbin).ptSpy() & 0x7F00) >> 8);
              track.first.setChargeValidPacked(
                  (iter->f_r() ? (sp->record(tbin).ptSpy() & 0x80) >> 8 : (sp->record(tbin).ptSpy() & 0x8000) >> 15));
            } else {
              track.first.m_rank = 0;
              track.first.setChargeValidPacked(0);
            }
            track.first.setFineHaloPacked(iter->halo());
 
            track.first.m_winner = iter->MS_id() & (1 << trkNumber);
 
            std::vector<CSCSP_MEblock> lcts = iter->LCTs();
            for (std::vector<CSCSP_MEblock>::const_iterator lct = lcts.begin(); lct != lcts.end(); lct++) {
              int station = (lct->spInput() > 6 ? (lct->spInput() - 1) / 3 : 1);
              int subsector = (lct->spInput() > 6 ? 0 : (lct->spInput() - 1) / 3 + 1);
              try {
                CSCDetId id =
                    mapping->detId(track.first.m_endcap, station, track.first.m_sector, subsector, lct->csc(), 0);
                track.second.insertDigi(id,
                                        CSCCorrelatedLCTDigi(0,
                                                             lct->vp(),
                                                             lct->quality(),
                                                             lct->wireGroup(),
                                                             (swapME1strips && lct->csc() <= 3 && station == 1 &&
                                                                      track.first.m_endcap == 2 && lct[0].strip() < 65
                                                                  ? 65 - lct[0].strip()
                                                                  : lct[0].strip()),
                                                             lct->pattern(),
                                                             lct->l_r(),
                                                             (lct->tbin() + (central_lct_bx - central_sp_bx)),
                                                             lct->link(),
                                                             lct->BXN(),
                                                             0,
                                                             lct->csc()));
              } catch (cms::Exception& e) {
                edm::LogInfo("CSCTFUnpacker|produce")
                    << e.what() << "Not adding track digi to collection in event" << sp->header().L1A()
                    << " (endcap=" << track.first.m_endcap << ",station=" << station
                    << ",sector=" << track.first.m_sector << ",subsector=" << subsector << ",cscid=" << lct->csc()
                    << ",spSlot=" << sp->header().slot() << ")";
              }
            }
 
            std::vector<CSCSP_MBblock> mbStubs = iter->dtStub();
            for (std::vector<CSCSP_MBblock>::const_iterator iter = mbStubs.begin(); iter != mbStubs.end(); iter++) {
              CSCDetId id = mapping->detId(track.first.m_endcap, 1, track.first.m_sector, iter->id(), 1, 0);
              track.second.insertDigi(id,
                                      CSCCorrelatedLCTDigi(iter->phi(),
                                                           iter->vq(),
                                                           iter->quality() + 100,
                                                           iter->cal(),
                                                           iter->flag(),
                                                           iter->bc0(),
                                                           iter->phi_bend(),
                                                           tbin + (central_lct_bx - central_sp_bx),
                                                           iter->id(),
                                                           iter->bxn(),
                                                           iter->timingError(),
                                                           iter->BXN()));
            }
 
            trackProduct->push_back(track);
          }
        }
        statusProduct->second.push_back(status);
      }
    } else {
      edm::LogError("CSCTFUnpacker|produce")
          << " problem of unpacking TF event: 0x" << std::hex << unpacking_status << std::dec << " code";
    }
 
    statusProduct->first = unpacking_status;
 
  }  //end of fed cycle
  e.put(std::move(dtProduct), "DT");
  e.put(std::move(LCTProduct));  // put processed lcts into the event.
  e.put(std::move(trackProduct));
  e.put(std::move(statusProduct));
}

References CSCSPRecord::AFs(), newFWLiteAna::base, CSCSPHeader::BXN(), CSCSPRecord::BXs(), CSCSPEvent::counters(), FEDRawData::data(), TauDecayModes::dec, CSCTriggerMapping::detId(), MillePedeFileConverter_cfg::e, CSCSPHeader::endcap(), makeMuonMisalignmentScenario::endcap, FEDRawDataCollection::FEDData(), PixelSLinkDataInputSource_cfi::fedid, CSCSPEvent::header(), triggerObjects_cff::id, createfilelist::int, CSCSPHeader::L1A(), CSCSPRecord::LCT(), MainPageGenerator::link, m_maxBX, m_minBX, mapping, FEDNumbering::MAXCSCTFFEDID, CSCSPRecord::mbStubs(), FEDNumbering::MINCSCTFFEDID, eostools::move(), CSCSPHeader::nTBINs(), CSCSPCounters::orbit_counter(), topSingleLeptonDQM_PU_cfi::pattern, CSCSPRecord::ptSpy(), quality, Raw_token, CSCSPEvent::record(), CSCSPHeader::sector(), CSCSPRecord::SEs(), FEDRawData::size(), CSCSPHeader::slot(), slot2sector, CSCSPRecord::SMs(), CSCTFEvent::SPs_fast(), relativeConstraints::station, mps_update::status, CSCSPHeader::status(), digitizers_cfi::strip, swapME1strips, tfEvent, HLT_FULL_cff::track, CSCSPCounters::track_counter(), CSCSPRecord::tracks(), tracks, CSCTFEvent::unpack(), and CSCSPRecord::VPs().

Member Data Documentation

m_maxBX

int CSCTFUnpacker::m_maxBX

private

Definition at line 24 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker(), and produce().

m_minBX

int CSCTFUnpacker::m_minBX

private

Definition at line 24 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker(), and produce().

mapping

CSCTriggerMapping* CSCTFUnpacker::mapping

private

Definition at line 27 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker(), produce(), progressbar.FormatLabel::update(), and ~CSCTFUnpacker().

producer

edm::InputTag CSCTFUnpacker::producer

private

Definition at line 36 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker().

Raw_token

edm::EDGetTokenT<FEDRawDataCollection> CSCTFUnpacker::Raw_token

private

Definition at line 38 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker(), and produce().

slot2sector

std::vector<int> CSCTFUnpacker::slot2sector

private

Definition at line 33 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker(), and produce().

swapME1strips

bool CSCTFUnpacker::swapME1strips

private

Definition at line 25 of file CSCTFUnpacker.h.

Referenced by CSCTFUnpacker(), and produce().

tfEvent

CSCTFEvent CSCTFUnpacker::tfEvent

private

Definition at line 29 of file CSCTFUnpacker.h.

Referenced by produce().