FitAmpPairCovariance Class Reference

#include <FitAmpPairCovariance.h>

Public Member Functions
	FitAmpPairCovariance (const FitAmpPairList *myList)
	FitAmpPairCovariance (const FitAmpPairCovariance &other)
	FitAmpPairCovariance (const FitAmpPairCovariance &other, const FitAmpPairList *newList)
virtual	~FitAmpPairCovariance ()
bool	add (const FitAmpPairCovariance &other)
bool	reset ()
bool	addLastEventFromList ()
bool	calculateAll ()
double	getFractionVariance (int i)
double	getFractionError (int i)
double	getFractionSumVariance ()
double	getFractionSumError ()
double	getIntegralVariance ()
double	getIntegralError ()
unsigned int	size () const
bool	isValid () const
bool	save (const std::string &asSubdirOf=".") const
bool	retrieve (const std::string &fromDirectory=".")
int	Nevents () const
FitAmpPairCovariance &	operator+= (const FitAmpPairCovariance &other)
FitAmpPairCovariance	operator+ (const FitAmpPairCovariance &other) const

Private Member Functions
std::string	realName (int listPosition) const
std::string	imagName (int listPosition) const
int	indexReal (int listPosition) const
int	indexImag (int listPosition) const
int	listPositionFromIndex (int idx) const
bool	make_dN_by_d2N ()
bool	make_dFraction_by_dN ()
bool	make_dFractionSum_by_dFraction ()
bool	make_dIntegral_by_dN ()
bool	resize ()
unsigned int	maxSize () const
bool	make2NCovariance ()
bool	makeNCovariance ()
bool	makeFractionCovariance ()
bool	makeFractionSumCovariance ()
bool	makeIntegralCovariance ()
bool	clearAll ()
bool	makeDirectory (const std::string &asSubdirOf=".") const
std::string	dirName (const std::string &asSubdirOf=".") const
std::string	matrix_x_fname (const std::string &asSubdirOf=".") const
std::string	matrix_xy_fname (const std::string &asSubdirOf=".") const
bool	save_sum_x (const std::string &inDirectory) const
bool	save_sum_xy (const std::string &inDirectory) const
bool	retrieve_sum_x (const std::string &inDirectory)
bool	retrieve_sum_xy (const std::string &inDirectory)
bool	isCompatibleWith (const FitAmpPairCovariance &other) const

Private Attributes
const FitAmpPairList *	_myList
unsigned int	_Nevents
unsigned int	_size
bool	_needToRecalculate
TMatrixT< float >	_sum_x
TMatrixTSym< float >	_sum_xy
TMatrixTSym< float >	_cov2N
TMatrixTSym< float >	_corr2N
TMatrixTSym< float >	_cov
TMatrixTSym< float >	_corr
TMatrixTSym< float >	_fractionCov
TMatrixTSym< float >	_fractionSumCov
TMatrixTSym< float >	_integralCov
TMatrixT< float >	_dN_by_d2N
TMatrixT< float >	_dFraction_by_dN
TMatrixT< float >	_dFractionSum_by_dFraction
TMatrixT< float >	_dIntegral_by_dN

Static Private Attributes
static unsigned int	_maxSize = 1500

Detailed Description

Definition at line 11 of file FitAmpPairCovariance.h.

Constructor & Destructor Documentation

FitAmpPairCovariance() [1/3]

FitAmpPairCovariance::FitAmpPairCovariance

(

const FitAmpPairList *

myList

)

Definition at line 18 of file FitAmpPairCovariance.cpp.

  : _myList(myList)
  , _Nevents(0)
  , _size()
  , _needToRecalculate(true)
  , _sum_x(1, 1)
  , _sum_xy(1)
  , _cov2N(2)
  , _corr2N(2)
  , _cov(1)
  , _corr(1)
  , _fractionCov(1)
  , _fractionSumCov(1)
  , _integralCov(1)
  , _dN_by_d2N(1, 2)
  , _dFraction_by_dN(1, 1)
  , _dFractionSum_by_dFraction(1, 1)
  , _dIntegral_by_dN(1, 1)
{  
  if(0 != _myList && _myList->size() >=1) resize();
}

FitAmpPairCovariance() [2/3]

FitAmpPairCovariance::FitAmpPairCovariance

(

const FitAmpPairCovariance &

other

)

Definition at line 39 of file FitAmpPairCovariance.cpp.

  : _myList(other._myList)
  , _Nevents(other._Nevents)
  , _size(other._size)
  , _needToRecalculate(true)
  , _sum_x(other._sum_x)
  , _sum_xy(other._sum_xy)
  , _cov2N(other._cov2N)
  , _corr2N(other._corr2N)
  , _cov(other._cov)
  , _corr(other._corr)
  , _fractionCov(other._fractionCov)
  , _fractionSumCov(other._fractionSumCov)
  , _integralCov(other._integralCov)
  , _dN_by_d2N(other._dN_by_d2N)
  , _dFraction_by_dN(other._dFraction_by_dN)
  , _dFractionSum_by_dFraction(other._dFractionSum_by_dFraction)
  , _dIntegral_by_dN(other._dIntegral_by_dN)
{
}

