GoldenPattern Class Reference

#include <GoldenPattern.h>

Inheritance diagram for GoldenPattern:

Public Types
typedef std::pair< int, bool >	layerResult
typedef boost::multi_array < short, 3 >	meanDistPhiArrayType
typedef boost::multi_array < PdfValueType, 3 >	pdfArrayType
typedef std::vector< int >	vector1D
typedef std::vector< vector1D >	vector2D
typedef std::vector< vector2D >	vector3D
Public Types inherited from GoldenPatternBase
typedef boost::multi_array < GoldenPatternResult, 2 >	resultsArrayType
typedef std::vector< int >	vector1D

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...
int	getDistPhiBitShift (unsigned int iLayer, unsigned int iRefLayer) const override
virtual const meanDistPhiArrayType &	getMeanDistPhi () const
const vector2D &	getMeanDistPhi () const
virtual const pdfArrayType &	getPdf () const
const vector3D &	getPdf () const
	GoldenPattern (const Key &aKey, unsigned int nLayers, unsigned int nRefLayers, unsigned int nPdfAddrBits)
	GoldenPattern (const Key &aKey, const OMTFConfiguration *omtfConfig)
	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, int refLayerPhiB=0) const override
int	meanDistPhiValue (unsigned int iLayer, unsigned int iRefLayer) const
void	normalise (unsigned int nPdfAddrBits)
PdfValueType	pdfValue (unsigned int iLayer, unsigned int iRefLayer, unsigned int iBin, int refLayerPhiB=0) const override
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) override
int	propagateRefPhi (int phiRef, int etaRef, unsigned int iRefLayer)
virtual void	reset ()
	Reset contents of all data vectors, keeping the vectors size. More...
void	reset ()
	Reset contents of all data vectors, keeping the vectors size. More...
void	setDistPhiBitShift (int value, unsigned int iLayer, unsigned int iRefLayer) override
virtual void	setMeanDistPhi (const meanDistPhiArrayType &aMeanDistPhi)
void	setMeanDistPhi (const vector2D &aMeanDistPhi)
void	setMeanDistPhiValue (int value, unsigned int iLayer, unsigned int iRefLayer, unsigned int paramIndex=0) override
virtual void	setPdf (pdfArrayType &aPdf)
void	setPdf (const vector3D &aPdf)
void	setPdfValue (PdfValueType value, unsigned int iLayer, unsigned int iRefLayer, unsigned int iBin, int refLayerPhiB=0) override
	~GoldenPattern () override
Public Member Functions inherited from GoldenPatternBase
virtual void	finalise (unsigned int procIndx)
	last step of the event processing, before sorting and ghost busting More...
const OMTFConfiguration *	getConfig () const
resultsArrayType &	getResults ()
	GoldenPatternBase (const Key &aKey)
	GoldenPatternBase (const Key &aKey, const OMTFConfiguration *omtfConfig)
virtual Key &	key ()
virtual StubResult	process1Layer1RefLayer (unsigned int iRefLayer, unsigned int iLayer, MuonStubPtrs1D layerStubs, const MuonStubPtr refStub)
virtual void	setConfig (const OMTFConfiguration *omtfConfig)
virtual	~GoldenPatternBase ()

Protected Attributes
boost::multi_array< short, 2 >	distPhiBitShift
meanDistPhiArrayType	meanDistPhi
pdfArrayType	pdfAllRef
Protected Attributes inherited from GoldenPatternBase
const OMTFConfiguration *	myOmtfConfig
resultsArrayType	results
Key	theKey
	Pattern kinematic identification (iEta,iPt,iCharge) More...

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)
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 boost::multi_array<short, 3> GoldenPattern::meanDistPhiArrayType

Definition at line 18 of file GoldenPattern.h.

typedef boost::multi_array<PdfValueType, 3> GoldenPattern::pdfArrayType

Definition at line 17 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) {}

GoldenPattern::GoldenPattern ( const Key &  aKey, unsigned int  nLayers, unsigned int  nRefLayers, unsigned int  nPdfAddrBits  )

inline

Definition at line 22 of file GoldenPattern.h.

