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