MuonResidualsTwoBin Class Reference

#include <MuonResidualsTwoBin.h>

Public Member Functions
double	antisym (int parNum)
void	correctBField ()
void	eraseNotSelectedResiduals ()
double	errorerror (int parNum)
void	fiducialCuts ()
void	fill (char charge, double *residual)
bool	fit (Alignable *ali)
void	fix (int parNum, bool value=true)
bool	fixed (int parNum)
double	loglikelihood ()
double	mean (int which, double truncate)
double	median (int which)
	MuonResidualsTwoBin (bool twoBin, MuonResidualsFitter *pos, MuonResidualsFitter *neg)
int	ndata ()
int	npar ()
long	numResidualsNeg () const
long	numResidualsPos () const
double	numsegments ()
double	plot (std::string name, TFileDirectory *dir, Alignable *ali)
void	plotsimple (std::string name, TFileDirectory *dir, int which, double multiplier)
void	plotweighted (std::string name, TFileDirectory *dir, int which, int whichredchi2, double multiplier)
void	read (FILE *file, int which=0)
int	residualsModel () const
std::vector< double * > ::const_iterator	residualsNeg_begin () const
std::vector< double * > ::const_iterator	residualsNeg_end () const
std::vector< bool >::const_iterator	residualsNeg_ok_begin () const
std::vector< bool >::const_iterator	residualsNeg_ok_end () const
std::vector< double * > ::const_iterator	residualsPos_begin () const
std::vector< double * > ::const_iterator	residualsPos_end () const
std::vector< bool >::const_iterator	residualsPos_ok_begin () const
std::vector< bool >::const_iterator	residualsPos_ok_end () const
void	selectPeakResiduals (double nsigma, int nvar, int *vars)
void	setPrintLevel (int printLevel) const
void	setStrategy (int strategy) const
double	stdev (int which, double truncate)
double	sumofweights ()
int	type () const
int	useRes () const
double	value (int parNum)
double	wmean (int which, int whichredchi2, double truncate)
void	write (FILE *file, int which=0)
	~MuonResidualsTwoBin ()

Protected Attributes
MuonResidualsFitter *	m_neg
MuonResidualsFitter *	m_pos
bool	m_twoBin

Detailed Description

Date:

2011/04/15 21:51:13

Revision:

1.10

Author

J. Pivarski - Texas A&M University pivar.nosp@m.ski@.nosp@m.physi.nosp@m.cs.t.nosp@m.amu.e.nosp@m.du

Definition at line 15 of file MuonResidualsTwoBin.h.

Constructor & Destructor Documentation

MuonResidualsTwoBin::MuonResidualsTwoBin ( bool  twoBin, MuonResidualsFitter *  pos, MuonResidualsFitter *  neg  )

inline

Definition at line 17 of file MuonResidualsTwoBin.h.

      : m_twoBin(twoBin), m_pos(pos), m_neg(neg){};

MuonResidualsTwoBin::~MuonResidualsTwoBin ( )

inline

Definition at line 19 of file MuonResidualsTwoBin.h.

References m_neg, and m_pos.

                         {
    if (m_pos != nullptr)
      delete m_pos;
    if (m_neg != nullptr)
      delete m_neg;
  };

Member Function Documentation

double MuonResidualsTwoBin::antisym ( int  parNum )

inline

Definition at line 78 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, and MuonResidualsFitter::value().

{ return (m_twoBin ? ((m_pos->value(parNum) - m_neg->value(parNum)) / 2.) : 0.); };

void MuonResidualsTwoBin::correctBField ( )

inline

Definition at line 240 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::correctBField(), and m_pos.

Referenced by MuonAlignmentFromReference::correctBField().

                       {
    m_pos->correctBField();
    //if (m_twoBin) m_neg->correctBField();
  };

void MuonResidualsTwoBin::eraseNotSelectedResiduals ( )

inline

Definition at line 247 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::eraseNotSelectedResiduals(), m_neg, m_pos, and m_twoBin.

