MuonsGrabber.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     RPCTrigger
// Class  :     MuonsGrabber
//
// Implementation:
//     <Notes on implementation>
//
// Original Author:
//         Created:  Thu Sep 17 14:21:01 CEST 2009
//
 
// system include files
 
// user include files
#include "L1Trigger/RPCTrigger/interface/MuonsGrabber.h"
#include "Utilities/Xerces/interface/Xerces.h"
#include <xercesc/util/XMLString.hpp>
#include <xercesc/dom/DOM.hpp>
#include <xercesc/framework/LocalFileFormatTarget.hpp>
#include <sstream>
#include <algorithm>
 
XERCES_CPP_NAMESPACE_USE
 
class XStr {
public:
  XStr(const char* const toTranscode) { fUnicodeForm = XMLString::transcode(toTranscode); }
 
  ~XStr() { XMLString::release(&fUnicodeForm); }
 
  const XMLCh* unicodeForm() const { return fUnicodeForm; }
 
private:
  XMLCh* fUnicodeForm;
};
 
#define X(str) XStr(str).unicodeForm()
 
//
// constants, enums and typedefs
//
 
//
// static data member definitions
//
 
MuonsGrabber& MuonsGrabber::Instance() {
  static MuonsGrabber grabber;
  return grabber;
}
 
//
// constructors and destructor
//
MuonsGrabber::MuonsGrabber() {
  try {
    cms::concurrency::xercesInitialize();
  } catch (const XMLException& toCatch) {
    throw std::string("Error during Xerces-c Initialization: " +
                      std::string(XMLString::transcode(toCatch.getMessage())));
  }
 
  m_dom = DOMImplementationRegistry::getDOMImplementation(X("Core"));
  if (m_dom == nullptr)
    throw cms::Exception("RPCMuonsGrabber") << "Cannot get DOM" << std::endl;
 
  m_doc = m_dom->createDocument(nullptr,              // root element namespace URI.
                                X("rpctDataStream"),  // root element name
                                nullptr);             // document type object (DTD).
 
  m_rootElem = m_doc->getDocumentElement();
 
  m_currEvent = nullptr;
}
 
// MuonsGrabber::MuonsGrabber(const MuonsGrabber& rhs)
// {
//    // do actual copying here;
// }
 
MuonsGrabber::~MuonsGrabber() {
  // save xmlfile
  XMLCh tempStr[100];
  XMLString::transcode("LS", tempStr, 99);
  DOMImplementation* impl = DOMImplementationRegistry::getDOMImplementation(tempStr);
  DOMLSSerializer* theSerializer = ((DOMImplementationLS*)impl)->createLSSerializer();
  DOMConfiguration* dc = theSerializer->getDomConfig();
  dc->setParameter(XMLUni::fgDOMWRTFormatPrettyPrint, true);
  DOMLSOutput* outputDesc = ((DOMImplementationLS*)impl)->createLSOutput();
  outputDesc->setEncoding(X("UTF-8"));
 
  XMLFormatTarget* myFormTarget = new LocalFileFormatTarget(X("testpulses.xml"));
  outputDesc->setByteStream(myFormTarget);
  DOMNode* xmlstylesheet =
      m_doc->createProcessingInstruction(X("xml-stylesheet"), X("type=\"text/xsl\"href=\"default.xsl\""));
 
  m_doc->insertBefore(xmlstylesheet, m_rootElem);
  theSerializer->write(m_doc, outputDesc);
 
  delete theSerializer;
  delete myFormTarget;
  m_doc->release();
  cms::concurrency::xercesTerminate();
}
 
void MuonsGrabber::startNewEvent(int event, int bx) {
  //<event bx="0" num="0">
  m_currEvent = m_doc->createElement(X("event"));
  m_currEvent->setAttribute(X("num"), X(IntToString(event).c_str()));
  m_currEvent->setAttribute(X("bx"), X(IntToString(bx).c_str()));
  m_rootElem->appendChild(m_currEvent);
 
  m_currentEvent = event;
  m_currentBX = bx;
}
 
void MuonsGrabber::addMuon(RPCTBMuon& mu, int lvl, int region, int hs, int index) {
  if (mu.getPtCode() > 0)
    m_muons.push_back(RPCMuonExtraStruct(lvl, region, hs, index, mu));
}
 
