TrigConf::JsonFileWriterL1 Class Reference

Loader of trigger configurations from Json files. More...

#include <JsonFileWriterL1.h>

Inheritance diagram for TrigConf::JsonFileWriterL1:

Collaboration diagram for TrigConf::JsonFileWriterL1:

Public Member Functions
	JsonFileWriterL1 ()
	Constructor.
virtual	~JsonFileWriterL1 () override=default
	Destructor.
bool	writeJsonFile (const std::string &filename, const L1Menu &l1menu) const
bool	writeJsonFile (const std::string &filename, const L1BunchGroupSet &l1bgs) const
bool	writeJsonFile (const std::string &filename, const L1PrescalesSet &l1bgs) const
bool	msgLvl (const MSGTC::Level lvl) const
	Test the output level.
MsgStreamTC &	msg () const
	The standard message stream.
MsgStreamTC &	msg (const MSGTC::Level lvl) const
	The standard message stream.
const std::string &	getName () const
	name accessor

Private Attributes
boost::thread_specific_ptr< MsgStreamTC >	m_msg_tls
	MsgStreamTC instance (a std::cout like with print-out levels)
std::string	m_name

Detailed Description

Loader of trigger configurations from Json files.

Definition at line 25 of file JsonFileWriterL1.h.

Constructor & Destructor Documentation

JsonFileWriterL1()

TrigConf::JsonFileWriterL1::JsonFileWriterL1

(

)

Constructor.

Definition at line 19 of file JsonFileWriterL1.cxx.

                                           : 
   TrigConfMessaging( "JsonFileWriterL1")
{}

~JsonFileWriterL1()

virtual TrigConf::JsonFileWriterL1::~JsonFileWriterL1 ( )

overridevirtualdefault

Destructor.

Member Function Documentation

getName()

const std::string & TrigConf::TrigConfMessaging::getName ( ) const

inlineinherited

name accessor

Returns

the name

Definition at line 101 of file TrigConfMessaging.h.

                                                           {
    return m_name;
  }

msg() [1/2]

MsgStreamTC & TrigConf::TrigConfMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 86 of file TrigConfMessaging.h.

  {
    MsgStreamTC* ms = m_msg_tls.get();
    if (!ms) {
      ms = new MsgStreamTC(m_name);
      m_msg_tls.reset(ms);
    }
    return *ms;
  }

msg() [2/2]

