L1TMicroGMTInputProducer Class Reference

#include <L1Trigger/L1TGlobalMuon/plugins/L1TMicroGMTInputProducer.cc>

Inheritance diagram for L1TMicroGMTInputProducer:

Public Member Functions
	L1TMicroGMTInputProducer (const edm::ParameterSet &)
	~L1TMicroGMTInputProducer () 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	beginLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override
void	beginRun (const edm::Run &, edm::EventSetup const &) override
int	convertToInt (std::string &bitstr) const
void	endLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override
void	endRun (const edm::Run &, edm::EventSetup const &) override
void	openFile ()
void	produce (edm::Event &, const edm::EventSetup &) override
void	skipHeader ()

Static Private Member Functions
static bool	cmpProc (const RegionalMuonCand &, const RegionalMuonCand &)

Private Attributes
int	m_currEvt
int	m_currType
bool	m_endOfBx
std::ifstream	m_filestream
std::string	m_fname
bool	m_lastMuInBx

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

Description: Takes txt-file input and produces barrel- / overlap- / forward TF muons

Implementation: [Notes on implementation]

Definition at line 46 of file L1TMicroGMTInputProducer.cc.

Constructor & Destructor Documentation

L1TMicroGMTInputProducer()

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

explicit

Definition at line 86 of file L1TMicroGMTInputProducer.cc.

References edm::ParameterSet::getParameter(), m_fname, openFile(), skipHeader(), and AlCaHLTBitMon_QueryRunRegistry::string.

    : m_endOfBx(false), m_currType(0), m_currEvt(0) {
  //register your products
  produces<RegionalMuonCandBxCollection>("BarrelTFMuons");
  produces<RegionalMuonCandBxCollection>("OverlapTFMuons");
  produces<RegionalMuonCandBxCollection>("ForwardTFMuons");
  produces<MuonCaloSumBxCollection>("TriggerTowerSums");

  //now do what ever other initialization is needed
  m_fname = iConfig.getParameter<std::string>("inputFileName");

  openFile();
  skipHeader();
}

~L1TMicroGMTInputProducer()

L1TMicroGMTInputProducer::~L1TMicroGMTInputProducer ( )

override

Definition at line 101 of file L1TMicroGMTInputProducer.cc.

References m_filestream.

                                                    {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
  m_filestream.close();
}

Member Function Documentation

beginLuminosityBlock()

void L1TMicroGMTInputProducer::beginLuminosityBlock ( const edm::LuminosityBlock &  , edm::EventSetup const &    )

overrideprivate

Definition at line 307 of file L1TMicroGMTInputProducer.cc.

{}

beginRun()

void L1TMicroGMTInputProducer::beginRun ( const edm::Run &  , edm::EventSetup const &    )

overrideprivate

Definition at line 301 of file L1TMicroGMTInputProducer.cc.

{}

cmpProc()

bool L1TMicroGMTInputProducer::cmpProc ( const RegionalMuonCand &  mu1, const RegionalMuonCand &  mu2  )

staticprivate

Definition at line 110 of file L1TMicroGMTInputProducer.cc.

References l1t::RegionalMuonCand::processor().

                                                                                               {
  return mu1.processor() < mu2.processor();
}

convertToInt()

int L1TMicroGMTInputProducer::convertToInt ( std::string &  bitstr ) const

private

Definition at line 130 of file L1TMicroGMTInputProducer.cc.

References HltBtagPostValidation_cff::c, and EgammaValidation_cff::num.

                                                                  {
  int num = 0;
  for (size_t cntr = 0; cntr < bitstr.size(); ++cntr) {
    char c = bitstr[cntr];
    num = (num << 1) |  // Shift the current set of bits to the left one bit
          (c - '0');    // Add in the current bit via a bitwise-or
  }
  return num;
}

endLuminosityBlock()

void L1TMicroGMTInputProducer::endLuminosityBlock ( const edm::LuminosityBlock &  , edm::EventSetup const &    )

overrideprivate

Definition at line 310 of file L1TMicroGMTInputProducer.cc.

{}

endRun()

void L1TMicroGMTInputProducer::endRun ( const edm::Run &  , edm::EventSetup const &    )

overrideprivate

Definition at line 304 of file L1TMicroGMTInputProducer.cc.

{}

fillDescriptions()

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

static

