Run2toRun3ConvertersL1.h File Reference

#include <string>

Include dependency graph for Run2toRun3ConvertersL1.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Namespaces
namespace	TXC
namespace	TrigConf
	Forward iterator to traverse the main components of the trigger configuration.

Functions
void	convertRun2L1MenuToRun3 (const TrigConf::CTPConfig *ctpConfig, const TXC::L1TopoMenu *topoMenu, const std::string &filename, bool writeTmpFile=false)
	Run 2 to Run 3 L1 menu converter Converts Run 2 L1 menu and L1Topo menu into run 3 format and writes out a json file.
void	convertRun2BunchGroupsToRun3 (const TrigConf::CTPConfig *ctpConfig, const std::string &filename, bool writeTmpFile=false)
	Run 2 to Run 3 bunchgroup converter.
void	convertRun2L1PrescalesToRun3 (const TrigConf::CTPConfig *ctpConfig, const std::string &filename, bool writeTmpFile=false)
	Run 2 to Run 3 L1 prescale converter.

Function Documentation

convertRun2BunchGroupsToRun3()

void convertRun2BunchGroupsToRun3	(	const TrigConf::CTPConfig *	ctpConfig,
		const std::string &	filename,
		bool	writeTmpFile = false )

Run 2 to Run 3 bunchgroup converter.

Parameters

ctpConfig	pointer to L1 menu (which contains the bunchgroups)
filename	name of the output json file

Definition at line 483 of file Run2toRun3ConvertersL1.cxx.

{
   const TrigConf::BunchGroupSet & bgs = ctpConfig->bunchGroupSet();
   if(bgs.bunchGroups().size()==0) {
      std::cout << "BunchgroupSet is empty, no file will be produced" << std::endl;
      return;
   }
 
   json bgset;
   bgset["filetype"] ="bunchgroupset";
   bgset["name"] = bgs.name();
 
   json jGroups;
   const std::vector<TrigConf::BunchGroup>& bgVec = bgs.bunchGroups();
   for (auto & group : bgVec) {
      auto ranges = toRanges(group.bunches());
      json jGroup;
      jGroup["name"] =group.name();
      jGroup["id"] =group.internalNumber();
      jGroup["info"] = std::to_string(group.bunches().size()) + " bunches, " + std::to_string(ranges.size()) + " groups";
 
      json jBCIDS = json::array_t{};
      for ( auto [start, end] : ranges) {
         jBCIDS += json{{"first", start}, {"length", 1+end-start}};
      }
      jGroup["bcids"] =  std::move(jBCIDS);
 
      jGroups[std::string("BGRP")+std::to_string(group.internalNumber())] = std::move(jGroup);
   }
   bgset["bunchGroups"] = std::move(jGroups);
 
   if(writeTmpFile) {
      std::ofstream outfile("tmp" + filename);
      outfile << std::setw(4) << bgset << std::endl;
      std::cout << "Wrote tmp" << filename << std::endl;
   }
   std::stringstream ss;
   ss << bgset;
   ptree top;
   boost::property_tree::read_json(ss, top);
   TrigConf::L1BunchGroupSet newbgs(std::move(top));
   TrigConf::JsonFileWriterL1 writer;
   writer.writeJsonFile(filename, newbgs);
}

convertRun2L1MenuToRun3()

void convertRun2L1MenuToRun3	(	const TrigConf::CTPConfig *	ctpConfig,
		const TXC::L1TopoMenu *	topoMenu,
		const std::string &	filename,
		bool	writeTmpFile )

Run 2 to Run 3 L1 menu converter Converts Run 2 L1 menu and L1Topo menu into run 3 format and writes out a json file.

Parameters

ctpConfig	pointer to L1 menu
topoMenu	pointer to L1Topo menu (can be empty)
filename	name of the output json file
writeTmpFile	writes out a temporary version of the json file (before filling the L1Menu structure)

Run 2 to Run 3 L1 menu converter Converts Run 2 L1 menu and L1Topo menu into run 3 format and writes out a json file.

Definition at line 88 of file Run2toRun3ConvertersL1.cxx.

