GoldenPattern Class Reference

#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

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

Definition at line 49 of file GoldenPattern.h.

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

Definition at line 46 of file GoldenPattern.h.

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

Definition at line 47 of file GoldenPattern.h.

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

Definition at line 48 of file GoldenPattern.h.

Constructor & Destructor Documentation

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

inline

Definition at line 54 of file GoldenPattern.h.

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

Member Function Documentation

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.

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

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

const vector2D& GoldenPattern::getMeanDistPhi ( ) const

inline

Definition at line 60 of file GoldenPattern.h.

Referenced by OMTFProcessor::averagePatterns().

{ return meanDistPhi; }

const vector3D& GoldenPattern::getPdf ( ) const

inline

Definition at line 62 of file GoldenPattern.h.

Referenced by OMTFProcessor::shiftGP().

{ return pdfAllRef; }

bool GoldenPattern::hasCounts

(

)

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

Definition at line 226 of file GoldenPattern.cc.

References meanDistPhi, and meanDistPhiCounts.

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

Key GoldenPattern::key ( ) const

inline

Definition at line 56 of file GoldenPattern.h.

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

{ return theKey; }

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

inline

Definition at line 66 of file GoldenPattern.h.

Referenced by XMLConfigWriter::writeGPData().

{ return meanDistPhi[iLayer][iRefLayer]; }

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.

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

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

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

inline

Definition at line 68 of file GoldenPattern.h.

References Key::operator<<, MillePedeFileConverter_cfg::out, and reset().

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

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

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.

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

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

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.

References meanDistPhi.

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

void GoldenPattern::reset

(

void

)

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

Definition at line 136 of file GoldenPattern.cc.

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

Referenced by OMTFPatternMaker::endJob(), and OMTFPatternMaker::writeMergedGPs().

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

void GoldenPattern::setMeanDistPhi ( const vector2D &  aMeanDistPhi )

inline

Definition at line 58 of file GoldenPattern.h.

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

{ meanDistPhi = aMeanDistPhi; }

void GoldenPattern::setPdf ( const vector3D &  aPdf )

inline

Definition at line 64 of file GoldenPattern.h.

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

{ pdfAllRef = aPdf; }

Friends And Related Function Documentation

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

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 hasCounts(), normalise(), operator<<(), process1Layer1RefLayer(), propagateRefPhi(), and reset().

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().

const OMTFConfiguration* GoldenPattern::myOmtfConfig

private

Definition at line 120 of file GoldenPattern.h.

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

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(), normalise(), operator<<(), process1Layer1RefLayer(), and reset().

Key GoldenPattern::theKey

private

Pattern kinematical identification (iEta,iPt,iCharge)

Definition at line 103 of file GoldenPattern.h.

Referenced by operator<<().