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
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

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<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

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

Implementation: [Notes on implementation]

Definition at line 47 of file L1TMicroGMTInputProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 88 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 ( )

override

Definition at line 105 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

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

overrideprivate

Definition at line 325 of file L1TMicroGMTInputProducer.cc.

{
}

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

overrideprivate

Definition at line 313 of file L1TMicroGMTInputProducer.cc.

{
}

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

staticprivate

Definition at line 117 of file L1TMicroGMTInputProducer.cc.

References l1t::RegionalMuonCand::processor().

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

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

private

Definition at line 143 of file L1TMicroGMTInputProducer.cc.

References EnergyCorrector::c, and pileupDistInMC::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;
}

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

overrideprivate

Definition at line 331 of file L1TMicroGMTInputProducer.cc.

{
}

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

overrideprivate

Definition at line 319 of file L1TMicroGMTInputProducer.cc.

{
}

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

static

Definition at line 337 of file L1TMicroGMTInputProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

                                                                                     {
  //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);
}

void L1TMicroGMTInputProducer::openFile ( )

private

Definition at line 123 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";
    }
  }
}

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

overrideprivate

Definition at line 158 of file L1TMicroGMTInputProducer.cc.

References l1t::bmtf, l1t::emtf_neg, l1t::emtf_pos, stringResolutionProvider_cfi::et, f, mps_fire::i, createfilelist::int, m_currEvt, m_currType, m_endOfBx, m_filestream, eostools::move(), RPCpg::mu, l1t::omtf_neg, l1t::omtf_pos, edm::Event::put(), l1t::MuonCaloSum::setEtaBits(), l1t::MuonCaloSum::setEtBits(), l1t::RegionalMuonCand::setHwEta(), l1t::RegionalMuonCand::setHwPhi(), l1t::RegionalMuonCand::setHwPt(), l1t::RegionalMuonCand::setHwQual(), l1t::RegionalMuonCand::setHwSign(), l1t::RegionalMuonCand::setHwSignValid(), l1t::MuonCaloSum::setIndex(), l1t::MuonCaloSum::setPhiBits(), l1t::RegionalMuonCand::setTFIdentifiers(), createPayload::skip, AlCaHLTBitMon_QueryRunRegistry::string, and 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++;

}

void L1TMicroGMTInputProducer::skipHeader ( )

private

Definition at line 134 of file L1TMicroGMTInputProducer.cc.

References m_filestream, AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

Referenced by L1TMicroGMTInputProducer().

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

Member Data Documentation

int L1TMicroGMTInputProducer::m_currEvt

private

Definition at line 73 of file L1TMicroGMTInputProducer.cc.

Referenced by produce().

int L1TMicroGMTInputProducer::m_currType

private

Definition at line 72 of file L1TMicroGMTInputProducer.cc.

Referenced by produce().

bool L1TMicroGMTInputProducer::m_endOfBx

private

Definition at line 70 of file L1TMicroGMTInputProducer.cc.

Referenced by produce().

std::ifstream L1TMicroGMTInputProducer::m_filestream

private

Definition at line 69 of file L1TMicroGMTInputProducer.cc.

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

std::string L1TMicroGMTInputProducer::m_fname

private

Definition at line 68 of file L1TMicroGMTInputProducer.cc.

Referenced by L1TMicroGMTInputProducer(), and openFile().

bool L1TMicroGMTInputProducer::m_lastMuInBx

private

Definition at line 71 of file L1TMicroGMTInputProducer.cc.