#include <PndSdsStripCorrelator.h>

Public Member Functions
	PndSdsStripCorrelator (int mode=0, double cut=0., double noise=0., double threshold=0.)

	~PndSdsStripCorrelator ()

void	Reset ()

void	Setup (int mode=0, double cut=0., double noise=0., double threshold=0.)

void	AddCluster (int moduleId, int side, int clusterIndex, double charge)

vector< pair< int, int > >	GetCorrelationList ()

vector< double >	GetProbList ()

vector< double >	GetSecondProbList ()

vector< int >	GetMultProbList ()

Private Member Functions
void	CalcChargeDifferenceCut ()

void	CalcLikelihoodAlgo ()

void	CalcAll ()

std::vector < PndSdsStripCorrelatorCombi >	getCombinations (std::map< int, std::map< int, PndSdsStripCorrelatorCand > > matrix, int cols, int rows)

std::map< int, std::map< int, PndSdsStripCorrelatorCand > >	getSubMatrix (std::map< int, std::map< int, PndSdsStripCorrelatorCand > > matrix, int cols, int rows, int pivotCol, int pivotRow)

Private Attributes
map< int, map< int, vector < pair< int, double > > > >	fClusterList

vector< pair< int, int > >	fCorrelationList

vector< double >	fCorrelationProbList

vector< double >	fSecondProbList

vector< int >	fMultProbList

double	fCut

double	fNoise

double	fThreshold

int	fMode

bool	fCalculated

Detailed Description

Definition at line 59 of file PndSdsStripCorrelator.h.

Constructor & Destructor Documentation

PndSdsStripCorrelator::PndSdsStripCorrelator	(	int	mode = `0`,
		double	cut = `0.`,
		double	noise = `0.`,
		double	threshold = `0.`
	)

Definition at line 19 of file PndSdsStripCorrelator.cxx.

                                                                                                  :
   fClusterList(),
   fCorrelationList(),
   fCorrelationProbList(),
   fSecondProbList(),
   fMultProbList(),
   fCut(cut),
   fNoise(noise),
   fThreshold(threshold),
   fMode(mode),
   fCalculated(false)
 {}

PndSdsStripCorrelator::~PndSdsStripCorrelator ( )

Definition at line 33 of file PndSdsStripCorrelator.cxx.

33 {}

Member Function Documentation

void PndSdsStripCorrelator::AddCluster	(	int	moduleId,
		int	side,
		int	clusterIndex,
		double	charge
	)

Definition at line 36 of file PndSdsStripCorrelator.cxx.

References fClusterList.

 {
   fClusterList[moduleId][side].push_back(make_pair(clusterIndex,charge));
 }

void PndSdsStripCorrelator::CalcAll ( )

private

Definition at line 117 of file PndSdsStripCorrelator.cxx.

References CalcChargeDifferenceCut(), and fCut.

Referenced by GetCorrelationList().

 {
   double oldcut = fCut;
   fCut=0.;
   CalcChargeDifferenceCut();
   fCut=oldcut;
 }

void PndSdsStripCorrelator::CalcChargeDifferenceCut ( )

private

Definition at line 126 of file PndSdsStripCorrelator.cxx.

References fabs(), fClusterList, fCorrelationList, fCorrelationProbList, fCut, fSecondProbList, and map.

Referenced by CalcAll(), and GetCorrelationList().

 {
   double dq=0;
   fCorrelationList.clear();
   for (map<int,map<int,vector<pair<int,double> > > >::iterator modIt=fClusterList.begin(); modIt!=fClusterList.end(); ++modIt) 
   {
     for (vector<pair<int,double> >::iterator topIt = modIt->second[0].begin(); topIt != modIt->second[0].end(); ++topIt) 
     {
       for (vector<pair<int,double> >::iterator botIt = modIt->second[1].begin(); botIt != modIt->second[1].end(); ++botIt) 
       {
         dq=fabs((*topIt).second - (*botIt).second);
         if (dq<fCut) 
         {
           fCorrelationList.push_back(make_pair((*topIt).first,(*botIt).first));
           fCorrelationProbList.push_back( (dq>5e4) ? 0 : (1 - 2e-5*dq) );
           fSecondProbList.push_back( -1 );
         }
       }
     }
   }
 }

