---

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 () ;
    
    void Fill (std::map<int,double>::const_iterator,
               std::map<int,double>::const_iterator,
               double pTk,
               double pSubtract,
               double sigma = 1.) ;

    void reset () ;
    void solve (int usingBlockSolver, double min, double max) ;
  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

/* TODOS and dubios
   ----------------

- fare l'inversione finale
- uso le CLHEP
- ritorno CLHEP o un vector<double>? forse la seconda e' meglio
- oppure una mappa con uindex, coeff? cosi' poi si possono incollare insieme
  tutte le mappe facilmente?
- CLHEP hanno qualche metodo utile per usare STL che non ho visto?
- rem se ho gia' invertito con una variabile: ha senso?

*/

---