FitAmpPairCovariance() [3/3]

FitAmpPairCovariance::FitAmpPairCovariance	(	const FitAmpPairCovariance &	other,
		const FitAmpPairList *	newList
	)

Definition at line 59 of file FitAmpPairCovariance.cpp.

  : _myList(newList)
  , _Nevents(other._Nevents)
  , _size(other._size)
  , _needToRecalculate(true)
  , _sum_x(other._sum_x)
  , _sum_xy(other._sum_xy)
  , _cov2N(other._cov2N)
  , _corr2N(other._corr2N)
  , _cov(other._cov)
  , _corr(other._corr)
  , _fractionCov(other._fractionCov)
  , _fractionSumCov(other._fractionSumCov)
  , _integralCov(other._integralCov)
  , _dN_by_d2N(other._dN_by_d2N)
  , _dFraction_by_dN(other._dFraction_by_dN)
  , _dFractionSum_by_dFraction(other._dFractionSum_by_dFraction)
  , _dIntegral_by_dN(other._dIntegral_by_dN)
{  
}

~FitAmpPairCovariance()

FitAmpPairCovariance::~FitAmpPairCovariance ( )

virtual

Definition at line 81 of file FitAmpPairCovariance.cpp.

                                           {
  clearAll();
}

Member Function Documentation

add()

bool FitAmpPairCovariance::add

(

const FitAmpPairCovariance &

other

)

Definition at line 99 of file FitAmpPairCovariance.cpp.

                                                               {
  if(0 == _myList || 0 == other._myList) {
    cout << "WARNING in FitAmpPairCovariance::add: "
     << "trying to add incompatible FitAmpPairCovariances"
     << "\n\t _myListPtr = " << _myList << " other's: " << other._myList
     << endl;
    return false;
  }
  if(! isCompatibleWith(other)){
     cout << "WARNING in FitAmpPairCovariance::add: "
     << "trying to add incompatible FitAmpPairCovariances"
     << endl;   
    return false;
  }
  _needToRecalculate = true;
  _Nevents += other._Nevents;
  _sum_x   += other._sum_x;
  _sum_xy  += other._sum_xy;

  return true;
}

addLastEventFromList()

bool FitAmpPairCovariance::addLastEventFromList

(

)

Definition at line 177 of file FitAmpPairCovariance.cpp.

                                               {
  bool dbThis=false;
  if(size() != _myList->size()) {
    if(dbThis) cout << "FitAmpPairCovariance::addLastEventFromList() size() " 
            << size() << " == " <<  _myList->size() << endl;
    if(! resize()) return false;
  }
  _needToRecalculate = true;
  _Nevents++;
  if(dbThis)cout << "FitAmpPairCovariance::addLastEventFromList()"
         << " adding event to cov matrix." << endl;

  if(! isValid()) return false;
  for(unsigned int i=0; i < _myList->size(); i++){
    if(dbThis) cout << " i= " << i << endl;
    double re_i =  _myList->at(i).lastEntry().real();
    double im_i =  _myList->at(i).lastEntry().imag();
    _sum_x(indexReal(i), 0) += re_i;
    _sum_x(indexImag(i), 0) += im_i;

    for(unsigned int j=0; j < _myList->size(); j++){
      double re_j =  _myList->at(j).lastEntry().real();
      double im_j =  _myList->at(j).lastEntry().imag();
 
      _sum_xy(indexReal(i), indexReal(j)) += re_i * re_j;
      _sum_xy(indexReal(i), indexImag(j)) += re_i * im_j;
      _sum_xy(indexImag(i), indexReal(j)) += im_i * re_j;
      _sum_xy(indexImag(i), indexImag(j)) += im_i * im_j;
    }
  }
  return true;
}

calculateAll()

bool FitAmpPairCovariance::calculateAll

(

)

Definition at line 210 of file FitAmpPairCovariance.cpp.

                                       {
  bool dbThis=false;
  if(! isValid()) return false;
  bool sc=true;
  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making dN_by_d2N" << endl;
  sc &= make_dN_by_d2N();

  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making dFraction_by_dN" << endl;
  sc &= make_dFraction_by_dN();

  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making dFractionSum_by_dFraction" << endl;
  sc &= make_dFractionSum_by_dFraction();

  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making dIntegral_by_dN" << endl;
  sc &= make_dIntegral_by_dN();

  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making 2NCovariance" << endl;
  sc &= make2NCovariance();
  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making NCovariance" << endl;
  sc &= makeNCovariance();
  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making FractionCovariance" << endl;
  sc &= makeFractionCovariance();
  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making FractionSumCovariance" << endl;
  sc &= makeFractionSumCovariance();
  if(dbThis) cout << "FitAmpPairCovariance::calculateAll"
          << " making IntegralCovariance" << endl;
  sc &= makeIntegralCovariance();

  

  //_needToRecalculate=false;
  return sc;
}