void PndSdsStripCorrelator::CalcLikelihoodAlgo ( )

private

Definition at line 149 of file PndSdsStripCorrelator.cxx.

References CLUSTER_MULT_CUTOFF, fabs(), fClusterList, fCorrelationList, fCorrelationProbList, fMultProbList, fSecondProbList, fThreshold, getCombinations(), map, sqrt(), and top.

Referenced by GetCorrelationList().

 {
   //per sensor module fill a combinations matrix, calculate likelihoods, assign b
   
   fCorrelationList.clear();
   for (map<int,map<int,vector<pair<int,double> > > >::iterator modIt=fClusterList.begin(); modIt!=fClusterList.end(); ++modIt) 
   {
     if (modIt->second[0].size()>CLUSTER_MULT_CUTOFF) continue; // skip modules with too busy events
     if (modIt->second[1].size()>CLUSTER_MULT_CUTOFF) continue; // skip modules with too busy events
                                                                //fill matrix
     std::map<int,std::map<int,PndSdsStripCorrelatorCand> > matrix;
     double top_total=0.,wert=0.;
     int nTop=0, nBot=0, iBot=0; //,iTop=-1,iBot=-1;
     for (vector<pair<int,double> >::iterator topIt = modIt->second[0].begin(); topIt != modIt->second[0].end(); ++topIt) 
     {
       for (vector<pair<int,double> >::iterator botIt = modIt->second[1].begin(); botIt != modIt->second[1].end(); ++botIt) 
       {
         wert = (*topIt).second - (*botIt).second;
         //        wert *= wert;
         //        top_total += wert;
         top_total += wert*wert;
         wert = erfc(fabs(wert)/(2.*fThreshold)/sqrt(2.)); // [R.K.31.7.'12]
         //wert = erfc(fabs(wert)/(6.*fNoise)/sqrt(2.)); // original
         matrix[nTop][iBot]=PndSdsStripCorrelatorCand((*topIt).first,(*botIt).first,(*topIt).second,(*botIt).second,wert);
         //printf("add to matrix: t=%i, b=%i, qt=%f, qb=%f, p=%f \n",(*topIt).first,(*botIt).first,(*topIt).second,(*botIt).second,wert);
         iBot++;
       }
       if(iBot>nBot)nBot=iBot;
       iBot=0;
       nTop++;
     }
     
     // Calculate combinations 
     std::vector<PndSdsStripCorrelatorCombi> combinations = getCombinations(matrix,nTop,nBot);
     
     if(combinations.size()==0) continue; //next module
     
     //select the lowest probability
     double klein=-1.,ganzklein=-1.;
     int who=-1, whoelse=-1;
     
     for(int like=0;like<(int)combinations.size();like++)
     {
       if(combinations[like].prob>klein)
       { 
         if(combinations[like].prob>ganzklein)
         {
           klein=ganzklein;
           whoelse=who;
           ganzklein=combinations[like].prob;
           who=like;
         }
         else
         {
           klein=combinations[like].prob;
           whoelse=like;
         }
       }
     }
     //    if(combinations[who].prob>fCut)
     //    { // index pairs for the output
     for (int kk=0; kk<(int)combinations[who].pairlist.size(); ++kk) 
     {
       fCorrelationList.push_back(make_pair(combinations[who].pairlist[kk].top,combinations[who].pairlist[kk].bot));
       //printf("add correlation [%i|%i] p=%f \n",combinations[who].pairlist[kk].top,combinations[who].pairlist[kk].bot,combinations[who].prob);
       fCorrelationProbList.push_back(combinations[who].prob);
       if(whoelse>=0) fSecondProbList.push_back(combinations[whoelse].prob);
       else fSecondProbList.push_back(-1.);
       fMultProbList.push_back(combinations.size());
     }
     //    }
   }
 }

