RPCEMap.h

Go to the documentation of this file.

#ifndef RPCEMap_H
#define RPCEMap_H
 
#include "CondFormats/Serialization/interface/Serializable.h"
 
#include "CondFormats/RPCObjects/interface/RPCReadOutMapping.h"
#include <map>
#include <vector>
#include <utility>
#include <string>
#include <iostream>
 
class RPCEMap {
public:
  RPCEMap(const std::string& version = "") : theVersion(version) {}
 
  virtual ~RPCEMap() {}
 
  std::string theVersion;
 
  struct dccItem {
    int theId;
    int nTBs;
 
    COND_SERIALIZABLE;
  };
  struct tbItem {
    int theNum;
    int nLinks;
 
    COND_SERIALIZABLE;
  };
  struct linkItem {
    int theTriggerBoardInputNumber;
    int nLBs;
 
    COND_SERIALIZABLE;
  };
  struct lbItem {
    bool theMaster;
    int theLinkBoardNumInLink;
    int theCode;
    int nFebs;
 
    lbItem() : theMaster(false), theLinkBoardNumInLink(0), theCode(0), nFebs(0){/* nop */};
 
    COND_SERIALIZABLE;
  };
  struct febItem {
    int theLinkBoardInputNum;
    int thePartition;
    int theChamber;
    int theAlgo;
 
    COND_SERIALIZABLE;
  };
 
  std::vector<dccItem> theDccs;
  std::vector<tbItem> theTBs;
  std::vector<linkItem> theLinks;
  std::vector<lbItem> theLBs;
  std::vector<febItem> theFebs;
 
  RPCReadOutMapping const* convert() const {
    RPCReadOutMapping* cabling = new RPCReadOutMapping(theVersion);
    int diskOffset = 4;
    int year = atoi(theVersion.substr(6, 4).c_str());
    int month = atoi(theVersion.substr(3, 2).c_str());
    if (year < 2012 || (year == 2012 && month < 11))
      diskOffset = 3;
    int lastTB = 0;
    int lastLink = 0;
    int lastLB = 0;
    int lastFeb = 0;
    for (unsigned int idcc = 0; idcc < theDccs.size(); idcc++) {
      DccSpec dcc(theDccs[idcc].theId);
      for (int itb = lastTB; itb < lastTB + theDccs[idcc].nTBs; itb++) {
        TriggerBoardSpec tb(theTBs[itb].theNum);
        for (int ilink = lastLink; ilink < lastLink + theTBs[itb].nLinks; ilink++) {
          LinkConnSpec lc(theLinks[ilink].theTriggerBoardInputNumber);
          for (int ilb = lastLB; ilb < lastLB + theLinks[ilink].nLBs; ilb++) {
            LinkBoardSpec lb(theLBs[ilb].theMaster, theLBs[ilb].theLinkBoardNumInLink, theLBs[ilb].theCode);
            for (int ifeb = lastFeb; ifeb < lastFeb + theLBs[ilb].nFebs; ifeb++) {
              int sector = (theFebs[ifeb].theChamber) % 100;
              char subsector = ((theFebs[ifeb].theChamber) / 100) % 10 - 2;
              char febZRadOrnt = ((theFebs[ifeb].theChamber) / 1000) % 5;
              char febZOrnt = ((theFebs[ifeb].theChamber) / 5000) % 2;
              char diskOrWheel = ((theFebs[ifeb].theChamber) / 10000) % 10 - diskOffset;
              char layer = ((theFebs[ifeb].theChamber) / 100000) % 10;
              char barrelOrEndcap = (theFebs[ifeb].theChamber) / 1000000;
              ChamberLocationSpec chamber = {
                  diskOrWheel, layer, sector, subsector, febZOrnt, febZRadOrnt, barrelOrEndcap};
              char cmsEtaPartition = (theFebs[ifeb].thePartition) / 1000;
              char positionInCmsEtaPartition = ((theFebs[ifeb].thePartition) % 1000) / 100;
              char localEtaPartition = ((theFebs[ifeb].thePartition) % 100) / 10;
              char positionInLocalEtaPartition = (theFebs[ifeb].thePartition) % 10;
              FebLocationSpec afeb = {
                  cmsEtaPartition, positionInCmsEtaPartition, localEtaPartition, positionInLocalEtaPartition};
              FebConnectorSpec febConnector(theFebs[ifeb].theLinkBoardInputNum, chamber, afeb);
              febConnector.addStrips(theFebs[ifeb].theAlgo);
              lb.add(febConnector);
              //              std::cout<<"End of FEB"<<std::endl;
            }
            lc.add(lb);
            lastFeb += theLBs[ilb].nFebs;
          }
          tb.add(lc);
          lastLB += theLinks[ilink].nLBs;
        }
        dcc.add(tb);
        lastLink += theTBs[itb].nLinks;
      }
      cabling->add(dcc);
      lastTB += theDccs[idcc].nTBs;
    }
    return cabling;
  };
 
private:
  COND_SERIALIZABLE;
};
 
#endif  // RPCEMap_H