reco::InverseCrosstalkMatrix Class Reference

#include <CrosstalkInversion.h>

Public Member Functions
	InverseCrosstalkMatrix (const unsigned N, const float x)
float	operator() (const unsigned i, const unsigned j) const

Static Public Member Functions
static std::vector< stats_t< float > >	unfold (const SiStripCluster &q, const float x)

Private Member Functions
float	element (const unsigned, const unsigned) const

Private Attributes
const double	denominator
const double	lambdaM
const double	lambdaP
const unsigned	N
const double	sq

Detailed Description

Definition at line 13 of file CrosstalkInversion.h.

Constructor & Destructor Documentation

InverseCrosstalkMatrix()

reco::InverseCrosstalkMatrix::InverseCrosstalkMatrix	(	const unsigned	N,
		const float	x
	)

Definition at line 45 of file CrosstalkInversion.cc.

References x.

      : N(x > 0 ? N : 0),
        sq(sqrt(-x * 4 + 1)),
        lambdaP(1 + (1 + sq) / (-x * 2)),
        lambdaM(1 + (1 - sq) / (-x * 2)),
        denominator(sq * (pow(lambdaP, N + 1) - pow(lambdaM, N + 1))) {}

Member Function Documentation

element()

float reco::InverseCrosstalkMatrix::element ( const unsigned  i, const unsigned  j  ) const

inlineprivate

Definition at line 56 of file CrosstalkInversion.cc.

References denominator, mps_fire::i, dqmiolumiharvest::j, lambdaM, lambdaP, N, and conifer::pow().

Referenced by operator()().

                                                                                       {
    return (pow(lambdaM, N + 1 - i) - pow(lambdaP, N + 1 - i)) * (pow(lambdaM, j) - pow(lambdaP, j)) / denominator;
  }

operator()()

float reco::InverseCrosstalkMatrix::operator()	(	const unsigned	i,
		const unsigned	j
	)		const

Definition at line 52 of file CrosstalkInversion.cc.

References denominator, element(), mps_fire::i, edm::isNotFinite(), dqmiolumiharvest::j, and N.

                                                                                   {
    return N == 0 || edm::isNotFinite(denominator) ? i == j : i >= j ? element(i, j) : element(j, i);
  }

unfold()

std::vector< stats_t< float > > reco::InverseCrosstalkMatrix::unfold ( const SiStripCluster &  q, const float  x  )

static

Definition at line 7 of file CrosstalkInversion.cc.

References A, SiStripCluster::amplitudes(), mps_fire::i, dqmiolumiharvest::j, N, submitPVResolutionJobs::q, SiStripMonitorCluster_cfi::q0, STATS, and x.

Referenced by StripCPEgeometric::localParameters().

                                                                                                     {
    auto const& q = clus.amplitudes();
    const stats_t<float> suppressed(-5, 100);
    const stats_t<float> saturated(254, 400);
#define STATS(value) ((value < 254) ? stats_t<float>(value) : saturated)

    const unsigned N = q.size();
    std::vector<stats_t<float> > Q(N + 2, stats_t<float>(0));
    Q[0] = Q[N + 1] = suppressed;

    if (N == 1)  //optimize N==1
      Q[1] = STATS(q[0]) / (1 - 2 * x);
    else if (N == 2) {  //optimize N==2
      const double A = 1 - 2 * x;
      const stats_t<float> q0 = STATS(q[0]);
      const stats_t<float> q1 = STATS(q[1]);
      Q[1] = (A * q0 - x * q1) / (A * A - x * x);
      Q[2] = (A * q1 - x * q0) / (A * A - x * x);
    } else {  //general case
      const InverseCrosstalkMatrix inverse(N, x);
      for (unsigned i = 0; i < (N + 1) / 2; i++) {
        for (unsigned j = i; j < N - i; j++) {
          const float Cij = inverse(i + 1, j + 1);
          Q[i + 1] += Cij * STATS(q[j]);
          if (i != j)
            Q[j + 1] += Cij * STATS(q[i]);
          if (N != i + j + 1) {
            Q[N - i] += Cij * STATS(q[N - j - 1]);
            if (i != j)
              Q[N - j] += Cij * STATS(q[N - i - 1]);
          }
        }
      }
    }
#undef STATS
    return Q;
  }

Member Data Documentation

denominator

const double reco::InverseCrosstalkMatrix::denominator

private

Definition at line 21 of file CrosstalkInversion.h.

Referenced by element(), and operator()().

lambdaM

const double reco::InverseCrosstalkMatrix::lambdaM

private

Definition at line 21 of file CrosstalkInversion.h.

Referenced by element().

lambdaP

const double reco::InverseCrosstalkMatrix::lambdaP

private

Definition at line 21 of file CrosstalkInversion.h.

Referenced by element().

N

const unsigned reco::InverseCrosstalkMatrix::N

private

Definition at line 20 of file CrosstalkInversion.h.

Referenced by element(), operator()(), and unfold().

sq

const double reco::InverseCrosstalkMatrix::sq

private

Definition at line 21 of file CrosstalkInversion.h.