Definition at line 313 of file L1TMicroGMTInputProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                          {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

openFile()

void L1TMicroGMTInputProducer::openFile ( )

private

Definition at line 114 of file L1TMicroGMTInputProducer.cc.

References Exception, m_filestream, and m_fname.

Referenced by L1TMicroGMTInputProducer().

                                        {
  if (!m_filestream.is_open()) {
    m_filestream.open(m_fname.c_str());
    if (!m_filestream.good()) {
      cms::Exception("FileOpenError") << "Failed to open input file";
    }
  }
}

produce()

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

overrideprivate

Definition at line 141 of file L1TMicroGMTInputProducer.cc.

References l1t::bmtf, l1t::emtf_neg, l1t::emtf_pos, EgHLTOffHistBins_cfi::et, f, mps_fire::i, hcalRecHitTable_cff::ieta, iEvent, createfilelist::int, hcalRecHitTable_cff::iphi, m_currEvt, m_currType, m_endOfBx, m_filestream, eostools::move(), amptDefaultParameters_cff::mu, l1t::omtf_neg, l1t::omtf_pos, l1t::MuonCaloSum::setEtaBits(), l1t::MuonCaloSum::setEtBits(), l1t::MuonCaloSum::setIndex(), l1t::MuonCaloSum::setPhiBits(), optionsL1T::skip, AlCaHLTBitMon_QueryRunRegistry::string, and createJobs::tmp.

                                                                                    {
  using namespace edm;

  std::unique_ptr<RegionalMuonCandBxCollection> barrelMuons(new RegionalMuonCandBxCollection());
  std::unique_ptr<RegionalMuonCandBxCollection> overlapMuons(new RegionalMuonCandBxCollection());
  std::unique_ptr<RegionalMuonCandBxCollection> endcapMuons(new RegionalMuonCandBxCollection());
  std::unique_ptr<MuonCaloSumBxCollection> towerSums(new MuonCaloSumBxCollection());

  RegionalMuonCand mu;
  MuonCaloSum tSum;
  m_endOfBx = false;
  int caloCounter = 0;
  std::vector<int> bar{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
  std::vector<int> ovl_neg{0, 0, 0, 0, 0, 0};
  std::vector<int> ovl_pos{0, 0, 0, 0, 0, 0};
  std::vector<int> fwd_neg{0, 0, 0, 0, 0, 0};
  std::vector<int> fwd_pos{0, 0, 0, 0, 0, 0};
  while (!m_endOfBx && !m_filestream.eof()) {
    std::string lineID;
    m_filestream >> lineID;
    std::string restOfLine;

    if (lineID == "BAR" || lineID == "OVL-" || lineID == "FWD-" || lineID == "OVL+" || lineID == "FWD+") {
      int tmp;
      m_filestream >> tmp;  // cable no
      // if (lineID == "BAR") tmp += 12;
      // if (lineID == "OVL-") tmp = (tmp-6)+24;
      // if (lineID == "OVL+") tmp = tmp + 6;
      // if (lineID == "FWD-") tmp = (tmp-6)+30;

      // mu.setLink(tmp);
      m_filestream >> tmp;
      mu.setHwPt(tmp);

      m_filestream >> tmp;

      int globalPhi = int(tmp * 0.560856864654333f);  // correction from txt file producer!
      int globalWedgePhi = (globalPhi + 24) % 576;    // this sets CMS phi = 0 to -15 deg
      int globalSectorPhi = (globalPhi - 24);         // this sets CMS phi = 0 to +15 deg
      if (globalSectorPhi < 0) {
        globalSectorPhi += 576;
      }

      // int globalMuonPhi = int(tmp*0.560856864654333f); // make sure scale is correct
      bool skip = false;
      if (lineID == "BAR") {
        int processor = globalWedgePhi / 48 + 1;
        int localPhi = globalWedgePhi % 48;
        mu.setTFIdentifiers(processor, tftype::bmtf);
        mu.setHwPhi(localPhi);
        bar[processor - 1]++;
        if (bar[processor - 1] > 3)
          skip = true;
      }
      if (lineID == "OVL-") {
        int processor = globalSectorPhi / 96 + 1;
        int localPhi = globalSectorPhi % 96;
        mu.setTFIdentifiers(processor, tftype::omtf_neg);
        mu.setHwPhi(localPhi);
        ovl_neg[processor - 1]++;
        if (ovl_neg[processor - 1] > 3)
          skip = true;
      }
      if (lineID == "OVL+") {
        int processor = globalSectorPhi / 96 + 1;
        int localPhi = globalSectorPhi % 96;
        mu.setTFIdentifiers(processor, tftype::omtf_pos);
        mu.setHwPhi(localPhi);
        ovl_pos[processor - 1]++;
        if (ovl_pos[processor - 1] > 3)
          skip = true;
      }
      if (lineID == "FWD-") {
        int processor = globalSectorPhi / 96 + 1;
        int localPhi = globalSectorPhi % 96;
        mu.setTFIdentifiers(processor, tftype::emtf_neg);
        mu.setHwPhi(localPhi);
        fwd_neg[processor - 1]++;
        if (fwd_neg[processor - 1] > 3)
          skip = true;
      }
      if (lineID == "FWD+") {
        int processor = globalSectorPhi / 96 + 1;
        int localPhi = globalSectorPhi % 96;
        mu.setTFIdentifiers(processor, tftype::emtf_pos);
        mu.setHwPhi(localPhi);
        fwd_pos[processor - 1]++;
        if (fwd_pos[processor - 1] > 3)
          skip = true;
      }

      m_filestream >> tmp;
      tmp = int(tmp * 0.9090909090f);
      mu.setHwEta(tmp);

      m_filestream >> tmp;
      mu.setHwSign(tmp);

      m_filestream >> tmp;
      mu.setHwSignValid(tmp);

      m_filestream >> tmp;
      mu.setHwQual(tmp);

      if (lineID == "BAR")
        m_currType = 0;
      if (lineID == "OVL-")
        m_currType = 1;
      if (lineID == "OVL+")
        m_currType = 2;
      if (lineID == "FWD-")
        m_currType = 3;
      if (lineID == "FWD+")
        m_currType = 4;

      if (m_currType == 0 && !skip)
        barrelMuons->push_back(0, mu);
      if ((m_currType == 1 || m_currType == 2) && !skip)
        overlapMuons->push_back(0, mu);
      if ((m_currType == 3 || m_currType == 4) && !skip)
        endcapMuons->push_back(0, mu);
    }

    if (lineID == "EVT" && m_currEvt != 0) {
      m_endOfBx = true;
    } else if (lineID == "EVT") {
      m_currEvt++;
    }

    if (lineID == "CALO") {
      for (int i = 0; i < 28; ++i) {
        int ieta = i;  //caloCounter%28;
        int iphi = caloCounter;
        int et;

        m_filestream >> et;
        tSum.setEtBits(et);
        tSum.setEtaBits(ieta);
        tSum.setPhiBits(iphi);
        tSum.setIndex(caloCounter * 28 + i);
        towerSums->push_back(0, tSum);
      }
      caloCounter++;
    }
    getline(m_filestream, restOfLine);
    //std::cout << restOfLine;
  }

  // std::sort(barrelMuons->begin(0), barrelMuons->end(0), L1TMicroGMTInputProducer::cmpProc);
  // std::sort(overlapMuons->begin(0), overlapMuons->end(0), L1TMicroGMTInputProducer::cmpProc);
  // std::sort(endcapMuons->begin(0), endcapMuons->end(0), L1TMicroGMTInputProducer::cmpProc);

  iEvent.put(std::move(barrelMuons), "BarrelTFMuons");
  iEvent.put(std::move(overlapMuons), "OverlapTFMuons");
  iEvent.put(std::move(endcapMuons), "ForwardTFMuons");
  iEvent.put(std::move(towerSums), "TriggerTowerSums");
  m_currEvt++;
}

skipHeader()

void L1TMicroGMTInputProducer::skipHeader ( )

private

Definition at line 123 of file L1TMicroGMTInputProducer.cc.

References m_filestream, AlCaHLTBitMon_QueryRunRegistry::string, and createJobs::tmp.

Referenced by L1TMicroGMTInputProducer().

                                          {
  while (m_filestream.peek() == '#') {
    std::string tmp;
    getline(m_filestream, tmp);
  }
}

Member Data Documentation

m_currEvt

int L1TMicroGMTInputProducer::m_currEvt

private

Definition at line 72 of file L1TMicroGMTInputProducer.cc.

Referenced by produce().

m_currType

int L1TMicroGMTInputProducer::m_currType

private

Definition at line 71 of file L1TMicroGMTInputProducer.cc.

Referenced by produce().

m_endOfBx

bool L1TMicroGMTInputProducer::m_endOfBx

private

Definition at line 69 of file L1TMicroGMTInputProducer.cc.

Referenced by produce().

m_filestream

std::ifstream L1TMicroGMTInputProducer::m_filestream

private

Definition at line 68 of file L1TMicroGMTInputProducer.cc.

Referenced by openFile(), produce(), skipHeader(), and ~L1TMicroGMTInputProducer().

m_fname

std::string L1TMicroGMTInputProducer::m_fname

private

Definition at line 67 of file L1TMicroGMTInputProducer.cc.

Referenced by L1TMicroGMTInputProducer(), and openFile().

m_lastMuInBx

bool L1TMicroGMTInputProducer::m_lastMuInBx

private

Definition at line 70 of file L1TMicroGMTInputProducer.cc.