L1TriggerConfig/DTTPGConfig/src/DTConfigManager.cc - master - githubmirror/cmssw

L1TriggerConfig/DTTPGConfig/src/DTConfigManager.cc
255 строк9 KB
Cms Build
Clang-Format
22 май 2019, 04:56
22 май 2019, 04:565d82455
//----------------------------------------------------------------------
//
//   Class: DTConfigManager
//
//   Description: DT Configuration manager includes config classes for every single chip
//
//
//   Author List:
//   C.Battilana
//
//   april 07 : SV DTConfigTrigUnit added
//   april 07 : CB Removed DTGeometry dependecies
//-----------------------------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------
#include "L1TriggerConfig/DTTPGConfig/interface/DTConfigManager.h"

//---------------
// C++ Headers --
//---------------
#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"

using namespace std;

//----------------
// Constructors --
//----------------

DTConfigManager::DTConfigManager() {}

//--------------
// Destructor --
//--------------

DTConfigManager::~DTConfigManager() {
  my_sectcollmap.clear();
  my_trigunitmap.clear();
  my_tsphimap.clear();
  my_tsthetamap.clear();
  my_tracomap.clear();
  my_btimap.clear();
}

//--------------
// Operations --
//--------------

const DTConfigBti* DTConfigManager::getDTConfigBti(DTBtiId btiid) const {
  DTChamberId chambid = btiid.SLId().chamberId();
  BtiMap::const_iterator biter1 = my_btimap.find(chambid);
  if (biter1 == my_btimap.end()) {
    std::cout << "DTConfigManager::getConfigBti : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
    return nullptr;
  }

  innerBtiMap::const_iterator biter2 = (*biter1).second.find(btiid);
  if (biter2 == (*biter1).second.end()) {
    std::cout << "DTConfigManager::getConfigBti : BTI (" << btiid.wheel() << "," << btiid.sector() << ","
              << btiid.station() << "," << btiid.superlayer() << "," << btiid.bti() << ") not found, return 0"
              << std::endl;
    return nullptr;
  }
  return &(*biter2).second;
}

const std::map<DTBtiId, DTConfigBti>& DTConfigManager::getDTConfigBtiMap(DTChamberId chambid) const {
  BtiMap::const_iterator biter = my_btimap.find(chambid);
  if (biter == my_btimap.end()) {
    std::cout << "DTConfigManager::getConfigBtiMap : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return a reference to the end of the map" << std::endl;
  }

  return (*biter).second;
}

const DTConfigTraco* DTConfigManager::getDTConfigTraco(DTTracoId tracoid) const {
  DTChamberId chambid = tracoid.ChamberId();
  TracoMap::const_iterator titer1 = my_tracomap.find(chambid);
  if (titer1 == my_tracomap.end()) {
    std::cout << "DTConfigManager::getConfigTraco : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
    return nullptr;
  }

  innerTracoMap::const_iterator titer2 = (*titer1).second.find(tracoid);
  if (titer2 == (*titer1).second.end()) {
    std::cout << "DTConfigManager::getConfigTraco : TRACO (" << tracoid.wheel() << "," << tracoid.sector() << ","
              << tracoid.station() << "," << tracoid.traco() << ") not found, return a reference to the end of the map"
              << std::endl;
    return nullptr;
  }
  return &(*titer2).second;
}

const std::map<DTTracoId, DTConfigTraco>& DTConfigManager::getDTConfigTracoMap(DTChamberId chambid) const {
  TracoMap::const_iterator titer = my_tracomap.find(chambid);
  if (titer == my_tracomap.end()) {
    std::cout << "DTConfigManager::getConfigTracoMap : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
  }

  return (*titer).second;
}

const DTConfigTSTheta* DTConfigManager::getDTConfigTSTheta(DTChamberId chambid) const {
  TSThetaMap::const_iterator thiter = my_tsthetamap.find(chambid);
  if (thiter == my_tsthetamap.end()) {
    std::cout << "DTConfigManager::getConfigTSTheta : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
    return nullptr;
  }

  return &(*thiter).second;
}

const DTConfigTSPhi* DTConfigManager::getDTConfigTSPhi(DTChamberId chambid) const {
  TSPhiMap::const_iterator phiter = my_tsphimap.find(chambid);
  if (phiter == my_tsphimap.end()) {
    std::cout << "DTConfigManager::getConfigTSPhi : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
    return nullptr;
  }

  return &(*phiter).second;
}

const DTConfigTrigUnit* DTConfigManager::getDTConfigTrigUnit(DTChamberId chambid) const {
  TrigUnitMap::const_iterator tuiter = my_trigunitmap.find(chambid);
  if (tuiter == my_trigunitmap.end()) {
    std::cout << "DTConfigManager::getConfigTrigUnit : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
    return nullptr;
  }

  return &(*tuiter).second;
}