clearAll()

bool FitAmpPairCovariance::clearAll ( )

private

Definition at line 573 of file FitAmpPairCovariance.cpp.

                                   {
  _Nevents=0;
  _needToRecalculate=true;
  resize();
  if(0 != this->size()){
    _sum_x  *= 0;
    _sum_xy *= 0;
  }
  return 0;
}

dirName()

std::string FitAmpPairCovariance::dirName ( const std::string &  asSubdirOf = "." ) const

private

Definition at line 625 of file FitAmpPairCovariance.cpp.

                                                                        {
  return asSubdirOf +"/" + "FitAmpPairCovariance";
}

getFractionError()

double FitAmpPairCovariance::getFractionError

(

int

i

)

Definition at line 501 of file FitAmpPairCovariance.cpp.

                                                  {
  double V = getFractionVariance(i);
  if(V < 0){
    cout << "ERROR in FitAmpPairCovariance::getFractionError"
     << " variance is " << V << " which is < 0."
     << " will return -9999."
     << endl;
    return -9999.0;
  }
  return sqrt(V);
}

getFractionSumError()

double FitAmpPairCovariance::getFractionSumError

(

)

Definition at line 518 of file FitAmpPairCovariance.cpp.

                                                {
  double V = getFractionSumVariance();
  if(V < 0){
    cout << "ERROR in FitAmpPairCovariance::getFractionSumError"
     << " variance is " << V << " which is < 0."
     << " will return 0."
     << endl;
    return 0;
  }
  return sqrt(V);
}

getFractionSumVariance()

double FitAmpPairCovariance::getFractionSumVariance

(

)

Definition at line 513 of file FitAmpPairCovariance.cpp.

                                                   {
  if(_needToRecalculate)  calculateAll();
  return _fractionSumCov(0,0);
}

getFractionVariance()

double FitAmpPairCovariance::getFractionVariance

(

int

i

)

Definition at line 495 of file FitAmpPairCovariance.cpp.

                                                     {
  if(_needToRecalculate)  calculateAll();
  if(i < 0 || i >= _fractionCov.GetNcols()) return -9999.0;
  if(i < 0 || i >= _fractionCov.GetNrows()) return -9999.0;
  return _fractionCov(i, i);
}

getIntegralError()

double FitAmpPairCovariance::getIntegralError

(

)

Definition at line 537 of file FitAmpPairCovariance.cpp.

                                             {
  double V = getIntegralVariance();
  if(V < 0){
    cout << "ERROR in FitAmpPairCovariance::getFractionError"
     << " variance is " << V << " which is < 0."
     << " will return -9999."
     << endl;
    return -9999.0;
  }
  return sqrt(V);
}

getIntegralVariance()

double FitAmpPairCovariance::getIntegralVariance

(

)

Definition at line 530 of file FitAmpPairCovariance.cpp.

                                                {
  bool dbThis=false;
  if(dbThis) cout << "need to recalculate ? " << _needToRecalculate << endl;
  if(_needToRecalculate)  calculateAll();
  if(_integralCov.GetNcols() < 1 || _integralCov.GetNrows() < 1) return -9999.0;
  return _integralCov(0,0);
}

imagName()

std::string FitAmpPairCovariance::imagName ( int  listPosition ) const

private

Definition at line 173 of file FitAmpPairCovariance.cpp.

                                                              {
  return _myList->at(listPosition).name() + "_imag";
}

indexImag()

int FitAmpPairCovariance::indexImag ( int  listPosition ) const

inlineprivate

Definition at line 36 of file FitAmpPairCovariance.h.

{return 2*listPosition + 1;}

indexReal()

int FitAmpPairCovariance::indexReal ( int  listPosition ) const

inlineprivate

Definition at line 35 of file FitAmpPairCovariance.h.

{return 2*listPosition;}

isCompatibleWith()

bool FitAmpPairCovariance::isCompatibleWith ( const FitAmpPairCovariance &  other ) const

private

Definition at line 84 of file FitAmpPairCovariance.cpp.

                                    {
  if(0 == _myList && 0 == other._myList) return true;
  if(0 == _myList || 0 == other._myList) return false;
  if(_myList->size() != other._myList->size()) return false;

  if(! _myList->isCompatibleWith(*(other._myList))) return false;

  // the following check should be redundant:
  for(unsigned int i=0; i < this->size(); i++){
    if(realName(i) != other.realName(i)) return false;
    if(imagName(i) != other.imagName(i)) return false;
  }
  return true;
}

isValid()

bool FitAmpPairCovariance::isValid

(

)

const

