MINT::MinuitParameterSet Class Reference

#include <MinuitParameterSet.h>

Inheritance diagram for MINT::MinuitParameterSet:

Public Types
typedef TMatrixTSym< double >	CovMatrix

Public Member Functions
	MinuitParameterSet ()
	MinuitParameterSet (const MinuitParameterSet &other)
MinuitParameterSet	getFloating ()
bool	add (IMinuitParameter *parPtr)
bool	unregister (IMinuitParameter *patPtr)
unsigned int	size () const
IMinuitParameter *	getParPtr (unsigned int i)
const IMinuitParameter *	getParPtr (unsigned int i) const
IMinuitParameter *	getParPtr (std::string name)
const IMinuitParameter *	getParPtr (std::string name) const
int	findParPtr (std::string name)
void	deleteListAndObjects ()
void	deleteListKeepObjects ()
TNtupleD *	makeNewNtpForOwner (TFile *&ntpFile) const
std::string	ntpNames () const
void	fillNtp (TFile *&ntpFile, TNtupleD *&ntp) const
void	print (std::ostream &os=std::cout) const
void	printVariable (std::ostream &os=std::cout) const
void	printResultVsInput (std::ostream &os=std::cout) const
bool	setCovMatrix (const CovMatrix &)
	Set the covariance matrix. More...
const CovMatrix &	covMatrix () const
	Get the covariance matrix. More...

Static Public Member Functions
static MinuitParameterSet *	getDefaultSet ()

Protected Member Functions
bool	addToEnd (IMinuitParameter *parPtr)
bool	setAllIndices ()

Static Protected Member Functions
static std::string	prtToNtpName (const std::string &s_in)

Protected Attributes
std::vector< IMinuitParameter * >	_parPtrList
CovMatrix	_covMatrix

Static Protected Attributes
static MinuitParameterSet *	_defaultMinuitParameterSet =0
static const char	ntpNameChars [] = { 'p', 'm', 's', '_', '_', '_', '_', '_', '_', '\0'}
static const char	prtNameChars [] = { '+', '-', '*', '>', ',', '(', ')', '[', ']', '\0'}

Detailed Description

Definition at line 20 of file MinuitParameterSet.h.

Member Typedef Documentation

CovMatrix

typedef TMatrixTSym<double> MINT::MinuitParameterSet::CovMatrix

Definition at line 22 of file MinuitParameterSet.h.

Constructor & Destructor Documentation

MinuitParameterSet() [1/2]

MinuitParameterSet::MinuitParameterSet

(

)

Definition at line 32 of file MinuitParameterSet.cpp.

{}

MinuitParameterSet() [2/2]

MinuitParameterSet::MinuitParameterSet

(

const MinuitParameterSet &

other

)

Definition at line 33 of file MinuitParameterSet.cpp.

  : _parPtrList(other._parPtrList)
{}

Member Function Documentation

add()

bool MinuitParameterSet::add

(

IMinuitParameter *

parPtr

)

Definition at line 67 of file MinuitParameterSet.cpp.

                                                    {
  return addToEnd(parPtr);
}

addToEnd()

bool MinuitParameterSet::addToEnd ( IMinuitParameter *  parPtr )

protected

Definition at line 47 of file MinuitParameterSet.cpp.

                                                         {
  bool dbThis=false;
  bool success=true;
  if(0 == parPtr) return false;
  if(dbThis){
    cout << "MinuitParameterSet::add "
     << " adding parPtr " << parPtr->name() << endl;
  }
  if(find(_parPtrList.begin(), _parPtrList.end(), parPtr) != _parPtrList.end()){
    if(dbThis){
      cout << "MinuitParameterSet::add "
       << " already have parPtr " << parPtr->name() << endl;
    }
    return true ;
  }
  _parPtrList.push_back(parPtr);
  success &= parPtr->setParSet(this);
  success &= parPtr->setParSetIndex(_parPtrList.size()-1);
  return success;
}

covMatrix()

const MinuitParameterSet::CovMatrix & MinuitParameterSet::covMatrix

(

)

