/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/SimCalorimetry/EcalSelectiveReadoutProducers/src/EcalSRCondTools.cc

Go to the documentation of this file.

//emacs settings:-*- mode: c++; c-basic-offset: 2; indent-tabs-mode: nil -*-
/*
 * $Id: EcalSRCondTools.cc,v 1.4 2012/06/26 18:53:24 eulisse Exp $
 *
 * author: Ph Gras. June, 2010
 */

#include "SimCalorimetry/EcalSelectiveReadoutProducers/interface/EcalSRCondTools.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"

#include "CondFormats/EcalObjects/interface/EcalSRSettings.h"
#include "CondFormats/DataRecord/interface/EcalSRSettingsRcd.h"
#include "CondFormats/EcalObjects/interface/EcalTPGPhysicsConst.h"
#include "CondFormats/DataRecord/interface/EcalTPGPhysicsConstRcd.h"

#include <string>
#include <fstream>
#include <algorithm>

constexpr int tccNum[12][12] = {
  /* EE- */
  {  36, 19, 20, 21, 22, 23, 18,  1,  2,  3,  4,  5}, //SRP 1
  {  24, 25, 26, 27, 28, 29,  6,  7,  8,  9, 10, 11}, //SRP 2
  {  30, 31, 32, 33, 34, 35, 12, 13, 14, 15, 16, 17}, //SRP 3
  /* EB- */
  { 54, 37, 38, 39, 40, 41, -1, -1, -1, -1, -1, -1},  //SRP 4
  { 42, 43, 44, 45, 46, 47, -1, -1, -1, -1, -1, -1},  //SRP 5
  { 48, 49, 50, 51, 52, 53, -1, -1, -1, -1, -1, -1},  //SRP 6
  /* EB+ */
  { 72, 55, 56, 57, 58, 59, -1, -1, -1, -1, -1, -1},  //SRP 7
  { 60, 61, 62, 63, 64, 65, -1, -1, -1, -1, -1, -1},  //SRP 8
  { 66, 67, 68, 69, 70, 71, -1, -1, -1, -1, -1, -1},  //SRP 9
  /* EE+ */
  {  90, 73, 74, 75, 76, 77,108, 91, 92, 93, 94, 95}, //SRP 10
  {  78, 79, 80, 81, 82, 83, 96, 97, 98, 99,100,101}, //SRP 11
  {  84, 85, 86, 87, 88, 89,102,103,104,105,106,107}   //SRP 12
};

constexpr int dccNum[12][12] = {
  {  1,  2,  3, -1, -1, -1, -1, -1, -1, -1, -1, -1},  //SRP 1 
  {  4,  5,  6, -1, -1, -1, -1, -1, -1, -1, -1, -1},  //SRP 2 
  {  7,  8,  9, -1, -1, -1, -1, -1, -1, -1, -1, -1},  //SRP 3 
  { 10, 11, 12, 13, 14, 15, -1, -1, -1, -1, -1, -1},  //SRP 4 
  { 16, 17, 18, 19, 20, 21, -1, -1, -1, -1, -1, -1},  //SRP 5 
  { 22, 23, 24, 25, 26, 27, -1, -1, -1, -1, -1, -1},  //SRP 6 
  { 28, 29, 30, 31, 32, 33, -1, -1, -1, -1, -1, -1},  //SRP 7 
  { 34, 35, 36, 37, 38, 39, -1, -1, -1, -1, -1, -1},  //SRP 8 
  { 40, 41, 42, 43, 44, 45, -1, -1, -1, -1, -1, -1},  //SRP 9 
  { 46, 47, 48, -1, -1, -1, -1, -1, -1, -1, -1, -1},  //SRP 10
  { 49, 50, 51, -1, -1, -1, -1, -1, -1, -1, -1, -1},  //SRP 11
  { 52, 53, 54, -1, -1, -1, -1, -1, -1, -1, -1, -1}  //SRP 12
};


using namespace std;

EcalSRCondTools::EcalSRCondTools(const edm::ParameterSet& ps):ps_(ps){
}


EcalSRCondTools::~EcalSRCondTools(){
}