Definition at line 549 of file FitAmpPairCovariance.cpp.

                                       {
  bool dbThis=false;
  static int numWarnings=0;
  const static int maxWarnings=10;
  if(dbThis){
    cout << " FitAmpPairCovariance::isValid() with: "
     << size() << " == " << _sum_x.GetNcols() << " ? " << endl;
  }
  if(size() > maxSize()){
    if(numWarnings < 10){
      cout << "WARNING: FitAmpPairCovariance::isValid returns false"
       << " because size (" << size() << ") is too big"
       << endl;
      cout << " Will print this " << maxWarnings 
       << " times (" << numWarnings << ")." << endl;
    }
    numWarnings++;
    return false;
  }
  if(size() <= 0)return false;
  //if((unsigned int) _sum_x.GetNcols() != 2*size()) return false;
  return true;
}

listPositionFromIndex()

int FitAmpPairCovariance::listPositionFromIndex ( int  idx ) const

inlineprivate

Definition at line 37 of file FitAmpPairCovariance.h.

{return idx/2;}

make2NCovariance()

bool FitAmpPairCovariance::make2NCovariance ( )

private

Definition at line 349 of file FitAmpPairCovariance.cpp.

                                           {
  bool dbThis=false;
  if(dbThis) cout << "making 2N covariance" << endl;

  double dN = _myList->numEvents();
  TMatrixT<float> mean_x(_myList->size()*2, 1);
  TMatrixTSym<float> mean_xy(_myList->size()*2);

  mean_x  = _sum_x;
  mean_x *= 1./dN;

  mean_xy  = _sum_xy;
  mean_xy *= 1./dN;
                
  for(unsigned int i=0; i < _myList->size()*2; i++){
    for(unsigned int j=0; j < _myList->size()*2; j++){
      _cov2N(i, j) = (mean_xy(i,j) -  mean_x(i,0)*mean_x(j,0));
      // above: cov of the distribution
      // now: .. of the mean
      _cov2N(i, j) /= dN;
      // now taking into account other factors that I muliply the result
      // with.
      int il = listPositionFromIndex(i);
      int jl = listPositionFromIndex(j);
      double sfden = _myList->at(il).weightSum() * _myList->at(jl).weightSum();
      if(sfden > 0){
    _cov2N(i, j) *= dN*dN/sfden;
      }
    }
  }
  for(unsigned int i=0; i < _myList->size()*2; i++){
    for(unsigned int j=0; j < _myList->size()*2; j++){
      double den = _cov2N(i,i)*_cov2N(j,j);
      if(den > 0){
    _corr2N(i,j) = _cov2N(i,j)/sqrt(den);
      }
    }
  }

  if(dbThis){
    cout << "After " << _Nevents << "events" << endl;
    cout << "mean matrix" << endl;
    mean_x.Print();
    cout << "Covariance Matrix 2N" << endl;
    _cov2N.Print();
    cout << "Correlation Matrix 2N" << endl;
    _corr2N.Print();
  }
  return true;

}

make_dFraction_by_dN()

bool FitAmpPairCovariance::make_dFraction_by_dN ( )

private

Definition at line 278 of file FitAmpPairCovariance.cpp.

                                               {
  bool dbThis=false;
  if(dbThis) cout << "FitAmpPairCovariance::make_dFraction_by_dN() called" 
          << endl;
  /*
    We are calculating trafo matrix for the following trafo:

    f_a = a/(a + b + c...)


    diagonal terms:

      d(f_a)/d(a) = 1/(a + b + c +...) - a/(a + b + c ...)^2

    off-diagonal terms:

      d(f_a)/d(b) = -a/(a + b + c ...)^2
   */

  double integ = _myList->integral();
  double integ2 = integ*integ;
  for(unsigned int i=0; i < _myList->size(); i++){
    double num = _myList->at(i).integral();
    for(unsigned int j=0; j < _myList->size(); j++){
      _dFraction_by_dN(i, j) = -num/integ2;
    }
    _dFraction_by_dN(i, i) += 1.0/integ;
  }
  return true;
}

make_dFractionSum_by_dFraction()

bool FitAmpPairCovariance::make_dFractionSum_by_dFraction ( )

private

Definition at line 326 of file FitAmpPairCovariance.cpp.

                                                         {
  bool dbThis=false;
  if(dbThis) cout << "FitAmpPairCovariance::make_dFractionSum_by_dFraction()"
          << " called" << endl;
  /*
    We are calculating trafo matrix for the trafo to the sum 
    of all terms entering the fraction calculation, this 
    excludes the terms for interference effect.

    It's just ones (those included) and zeros (for those not included)
   */

  for(unsigned int i=0; i < _myList->size(); i++){
    if(_myList->at(i).isSingleAmp()){
      _dFractionSum_by_dFraction(0, i) = 1.0;
    }else{
      _dFractionSum_by_dFraction(0, i) = 0.0;
    }
  }
  return true;
}

make_dIntegral_by_dN()

bool FitAmpPairCovariance::make_dIntegral_by_dN ( )

private

Definition at line 309 of file FitAmpPairCovariance.cpp.

                                               {
  bool dbThis=false;
  if(dbThis) cout << "FitAmpPairCovariance::make_dIntegral_by_dN() called" 
          << endl;
  /*
    We are calculating trafo matrix for the following trafo:

    sum = a + b + c...

   */

  for(unsigned int i=0; i < _myList->size(); i++){
    _dIntegral_by_dN(0, i) = 1;
  }
  return true;
}