Referenced by MuonAlignmentFromReference::eraseNotSelectedResiduals().

                                   {
    if (m_twoBin) {
      m_pos->eraseNotSelectedResiduals();
      m_neg->eraseNotSelectedResiduals();
    } else {
      m_pos->eraseNotSelectedResiduals();
    }
  }

double MuonResidualsTwoBin::errorerror ( int  parNum )

inline

Definition at line 74 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::errorerror(), m_neg, m_pos, m_twoBin, funct::pow(), and mathSSE::sqrt().

                                {
    return (m_twoBin ? (sqrt(pow(m_pos->errorerror(parNum), 2.) + pow(m_neg->errorerror(parNum), 2.)) / 2.)
                     : m_pos->errorerror(parNum));
  };

void MuonResidualsTwoBin::fiducialCuts ( )

inline

Definition at line 245 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::fiducialCuts(), and m_pos.

Referenced by MuonAlignmentFromReference::fiducialCuts().

{ m_pos->fiducialCuts(); };

void MuonResidualsTwoBin::fill ( char  charge, double *  residual  )

inline

Definition at line 63 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::fill(), m_neg, m_pos, and m_twoBin.

                                           {
    if (!m_twoBin || charge > 0)
      m_pos->fill(residual);
    else
      m_neg->fill(residual);
  };

bool MuonResidualsTwoBin::fit ( Alignable *  ali )

inline

Definition at line 70 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::fit(), m_neg, m_pos, and m_twoBin.

Referenced by trackingPlots.Iteration::modules().

{ return (m_twoBin ? (m_pos->fit(ali) && m_neg->fit(ali)) : m_pos->fit(ali)); };

void MuonResidualsTwoBin::fix ( int  parNum, bool  value = true  )

inline

Definition at line 46 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::fix(), m_neg, m_pos, and value().

                                          {
    m_pos->fix(parNum, value);
    m_neg->fix(parNum, value);
  };

bool MuonResidualsTwoBin::fixed ( int  parNum )

inline

Definition at line 51 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::fixed(), m_neg, and m_pos.

{ return m_pos->fixed(parNum) && m_neg->fixed(parNum); };

double MuonResidualsTwoBin::loglikelihood ( )

inline

Definition at line 79 of file MuonResidualsTwoBin.h.

References MuonResidualsFitter::loglikelihood(), m_neg, m_pos, and m_twoBin.

                         {
    return (m_twoBin ? (m_pos->loglikelihood() + m_neg->loglikelihood()) : m_pos->loglikelihood());
  };

double MuonResidualsTwoBin::mean ( int  which, double  truncate  )

inline

Definition at line 134 of file MuonResidualsTwoBin.h.

References m_twoBin, dqmiodumpmetadata::n, alignCSCRings::r, residualsNeg_begin(), residualsNeg_end(), residualsPos_begin(), residualsPos_end(), value(), and eostools::which().

                                          {
    double sum = 0.;
    double n = 0.;
    for (std::vector<double *>::const_iterator r = residualsPos_begin(); r != residualsPos_end(); ++r) {
      double value = (*r)[which];
      if (fabs(value) < truncate) {
        sum += value;
        n += 1.;
      }
    }
    if (m_twoBin) {
      for (std::vector<double *>::const_iterator r = residualsNeg_begin(); r != residualsNeg_end(); ++r) {
        double value = (*r)[which];
        if (fabs(value) < truncate) {
          sum += value;
          n += 1.;
        }
      }
    }
    return sum / n;
  };

double MuonResidualsTwoBin::median ( int  which )

inline

Definition at line 119 of file MuonResidualsTwoBin.h.