References reset().

      : GoldenPatternBase(aKey),
        pdfAllRef(boost::extents[nLayers][nRefLayers][1 << nPdfAddrBits]),
        meanDistPhi(boost::extents[nLayers][nRefLayers][2]),
        distPhiBitShift(boost::extents[nLayers][nRefLayers]) {
    reset();
  }

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

inline

Definition at line 30 of file GoldenPattern.h.

References reset().

      : GoldenPatternBase(aKey, omtfConfig),
        pdfAllRef(boost::extents[omtfConfig->nLayers()][omtfConfig->nRefLayers()][omtfConfig->nPdfBins()]),
        meanDistPhi(boost::extents[omtfConfig->nLayers()][omtfConfig->nRefLayers()][2]),
        distPhiBitShift(boost::extents[omtfConfig->nLayers()][omtfConfig->nRefLayers()]) {
    reset();
  }

GoldenPattern::~GoldenPattern ( )

inlineoverride

Definition at line 38 of file GoldenPattern.h.

{};

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

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

inlineoverridevirtual

Implements GoldenPatternBase.

Definition at line 72 of file GoldenPattern.h.

References distPhiBitShift.

Referenced by XMLConfigWriter::writeGPData().

                                                                                     {
    return distPhiBitShift[iLayer][iRefLayer];
  }

virtual const meanDistPhiArrayType& GoldenPattern::getMeanDistPhi ( ) const

inlinevirtual

Definition at line 42 of file GoldenPattern.h.

References meanDistPhi.

{ return meanDistPhi; }

const vector2D& GoldenPattern::getMeanDistPhi ( ) const

inline

Definition at line 60 of file GoldenPattern.h.

References meanDistPhi.

Referenced by OMTFProcessor< GoldenPatternType >::averagePatterns(), and XMLConfigWriter::writeGPData().

{ return meanDistPhi; }

virtual const pdfArrayType& GoldenPattern::getPdf ( ) const

inlinevirtual

Definition at line 44 of file GoldenPattern.h.

References pdfAllRef.

{ return pdfAllRef; }

const vector3D& GoldenPattern::getPdf ( ) const

inline

Definition at line 62 of file GoldenPattern.h.

References pdfAllRef.

Referenced by OMTFProcessor< GoldenPatternType >::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.

References theKey.

Referenced by OMTFProcessor< GoldenPatternType >::addGP(), PatternOptimizerBase::observeEventEnd(), XMLConfigWriter::writeGPData(), and OMTFPatternMaker::writeMergedGPs().

{ return theKey; }

int GoldenPattern::meanDistPhiValue ( unsigned int  iLayer, unsigned int  iRefLayer, int  refLayerPhiB = 0  ) const

overridevirtual

Implements GoldenPatternBase.

Definition at line 9 of file GoldenPattern.cc.

References meanDistPhi, myOmtfConfig, and OMTFConfiguration::nPdfAddrBits().

                                                                                                       {
  //return meanDistPhi[iLayer][iRefLayer][0];
  return (((meanDistPhi[iLayer][iRefLayer][1] * refLayerPhiB) >> myOmtfConfig->nPdfAddrBits()) +
          meanDistPhi[iLayer][iRefLayer][0]);
  //assumes that the meanDistPhi[1] is float alpha from the fit to the phiB-phi distribution multiplied by 2^myOmtfConfig->nPdfAddrBits()
}

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

inline

Definition at line 66 of file GoldenPattern.h.

References meanDistPhi.

Referenced by PatternGenerator::updateStat(), and 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(), 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];
      }
    }
  }
}

PdfValueType GoldenPattern::pdfValue ( unsigned int  iLayer, unsigned int  iRefLayer, unsigned int  iBin, int  refLayerPhiB = 0  ) const

inlineoverridevirtual

Implements GoldenPatternBase.

Definition at line 50 of file GoldenPattern.h.

References pdfAllRef.

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

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

inline

Definition at line 68 of file GoldenPattern.h.

References pdfAllRef.

Referenced by OMTFProcessor< GoldenPatternType >::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  )

overridevirtual

Propagate phi from given reference layer to MB2 or ME2 expressed in integer MicroGMT scale: 1.1/2.61*240 = 101