make_dN_by_d2N()

bool FitAmpPairCovariance::make_dN_by_d2N ( )

private

Definition at line 251 of file FitAmpPairCovariance.cpp.

                                         {
  bool dbThis=false;
  if(dbThis) cout << "FitAmpPairCovariance::make_dN_by_d2N() called" << endl;

  /*
    each FitAmpPair is evaluated as the real part of
    the product of two complex numbers.
    The product of fit parameters: oneOrTwo * (x_1 + iy_1)*(x_2 - iy_2)
    The integral: a + ib
    So we want to transfrom from (a, b) to c = xa - yb.

    Here we do the covariance due to the statistical uncertainty
    in the MC integration, not the uncertainy in the fit parameters x,y.

    So we calculate d(c) by d(a,b) = (x, -y) 
   */
  for(unsigned int i=0; i < _myList->size(); i++){
    complex<double> fpv = 
      _myList->at(i).fitParValue()
      * (double)(_myList->at(i).oneOrTwo());

    _dN_by_d2N(i, indexReal(i)) =   fpv.real();
    _dN_by_d2N(i, indexImag(i)) = - fpv.imag();
  }
  return true;
}

makeDirectory()

bool FitAmpPairCovariance::makeDirectory ( const std::string &  asSubdirOf = "." ) const

private

Definition at line 590 of file FitAmpPairCovariance.cpp.

                                                                        {
  /*
    A mode is created from or'd permission bit masks defined
     in <sys/stat.h>:
           #define S_IRWXU 0000700     RWX mask for owner 
           #define S_IRUSR 0000400     R for owner 
           #define S_IWUSR 0000200     W for owner 
           #define S_IXUSR 0000100     X for owner 

           #define S_IRWXG 0000070     RWX mask for group 
           #define S_IRGRP 0000040     R for group 
           #define S_IWGRP 0000020     W for group 
           #define S_IXGRP 0000010     X for group 

           #define S_IRWXO 0000007     RWX mask for other 
           #define S_IROTH 0000004     R for other 
           #define S_IWOTH 0000002     W for other 
           #define S_IXOTH 0000001     X for other 

           #define S_ISUID 0004000     set user id on execution 
           #define S_ISGID 0002000     set group id on execution 
           #define S_ISVTX 0001000     save swapped text even after use
   */

  mode_t mode = S_IRWXU | S_IRWXG | S_IRWXO 
              | S_ISUID | S_ISGID; 
  // see above for meaning. I want everybody to be allowed to read/write/exec.
  // Not sure about the last two bits.

  int zeroForSuccess = 0;
  zeroForSuccess |= mkdir( (asSubdirOf ).c_str(), mode );
  zeroForSuccess |= mkdir( (dirName(asSubdirOf) ).c_str(), mode );
  return (0 == zeroForSuccess);
}

makeFractionCovariance()

bool FitAmpPairCovariance::makeFractionCovariance ( )

private

Definition at line 441 of file FitAmpPairCovariance.cpp.

                                                 {
  bool dbThis=false;

  TMatrixT<float> dFraction_by_dN_Transpose(size(), size());
  dFraction_by_dN_Transpose.Transpose(_dFraction_by_dN);

  _fractionCov = makeTMatrixTSym(_dFraction_by_dN * _cov * 
                 dFraction_by_dN_Transpose);

  if(dbThis){
    cout << "Correlation Matrix Fractions" << endl;
    _fractionCov.Print();
  }

  return true;
}

makeFractionSumCovariance()

bool FitAmpPairCovariance::makeFractionSumCovariance ( )

private

Definition at line 457 of file FitAmpPairCovariance.cpp.

                                                    {
  bool dbThis=false;

  TMatrixT<float> dFractionSum_by_dFraction_Transpose(size(), 1);
  dFractionSum_by_dFraction_Transpose.Transpose(_dFractionSum_by_dFraction);

  _fractionSumCov = makeTMatrixTSym(_dFractionSum_by_dFraction *
                    _fractionCov * 
                    dFractionSum_by_dFraction_Transpose);

  if(dbThis){
    cout << "After " << _Nevents << "events" << endl;
    cout << "Covariance Matrix for FractionSum" << endl;
    _fractionSumCov.Print();
   }

  return true;
}

makeIntegralCovariance()

bool FitAmpPairCovariance::makeIntegralCovariance ( )

private

Definition at line 476 of file FitAmpPairCovariance.cpp.

                                                 {
  bool dbThis=false;

  TMatrixT<float> dIntegral_by_dN_Transpose(size(), 1);
  dIntegral_by_dN_Transpose.Transpose(_dIntegral_by_dN);

  _integralCov = makeTMatrixTSym(_dIntegral_by_dN * 
                 _cov 
                 * dIntegral_by_dN_Transpose);

  if(dbThis){
    cout << "Correlation Matrix Integral" << endl;
    _fractionCov.Print();
  }

  return true;
}