std::vector< PndSdsStripCorrelatorCombi > PndSdsStripCorrelator::getCombinations	(	std::map< int, std::map< int, PndSdsStripCorrelatorCand > >	matrix,
		int	cols,
		int	rows
	)

private

Definition at line 226 of file PndSdsStripCorrelator.cxx.

References cols, combi, fabs(), fThreshold, getSubMatrix(), i, PndSdsStripCorrelatorCombi::pairlist, PndSdsStripCorrelatorCombi::prob, sqrt(), and top.

Referenced by CalcLikelihoodAlgo().

 {
   //printf("getCombinations: Matrix(%i,%i) \n",cols,rows);
   std::vector<PndSdsStripCorrelatorCombi> combinations;
   if (cols==1 && rows==1) {
     PndSdsStripCorrelatorCombi combi;
     combi.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[0][0]));
     combi.prob=matrix[0][0].prob;
     combinations.push_back(combi);
     return combinations;
   }
   if (cols==1) {
     PndSdsStripCorrelatorCombi combi;
     double q_t[rows];
     double bsum=0.;
     for (int i=0; i< rows; i++) bsum+=matrix[0][i].q_bot;
     for (int i=0; i< rows; i++) q_t[i]=matrix[0][0].q_top - bsum + matrix[0][i].q_bot;
     //double tot=0.; //unused?
     //double diag=0.; //unused?
     for (int i=0; i< rows; i++)
     {
       combi.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[0][i].top,matrix[0][i].bot,q_t[i],matrix[0][i].q_bot,matrix[0][i].prob));
       combi.prob*=erfc(fabs(q_t[i]-matrix[0][i].q_bot)/(2.*fThreshold*sqrt(2.)));// [R.K.31.7.'12]
       //combi.prob*=erfc(fabs(q_t[i]-matrix[0][i].q_bot)/(6.*fNoise*sqrt(2.)));//original
     }
     double prob_full=combi.prob;
     combinations.push_back(combi);
     //double _maxprob = combi.prob; //unused?
     //int _maxcombi = 0; //unused?
     for (int i=0; i<rows; ++i)
     { // Recursive call
       std::vector<PndSdsStripCorrelatorCombi> combilist_reduced=getCombinations(getSubMatrix(matrix,1,rows,-1,i),1,rows-1);
       for (int ii=0; ii<(int)combilist_reduced.size(); ii++) 
       { 
         combilist_reduced[ii].prob*=(1.-prob_full);
       }
       combinations.insert(combinations.end(),combilist_reduced.begin(),combilist_reduced.end());
     }
     return combinations;
   }
   if (rows==1) {
     PndSdsStripCorrelatorCombi combi;
     double q_b[cols];
     double tsum=0.;
     for (int i=0; i< cols; i++) tsum+=matrix[i][0].q_top;
     for (int i=0; i< cols; i++) q_b[i]=matrix[0][0].q_bot - tsum + matrix[i][0].q_top;
     //double tot=0.; //unused?
     //double diag=0.; //unused?
     for (int i=0; i< cols; i++)
     {
       combi.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[i][0].top,matrix[i][0].bot,matrix[i][0].q_top,q_b[i],matrix[i][0].prob));
       combi.prob*=erfc(fabs(matrix[i][0].q_top-q_b[i])/(2.*fThreshold*sqrt(2.)));// [R.K.31.7.'12]
       //combi.prob*=erfc(fabs(matrix[i][0].q_top-q_b[i])/(6.*fNoise*sqrt(2.)));//original
     }
     double prob_full=combi.prob;
     combinations.push_back(combi);
     //double _maxprob = combi.prob; //unused?
     //int _maxcombi = 0; //unused?
     for (int i=0; i<cols; ++i)
     { // Recursive call
       std::vector<PndSdsStripCorrelatorCombi> combilist_reduced=getCombinations(getSubMatrix(matrix,cols,1,i,-1),cols-1,1);
       for (int ii=0; ii<(int)combilist_reduced.size(); ii++)
       {
         combilist_reduced[ii].prob*=(1.-prob_full);
       }
       combinations.insert(combinations.end(),combilist_reduced.begin(),combilist_reduced.end());
     }
     return combinations;
   }
   // treat 2x2 matrix as separate case since the nr. of independent combinations is just half of the nr. of possible combinations
   if (rows==2 && cols==2) {
     PndSdsStripCorrelatorCombi combi1;
     PndSdsStripCorrelatorCombi combi2;
     combi1.prob=1.;
     combi2.prob=1.;
     combi1.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[0][0]));
     combi1.prob*=matrix[0][0].prob;
     combi1.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[1][1]));
     combi1.prob*=matrix[1][1].prob;
     combi2.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[1][0]));
     combi2.prob*=matrix[1][0].prob;
     combi2.pairlist.push_back(PndSdsStripCorrelatorCand(matrix[0][1]));
     combi2.prob*=matrix[0][1].prob;
     combinations.push_back(combi1);
     combinations.push_back(combi2);
     return combinations;
   }
   for (int i=0; i<cols; i++) {
     for (int j=0; j<rows; j++)
     { // Recursive call
       std::vector<PndSdsStripCorrelatorCombi> combilist=getCombinations(getSubMatrix(matrix,cols,rows,i,j),cols-1,rows-1);
       for (int k=0; k<(int)combilist.size(); k++) {
         combilist[k].pairlist.push_back(PndSdsStripCorrelatorCand(matrix[i][j]));
         combilist[k].prob*=matrix[i][j].prob;
         combinations.push_back(PndSdsStripCorrelatorCombi(combilist[k].pairlist,combilist[k].prob));
       }
     }
   }
   return combinations;
 }

