#include <GoldenPattern.h>

Public Types
typedef std::pair< int, bool >	layerResult

typedef std::vector< int >	vector1D

typedef std::vector< vector1D >	vector2D

typedef std::vector< vector2D >	vector3D

Public Member Functions
void	addCount (unsigned int iRefLayer, unsigned int iLayer, const int refPhi, const OMTFinput::vector1D &layerHits)
	Add a single count to the relevant pdf bin in three dimensions. More...

const vector2D &	getMeanDistPhi () const

const vector3D &	getPdf () const

	GoldenPattern (const Key &aKey, const OMTFConfiguration *omtfConfig)

bool	hasCounts ()
	Check if the GP has any counts in any of referecne layers;. More...

Key	key () const

int	meanDistPhiValue (unsigned int iLayer, unsigned int iRefLayer) const

void	normalise (unsigned int nPdfAddrBits)

int	pdfValue (unsigned int iLayer, unsigned int iRefLayer, unsigned int iBin) const

GoldenPattern::layerResult	process1Layer1RefLayer (unsigned int iRefLayer, unsigned int iLayer, const int refPhi, const OMTFinput::vector1D &layerHits)

int	propagateRefPhi (int phiRef, int etaRef, unsigned int iRefLayer)

void	reset ()
	Reset contents of all data vectors, keeping the vectors size. More...

void	setMeanDistPhi (const vector2D &aMeanDistPhi)

void	setPdf (const vector3D &aPdf)

Private Attributes
vector2D	meanDistPhi

vector2D	meanDistPhiCounts

const OMTFConfiguration *	myOmtfConfig

vector3D	pdfAllRef

Key	theKey
	Pattern kinematical identification (iEta,iPt,iCharge) More...

Friends
std::ostream &	operator<< (std::ostream &out, const GoldenPattern &aPattern)

Detailed Description

Definition at line 44 of file GoldenPattern.h.

Member Typedef Documentation

◆ layerResult

typedef std::pair<int, bool> GoldenPattern::layerResult

Definition at line 49 of file GoldenPattern.h.

◆ vector1D

typedef std::vector<int> GoldenPattern::vector1D

Definition at line 46 of file GoldenPattern.h.

◆ vector2D

typedef std::vector<vector1D> GoldenPattern::vector2D

Definition at line 47 of file GoldenPattern.h.

◆ vector3D

typedef std::vector<vector2D> GoldenPattern::vector3D

Definition at line 48 of file GoldenPattern.h.

Constructor & Destructor Documentation

◆ GoldenPattern()

GoldenPattern::GoldenPattern	(	const Key &	aKey,
		const OMTFConfiguration *	omtfConfig
	)

inline

Definition at line 54 of file GoldenPattern.h.

54 : theKey(aKey), myOmtfConfig(omtfConfig) {}

Member Function Documentation

◆ addCount()

void GoldenPattern::addCount	(	unsigned int	iRefLayer,
		unsigned int	iLayer,
		const int	refPhi,
		const OMTFinput::vector1D &	layerHits
	)

Add a single count to the relevant pdf bin in three dimensions.

For making the patterns take events with a single hit in each layer

Shift phiDist so it is in +-Pi range

Shift phidist, so 0 is at the middle of the range

Check if phiDist is within pdf range in -64 +63 U2 code Find more elegant way to check this.

Definition at line 50 of file GoldenPattern.cc.

                                                                  {
   int nHitsInLayer = 0;
   int phiDist = exp2(myOmtfConfig->nPdfAddrBits());
   for (auto itHit : layerHits) {
     if (itHit >= (int)myOmtfConfig->nPhiBins())
       continue;
     if (abs(itHit - phiRefHit) < phiDist)
       phiDist = itHit - phiRefHit;
     ++nHitsInLayer;
   }
   if (nHitsInLayer > 1 || nHitsInLayer == 0)
     return;
  
   if (phiDist >= (int)myOmtfConfig->nPhiBins() / 2)
     phiDist -= (int)myOmtfConfig->nPhiBins();
   if (phiDist <= -(int)myOmtfConfig->nPhiBins() / 2)
     phiDist += (int)myOmtfConfig->nPhiBins();
  
   int phiDistShift = phiDist + exp2(myOmtfConfig->nPdfAddrBits() - 1);
  
   if (phiDistShift < 0 || phiDistShift > exp2(myOmtfConfig->nPdfAddrBits()) - 1) {
     return;
   }
  
   if ((int)iLayer == myOmtfConfig->getRefToLogicNumber()[iRefLayer])
     ++meanDistPhiCounts[iLayer][iRefLayer];
   ++pdfAllRef[iLayer][iRefLayer][phiDistShift];
 }