makeNCovariance()

bool FitAmpPairCovariance::makeNCovariance ( )

private

Definition at line 401 of file FitAmpPairCovariance.cpp.

                                          {
  bool dbThis=false;

  bool success = true;
  //  success &= make2NCovariance();
  //  success &= make_dN_by_d2N();

  TMatrixT<float> dN_by_d2N_Transpose(2*size(), size());
  dN_by_d2N_Transpose.Transpose(_dN_by_d2N);
  if(dbThis){
    cout << "FitAmpPairCovariance::makeNCovariance()"
     << " _dN_by_d2N(): " 
     << endl;
    _dN_by_d2N.Print();
    cout << "and its transpose:" << endl;
    dN_by_d2N_Transpose.Print();
  }

  _cov = makeTMatrixTSym(_dN_by_d2N * _cov2N * dN_by_d2N_Transpose);

  for(unsigned int i=0; i < _myList->size(); i++){
    if(_cov(i, i) <= 0) continue;
    for(unsigned int j=0; j < _myList->size(); j++){
      double den = _cov(i,i)*_cov(j,j);
      if(den > 0){
    _corr(i,j) = _cov(i,j)/sqrt(den);
      }
    }
  }

  if(dbThis){
    cout << "After " << _Nevents << "events" << endl;
    cout << "Covariance Matrix N" << endl;
    _cov.Print();
    cout << "Correlation Matrix N" << endl;
    _corr.Print();
  }
  return success;
}

matrix_x_fname()

std::string FitAmpPairCovariance::matrix_x_fname ( const std::string &  asSubdirOf = "." ) const

private

Definition at line 628 of file FitAmpPairCovariance.cpp.

                                                                               {
  return dirName(asSubdirOf) + "/" + "matrix_x.txt";
}

matrix_xy_fname()

std::string FitAmpPairCovariance::matrix_xy_fname ( const std::string &  asSubdirOf = "." ) const

private

Definition at line 631 of file FitAmpPairCovariance.cpp.

                                                                                {
  return dirName(asSubdirOf) + "/" + "matrix_xy.txt";
}

maxSize()

unsigned int FitAmpPairCovariance::maxSize ( ) const

inlineprivate

Definition at line 45 of file FitAmpPairCovariance.h.

{return _maxSize;}

Nevents()

int FitAmpPairCovariance::Nevents ( ) const

inline

Definition at line 94 of file FitAmpPairCovariance.h.

{return _Nevents;}

operator+()

FitAmpPairCovariance FitAmpPairCovariance::operator+

(

const FitAmpPairCovariance &

other

)

const

Definition at line 814 of file FitAmpPairCovariance.cpp.

                                                                                           {
  FitAmpPairCovariance returnVal(*this);
  return (returnVal += other);
}

operator+=()

FitAmpPairCovariance & FitAmpPairCovariance::operator+=

(

const FitAmpPairCovariance &

other

)

Definition at line 810 of file FitAmpPairCovariance.cpp.

                                                                                       {
  this->add(other);
  return *this;
}

realName()

std::string FitAmpPairCovariance::realName ( int  listPosition ) const

private

Definition at line 170 of file FitAmpPairCovariance.cpp.

                                                              {
  return _myList->at(listPosition).name() + "_real";
}

reset()

bool FitAmpPairCovariance::reset

(

)

Definition at line 583 of file FitAmpPairCovariance.cpp.

                                {
  clearAll();
  resize();
  return true;
}

resize()

bool FitAmpPairCovariance::resize ( )

private

Definition at line 124 of file FitAmpPairCovariance.cpp.

                                 {
  bool dbThis=false;
  if(dbThis){
    cout << "FitAmpPairCovariance::resize()"
     << " resizing from " << size()
     << " to " << _myList->size() 
     << " with maxSize() = " << maxSize()
     << endl;
  }

  _needToRecalculate=true;

  _size = _myList->size();

  if(size() > maxSize()){
    cout << "FitAmpPairCovariance::resize size() = " << size()
     << " that's too big! won't do it!" << endl;
    return false;
  }

  _sum_x. ResizeTo(size()*2, 1);
  _sum_xy.ResizeTo(size()*2, size()*2);
  _cov2N. ResizeTo(size()*2, size()*2);
  _corr2N.ResizeTo(size()*2, size()*2);

  _cov.   ResizeTo(size(), size());
  _corr.  ResizeTo(size(), size());

  _fractionCov.     ResizeTo(size(), size());
  _fractionSumCov.  ResizeTo(1,1);
  _integralCov.     ResizeTo(1,1);

  _dN_by_d2N.       ResizeTo(size(), size()*2);
  _dFraction_by_dN. ResizeTo(size(), size());

  _dFractionSum_by_dFraction.ResizeTo(1, size());

  _dIntegral_by_dN. ResizeTo(1, size());

  if(dbThis){
    cout << "FitAmpPairCovariance::resize() to "
     << size() << " done." << endl;
  }
  return true;
}

retrieve()