References m_twoBin, alignCSCRings::r, residualsNeg_begin(), residualsNeg_end(), residualsPos_begin(), and residualsPos_end().

                           {
    std::vector<double> residuals;
    for (std::vector<double *>::const_iterator r = residualsPos_begin(); r != residualsPos_end(); ++r) {
      residuals.push_back((*r)[which]);
    }
    if (m_twoBin) {
      for (std::vector<double *>::const_iterator r = residualsNeg_begin(); r != residualsNeg_end(); ++r) {
        residuals.push_back((*r)[which]);
      }
    }
    std::sort(residuals.begin(), residuals.end());
    int length = residuals.size();
    return residuals[length / 2];
  };

int MuonResidualsTwoBin::ndata ( )

inline

Definition at line 36 of file MuonResidualsTwoBin.h.

References cms::cuda::assert(), m_neg, m_pos, and MuonResidualsFitter::ndata().

              {
    assert(m_pos->ndata() == m_neg->ndata());
    return m_pos->ndata();
  };

int MuonResidualsTwoBin::npar ( )

inline

Definition at line 32 of file MuonResidualsTwoBin.h.

References cms::cuda::assert(), m_neg, m_pos, and MuonResidualsFitter::npar().

             {
    assert(m_pos->npar() == m_neg->npar());
    return m_pos->npar();
  };

long MuonResidualsTwoBin::numResidualsNeg ( ) const

inline

Definition at line 31 of file MuonResidualsTwoBin.h.

References m_neg, and MuonResidualsFitter::numResiduals().

{ return m_neg->numResiduals(); };

long MuonResidualsTwoBin::numResidualsPos ( ) const

inline

Definition at line 30 of file MuonResidualsTwoBin.h.

References m_pos, and MuonResidualsFitter::numResiduals().

{ return m_pos->numResiduals(); };

double MuonResidualsTwoBin::numsegments ( )

inline

Definition at line 82 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, and MuonResidualsFitter::numsegments().

{ return (m_twoBin ? (m_pos->numsegments() + m_neg->numsegments()) : m_pos->numsegments()); };

double MuonResidualsTwoBin::plot ( std::string  name, TFileDirectory *  dir, Alignable *  ali  )

inline

