CondFormats/PhysicsToolsObjects/src/MVAComputer.cc - master - githubmirror/cmssw

CondFormats/PhysicsToolsObjects/src/MVAComputer.cc
192 строки7 KB
Matti Kortelainen
Make PhysicsTools::Calibration::MVAComputer::getProcessors() return const pointers
12 ноя 2021, 04:52
12 ноя 2021, 04:52bb0e6ee
// -*- C++ -*-
//
// Package:     PhysicsToolsObjects
// Class  :     MVAComputer
//

// Implementation:
//     getProcessor() and addProcessor() methods to add processors to
//     the discriminator computer calibration object. POOL doesn't support
//     polymorph pointers, so this is implemented using multiple containers
//     for each possible sub-class and an index array from which the
//     array of pointers can be reconstructed.
//
// Author:      Christophe Saout
// Created:     Sat Apr 24 15:18 CEST 2007
//
#include <algorithm>
#include <memory>

#include <typeinfo>
#include <iostream>
#include <cstring>
#include <cstddef>

#include <atomic>

#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Utilities/interface/TypeID.h"

#include "CondFormats/PhysicsToolsObjects/interface/MVAComputer.h"

namespace PhysicsTools {
  namespace Calibration {

    std::string VarProcessor::getInstanceName() const {
      static const char prefix[] = "PhysicsTools::Calibration::";
      edm::TypeID typeID(typeid(*this));
      const std::string &type(typeID.className());
      if (type.size() <= sizeof prefix - 1 || type.substr(0, sizeof prefix - 1) != prefix)
        throw cms::Exception("MVAComputerCalibration")
            << "getInstanceName failed for " << typeid(*this).name() << "." << std::endl;

      return type.substr(sizeof prefix - 1);
    }

    std::unique_ptr<VarProcessor> VarProcessor::clone() const { return (std::make_unique<VarProcessor>(*this)); }

    std::unique_ptr<VarProcessor> ProcOptional::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcOptional(*this)));
    }

    std::unique_ptr<VarProcessor> ProcCount::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcCount(*this)));
    }

    std::unique_ptr<VarProcessor> ProcClassed::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcClassed(*this)));
    }

    std::unique_ptr<VarProcessor> ProcSplitter::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcSplitter(*this)));
    }

    std::unique_ptr<VarProcessor> ProcForeach::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcForeach(*this)));
    }

    std::unique_ptr<VarProcessor> ProcSort::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcSort(*this)));
    }

    std::unique_ptr<VarProcessor> ProcCategory::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcCategory(*this)));
    }

    std::unique_ptr<VarProcessor> ProcNormalize::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcNormalize(*this)));
    }

    std::unique_ptr<VarProcessor> ProcLikelihood::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcLikelihood(*this)));
    }

    std::unique_ptr<VarProcessor> ProcLinear::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcLinear(*this)));
    }

    std::unique_ptr<VarProcessor> ProcMultiply::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcMultiply(*this)));
    }

    std::unique_ptr<VarProcessor> ProcMatrix::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcMatrix(*this)));
    }

    std::unique_ptr<VarProcessor> ProcExternal::clone() const {
      return (std::unique_ptr<VarProcessor>(new ProcExternal(*this)));
    }

    std::unique_ptr<VarProcessor> ProcMLP::clone() const { return (std::unique_ptr<VarProcessor>(new ProcMLP(*this))); }

    std::string ProcExternal::getInstanceName() const { return method; }

    static MVAComputer::CacheId getNextMVAComputerCacheId() {
      static std::atomic<MVAComputer::CacheId> nextCacheId{0};

      return ++nextCacheId;
    }

    MVAComputer::MVAComputer() : cacheId(getNextMVAComputerCacheId()) {}

    MVAComputer::MVAComputer(const MVAComputer &orig)
        : inputSet(orig.inputSet), output(orig.output), cacheId(orig.cacheId) {
      for (std::vector<VarProcessor *>::const_iterator iter = orig.processors.begin(); iter != orig.processors.end();
           ++iter)
        addProcessor(*iter);
    }

    MVAComputer::~MVAComputer() {
      for (std::vector<VarProcessor *>::iterator iter = processors.begin(); iter != processors.end(); ++iter)
        delete *iter;
      processors.clear();
    }

    MVAComputer &MVAComputer::operator=(const MVAComputer &orig) {
      inputSet = orig.inputSet;
      output = orig.output;
      cacheId = orig.cacheId;

      for (std::vector<VarProcessor *>::iterator iter = processors.begin(); iter != processors.end(); ++iter)
        delete *iter;
      processors.clear();

      for (std::vector<VarProcessor *>::const_iterator iter = orig.processors.begin(); iter != orig.processors.end();
           ++iter)
        addProcessor(*iter);

      return *this;
    }

    std::vector<VarProcessor *> MVAComputer::getProcessors() { return processors; }

    std::vector<const VarProcessor *> MVAComputer::getProcessors() const {
      std::vector<const VarProcessor *> ret(processors.size());
      std::copy(processors.begin(), processors.end(), ret.begin());
      return ret;
    }

    void MVAComputer::addProcessor(const VarProcessor *proc) {
      cacheId = getNextMVAComputerCacheId();
      processors.push_back(proc->clone().release());
    }

    static MVAComputerContainer::CacheId getNextMVAComputerContainerCacheId() {
      static MVAComputerContainer::CacheId nextCacheId = 0;
      return ++nextCacheId;
    }

    MVAComputerContainer::MVAComputerContainer() : cacheId(getNextMVAComputerContainerCacheId()) {}

    MVAComputer &MVAComputerContainer::add(const std::string &label) {
      cacheId = getNextMVAComputerContainerCacheId();

      entries.push_back(std::make_pair(label, MVAComputer()));
      return entries.back().second;
    }

    const MVAComputer &MVAComputerContainer::find(const std::string &label) const {
      std::vector<Entry>::const_iterator pos =
          std::find_if(entries.begin(), entries.end(), [&label](const MVAComputerContainer::Entry &entry) {
            return entry.first == label;
          });

      if (pos == entries.end())
        throw cms::Exception("MVAComputerCalibration")
            << "Calibration record " << label << " not found in MVAComputerContainer." << std::endl;

      return pos->second;
    }

    bool MVAComputerContainer::contains(const std::string &label) const {
      std::vector<Entry>::const_iterator pos =
          std::find_if(entries.begin(), entries.end(), [&label](const MVAComputerContainer::Entry &entry) {
            return entry.first == label;
          });
      if (pos == entries.end())
        return false;
      return true;
    }

  }  // namespace Calibration
}  // namespace PhysicsTools