IMACalibBlock Class Reference

#include <IMACalibBlock.h>

Inheritance diagram for IMACalibBlock:

Public Member Functions
void	Fill (std::map< int, double >::const_iterator, std::map< int, double >::const_iterator, double pTk, double pSubtract, double sigma=1.) override
	insert an entry More...
	IMACalibBlock (const int)
	ctor More...
void	reset () override
	reset the chi2 matrices More...
int	solve (int usingBlockSolver, double min, double max) override
	solve the chi2 linear system More...
	~IMACalibBlock () override
	dtor More...
Public Member Functions inherited from VEcalCalibBlock
float	at (const unsigned int index)
	retrieve the coefficients More...
	VEcalCalibBlock (int numberOfElements)
	ctor More...
virtual	~VEcalCalibBlock ()
	dtor More...

Private Member Functions
void	complete ()
	complete the triangolar chi2 matrix to a sym one More...
int	evalX2Size ()
	give the size of a chi2 matrix More...
void	fillMap (const CLHEP::HepVector &result)
	fill the coefficients map from the CLHEP vector solution More...
void	riempiMtr (const std::vector< double > &piena, CLHEP::HepMatrix &vuota)
	copy a vector into a CLHEP object More...
void	riempiVtr (const std::vector< double > &pieno, CLHEP::HepVector &vuoto)
	copy a vector into a CLHEP object More...

Private Attributes
std::vector< double >	m_kaliMatrix
	matrix for the chi2 inversion More...
std::vector< double >	m_kaliVector
	vector for the chi2 inversion More...

Additional Inherited Members
Protected Attributes inherited from VEcalCalibBlock
std::map< unsigned int, float >	m_coefficients
	map of coefficients More...
unsigned int	m_numberOfElements
	The only parameter! More...

Detailed Description

Definition at line 22 of file IMACalibBlock.h.

Constructor & Destructor Documentation

IMACalibBlock()

IMACalibBlock::IMACalibBlock

(

const int

numberOfElements

)

ctor

Definition at line 14 of file IMACalibBlock.cc.

    : VEcalCalibBlock(numberOfElements), m_kaliVector(m_numberOfElements), m_kaliMatrix(evalX2Size()) {
  reset();
}

References reset().

~IMACalibBlock()

IMACalibBlock::~IMACalibBlock ( )

override

dtor

Definition at line 21 of file IMACalibBlock.cc.

{}

Member Function Documentation

complete()

void IMACalibBlock::complete ( )

private

complete the triangolar chi2 matrix to a sym one

Definition at line 56 of file IMACalibBlock.cc.

                             {
  int bef;
  int aft;
  for (unsigned int i = 0; i < m_numberOfElements; ++i) {
    for (unsigned int j = i + 1; j < m_numberOfElements; ++j) {
      bef = (i * m_numberOfElements + j);
      aft = (j * m_numberOfElements + i);
      m_kaliMatrix.at(aft) = m_kaliMatrix.at(bef);
    }  //LP second loop over xtals
  }    //LP first loop over xtals
 
  return;
}

References mps_fire::i, dqmiolumiharvest::j, m_kaliMatrix, and VEcalCalibBlock::m_numberOfElements.

Referenced by solve().

evalX2Size()

int IMACalibBlock::evalX2Size ( )

inlineprivate

give the size of a chi2 matrix

Definition at line 108 of file IMACalibBlock.cc.

{ return m_numberOfElements * m_numberOfElements; }

References VEcalCalibBlock::m_numberOfElements.

Fill()

void IMACalibBlock::Fill ( std::map< int, double >::const_iterator  MapBegin, std::map< int, double >::const_iterator  MapEnd, double  pTk, double  pSubtract, double  sigma = 1.  )

overridevirtual

insert an entry

Implements VEcalCalibBlock.

Definition at line 25 of file IMACalibBlock.cc.

                                       {
  //  std::cerr<<"\n\nfaccio il fill!\n";
  double inverror = 1. / sigma;
  //LP fist loop over energies
  for (std::map<int, double>::const_iterator itMap1 = MapBegin; itMap1 != MapEnd; ++itMap1) {
    //      std::cerr<<"numero "<<itMap1->first<<" vale "<<itMap1->second<<"\t";
    for (std::map<int, double>::const_iterator itMap2 = itMap1; itMap2 != MapEnd; ++itMap2) {
      //LP calculate the chi square value
      double dummy = itMap1->second * itMap2->second;
      dummy *= inverror;
      //LP fill the calib matrix
      m_kaliMatrix.at(itMap1->first * m_numberOfElements + itMap2->first) += dummy;
    }  //LP second loop over xtals
 
    //LP calculate the vector value
    double dummy = pTk;
    dummy -= pSubtract;  //LP borders
    dummy *= itMap1->second;
    dummy *= inverror;
    //LP fill the calib vector
    m_kaliVector.at(itMap1->first) += dummy;
  }  //LP first loop over energies
  return;
}

References m_kaliMatrix, m_kaliVector, and VEcalCalibBlock::m_numberOfElements.

fillMap()

void IMACalibBlock::fillMap ( const CLHEP::HepVector &  result )

private

fill the coefficients map from the CLHEP vector solution

Definition at line 157 of file IMACalibBlock.cc.

                                                        {
  for (unsigned int i = 0; i < m_numberOfElements; ++i) {
    m_coefficients[i] = result[i];
  }
 
  return;
}