void
EcalSRCondTools::analyze(const edm::Event& event, const edm::EventSetup& es){
  static bool done = false;
  if(done) return;
  EcalSRSettings* sr = new EcalSRSettings;

  string mode = ps_.getParameter<string>("mode");
  bool iomode;
  const bool iomode_read = false;
  const bool iomode_write = true;

  if(mode == "online_config" || mode == "combine_config"){
    iomode =iomode_write;
    string fname = ps_.getParameter<string>("onlineSrpConfigFile");
    ifstream f(fname.c_str());
    if(!f.good()){
      throw cms::Exception("EcalSRCondTools") << "Failed to open file " << fname;
    }
    importSrpConfigFile(*sr, f, true);
  }

  if(mode=="python_config" || mode=="combine_config"){
    iomode = iomode_write;
    importParameterSet(*sr, ps_);
  }

  if(mode=="read"){
    iomode = iomode_read;
  }

  if(!(mode=="python_config" || mode == "online_config" || mode == "combine_config" || (mode=="read"))){
    throw cms::Exception("Config") << "Invalid value," << mode << ",  for parameter mode. "
                                   << "Valid values: online_config, python_config, combine_config, read";
  }

  if(iomode==iomode_write){
    sr->bxGlobalOffset_ = ps_.getParameter<int>("bxGlobalOffset");
    sr->automaticSrpSelect_ = ps_.getParameter<int>("automaticSrpSelect");
    sr->automaticMasks_ = ps_.getParameter<int>("automaticMasks");
    
    edm::Service<cond::service::PoolDBOutputService> db;
    if( !db.isAvailable() ){
      throw cms::Exception("CondDBAccess") << "Failed to connect to PoolDBOutputService\n";
    }
    //fillup DB
    //create new infinite IOV
    cond::Time_t firstSinceTime = db->beginOfTime();
    db->writeOne(sr,firstSinceTime,"EcalSRSettingsRcd");
    done = true;
  } else {//read mode
    edm::ESHandle<EcalSRSettings> hSr;
    es.get<EcalSRSettingsRcd>().get(hSr);
    if(!hSr.isValid()){
      cout << "EcalSRSettings record not found. Check the Cond DB Global tag.\n";
    } else{
      const EcalSRSettings* sr = hSr.product();
      cout << "ECAL Seletive readout settings:\n";
      cout << *sr << "\n" << endl;
    }
    
    //trigger tower thresholds (from FENIX configuration):
    edm::ESHandle<EcalTPGPhysicsConst> hTp;
    es.get<EcalTPGPhysicsConstRcd>().get(hTp);
    if(!hTp.isValid()){
      cout << "EcalTPGPhysicsConst record not found. Check the Cond DB Global tag.\n";
    } else{
      const EcalTPGPhysicsConst * tp = hTp.product();
      const EcalTPGPhysicsConstMap mymap = tp->getMap();
      if(mymap.size()!=2){
        cout << "Error: TPG physics record is of unexpected size: "
             << mymap.size()  << " elements instead of two (one for EB, one for EE)\n";
      } else{
        EcalTPGPhysicsConstMap::const_iterator it=mymap.begin();
        cout << "----------------------------------------------------------------------\n"
          "Trigger tower Et thresholds extracted from TPG configuration \n"
          "(EcalSRCondTools modules supports only read mode for these parameters):\n\n";
        cout<< "EB: " << "LT = " << it->second.ttf_threshold_Low << " GeV "
            << "HT = "<< it->second.ttf_threshold_High << " GeV\n";
        ++it;
        cout<< "EE: " << "LT = " << it->second.ttf_threshold_Low << " GeV "
            << "HT = "<< it->second.ttf_threshold_High << " GeV\n";
        
      }
    }
  }
}