Implements GoldenPatternBase.

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

virtual

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

Implements GoldenPatternBase.

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

virtual void GoldenPattern::reset ( )

virtual

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

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(), GoldenPattern(), 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::setDistPhiBitShift ( int  value, unsigned int  iLayer, unsigned int  iRefLayer  )

inlineoverridevirtual

Implements GoldenPatternBase.

Definition at line 76 of file GoldenPattern.h.

References distPhiBitShift, and relativeConstraints::value.

                                                                                           {
    distPhiBitShift[iLayer][iRefLayer] = value;
  }

virtual void GoldenPattern::setMeanDistPhi ( const meanDistPhiArrayType &  aMeanDistPhi )

inlinevirtual

Definition at line 40 of file GoldenPattern.h.

References meanDistPhi.

{ meanDistPhi = aMeanDistPhi; }

void GoldenPattern::setMeanDistPhi ( const vector2D &  aMeanDistPhi )

inline

Definition at line 58 of file GoldenPattern.h.

References meanDistPhi.

Referenced by OMTFProcessor< GoldenPatternType >::averagePatterns(), and OMTFProcessor< GoldenPatternType >::configure().

{ meanDistPhi = aMeanDistPhi; }

void GoldenPattern::setMeanDistPhiValue ( int  value, unsigned int  iLayer, unsigned int  iRefLayer, unsigned int  paramIndex = 0  )

inlineoverridevirtual

Implements GoldenPatternBase.

Definition at line 57 of file GoldenPattern.h.

References meanDistPhi, and relativeConstraints::value.

                                                                 {
    meanDistPhi[iLayer][iRefLayer][paramIndex] = value;
  }

virtual void GoldenPattern::setPdf ( pdfArrayType &  aPdf )

inlinevirtual

Definition at line 46 of file GoldenPattern.h.

References pdfAllRef.

{ pdfAllRef = aPdf; }

void GoldenPattern::setPdf ( const vector3D &  aPdf )

inline

Definition at line 64 of file GoldenPattern.h.

References pdfAllRef.

Referenced by OMTFProcessor< GoldenPatternType >::configure(), and OMTFProcessor< GoldenPatternType >::shiftGP().

{ pdfAllRef = aPdf; }

void GoldenPattern::setPdfValue ( PdfValueType  value, unsigned int  iLayer, unsigned int  iRefLayer, unsigned int  iBin, int  refLayerPhiB = 0  )

inlineoverridevirtual

Implements GoldenPatternBase.

Definition at line 64 of file GoldenPattern.h.

References pdfAllRef, and relativeConstraints::value.

                                                  {
    pdfAllRef[iLayer][iRefLayer][iBin] = value;
  }

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

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

boost::multi_array<short, 2> GoldenPattern::distPhiBitShift

protected

distPhi resolution can be reduced to reduce the number of bit on the LUT input distPhi = distPhi<<distPhiBitShift[layer]refLayer First index: measurement layer number Second index: refLayer number

Definition at line 106 of file GoldenPattern.h.

Referenced by getDistPhiBitShift(), and setDistPhiBitShift().

meanDistPhiArrayType GoldenPattern::meanDistPhi

protected

Mean positions in each layer First index: measurement layer number Second index: refLayer number Third index: index = 0 - a0, index = 1 - a1 for the linear fit meanDistPhi = a0 + a1 * phi_b

Definition at line 100 of file GoldenPattern.h.

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(), setMeanDistPhi(), and setMeanDistPhiValue().

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

pdfArrayType GoldenPattern::pdfAllRef

protected

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

Definition at line 94 of file GoldenPattern.h.

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(), GoldenPatternWithStat::iniStatisitics(), normalise(), operator<<(), pdfValue(), process1Layer1RefLayer(), reset(), setPdf(), and setPdfValue().

Key GoldenPattern::theKey

private

Pattern kinematical identification (iEta,iPt,iCharge)

Definition at line 103 of file GoldenPattern.h.

Referenced by key(), operator<<(), GoldenPatternWithStat::setKeyNumber(), and GoldenPatternWithStat::setKeyPt().