MsgStreamTC & TrigConf::TrigConfMessaging::msg ( const MSGTC::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 96 of file TrigConfMessaging.h.

  {
    return msg() << lvl;
  }

msgLvl()

bool TrigConf::TrigConfMessaging::msgLvl ( const MSGTC::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicting if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 75 of file TrigConfMessaging.h.

  {
    if (msg().level() <= lvl) {
      msg() << lvl;
      return true;
    }
    else {
      return false;
    }
  }

writeJsonFile() [1/3]

bool TrigConf::JsonFileWriterL1::writeJsonFile	(	const std::string &	filename,
		const L1BunchGroupSet &	l1bgs ) const

Definition at line 860 of file JsonFileWriterL1.cxx.

                                                                                                               {
 
   json j;
   j["filetype"] = "bunchgroupset";
   j["name"]  = l1bgs.name();
 
   json groups;
   for (size_t i = 0 ; i< l1bgs.size(); ++i) {
      const auto & group = l1bgs.getBunchGroup(i);
      json jgroup({});
      jgroup["name"] = group->name();
      jgroup["id"] = group->id();
      jgroup["info"] = std::to_string(group->size()) + " bunches, " + std::to_string(group->nGroups()) + " groups";
      json trains = json::array();
      for (auto [first, len]: group->trains()) {
         json train({});
         train["first"] = first;
         train["length"] = len;
         trains.push_back(train);
      }
      jgroup["bcids"] = std::move(trains);
      groups["BGRP"+std::to_string(group->id())] = std::move(jgroup);
   }
   j["bunchGroups"] = std::move(groups);
   std::ofstream outfile(filename);
   outfile << std::setw(4) << j << std::endl;
   TRG_MSG_INFO("Saved file " << filename);
   return true;
}

writeJsonFile() [2/3]

bool TrigConf::JsonFileWriterL1::writeJsonFile	(	const std::string &	filename,
		const L1Menu &	l1menu ) const

Definition at line 25 of file JsonFileWriterL1.cxx.

{
 
   json items = json::object_t{};
   for( auto & item : l1menu ) {
      json jItem({});
      jItem["name"] = item.name();
      jItem["ctpid"] = item.ctpId();
      jItem["definition"] = item.definition();
      jItem["monitor"] = item.monitor();
      jItem["partition"] = item.partition();
      jItem["triggerType"] = item.triggerType();
      jItem["bunchgroups"] = json(item.bunchgroups());
      if(auto legacy = item.legacy() )
         jItem["legacy"] = *legacy;
      items[item.name()] = jItem;
   };
 
   json thresholds = json::object_t{};
   for(const std::string & thrType : l1menu.thresholdTypes()) {
      json jThresholsByType({});
      // first the thresholds of this type
      for(auto & thr : l1menu.thresholds(thrType)) {
         json jThr({});
         jThr["mapping"] = thr->mapping();
 
         if(thr->hasChild("sectors")) { // for MBTS_A and MBTS_C
            std::vector<std::string> sectors; 
            for(auto & s : thr->getList("sectors")) {
               sectors.push_back(s.getValue());
            }
            jThr["sectors"] = sectors;
         } 
         if(thr->hasAttribute("voltage")) { // for MBTSII
            jThr["voltage"] = thr->getAttribute<float>("voltage");
         }
 
         // value
         try {
            if(thr->hasAttribute("value")) {
               auto & caloThr = dynamic_cast<const TrigConf::L1Threshold_Calo &>(*thr); // for MBTSII
               jThr["value"] = int(caloThr.thrValue());
            }
         } catch(std::bad_cast&) {};
 
         // MU
         try {
            const auto & muThr = dynamic_cast<const TrigConf::L1Threshold_MU &>(*thr);
            jThr["baThr"] = muThr.ptBarrel();
            jThr["ecThr"] = muThr.ptEndcap();
            jThr["fwThr"] = muThr.ptForward(); 
            jThr["baIdx"] = muThr.idxBarrel();
            jThr["ecIdx"] = muThr.idxEndcap();
            jThr["fwIdx"] = muThr.idxForward();
            jThr["region"] = muThr.region();
            jThr["tgcFlags"] = muThr.tgcFlags();
            jThr["rpcFlags"] = muThr.rpcFlags();
            if(const std::string & roiExcl = muThr.rpcExclROIList(); !roiExcl.empty()) {
               jThr["rpcExclROIList"] = roiExcl;
            }
         } catch(std::bad_cast&) {};
 
         // TAU
         try {
            const auto & tauThr = dynamic_cast<const TrigConf::L1Threshold_TAU &>(*thr);
            std::string isobits = "00000";
            auto isomask = tauThr.isolationMask();
            for(size_t b=0; b<5; ++b) {
               if(isomask & (1<<b)) { isobits[4-b] = '1'; }
            }
            jThr["isobits"] = isobits;
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : tauThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["phimin"] = 0; // never used, so not read 
               jRV["phimax"] = 64; // never used, so not read
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // EM
         try {
            const auto & EMThr = dynamic_cast<const TrigConf::L1Threshold_EM &>(*thr);
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : EMThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["phimin"] = 0; // never used, so not read 
               jRV["phimax"] = 64; // never used, so not read
               std::string isobits = "00000";
               auto isomask = EMThr.isolationMask(rv.etaMin());
               for(size_t b=0; b<5; ++b) {
                  if(isomask & (1<<b)) { isobits[4-b] = '1'; }
               }
               jRV["isobits"] = isobits;
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // JET
         try {
            const auto & JThr = dynamic_cast<const TrigConf::L1Threshold_JET &>(*thr);
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : JThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["phimin"] = 0; // never used, so not read 
               jRV["phimax"] = 64; // never used, so not read
               jRV["window"] = JThr.window(0);
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // TE
         try {
            const auto & teThr = dynamic_cast<const TrigConf::L1Threshold_TE &>(*thr);
            //jThr["thrValues"] = json::array_t({});
            for(auto & rv : teThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // ZB
         try {
            const auto & zbThr = dynamic_cast<const TrigConf::L1Threshold_ZB &>(*thr);
            jThr["seed"] = zbThr.seed();
            jThr["seedBcdelay"] = zbThr.seedBcdelay();
            jThr["seedMultiplicity"] = zbThr.seedMultiplicity();
         } catch(std::bad_cast&) {};
 
         // ZBTopo
         try {
            const auto & zbTopoThr = dynamic_cast<const TrigConf::L1Threshold_ZBTopo &>(*thr);
            jThr["mask0"] = zbTopoThr.mask0();
            jThr["mask0"] = zbTopoThr.mask0();
            jThr["mask0"] = zbTopoThr.mask0();
            jThr["mask0"] = zbTopoThr.mask0();
            jThr["mask0"] = zbTopoThr.mask0();
            jThr["mask0"] = zbTopoThr.mask0();
            jThr["seedBcdelay"] = zbTopoThr.seedBcdelay();
            
         } catch(std::bad_cast&) {};
         
         // eEM
         try {
            const auto & eEMThr = dynamic_cast<const TrigConf::L1Threshold_eEM &>(*thr);
            jThr["reta"] = TrigConf::Selection::wpToString(eEMThr.reta());
            jThr["rhad"] = TrigConf::Selection::wpToString(eEMThr.rhad());
            jThr["wstot"] = TrigConf::Selection::wpToString(eEMThr.wstot());
            jThr["thrValues"] = json::value_t::array;
            for(auto & rv : eEMThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // jEM
         try {
            const auto & jEMThr = dynamic_cast<const TrigConf::L1Threshold_jEM &>(*thr);
            jThr["iso"] = TrigConf::Selection::wpToString(jEMThr.iso());
            jThr["frac"] = TrigConf::Selection::wpToString(jEMThr.frac());
            jThr["frac2"] = TrigConf::Selection::wpToString(jEMThr.frac2());
            jThr["thrValues"] = json::value_t::array;
            for(auto & rv : jEMThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // eTAU
         try {
            const auto & eTAUThr = dynamic_cast<const TrigConf::L1Threshold_eTAU &>(*thr);
            jThr["rCore"] = TrigConf::Selection::wpToString(eTAUThr.rCore());
            jThr["rHad"] = TrigConf::Selection::wpToString(eTAUThr.rHad());
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : eTAUThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["phimin"] = 0; // never used, so not read 
               jRV["phimax"] = 64; // never used, so not read
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
         
         // jTAU
         try {
            const auto & jTAUThr = dynamic_cast<const TrigConf::L1Threshold_jTAU &>(*thr);
            jThr["isolation"] = TrigConf::Selection::wpToString(jTAUThr.isolation());
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : jTAUThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["phimin"] = 0; // never used, so not read 
               jRV["phimax"] = 64; // never used, so not read
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};         
 
         // cTAU
         try {
            const auto & cTAUThr = dynamic_cast<const TrigConf::L1Threshold_cTAU &>(*thr);
            jThr["isolation"] = TrigConf::Selection::wpToString(cTAUThr.isolation());
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : cTAUThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["phimin"] = 0; // never used, so not read 
               jRV["phimax"] = 64; // never used, so not read
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // jJ
         try {
            const auto & jJThr = dynamic_cast<const TrigConf::L1Threshold_jJ &>(*thr);
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : jJThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // jLJ
         try {
            const auto & jLJThr = dynamic_cast<const TrigConf::L1Threshold_jLJ &>(*thr);
            jThr["thrValues"] = json::array_t({});
            for(auto & rv : jLJThr.thrValues()) {
               json jRV({});
               jRV["value"] = static_cast<unsigned int>(rv.value());
               jRV["etamin"] = rv.etaMin();
               jRV["etamax"] = rv.etaMax();
               jRV["priority"] = rv.priority();
               jThr["thrValues"] += jRV;
            }
         } catch(std::bad_cast&) {};
 
         // gJ
         try {
            const auto & gJThr = dynamic_cast<const TrigConf::L1Threshold_gJ &>(*thr);
            jThr["value"] = int(gJThr.thrValue());
         } catch(std::bad_cast&) {};
 
         // gLJ
         try {
            const auto & gLJThr = dynamic_cast<const TrigConf::L1Threshold_gLJ &>(*thr);
            jThr["value"] = int(gLJThr.thrValue());
         } catch(std::bad_cast&) {};         
 
         // jXE
         try {
            const auto & jXEThr = dynamic_cast<const TrigConf::L1Threshold_jXE &>(*thr);
            jThr["value"] = int(jXEThr.thrValue());
         } catch(std::bad_cast&) {};
 
         // jTE
         try {
            const auto & jTEThr = dynamic_cast<const TrigConf::L1Threshold_jTE &>(*thr);
            jThr["value"] = int(jTEThr.thrValue());
         } catch(std::bad_cast&) {};         
         
         // gXE
         try {
            const auto & gXEThr = dynamic_cast<const TrigConf::L1Threshold_gXE &>(*thr);
            jThr["value"] = int(gXEThr.thrValue());
         } catch(std::bad_cast&) {};         
 
         // gTE
         try {
            const auto & gTEThr = dynamic_cast<const TrigConf::L1Threshold_gTE &>(*thr);
            jThr["value"] = int(gTEThr.thrValue());
         } catch(std::bad_cast&) {};
 
         jThresholsByType[thr->name()] = std::move(jThr);
      };
      json jThrType({});
      if(thrType != "internal") {
         jThrType["thresholds"] = std::move(jThresholsByType);
      }
      jThrType["type"] = thrType;
 
      // extra info 
      auto & extraInfo = l1menu.thrExtraInfo().thrExtraInfo(thrType);
      if(extraInfo.hasAttribute("ptMinToTopo")) { // for MBTSII
         jThrType["ptMinToTopo"] = extraInfo.getAttribute<unsigned int>("ptMinToTopo");
      }
      if(extraInfo.hasAttribute("resolutionMeV")) { // for MBTSII
         jThrType["resolutionMeV"] = extraInfo.getAttribute<unsigned int>("resolutionMeV");
      }
 
      // extra info using type specific accessors
      if(thrType == "internal") {
         jThrType["names"] = json::array_t({});
         for(auto & thr : l1menu.thresholds(thrType)) {
            jThrType["names"] += thr->name();
         }
         auto & extInfo = l1menu.thrExtraInfo().thrExtraInfo("internal");
         if(auto randoms = extInfo.getExtraInfo("randoms")) {
            for(size_t rc=0; rc<4; ++rc) {
               jThrType["randoms"]["RNDM" + std::to_string(rc)]["cut"] = 
                  randoms->get().getAttribute<unsigned int>("RNDM" + std::to_string(rc)+".cut");
            }
         }
      }
 
      if(thrType == "MU") {
         auto & muinfo = l1menu.thrExtraInfo().MU();
         jThrType["exclusionLists"] = json::value_t::object;
         for(auto & listName : muinfo.exclusionListNames()) {
            jThrType["exclusionLists"][listName] = json::array_t({});
            for(auto & x : muinfo.exclusionList(listName)) {
               jThrType["exclusionLists"][listName] += json({ {"sectorName", x.first}, {"rois", x.second}});
            }
         }
         for(auto & rpcPt : muinfo.knownRpcPtValues()) {
            jThrType["roads"]["rpc"][std::to_string(rpcPt)] = muinfo.rpcIdxForPt(rpcPt);
         }
         for(auto & tgcPt : muinfo.knownTgcPtValues()) {
            jThrType["roads"]["tgc"][std::to_string(tgcPt)] = muinfo.tgcIdxForPt(tgcPt);
         }
      }
 
      if(thrType == "EM" && l1menu.run()>1) {
         auto & eminfo = l1menu.thrExtraInfo().EM();
         for(const std::string isoSet : { "EMIsoForEMthr", "HAIsoForEMthr" }) {
            jThrType["isolation"][isoSet]["thrtype"] = isoSet;
            jThrType["isolation"][isoSet]["Parametrization"] = json::array_t({});
            for(size_t bit = 1; bit<=5; ++bit) {
               auto & iso = eminfo.isolation(isoSet,bit);
               json jIso({});
               jIso["etamax"] = iso.etamax();
               jIso["etamin"] = iso.etamin();
               jIso["isobit"] = iso.isobit();
               jIso["mincut"] = iso.mincut();
               jIso["offset"] = iso.offset();
               jIso["priority"] = iso.priority();
               jIso["slope"] = iso.slope();
               jIso["upperlimit"] = iso.upperlimit();
               jThrType["isolation"][isoSet]["Parametrization"] += jIso;
            }
         }
      }
 
      if(thrType == "TAU" && l1menu.run()>1) {
         auto & tauinfo = l1menu.thrExtraInfo().TAU();
         const std::string isoSet{ "EMIsoForTAUthr" };
         jThrType["isolation"][isoSet]["thrtype"] = isoSet;
         jThrType["isolation"][isoSet]["Parametrization"] = json::array_t({});
         for(size_t bit = 1; bit<=5; ++bit) {
            auto & iso = tauinfo.isolation(isoSet,bit);
            json jIso({});
            jIso["etamax"] = iso.etamax();
            jIso["etamin"] = iso.etamin();
            jIso["isobit"] = iso.isobit();
            jIso["mincut"] = iso.mincut();
            jIso["offset"] = iso.offset();
            jIso["priority"] = iso.priority();
            jIso["slope"] = iso.slope();
            jIso["upperlimit"] = iso.upperlimit();
            jThrType["isolation"][isoSet]["Parametrization"] += jIso;
         }
      }
 
      if(thrType == "JET") {
         auto & jetinfo = l1menu.thrExtraInfo().JET();
         jThrType["ptMinToTopoSmallWindow"] = (int)jetinfo.ptMinToTopoSmallWindow();
         jThrType["ptMinToTopoLargeWindow"] = (int)jetinfo.ptMinToTopoLargeWindow();
      }
 
      if(thrType == "XS") {
         auto & xsinfo = l1menu.thrExtraInfo().XS();
         jThrType["significance"]["xeMin"] = xsinfo.xeMin();
         jThrType["significance"]["xeMax"] = xsinfo.xeMax();
         jThrType["significance"]["teSqrtMin"] = xsinfo.teSqrtMin();
         jThrType["significance"]["teSqrtMax"] = xsinfo.teSqrtMax();
         jThrType["significance"]["xsSigmaScale"] = xsinfo.xsSigmaScale();
         jThrType["significance"]["xsSigmaOffset"] = xsinfo.xsSigmaOffset();
      }
 
      if(thrType == "eEM") {
         auto & eeminfo = l1menu.thrExtraInfo().eEM();
         jThrType["maxEt"] = (int)eeminfo.maxEt();
         for( auto wp : {TrigConf::Selection::WP::LOOSE, TrigConf::Selection::WP::MEDIUM, TrigConf::Selection::WP::TIGHT} ) {
            const auto & wpstr = TrigConf::Selection::wpToString(wp);
            jThrType["workingPoints"][wpstr] = json::array_t({});
            for(auto & iso : eeminfo.isolation(wp)) {
               json jWPIso({});
               std::stringstream stream;
               stream << std::fixed << std::setprecision(3) << iso.value().reta_d();
               jWPIso["reta"] = std::stod(stream.str());
               jWPIso["reta_fw"] = iso.value().reta_fw();
               stream.str("");
               stream.clear();
               stream << std::fixed << std::setprecision(3) << iso.value().rhad_d();
               jWPIso["rhad"] = std::stod(stream.str());
               jWPIso["rhad_fw"] = iso.value().rhad_fw();
               stream.str(""); 
               stream.clear();
               stream << std::fixed << std::setprecision(3) << iso.value().wstot_d();               
               jWPIso["wstot"] = std::stod(stream.str());
               jWPIso["wstot_fw"] = iso.value().wstot_fw();
               jWPIso["etamin"] = iso.etaMin();
               jWPIso["etamax"] = iso.etaMax();
               jWPIso["priority"] = iso.priority();
               jThrType["workingPoints"][wpstr] += jWPIso;
            }
         }
      }
 
      if(thrType == "jEM") {
         auto & jeminfo = l1menu.thrExtraInfo().jEM();
         jThrType["maxEt"] = (int)jeminfo.maxEt();
         jThrType["ptMinToTopo1"] = (int)jeminfo.ptMinToTopo("1A");
         jThrType["ptMinToTopo2"] = (int)jeminfo.ptMinToTopo("2A");
         jThrType["ptMinToTopo3"] = (int)jeminfo.ptMinToTopo("3A");
         jThrType["ptMinxTOB1"] = (int)jeminfo.ptMinxTOB("1A");
         jThrType["ptMinxTOB2"] = (int)jeminfo.ptMinxTOB("2A");
         jThrType["ptMinxTOB3"] = (int)jeminfo.ptMinxTOB("3A");
         for( auto wp : {TrigConf::Selection::WP::LOOSE, TrigConf::Selection::WP::MEDIUM, TrigConf::Selection::WP::TIGHT} ) {
            const auto & wpstr = TrigConf::Selection::wpToString(wp);
            jThrType["workingPoints"][wpstr] = json::array_t({});
            for(auto & iso : jeminfo.isolation(wp)) {
               json jWPIso({});
               std::stringstream stream;
               stream << std::fixed << std::setprecision(3) << iso.value().iso_d();
               jWPIso["iso"] = std::stod(stream.str());
               jWPIso["iso_fw"] = iso.value().iso_fw();
               stream.str("");
               stream.clear();
               stream << std::fixed << std::setprecision(3) << iso.value().frac_d();
               jWPIso["frac"] = std::stod(stream.str());
               jWPIso["frac_fw"] = iso.value().frac_fw();
               stream.str("");
               stream.clear();
               stream << std::fixed << std::setprecision(3) << iso.value().frac2_d();
               jWPIso["frac2"] = std::stod(stream.str());
               jWPIso["frac2_fw"] = iso.value().frac2_fw();
               jWPIso["etamin"] = iso.etaMin();
               jWPIso["etamax"] = iso.etaMax();
               jWPIso["priority"] = iso.priority();
               jThrType["workingPoints"][wpstr] += jWPIso;
            }
         }
      }
 
      if(thrType == "eTAU") {
         auto & eeminfo = l1menu.thrExtraInfo().eTAU();
         jThrType["minIsoEt"] = (int)eeminfo.minIsoEt();
         jThrType["maxEt"] = (int)eeminfo.maxEt();
         jThrType["algoVersion"] = (int)eeminfo.algoVersion();
         for( auto wp : {TrigConf::Selection::WP::LOOSE, TrigConf::Selection::WP::MEDIUM, TrigConf::Selection::WP::TIGHT, 
                         TrigConf::Selection::WP::HADLOOSE, TrigConf::Selection::WP::HADMEDIUM, TrigConf::Selection::WP::HADTIGHT} ) {
            const auto & wpstr = TrigConf::Selection::wpToString(wp);
            jThrType["workingPoints"][wpstr] = json::array_t({});
            for(auto & iso : eeminfo.isolation(wp)) {
               json jWPIso({});
               std::stringstream stream;
               stream << std::fixed << std::setprecision(3) << iso.value().rCore_d();
               jWPIso["rCore"] = std::stod(stream.str());
               jWPIso["rCore_fw"] = iso.value().rCore_fw();
               stream.str("");
               stream.clear();
               stream << std::fixed << std::setprecision(3) << iso.value().rHad_d();
               jWPIso["rHad"] = std::stod(stream.str());
               jWPIso["rHad_fw"] = iso.value().rHad_fw();
               jThrType["workingPoints"][wpstr] += jWPIso;
            }
         }
      }
 
      if(thrType == "jTAU") {
         auto & jtauinfo = l1menu.thrExtraInfo().jTAU();
         jThrType["maxEt"] = (int)jtauinfo.maxEt();
         jThrType["ptMinToTopo1"] = (int)jtauinfo.ptMinToTopo("1A");
         jThrType["ptMinToTopo2"] = (int)jtauinfo.ptMinToTopo("2A");
         jThrType["ptMinToTopo3"] = (int)jtauinfo.ptMinToTopo("3A");
         jThrType["ptMinxTOB1"] = (int)jtauinfo.ptMinxTOB("1A");
         jThrType["ptMinxTOB2"] = (int)jtauinfo.ptMinxTOB("2A");
         jThrType["ptMinxTOB3"] = (int)jtauinfo.ptMinxTOB("3A");
         for( auto wp : {TrigConf::Selection::WP::LOOSE, TrigConf::Selection::WP::MEDIUM, TrigConf::Selection::WP::TIGHT} ) {
            const auto & wpstr = TrigConf::Selection::wpToString(wp);
            jThrType["workingPoints"][wpstr] = json::array_t({});
            for(auto & iso : jtauinfo.isolation(wp)) {
               json jWPIso({});
               std::stringstream stream;
               stream << std::fixed << std::setprecision(3) << iso.value().isolation_d();
               jWPIso["isolation"] = std::stod(stream.str());
               jWPIso["isolation_fw"] = iso.value().isolation_fw();
               jThrType["workingPoints"][wpstr] += jWPIso;
            }
         }
      }
 
      if(thrType == "cTAU") {
         auto & eeminfo = l1menu.thrExtraInfo().cTAU();
         for( auto wp : {TrigConf::Selection::WP::LOOSE, TrigConf::Selection::WP::MEDIUM, TrigConf::Selection::WP::TIGHT} ) {
            const auto & wpstr = TrigConf::Selection::wpToString(wp);
            jThrType["workingPoints"][wpstr] = json::array_t({});
            for(auto & iso : eeminfo.isolation(wp)) {
               json jWPIso({});
               std::stringstream stream;
               stream << std::fixed << std::setprecision(3) << iso.value().isolation_d();
               jWPIso["isolation"] = std::stod(stream.str());
               jWPIso["isolation_fw"] = iso.value().isolation_fw();
           stream.str("");
               stream << std::fixed << std::setprecision(3) << iso.value().isolation_jTAUCoreScale_d();
               jWPIso["isolation_jTAUCoreScale"] = std::stod(stream.str());
               jWPIso["isolation_jTAUCoreScale_fw"] = iso.value().isolation_jTAUCoreScale_fw();
           stream.str("");
               stream << std::fixed << std::setprecision(3) << iso.value().eTAU_rCoreMin_WP_d();
               jWPIso["eTAU_rCoreMin"] = std::stod(stream.str());
               jWPIso["eTAU_rCoreMin_WP_fw"] = iso.value().eTAU_rCoreMin_WP_fw();
           stream.str("");
               stream << std::fixed << std::setprecision(3) << iso.value().eTAU_rHadMin_WP_d();
               jWPIso["eTAU_rHadMin"] = std::stod(stream.str());
               jWPIso["eTAU_rHadMin_WP_fw"] = iso.value().eTAU_rHadMin_WP_fw();
               jThrType["workingPoints"][wpstr] += jWPIso;
            }
         }
      }
 
      if(thrType == "jJ") {
         auto & jjinfo = l1menu.thrExtraInfo().jJ();
         jThrType["ptMinToTopo1"] = (int)jjinfo.ptMinToTopo("1A");
         jThrType["ptMinToTopo2"] = (int)jjinfo.ptMinToTopo("2A");
         jThrType["ptMinToTopo3"] = (int)jjinfo.ptMinToTopo("3A");
         jThrType["ptMinxTOB1"] = (int)jjinfo.ptMinxTOB("1A");
         jThrType["ptMinxTOB2"] = (int)jjinfo.ptMinxTOB("2A");
         jThrType["ptMinxTOB3"] = (int)jjinfo.ptMinxTOB("3A");
      }
 
      if(thrType == "jLJ") {
         auto & jljinfo = l1menu.thrExtraInfo().jLJ();
         jThrType["ptMinToTopo1"] = (int)jljinfo.ptMinToTopo("1A");
         jThrType["ptMinToTopo2"] = (int)jljinfo.ptMinToTopo("2A");
         jThrType["ptMinToTopo3"] = (int)jljinfo.ptMinToTopo("3A");
         jThrType["ptMinxTOB1"] = (int)jljinfo.ptMinxTOB("1A");
         jThrType["ptMinxTOB2"] = (int)jljinfo.ptMinxTOB("2A");
         jThrType["ptMinxTOB3"] = (int)jljinfo.ptMinxTOB("3A");
      }
 
      if(thrType == "gJ") {
         auto & gjinfo = l1menu.thrExtraInfo().gJ();
         jThrType["ptMinToTopo1"] = (int)gjinfo.ptMinToTopo(1);
         jThrType["ptMinToTopo2"] = (int)gjinfo.ptMinToTopo(2);
      }
 
      if(thrType == "gLJ") {
         auto & gljinfo = l1menu.thrExtraInfo().gLJ();
         jThrType["ptMinToTopo1"] = (int)gljinfo.ptMinToTopo(1);
         jThrType["ptMinToTopo2"] = (int)gljinfo.ptMinToTopo(2);
         jThrType["seedThrA"] = (int)gljinfo.seedThr('A');
         jThrType["seedThrB"] = (int)gljinfo.seedThr('B');
         jThrType["seedThrC"] = (int)gljinfo.seedThr('C');
         std::stringstream stream;
         stream << std::fixed << std::setprecision(3) << gljinfo.rhoTowerMin('A');
         jThrType["rhoTowerMinA"] = std::stod(stream.str());
         stream.str("");
         stream.clear(); 
         stream << std::fixed << std::setprecision(3) << gljinfo.rhoTowerMin('B');
         jThrType["rhoTowerMinB"] = std::stod(stream.str());
         stream.str("");
         stream.clear();      
         stream << std::fixed << std::setprecision(3) << gljinfo.rhoTowerMin('C');
         jThrType["rhoTowerMinC"] = std::stod(stream.str());
         stream.str("");
         stream.clear(); 
         stream << std::fixed << std::setprecision(3) << gljinfo.rhoTowerMax('A');
         jThrType["rhoTowerMaxA"] = std::stod(stream.str());
         stream.str("");
         stream.clear();      
         stream << std::fixed << std::setprecision(3) << gljinfo.rhoTowerMax('B');
         jThrType["rhoTowerMaxB"] = std::stod(stream.str());
         stream.str("");
         stream.clear();
         stream << std::fixed << std::setprecision(3) << gljinfo.rhoTowerMax('C');
         jThrType["rhoTowerMaxC"] = std::stod(stream.str());
      }
 
      if(thrType == "jXE") {
          // nothing to do for now...
      }
 
      if(thrType == "jTE") {
         auto & ei = l1menu.thrExtraInfo().jTE();
         jThrType["etaBoundary1"] = ei.etaBoundary("1A");
         jThrType["etaBoundary1_fw"] = ei.etaBoundary_fw("1A");
         jThrType["etaBoundary2"] = ei.etaBoundary("2A");
         jThrType["etaBoundary2_fw"] = ei.etaBoundary_fw("2A");
         jThrType["etaBoundary3"] = ei.etaBoundary("3A");
         jThrType["etaBoundary3_fw"] = ei.etaBoundary_fw("3A");   
      }
 
      if(thrType == "gXE") {
         auto & ei = l1menu.thrExtraInfo().gXE();
         jThrType["seedThrA"] = (int)ei.seedThr('A');
         jThrType["seedThrB"] = (int)ei.seedThr('B');
         jThrType["seedThrC"] = (int)ei.seedThr('C');
         jThrType["XERHO_sigmaPosA"] = ei.XERHO_param('A',true);
         jThrType["XERHO_sigmaPosB"] = ei.XERHO_param('B',true);
         jThrType["XERHO_sigmaPosC"] = ei.XERHO_param('C',true);
         jThrType["XERHO_sigmaNegA"] = ei.XERHO_param('A',false);
         jThrType["XERHO_sigmaNegB"] = ei.XERHO_param('B',false);
         jThrType["XERHO_sigmaNegC"] = ei.XERHO_param('C',false);
         jThrType["XEJWOJ_a_A"] = ei.JWOJ_param('A','a');
         jThrType["XEJWOJ_a_B"] = ei.JWOJ_param('B','a');
         jThrType["XEJWOJ_a_C"] = ei.JWOJ_param('C','a');
         jThrType["XEJWOJ_b_A"] = ei.JWOJ_param('A','b');
         jThrType["XEJWOJ_b_B"] = ei.JWOJ_param('B','b');
         jThrType["XEJWOJ_b_C"] = ei.JWOJ_param('C','b');
         jThrType["XEJWOJ_c_A"] = ei.JWOJ_param('A','c');
         jThrType["XEJWOJ_c_B"] = ei.JWOJ_param('B','c');
         jThrType["XEJWOJ_c_C"] = ei.JWOJ_param('C','c');
         jThrType["XENOISECUT_noiseCutThrA"] = ei.noiseCutThr('A');
         jThrType["XENOISECUT_noiseCutThrB"] = ei.noiseCutThr('B');
         jThrType["XENOISECUT_noiseCutThrC"] = ei.noiseCutThr('C');
      }
 
      if(thrType == "gTE") {
          // nothing to do for now...
      }
 
      std::vector<std::string> legacyCalo = {"EM", "JET", "TAU", "XE", "TE", "XS", "ZB", "JB", "JF", "JE", "R2TOPO"};
      if( std::any_of(begin(legacyCalo), end(legacyCalo), [&thrType](const std::string &c) { return c==thrType; }) ) {
         thresholds["legacyCalo"][thrType] = std::move(jThrType);
      } else {
         thresholds[thrType] = std::move(jThrType);
      }
   };
 
   json boards = json::object_t{};
   for( auto & bname : l1menu.boardNames() ) {
      auto & bdef = l1menu.board(bname);
      boards[bname] = json{ {"connectors", bdef.connectorNames()}, {"type", bdef.type()} };
      if(bdef.legacy())
         boards[bname]["legacy"] = true;
   };
 
   json connectors = json::object_t{};
   for( auto & cname : l1menu.connectorNames() ) {
      auto jConn = json{};
      auto & cdef = l1menu.connector(cname);
      jConn["type"] = cdef.type();
      if(cdef.legacy())
         jConn["legacy"] = true;
      jConn["triggerlines"] = json::array_t();
      if(cdef.maxClock() == 2){
          if(cdef.maxFpga() == 2){
              // legacy topo, TOPO2, TOPO3 and muctpi
              for(size_t fpga = 0; fpga<cdef.maxFpga(); ++fpga) {
                  for(size_t clock = 0; clock<cdef.maxClock(); ++clock) {
                      for(auto & tl : cdef.triggerLines(fpga, clock)) {
                         jConn["triggerlines"] += json({ {"name", tl.name()}, {"nbits",tl.nbits()}, {"startbit", tl.startbit()}, {"flatindex", tl.flatindex()}, {"fpga", tl.fpga()}, {"clock", tl.clock()}, });
                      }
                  }
              }              
          } else {
              // AlfaCpti and merger board
              for(size_t fpga = 0; fpga<cdef.maxFpga(); ++fpga) {
                  for(size_t clock = 0; clock<cdef.maxClock(); ++clock) {
                      for(auto & tl : cdef.triggerLines(fpga, clock)) {
                         jConn["triggerlines"] += json({ {"name", tl.name()}, {"nbits",tl.nbits()}, {"startbit", tl.startbit()}, {"flatindex", tl.flatindex()}, {"clock", tl.clock()}, });
                      }
                  }
              }
          }
      } else {
          for(auto & tl : cdef.triggerLines()) {
             jConn["triggerlines"] += json({ {"name", tl.name()}, {"nbits",tl.nbits()}, {"startbit", tl.startbit()}, {"flatindex", tl.flatindex()} });
          }
      }
      connectors[cname] = std::move(jConn);
   }
 
   json ctp = json::object_t{};
   {
      for(size_t slot=7; slot<=9; ++slot) {
         std::string sName = "slot" + std::to_string(slot);
         ctp["inputs"]["ctpin"][sName] = json({});
         for(size_t conn=0; conn<=3; ++conn) {
            ctp["inputs"]["ctpin"][sName]["connector" + std::to_string(conn)] = l1menu.ctp().ctpin(slot,conn);
         }
      }
      ctp["inputs"]["electrical"] = json::object_t{};
      for(size_t conn=0; conn<3; ++conn) {
         if(l1menu.ctp().electrical(conn)=="")
            continue;
         ctp["inputs"]["electrical"]["connector" + std::to_string(conn)] = l1menu.ctp().electrical(conn);
      }
      ctp["inputs"]["optical"] = json::object_t{};
      for(size_t conn=0; conn<12; ++conn) {
         if(l1menu.ctp().optical(conn)=="")
            continue;
         ctp["inputs"]["optical"]["connector" + std::to_string(conn)] = l1menu.ctp().optical(conn);
      }
      json ctpmon = json::object_t{};
      for(auto & mon : l1menu.ctp().ctpMon()) {
         ctpmon[mon.first] = json({{"multiplicity",mon.second.first},{"thr",mon.second.second}});
      }
      json ctpin = json::object_t{};
      for(auto & mon : l1menu.ctp().ctpinMon()) {
         ctpin[mon.first] = json({{"multiplicity", mon.second.first},{"thr", mon.second.second}});
      }
      ctp["monitoring"] = json({{"ctpmon", ctpmon}, {"ctpin", ctpin}});
   }
 
   json jtopo = json::object_t{};
   if(l1menu.run()>1) {
      std::map<L1TopoAlgorithm::AlgorithmType,std::string> algTypeNames = {
         {L1TopoAlgorithm::AlgorithmType::SORTING, "sortingAlgorithms"},
         {L1TopoAlgorithm::AlgorithmType::DECISION, "decisionAlgorithms"},
         {L1TopoAlgorithm::AlgorithmType::MULTIPLICITY, "multiplicityAlgorithms"}
      };
      auto topoCategories = l1menu.isRun2() ? std::vector<std::string> {"R2TOPO"} : std::vector<std::string> {"TOPO", "MUTOPO", "MULTTOPO", "R2TOPO"};
      for( const std::string& topoCat : topoCategories ) {
         for(auto & algName : l1menu.topoAlgorithmNames(topoCat)) {
            json jalg = json::object_t({});
            auto & alg = l1menu.algorithm(algName,topoCat);
            jalg["klass"]     = alg.klass();
            // input
            if(alg.type()==L1TopoAlgorithm::AlgorithmType::MULTIPLICITY) {
               if(alg.inputs().size()>0) {
                  jalg["input"] = alg.inputs()[0];
               } else {
                  jalg["input"] = nullptr;
               }
            } else if(alg.type()==L1TopoAlgorithm::AlgorithmType::SORTING) {
               jalg["input"] = alg.inputs()[0];
            } else {
               jalg["input"] = alg.inputs();
            }
            // output
            if( alg.type() == L1TopoAlgorithm::AlgorithmType::DECISION ) {
               jalg["output"] = alg.outputs();
            } else {
               jalg["output"] = alg.outputs()[0];
            }
            // generic parameters
            if(topoCat == "MULTTOPO") {
               jalg["nbits"] = alg.getAttribute<unsigned int>("nbits");
               jalg["threshold"] = alg.getAttribute("threshold");
               if(alg.klass()=="EnergyThreshold"){
                  for(const L1TopoAlgorithm::VariableParameter & vpar : alg.parameters()) { 
                     if(vpar.name()=="flavour") jalg["flavour"] = XEFlavour::flavourIntToStr(vpar.value());
                  }
               }
            } else {
               const auto & ds = alg.generics();
               for(const auto & gpname : ds.getKeys()) {
                  auto gp = ds.getObject(gpname);
                  {
                     if(gp.hasAttribute("position")) {
                        jalg["fixedParameters"]["generics"][gpname]["position"] = gp.getAttribute<unsigned int>("position");
                     }
                     std::string pval = alg.genericParameter(gpname);
                     try{
                        int pvali = std::stoi(pval);
                        jalg["fixedParameters"]["generics"][gpname]["value"] = pvali;
                     } catch(std::invalid_argument &) {
                        jalg["fixedParameters"]["generics"][gpname]["value"] = pval;
                     }
                  }
               }
            }
            // variable parameters
            if(topoCat != "MULTTOPO") {
               jalg["variableParameters"] = json::array_t({});
               for(const L1TopoAlgorithm::VariableParameter & vpar : alg.parameters()) {
                  json jVPar({});
                  jVPar["name"] = vpar.name();
                  jVPar["value"] = vpar.value();
                  if(auto sel = vpar.selection_optional()) { jVPar["selection"] = *sel; }
                  jalg["variableParameters"] += jVPar;
               }
            }
            jtopo[topoCat][algTypeNames[alg.type()]][algName] = std::move(jalg);
         }
      }
   }
 
   json j = json::object_t{};
   j["filetype"] = "l1menu";
   j["run"] = l1menu.run();
   j["name"] = l1menu.name();
   j["items"] = std::move(items);
   j["thresholds"] = std::move(thresholds);
   j["topoAlgorithms"] = std::move(jtopo);
   j["boards"] = std::move(boards);
   j["connectors"] = std::move(connectors);
   j["ctp"] = std::move(ctp);
 
   std::ofstream outfile(filename);
   outfile << std::setw(4) << j << std::endl;
 
   TRG_MSG_INFO("Saved file " << filename);
   return true;
}

writeJsonFile() [3/3]

bool TrigConf::JsonFileWriterL1::writeJsonFile	(	const std::string &	filename,
		const L1PrescalesSet &	l1bgs ) const

Definition at line 890 of file JsonFileWriterL1.cxx.

                                                                                                                  {
   json j;
   j["filetype"] = "l1prescale";
   j["name"]  = l1ps.name();
   json cuts;
   for ( const auto & [itemName, ps]: l1ps.prescales()){
      json cut;
      cut["cut"] = ps.cut;
      cut["enabled"] = ps.enabled;
      cut["info"] = "prescale: " + std::to_string(ps.prescale);
      cuts[itemName] = std::move(cut);
   }
   j["cutValues"] = std::move(cuts);
   std::ofstream outfile(filename);
   outfile << std::setw(4) << j << std::endl;
   TRG_MSG_INFO("Saved file " << filename);
   return true;
}

Member Data Documentation

m_msg_tls

boost::thread_specific_ptr<MsgStreamTC> TrigConf::TrigConfMessaging::m_msg_tls

mutableprivateinherited

MsgStreamTC instance (a std::cout like with print-out levels)

Definition at line 71 of file TrigConfMessaging.h.

m_name

std::string TrigConf::TrigConfMessaging::m_name

privateinherited

Definition at line 72 of file TrigConfMessaging.h.

---

The documentation for this class was generated from the following files:

• JsonFileWriterL1.h
• JsonFileWriterL1.cxx