BackgroundFunction.h

Go to the documentation of this file.

#ifndef BackgroundFunction_h
#define BackgroundFunction_h

#include <fstream>
#include <sstream>
#include "MuonAnalysis/MomentumScaleCalibration/interface/BaseFunction.h"
#include "MuonAnalysis/MomentumScaleCalibration/interface/Functions.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"

class BackgroundFunction : public BaseFunction {
public:
  BackgroundFunction(TString identifier) {
    identifier.Prepend("MuonAnalysis/MomentumScaleCalibration/data/");
    identifier.Append(".txt");
    edm::FileInPath fileWithFullPath(identifier.Data());
    readParameters(fileWithFullPath.fullPath());

    std::vector<int>::const_iterator idIt = functionId_.begin();
    for (; idIt != functionId_.end(); ++idIt)
      std::cout << "idIt = " << *idIt << std::endl;
  }
  BackgroundFunction(const MuScleFitDBobject* dbObject) : BaseFunction(dbObject) {
    std::vector<int>::const_iterator id = functionId_.begin();
    for (; id != functionId_.end(); ++id) {
      // TODO: fix the values for the lower and upper limits
      backgroundFunctionVec_.push_back(backgroundFunctionService(*id, 0., 200.));
    }
    // Fill the arrays that will be used when calling the correction function.
    convertToArrays(backgroundFunction_, backgroundFunctionVec_);
  }

  ~BackgroundFunction() {
    if (parArray_ != nullptr) {
      for (unsigned int i = 0; i < functionId_.size(); ++i) {
        delete[] parArray_[i];
        delete backgroundFunction_[i];
      }
      delete[] parArray_;
      delete[] backgroundFunction_;
    }
  }
  backgroundFunctionBase* function(const unsigned int i) {
    if (backgroundFunctionVec_.size() > i)
      return backgroundFunction_[i];
    else
      return nullptr;
  }

protected:
  void readParameters(TString fileName);

  backgroundFunctionBase** backgroundFunction_;
  std::vector<backgroundFunctionBase*> backgroundFunctionVec_;
};

#endif  // BackgroundFunction_h