References funct::abs(), OMTFConfiguration::getRefToLogicNumber(), meanDistPhiCounts, myOmtfConfig, OMTFConfiguration::nPdfAddrBits(), OMTFConfiguration::nPhiBins(), and pdfAllRef.

◆ getMeanDistPhi()

const vector2D& GoldenPattern::getMeanDistPhi ( ) const

inline

Definition at line 60 of file GoldenPattern.h.

60 { return meanDistPhi; }

References meanDistPhi.

Referenced by OMTFProcessor::averagePatterns().

◆ getPdf()

const vector3D& GoldenPattern::getPdf ( ) const

inline

Definition at line 62 of file GoldenPattern.h.

62 { return pdfAllRef; }

References pdfAllRef.

Referenced by OMTFProcessor::shiftGP().

◆ hasCounts()

bool GoldenPattern::hasCounts ( )

Check if the GP has any counts in any of referecne layers;.

Definition at line 226 of file GoldenPattern.cc.

                               {
   for (unsigned int iRefLayer = 0; iRefLayer < meanDistPhi[0].size(); ++iRefLayer) {
     for (unsigned int iLayer = 0; iLayer < meanDistPhi.size(); ++iLayer) {
       if (!meanDistPhiCounts.empty() && meanDistPhiCounts[iLayer][iRefLayer])
         return true;
     }
   }
   return false;
 }

References meanDistPhi, and meanDistPhiCounts.

◆ key()

Key GoldenPattern::key ( ) const

inline

Definition at line 56 of file GoldenPattern.h.

56 { return theKey; }

References theKey.

Referenced by OMTFProcessor::addGP(), XMLConfigWriter::writeGPData(), and OMTFPatternMaker::writeMergedGPs().

◆ meanDistPhiValue()

int GoldenPattern::meanDistPhiValue	(	unsigned int	iLayer,
		unsigned int	iRefLayer
	)		const

inline

Definition at line 66 of file GoldenPattern.h.

66 { return meanDistPhi[iLayer][iRefLayer]; }

References meanDistPhi.

Referenced by XMLConfigWriter::writeGPData().

◆ normalise()

void GoldenPattern::normalise ( unsigned int nPdfAddrBits )

Normalise event counts in mean dist phi, and pdf vectors to get the real values of meand dist phi and probability. The pdf width is passed to this method, since the width stored in configuration is extended during the pattern making phase.

Mean dist phi

Probabilities. Normalise and change from float to integer values

If there are only a few counts in given measurement layer, set pdf value to 0

Digitisation Values remapped 0->std::pow(2,nPdfValBits) minPlog->0

Make sure digitised value is saved using nBitsVal bits

Shift pdf index by meanDistPhi