void MuonsGrabber::writeDataForRelativeBX(int bx) {
  if (m_muons.empty())
    return;
 
  //<bxData num="11">
  DOMElement* currRelBx = m_doc->createElement(X("bxData"));
  currRelBx->setAttribute(X("num"), X(IntToString(bx).c_str()));
  m_currEvent->appendChild(currRelBx);
 
  //std::cout << "Writing out for relBx "   << bx << std::endl;
 
  //   <hs num="1" be="0">
  // levels
  // 0 - PAC
  // 1 - tbgs
  // 2 - tcgs
  // 3 - hs
  // 4 fs
 
  std::sort(m_muons.begin(), m_muons.end(), RPCMuonExtraStruct::lvlCompare);
  for (int tcNum = 0; tcNum <= 11; ++tcNum) {
    DOMElement* tc = nullptr;
    DOMElement* tcgs = nullptr;
    for (int tbNum = 0; tbNum <= 10; ++tbNum) {  // check actual range, probably till 9 total
      DOMElement* tb = nullptr;
      DOMElement* tbgs = nullptr;
      for (int PAC = 0; PAC <= 4; ++PAC) {  // same here
 
        DOMElement* pac = nullptr;
        // for (int segment = 0; segment <= 11; ++segment ) {
        std::vector<RPCMuonExtraStruct>::iterator it = m_muons.begin();
        while (it != m_muons.end()) {
          int muSegment = it->_mu.getLogSegment();
          //if
 
          int muTBno = m_trigConfig->getTBNum(it->_mu.getConeCrdnts());
          int muPACChipNo = m_trigConfig->getTowerNumOnTb(it->_mu.getConeCrdnts());
          int muTC = m_trigConfig->getTCNum(it->_mu.getConeCrdnts());
 
          if (!((int(it->_level) == 0 && tbNum == muTBno && muTC == tcNum && PAC == muPACChipNo) ||
                (int(it->_level) == 1 && tbNum == muTBno && muTC == tcNum) || (int(it->_level) == 2 && muTC == tcNum)))
 
          {
            ++it;
            continue;
          }
          // std::cout << int(it->_level) << int(it->_region) << int(it->_hsHalf)
          //         << " " << int(it->_index)
          //       << " " << it->_mu.printDebugInfo(2) << std::endl;
 
          if (tc == nullptr) {
            tc = m_doc->createElement(X("tc"));
            currRelBx->appendChild(tc);
            tc->setAttribute(X("num"), X(IntToString(tcNum).c_str()));
 
            tcgs = m_doc->createElement(X("tcgs"));
            tc->appendChild(tcgs);
          }
          if (tb == nullptr && int(it->_level) <= 1) {
            tb = m_doc->createElement(X("tb"));
            tc->appendChild(tb);
            tb->setAttribute(X("num"), X(IntToString(tbNum).c_str()));
 
            tbgs = m_doc->createElement(X("tbgs"));
            tb->appendChild(tbgs);
          }
 
          if (pac == nullptr && int(it->_level) == 0) {
            pac = m_doc->createElement(X("pac"));
            tb->appendChild(pac);
            pac->setAttribute(X("num"), X(IntToString(muPACChipNo).c_str()));
          }
 
          DOMElement* mu = m_doc->createElement(X("mu"));
          mu->setAttribute(X("pt"), X(IntToString(int(it->_mu.getPtCode())).c_str()));
          mu->setAttribute(X("qual"), X(IntToString(int(it->_mu.getQuality())).c_str()));
          mu->setAttribute(X("sign"), X(IntToString(int(it->_mu.getSign())).c_str()));
 
          if (int(it->_level) == 0) {
            mu->setAttribute(X("num"), X(IntToString(muSegment).c_str()));
            pac->appendChild(mu);
          } else {
            mu->setAttribute(X("num"), X(IntToString(int(it->_index)).c_str()));
            mu->setAttribute(X("phi"), X(IntToString(int(it->_mu.getPhiAddr())).c_str()));
            mu->setAttribute(X("eta"), X(IntToString(int(it->_mu.getEtaAddr())).c_str()));
            mu->setAttribute(X("gbD"), X(IntToString(int(it->_mu.getGBData())).c_str()));
            if (int(it->_level) == 1) {
              tbgs->appendChild(mu);
            } else if (int(it->_level) == 2) {
              tcgs->appendChild(mu);
            } else {
              throw cms::Exception("RPCMuonsGrabber") << "xx Unexpected level" << std::endl;
            }
          }
 
          it = m_muons.erase(it);
 
        }  // muons iter
        // } // segment
      }  // PAC
    }    // TB
  }      // TC
 
  for (int level = 3; level <= 4; ++level) {
    for (int half = 0; half <= 1; ++half) {
      for (int be = 0; be <= 1; ++be) {  // brl/endcap
 
        std::vector<RPCMuonExtraStruct>::iterator it = m_muons.begin();
        DOMElement* hs = nullptr;
        while (it != m_muons.end()) {
          if ((int(it->_level) != level) || int(it->_hsHalf) != half || int(it->_region) != be) {
            ++it;
            continue;
          }
 
          if (hs == nullptr) {
            if (level == 3) {
              hs = m_doc->createElement(X("hs"));
              hs->setAttribute(X("num"), X(IntToString(half).c_str()));
            } else if (level == 4) {
              hs = m_doc->createElement(X("fs"));
            } else {  // shoudlnt get here
              throw cms::Exception("RPCMuonsGrabber") << "Problem writing out muons - lvl " << level << std::endl;
            }
            hs->setAttribute(X("be"), X(IntToString(be).c_str()));
            currRelBx->appendChild(hs);
          }
 
          DOMElement* mu = m_doc->createElement(X("mu"));
          hs->appendChild(mu);
          mu->setAttribute(X("num"), X(IntToString(int(it->_index)).c_str()));
          mu->setAttribute(X("pt"), X(IntToString(int(it->_mu.getPtCode())).c_str()));
          mu->setAttribute(X("qual"), X(IntToString(int(it->_mu.getQuality())).c_str()));
          mu->setAttribute(X("sign"), X(IntToString(int(it->_mu.getSign())).c_str()));
          mu->setAttribute(X("phi"), X(IntToString(int(it->_mu.getPhiAddr())).c_str()));
          mu->setAttribute(X("eta"), X(IntToString(int(it->_mu.getEtaAddr())).c_str()));
          mu->setAttribute(X("gbD"), X(IntToString(int(it->_mu.getGBData())).c_str()));
 
          //std::cout << int(it->_level) << int(it->_region) << int(it->_hsHalf)
          //   << " " << int(it->_index)
          //   << " " << it->_mu.printDebugInfo(2) << std::endl;
 
          it = m_muons.erase(it);
 
        }  // muons iter
      }    // be iter
    }      //half iteration
  }        // lvl iteration
 
  if (!m_muons.empty()) {
    throw cms::Exception("RPCMuonsGrabber") << " There are still some muons in muons vec" << std::endl;
  }
}
 
std::string MuonsGrabber::IntToString(int i) {
  std::stringstream ss;
  ss << i;
 
  return ss.str();
}