vector< pair< int, int > > PndSdsStripCorrelator::GetCorrelationList ( )

Definition at line 63 of file PndSdsStripCorrelator.cxx.

References CalcAll(), CalcChargeDifferenceCut(), CalcLikelihoodAlgo(), fCalculated, fCorrelationList, and fMode.

 {
   if(!fCalculated)
   {
     switch (fMode) {
       case 0:
         CalcChargeDifferenceCut();
         break;
       case 1:
         CalcLikelihoodAlgo();
         break;
       default:
         CalcAll();
         break;
     }
     fCalculated=true;
   }
   return fCorrelationList;
 }

vector< int > PndSdsStripCorrelator::GetMultProbList ( )

Definition at line 106 of file PndSdsStripCorrelator.cxx.

References fCalculated, fCorrelationProbList, fMultProbList, and fSecondProbList.

 {
   if(!fCalculated){
     fCorrelationProbList.clear(); // empty on invalid state, just in case
     fSecondProbList.clear(); // empty on invalid state, just in case
     fMultProbList.clear(); // empty on invalid state, just in case
   }
   return fMultProbList;
 }

vector< double > PndSdsStripCorrelator::GetProbList ( )

Definition at line 84 of file PndSdsStripCorrelator.cxx.

References fCalculated, fCorrelationProbList, fMultProbList, and fSecondProbList.

 {
   if(!fCalculated){
     fCorrelationProbList.clear(); // empty on invalid state, just in case
     fSecondProbList.clear(); // empty on invalid state, just in case
     fMultProbList.clear(); // empty on invalid state, just in case
   }
   return fCorrelationProbList;
 }

vector< double > PndSdsStripCorrelator::GetSecondProbList ( )

Definition at line 95 of file PndSdsStripCorrelator.cxx.

References fCalculated, fCorrelationProbList, fMultProbList, and fSecondProbList.

 {
   if(!fCalculated){
     fCorrelationProbList.clear(); // empty on invalid state, just in case
     fSecondProbList.clear(); // empty on invalid state, just in case
     fMultProbList.clear(); // empty on invalid state, just in case
   }
   return fSecondProbList;
 }

std::map< int, std::map< int, PndSdsStripCorrelatorCand > > PndSdsStripCorrelator::getSubMatrix	(	std::map< int, std::map< int, PndSdsStripCorrelatorCand > >	matrix,
		int	cols,
		int	rows,
		int	pivotCol,
		int	pivotRow
	)