Definition at line 152 of file GoldenPattern.cc.

                                                        {
   for (unsigned int iRefLayer = 0; iRefLayer < pdfAllRef[0].size(); ++iRefLayer) {
     for (unsigned int iLayer = 0; iLayer < pdfAllRef.size(); ++iLayer) {
       for (unsigned int iPdf = 0; iPdf < pdfAllRef[iLayer][iRefLayer].size(); ++iPdf) {
         float pVal = log((float)pdfAllRef[iLayer][iRefLayer][iPdf] /
                          meanDistPhiCounts[myOmtfConfig->getRefToLogicNumber()[iRefLayer]][iRefLayer]);
         if (pVal < log(myOmtfConfig->minPdfVal()))
           continue;
         meanDistPhi[iLayer][iRefLayer] +=
             (iPdf - exp2(myOmtfConfig->nPdfAddrBits() - 1)) * pdfAllRef[iLayer][iRefLayer][iPdf];
         if ((int)iLayer != myOmtfConfig->getRefToLogicNumber()[iRefLayer])
           meanDistPhiCounts[iLayer][iRefLayer] += pdfAllRef[iLayer][iRefLayer][iPdf];
       }
     }
   }
  
   for (unsigned int iRefLayer = 0; iRefLayer < meanDistPhi[0].size(); ++iRefLayer) {
     for (unsigned int iLayer = 0; iLayer < meanDistPhi.size(); ++iLayer) {
       if (!meanDistPhiCounts.empty() && meanDistPhiCounts[iLayer][iRefLayer]) {
         if (meanDistPhiCounts[iLayer][iRefLayer] < 1000)
           meanDistPhi[iLayer][iRefLayer] = 0;
         else
           meanDistPhi[iLayer][iRefLayer] =
               rint((float)meanDistPhi[iLayer][iRefLayer] / meanDistPhiCounts[iLayer][iRefLayer]);
       }
     }
   }
   const float minPlog = log(myOmtfConfig->minPdfVal());
   const unsigned int nPdfValBits = myOmtfConfig->nPdfValBits();
   float pVal;
   int digitisedVal, truncatedValue;
   for (unsigned int iRefLayer = 0; iRefLayer < pdfAllRef[0].size(); ++iRefLayer) {
     for (unsigned int iLayer = 0; iLayer < pdfAllRef.size(); ++iLayer) {
       for (unsigned int iPdf = 0; iPdf < pdfAllRef[iLayer][iRefLayer].size(); ++iPdf) {
         if (!meanDistPhiCounts[myOmtfConfig->getRefToLogicNumber()[iRefLayer]][iRefLayer] ||
             !pdfAllRef[iLayer][iRefLayer][iPdf])
           continue;
         pVal = log((float)pdfAllRef[iLayer][iRefLayer][iPdf] /
                    meanDistPhiCounts[myOmtfConfig->getRefToLogicNumber()[iRefLayer]][iRefLayer]);
         if ((pVal < minPlog || meanDistPhiCounts[iLayer][iRefLayer] < 1000)) {
           pdfAllRef[iLayer][iRefLayer][iPdf] = 0;
           continue;
         }
         digitisedVal = rint((std::pow(2, nPdfValBits) - 1) - (pVal / minPlog) * (std::pow(2, nPdfValBits) - 1));
         truncatedValue = 0 | (digitisedVal & ((int)pow(2, nPdfValBits) - 1));
         pdfAllRef[iLayer][iRefLayer][iPdf] = truncatedValue;
       }
     }
   }
  
   vector3D pdfAllRefTmp = pdfAllRef;
   for (unsigned int iRefLayer = 0; iRefLayer < pdfAllRef[0].size(); ++iRefLayer) {
     for (unsigned int iLayer = 0; iLayer < pdfAllRef.size(); ++iLayer) {
       for (unsigned int iPdf = 0; iPdf < pdfAllRef[iLayer][iRefLayer].size(); ++iPdf) {
         pdfAllRef[iLayer][iRefLayer][iPdf] = 0;
         int index =
             iPdf - exp2(myOmtfConfig->nPdfAddrBits() - 1) - meanDistPhi[iLayer][iRefLayer] + exp2(nPdfAddrBits - 1);
         if (index < 0 || index > exp2(nPdfAddrBits) - 1)
           continue;
         pdfAllRef[iLayer][iRefLayer][index] = pdfAllRefTmp[iLayer][iRefLayer][iPdf];
       }
     }
   }
 }

References OMTFConfiguration::getRefToLogicNumber(), createfilelist::int, dqm-mbProfile::log, meanDistPhi, meanDistPhiCounts, OMTFConfiguration::minPdfVal(), myOmtfConfig, OMTFConfiguration::nPdfAddrBits(), OMTFConfiguration::nPdfValBits(), pdfAllRef, and funct::pow().

◆ pdfValue()

int GoldenPattern::pdfValue	(	unsigned int	iLayer,
		unsigned int	iRefLayer,
		unsigned int	iBin
	)		const

inline

Definition at line 68 of file GoldenPattern.h.

                                                                                      {
     return pdfAllRef[iLayer][iRefLayer][iBin];
   }

References pdfAllRef.

Referenced by OMTFProcessor::shiftGP(), and XMLConfigWriter::writeGPData().

◆ process1Layer1RefLayer()

GoldenPattern::layerResult GoldenPattern::process1Layer1RefLayer	(	unsigned int	iRefLayer,
		unsigned int	iLayer,
		const int	refPhi,
		const OMTFinput::vector1D &	layerHits
	)

Process single measurement layer with a single ref layer Method should be thread safe

Select hit closest to the mean of probability distribution in given layer

