TileFilterResult Class Reference

Auxiliary class for TileRawChannelMakerManyAmps. More...

#include <TileFilterResult.h>

Collaboration diagram for TileFilterResult:

Public Member Functions
	TileFilterResult (std::vector< float > &dig, double sigma)
	~TileFilterResult ()
double	getSigDig () const
CLHEP::HepVector &	getDigRef ()
std::vector< int > &	getVcrossRef ()
int &	getFitIndexRef ()
int &	getNparRef ()
CLHEP::HepVector &	getParamRef ()
CLHEP::HepVector &	getErrRef ()
CLHEP::HepVector &	getResidRef ()
double &	getChi2Ref ()
int	addCross (int kcrIndex)
int	dropCross (int kcrIndex)
void	printFitParam ()
double	getInTime (double &, double &, double &, double &, double &)
void	snapShot (int imode)

Private Attributes
bool	m_debug
CLHEP::HepVector	m_digits
double	m_sigDigit
int	m_nParam
int	m_nPileup
int	m_iFitIndex
std::vector< int >	m_vCross
CLHEP::HepVector	m_fitParam
CLHEP::HepVector	m_fitErr
CLHEP::HepVector	m_residuals
double	m_chi2

Detailed Description

Auxiliary class for TileRawChannelMakerManyAmps.

Definition at line 21 of file TileFilterResult.h.

Constructor & Destructor Documentation

TileFilterResult()

TileFilterResult::TileFilterResult	(	std::vector< float > &	dig,
		double	sigma )

Definition at line 18 of file TileFilterResult.cxx.

                                                                    {
  m_debug = false;
  int ND = dig.size();
  CLHEP::HepVector digtem(ND);
  m_digits = digtem;
  for (int id = 0; id < ND; id++) {
    m_digits[id] = dig[id];
  }
  //int nrow = digits.num_row();
  //int nrowtem = digtem.num_row();
  m_sigDigit = sig;
  // Initialize some of the variables that will be set later.
  m_nParam = 1;
  m_nPileup = 0;
  m_vCross.clear();
  m_chi2 = 999.;
  m_iFitIndex = -1;
}

~TileFilterResult()

TileFilterResult::~TileFilterResult

(

)

Definition at line 37 of file TileFilterResult.cxx.

                                    {
  return;
}

Member Function Documentation

addCross()

int TileFilterResult::addCross

(

int

kcrIndex

)

Definition at line 133 of file TileFilterResult.cxx.

                                           {
 
  int iret = -1;
 
  //Check that kcrIndex is not already in list.
  bool ldup = false;
  if (m_nParam > 1) {
    int Namp = m_nParam - 1;
    for (int icr = 0; icr < Namp; icr++) {
      if (kcrIndex == m_vCross[icr]) {
        ldup = true;
        if (m_debug) {
          std::cout << " TileFilterResult.addCross: kcrIndex=" << kcrIndex << " is already in crossing list: Kcross =";
          for (int j = 0; j < Namp; j++) {
            std::cout << " " << m_vCross[j];
          }
          std::cout << std::endl;
        }
        if (ldup) break;
      }
    } // end for loop
  } // end "if(Nparam>1)"
  if (ldup) {
    iret = 1;
    return iret;
  }
 
  // Add the new crossing.
  iret = 0;
  m_vCross.push_back(kcrIndex);
  std::sort(m_vCross.begin(), m_vCross.end());
  m_nParam = m_nParam + 1;
 
  // Since we have a new Vcross configuration, iFitIndex is not yet defined for it.
  m_iFitIndex = -1;
 
  if (m_debug) {
    int Namp = m_nParam - 1;
    std::cout << " TileFilterResult.addCross.  Exit with Nparam=" << m_nParam << ", Vcross=";
    for (int icr = 0; icr < Namp; icr++) {
      std::cout << " " << std::setw(3) << m_vCross[icr];
    }
    std::cout << std::endl;
  }
 
  return iret;
}

dropCross()

int TileFilterResult::dropCross

(

int

kcrIndex

)

Definition at line 181 of file TileFilterResult.cxx.

                                         {
  // Drop the crossing stored in amplitude # iamp.
  int iret = -1;
  int Namp = m_nParam - 1;
  for (int iamp = 1; iamp < Namp; iamp++) {
    if (m_vCross[iamp] == idrop) {
      // Erase the crossing.
      iret = 0;
      m_vCross.erase(m_vCross.begin() + iamp);
      //  std::sort(Vcross.begin(), Vcross.end() );
      m_nParam = m_nParam - 1;
    }
  }
  if (iret != 0) {
    std::cout << "error in TileFilterResult.dropCross:  idrop=" << idrop << "but vcross =";
    for (int iamp = 0; iamp < Namp; iamp++) {
      std::cout << " " << m_vCross[iamp];
    }
    std::cout << std::endl;
    return iret;
  }
 
  // Since we have a new Vcross configuration, iFitIndex is not yet defined for it.
  m_iFitIndex = -1;
 
  if (m_debug) {
    int Namp = m_nParam - 1;
    std::cout << " TileFilterResult.dropCross.  Exit with Nparam=" << m_nParam << ", Vcross=";
    for (int icr = 0; icr < Namp; icr++) {
      std::cout << " " << std::setw(3) << m_vCross[icr];
    }
    std::cout << std::endl;
  }
 
  return iret;
}

getChi2Ref()

double & TileFilterResult::getChi2Ref

(

)

Definition at line 74 of file TileFilterResult.cxx.

                                     {
  return m_chi2;
}

getDigRef()

CLHEP::HepVector & TileFilterResult::getDigRef