{
 
   if( !ctpConfig ) {
      std::cout << "No CTPConfig, no L1Menu file will be produced" << std::endl;
      return;
   }
 
   unsigned int run = ctpConfig->ctpVersion() <=3 ? 1 : 2;
 
   // items
   json items = json::object_t{};
   for( const TrigConf::TriggerItem * sourceItem: ctpConfig->menu().items() ) {
      json item({});
      item["name"] = sourceItem->name();
      item["legacy"] = true;
      item["ctpid"] = sourceItem->ctpId();
      auto [definition, bunchgroups] = decodeItemDefinition(sourceItem);
      item["definition"] = definition;
      item["bunchgroups"] = bunchgroups;
      item["triggerType"] = TrigConf::uint2bin(sourceItem->triggerType(),8);
      item["partition"] = sourceItem->partition();
      std::string smon("LF:");
      smon += (sourceItem->monitor() & 0x04 ? '1' : '0');
      smon += (sourceItem->monitor() & 0x02 ? '1' : '0');
      smon += (sourceItem->monitor() & 0x01 ? '1' : '0');
      smon += "|HF:";
      smon += (sourceItem->monitor() & 0x20 ? '1' : '0');
      smon += (sourceItem->monitor() & 0x10 ? '1' : '0');
      smon += (sourceItem->monitor() & 0x08 ? '1' : '0');
      item["monitor"] = smon;
      items[sourceItem->name()] = item;
   };
 
   // thresholds
   json thresholds = json::object_t{};
   std::unordered_set<std::string> legacyThrTypes{"EM", "TAU", "JET", "JB", "JE", "JF", "XE", "XS", "TE", "ZB", "R2TOPO"};
   for( const TrigConf::TriggerThreshold * sourceThr : ctpConfig->menu().thresholdVector()) {
      std::string thrType = sourceThr->type();
      if(thrType=="BGRP" || thrType=="RNDM") {
         continue;
      }
      if(thrType=="MUON") {
         thrType = "MU";
      } 
      else if(thrType=="TOPO") {
         thrType = "R2TOPO";
      }
 
      bool isLegacyCalo = legacyThrTypes.count(thrType) > 0;
      json &jThisType = isLegacyCalo ? thresholds["legacyCalo"][thrType] : thresholds[thrType];
      bool firstOfItsType = jThisType.empty();
      if (firstOfItsType)
      {
         jThisType["type"] = thrType;
      }
      json thr;
      size_t mapping = sourceThr->mapping();
      thr["mapping"] = mapping;
      if(thrType=="MU") {
         if (firstOfItsType)
         {
            jThisType["type"] = "MU";
            jThisType["exclusionLists"] = json::object_t{};
            jThisType["roads"]["rpc"] = json::object_t{};
            jThisType["roads"]["tgc"] = json::object_t{};
         }
         int ptCut = sourceThr->triggerThresholdValue(0,0)->ptCutCount();
         thr["baThr"] = ptCut;
         thr["ecThr"] = ptCut;
         thr["fwThr"] = ptCut;
         thr["baIdx"] = mapping;
         thr["ecIdx"] = mapping;
         thr["fwIdx"] = mapping;
         thr["tgcFlags"] = "";
         thr["rpcFlags"] = "";
         thr["region"] = "ALL";
         jThisType["roads"]["rpc"][std::to_string(ptCut)] = mapping;
         jThisType["roads"]["tgc"][std::to_string(ptCut)] = mapping;
      }
      else if (thrType == "EM")
      {
         thr["thrValues"] = json::array_t{};
         for (const TrigConf::TriggerThresholdValue *tv : sourceThr->thresholdValueVector())
         {
            auto cl = dynamic_cast<const TrigConf::ClusterThresholdValue *>(tv);
            if (not cl) continue;
            json jtv;
            jtv["value"] = static_cast<unsigned int>(tv->ptcut());
            std::string isobits = "00000";
            auto isomask = cl->isolationMask();
            for (size_t b = 0; b < 5; ++b)
            {
               if (isomask & (1 << b))
               {
                  isobits[4 - b] = '1';
               }
            }
            jtv["isobits"] = isobits;
            jtv["etamin"] = tv->etamin();
            jtv["etamax"] = tv->etamax();
            jtv["phimin"] = tv->phimin();
            jtv["phimax"] = tv->phimax();
            jtv["priority"] = static_cast<unsigned int>(tv->priority());
            thr["thrValues"] += jtv;
         }
      }
      else if (thrType == "TAU")
      {
         auto cl = dynamic_cast<const TrigConf::ClusterThresholdValue*>(sourceThr->triggerThresholdValue(0,0));
         if (not cl) continue;
         int ptCut = (int)cl->ptcut();
         thr["value"] = ptCut;
         std::string isobits = "00000";
         auto isomask = cl->isolationMask();
         for (size_t b = 0; b < 5; ++b)
         {
            if (isomask & (1 << b))
            {
               isobits[4 - b] = '1';
            }
         }
         thr["isobits"] = isobits;
      }
      else if (thrType == "JET" || thrType == "JB" ||  thrType == "JF")
      {
         thr["thrValues"] = json::array_t{};
         for (const TrigConf::TriggerThresholdValue *tv : sourceThr->thresholdValueVector())
         {
            auto jetThrVal = dynamic_cast<const TrigConf::JetThresholdValue *>(sourceThr->triggerThresholdValue(0, 0));
            if (not jetThrVal) continue;
            json jtv;
            jtv["value"] = static_cast<unsigned int>(tv->ptcut());
            jtv["etamin"] = tv->etamin();
            jtv["etamax"] = tv->etamax();
            jtv["phimin"] = tv->phimin();
            jtv["phimax"] = tv->phimax();
            jtv["window"] = jetThrVal->windowSize() == TrigConf::JetWindowSize::SMALL ? 4 : 8; // 4 and 8 are the two possible jet window sizes
            jtv["priority"] = static_cast<unsigned int>(tv->priority());
            thr["thrValues"] += jtv;
         }
      }
      else if (thrType == "XE" || thrType == "TE" ||  thrType == "XS" || thrType == "JE")
      {
         thr["value"] = static_cast<unsigned int>(sourceThr->thresholdValueVector().at(0)->ptcut());
      }
      else if (thrType == "ZB")
      {
         thr["seed"] = sourceThr->zbSeedingThresholdName();
         thr["seedMultiplicity"] = sourceThr->zbSeedingThresholdMulti();
         thr["seedBcdelay"] = sourceThr->bcDelay();
      }
      jThisType["thresholds"][sourceThr->name()] = thr;
   }
 
   // extra info for thresholds
   const TrigConf::CaloInfo& ci = ctpConfig->menu().caloInfo();
   thresholds["legacyCalo"]["EM"]["ptMinToTopo"] = ci.minTobEM().ptmin;
   thresholds["legacyCalo"]["EM"]["resolutionMeV"] = (int)(1000/ci.globalEmScale());
   thresholds["legacyCalo"]["TAU"]["ptMinToTopo"] = ci.minTobTau().ptmin;
   thresholds["legacyCalo"]["JET"]["ptMinToTopoSmallWindow"] = ci.minTobJetSmall().ptmin;
   thresholds["legacyCalo"]["JET"]["ptMinToTopoLargeWindow"] = ci.minTobJetLarge().ptmin;
   if(run==2) {
      json isoHAforEM{ {"thrtype", "HAIsoForEMthr"}, {"Parametrization", json::array_t{}} };
      json isoEMforEM{ {"thrtype", "EMIsoForEMthr"}, {"Parametrization", json::array_t{}} };
      json isoEMforTAU{ {"thrtype", "EMIsoForTAUthr"}, {"Parametrization", json::array_t{}} };
      for(const TrigConf::IsolationParam & iso : ci.isolationHAIsoForEMthr()) {
         json p{ {"etamax", iso.etamax()}, {"etamin", iso.etamin()}, {"isobit", iso.isobit()}, {"mincut", iso.mincut()}, 
               {"offset", iso.offset()}, {"priority", iso.priority()}, {"slope", iso.slope()}, {"upperlimit", iso.upperlimit()} };
         isoHAforEM["Parametrization"] += p;
      }
      for(const TrigConf::IsolationParam & iso : ci.isolationEMIsoForEMthr()) {
         json p{ {"etamax", iso.etamax()}, {"etamin", iso.etamin()}, {"isobit", iso.isobit()}, {"mincut", iso.mincut()}, 
               {"offset", iso.offset()}, {"priority", iso.priority()}, {"slope", iso.slope()}, {"upperlimit", iso.upperlimit()} };
         isoEMforEM["Parametrization"] += p;
      }
      for(const TrigConf::IsolationParam & iso : ci.isolationEMIsoForTAUthr()) {
         json p{ {"etamax", iso.etamax()}, {"etamin", iso.etamin()}, {"isobit", iso.isobit()}, {"mincut", iso.mincut()}, 
               {"offset", iso.offset()}, {"priority", iso.priority()}, {"slope", iso.slope()}, {"upperlimit", iso.upperlimit()} };
         isoEMforTAU["Parametrization"] += p;
      }
      thresholds["legacyCalo"]["EM"]["isolation"]["HAIsoForEMthr"] =  std::move(isoHAforEM);
      thresholds["legacyCalo"]["EM"]["isolation"]["EMIsoForEMthr"] =  std::move(isoEMforEM);
      thresholds["legacyCalo"]["TAU"]["isolation"]["EMIsoForTAUthr"] =  std::move(isoEMforTAU);
   }
   const TrigConf::METSigParam &xs = ci.metSigParam();
   thresholds["legacyCalo"]["XS"]["significance"] = json::object_t{
       {"xsSigmaScale", xs.xsSigmaScale()}, {"xsSigmaOffset", xs.xsSigmaOffset()}, {"xeMin", xs.xeMin()}, {"xeMax", xs.xeMax()}, {"teSqrtMin", xs.teSqrtMin()}, {"teSqrtMax", xs.teSqrtMax()}};
 
   // boards
   json boards = json::object_t{};
   boards["Ctpin7"] = json::object_t{ {"type", "CTPIN"}, {"legacy", true},
                                      {"connectors", std::vector<std::string>{"EM1", "EM2", "TAU1", "TAU2"}} };
   boards["Ctpin8"] = json::object_t{ {"type", "CTPIN"}, {"legacy", true},
                                      {"connectors", std::vector<std::string>{"JET1", "JET2", "EN1", "EN2"}} };
   boards["Ctpin9"] = json::object_t{ {"type", "CTPIN"}, {"legacy", true},
                                      {"connectors", std::vector<std::string>{"MUCTPI", "CTPCAL", "NIM1", "NIM2"}} };
   boards["LegacyTopo0"] = json::object_t{ {"type", "TOPO"}, {"legacy", true},
                                           {"connectors", std::vector<std::string>{"LegacyTopo0"}} };
   boards["LegacyTopo1"] = json::object_t{ {"type", "TOPO"}, {"legacy", true},
                                           {"connectors", std::vector<std::string>{"LegacyTopo1"}} };
 
   // connectors
   json connectors = json::object_t{};
   std::map<std::string,std::vector<const TrigConf::TriggerThreshold*>> triggerlinesMap;
   for( const TrigConf::TriggerThreshold * sourceThr : ctpConfig->menu().thresholdVector()) {
      if(sourceThr->isInternal()) {
         continue;
      }
      std::string cableName = sourceThr->cableName();
      if(cableName=="ALFA") {
         cableName = "AlfaCtpin";
      }
      else if(cableName=="TOPO1") {
         cableName = "LegacyTopo0";
      }
      else if(cableName=="TOPO2") {
         cableName = "LegacyTopo1";
      }
      triggerlinesMap[cableName].push_back(sourceThr);
   }
   for( auto & [type, triggerlines] : triggerlinesMap) {
      std::sort(std::begin(triggerlines), std::end(triggerlines),
                [](const TrigConf::TriggerThreshold *a, const TrigConf::TriggerThreshold *b) { return a->cableStart() < b->cableStart(); });
      for( const TrigConf::TriggerThreshold * thr : triggerlines) {
         if(!connectors.contains(type)) {
            if(thr->input()=="ctpcore") {
               connectors[type]["type"] =  "electrical";
               if(type=="AlfaCtpin") {
                  connectors[type]["triggerlines"]["clock0"] = json::array_t{};
                  connectors[type]["triggerlines"]["clock1"] = json::array_t{};
               } else {
                  connectors[type]["triggerlines"]["fpga0"]["clock0"] = json::array_t{};
                  connectors[type]["triggerlines"]["fpga0"]["clock1"] = json::array_t{};
                  connectors[type]["triggerlines"]["fpga1"]["clock0"] = json::array_t{};
                  connectors[type]["triggerlines"]["fpga1"]["clock1"] = json::array_t{};
               }
            } else if(thr->input()=="ctpin") {
               connectors[type]["type"] =  "ctpin";
               connectors[type]["triggerlines"] = json::array_t{};
            } else {
               throw std::runtime_error("Unknown connector type" + thr->input());
            }
            connectors[type]["legacy"] = true;
         }
         size_t start = thr->cableStart();
         size_t nbits = thr->cableEnd() - thr->cableStart() + 1;
         if(thr->input()=="ctpcore") {
            unsigned int clock = thr->clock();
            if(type=="AlfaCtpin") {
               json tl{{"name", thr->name()}, {"startbit", start}, {"nbits", nbits}};
               connectors[type]["triggerlines"]["clock" + std::to_string(clock)] += tl;
            } else {
               size_t fpga = 0;
               if(start>=16) {
                  start -= 16;
                  fpga++;
               }
               json tl{{"name", "R2TOPO_" + thr->name()}, {"startbit", start}, {"nbits", nbits}};
               connectors[type]["triggerlines"]["fpga" + std::to_string(fpga)]["clock" + std::to_string(clock)] += tl;
            }
         } else {
            json tl{{"name", thr->name()}, {"startbit", start}, {"nbits", nbits}};
            connectors[type]["triggerlines"] += tl;
         }
      }
   }
 
   // algorithms
   json::object_t decAlgos;
   json::object_t sortAlgos;
   for(const TXC::L1TopoConfigAlg & alg : topoMenu->getL1TopoConfigAlgs()) {
      json jAlg;
      jAlg["algId"] = alg.algoID();
      jAlg["klass"] = alg.type();
      jAlg["fixedParameters"]["generics"] = json::object_t{};
      size_t pos{0};
      for(const TXC::FixedParameter & fixP : alg.getFixedParameters()) {
         try
         {
            auto vInt = std::stoi(fixP.value);
            jAlg["fixedParameters"]["generics"][fixP.name] = json::object_t{{"value", vInt}, {"position", pos++}};
         }
         catch (std::invalid_argument &)
         {
            jAlg["fixedParameters"]["generics"][fixP.name] = json::object_t{{"value", fixP.value}, {"position", pos++}};
         }
      }
      jAlg["variableParameters"] = json::array_t{};
      for(const TXC::RegisterParameter & regP : alg.getParameters()) {
         // Work around overflow in the database...
         if (regP.name == "MaxMSqr") {
           unsigned val = std::stoul(regP.value);
           if (val >= 1u<<31) {
             val = 999;
           }
           jAlg["variableParameters"] += json::object_t{{"name", regP.name}, {"selection", regP.selection}, {"value", val}};
         }
         else {
           jAlg["variableParameters"] += json::object_t{{"name", regP.name}, {"selection", regP.selection}, {"value", std::stoi(regP.value)}};
         }
      }
      if(alg.isSortAlg()) {
         jAlg["input"] = alg.getInputNames()[0];
         jAlg["output"] = alg.output();
         sortAlgos[alg.name()] =  std::move(jAlg);
      } else if(alg.isDecAlg()) {
         jAlg["input"] = alg.getInputNames();
         jAlg["output"] = alg.getOutputNames();
         decAlgos[alg.name()] =  std::move(jAlg);
      }
   }
   json topo = json::object_t{};
   if(run==2) {
      topo["R2TOPO"]["decisionAlgorithms"] = decAlgos;
      topo["R2TOPO"]["sortingAlgorithms"] = sortAlgos;
   }
 
   // ctp
   json ctp = json::object_t{};
   ctp["inputs"]["ctpin"]["slot7"] = json::object_t{{"connector0", "EM1"}, {"connector1", "EM2"}, 
                                                    {"connector2", "TAU1"}, {"connector3", "TAU2"}};
   ctp["inputs"]["ctpin"]["slot8"] = json::object_t{{"connector0", "JET1"}, {"connector1", "JET2"}, 
                                                    {"connector2", "EN1"}, {"connector3", "EN2"}};
   ctp["inputs"]["ctpin"]["slot9"] = json::object_t{{"connector0", "MUCTPI"}, {"connector1", "CTPCAL"}, 
                                                    {"connector2", "NIM1"}, {"connector3", "NIM2"}};
   if(run==2) {
      ctp["inputs"]["electrical"] = json::object_t{{"connector0", "AlfaCtpin"}, {"connector1", "LegacyTopo0"}, {"connector2", "LegacyTopo1"}};
   } else {
      ctp["inputs"]["electrical"] = json::object_t{};
   }
   ctp["inputs"]["optical"] = json::object_t{};
   ctp["monitoring"]["ctpmon"] = json::object_t{};
 
   // putting the menu together
   json menu = json::object_t{};
   menu["filetype"] = "l1menu";
   menu["run"] = run;
   menu["name"] = ctpConfig->name();
   menu["items"] =  std::move(items);
   menu["thresholds"] =  std::move(thresholds);
   menu["topoAlgorithms"] =  std::move(topo);
   menu["boards"] =  std::move(boards);
   menu["connectors"] =  std::move(connectors);
   menu["ctp"] =  std::move(ctp);
 
 
   if(writeTmpFile) {
      std::ofstream outfile("tmp" + filename);
      outfile << std::setw(4) << menu << std::endl;
      std::cout << "Wrote tmp" << filename << std::endl;
   }
   std::stringstream ss;
   ss << menu;
   ptree top;
   boost::property_tree::read_json(ss, top);
   TrigConf::L1Menu l1menu(std::move(top));
   TrigConf::JsonFileWriterL1 writer;
   writer.writeJsonFile(filename, l1menu);
}