Check if phiDist is within pdf range -63 +63

Shift phidist, so 0 is at the middle of the range

Definition at line 13 of file GoldenPattern.cc.

                                                                                                      {
   GoldenPattern::layerResult aResult;
  
   int phiMean = meanDistPhi[iLayer][iRefLayer];
   int phiDist = exp2(myOmtfConfig->nPdfAddrBits());
   for (auto itHit : layerHits) {
     if (itHit >= (int)myOmtfConfig->nPhiBins())
       continue;
     if (abs(itHit - phiMean - phiRefHit) < abs(phiDist))
       phiDist = itHit - phiMean - phiRefHit;
   }
  
   if (abs(phiDist) > (exp2(myOmtfConfig->nPdfAddrBits() - 1) - 1))
     return aResult;
  
   phiDist += exp2(myOmtfConfig->nPdfAddrBits() - 1);
  
   int pdfVal = pdfAllRef[iLayer][iRefLayer][phiDist];
  
   return GoldenPattern::layerResult(pdfVal, pdfVal > 0);
 }

References funct::abs(), meanDistPhi, myOmtfConfig, OMTFConfiguration::nPdfAddrBits(), OMTFConfiguration::nPhiBins(), and pdfAllRef.

◆ propagateRefPhi()

int GoldenPattern::propagateRefPhi	(	int	phiRef,
		int	etaRef,
		unsigned int	iRefLayer
	)

Propagate phi from given reference layer to MB2 or ME2 ME2 is used if eta of reference hit is larger than 1.1 expressed in ingerer MicroGMT scale: 1.1/2.61*240 = 101

Definition at line 43 of file GoldenPattern.cc.

                                                                                  {
   unsigned int iLayer = 2;  //MB2
   //if(etaRef>101) iLayer = 7;//RE2
   return phiRef + meanDistPhi[iLayer][iRefLayer];
 }

References meanDistPhi.

◆ reset()

void GoldenPattern::reset ( void )

Reset contents of all data vectors, keeping the vectors size.

Definition at line 136 of file GoldenPattern.cc.

                           {
   GoldenPattern::vector1D meanDistPhi1D(myOmtfConfig->nRefLayers());
   GoldenPattern::vector2D meanDistPhi2D(myOmtfConfig->nLayers());
   meanDistPhi2D.assign(myOmtfConfig->nLayers(), meanDistPhi1D);
   meanDistPhi = meanDistPhi2D;
   meanDistPhiCounts = meanDistPhi2D;
  
   GoldenPattern::vector1D pdf1D(exp2(myOmtfConfig->nPdfAddrBits()));
   GoldenPattern::vector2D pdf2D(myOmtfConfig->nRefLayers());
   GoldenPattern::vector3D pdf3D(myOmtfConfig->nLayers());
  
   pdf2D.assign(myOmtfConfig->nRefLayers(), pdf1D);
   pdfAllRef.assign(myOmtfConfig->nLayers(), pdf2D);
 }

References meanDistPhi, meanDistPhiCounts, myOmtfConfig, OMTFConfiguration::nLayers(), OMTFConfiguration::nPdfAddrBits(), OMTFConfiguration::nRefLayers(), and pdfAllRef.

Referenced by OMTFPatternMaker::endJob().

◆ setMeanDistPhi()

void GoldenPattern::setMeanDistPhi ( const vector2D & aMeanDistPhi )

inline

Definition at line 58 of file GoldenPattern.h.

58 { meanDistPhi = aMeanDistPhi; }

References meanDistPhi.

Referenced by OMTFProcessor::averagePatterns(), and OMTFProcessor::configure().

◆ setPdf()

void GoldenPattern::setPdf ( const vector3D & aPdf )

inline

Definition at line 64 of file GoldenPattern.h.

64 { pdfAllRef = aPdf; }

References pdfAllRef.

Referenced by OMTFProcessor::configure(), and OMTFProcessor::shiftGP().

Friends And Related Function Documentation

◆ operator<<

std::ostream& operator<<	(	std::ostream &	out,
		const GoldenPattern &	aPattern
	)

friend