(

)

Definition at line 46 of file TileFilterResult.cxx.

                                            {
  return m_digits;
}

getErrRef()

CLHEP::HepVector & TileFilterResult::getErrRef

(

)

Definition at line 66 of file TileFilterResult.cxx.

                                            {
  return m_fitErr;
}

getFitIndexRef()

int & TileFilterResult::getFitIndexRef

(

)

Definition at line 54 of file TileFilterResult.cxx.

                                      {
  return m_iFitIndex;
}

getInTime()

double TileFilterResult::getInTime	(	double &	amp,
		double &	err,
		double &	ped,
		double &	chi2,
		double &	t )

Definition at line 93 of file TileFilterResult.cxx.

                                                                                                 {
  amp = m_fitParam[1];
  err = m_fitErr[1];
  ped = m_fitParam[0];
  chi2 = m_chi2;
  t = 0.;
 
  return m_chi2;
}

getNparRef()

int & TileFilterResult::getNparRef

(

)

Definition at line 58 of file TileFilterResult.cxx.

                                  {
  return m_nParam;
}

getParamRef()

CLHEP::HepVector & TileFilterResult::getParamRef

(

)

Definition at line 62 of file TileFilterResult.cxx.

                                              {
  return m_fitParam;
}

getResidRef()

CLHEP::HepVector & TileFilterResult::getResidRef

(

)

Definition at line 70 of file TileFilterResult.cxx.

                                              {
  return m_residuals;
}

getSigDig()

double TileFilterResult::getSigDig

(

)

const

Definition at line 41 of file TileFilterResult.cxx.

                                         {
  return m_sigDigit;
}

getVcrossRef()

std::vector< int > & TileFilterResult::getVcrossRef

(

)

Definition at line 50 of file TileFilterResult.cxx.

                                               {
  return m_vCross;
}

printFitParam()

void TileFilterResult::printFitParam

(

)

Definition at line 78 of file TileFilterResult.cxx.

                                     {
  boost::io::ios_base_all_saver coutsave(std::cout);
  std::cout << " Print fitted param from TileFilterResult:  Nparam=" << m_nParam << ", chisq=" << m_chi2 << std::endl;
  for (int ipar = 0; ipar < m_nParam; ipar++) {
    if (ipar == 0) {
      std::cout << " i=" << ipar << ", kcr=P" << ", A=" << std::setw(5) << std::setprecision(2) << m_fitParam[ipar]
          << " +-" << m_fitErr[ipar] << std::endl;
    } else {
      std::cout << " i=" << ipar << ", kcr=" << m_vCross[ipar - 1] << ", A=" << m_fitParam[ipar] << " +-"
          << m_fitErr[ipar] << std::endl;
    }
  }
  return;
}

snapShot()

void TileFilterResult::snapShot

(

int

imode

)

Definition at line 103 of file TileFilterResult.cxx.

                                         {
  boost::io::ios_base_all_saver coutsave(std::cout);
  // This print a short snapshot of the FilterResult state.
  std::cout << " SnapShot: imode=" << imode << ".  Nparam=" << m_nParam << ", chisq=" << m_chi2 << ", iFitIndex"
      << m_iFitIndex << ", Vcross=";
  int Namp = m_nParam - 1;
  for (int jamp = 0; jamp < Namp; jamp++) {
    std::cout << " " << m_vCross[jamp];
  }
  std::cout << std::endl;
  if (m_iFitIndex < 0) return;
 
  if (imode > 0) {
    std::cout << "   FitParam=";
    for (int ipar = 0; ipar < m_nParam; ipar++) {
      std::cout << std::setw(5) << std::setprecision(1) << m_fitParam[ipar] << "+-" << m_fitErr[ipar];
      if (ipar < m_nParam - 1) std::cout << ", ";
    }
    std::cout << std::endl;
  }
  if (imode > 1) {
    int Ndig = m_digits.num_row();
    std::cout << "   Residuals=";
    for (int idig = 0; idig < Ndig; idig++) {
      std::cout << " " << std::setw(4) << std::setprecision(3) << m_residuals[idig];
    }
    std::cout << std::endl;
  }
}

Member Data Documentation

m_chi2

double TileFilterResult::m_chi2

private

Definition at line 62 of file TileFilterResult.h.

m_debug

bool TileFilterResult::m_debug

private

Definition at line 48 of file TileFilterResult.h.

m_digits

CLHEP::HepVector TileFilterResult::m_digits

private

Definition at line 50 of file TileFilterResult.h.

m_fitErr

CLHEP::HepVector TileFilterResult::m_fitErr

private

Definition at line 59 of file TileFilterResult.h.

m_fitParam

CLHEP::HepVector TileFilterResult::m_fitParam

private

Definition at line 58 of file TileFilterResult.h.

m_iFitIndex

int TileFilterResult::m_iFitIndex

private

Definition at line 55 of file TileFilterResult.h.

m_nParam

int TileFilterResult::m_nParam

private

Definition at line 53 of file TileFilterResult.h.

m_nPileup

int TileFilterResult::m_nPileup

private

Definition at line 54 of file TileFilterResult.h.

m_residuals

CLHEP::HepVector TileFilterResult::m_residuals

private

Definition at line 61 of file TileFilterResult.h.

m_sigDigit

double TileFilterResult::m_sigDigit

private

Definition at line 51 of file TileFilterResult.h.

m_vCross

std::vector<int> TileFilterResult::m_vCross

private

Definition at line 56 of file TileFilterResult.h.

---

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

• TileFilterResult.h
• TileFilterResult.cxx