LUCID_DigiSettings.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LUCID_DigiTools/LUCID_DigiSettings.h"
#include "CLHEP/Units/PhysicalConstants.h"
#include "CLHEP/Units/SystemOfUnits.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/MsgStream.h"
#include <iostream>
#include <utility>
 
 
LUCID_DigiSettings::LUCID_DigiSettings(): 
  m_msgSvc(nullptr),
  m_parValueNotSetByUserDouble(10.0e39),
  m_parValueNotSetByUserInt   (32767) {
  
  SetDefaultValues();
}
 
void LUCID_DigiSettings::Initialize(IMessageSvc* msgSvc) {
  
  m_msgSvc = msgSvc;
  
  MsgStream log(m_msgSvc, "LUCID_DigiSettings::initialize");
  
  log << MSG::INFO << "Initializing" << endmsg;
  
  OverwriteDigiParValues();
}
 
void LUCID_DigiSettings::SetDefaultValues() {
  
  MsgStream log(m_msgSvc, "LUCID_DigiSettings::SetDefaultValues");
 
  DefNewParameterInt   ("Number of LUCID tubes"          , "numTubes"         , &m_numTubes         , 1,  100);
  DefNewParameterInt   ("QDC channels per photo-electron", "qdcChannelsPerPE" , &m_qdcChannelsPerPE , 1 , 100);
  DefNewParameterDouble("QDC FED NoiseFactor"            , "qdcFedNoiseFactor", &m_qdcFedNoiseFactor, 1., 50.);
  DefNewParameterDouble("TDC PMT NoiseFactor"            , "tdcPmtNoiseFactor", &m_tdcPmtNoiseFactor, 1., 50.);
  DefNewParameterDouble("TDC FED NoiseFactor"            , "tdcFedNoiseFactor", &m_tdcFedNoiseFactor, 1., 50.);
  
  m_numTubes          = 40;
  m_qdcChannelsPerPE  = 10;
  m_qdcFedNoiseFactor = 3.0;
  m_tdcPmtNoiseFactor = 2.8;
  m_tdcFedNoiseFactor = 0.8;
}
 
void LUCID_DigiSettings::OverwriteDigiParProperties(Gaudi::Algorithm* alg) {
  
  std::map <std::string, parDouble>::iterator doubleMapIt  = m_doubleMap.begin();
  std::map <std::string, parInt>::iterator    intMapIt     = m_intMap.begin();
  
  for (; doubleMapIt != m_doubleMap.end(); ++doubleMapIt) alg->declareProperty(std::string("Overwrite_") + doubleMapIt->first, doubleMapIt->second.overwriteVal); 
  for (; intMapIt    != m_intMap.end()   ; ++intMapIt   ) alg->declareProperty(std::string("Overwrite_") + intMapIt->first   , intMapIt->second.overwriteVal);   
}
 
void LUCID_DigiSettings::OverwriteDigiParValues() {
  
  MsgStream log(m_msgSvc, "LUCID_DigiSettings");
 
  log << MSG::INFO << " OverwriteDigiParValues " << endmsg;
  
  std::map <std::string, parDouble>::iterator doubleMapIt = m_doubleMap.begin();
  std::map <std::string, parInt>::iterator    intMapIt    = m_intMap.begin();
 
  if (doubleMapIt->second.overwriteVal != m_parValueNotSetByUserDouble) {
    
    for (; doubleMapIt != m_doubleMap.end(); ++doubleMapIt) {
      
      if (doubleMapIt->second.overwriteVal > doubleMapIt->second.low && 
      doubleMapIt->second.overwriteVal < doubleMapIt->second.high) {
    
    *(doubleMapIt->second.par) = doubleMapIt->second.overwriteVal;
    
    log << MSG::INFO << " Overwriting: "
        << doubleMapIt->second.parDescription << " to : " 
        << doubleMapIt->second.overwriteVal   << endmsg;
      } 
    }
  }
  
  if (intMapIt->second.overwriteVal != m_parValueNotSetByUserInt) {
    
    for (; intMapIt != m_intMap.end(); ++intMapIt) {
      
      if (intMapIt->second.overwriteVal > intMapIt->second.low &&
      intMapIt->second.overwriteVal < intMapIt->second.high) {
    
    *(intMapIt->second.par) = intMapIt->second.overwriteVal;
    
    log << MSG::INFO << " Overwriting: " 
        << intMapIt->second.parDescription << " to : " 
        << intMapIt->second.overwriteVal   << endmsg;
      } 
    }
  }
}
 