Definition at line 89 of file GoldenPattern.cc.

                                                                        {
   out << "GoldenPattern " << aPattern.theKey << std::endl;
   out << "Number of reference layers: " << aPattern.meanDistPhi[0].size()
       << ", number of measurement layers: " << aPattern.pdfAllRef.size() << std::endl;
  
   if (aPattern.meanDistPhi.empty())
     return out;
   if (aPattern.pdfAllRef.empty())
     return out;
  
   out << "Mean dist phi per layer:" << std::endl;
   for (unsigned int iRefLayer = 0; iRefLayer < aPattern.meanDistPhi[0].size(); ++iRefLayer) {
     out << "Ref layer: " << iRefLayer << " (";
     for (unsigned int iLayer = 0; iLayer < aPattern.meanDistPhi.size(); ++iLayer) {
       out << std::setw(3) << aPattern.meanDistPhi[iLayer][iRefLayer] << "\t";
     }
     out << ")" << std::endl;
   }
  
   if (!aPattern.meanDistPhiCounts.empty()) {
     out << "Counts number per layer:" << std::endl;
     for (unsigned int iRefLayer = 0; iRefLayer < aPattern.meanDistPhi[0].size(); ++iRefLayer) {
       out << "Ref layer: " << iRefLayer << " (";
       for (unsigned int iLayer = 0; iLayer < aPattern.meanDistPhi.size(); ++iLayer) {
         out << aPattern.meanDistPhiCounts[iLayer][iRefLayer] << "\t";
       }
       out << ")" << std::endl;
     }
   }
  
   unsigned int nPdfAddrBits = 7;
   out << "PDF per layer:" << std::endl;
   for (unsigned int iRefLayer = 0; iRefLayer < aPattern.pdfAllRef[0].size(); ++iRefLayer) {
     out << "Ref layer: " << iRefLayer;
     for (unsigned int iLayer = 0; iLayer < aPattern.pdfAllRef.size(); ++iLayer) {
       out << ", measurement layer: " << iLayer << std::endl;
       for (unsigned int iPdf = 0; iPdf < exp2(nPdfAddrBits); ++iPdf) {
         out << std::setw(2) << aPattern.pdfAllRef[iLayer][iRefLayer][iPdf] << " ";
       }
       out << std::endl;
     }
   }
  
   return out;
 }

Member Data Documentation

◆ meanDistPhi

vector2D GoldenPattern::meanDistPhi

private

Mean positions in each layer First index: measurement layer number Second index: refLayer number

Definition at line 114 of file GoldenPattern.h.

Referenced by getMeanDistPhi(), hasCounts(), meanDistPhiValue(), normalise(), operator<<(), process1Layer1RefLayer(), propagateRefPhi(), reset(), and setMeanDistPhi().

◆ meanDistPhiCounts

vector2D GoldenPattern::meanDistPhiCounts

private

Vector holding number of counts. Used for making the patterns

Definition at line 118 of file GoldenPattern.h.

Referenced by addCount(), hasCounts(), normalise(), operator<<(), and reset().

◆ myOmtfConfig

const OMTFConfiguration* GoldenPattern::myOmtfConfig

private

Definition at line 120 of file GoldenPattern.h.

Referenced by addCount(), normalise(), process1Layer1RefLayer(), and reset().

◆ pdfAllRef

vector3D GoldenPattern::pdfAllRef

private

Distributions for all reference layers First index: measurement layer number Second index: refLayer number Third index: pdf bin number within layer

Definition at line 109 of file GoldenPattern.h.

Referenced by addCount(), getPdf(), normalise(), operator<<(), pdfValue(), process1Layer1RefLayer(), reset(), and setPdf().

◆ theKey

Key GoldenPattern::theKey

private

Pattern kinematical identification (iEta,iPt,iCharge)

Definition at line 103 of file GoldenPattern.h.

Referenced by key(), and operator<<().

Public Types

Public Member Functions

Private Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ layerResult

◆ vector1D

◆ vector2D

◆ vector3D

Constructor & Destructor Documentation

◆ GoldenPattern()

Member Function Documentation

◆ addCount()

◆ getMeanDistPhi()

◆ getPdf()

◆ hasCounts()

◆ key()

◆ meanDistPhiValue()

◆ normalise()

◆ pdfValue()

◆ process1Layer1RefLayer()

◆ propagateRefPhi()

◆ reset()

◆ setMeanDistPhi()

◆ setPdf()

Friends And Related Function Documentation

◆ operator<<

Member Data Documentation

◆ meanDistPhi

◆ meanDistPhiCounts

◆ myOmtfConfig

◆ pdfAllRef

◆ theKey