const

Get the covariance matrix.

Definition at line 255 of file MinuitParameterSet.cpp.

                                                                       {
  return _covMatrix ;
}

deleteListAndObjects()

void MinuitParameterSet::deleteListAndObjects

(

)

Definition at line 110 of file MinuitParameterSet.cpp.

                                             {
  for(std::vector<IMinuitParameter*>::iterator it = _parPtrList.begin();
      it != _parPtrList.end(); it++){
    delete (*it);
  }
  _parPtrList.clear();
}

deleteListKeepObjects()

void MinuitParameterSet::deleteListKeepObjects

(

)

Definition at line 118 of file MinuitParameterSet.cpp.

                                              {
  _parPtrList.clear();
}

fillNtp()

void MinuitParameterSet::fillNtp	(	TFile *&	ntpFile,
		TNtupleD *&	ntp
	)		const

Definition at line 179 of file MinuitParameterSet.cpp.

                                                                     {
  bool dbThis=false;
  if(dbThis){
    cout << "hello from MinuitParameterSet::fillNtp " << endl;
    cout << " you called me with: " << ntp << endl;
  }
  
  if(0 == ntp) ntp = makeNewNtpForOwner(ntpFile);
  if(dbThis) cout << "got an ntuple ptr: " << ntp << endl;
  Double_t* array = new Double_t[size()*4];
  if(dbThis) cout << "made array" << endl;
  int n=0;
  for(unsigned int i=0; i < size(); i++){
    if(0 == getParPtr(i)) continue;
    if(0 != getParPtr(i)->iFixInit()) continue;
    array[n++] = ((FitParameter*)getParPtr(i))->blindedMean();
    array[n++] = getParPtr(i)->meanInit();
    array[n++] = getParPtr(i)->err();
    Double_t pull=-9999.0;
    if(getParPtr(i)->err() > 0){
      pull = (((FitParameter*)getParPtr(i))->blindedMean() - getParPtr(i)->meanInit())/getParPtr(i)->err();
    }
    array[n++] = pull;
  }
  if(dbThis){
    cout << "checking how ntuple behaves with GetEntries()" << endl;
    cout << " ntp->GetEntries() = " << ntp->GetEntries() << endl;
  }
  if(dbThis)cout << "filled the array, now putting this into the ntp" << endl;
  ntp->Fill(array);
  if(dbThis)cout << "filled into the ntuple, now deleting array" << endl;
  delete[] array;
  if(dbThis)cout << "done. returning" << endl;
  return;
}

findParPtr()

int MINT::MinuitParameterSet::findParPtr ( std::string  name )

inline

Definition at line 67 of file MinuitParameterSet.h.

                                {
      for(unsigned int i=0; i <this->size(); i++){
          if(0 == this->getParPtr(i)) continue;
          if(A_is_in_B(name,this->getParPtr(i)->name()))return i;
      }
      return -1;
  }

getDefaultSet()

MinuitParameterSet * MinuitParameterSet::getDefaultSet ( )

static

Definition at line 25 of file MinuitParameterSet.cpp.

                                                     {
  if(0 == _defaultMinuitParameterSet){
    _defaultMinuitParameterSet = new MinuitParameterSet();
  }
  return _defaultMinuitParameterSet;
}

getFloating()

MinuitParameterSet MinuitParameterSet::getFloating

(

)

Definition at line 37 of file MinuitParameterSet.cpp.

                                                  {
  MinuitParameterSet floating;
  for(unsigned int i=0; i< size(); i++){
    if(0 != getParPtr(i)){
      floating.add(getParPtr(i));
    }
  }
  return floating;
}

getParPtr() [1/4]

IMinuitParameter * MinuitParameterSet::getParPtr

(

unsigned int

i

)

Definition at line 100 of file MinuitParameterSet.cpp.

                                                             {
  if(i >= _parPtrList.size()) return 0;
  return _parPtrList[i];
}

getParPtr() [2/4]

const IMinuitParameter * MinuitParameterSet::getParPtr

(

unsigned int

i

)

const