bool FitAmpPairCovariance::retrieve

(

const std::string &

fromDirectory = "."

)

Definition at line 711 of file FitAmpPairCovariance.cpp.

                                                               {
  if(! resize()) return false;
  if(! retrieve_sum_x(inDirectory)) return false;
  if(! retrieve_sum_xy(inDirectory)) return false;
  _needToRecalculate=true;
  return true;
}

retrieve_sum_x()

bool FitAmpPairCovariance::retrieve_sum_x ( const std::string &  inDirectory )

private

Definition at line 719 of file FitAmpPairCovariance.cpp.

                                                                      {
  bool dbThis=false;
  if(size() > maxSize()) return false;
  ifstream is(matrix_x_fname(inDirectory).c_str());
  if(is.bad()) return false;

  if(dbThis) cout << "opened file: "
          << "\n" << matrix_x_fname(inDirectory)
          << endl;
  std::string st1;
  double val;
  std::map<string, double> map_x;

  int counter=0;
  while(is){
    if(dbThis) cout << "reading " << counter << ". line" << endl;
    if(0 == counter){
      is >> st1 >> _Nevents;
      if(dbThis){
    cout << "just read " << st1 << " " << _Nevents << endl;
      }
     }else{
      is >> st1 >> val;
      map_x[st1] = val;
      if(dbThis){
    cout << "just read " << st1 << " " << val << endl;
      }
    }
    counter++;
  }
  is.close();
  if(dbThis)cout << "retrieve_sum_x, filling _sum_x for size " 
         << size() << endl;
  for(unsigned int i=0; i < size(); i++){
    _sum_x(indexReal(i), 0) = map_x[realName(i)];
    _sum_x(indexImag(i), 0) = map_x[imagName(i)];
   if(dbThis){
      cout << "Tried to retrieve from map entry for " << realName(i)
       << " \n got: " << map_x[realName(i)] << endl;
    }
 
  }
  if(dbThis){
    cout << " did so, look at what I got for sum_x" << endl;
    _sum_x.Print();
  }
  _needToRecalculate=true;
  return true;
}

retrieve_sum_xy()

bool FitAmpPairCovariance::retrieve_sum_xy ( const std::string &  inDirectory )

private

Definition at line 768 of file FitAmpPairCovariance.cpp.

                                                                       {
  bool dbThis=false;
  if(size() > maxSize()) return false;
  ifstream is(matrix_xy_fname(inDirectory).c_str());
  if(is.bad()) return false;
  if(dbThis) cout << "opened file: "
          << "\n" << matrix_xy_fname(inDirectory)
          << endl;
  std::string st1, st2;
  double val;
  map< string, map<string, double> > map_xy;

  int counter=0;
  while(is){
    if(0 == counter){
      is >> st1 >> _Nevents;
    }else{
      is >> st1 >> st2 >> val;
      map_xy[st1][st2] = val;
      if(dbThis) cout << "just read " << counter << ". line: "
              << st1 << " " << st2 << " " << val << endl;
    }
    counter++;
  }

  is.close();
  for(unsigned int i=0; i < size(); i++){
    for(unsigned int j=0; j < size(); j++){
      _sum_xy(indexReal(i), indexReal(j)) = map_xy[realName(i)][realName(j)];
      _sum_xy(indexReal(i), indexImag(j)) = map_xy[realName(i)][imagName(j)];
      _sum_xy(indexImag(i), indexReal(j)) = map_xy[imagName(i)][realName(j)];
      _sum_xy(indexImag(i), indexImag(j)) = map_xy[imagName(i)][imagName(j)];
    }
  }
  if(dbThis){
    cout << "got _sum_xy: " << endl;
    _sum_xy.Print();
  }
  _needToRecalculate=true;
  return true;
}

save()

bool FitAmpPairCovariance::save

(

const std::string &

asSubdirOf = "."

)

const

Definition at line 635 of file FitAmpPairCovariance.cpp.

                                                                {
  bool dbThis=false;
  if(! isValid()){
    cout << "WARNING FitAmpPairCovariance::save(" << inDirectory <<")"
     << " I'm not valid, I won't save anything." << endl;
    return false;
  }
  if(dbThis){
    cout << "saving FitAmpPairCovariance in " << inDirectory << "." << endl;
  }
  makeDirectory(inDirectory);
  bool sc=true;
  sc &= save_sum_x(inDirectory);
  sc &= save_sum_xy(inDirectory);

  return true;
}

save_sum_x()

bool FitAmpPairCovariance::save_sum_x ( const std::string &  inDirectory ) const

private

Definition at line 653 of file FitAmpPairCovariance.cpp.

                                                                       {
  bool dbThis=false;
  if(! isValid()){
    cout << "WARNING FitAmpPairCovariance::save_sum_x(" << inDirectory <<")"
     << " I'm not valid, I won't save anything." << endl;
    return false;
  }
  ofstream os(matrix_x_fname(inDirectory).c_str());
  os << "N " << _Nevents << endl;
  if(dbThis){
    cout << "about to save: ";
    _sum_x.Print();
  }
  for(unsigned int i=0; i < _myList->size(); i++){

    os << setprecision(20) << realName(i) << " "
       << _sum_x(indexReal(i), 0) << endl;

    os << setprecision(20) << imagName(i) << " "
       << _sum_x(indexImag(i), 0) << endl;
  }
  os.close();
  return true;
}

