#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.

       : 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.

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

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

Referenced by operator()().

◆ operator()()

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

Definition at line 52 of file CrosstalkInversion.cc.

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

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

◆ unfold()

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

static

Definition at line 7 of file CrosstalkInversion.cc.

                                                                                                      {
     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;
   }