void LUCID_DigiSettings::Print() const {
  
  MsgStream log(m_msgSvc, "LUCID_DigiSettings::Print");
  
  std::map <std::string, parDouble>::const_iterator doubleMapIt = m_doubleMap.begin();
  std::map <std::string, parInt>::const_iterator    intMapIt    = m_intMap.begin();
  
  log << MSG::INFO << "===================== Doubles ======================" << endmsg;
 
  for (; doubleMapIt != m_doubleMap.end(); ++doubleMapIt) 
    log << MSG::INFO << doubleMapIt->second.parDescription << ": " << *(doubleMapIt->second.par) <<endmsg; 
  
  log << MSG::INFO << "===================== Integers ======================" << endmsg;
    
  for (; intMapIt != m_intMap.end(); ++intMapIt)
    log << MSG::INFO  << intMapIt->second.parDescription << ": " << *(intMapIt->second.par) << endmsg;
 
  log << MSG::INFO << "=====================================================" << endmsg;
}
 
void LUCID_DigiSettings::SetDigiParDouble(const std::string& parname, double parval) {
  
  if (!m_doubleMap.count(parname))
    if (!m_intMap.count(parname)) {
      MsgStream log(m_msgSvc, "LUCID_DigiSettings::SetDigiParDouble");
      log << MSG::FATAL << " LUCID_DigiSettings: parameter does not exist: "  << parname << endmsg; 
      std::abort();
    }
  
  m_doubleMap[parname].par = &parval;
}
 
void LUCID_DigiSettings::SetDigiParInt(const std::string& parname, int parval) {
  
  if (!m_doubleMap.count(parname))
    if (!m_intMap.count(parname)) {
      MsgStream log(m_msgSvc, "LUCID_DigiSettings::SetDigiParInt");
      log << MSG::FATAL << " LUCID_DigiSettings: parameter does not exist: "  << parname << endmsg;
      std::abort();
    }
  
  m_intMap[parname].par = &parval;
}
 
double LUCID_DigiSettings::GetDigiParDouble(const std::string& parname) {
  
  if (!m_doubleMap.count(parname))
    if (!m_intMap.count(parname)) {
      MsgStream log(m_msgSvc, "LUCID_DigiSettings::GetDigiParDouble");
      log << MSG::FATAL << " LUCID_DigiSettings: parameter does not exist: "  << parname << endmsg;
      std::abort();
    }
  
  return *m_doubleMap[parname].par;
}
 
int LUCID_DigiSettings::GetDigiParInt(const std::string& parname) {
  
  if (!m_doubleMap.count(parname))
    if (!m_intMap.count(parname)) {
      MsgStream log(m_msgSvc, "LUCID_DigiSettings::GetDigiParInt");
      log << MSG::FATAL << " LUCID_DigiSettings: parameter does not exist: "  << parname << endmsg;
      std::abort();
    }
  
  return *m_intMap[parname].par;
}
 
void LUCID_DigiSettings::DefNewParameterDouble(std::string parDescription,
                           const std::string& parname,
                           double* par,
                           double  low,
                           double  high) { 
  
  if (low > high)                 {
    MsgStream log(m_msgSvc, "LUCID_DigiSettings::DefNewParameterDouble");
    log << MSG::FATAL << " Swap the boundaries on the range "       << parname << endmsg;
    std::abort();
  }
  if (m_doubleMap.count(parname)) {
    MsgStream log(m_msgSvc, "LUCID_DigiSettings::DefNewParameterDouble");
    log << MSG::FATAL << " More than one parameter with that name " << parname << endmsg;
    std::abort();
  }
  if (m_intMap.count   (parname)) {
    MsgStream log(m_msgSvc, "LUCID_DigiSettings::DefNewParameterDouble");
    log << MSG::FATAL << " More than one parameter with that name " << parname << endmsg;
    std::abort();
  }
  
  parDouble parD;
  
  parD.parDescription = std::move(parDescription);
  parD.par            = par;
  parD.low            = low;
  parD.high           = high;
  parD.overwriteVal   = m_parValueNotSetByUserDouble;
  
  m_doubleMap[parname] = parD;
}
 
void LUCID_DigiSettings::DefNewParameterInt(std::string parDescription,
                        const std::string& parname,
                        int* par,
                        int  low,
                        int  high) {
  
  if (low > high)                 {
    MsgStream log(m_msgSvc, "LUCID_DigiSettings::DefNewParameterInt");
    log << MSG::FATAL << " Swap the boundaries on the range "       << parname << endmsg;
    std::abort();
  }
  if (m_doubleMap.count(parname)) {
    MsgStream log(m_msgSvc, "LUCID_DigiSettings::DefNewParameterInt");
    log << MSG::FATAL << " More than one parameter with that name " << parname << endmsg;
    std::abort();
  }
  if (m_intMap.count   (parname)) {
    MsgStream log(m_msgSvc, "LUCID_DigiSettings::DefNewParameterInt");
    log << MSG::FATAL << " More than one parameter with that name " << parname << endmsg;
    std::abort();
  }
  
  parInt parI;
  
  parI.parDescription = std::move(parDescription);
  parI.par            = par;
  parI.low            = low;
  parI.high           = high;
  parI.overwriteVal   = m_parValueNotSetByUserInt;
  
  m_intMap[parname] = parI;
}