void EcalSRCondTools::importParameterSet(EcalSRSettings& sr, const edm::ParameterSet& ps){
  sr.deltaPhi_.resize(1);
  sr.deltaPhi_[0] = ps.getParameter<int >("deltaPhi");
  sr.deltaEta_.resize(1);
  sr.deltaEta_[0] = ps.getParameter<int >("deltaEta");
  sr.ecalDccZs1stSample_.resize(1);
  sr.ecalDccZs1stSample_[0] = ps.getParameter<int >("ecalDccZs1stSample");
  sr.ebDccAdcToGeV_ = ps.getParameter<double >("ebDccAdcToGeV");
  sr.eeDccAdcToGeV_ = ps.getParameter<double >("eeDccAdcToGeV");
  sr.dccNormalizedWeights_.resize(1);
  const std::vector<double>& w = ps.getParameter<std::vector<double>  >("dccNormalizedWeights");
  sr.dccNormalizedWeights_[0].resize(w.size());
  std::copy(w.begin(), w.end(), sr.dccNormalizedWeights_[0].begin());
  sr.symetricZS_.resize(1);
  sr.symetricZS_[0] = ps.getParameter<bool >("symetricZS");
  sr.srpLowInterestChannelZS_.resize(2);
  const int eb = 0;
  const int ee = 1;
  sr.srpLowInterestChannelZS_[eb] = ps.getParameter<double >("srpBarrelLowInterestChannelZS");
  sr.srpLowInterestChannelZS_[ee] = ps.getParameter<double >("srpEndcapLowInterestChannelZS");
  sr.srpHighInterestChannelZS_.resize(2);
  sr.srpHighInterestChannelZS_[eb] = ps.getParameter<double >("srpBarrelHighInterestChannelZS");
  sr.srpHighInterestChannelZS_[ee] = ps.getParameter<double >("srpEndcapHighInterestChannelZS");
  //sr.trigPrimBypass_.resize(1);
  //sr.trigPrimBypass_[0] = ps.getParameter<bool >("trigPrimBypass");
  //sr.trigPrimBypassMode_.resize(1);
  //sr.trigPrimBypassMode_[0] = ps.getParameter<int >("trigPrimBypassMode");
  //sr.trigPrimBypassLTH_.resize(1);
  //sr.trigPrimBypassLTH_[0] = ps.getParameter<double >("trigPrimBypassLTH");
  //sr.trigPrimBypassHTH_.resize(1);
  //sr.trigPrimBypassHTH_[0] = ps.getParameter<double >("trigPrimBypassHTH");
  //sr.trigPrimBypassWithPeakFinder_.resize(1);
  //sr.trigPrimBypassWithPeakFinder_[0] = ps.getParameter<bool >("trigPrimBypassWithPeakFinder");
  //sr.defaultTtf_.resize(1);
  //sr.defaultTtf_[0] = ps.getParameter<int >("defaultTtf");
  sr.actions_ = ps.getParameter<std::vector<int> >("actions");
}