Definition at line 105 of file MinuitParameterSet.cpp.

                                                                         {
  if(i >= _parPtrList.size()) return 0;
  return _parPtrList[i];
}

getParPtr() [3/4]

IMinuitParameter* MINT::MinuitParameterSet::getParPtr ( std::string  name )

inline

Definition at line 52 of file MinuitParameterSet.h.

                                             {
      for(unsigned int i=0; i <this->size(); i++){
          if(0 == this->getParPtr(i)) continue;
          if(A_is_in_B(name,this->getParPtr(i)->name()))return this->getParPtr(i);
      }
      return 0;
  }

getParPtr() [4/4]

const IMinuitParameter* MINT::MinuitParameterSet::getParPtr ( std::string  name ) const

inline

Definition at line 59 of file MinuitParameterSet.h.

                                                          {
      for(unsigned int i=0; i <this->size(); i++){
          if(0 == this->getParPtr(i)) continue;
          if(A_is_in_B(name,this->getParPtr(i)->name()))return this->getParPtr(i);
      }
      return 0;
  }

makeNewNtpForOwner()

TNtupleD * MinuitParameterSet::makeNewNtpForOwner

(

TFile *&

ntpFile

)

const

Definition at line 122 of file MinuitParameterSet.cpp.

                                                                     {
  bool dbThis=false;
  if(dbThis){
    cout << "Hello from MinuitParameterSet::makeNewNtpForOwner" << endl;
  }
  gDirectory->cd();
  if(0 == ntpFile) ntpFile = new TFile("MinuitParameterSetFile.root"
                       , "RECREATE");
  ntpFile->cd();
  TNtupleD* ntp = new TNtupleD("MinuitParameterSetNtp"
                   , "MinuitParameterSetNtp"
                   , ntpNames().c_str());
  return ntp;
}

ntpNames()

std::string MinuitParameterSet::ntpNames

(

)

const

Definition at line 152 of file MinuitParameterSet.cpp.

                                            {
  bool dbThis=false;
  std::string str="";
  int n=0;
  for(unsigned int i=0; i < size(); i++){
    if(0 == getParPtr(i)) continue;
    if(0 != getParPtr(i)->iFixInit()) continue;
    std::string name = //"p" + anythingToString(i) + "_" + 
    prtToNtpName(getParPtr(i)->name());
    name = ((string) (TString(name).ReplaceAll("GS","")));
    name = ((string) (TString(name).ReplaceAll("Lass","")));
    name = ((string) (TString(name).ReplaceAll("RhoOmega","")));
    if(str != "") str += ":";
    str += (name + "_mean" + ":"); n++;
    str += (name + "_init" + ":"); n++;
    str += (name + "_err:");       n++;
    str += (name + "_pull");       n++;
  }
  if(dbThis){
    cout << "MinuitParameterSet::ntpNames():"
     << " made " << n << " names " << str << endl;
    cout << " this becomes: prtToNtpName(str) " 
     << prtToNtpName(str) << endl;
  }
  //return str;
  return prtToNtpName(str);
}

print()

void MinuitParameterSet::print

(

std::ostream &

os = std::cout

)

const

Definition at line 216 of file MinuitParameterSet.cpp.

                                                  {
  for(unsigned int i=0; i < size(); i++){
    if(0 == getParPtr(i)) continue;
    os << '\n';
    getParPtr(i)->print(os);
  }
}

printResultVsInput()

void MinuitParameterSet::printResultVsInput

(

std::ostream &

os = std::cout

)

const

Definition at line 231 of file MinuitParameterSet.cpp.

                                                               {
  os << "\n ======== MinuitParameterSet::(FitResult - StartValue)/error ========\n";
  for(unsigned int i=0; i < this->size(); i++){
    if(0 == getParPtr(i)) continue;
    if(getParPtr(i)->iFixInit() == 1) continue;
    os << '\n';
    getParPtr(i)->printResultVsInput(os);
  }
  os << "\n ===================================================================="
     << endl;
}

printVariable()

void MinuitParameterSet::printVariable

