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 132 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 180 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 90 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.

                                     {
  std::cout << std::format(" Print fitted param from TileFilterResult:  Nparam={}, chisq={}\n", m_nParam, m_chi2);
  for (int ipar = 0; ipar < m_nParam; ipar++) {
    if (ipar == 0) {
      std::cout << std::format(" i={}, kcr=P, A={:>5.2g} +-{}\n", ipar, m_fitParam[ipar], m_fitErr[ipar]);
    } else {
      std::cout << std::format(" i={}, kcr={}, A={} +-{}\n", ipar, m_vCross[ipar - 1], m_fitParam[ipar], m_fitErr[ipar]);
    }
  }
  return;
}

◆ snapShot()

void TileFilterResult::snapShot ( int imode )

Definition at line 100 of file TileFilterResult.cxx.

                                         {
  // This print a short snapshot of the FilterResult state.
  std::cout << std::format(" SnapShot: imode={}.  Nparam={}, chisq={}, iFitIndex{}, Vcross=",
                           imode,
                           m_nParam,
                           m_chi2,
                           m_iFitIndex);
  int Namp = m_nParam - 1;
  for (int jamp = 0; jamp < Namp; jamp++) {
    std::cout << std::format(" {}", 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::format("{:>5.1g}+-{}", 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::format(" {:>4.3g}", 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:

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TileFilterResult()

◆ ~TileFilterResult()

Member Function Documentation

◆ addCross()

◆ dropCross()

◆ getChi2Ref()

◆ getDigRef()

◆ getErrRef()

◆ getFitIndexRef()

◆ getInTime()

◆ getNparRef()

◆ getParamRef()

◆ getResidRef()

◆ getSigDig()

◆ getVcrossRef()

◆ printFitParam()

◆ snapShot()

Member Data Documentation

◆ m_chi2

◆ m_debug

◆ m_digits

◆ m_fitErr

◆ m_fitParam

◆ m_iFitIndex

◆ m_nParam

◆ m_nPileup

◆ m_residuals

◆ m_sigDigit

◆ m_vCross