save_sum_xy()

bool FitAmpPairCovariance::save_sum_xy ( const std::string &  inDirectory ) const

private

Definition at line 678 of file FitAmpPairCovariance.cpp.

                                                                        {
  bool dbThis=false;
  if(! isValid()){
    cout << "WARNING FitAmpPairCovariance::save_sum_xy(" << inDirectory <<")"
     << " I'm not valid, I won't save anything." << endl;
    return false;
  }
  ofstream os(matrix_xy_fname(inDirectory).c_str());
  os << "N " << _Nevents << endl;
  if(dbThis){
    cout << "about to save: ";
    _sum_xy.Print();
  }
  for(unsigned int i=0; i < _myList->size(); i++){
    for(unsigned int j=0; j < _myList->size(); j++){

      os << setprecision(20)<< realName(i) << " " << realName(j) << " "
     << _sum_xy(indexReal(i), indexReal(j)) << endl;

      os << setprecision(20)<< realName(i) << " " << imagName(j) << " "
     << _sum_xy(indexReal(i), indexImag(j)) << endl;

      os << setprecision(20)<< imagName(i) << " " << realName(j) << " "
     << _sum_xy(indexImag(i), indexReal(j)) << endl;

      os << setprecision(20)<< imagName(i) << " " << imagName(j) << " "
     << _sum_xy(indexImag(i), indexImag(j)) << endl;
    }
  }
  os.close();
  return true;
}

size()

unsigned int FitAmpPairCovariance::size

(

)

const

Definition at line 121 of file FitAmpPairCovariance.cpp.

                                             {
  return _size;
}

Member Data Documentation

_corr

TMatrixTSym<float> FitAmpPairCovariance::_corr

private

Definition at line 22 of file FitAmpPairCovariance.h.

_corr2N

TMatrixTSym<float> FitAmpPairCovariance::_corr2N

private

Definition at line 22 of file FitAmpPairCovariance.h.

_cov

TMatrixTSym<float> FitAmpPairCovariance::_cov

private

Definition at line 22 of file FitAmpPairCovariance.h.

_cov2N

TMatrixTSym<float> FitAmpPairCovariance::_cov2N

private

Definition at line 22 of file FitAmpPairCovariance.h.

_dFraction_by_dN

TMatrixT<float> FitAmpPairCovariance::_dFraction_by_dN

private

Definition at line 28 of file FitAmpPairCovariance.h.

_dFractionSum_by_dFraction

TMatrixT<float> FitAmpPairCovariance::_dFractionSum_by_dFraction

private

Definition at line 29 of file FitAmpPairCovariance.h.

_dIntegral_by_dN

TMatrixT<float> FitAmpPairCovariance::_dIntegral_by_dN

private

Definition at line 30 of file FitAmpPairCovariance.h.

_dN_by_d2N

TMatrixT<float> FitAmpPairCovariance::_dN_by_d2N

private

Definition at line 27 of file FitAmpPairCovariance.h.

_fractionCov

TMatrixTSym<float> FitAmpPairCovariance::_fractionCov

private

Definition at line 23 of file FitAmpPairCovariance.h.

_fractionSumCov

TMatrixTSym<float> FitAmpPairCovariance::_fractionSumCov

private

Definition at line 24 of file FitAmpPairCovariance.h.

_integralCov

TMatrixTSym<float> FitAmpPairCovariance::_integralCov

private

Definition at line 25 of file FitAmpPairCovariance.h.

_maxSize

unsigned int FitAmpPairCovariance::_maxSize = 1500

staticprivate

Definition at line 12 of file FitAmpPairCovariance.h.

_myList

const FitAmpPairList* FitAmpPairCovariance::_myList

private

Definition at line 14 of file FitAmpPairCovariance.h.

_needToRecalculate

bool FitAmpPairCovariance::_needToRecalculate

private

Definition at line 17 of file FitAmpPairCovariance.h.

_Nevents

unsigned int FitAmpPairCovariance::_Nevents

private

Definition at line 15 of file FitAmpPairCovariance.h.

_size

unsigned int FitAmpPairCovariance::_size

private

Definition at line 16 of file FitAmpPairCovariance.h.

_sum_x

TMatrixT<float> FitAmpPairCovariance::_sum_x

private

Definition at line 19 of file FitAmpPairCovariance.h.

_sum_xy

TMatrixTSym<float> FitAmpPairCovariance::_sum_xy

private

Definition at line 20 of file FitAmpPairCovariance.h.

---

The documentation for this class was generated from the following files:

• Mint/FitAmpPairCovariance.h
• src/Mojito/FitAmplitude/FitAmpPairCovariance.cpp