convertRun2L1PrescalesToRun3()

void convertRun2L1PrescalesToRun3	(	const TrigConf::CTPConfig *	ctpConfig,
		const std::string &	filename,
		bool	writeTmpFile )

Run 2 to Run 3 L1 prescale converter.

Parameters

ctpConfig	pointer to L1 menu (which contains the prescales)
filename	name of the output json file

Run 2 to Run 3 L1 prescale converter.

Definition at line 532 of file Run2toRun3ConvertersL1.cxx.

                                                                                                               {
   const TrigConf::PrescaleSet & l1pss = ctpConfig->prescaleSet();
   if(l1pss.isNull()) {
      std::cout << "L1PrescaleSet is empty, no file will be produced" << std::endl;
      return;
   }
 
   json psset;
   psset["filetype"] = "l1prescale";
   psset["name"] = l1pss.name();
   json jCuts;
   for (size_t id = 0; id < l1pss.prescales_float().size(); ++id) {
      json jCut;
      auto itemPtr = ctpConfig->menu().item(id);
      if ( itemPtr != nullptr ) {
         int32_t cut = l1pss.cuts().at(id);
         jCut["cut"] = abs(cut);
         jCut["enabled"] = cut > 0;
         double ps = static_cast<double>(0xFFFFFF) / ( 0x1000000 - cut );
         jCut["info"] = "prescale: "+std::to_string(ps);
         jCuts[itemPtr->name()] = std::move(jCut);
      }
   }
   psset["cutValues"] = std::move(jCuts);
 
   if(writeTmpFile) {
      std::ofstream outfile("tmp" + filename);
      outfile << std::setw(4) << psset << std::endl;
      std::cout << "Wrote tmp" << filename << std::endl;
   }
   std::stringstream ss;
   ss << psset;
   ptree top;
   boost::property_tree::read_json(ss, top);
   TrigConf::L1PrescalesSet ps(std::move(top));
   TrigConf::JsonFileWriterL1 writer;
   writer.writeJsonFile(filename, ps);
}