(

std::ostream &

os = std::cout

)

const

Definition at line 223 of file MinuitParameterSet.cpp.

                                                          {
  for(unsigned int i=0; i < size(); i++){
    if(0 == getParPtr(i)) continue;
    if(getParPtr(i)->iFixInit() == 1) continue;
    os << '\n';
    getParPtr(i)->print(os);
  }
}

prtToNtpName()

std::string MinuitParameterSet::prtToNtpName ( const std::string &  s_in )

staticprotected

Definition at line 137 of file MinuitParameterSet.cpp.

                                                               {
  std::string s(s_in);
  for(unsigned int i=1; i < s.size(); i++){
    // cout << "... " << s[i-1] << s[i];
    if(s[i-1] == '-' && s[i] == '>'){
      s[i-1] = 't'; s[i] = 'o';
    }
  }

  for(int i=0; prtNameChars[i] != '\0'; i++){
    replace(s.begin(), s.end(), prtNameChars[i], ntpNameChars[i]);
  }
  return s;
}

setAllIndices()

bool MinuitParameterSet::setAllIndices ( )

protected

Definition at line 88 of file MinuitParameterSet.cpp.

                                      {
  bool success=true;
  for(unsigned int i=0; i < _parPtrList.size(); i++){
    success &= _parPtrList[i]->setParSetIndex(i);
  }
  return success;
}

setCovMatrix()

bool MinuitParameterSet::setCovMatrix

(

const CovMatrix &

covMatrix

)

Set the covariance matrix.

Definition at line 243 of file MinuitParameterSet.cpp.

                                                                {
  if(covMatrix.GetNcols() != size() || covMatrix.GetNrows() != size()){
    cout << "WARNING: MiniutParameterSet::setCovMatrix: given matrix has the wrong size!" << endl ;
    cout << "N. parameters: " << size() << ", matrix is " << covMatrix.GetNcols() << "x" 
     << covMatrix.GetNrows() << endl ;
    return false ;
  }
  _covMatrix.ResizeTo(covMatrix) ;
  _covMatrix = covMatrix ;
  return true ;
}

size()

unsigned int MinuitParameterSet::size

(

)

const

Definition at line 96 of file MinuitParameterSet.cpp.

                                          {
  return _parPtrList.size();
}

unregister()

bool MinuitParameterSet::unregister

(

IMinuitParameter *

patPtr

)

Definition at line 71 of file MinuitParameterSet.cpp.

                                                           {

  if(_parPtrList.end() == find(_parPtrList.begin(), _parPtrList.end(), parPtr)){
    cout << "WARNING in MinuitParameterSet::unregister"
     << " parPtr you want to unregister is not part of this list!"
     << endl;
    return false;
  }

  parPtr->setParSetIndex(-9999);
  _parPtrList.erase(remove(_parPtrList.begin(), _parPtrList.end(), parPtr)
            , _parPtrList.end()
            );
  setAllIndices();
  return true;
}

Member Data Documentation

_covMatrix

CovMatrix MINT::MinuitParameterSet::_covMatrix

protected

Definition at line 34 of file MinuitParameterSet.h.

_defaultMinuitParameterSet

MinuitParameterSet * MinuitParameterSet::_defaultMinuitParameterSet =0

staticprotected

Definition at line 25 of file MinuitParameterSet.h.

_parPtrList

std::vector<IMinuitParameter*> MINT::MinuitParameterSet::_parPtrList

protected

Definition at line 24 of file MinuitParameterSet.h.

ntpNameChars

const char MinuitParameterSet::ntpNameChars = { 'p', 'm', 's', '_', '_', '_', '_', '_', '_', '\0'}

staticprotected

Definition at line 27 of file MinuitParameterSet.h.

prtNameChars

const char MinuitParameterSet::prtNameChars = { '+', '-', '*', '>', ',', '(', ')', '[', ']', '\0'}

staticprotected

Definition at line 28 of file MinuitParameterSet.h.

---

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

• Mint/MinuitParameterSet.h
• src/Mint/Minimiser/MinuitParameterSet.cpp