void EcalSRCondTools::importSrpConfigFile(EcalSRSettings& sr, std::istream& f, bool d){
  //initialize vectors:
  sr.deltaEta_ = vector<int>(1,0);
  sr.deltaPhi_ = vector<int>(1,0);
  sr.actions_ = vector<int>(4, 0);
  sr.tccMasksFromConfig_ = vector<short>(EcalSRSettings::nTccs_, 0);
  sr.srpMasksFromConfig_ = vector<vector<short> >(EcalSRSettings::nSrps_, vector<short>(8, 0));
  sr.dccMasks_ = vector<short>(EcalSRSettings::nDccs_);
  sr.srfMasks_ = vector<short>(EcalSRSettings::nDccs_);
  sr.substitutionSrfs_ = vector<vector<short> >(EcalSRSettings::nSrps_, vector<short>(68,0));
  sr.testerTccEmuSrpIds_ = vector<int>(EcalSRSettings::nSrps_, 0);
  sr.testerSrpEmuSrpIds_ = vector<int>(EcalSRSettings::nSrps_, 0);
  sr.testerDccTestSrpIds_ = vector<int>(EcalSRSettings::nSrps_, 0);
  sr.testerSrpTestSrpIds_ = vector<int>(EcalSRSettings::nSrps_, 0);
  sr.bxOffsets_ = vector<short>(EcalSRSettings::nSrps_, 0);
  sr.automaticMasks_ = 0;
  sr.automaticSrpSelect_ = 0;
  
  //string line;
  int iLine = 0;
  int iValueSet = -1;
  const int nValueSets = 6*EcalSRSettings::nSrps_+9;
  string line;
  stringstream sErr("");
  while(!f.eof() && sErr.str().empty()){
    getline(f, line);
    ++iLine;
    line = trim(line);
    if(line[0] == '#' || line.empty()){//comment line and empty line to ignore
      continue;
    } else{
      ++iValueSet;
    }
    if(iValueSet>=nValueSets) break;
    uint32_t value;
    string sValue;
    int pos = 0;
    int iCh = 0;
    int nChs[nValueSets] = {
      //TCC masks: 0-11
      12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
      //SRP masks: 12-23
      8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,  8,
      //DCC masks: 24-35
      12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
      //SRF Masks: 36-47
      6,  6,  6,  6,  6,  6,  6,  6,  6,  6,  6,  6,
      //substitution SRFs: 48-59
      68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
      //Tester card to emulate or test: 60-71
      4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,  4,
      //Bx offsets: 72
      12,
      //algo type: 73
      1,
      //action flags: 74
      4,
      //pattern file directory: 75
      1,
      //VME slots: 76
      12,
      //card types: 77
      12,
      //config Mode
      1,
      //VME Interface card
      1,
      //Spy Mode
      12,
    };

    while(((sValue = tokenize(line, " \t", pos))!=string(""))
          && (iCh<nChs[iValueSet]) && sErr.str().empty()){
      value = strtoul(sValue.c_str(), 0, 0);
      const int iSrp = iValueSet%EcalSRSettings::nSrps_;
      if(iValueSet<12){//TCC
        assert((unsigned)iSrp < sizeof(tccNum) / sizeof(tccNum[0]));
        assert((unsigned)iCh < sizeof(tccNum[0]) / sizeof(tccNum[0][0]));
        int tcc = tccNum[iSrp][iCh];
        if(tcc>=0) {
          if(d) cout << "tccMasksFromConfig_[" << tcc << "] <- "
                     << value << "\n";
          sr.tccMasksFromConfig_[tcc-1] = value;
        }
      } else if(iValueSet<24){//SRP-SRP
        if(d) cout << "srpMasks_[" << iSrp << "][" << iCh << "] <- "
                   << value << "\n";
        sr.srpMasksFromConfig_[iSrp][iCh] = value;
      } else if(iValueSet<36) {//DCC output
        assert((unsigned)iSrp < sizeof(dccNum) / sizeof(dccNum[0]));
        assert((unsigned)iCh < sizeof(dccNum[0]) / sizeof(dccNum[0][0]));
        int dcc = dccNum[iSrp][iCh];
        if(dcc > 0){
          assert((unsigned)(dcc-1) < sr.dccMasks_.size());
          if(d) cout << "dccMasks_[" << (dcc-1) << "] <- "
                     << value << "\n";
          sr.dccMasks_[dcc-1] = value;
        }
      } else if(iValueSet<48){//SRF masks
        assert((unsigned)iSrp < sizeof(dccNum) / sizeof(dccNum[0]));
        assert((unsigned)iCh < sizeof(dccNum[0]) / sizeof(dccNum[0][0]));
        int dcc = dccNum[iSrp][iCh];
        if(dcc > 0){
          if(d) cout << "srfMasks_[" << (dcc-1) << "] <- "
                     << value << "\n";
          assert((unsigned)(dcc-1) < sr.srfMasks_.size());
          sr.srfMasks_[dcc-1] = value;
        }
      } else if(iValueSet<60){//substiution SRFs
        assert((unsigned)iSrp < sr.substitutionSrfs_.size());
        assert((unsigned)iCh <  sr.substitutionSrfs_[0].size());
        if(d) cout << "substitutionMasks_[" << iSrp << "][" << iCh << "] <- "
                   << value << "\n";
        sr.substitutionSrfs_[iSrp][iCh] = value;
      } else if(iValueSet<72){//Tester card config
        switch(iCh){
        case 0:
          assert((unsigned)iSrp < sr.testerTccEmuSrpIds_.size());
          if(d) cout << "testerTccEmuSrpIds_[" << iSrp << "] <- "
                     << value << "\n";
          sr.testerTccEmuSrpIds_[iSrp] = value;
          break;
        case 1:
          assert((unsigned)iSrp < sr.testerSrpEmuSrpIds_.size());
          if(d) cout << "testerSrpEmuSrpIds_[" << iSrp << "] <- "
                     << value << "\n";
          sr.testerSrpEmuSrpIds_[iSrp] = value;
          break;
        case 2:
          assert((unsigned)iSrp < sr.testerDccTestSrpIds_.size());
          if(d) cout << "testerDccTestSrpIds_[" << iSrp << "] <- "
                     << value << "\n";
          sr.testerDccTestSrpIds_[iSrp] = value;
          break;
        case 3:
          assert((unsigned)iSrp < sr.testerSrpTestSrpIds_.size());
          if(d) cout << "testerSrpTestSrpIds_[" << iSrp << "] <- "
                     << value << "\n";
          sr.testerSrpTestSrpIds_[iSrp] = value;
          break;
        default:
          sErr << "Syntax error in SRP system configuration "
               << " line " << iLine << ".";
        }
      } else if(iValueSet<73){//bx offsets
        assert((unsigned)iCh < sr.bxOffsets_.size());
        if(d) cout << "bxOffset_[" << iCh << "] <- "
                   << value << "\n";
        sr.bxOffsets_[iCh] = value;
      } else if(iValueSet<74){//algo type
        int algo = value;
        switch(algo){
        case 0:
          sr.deltaEta_[0] = sr.deltaPhi_[0] = 1;
          break;
        case 1:
          sr.deltaEta_[0] = sr.deltaPhi_[0]  = 2;
          break;
        default:
          throw cms::Exception("OutOfRange") << "Value of parameter algo ," << algo
                                             << ", is invalid. Valid values are 0 and 1.";
        }
        if(d) cout << "deltaEta_[0] <- " << sr.deltaEta_[0] << "\t"
                   << "deltaPhi_[0] <- " << sr.deltaPhi_[0] << "\n";
      } else if(iValueSet<75){//action flags
        assert((unsigned)iCh < sr.actions_.size());
        if(d) cout << "actions_[" << iCh << "] <- "
                   << value << "\n";
        sr.actions_[iCh] = value;
      } else if(iValueSet<76){//pattern file directory
//      emuDir_ = sValue;
//      if(d) cout << "emuDir_ <= "
//                 << value << "\n";
      } else if(iValueSet<77){//VME slots
        //      slotIds_[iCh] = value;
  //    if(d) cout << "slotIds_[" << iCh << "] <= "
        //                 << value << "\n";
      } else if(iValueSet<78){//card types
        //      cardTypes_[iCh] = sValue[0];
        //      if(d) cout << "cardTypes_[" << iCh << "] <= "
        //                 << value << "\n";
      } else if (iValueSet<79){//config mode
        //TODO validity check on value
        //      configMode_ = (ConfigMode)value;
        //      if(d) cout << "config mode <= " << value << "\n";
      } else if (iValueSet<80){//VME I/F
        //TODO validity check on value
        //     vmeInterface_ = (Vme::type_t)value;
        //if(d) cout << "Vme Interface code <= " << value << "\n";
      } else if (iValueSet<81){//Spy Mode
        //TODO validity check on value
        //      spyMode_[iCh] = value & 0x3;
        //      if(d) cout << "Spy mode <= " << value << "\n";
      } else{//should never be reached!
        assert(false);
      }
      ++iCh;
    }
    if(iCh!=nChs[iValueSet]){//error
      sErr << "Syntax error in imported SRP system configuration file "
        /*<< filename <<*/ " line " << iLine << ".";
    }
  }
  if(sErr.str().empty() && iValueSet!=(nValueSets-1)){//error
    sErr << "Syntax Error in imported SRP system configuration file "
      /*<< filename <<*/ " line " << iLine << ".";
  }
  if(sErr.str().size()!=0) throw cms::Exception("SyntaxError") << sErr.str();
}

double EcalSRCondTools::normalizeWeights(int hwWeight){
  //Fix sign bit in case only the 12 least significant bits of hwWeight were set
  //(hardware reprensentation uses only 12 bits)
  if(hwWeight & (1<<11)) hwWeight |= ~0xEFF;
  return hwWeight/1024.;
}

string EcalSRCondTools::tokenize(const string& s, const string& delim, int& pos){
  if(pos<0) return "";
  int pos0 = pos;
  int len = s.size();
  //eats delimeters at beginning of the string
  while(pos0<len && find(delim.begin(), delim.end(), s[pos0])!=delim.end()){
    ++pos0;
  }
  if(pos0==len) return "";
  pos = s.find_first_of(delim, pos0);
  return s.substr(pos0, (pos>0?pos:len)-pos0);
}

std::string EcalSRCondTools::trim(std::string s){
  std::string::size_type pos0 = s.find_first_not_of(" \t");
  if(pos0==string::npos){
    pos0=0;
  }
  string::size_type pos1 = s.find_last_not_of(" \t") + 1;
  if(pos1==string::npos){
    pos1 = pos0;
  }
  return s.substr(pos0, pos1-pos0);
}