private

Definition at line 329 of file PndSdsStripCorrelator.cxx.

References cols, i, and jj.

Referenced by getCombinations().

 {
   std::map<int,std::map<int,PndSdsStripCorrelatorCand> > result;
   int ii=0;
   for (int i=0; i<cols;i++)
   {
     if (i==pivotCol) continue;
     int jj=0;
     for (int j=0; j<rows; j++)
     {
       result[ii][jj]=matrix[i][j];
       if (j==pivotRow) continue;
       jj++;
     }
     ii++;
   }
   return result;
 }

void PndSdsStripCorrelator::Reset ( )

Definition at line 42 of file PndSdsStripCorrelator.cxx.

References fCalculated, fClusterList, fCorrelationList, fCorrelationProbList, fMultProbList, and fSecondProbList.

 {
   fClusterList.clear();
   fCorrelationList.clear();
   fCorrelationProbList.clear();
   fSecondProbList.clear();
   fMultProbList.clear();
   fCalculated=false;
 }

void PndSdsStripCorrelator::Setup	(	int	mode = `0`,
		double	cut = `0.`,
		double	noise = `0.`,
		double	threshold = `0.`
	)

Definition at line 53 of file PndSdsStripCorrelator.cxx.

References cut, fCalculated, fCut, fMode, fNoise, fThreshold, mode, noise, and threshold.

 {
   fMode=mode;
   fCut=cut;
   fNoise=noise;
   fThreshold=threshold;
   fCalculated=false;
 }

Member Data Documentation

bool PndSdsStripCorrelator::fCalculated

private

Definition at line 95 of file PndSdsStripCorrelator.h.

Referenced by GetCorrelationList(), GetMultProbList(), GetProbList(), GetSecondProbList(), Reset(), and Setup().

map<int,map<int,vector<pair<int,double> > > > PndSdsStripCorrelator::fClusterList

private

Definition at line 85 of file PndSdsStripCorrelator.h.

Referenced by AddCluster(), CalcChargeDifferenceCut(), CalcLikelihoodAlgo(), and Reset().

vector<pair<int,int> > PndSdsStripCorrelator::fCorrelationList

private

Definition at line 86 of file PndSdsStripCorrelator.h.

Referenced by CalcChargeDifferenceCut(), CalcLikelihoodAlgo(), GetCorrelationList(), and Reset().

vector<double> PndSdsStripCorrelator::fCorrelationProbList

private

Definition at line 87 of file PndSdsStripCorrelator.h.

Referenced by CalcChargeDifferenceCut(), CalcLikelihoodAlgo(), GetMultProbList(), GetProbList(), GetSecondProbList(), and Reset().

double PndSdsStripCorrelator::fCut

private

Definition at line 91 of file PndSdsStripCorrelator.h.

Referenced by CalcAll(), CalcChargeDifferenceCut(), and Setup().

int PndSdsStripCorrelator::fMode

private

Definition at line 94 of file PndSdsStripCorrelator.h.

Referenced by GetCorrelationList(), and Setup().

vector<int> PndSdsStripCorrelator::fMultProbList

private

Definition at line 89 of file PndSdsStripCorrelator.h.

Referenced by CalcLikelihoodAlgo(), GetMultProbList(), GetProbList(), GetSecondProbList(), and Reset().

double PndSdsStripCorrelator::fNoise

private

Definition at line 92 of file PndSdsStripCorrelator.h.

Referenced by Setup().

vector<double> PndSdsStripCorrelator::fSecondProbList

private

Definition at line 88 of file PndSdsStripCorrelator.h.

Referenced by CalcChargeDifferenceCut(), CalcLikelihoodAlgo(), GetMultProbList(), GetProbList(), GetSecondProbList(), and Reset().

double PndSdsStripCorrelator::fThreshold

private

Definition at line 93 of file PndSdsStripCorrelator.h.

Referenced by CalcLikelihoodAlgo(), getCombinations(), and Setup().

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation