IMACalibBlock.h

Go to the documentation of this file.

#ifndef __CINT__
#ifndef EcalCalibBlock_H
#define EcalCalibBlock_H

#include <map>
#include <string>
#include <vector>

#include "CLHEP/Matrix/GenMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Random/RandGaussQ.h"

#include "Calibration/Tools/interface/InvMatrixCommonDefs.h"
#include "Calibration/EcalCalibAlgos/interface/VEcalCalibBlock.h"

class IMACalibBlock : public VEcalCalibBlock {
public:
  IMACalibBlock(const int);
  ~IMACalibBlock() override;

  void Fill(std::map<int, double>::const_iterator,
            std::map<int, double>::const_iterator,
            double pTk,
            double pSubtract,
            double sigma = 1.) override;

  void reset() override;
  int solve(int usingBlockSolver, double min, double max) override;

private:
  int evalX2Size();
  void complete();
  void riempiMtr(const std::vector<double>& piena, CLHEP::HepMatrix& vuota);
  void riempiVtr(const std::vector<double>& pieno, CLHEP::HepVector& vuoto);
  void fillMap(const CLHEP::HepVector& result);

private:
  std::vector<double> m_kaliVector;
  std::vector<double> m_kaliMatrix;
};

#endif
#endif