Definition at line 88 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, convertSQLitetoXML_cfg::output, MuonResidualsFitter::plot(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                   {
    if (m_twoBin) {
      std::string namePos = name + std::string("Pos");
      std::string nameNeg = name + std::string("Neg");
      double output = 0.;
      output += m_pos->plot(namePos, dir, ali);
      output += m_neg->plot(nameNeg, dir, ali);
      return output;
    } else {
      return m_pos->plot(name, dir, ali);
    }
  };

void MuonResidualsTwoBin::plotsimple ( std::string  name, TFileDirectory *  dir, int  which, double  multiplier  )

inline

Definition at line 209 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, MuonResidualsFitter::plotsimple(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                     {
    if (m_twoBin) {
      std::string namePos = name + std::string("Pos");
      std::string nameNeg = name + std::string("Neg");
      m_pos->plotsimple(namePos, dir, which, multiplier);
      m_neg->plotsimple(nameNeg, dir, which, multiplier);
    } else {
      m_pos->plotsimple(name, dir, which, multiplier);
    }
  };

void MuonResidualsTwoBin::plotweighted ( std::string  name, TFileDirectory *  dir, int  which, int  whichredchi2, double  multiplier  )

inline

Definition at line 220 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, MuonResidualsFitter::plotweighted(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                         {
    if (m_twoBin) {
      std::string namePos = name + std::string("Pos");
      std::string nameNeg = name + std::string("Neg");
      m_pos->plotweighted(namePos, dir, which, whichredchi2, multiplier);
      m_neg->plotweighted(nameNeg, dir, which, whichredchi2, multiplier);
    } else {
      m_pos->plotweighted(name, dir, which, whichredchi2, multiplier);
    }
  };

void MuonResidualsTwoBin::read ( FILE *  file, int  which = 0  )

inline

Definition at line 110 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, MuonResidualsFitter::read(), and eostools::which().

Referenced by edmIntegrityCheck.PublishToFileSystem::get(), and MuonAlignmentFromReference::readTmpFiles().

                                       {
    if (m_twoBin) {
      m_pos->read(file, 2 * which);
      m_neg->read(file, 2 * which + 1);
    } else {
      m_pos->read(file, which);
    }
  };

int MuonResidualsTwoBin::residualsModel ( ) const

inline

Definition at line 26 of file MuonResidualsTwoBin.h.

References cms::cuda::assert(), m_neg, m_pos, and MuonResidualsFitter::residualsModel().

                             {
    assert(m_pos->residualsModel() == m_neg->residualsModel());
    return m_pos->residualsModel();
  };

std::vector<double *>::const_iterator MuonResidualsTwoBin::residualsNeg_begin ( ) const

inline

Definition at line 258 of file MuonResidualsTwoBin.h.

References m_neg, and MuonResidualsFitter::residuals_begin().

Referenced by mean(), median(), stdev(), and wmean().

{ return m_neg->residuals_begin(); };

std::vector<double *>::const_iterator MuonResidualsTwoBin::residualsNeg_end ( ) const

inline

Definition at line 259 of file MuonResidualsTwoBin.h.

References m_neg, and MuonResidualsFitter::residuals_end().

Referenced by mean(), median(), stdev(), and wmean().

{ return m_neg->residuals_end(); };

std::vector<bool>::const_iterator MuonResidualsTwoBin::residualsNeg_ok_begin ( ) const

inline

Definition at line 263 of file MuonResidualsTwoBin.h.

References m_neg, and MuonResidualsFitter::selectedResidualsFlags().

{ return m_neg->selectedResidualsFlags().begin(); };

std::vector<bool>::const_iterator MuonResidualsTwoBin::residualsNeg_ok_end ( ) const

inline

Definition at line 264 of file MuonResidualsTwoBin.h.

References m_neg, and MuonResidualsFitter::selectedResidualsFlags().

{ return m_neg->selectedResidualsFlags().end(); };

std::vector<double *>::const_iterator MuonResidualsTwoBin::residualsPos_begin ( ) const

inline

Definition at line 256 of file MuonResidualsTwoBin.h.

References m_pos, and MuonResidualsFitter::residuals_begin().

Referenced by MuonAlignmentFromReference::fillNtuple(), mean(), median(), stdev(), and wmean().

{ return m_pos->residuals_begin(); };

std::vector<double *>::const_iterator MuonResidualsTwoBin::residualsPos_end ( ) const

inline

Definition at line 257 of file MuonResidualsTwoBin.h.

References m_pos, and MuonResidualsFitter::residuals_end().

Referenced by MuonAlignmentFromReference::fillNtuple(), mean(), median(), stdev(), and wmean().

{ return m_pos->residuals_end(); };

std::vector<bool>::const_iterator MuonResidualsTwoBin::residualsPos_ok_begin ( ) const

inline

Definition at line 261 of file MuonResidualsTwoBin.h.

References m_pos, and MuonResidualsFitter::selectedResidualsFlags().

Referenced by MuonAlignmentFromReference::fillNtuple().

{ return m_pos->selectedResidualsFlags().begin(); };

std::vector<bool>::const_iterator MuonResidualsTwoBin::residualsPos_ok_end ( ) const

inline

Definition at line 262 of file MuonResidualsTwoBin.h.

References m_pos, and MuonResidualsFitter::selectedResidualsFlags().

{ return m_pos->selectedResidualsFlags().end(); };

void MuonResidualsTwoBin::selectPeakResiduals ( double  nsigma, int  nvar, int *  vars  )

inline

Definition at line 231 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, and MuonResidualsFitter::selectPeakResiduals().

Referenced by MuonAlignmentFromReference::selectResidualsPeaks().

                                                               {
    if (m_twoBin) {
      m_pos->selectPeakResiduals(nsigma, nvar, vars);
      m_neg->selectPeakResiduals(nsigma, nvar, vars);
    } else {
      m_pos->selectPeakResiduals(nsigma, nvar, vars);
    }
  }

void MuonResidualsTwoBin::setPrintLevel ( int  printLevel ) const

inline

Definition at line 53 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, and MuonResidualsFitter::setPrintLevel().

                                           {
    m_pos->setPrintLevel(printLevel);
    m_neg->setPrintLevel(printLevel);
  }

void MuonResidualsTwoBin::setStrategy ( int  strategy ) const

inline

Definition at line 58 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, and MuonResidualsFitter::setStrategy().

                                       {
    m_pos->setStrategy(strategy);
    m_neg->setStrategy(strategy);
  }

double MuonResidualsTwoBin::stdev ( int  which, double  truncate  )

inline

Definition at line 184 of file MuonResidualsTwoBin.h.

References m_twoBin, dqmiodumpmetadata::n, funct::pow(), alignCSCRings::r, residualsNeg_begin(), residualsNeg_end(), residualsPos_begin(), residualsPos_end(), mathSSE::sqrt(), value(), and eostools::which().

                                           {
    double sum2 = 0.;
    double sum = 0.;
    double n = 0.;
    for (std::vector<double *>::const_iterator r = residualsPos_begin(); r != residualsPos_end(); ++r) {
      double value = (*r)[which];
      if (fabs(value) < truncate) {
        sum2 += value * value;
        sum += value;
        n += 1.;
      }
    }
    if (m_twoBin) {
      for (std::vector<double *>::const_iterator r = residualsNeg_begin(); r != residualsNeg_end(); ++r) {
        double value = (*r)[which];
        if (fabs(value) < truncate) {
          sum2 += value * value;
          sum += value;
          n += 1.;
        }
      }
    }
    return sqrt(sum2 / n - pow(sum / n, 2));
  };

double MuonResidualsTwoBin::sumofweights ( )

inline

Definition at line 83 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, and MuonResidualsFitter::sumofweights().

                        {
    return (m_twoBin ? (m_pos->sumofweights() + m_neg->sumofweights()) : m_pos->sumofweights());
  };

int MuonResidualsTwoBin::type ( ) const

inline

Definition at line 40 of file MuonResidualsTwoBin.h.

References cms::cuda::assert(), m_neg, m_pos, and MuonResidualsFitter::type().

Referenced by MuonAlignmentFromReference::fillNtuple(), and MuonAlignmentFromReference::selectResidualsPeaks().

                   {
    assert(m_pos->type() == m_neg->type());
    return m_pos->type();
  };

int MuonResidualsTwoBin::useRes ( ) const

inline

Definition at line 44 of file MuonResidualsTwoBin.h.

References m_pos, and MuonResidualsFitter::useRes().

Referenced by MuonAlignmentFromReference::selectResidualsPeaks().

{ return m_pos->useRes(); };

double MuonResidualsTwoBin::value ( int  parNum )

inline

Definition at line 71 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, and MuonResidualsFitter::value().

Referenced by Types.int32::__nonzero__(), Types.uint32::__nonzero__(), Types.int64::__nonzero__(), Types.uint64::__nonzero__(), Types.double::__nonzero__(), Types.bool::__nonzero__(), Types.string::__nonzero__(), average.Average::average(), Types.string::configValue(), Types.FileInPath::configValue(), Mixins.UsingBlock::dumpPython(), fix(), Mixins.UsingBlock::insertInto(), Types.int32::insertInto(), Types.uint32::insertInto(), Types.int64::insertInto(), Types.uint64::insertInto(), Types.double::insertInto(), Types.bool::insertInto(), Types.string::insertInto(), Types.FileInPath::insertInto(), Types.vint32::insertInto(), Types.vuint32::insertInto(), Types.vint64::insertInto(), Types.vuint64::insertInto(), Types.vdouble::insertInto(), Types.vbool::insertInto(), Types.vstring::insertInto(), mean(), stdev(), and wmean().

                           {
    return (m_twoBin ? ((m_pos->value(parNum) + m_neg->value(parNum)) / 2.) : m_pos->value(parNum));
  };

double MuonResidualsTwoBin::wmean ( int  which, int  whichredchi2, double  truncate  )

inline

Definition at line 156 of file MuonResidualsTwoBin.h.

References m_twoBin, dqmiodumpmetadata::n, alignCSCRings::r, residualsNeg_begin(), residualsNeg_end(), residualsPos_begin(), residualsPos_end(), value(), histoStyle::weight, and eostools::which().

                                                             {
    double sum = 0.;
    double n = 0.;
    for (std::vector<double *>::const_iterator r = residualsPos_begin(); r != residualsPos_end(); ++r) {
      double value = (*r)[which];
      if (fabs(value) < truncate) {
        double weight = 1. / (*r)[whichredchi2];
        if (TMath::Prob(1. / weight * 12, 12) < 0.99) {
          sum += weight * value;
          n += weight;
        }
      }
    }
    if (m_twoBin) {
      for (std::vector<double *>::const_iterator r = residualsNeg_begin(); r != residualsNeg_end(); ++r) {
        double value = (*r)[which];
        if (fabs(value) < truncate) {
          double weight = 1. / (*r)[whichredchi2];
          if (TMath::Prob(1. / weight * 12, 12) < 0.99) {
            sum += weight * value;
            n += weight;
          }
        }
      }
    }
    return sum / n;
  };

void MuonResidualsTwoBin::write ( FILE *  file, int  which = 0  )

inline

Definition at line 102 of file MuonResidualsTwoBin.h.

References m_neg, m_pos, m_twoBin, eostools::which(), and MuonResidualsFitter::write().

Referenced by pkg.AbstractPkg::generate(), querying.connection::write_and_commit(), and MuonAlignmentFromReference::writeTmpFiles().

                                        {
    if (m_twoBin) {
      m_pos->write(file, 2 * which);
      m_neg->write(file, 2 * which + 1);
    } else {
      m_pos->write(file, which);
    }
  };

Member Data Documentation

MuonResidualsFitter * MuonResidualsTwoBin::m_neg

protected

Definition at line 268 of file MuonResidualsTwoBin.h.

Referenced by antisym(), eraseNotSelectedResiduals(), errorerror(), fill(), fit(), fix(), fixed(), loglikelihood(), ndata(), npar(), numResidualsNeg(), numsegments(), plot(), plotsimple(), plotweighted(), read(), residualsModel(), residualsNeg_begin(), residualsNeg_end(), residualsNeg_ok_begin(), residualsNeg_ok_end(), selectPeakResiduals(), setPrintLevel(), setStrategy(), sumofweights(), type(), value(), write(), and ~MuonResidualsTwoBin().

MuonResidualsFitter* MuonResidualsTwoBin::m_pos

protected

Definition at line 268 of file MuonResidualsTwoBin.h.

Referenced by antisym(), correctBField(), eraseNotSelectedResiduals(), errorerror(), fiducialCuts(), fill(), fit(), fix(), fixed(), loglikelihood(), ndata(), npar(), numResidualsPos(), numsegments(), plot(), plotsimple(), plotweighted(), read(), residualsModel(), residualsPos_begin(), residualsPos_end(), residualsPos_ok_begin(), residualsPos_ok_end(), selectPeakResiduals(), setPrintLevel(), setStrategy(), sumofweights(), type(), useRes(), value(), write(), and ~MuonResidualsTwoBin().

bool MuonResidualsTwoBin::m_twoBin

protected

Definition at line 264 of file MuonResidualsTwoBin.h.

Referenced by antisym(), eraseNotSelectedResiduals(), errorerror(), fill(), fit(), loglikelihood(), mean(), median(), numsegments(), plot(), plotsimple(), plotweighted(), read(), selectPeakResiduals(), stdev(), sumofweights(), value(), wmean(), and write().