References mps_fire::i, VEcalCalibBlock::m_coefficients, VEcalCalibBlock::m_numberOfElements, and mps_fire::result.

Referenced by solve().

reset()

void IMACalibBlock::reset ( void  )

overridevirtual

reset the chi2 matrices

Implements VEcalCalibBlock.

Definition at line 144 of file IMACalibBlock.cc.

                          {
  for (std::vector<double>::iterator vecIt = m_kaliVector.begin(); vecIt != m_kaliVector.end(); ++vecIt) {
    *vecIt = 0.;
  }
 
  for (std::vector<double>::iterator vecIt = m_kaliMatrix.begin(); vecIt != m_kaliMatrix.end(); ++vecIt) {
    *vecIt = 0.;
  }
}

References m_kaliMatrix, and m_kaliVector.

Referenced by MatrixReader.MatrixReader::__init__(), IMACalibBlock(), data_sources.json_list::next(), and MatrixReader.MatrixReader::showRaw().

riempiMtr()

void IMACalibBlock::riempiMtr ( const std::vector< double > &  piena, CLHEP::HepMatrix &  vuota  )

private

copy a vector into a CLHEP object

Definition at line 112 of file IMACalibBlock.cc.

                                                                                   {
  unsigned int max = m_numberOfElements;
 
  assert(piena.size() == max * max);
  assert(vuota.num_row() == int(max));
  assert(vuota.num_col() == int(max));
  for (unsigned int i = 0; i < max; ++i)
    for (unsigned int j = 0; j < max; ++j)
      if (edm::isNotFinite(piena[i * max + j]))
        vuota[i][j] = 0.;
      else
        vuota[i][j] = piena[i * max + j];
 
  return;
}

References cms::cuda::assert(), mps_fire::i, createfilelist::int, edm::isNotFinite(), dqmiolumiharvest::j, VEcalCalibBlock::m_numberOfElements, and SiStripPI::max.

Referenced by solve().

riempiVtr()

void IMACalibBlock::riempiVtr ( const std::vector< double > &  pieno, CLHEP::HepVector &  vuoto  )

private

copy a vector into a CLHEP object

Definition at line 130 of file IMACalibBlock.cc.

                                                                                   {
  int max = m_numberOfElements;
  assert(vuoto.num_row() == max);
  for (int i = 0; i < max; ++i)
    if (edm::isNotFinite(pieno[i]))
      vuoto[i] = 0.;
    else
      vuoto[i] = pieno[i];
 
  return;
}

References cms::cuda::assert(), mps_fire::i, edm::isNotFinite(), VEcalCalibBlock::m_numberOfElements, and SiStripPI::max.

Referenced by solve().

solve()

int IMACalibBlock::solve ( int  usingBlockSolver, double  min, double  max  )

overridevirtual

solve the chi2 linear system

Implements VEcalCalibBlock.

Definition at line 72 of file IMACalibBlock.cc.

                                                                     {
  complete();
  // TH1F vettore ("vettore","vettore",10,-0.1,9.9);
  // TH1F matrice ("matrice","matrice",100,-0.1,99.9);
  CLHEP::HepMatrix kaliMatrix(m_numberOfElements, m_numberOfElements);
  // for (std::vector<double>::iterator it = m_kaliVector.begin();
  //         it!= m_kaliVector.end();++it)
  //     vettore.Fill(it-m_kaliVector.begin(),*it);
  riempiMtr(m_kaliMatrix, kaliMatrix);
  // for (std::vector<double>::iterator it = m_kaliMatrix.begin();
  //         it!= m_kaliMatrix.end();++it)
  //     matrice.Fill(it-m_kaliMatrix.begin(),*it);
  //  TFile f ("fileInteressante.root","RECREATE");
  //  vettore.Write();
  // matrice.Write();
  CLHEP::HepVector kaliVector(m_numberOfElements);
  riempiVtr(m_kaliVector, kaliVector);
  //PG linear system solution
  CLHEP::HepVector result = CLHEP::solve(kaliMatrix, kaliVector);
  for (int i = 0; i < kaliVector.num_row(); ++i)
    if (result.normsq() < min * kaliVector.num_row() || result.normsq() > max * kaliVector.num_row()) {
      if (usingBlockSolver) {
        edm::LogWarning("IML") << "using  blocSlover " << std::endl;
        BlockSolver()(kaliMatrix, kaliVector, result);
      } else {
        edm::LogWarning("IML") << "coeff out of range " << std::endl;
        for (int i = 0; i < kaliVector.num_row(); ++i)
          result[i] = 1.;
      }
    }
  fillMap(result);
  return 0;
}

References complete(), fillMap(), mps_fire::i, m_kaliMatrix, m_kaliVector, VEcalCalibBlock::m_numberOfElements, SiStripPI::max, min(), mps_fire::result, riempiMtr(), and riempiVtr().

Member Data Documentation

m_kaliMatrix

std::vector<double> IMACalibBlock::m_kaliMatrix

private

matrix for the chi2 inversion

Definition at line 57 of file IMACalibBlock.h.

Referenced by complete(), Fill(), reset(), and solve().

m_kaliVector

std::vector<double> IMACalibBlock::m_kaliVector

private

vector for the chi2 inversion

Definition at line 55 of file IMACalibBlock.h.

Referenced by Fill(), reset(), and solve().