const DTConfigLUTs* DTConfigManager::getDTConfigLUTs(DTChamberId chambid) const {
  LUTMap::const_iterator lutiter = my_lutmap.find(chambid);
  if (lutiter == my_lutmap.end()) {
    std::cout << "DTConfigManager::getConfigLUTs : Chamber (" << chambid.wheel() << "," << chambid.sector() << ","
              << chambid.station() << ") not found, return 0" << std::endl;
    return nullptr;
  }

  return &(*lutiter).second;
}

const DTConfigSectColl* DTConfigManager::getDTConfigSectColl(DTSectCollId scid) const {
  SectCollMap::const_iterator sciter = my_sectcollmap.find(scid);
  if (sciter == my_sectcollmap.end()) {
    std::cout << "DTConfigManager::getConfigSectColl : SectorCollector (" << scid.wheel() << "," << scid.sector()
              << ") not found, return 0" << std::endl;
    return nullptr;
  }

  return &(*sciter).second;
}

const DTConfigPedestals* DTConfigManager::getDTConfigPedestals() const { return &my_pedestals; }

int DTConfigManager::getBXOffset() const {
  int ST = static_cast<int>(getDTConfigBti(DTBtiId(1, 1, 1, 1, 1))->ST());
  return (ST / 2 + ST % 2);
}

void DTConfigManager::setDTConfigBti(DTBtiId btiid, DTConfigBti conf) {
  DTChamberId chambid = btiid.SLId().chamberId();
  my_btimap[chambid][btiid] = conf;
}

void DTConfigManager::setDTConfigTraco(DTTracoId tracoid, DTConfigTraco conf) {
  DTChamberId chambid = tracoid.ChamberId();
  my_tracomap[chambid][tracoid] = conf;
}

void DTConfigManager::dumpLUTParam(DTChamberId& chambid) const {
  // open txt file
  string name = "Lut_from_param";
  name += ".txt";

  ofstream fout;
  fout.open(name.c_str(), ofstream::app);

  // get wheel, station, sector from chamber
  int wh = chambid.wheel();
  int st = chambid.station();
  int se = chambid.sector();

  //cout << "Dumping lut command for wh " << wh << " st " << st << " se " << se << endl;

  fout << wh;
  fout << "\t" << st;
  fout << "\t" << se;

  // get parameters from configuration
  // get DTConfigLUTs for this chamber
  const DTConfigLUTs* _confLUTs = getDTConfigLUTs(chambid);
  short int btic = getDTConfigTraco(DTTracoId(wh, st, se, 1))->BTIC();
  float d = _confLUTs->D();
  float xcn = _confLUTs->Xcn();
  //fout << "\td\t" << d << "\txcn\t" << xcn << "\t";
  //fout << "btic\t" << btic << "\t";

  // *** dump TRACO LUT command
  fout << "\tA8";
  short int Low_byte = (btic & 0x00FF);  // output in hex bytes format with zero padding
  short int High_byte = (btic >> 8 & 0x00FF);
  fout << setw(2) << setfill('0') << hex << High_byte << setw(2) << setfill('0') << Low_byte;

  // convert parameters from IEE32 float to DSP float format
  short int DSPmantissa = 0;
  short int DSPexp = 0;

  // d parameter conversion and dump
  _confLUTs->IEEE32toDSP(d, DSPmantissa, DSPexp);
  Low_byte = (DSPmantissa & 0x00FF);  // output in hex bytes format with zero padding
  High_byte = (DSPmantissa >> 8 & 0x00FF);
  fout << setw(2) << setfill('0') << hex << High_byte << setw(2) << setfill('0') << Low_byte;
  Low_byte = (DSPexp & 0x00FF);
  High_byte = (DSPexp >> 8 & 0x00FF);
  fout << setw(2) << setfill('0') << High_byte << setw(2) << setfill('0') << Low_byte;

  // xnc parameter conversion and dump
  DSPmantissa = 0;
  DSPexp = 0;
  _confLUTs->IEEE32toDSP(xcn, DSPmantissa, DSPexp);
  Low_byte = (DSPmantissa & 0x00FF);  // output in hex bytes format with zero padding
  High_byte = (DSPmantissa >> 8 & 0x00FF);
  fout << setw(2) << setfill('0') << hex << High_byte << setw(2) << setfill('0') << Low_byte;
  Low_byte = (DSPexp & 0x00FF);
  High_byte = (DSPexp >> 8 & 0x00FF);
  fout << setw(2) << setfill('0') << High_byte << setw(2) << setfill('0') << Low_byte;

  // sign bits
  short int xcn_sign = _confLUTs->Wheel();
  Low_byte = (xcn_sign & 0x00FF);  // output in hex bytes format with zero padding
  High_byte = (xcn_sign >> 8 & 0x00FF);
  fout << setw(2) << setfill('0') << hex << High_byte << setw(2) << setfill('0') << Low_byte << dec << "\n";

  fout.close();

  return;
}