RecoLocalTracker/SiStripZeroSuppression/src/SiStripRawProcessingFactory.cc - master - githubmirror/cmssw

RecoLocalTracker/SiStripZeroSuppression/src/SiStripRawProcessingFactory.cc
104 строки5 KB
Dan Riley
Fix logic error introduced in SiStripZeroSuppression/src/SiStripRawProcessingFactory.cc, convert some assertions to exceptions
20 мар 2023, 16:29
20 мар 2023, 16:2982d27f0
#include <memory>

#include "RecoLocalTracker/SiStripZeroSuppression/interface/SiStripRawProcessingFactory.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/SiStripRawProcessingAlgorithms.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/SiStripPedestalsSubtractor.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/SiStripFedZeroSuppression.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/MedianCMNSubtractor.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/PercentileCMNSubtractor.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/IteratedMedianCMNSubtractor.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/FastLinearCMNSubtractor.h"
#include "RecoLocalTracker/SiStripZeroSuppression/interface/TT6CMNSubtractor.h"

std::unique_ptr<SiStripRawProcessingAlgorithms> SiStripRawProcessingFactory::create(const edm::ParameterSet& conf,
                                                                                    edm::ConsumesCollector iC) {
  return std::unique_ptr<SiStripRawProcessingAlgorithms>(
      new SiStripRawProcessingAlgorithms(iC,
                                         create_SubtractorPed(conf, iC),
                                         create_SubtractorCMN(conf, iC),
                                         create_Suppressor(conf, iC),
                                         create_Restorer(conf, iC),
                                         conf.getParameter<bool>("doAPVRestore"),
                                         conf.getParameter<bool>("useCMMeanMap")));
}

std::unique_ptr<SiStripPedestalsSubtractor> SiStripRawProcessingFactory::create_SubtractorPed(
    const edm::ParameterSet& conf, edm::ConsumesCollector iC) {
  return std::unique_ptr<SiStripPedestalsSubtractor>(
      new SiStripPedestalsSubtractor(conf.getParameter<bool>("PedestalSubtractionFedMode"), iC));
}

std::unique_ptr<SiStripCommonModeNoiseSubtractor> SiStripRawProcessingFactory::create_SubtractorCMN(
    const edm::ParameterSet& conf, edm::ConsumesCollector iC) {
  const std::string mode = conf.getParameter<std::string>("CommonModeNoiseSubtractionMode");

  if (mode == "Median")
    return std::unique_ptr<SiStripCommonModeNoiseSubtractor>(new MedianCMNSubtractor());

  if (mode == "Percentile") {
    return std::unique_ptr<SiStripCommonModeNoiseSubtractor>(
        new PercentileCMNSubtractor(conf.getParameter<double>("Percentile")));
  }

  if (mode == "IteratedMedian") {
    return std::unique_ptr<SiStripCommonModeNoiseSubtractor>(new IteratedMedianCMNSubtractor(
        conf.getParameter<double>("CutToAvoidSignal"), conf.getParameter<int>("Iterations"), iC));
  }

  if (mode == "FastLinear")
    return std::unique_ptr<SiStripCommonModeNoiseSubtractor>(new FastLinearCMNSubtractor());

  if (mode == "TT6") {
    return std::unique_ptr<SiStripCommonModeNoiseSubtractor>(
        new TT6CMNSubtractor(conf.getParameter<double>("CutToAvoidSignal"), iC));
  }

  edm::LogError("SiStripRawProcessingFactory::create_SubtractorCMN")
      << "Unregistered Algorithm: " << mode << ". Use one of {Median, Percentile, IteratedMedian, FastLinear, TT6}";
  return std::unique_ptr<SiStripCommonModeNoiseSubtractor>(new MedianCMNSubtractor());
}

std::unique_ptr<SiStripFedZeroSuppression> SiStripRawProcessingFactory::create_Suppressor(const edm::ParameterSet& conf,
                                                                                          edm::ConsumesCollector iC) {
  const uint32_t mode = conf.getParameter<uint32_t>("SiStripFedZeroSuppressionMode");
  const bool trunc = conf.getParameter<bool>("TruncateInSuppressor");
  const bool trunc10bits = conf.getParameter<bool>("Use10bitsTruncation");
  switch (mode) {
    case 1:
    case 2:
    case 3:
    case 4:
      return std::make_unique<SiStripFedZeroSuppression>(mode, &iC, trunc, trunc10bits);
    default:
      edm::LogError("SiStripRawProcessingFactory::createSuppressor")
          << "Unregistered mode: " << mode << ". Use one of {1,2,3,4}.";
      return std::make_unique<SiStripFedZeroSuppression>(4, &iC, true, trunc10bits);
  }
}

std::unique_ptr<SiStripAPVRestorer> SiStripRawProcessingFactory::create_Restorer(const edm::ParameterSet& conf,
                                                                                 edm::ConsumesCollector iC) {
  if (!conf.exists("APVRestoreMode")) {
    return std::unique_ptr<SiStripAPVRestorer>(nullptr);
  } else {
    return std::unique_ptr<SiStripAPVRestorer>(new SiStripAPVRestorer(conf, iC));
  }
}

void SiStripRawProcessingFactory::fillDescriptions(edm::ParameterSetDescription& algorithms) {
  algorithms.add<std::string>("CommonModeNoiseSubtractionMode", "Median");
  algorithms.add("useCMMeanMap", false);
  algorithms.add("TruncateInSuppressor", true);
  algorithms.add("doAPVRestore", false);
  algorithms.add("SiStripFedZeroSuppressionMode", 4U);
  algorithms.add("PedestalSubtractionFedMode", true);
  algorithms.add("Use10bitsTruncation", false);
  algorithms.add("Percentile", 25.0);
  algorithms.add("CutToAvoidSignal", 2.0);
  algorithms.add("Iterations", 3);
  SiStripAPVRestorer::fillDescriptions(algorithms);
}