OMTFProcessor< GoldenPatternType > Class Template Reference

#include <OMTFProcessor.h>

Inheritance diagram for OMTFProcessor< GoldenPatternType >:

Public Types
typedef std::map< Key, OMTFResult >	resultsMap

Public Member Functions
void	averagePatterns (int charge)
bool	configure (const OMTFConfiguration *omtfParams, const L1TMuonOverlapParams *omtfPatterns)
	Fill GP map with patterns from CondFormats object. More...
int	extrapolateDtPhiB (const MuonStubPtr &refStub, const MuonStubPtr &targetStub, unsigned int targetLayer, const OMTFConfiguration *omtfConfig)
int	extrapolateDtPhiBFixedPoint (const int &refLogicLayer, const int &refPhi, const int &refPhiB, unsigned int targetLayer, const int &targetStubPhi, const int &targetStubQuality, const int &targetStubEta, const int &targetStubR, const OMTFConfiguration *omtfConfig)
int	extrapolateDtPhiBFloatPoint (const int &refLogicLayer, const int &refPhi, const int &refPhiB, unsigned int targetLayer, const int &targetStubPhi, const int &targetStubQuality, const int &targetStubEta, const int &targetStubR, const OMTFConfiguration *omtfConfig)
	Fill GP vec with patterns from CondFormats object. More...
void	fillCounts (unsigned int iProcessor, const OMTFinput &aInput, const SimTrack *aSimMuon)
std::vector< l1t::RegionalMuonCand >	getFinalcandidates (unsigned int iProcessor, l1t::tftype mtfType, const AlgoMuons &algoCands) override
const std::map< Key, GoldenPattern * > &	getPatterns () const
	Return map of GoldenPatterns. More...
AlgoMuons	ghostBust (AlgoMuons refHitCands, int charge=0) override
void	loadExtrapolFactors (const std::string &filename)
	OMTFProcessor ()
	OMTFProcessor (OMTFConfiguration *omtfConfig, const edm::ParameterSet &edmCfg, edm::EventSetup const &evSetup, const L1TMuonOverlapParams *omtfPatterns)
	OMTFProcessor (OMTFConfiguration *omtfConfig, const edm::ParameterSet &edmCfg, edm::EventSetup const &evSetup, GoldenPatternVec< GoldenPatternType > &&gps)
void	printInfo () const override
const std::vector< OMTFProcessor::resultsMap > &	processInput (unsigned int iProcessor, const OMTFinput &aInput)
void	processInput (unsigned int iProcessor, l1t::tftype mtfType, const OMTFinput &aInput, std::vector< std::unique_ptr< IOMTFEmulationObserver > > &observers) override
std::vector< l1t::RegionalMuonCand >	run (unsigned int iProcessor, l1t::tftype mtfType, int bx, OMTFinputMaker *inputMaker, std::vector< std::unique_ptr< IOMTFEmulationObserver > > &observers) override
void	saveExtrapolFactors ()
void	setGhostBuster (IGhostBuster *ghostBuster) override
	allows to use other IGhostBuster implementation than the default one More...
virtual void	setPtAssignment (PtAssignmentBase *ptAssignment)
virtual void	setSorter (SorterBase< GoldenPatternType > *sorter)
	allows to use other sorter implementation than the default one More...
AlgoMuons	sortResults (unsigned int iProcessor, l1t::tftype mtfType, int charge=0) override
	~OMTFProcessor ()
	~OMTFProcessor () override
Public Member Functions inherited from ProcessorBase< GoldenPatternType >
virtual void	addGP (GoldenPatternType *aGP)
	Add GoldenPattern to pattern vec. More...
virtual bool	configure (OMTFConfiguration *omtfParams, const L1TMuonOverlapParams *omtfPatterns)
	Fill GP vec with patterns from CondFormats object. More...
const std::vector< OMTFConfiguration::PatternPt > &	getPatternPtRange () const
virtual GoldenPatternVec< GoldenPatternType > &	getPatterns ()
	Return vector of GoldenPatterns. More...
virtual void	initPatternPtRange (bool firstPatFrom0)
	ProcessorBase (OMTFConfiguration *omtfConfig, const L1TMuonOverlapParams *omtfPatterns)
	ProcessorBase (OMTFConfiguration *omtfConfig, GoldenPatternVec< GoldenPatternType > &&gps)
virtual	~ProcessorBase ()
Public Member Functions inherited from IProcessorEmulator
virtual	~IProcessorEmulator ()

Private Member Functions
bool	addGP (GoldenPattern *aGP)
bool	checkHitPatternValidity (unsigned int hits) override
void	fillInputRange (unsigned int iProcessor, unsigned int iCone, const OMTFinput &aInput)
void	fillInputRange (unsigned int iProcessor, unsigned int iCone, unsigned int iRefLayer, unsigned int iHit)
virtual void	init (const edm::ParameterSet &edmCfg, edm::EventSetup const &evSetup)
void	resetConfiguration ()
	Reset all configuration parameters. More...
OMTFinput::vector1D	restrictInput (unsigned int iProcessor, unsigned int iCone, unsigned int iLayer, const OMTFinput::vector1D &layerHits)
void	shiftGP (GoldenPattern *aGP, const GoldenPattern::vector2D &meanDistPhiNew, const GoldenPattern::vector2D &meanDistPhiOld)

Private Attributes
std::vector< std::vector< std::map< int, double > > >	extrapolFactors
std::vector< std::vector< std::map< int, int > > >	extrapolFactorsNorm
int	extrapolMultiplier = 128
std::unique_ptr< IGhostBuster >	ghostBuster
const OMTFConfiguration *	myOmtfConfig
std::vector< OMTFProcessor::resultsMap >	myResults
PtAssignmentBase *	ptAssignment = nullptr
std::unique_ptr< SorterBase< GoldenPatternType > >	sorter
std::map< Key, GoldenPattern * >	theGPs
	Map holding Golden Patterns. More...
bool	useEndcapStubsRInExtr = false
bool	useFloatingPointExtrapolation = false
bool	useStubQualInExtr = false

Additional Inherited Members
Protected Member Functions inherited from ProcessorBase< GoldenPatternType >
virtual MuonStubPtrs1D	restrictInput (unsigned int iProcessor, unsigned int iCone, unsigned int iLayer, const OMTFinput &input)
Protected Attributes inherited from ProcessorBase< GoldenPatternType >
const OMTFConfiguration *	myOmtfConfig
std::vector< OMTFConfiguration::PatternPt >	patternPts
GoldenPatternVec< GoldenPatternType >	theGPs
	vector holding Golden Patterns More...

Detailed Description

template<class GoldenPatternType>
class OMTFProcessor< GoldenPatternType >

Definition at line 19 of file OMTFProcessor.h.

Member Typedef Documentation

resultsMap

template<class GoldenPatternType >

typedef std::map<Key, OMTFResult> OMTFProcessor< GoldenPatternType >::resultsMap

Definition at line 21 of file OMTFProcessor.h.

Constructor & Destructor Documentation

OMTFProcessor() [1/3]

template<class GoldenPatternType >

OMTFProcessor< GoldenPatternType >::OMTFProcessor ( )

inline

Definition at line 23 of file OMTFProcessor.h.

{};

~OMTFProcessor() [1/2]

template<class GoldenPatternType >

OMTFProcessor< GoldenPatternType >::~OMTFProcessor

(

)

Definition at line 21 of file OMTFProcessor.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and OMTFProcessor< GoldenPatternType >::theGPs.

                              {
  for (auto it : theGPs)
    delete it.second;
}

OMTFProcessor() [2/3]

template<class GoldenPatternType >

OMTFProcessor< GoldenPatternType >::OMTFProcessor	(	OMTFConfiguration *	omtfConfig,
		const edm::ParameterSet &	edmCfg,
		edm::EventSetup const &	evSetup,
		const L1TMuonOverlapParams *	omtfPatterns
	)

Definition at line 35 of file OMTFProcessor.cc.

References OMTFProcessor< GoldenPatternType >::init().

    : ProcessorBase<GoldenPatternType>(omtfConfig, omtfPatterns) {
  init(edmCfg, evSetup);
};

OMTFProcessor() [3/3]

template<class GoldenPatternType >

OMTFProcessor< GoldenPatternType >::OMTFProcessor	(	OMTFConfiguration *	omtfConfig,
		const edm::ParameterSet &	edmCfg,
		edm::EventSetup const &	evSetup,
		GoldenPatternVec< GoldenPatternType > &&	gps
	)

Definition at line 44 of file OMTFProcessor.cc.

References OMTFProcessor< GoldenPatternType >::init().

    : ProcessorBase<GoldenPatternType>(omtfConfig, std::forward<GoldenPatternVec<GoldenPatternType> >(gps)) {
  init(edmCfg, evSetup);
};

~OMTFProcessor() [2/2]

template<class GoldenPatternType >

OMTFProcessor< GoldenPatternType >::~OMTFProcessor ( )

override

Member Function Documentation

addGP()

template<class GoldenPatternType >

bool OMTFProcessor< GoldenPatternType >::addGP ( GoldenPattern *  aGP )

private

Add GoldenPattern to pattern map. If GP key already exists in map, a new entry is ignored

Definition at line 95 of file OMTFProcessor.cc.

References OMTFResult::configure(), Exception, GoldenPattern::key(), OMTFProcessor< GoldenPatternType >::myOmtfConfig, OMTFProcessor< GoldenPatternType >::myResults, and OMTFProcessor< GoldenPatternType >::theGPs.

Referenced by OMTFProcessor< GoldenPatternType >::configure().

                                            {
  if (theGPs.find(aGP->key()) != theGPs.end()) {
    throw cms::Exception("Corrupted Golden Patterns data")
        << "OMTFProcessor::addGP(...) "
        << " Reading two Golden Patterns with the same key: " << aGP->key() << std::endl;
  } else
    theGPs[aGP->key()] = aGP;

  for (auto &itRegion : myResults) {
    OMTFResult aResult;
    aResult.configure(myOmtfConfig);
    itRegion[aGP->key()] = aResult;
  }

  return true;
}

averagePatterns()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::averagePatterns

(

int

charge

)

Average patterns. Use same meanDistPhi for two patterns neighboring in pt code. Averaging is made saparately fo each charge

Previously pt codes were going by steps of 1, now this is not the case

Definition at line 113 of file OMTFProcessor.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, GoldenPattern::getMeanDistPhi(), ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), OMTFProcessor< GoldenPatternType >::myOmtfConfig, OMTFConfiguration::nLayers(), OMTFConfiguration::nRefLayers(), GoldenPattern::setMeanDistPhi(), OMTFProcessor< GoldenPatternType >::shiftGP(), OMTFProcessor< GoldenPatternType >::theGPs, and Key::thePtCode.

Referenced by OMTFPatternMaker::endJob().

                                              {
  Key aKey(0, 9, charge);

  while (theGPs.find(aKey) != theGPs.end()) {
    GoldenPattern *aGP1 = theGPs.find(aKey)->second;
    GoldenPattern *aGP2 = aGP1;
    GoldenPattern *aGP3 = aGP1;
    GoldenPattern *aGP4 = aGP1;

    ++aKey.thePtCode;
    while (theGPs.find(aKey) == theGPs.end() && aKey.thePtCode <= 401)
      ++aKey.thePtCode;
    if (aKey.thePtCode <= 401 && theGPs.find(aKey) != theGPs.end())
      aGP2 = theGPs.find(aKey)->second;

    if (aKey.thePtCode > 71) {
      ++aKey.thePtCode;
      while (theGPs.find(aKey) == theGPs.end() && aKey.thePtCode <= 401)
        ++aKey.thePtCode;
      if (aKey.thePtCode <= 401 && theGPs.find(aKey) != theGPs.end())
        aGP3 = theGPs.find(aKey)->second;

      ++aKey.thePtCode;
      while (theGPs.find(aKey) == theGPs.end() && aKey.thePtCode <= 401)
        ++aKey.thePtCode;
      if (aKey.thePtCode <= 401 && theGPs.find(aKey) != theGPs.end())
        aGP4 = theGPs.find(aKey)->second;
    } else {
      aGP3 = aGP1;
      aGP4 = aGP2;
    }
    //HACK. Have to clean this up.
    ++aKey.thePtCode;
    while (theGPs.find(aKey) == theGPs.end() && aKey.thePtCode <= 401)
      ++aKey.thePtCode;

    GoldenPattern::vector2D meanDistPhi = aGP1->getMeanDistPhi();

    GoldenPattern::vector2D meanDistPhi1 = aGP1->getMeanDistPhi();
    GoldenPattern::vector2D meanDistPhi2 = aGP2->getMeanDistPhi();
    GoldenPattern::vector2D meanDistPhi3 = aGP3->getMeanDistPhi();
    GoldenPattern::vector2D meanDistPhi4 = aGP4->getMeanDistPhi();

    for (unsigned int iLayer = 0; iLayer < myOmtfConfig->nLayers(); ++iLayer) {
      for (unsigned int iRefLayer = 0; iRefLayer < myOmtfConfig->nRefLayers(); ++iRefLayer) {
        meanDistPhi[iLayer][iRefLayer] += meanDistPhi2[iLayer][iRefLayer];
        meanDistPhi[iLayer][iRefLayer] += meanDistPhi3[iLayer][iRefLayer];
        meanDistPhi[iLayer][iRefLayer] += meanDistPhi4[iLayer][iRefLayer];
        meanDistPhi[iLayer][iRefLayer] /= 4;
      }
    }

    aGP1->setMeanDistPhi(meanDistPhi);
    aGP2->setMeanDistPhi(meanDistPhi);

    shiftGP(aGP1, meanDistPhi, meanDistPhi1);
    shiftGP(aGP2, meanDistPhi, meanDistPhi2);
    if (aGP3 != aGP1 && aGP4 != aGP2) {
      aGP3->setMeanDistPhi(meanDistPhi);
      aGP4->setMeanDistPhi(meanDistPhi);
      shiftGP(aGP3, meanDistPhi, meanDistPhi3);
      shiftGP(aGP4, meanDistPhi, meanDistPhi4);
    }
  }
}

checkHitPatternValidity()

template<class GoldenPatternType >

bool OMTFProcessor< GoldenPatternType >::checkHitPatternValidity ( unsigned int  hits )

overrideprivatevirtual

Check if the hit pattern of given OMTF candite is not on the list of invalid hit patterns. Invalid hit patterns provode very little to efficiency, but gives high contribution to rate. Candidate with invalid hit patterns is assigned quality=0. Currently the list of invalid patterns is hardcoded. This has to be read from configuration.

FIXME: read the list from configuration so this can be controlled at runtime.

Implements IProcessorEmulator.

Definition at line 310 of file OMTFProcessor.cc.

References hfClusterShapes_cfi::hits.

                                                                                {
  std::vector<unsigned int> badPatterns = {
      99840, 34304, 3075, 36928, 12300, 98816, 98944, 33408, 66688, 66176, 7171, 20528, 33856, 35840, 4156, 34880};

  /*
99840 01100001 1000 000000      011000011000000000
34304 00100001 1000 000000      001000011000000000
 3075 00000011 0000 000011      000000110000000011
36928 00100100 0001 000000      001001000001000000
12300 00001100 0000 001100      000011000000001100
98816 01100000 1000 000000      011000001000000000
98944 01100000 1010 000000      011000001010000000
33408 00100000 1010 000000      001000001010000000
66688 01000001 0010 000000      010000010010000000
66176 01000000 1010 000000      010000001010000000
 7171 00000111 0000 000011      000001110000000011
20528 00010100 0000 110000      000101000000110000
33856 00100001 0001 000000      001000010001000000
35840 00100011 0000 000000      001000110000000000
 4156 00000100 0000 111100      000001000000111100
34880 00100010 0001 000000      001000100001000000
   */
  for (auto aHitPattern : badPatterns) {
    if (hits == aHitPattern)
      return false;
  }

  return true;
}

configure()

template<class GoldenPatternType >

bool OMTFProcessor< GoldenPatternType >::configure	(	const OMTFConfiguration *	omtfParams,
		const L1TMuonOverlapParams *	omtfPatterns
	)

Fill GP map with patterns from CondFormats object.

Mean dist phi data

Pdf data

Definition at line 35 of file OMTFProcessor.cc.

References OMTFProcessor< GoldenPatternType >::addGP(), L1TMuonOverlapParams::chargeLUT(), l1t::LUT::data(), Phase2L1GMT::etaLUT, L1TMuonOverlapParams::etaLUT(), l1tPhase2CaloJetEmulator_cfi::iEta, createfilelist::int, L1TMuonOverlapParams::meanDistPhiLUT(), OMTFProcessor< GoldenPatternType >::myOmtfConfig, OMTFProcessor< GoldenPatternType >::myResults, OMTFConfiguration::nGoldenPatterns(), OMTFConfiguration::nLayers(), OMTFConfiguration::nPdfAddrBits(), l1t::LUT::nrBitsData(), OMTFConfiguration::nRefLayers(), OMTFConfiguration::nTestRefHits(), L1TMuonOverlapParams::pdfLUT(), Phase2L1GMT::ptLUT, L1TMuonOverlapParams::ptLUT(), OMTFProcessor< GoldenPatternType >::resetConfiguration(), GoldenPattern::setMeanDistPhi(), and GoldenPattern::setPdf().

Referenced by OMTFPatternMaker::beginRun(), and OMTFReconstruction::beginRun().

                                                                                                           {
  resetConfiguration();

  myOmtfConfig = omtfConfig;

  myResults.assign(myOmtfConfig->nTestRefHits(), OMTFProcessor::resultsMap());

  const l1t::LUT *chargeLUT = omtfPatterns->chargeLUT();
  const l1t::LUT *etaLUT = omtfPatterns->etaLUT();
  const l1t::LUT *ptLUT = omtfPatterns->ptLUT();
  const l1t::LUT *pdfLUT = omtfPatterns->pdfLUT();
  const l1t::LUT *meanDistPhiLUT = omtfPatterns->meanDistPhiLUT();

  unsigned int nGPs = myOmtfConfig->nGoldenPatterns();
  unsigned int address = 0;
  unsigned int iEta, iPt;
  int iCharge;
  for (unsigned int iGP = 0; iGP < nGPs; ++iGP) {
    address = iGP;
    iEta = etaLUT->data(address);
    iCharge = chargeLUT->data(address) == 0 ? -1 : 1;
    iPt = ptLUT->data(address);

    GoldenPattern::vector2D meanDistPhi2D(myOmtfConfig->nLayers());
    GoldenPattern::vector1D pdf1D(exp2(myOmtfConfig->nPdfAddrBits()));
    GoldenPattern::vector3D pdf3D(myOmtfConfig->nLayers());
    GoldenPattern::vector2D pdf2D(myOmtfConfig->nRefLayers());
    for (unsigned int iLayer = 0; iLayer < myOmtfConfig->nLayers(); ++iLayer) {
      GoldenPattern::vector1D meanDistPhi1D(myOmtfConfig->nRefLayers());
      for (unsigned int iRefLayer = 0; iRefLayer < myOmtfConfig->nRefLayers(); ++iRefLayer) {
        address = iRefLayer + iLayer * myOmtfConfig->nRefLayers() +
                  iGP * (myOmtfConfig->nRefLayers() * myOmtfConfig->nLayers());
        meanDistPhi1D[iRefLayer] = meanDistPhiLUT->data(address) - (1 << (meanDistPhiLUT->nrBitsData() - 1));
      }
      meanDistPhi2D[iLayer] = meanDistPhi1D;
      for (unsigned int iRefLayer = 0; iRefLayer < myOmtfConfig->nRefLayers(); ++iRefLayer) {
        pdf1D.assign(1 << myOmtfConfig->nPdfAddrBits(), 0);
        for (unsigned int iPdf = 0; iPdf < (unsigned int)(1 << myOmtfConfig->nPdfAddrBits()); ++iPdf) {
          address = iPdf + iRefLayer * (1 << myOmtfConfig->nPdfAddrBits()) +
                    iLayer * myOmtfConfig->nRefLayers() * (1 << myOmtfConfig->nPdfAddrBits()) +
                    iGP * myOmtfConfig->nLayers() * myOmtfConfig->nRefLayers() * (1 << myOmtfConfig->nPdfAddrBits());
          pdf1D[iPdf] = pdfLUT->data(address);
        }
        pdf2D[iRefLayer] = pdf1D;
      }
      pdf3D[iLayer] = pdf2D;
    }
    Key aKey(iEta, iPt, iCharge, iGP);

    GoldenPattern *aGP = new GoldenPattern(aKey, myOmtfConfig);
    aGP->setMeanDistPhi(meanDistPhi2D);
    aGP->setPdf(pdf3D);
    addGP(aGP);
  }
  return true;
}

extrapolateDtPhiB()

template<class GoldenPatternType >

int OMTFProcessor< GoldenPatternType >::extrapolateDtPhiB	(	const MuonStubPtr &	refStub,
		const MuonStubPtr &	targetStub,
		unsigned int	targetLayer,
		const OMTFConfiguration *	omtfConfig
	)

Definition at line 564 of file OMTFProcessor.cc.

References OMTFProcessor< GoldenPatternType >::extrapolateDtPhiBFixedPoint(), OMTFProcessor< GoldenPatternType >::extrapolateDtPhiBFloatPoint(), and simOmtfPhase2Digis_cfi::useFloatingPointExtrapolation.

                                                                                             {
  if (useFloatingPointExtrapolation)
    return OMTFProcessor<GoldenPatternType>::extrapolateDtPhiBFloatPoint(refStub->logicLayer,
                                                                         refStub->phiHw,
                                                                         refStub->phiBHw,
                                                                         targetLayer,
                                                                         targetStub->phiHw,
                                                                         targetStub->qualityHw,
                                                                         targetStub->etaHw,
                                                                         targetStub->r,
                                                                         omtfConfig);
  return OMTFProcessor<GoldenPatternType>::extrapolateDtPhiBFixedPoint(refStub->logicLayer,
                                                                       refStub->phiHw,
                                                                       refStub->phiBHw,
                                                                       targetLayer,
                                                                       targetStub->phiHw,
                                                                       targetStub->qualityHw,
                                                                       targetStub->etaHw,
                                                                       targetStub->r,
                                                                       omtfConfig);
}

extrapolateDtPhiBFixedPoint()

template<class GoldenPatternType >

int OMTFProcessor< GoldenPatternType >::extrapolateDtPhiBFixedPoint	(	const int &	refLogicLayer,
		const int &	refPhi,
		const int &	refPhiB,
		unsigned int	targetLayer,
		const int &	targetStubPhi,
		const int &	targetStubQuality,
		const int &	targetStubEta,
		const int &	targetStubR,
		const OMTFConfiguration *	omtfConfig
	)

Definition at line 505 of file OMTFProcessor.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, LogTrace, scaleFactor, simOmtfPhase2Digis_cfi::useEndcapStubsRInExtr, and simOmtfPhase2Digis_cfi::useStubQualInExtr.

Referenced by OMTFProcessor< GoldenPatternType >::extrapolateDtPhiB().

                                                                                                       {
  int phiExtr = 0;  //delta phi extrapolated

  int reflLayerIndex = refLogicLayer == 0 ? 0 : 1;
  int extrFactor = 0;

  if (targetLayer == 0 || targetLayer == 2 || targetLayer == 4) {
    if (useStubQualInExtr)
      extrFactor = extrapolFactors[reflLayerIndex][targetLayer][targetStubQuality];
    else
      extrFactor = extrapolFactors[reflLayerIndex][targetLayer][0];
  } else if (targetLayer == 1 || targetLayer == 3 || targetLayer == 5) {
    int deltaPhi = targetStubPhi - refPhi;  //here targetStubPhi is phi, not phiB

    int scaleFactor = this->myOmtfConfig->omtfPhiUnit() * this->myOmtfConfig->dtPhiBUnitsRad() * 512;
    //= 305 for phase-1, 512 is multiplier so that scaleFactor is non-zero integer

    deltaPhi = (deltaPhi * scaleFactor) / 512;  //here deltaPhi is converted to the phi_b hw scale

    phiExtr = refPhiB - deltaPhi;  //phiExtr is also in phi_b hw scale
    //LogTrace("l1tOmtfEventPrint") <<__FUNCTION__<<":"<<__LINE__<<" deltaPhi "<<deltaPhi<<" phiExtr "<<phiExtr<<std::endl;

  } else if (targetLayer >= 10 && targetLayer <= 14) {
    extrFactor = extrapolFactors[reflLayerIndex][targetLayer][0];
  } else if ((targetLayer >= 6 && targetLayer <= 9) || (targetLayer >= 15 && targetLayer <= 17)) {
    if (useEndcapStubsRInExtr) {
      //if given abs(targetStubEta) value is not present in the map, it is added with default value of 0
      //so it should be good. The only problem is that the map can grow...
      //TODO change to targetStubR when it is implemented in the FW
      extrFactor = extrapolFactors[reflLayerIndex][targetLayer][abs(targetStubEta)];
      //extrFactor = extrapolFactors[reflLayerIndex][targetLayer][abs(targetStubR)];
    } else {
      extrFactor = extrapolFactors[reflLayerIndex][targetLayer][0];
    }
  }

  if (this->myOmtfConfig->isBendingLayer(targetLayer) == false) {
    phiExtr = extrFactor * refPhiB / extrapolMultiplier;
  }

  LogTrace("l1tOmtfEventPrint") << "\n"
                                << __FUNCTION__ << ":" << __LINE__ << " refLogicLayer " << refLogicLayer
                                << " targetLayer " << targetLayer << std::endl;
  LogTrace("l1tOmtfEventPrint") << "refPhi " << refPhi << " refPhiB " << refPhiB << " targetStubPhi " << targetStubPhi
                                << " targetStubQuality " << targetStubQuality << " targetStubEta " << targetStubEta
                                << " extrFactor " << extrFactor << " phiExtr " << phiExtr << std::endl;

  return phiExtr;
}

extrapolateDtPhiBFloatPoint()

template<class GoldenPatternType >

int OMTFProcessor< GoldenPatternType >::extrapolateDtPhiBFloatPoint	(	const int &	refLogicLayer,
		const int &	refPhi,
		const int &	refPhiB,
		unsigned int	targetLayer,
		const int &	targetStubPhi,
		const int &	targetStubQuality,
		const int &	targetStubEta,
		const int &	targetStubR,
		const OMTFConfiguration *	omtfConfig
	)

Fill GP vec with patterns from CondFormats object.

Definition at line 349 of file OMTFProcessor.cc.

References funct::abs(), ztail::d, SiPixelRawToDigiRegional_cfi::deltaPhi, ProcConfigurationBase::dtPhiBUnitsRad(), nano_mu_digi_cff::float, LogTrace, OMTFConfiguration::omtfPhiUnit(), funct::tan(), simOmtfPhase2Digis_cfi::useEndcapStubsRInExtr, and simOmtfPhase2Digis_cfi::useStubQualInExtr.

Referenced by OMTFProcessor< GoldenPatternType >::extrapolateDtPhiB().

                                                                                                       {
  LogTrace("l1tOmtfEventPrint") << "\n"
                                << __FUNCTION__ << ":" << __LINE__ << " refLogicLayer " << refLogicLayer
                                << " targetLayer " << targetLayer << std::endl;
  LogTrace("l1tOmtfEventPrint") << "refPhi " << refPhi << " refPhiB " << refPhiB << " targetStubPhi " << targetStubPhi
                                << " targetStubQuality " << targetStubQuality << std::endl;

  int phiExtr = 0;  //delta phi extrapolated

  float rRefLayer = 431.133;  //[cm], MB1 i.e. refLogicLayer = 0
  if (refLogicLayer == 2)
    rRefLayer = 512.401;  //MB2
  else if (refLogicLayer != 0) {
    return 0;
    //throw cms::Exception("OMTFProcessor<GoldenPatternType>::extrapolateDtPhiB: wrong refStubLogicLayer " + std::to_string(refLogicLayer) );
  }

  int reflLayerIndex = refLogicLayer == 0 ? 0 : 1;

  if (targetLayer == 0 || targetLayer == 2 || targetLayer == 4 || (targetLayer >= 10 && targetLayer <= 14)) {
    //all units are cm. Values from the CMS geometry
    float rTargetLayer = 512.401;  //MB2

    if (targetLayer == 0)
      rTargetLayer = 431.133;  //MB1
    else if (targetLayer == 4)
      rTargetLayer = 617.946;  //MB3

    else if (targetLayer == 10)
      rTargetLayer = 413.675;  //RB1in
    else if (targetLayer == 11)
      rTargetLayer = 448.675;  //RB1out
    else if (targetLayer == 12)
      rTargetLayer = 494.975;  //RB2in
    else if (targetLayer == 13)
      rTargetLayer = 529.975;  //RB2out
    else if (targetLayer == 14)
      rTargetLayer = 602.150;  //RB3

    if (useStubQualInExtr) {
      if (targetLayer == 0 || targetLayer == 2 || targetLayer == 4) {
        if (targetStubQuality == 2 || targetStubQuality == 0)
          rTargetLayer = rTargetLayer - 23.5 / 2;  //inner superlayer
        else if (targetStubQuality == 3 || targetStubQuality == 1)
          rTargetLayer = rTargetLayer + 23.5 / 2;  //outer superlayer
      }
    }

    float d = rTargetLayer - rRefLayer;
    //formula in the form as in the slides explaining the extrapolation algorithm
    //float deltaPhiExtr = d/rTargetLayer * refPhiB / omtfConfig->dtPhiBUnitsRad(); //[rad]
    //phiExtr = round(deltaPhiExtr / omtfConfig->omtfPhiUnit()); //[halfStrip]

    //formula with approximation, used to calculate extrFactor
    float extrFactor = d / rTargetLayer / omtfConfig->dtPhiBUnitsRad() / omtfConfig->omtfPhiUnit();
    phiExtr = extrFactor * (float)refPhiB;  //[halfStrip]

    //formula without approximation
    float deltaPhiExtr = atan(d / rTargetLayer * tan(refPhiB / omtfConfig->dtPhiBUnitsRad()));  //[rad]
    phiExtr = round(deltaPhiExtr / omtfConfig->omtfPhiUnit());                                  //[halfStrip]

    if (useStubQualInExtr & (targetLayer == 0 || targetLayer == 2 || targetLayer == 4)) {
      extrapolFactors[reflLayerIndex][targetLayer][targetStubQuality] = extrFactor;
      extrapolFactorsNorm[reflLayerIndex][targetLayer][targetStubQuality] = 1;
    } else {
      extrapolFactors[reflLayerIndex][targetLayer][0] = extrFactor;
      extrapolFactorsNorm[reflLayerIndex][targetLayer][0] = 1;
    }

    //LogTrace("l1tOmtfEventPrint") <<__FUNCTION__<<":"<<__LINE__<<" deltaPhiExtr "<<deltaPhiExtr<<" phiExtr "<<phiExtr<<std::endl;

    LogTrace("l1tOmtfEventPrint") << "\n"
                                  << __FUNCTION__ << ":" << __LINE__ << " refLogicLayer " << refLogicLayer
                                  << " targetLayer " << std::setw(2) << targetLayer << " targetStubQuality "
                                  << targetStubQuality << " extrFactor " << extrFactor << std::endl;

    LogTrace("l1tOmtfEventPrint") << __FUNCTION__ << ":" << __LINE__ << " refPhiB " << refPhiB << " phiExtr " << phiExtr
                                  << std::endl;

  } else if (targetLayer == 1 || targetLayer == 3 || targetLayer == 5) {
    int deltaPhi = targetStubPhi - refPhi;  //[halfStrip]

    //deltaPhi is here in phi_b hw scale
    deltaPhi = round(deltaPhi * omtfConfig->omtfPhiUnit() * omtfConfig->dtPhiBUnitsRad());

    phiExtr = refPhiB - deltaPhi;  //phiExtr is also in phi_b hw scale
    LogTrace("l1tOmtfEventPrint") << __FUNCTION__ << ":" << __LINE__ << " deltaPhi " << deltaPhi << " phiExtr "
                                  << phiExtr << std::endl;
  } else if ((targetLayer >= 6 && targetLayer <= 9) || (targetLayer >= 15 && targetLayer <= 17)) {
    //if true, for the CSC and endcap RPC the R is taken from the hit coordinates

    float rME = targetStubR;
    if (!useEndcapStubsRInExtr) {
      //all units are cm. This are the average R values for a given chamber (more or less middle of the chamber, but taking into account the OMTF eta range)
      if (targetLayer == 6 || targetLayer == 15)  //ME1/3, RE1/3,
        rME = 600.;
      else if (targetLayer == 7 || targetLayer == 15) {  //ME2/2, RE2/3,
        if (refLogicLayer == 0)
          rME = 600.;
        else
          rME = 640.;
      } else if (targetLayer == 8 || rME == 16) {  //ME3/2, RE3/3,
        if (refLogicLayer == 0)
          rME = 620.;
        else
          rME = 680.;
      } else if (targetLayer == 9) {
        rME = 460.;  //for the refLogicLayer = 1. refLogicLayer = 2 is impossible
      }
    }

    float d = rME - rRefLayer;
    //formula in the form as in the slides explaining the extrapolation algorithm
    //float deltaPhiExtr = d / rME * refPhiB / omtfConfig->dtPhiBUnitsRad();  //[rad]
    //phiExtr = round(deltaPhiExtr / omtfConfig->omtfPhiUnit()); //[halfStrip]

    //formula with approximation, used to calculate extrFactor
    float extrFactor = d / rME / omtfConfig->dtPhiBUnitsRad() / omtfConfig->omtfPhiUnit();
    phiExtr = extrFactor * refPhiB;  //[halfStrip]

    //formula without approximation
    float deltaPhiExtr = atan(d / rME * tan(refPhiB / omtfConfig->dtPhiBUnitsRad()));  //[rad]
    phiExtr = round(deltaPhiExtr / omtfConfig->omtfPhiUnit());                         //[halfStrip]

    if (useEndcapStubsRInExtr) {
      //extrapolFactors[reflLayerIndex][targetLayer][std::abs(targetStubEta)] += extrFactor;
      //extrapolFactorsNorm[reflLayerIndex][targetLayer][std::abs(targetStubEta)]++;
      extrapolFactors[reflLayerIndex][targetLayer][std::abs(rME)] += extrFactor;
      extrapolFactorsNorm[reflLayerIndex][targetLayer][std::abs(rME)]++;
      //extrapolFactors[reflLayerIndex][targetLayer][0] += extrFactor;
      //extrapolFactorsNorm[reflLayerIndex][targetLayer][0]++;
    } else {
      extrapolFactors[reflLayerIndex][targetLayer][0] = extrFactor;
      extrapolFactorsNorm[reflLayerIndex][targetLayer][0] = 1;
    }
    LogTrace("l1tOmtfEventPrint") << "\n"
                                  << __FUNCTION__ << ":" << __LINE__ << " refLogicLayer " << refLogicLayer
                                  << " targetLayer " << std::setw(2) << targetLayer << " targetStubR " << targetStubR
                                  << " targetStubEta " << targetStubEta << " extrFactor "
                                  << " rRefLayer " << rRefLayer << " d " << d << " deltaPhiExtr " << deltaPhiExtr
                                  << " phiExtr " << phiExtr << std::endl;
  }
  //TODO restrict the range of the phiExtr and refPhiB !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  return phiExtr;
}

fillCounts()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::fillCounts	(	unsigned int	iProcessor,
		const OMTFinput &	aInput,
		const SimTrack *	aSimMuon
	)

Fill counts for a GoldenPattern of this processor unit. Pattern key is selcted according to the SimTrack parameters.

Stupid conersion. Have to go through PAC pt scale, as we later shift resulting pt code by +1
Number of reference hits to be checked.

Definition at line 286 of file OMTFProcessor.cc.

References funct::abs(), OMTFConfiguration::getBendingLayers(), OMTFinput::getLayerData(), OMTFinput::getRefHits(), OMTFConfiguration::getRefHitsDefs(), OMTFConfiguration::getRefToLogicNumber(), ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), RefHitDef::iInput, RPCConst::iptFromPt(), RefHitDef::iRefLayer, RefHitDef::iRegion, CoreSimTrack::momentum(), OMTFProcessor< GoldenPatternType >::myOmtfConfig, OMTFConfiguration::nLayers(), OMTFConfiguration::nRefHits(), RPCConst::ptFromIpt(), OMTFProcessor< GoldenPatternType >::restrictInput(), AlignmentTrackSelector_cfi::theCharge, OMTFProcessor< GoldenPatternType >::theGPs, and CoreSimTrack::type().

Referenced by OMTFPatternMaker::analyze().

                                                                                                         {
  int theCharge = (abs(aSimMuon->type()) == 13) ? aSimMuon->type() / -13 : 0;
  unsigned int iPt = RPCConst::iptFromPt(aSimMuon->momentum().pt());
  iPt += 1;
  if (iPt > 31)
    iPt = 200 * 2 + 1;
  else
    iPt = RPCConst::ptFromIpt(iPt) * 2.0 +
          1;  //MicroGMT has 0.5 GeV step size, with lower bin edge  (uGMT_pt_code - 1)*step_size

  std::bitset<128> refHitsBits = aInput.getRefHits(iProcessor);
  if (refHitsBits.none())
    return;

  std::ostringstream myStr;
  myStr << "iProcessor: " << iProcessor << std::endl;
  myStr << "Input: ------------" << std::endl;
  myStr << aInput << std::endl;
  edm::LogInfo("OMTF processor") << myStr.str();

  for (unsigned int iLayer = 0; iLayer < myOmtfConfig->nLayers(); ++iLayer) {
    const OMTFinput::vector1D &layerHits = aInput.getLayerData(iLayer);
    if (layerHits.empty())
      continue;
    for (unsigned int iRefHit = 0; iRefHit < myOmtfConfig->nRefHits(); ++iRefHit) {
      if (!refHitsBits[iRefHit])
        continue;
      const RefHitDef &aRefHitDef = myOmtfConfig->getRefHitsDefs()[iProcessor][iRefHit];
      int phiRef = aInput.getLayerData(myOmtfConfig->getRefToLogicNumber()[aRefHitDef.iRefLayer])[aRefHitDef.iInput];
      unsigned int iRegion = aRefHitDef.iRegion;
      if (myOmtfConfig->getBendingLayers().count(iLayer))
        phiRef = 0;
      const OMTFinput::vector1D restrictedLayerHits = restrictInput(iProcessor, iRegion, iLayer, layerHits);
      for (auto itGP : theGPs) {
        if (itGP.first.theCharge != theCharge)
          continue;
        if (itGP.first.thePtCode != iPt)
          continue;
        itGP.second->addCount(aRefHitDef.iRefLayer, iLayer, phiRef, restrictedLayerHits);
      }
    }
  }
}

fillInputRange() [1/2]

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::fillInputRange ( unsigned int  iProcessor, unsigned int  iCone, const OMTFinput &  aInput  )

private

Fill map of used inputs. FIXME: using hack from OMTFConfiguration

fillInputRange() [2/2]

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::fillInputRange ( unsigned int  iProcessor, unsigned int  iCone, unsigned int  iRefLayer, unsigned int  iHit  )

private

getFinalcandidates()

template<class GoldenPatternType >

std::vector< l1t::RegionalMuonCand > OMTFProcessor< GoldenPatternType >::getFinalcandidates ( unsigned int  iProcessor, l1t::tftype  mtfType, const AlgoMuons &  algoCands  )

overridevirtual

Implements IProcessorEmulator.

Definition at line 101 of file OMTFProcessor.cc.

References funct::abs(), l1t::RegionalMuonCand::hwPt(), l1t::RegionalMuonCand::hwPtUnconstrained(), l1t::omtf_pos, quality, mps_fire::result, l1t::RegionalMuonCand::setHwEta(), l1t::RegionalMuonCand::setHwPhi(), l1t::RegionalMuonCand::setHwPt(), l1t::RegionalMuonCand::setHwPtUnconstrained(), l1t::RegionalMuonCand::setHwQual(), l1t::RegionalMuonCand::setHwSign(), l1t::RegionalMuonCand::setHwSignValid(), l1t::RegionalMuonCand::setTFIdentifiers(), and l1t::RegionalMuonCand::setTrackAddress().

                                                                                                                    {
  std::vector<l1t::RegionalMuonCand> result;

  for (auto& myCand : algoCands) {
    l1t::RegionalMuonCand candidate;

    //the charge is only for the constrained measurement. The constrained measurement is always defined for a valid candidate
    if (ptAssignment) {
      if (myCand->getPdfSumConstr() > 0 && myCand->getFiredLayerCntConstr() >= 3)
        candidate.setHwPt(myCand->getPtNNConstr());
      else if (myCand->getPtUnconstr() > 0)
        candidate.setHwPt(1);
      else
        candidate.setHwPt(0);

      candidate.setHwSign(myCand->getChargeNNConstr() < 0 ? 1 : 0);
    } else {
      if (myCand->getPdfSumConstr() > 0 && myCand->getFiredLayerCntConstr() >= 3)
        candidate.setHwPt(myCand->getPtConstr());
      else if (myCand->getPtUnconstr() > 0)
        //if myCand->getPdfSumConstr() == 0, the myCand->getPtConstr() might not be 0, see the end of GhostBusterPreferRefDt::select
        //but 0 means empty candidate, 1 means pt=0, therefore here we set HwPt to 1, as the PtUnconstr > 0
        candidate.setHwPt(1);
      else
        candidate.setHwPt(0);

      candidate.setHwSign(myCand->getChargeConstr() < 0 ? 1 : 0);
    }

    if (mtfType == l1t::omtf_pos)
      candidate.setHwEta(myCand->getEtaHw());
    else
      candidate.setHwEta((-1) * myCand->getEtaHw());

    int phiValue = myCand->getPhi();
    if (phiValue >= int(this->myOmtfConfig->nPhiBins()))
      phiValue -= this->myOmtfConfig->nPhiBins();
    phiValue = this->myOmtfConfig->procPhiToGmtPhi(phiValue);
    candidate.setHwPhi(phiValue);

    candidate.setHwSignValid(1);

    if (myCand->getPtUnconstr() >= 0) {  //empty PtUnconstrained is -1, maybe should be corrected on the source
      //the upt has different hardware scale than the pt, the upt unit is 1 GeV
      candidate.setHwPtUnconstrained(myCand->getPtUnconstr());
    } else
      candidate.setHwPtUnconstrained(0);

    unsigned int quality = 12;
    if (this->myOmtfConfig->fwVersion() <= 6)
      quality = checkHitPatternValidity(myCand->getFiredLayerBits()) ? 0 | (1 << 2) | (1 << 3) : 0 | (1 << 2);  //12 : 4

    if (abs(myCand->getEtaHw()) == 115 &&  //115 is eta 1.25                        rrrrrrrrccccdddddd
        (static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000001110000000").to_ulong() ||
         static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000001110000000").to_ulong() ||
         static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000000110000000").to_ulong() ||
         static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000001100000000").to_ulong() ||
         static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000001010000000").to_ulong())) {
      if (this->myOmtfConfig->fwVersion() <= 6)
        quality = 4;
      else
        quality = 1;
    }

    if (this->myOmtfConfig->fwVersion() >= 5 && this->myOmtfConfig->fwVersion() <= 6) {
      if (static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010000000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000100000000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000000000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000000000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000000000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000100000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010000000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000100000000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000000000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000000000110000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000000000110000").to_ulong())
        quality = 1;
    } else if (this->myOmtfConfig->fwVersion() >= 8) {  //TODO fix the fwVersion     rrrrrrrrccccdddddd
      if (static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000110000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000110000000001").to_ulong() ||

          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000011000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000011000000000100").to_ulong() ||

          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000011000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000010000000001").to_ulong())
        quality = 1;
      else if (
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010001000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000011000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000011000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000011100000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100001000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100100000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000110100000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000111000000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000111000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000111000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000001000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001010000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001010000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001010000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100000000111").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100001000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001110000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001110000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010000000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010010000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010010000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010010000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010100000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010100000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000011110000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000011110000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000101000000010101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000010000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000011000000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000011000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000100000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000110000000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001001000000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001001100000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001010000000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000000010000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000000011000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000010000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000100000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000011000000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000110000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000110000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000011000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000011000000000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000010010000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001001000001000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000100000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001110000000111").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000110001000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001110000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000000001000100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000110001000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000000000101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001010000001000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001100000001000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("100000010000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000010010000000").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000010100000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000110000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001000000001100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000000000000111101").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000001100000110001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000100000000010100").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000100000000011").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("001000110000000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("010000100010000001").to_ulong() ||
          static_cast<unsigned int>(myCand->getFiredLayerBits()) == std::bitset<18>("000100000000110000").to_ulong())
        quality = 8;
    }  //  if (abs(myCand->getEta()) == 121) quality = 4;
    if (abs(myCand->getEtaHw()) >= 121)
      quality = 0;  // changed from 4 on request from HI

    candidate.setHwQual(quality);

    std::map<int, int> trackAddr;
    trackAddr[0] = myCand->getFiredLayerBits();
    //TODO in the hardware, the uPt is sent to the uGMT at the trackAddr = (uPt << 18) + trackAddr;
    //check if it matters if it needs to be here as well
    trackAddr[1] = myCand->getRefLayer();
    trackAddr[2] = myCand->getDisc();
    if (candidate.hwPt() > 0 || candidate.hwPtUnconstrained() > 0) {
      candidate.setTrackAddress(trackAddr);
      candidate.setTFIdentifiers(iProcessor, mtfType);
      result.push_back(candidate);
    }
  }
  return result;
}

getPatterns()

template<class GoldenPatternType >

const std::map< Key, GoldenPattern * > & OMTFProcessor< GoldenPatternType >::getPatterns

(

)

const

Return map of GoldenPatterns.

Definition at line 206 of file OMTFProcessor.cc.

References OMTFProcessor< GoldenPatternType >::theGPs.

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

{ return theGPs; }

ghostBust()

template<class GoldenPatternType >

AlgoMuons OMTFProcessor< GoldenPatternType >::ghostBust ( AlgoMuons  refHitCands, int  charge = 0  )

inlineoverridevirtual

Implements IProcessorEmulator.

Definition at line 96 of file OMTFProcessor.h.

References ALCARECOTkAlJpsiMuMu_cff::charge, and OMTFProcessor< GoldenPatternType >::ghostBuster.

                                                                      {
    return ghostBuster->select(refHitCands, charge);
  }

init()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::init ( const edm::ParameterSet &  edmCfg, edm::EventSetup const &  evSetup  )

privatevirtual

Definition at line 59 of file OMTFProcessor.cc.

References edm::ParameterSet::exists(), simOmtfDigis_cfi::extrapolFactorsFilename, contentValuesFiles::fullPath, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, simOmtfPhase2Digis_cfi::useEndcapStubsRInExtr, simOmtfPhase2Digis_cfi::useFloatingPointExtrapolation, simOmtfPhase2Digis_cfi::useStubQualInExtr, and trackerHitRTTI::vector.

Referenced by OMTFProcessor< GoldenPatternType >::OMTFProcessor().

                                                                                                       {
  setSorter(new OMTFSorter<GoldenPatternType>(this->myOmtfConfig->getSorterType()));
  //initialize with the default sorter

  if (this->myOmtfConfig->getGhostBusterType() == "GhostBusterPreferRefDt" ||
      this->myOmtfConfig->getGhostBusterType() == "byLLH" || this->myOmtfConfig->getGhostBusterType() == "byFPLLH" ||
      this->myOmtfConfig->getGhostBusterType() == "byRefLayer") {
    setGhostBuster(new GhostBusterPreferRefDt(this->myOmtfConfig));
    edm::LogVerbatim("OMTFReconstruction") << "setting " << this->myOmtfConfig->getGhostBusterType() << std::endl;
  } else {
    setGhostBuster(new GhostBuster(this->myOmtfConfig));  //initialize with the default sorter
    edm::LogVerbatim("OMTFReconstruction") << "setting GhostBuster" << std::endl;
  }

  edm::LogVerbatim("OMTFReconstruction") << "fwVersion 0x" << hex << this->myOmtfConfig->fwVersion() << std::endl;

  useStubQualInExtr = this->myOmtfConfig->useStubQualInExtr();
  useEndcapStubsRInExtr = this->myOmtfConfig->useEndcapStubsRInExtr();

  if (edmCfg.exists("useFloatingPointExtrapolation"))
    useFloatingPointExtrapolation = edmCfg.getParameter<bool>("useFloatingPointExtrapolation");

  std::string extrapolFactorsFilename;
  if (edmCfg.exists("extrapolFactorsFilename")) {
    extrapolFactorsFilename = edmCfg.getParameter<edm::FileInPath>("extrapolFactorsFilename").fullPath();
  }

  if (this->myOmtfConfig->usePhiBExtrapolationMB1() || this->myOmtfConfig->usePhiBExtrapolationMB2()) {
    extrapolFactors.resize(2, std::vector<std::map<int, double> >(this->myOmtfConfig->nLayers()));
    extrapolFactorsNorm.resize(2, std::vector<std::map<int, int> >(this->myOmtfConfig->nLayers()));

    //when useFloatingPointExtrapolation is true the extrapolFactors are not used,
    //all calculations are done in the extrapolateDtPhiBFloatPoint
    if (!extrapolFactorsFilename.empty() && !useFloatingPointExtrapolation)
      loadExtrapolFactors(extrapolFactorsFilename);
  }

  edm::LogVerbatim("OMTFReconstruction") << "useFloatingPointExtrapolation " << useFloatingPointExtrapolation
                                         << std::endl;
}

loadExtrapolFactors()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::loadExtrapolFactors

(

const std::string &

filename

)

Definition at line 893 of file OMTFProcessor.cc.

References corrVsCorr::filename, submitPVResolutionJobs::key, LogTrace, AlCaHLTBitMon_QueryRunRegistry::string, compare::tree, and relativeConstraints::value.

                                                                                    {
  boost::property_tree::ptree tree;

  boost::property_tree::read_xml(filename, tree);

  edm::LogVerbatim("OMTFReconstruction") << "loadExtrapolFactors from file " << filename << std::endl;

  extrapolMultiplier = tree.get<int>("ExtrapolationFactors.<xmlattr>.multiplier");
  edm::LogVerbatim("OMTFReconstruction") << "extrapolMultiplier " << extrapolMultiplier << std::endl;

  auto& lutNodes = tree.get_child("ExtrapolationFactors");
  for (boost::property_tree::ptree::value_type& lutNode : lutNodes) {
    if (lutNode.first == "Lut") {
      int iRefLayer = lutNode.second.get<int>("<xmlattr>.RefLayer");
      int iLayer = lutNode.second.get<int>("<xmlattr>.Layer");
      std::string keyType = lutNode.second.get<std::string>("<xmlattr>.KeyType");

      LogTrace("OMTFReconstruction") << "iRefLayer " << iRefLayer << " iLayer " << iLayer << " keyType " << keyType
                                     << std::endl;

      auto& valueNodes = lutNode.second;
      for (boost::property_tree::ptree::value_type& valueNode : valueNodes) {
        if (valueNode.first == "LutVal") {
          int key = valueNode.second.get<int>("<xmlattr>.key");
          float value = valueNode.second.get<float>("<xmlattr>.value");
          extrapolFactors.at(iRefLayer).at(iLayer)[key] = value;
          LogTrace("OMTFReconstruction") << "key " << key << " value " << value << std::endl;
        }
      }
    }
  }
}

printInfo()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::printInfo ( ) const

overridevirtual

Implements IProcessorEmulator.

Definition at line 834 of file OMTFProcessor.cc.

References ProcessorBase< GoldenPatternType >::printInfo().

                                                       {
  edm::LogVerbatim("OMTFReconstruction") << __PRETTY_FUNCTION__ << std::endl;

  ProcessorBase<GoldenPatternType>::printInfo();
}

processInput() [1/2]

template<class GoldenPatternType >

const std::vector< OMTFProcessor::resultsMap > & OMTFProcessor< GoldenPatternType >::processInput	(	unsigned int	iProcessor,
		const OMTFinput &	aInput
	)

Process input data from a single event Input data is represented by hits in logic layers expressed in local coordinates Vector index: logic region number Map key: GoldenPattern key

Number of reference hits to be checked.

Definition at line 209 of file OMTFProcessor.cc.

References OMTFConfiguration::getBendingLayers(), OMTFinput::getLayerData(), OMTFinput::getRefHits(), OMTFConfiguration::getRefHitsDefs(), OMTFConfiguration::getRefToLogicNumber(), ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), RefHitDef::iInput, RefHitDef::iRefLayer, RefHitDef::iRegion, OMTFProcessor< GoldenPatternType >::myOmtfConfig, OMTFProcessor< GoldenPatternType >::myResults, OMTFConfiguration::nLayers(), OMTFConfiguration::nRefHits(), OMTFConfiguration::nTestRefHits(), OMTFProcessor< GoldenPatternType >::restrictInput(), and OMTFProcessor< GoldenPatternType >::theGPs.

Referenced by OMTFReconstruction::getProcessorCandidates().

                                                                                                   {
  for (auto &itRegion : myResults)
    for (auto &itKey : itRegion)
      itKey.second.clear();

  std::bitset<128> refHitsBits = aInput.getRefHits(iProcessor);
  if (refHitsBits.none())
    return myResults;

  for (unsigned int iLayer = 0; iLayer < myOmtfConfig->nLayers(); ++iLayer) {
    const OMTFinput::vector1D &layerHits = aInput.getLayerData(iLayer);
    if (layerHits.empty())
      continue;
    unsigned int nTestedRefHits = myOmtfConfig->nTestRefHits();
    for (unsigned int iRefHit = 0; iRefHit < myOmtfConfig->nRefHits(); ++iRefHit) {
      if (!refHitsBits[iRefHit])
        continue;
      if (nTestedRefHits-- == 0)
        break;
      const RefHitDef &aRefHitDef = myOmtfConfig->getRefHitsDefs()[iProcessor][iRefHit];

      int phiRef = aInput.getLayerData(myOmtfConfig->getRefToLogicNumber()[aRefHitDef.iRefLayer])[aRefHitDef.iInput];
      int etaRef =
          aInput.getLayerData(myOmtfConfig->getRefToLogicNumber()[aRefHitDef.iRefLayer], true)[aRefHitDef.iInput];
      unsigned int iRegion = aRefHitDef.iRegion;

      if (myOmtfConfig->getBendingLayers().count(iLayer))
        phiRef = 0;
      const OMTFinput::vector1D restrictedLayerHits = restrictInput(iProcessor, iRegion, iLayer, layerHits);
      for (auto itGP : theGPs) {
        GoldenPattern::layerResult aLayerResult =
            itGP.second->process1Layer1RefLayer(aRefHitDef.iRefLayer, iLayer, phiRef, restrictedLayerHits);
        int phiRefSt2 = itGP.second->propagateRefPhi(phiRef, etaRef, aRefHitDef.iRefLayer);
        myResults[myOmtfConfig->nTestRefHits() - nTestedRefHits - 1][itGP.second->key()].setRefPhiRHits(
            aRefHitDef.iRefLayer, phiRef);
        myResults[myOmtfConfig->nTestRefHits() - nTestedRefHits - 1][itGP.second->key()].addResult(
            aRefHitDef.iRefLayer, iLayer, aLayerResult.first, phiRefSt2, etaRef);
      }
    }
  }
  for (auto &itRefHit : myResults)
    for (auto &itKey : itRefHit)
      itKey.second.finalise();

  std::ostringstream myStr;
  myStr << "iProcessor: " << iProcessor << std::endl;
  myStr << "Input: ------------" << std::endl;
  myStr << aInput << std::endl;
  edm::LogInfo("OMTF processor") << myStr.str();

  return myResults;
}

processInput() [2/2]

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::processInput ( unsigned int  iProcessor, l1t::tftype  mtfType, const OMTFinput &  aInput, std::vector< std::unique_ptr< IOMTFEmulationObserver > > &  observers  )

overridevirtual

Process input data from a single event Input data is represented by hits in logic layers expressed in local coordinates

Implements IProcessorEmulator.

Definition at line 592 of file OMTFProcessor.cc.

References OMTFinput::getMuonStub(), OMTFinput::getRefHits(), RefHitDef::iInput, RefHitDef::iRefLayer, RefHitDef::iRegion, LogTrace, mps_fire::result, and simOmtfPhase2Digis_cfi::useFloatingPointExtrapolation.

                                                                                                                  {
  unsigned int procIndx = this->myOmtfConfig->getProcIndx(iProcessor, mtfType);
  for (auto& itGP : this->theGPs) {
    for (auto& result : itGP->getResults()[procIndx]) {
      result.reset();
    }
  }

  LogTrace("l1tOmtfEventPrint") << __FUNCTION__ << "\n"
                                << __LINE__ << " iProcessor " << iProcessor << " mtfType " << mtfType << " procIndx "
                                << procIndx << " ----------------------" << std::endl;
  std::vector<const RefHitDef*> refHitDefs;

  {
    auto refHitsBits = aInput.getRefHits(iProcessor);
    if (refHitsBits.none())
      return;  // myResults;

    //loop over all possible refHits, e.g. 128
    for (unsigned int iRefHit = 0; iRefHit < this->myOmtfConfig->nRefHits(); ++iRefHit) {
      if (!refHitsBits[iRefHit])
        continue;

      refHitDefs.push_back(&(this->myOmtfConfig->getRefHitsDefs()[iProcessor][iRefHit]));

      if (refHitDefs.size() == this->myOmtfConfig->nTestRefHits())
        break;
    }
  }

  boost::property_tree::ptree procDataTree;
  LogTrace("l1tOmtfEventPrint") << __FUNCTION__ << " " << __LINE__;
  for (unsigned int iLayer = 0; iLayer < this->myOmtfConfig->nLayers(); ++iLayer) {
    //debug
    /*for(auto& h : layerHits) {
      if(h != 5400)
        LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<" "<<__LINE__<<" iLayer "<<iLayer<<" layerHit "<<h<<std::endl;
    }*/

    for (unsigned int iRefHit = 0; iRefHit < refHitDefs.size(); iRefHit++) {
      const RefHitDef& aRefHitDef = *(refHitDefs[iRefHit]);

      unsigned int refLayerLogicNum = this->myOmtfConfig->getRefToLogicNumber()[aRefHitDef.iRefLayer];
      const MuonStubPtr refStub = aInput.getMuonStub(refLayerLogicNum, aRefHitDef.iInput);
      //int etaRef = refStub->etaHw;

      unsigned int iRegion = aRefHitDef.iRegion;

      MuonStubPtrs1D restrictedLayerStubs = this->restrictInput(iProcessor, iRegion, iLayer, aInput);

      //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<" "<<__LINE__<<" iLayer "<<iLayer<<" iRefLayer "<<aRefHitDef.iRefLayer<<std::endl;
      //LogTrace("l1tOmtfEventPrint")<<"iLayer "<<iLayer<<" iRefHit "<<iRefHit;
      //LogTrace("l1tOmtfEventPrint")<<" nTestedRefHits "<<nTestedRefHits<<" aRefHitDef "<<aRefHitDef<<std::endl;

      std::vector<int> extrapolatedPhi(restrictedLayerStubs.size(), 0);

      //TODO make sure the that the iRefLayer numbers used here corresponds to this in the hwToLogicLayer_0x000X.xml
      if ((this->myOmtfConfig->usePhiBExtrapolationMB1() && aRefHitDef.iRefLayer == 0) ||
          (this->myOmtfConfig->usePhiBExtrapolationMB2() && aRefHitDef.iRefLayer == 2)) {
        if ((iLayer != refLayerLogicNum) && (iLayer != refLayerLogicNum + 1)) {
          unsigned int iStub = 0;
          for (auto& targetStub : restrictedLayerStubs) {
            if (targetStub) {
              extrapolatedPhi[iStub] = extrapolateDtPhiB(refStub, targetStub, iLayer, this->myOmtfConfig);

              LogTrace("l1tOmtfEventPrint")
                  << "\n"
                  << __FUNCTION__ << ":" << __LINE__ << " extrapolating from layer " << refLayerLogicNum
                  << " - iRefLayer " << aRefHitDef.iRefLayer << " to layer " << iLayer << " stub " << targetStub
                  << " value " << extrapolatedPhi[iStub] << std::endl;

              if (this->myOmtfConfig->getDumpResultToXML()) {
                auto& extrapolatedPhiTree = procDataTree.add_child("extrapolatedPhi", boost::property_tree::ptree());
                extrapolatedPhiTree.add("<xmlattr>.refLayer", refLayerLogicNum);
                extrapolatedPhiTree.add("<xmlattr>.layer", iLayer);
                extrapolatedPhiTree.add("<xmlattr>.refPhiBHw", refStub->phiBHw);
                extrapolatedPhiTree.add("<xmlattr>.iStub", iStub);
                extrapolatedPhiTree.add("<xmlattr>.qualityHw", targetStub->qualityHw);
                extrapolatedPhiTree.add("<xmlattr>.etaHw", targetStub->etaHw);
                extrapolatedPhiTree.add("<xmlattr>.phiExtr", extrapolatedPhi[iStub]);

                if (this->myOmtfConfig->isBendingLayer(iLayer))
                  extrapolatedPhiTree.add("<xmlattr>.dist_phi", targetStub->phiBHw - extrapolatedPhi[iStub]);
                else
                  extrapolatedPhiTree.add("<xmlattr>.dist_phi", targetStub->phiHw - extrapolatedPhi[iStub]);
              }
            }
            iStub++;
          }
        }
      }

      for (auto& itGP : this->theGPs) {
        if (itGP->key().thePt == 0)  //empty pattern
          continue;

        StubResult stubResult =
            itGP->process1Layer1RefLayer(aRefHitDef.iRefLayer, iLayer, restrictedLayerStubs, extrapolatedPhi, refStub);

        /* LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__
                                     <<" layerResult: valid"<<stubResult.getValid()
                                     <<" pdfVal "<<stubResult.getPdfVal()
                                     <<std::endl;*/

        itGP->getResults()[procIndx][iRefHit].setStubResult(iLayer, stubResult);
      }
    }
  }

  for (unsigned int iRefHit = 0; iRefHit < refHitDefs.size(); iRefHit++) {
    const RefHitDef& aRefHitDef = *(refHitDefs[iRefHit]);

    unsigned int refLayerLogicNum = this->myOmtfConfig->getRefToLogicNumber()[aRefHitDef.iRefLayer];
    const MuonStubPtr refStub = aInput.getMuonStub(refLayerLogicNum, aRefHitDef.iInput);

    int phiRef = refStub->phiHw;
    int etaRef = refStub->etaHw;

    //calculating the phiExtrp in the case the RefLayer is MB1, to include it in the  candidate phi of candidate
    int phiExtrp = 0;
    if ((this->myOmtfConfig->usePhiBExtrapolationMB1() && aRefHitDef.iRefLayer == 0)) {
      //||(this->myOmtfConfig->getUsePhiBExtrapolationMB2() && aRefHitDef.iRefLayer == 2) ) {  //the extrapolation from the layer 2 to the layer 2 has no sense, so phiExtrp is 0
      LogTrace("l1tOmtfEventPrint") << "\n"
                                    << __FUNCTION__ << ":" << __LINE__
                                    << "extrapolating ref hit to get the phi of the candidate" << std::endl;
      if (useFloatingPointExtrapolation)
        phiExtrp = extrapolateDtPhiBFloatPoint(
            aRefHitDef.iRefLayer, phiRef, refStub->phiBHw, 2, 0, 6, 0, 0, this->myOmtfConfig);
      else
        phiExtrp = extrapolateDtPhiBFixedPoint(
            aRefHitDef.iRefLayer, phiRef, refStub->phiBHw, 2, 0, 6, 0, 0, this->myOmtfConfig);
    }

    for (auto& itGP : this->theGPs) {
      if (itGP->key().thePt == 0)  //empty pattern
        continue;

      int phiRefSt2 = itGP->propagateRefPhi(phiRef + phiExtrp, etaRef, aRefHitDef.iRefLayer);
      itGP->getResults()[procIndx][iRefHit].set(aRefHitDef.iRefLayer, phiRefSt2, etaRef, phiRef);
    }
  }

  {
    for (auto& itGP : this->theGPs) {
      itGP->finalise(procIndx);
      //debug
      /*for(unsigned int iRefHit = 0; iRefHit < itGP->getResults()[procIndx].size(); ++iRefHit) {
        if(itGP->getResults()[procIndx][iRefHit].isValid()) {
          LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<"__LINE__"<<itGP->getResults()[procIndx][iRefHit]<<std::endl;
        }
      }*/
    }
  }

  for (auto& obs : observers)
    obs->addProcesorData("extrapolation", procDataTree);

  return;
}

resetConfiguration()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::resetConfiguration ( )

privatevirtual

Reset all configuration parameters.

Reimplemented from ProcessorBase< GoldenPatternType >.

Definition at line 27 of file OMTFProcessor.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, OMTFProcessor< GoldenPatternType >::myResults, and OMTFProcessor< GoldenPatternType >::theGPs.

Referenced by OMTFProcessor< GoldenPatternType >::configure().

                                       {
  myResults.clear();
  for (auto it : theGPs)
    delete it.second;
  theGPs.clear();
}

restrictInput()

template<class GoldenPatternType >

OMTFinput::vector1D OMTFProcessor< GoldenPatternType >::restrictInput ( unsigned int  iProcessor, unsigned int  iCone, unsigned int  iLayer, const OMTFinput::vector1D &  layerHits  )

private

Remove hits whis are outside input range for given processor and cone

Definition at line 269 of file OMTFProcessor.cc.

References OMTFConfiguration::getConnections(), OMTFProcessor< GoldenPatternType >::myOmtfConfig, and OMTFConfiguration::nPhiBins().

Referenced by OMTFProcessor< GoldenPatternType >::fillCounts(), and OMTFProcessor< GoldenPatternType >::processInput().

                                                                                     {
  OMTFinput::vector1D myHits = layerHits;

  unsigned int iStart = myOmtfConfig->getConnections()[iProcessor][iRegion][iLayer].first;
  unsigned int iEnd = iStart + myOmtfConfig->getConnections()[iProcessor][iRegion][iLayer].second - 1;

  for (unsigned int iInput = 0; iInput < 14; ++iInput) {
    if (iInput < iStart || iInput > iEnd)
      myHits[iInput] = myOmtfConfig->nPhiBins();
  }
  return myHits;
}

run()

template<class GoldenPatternType >

std::vector< l1t::RegionalMuonCand > OMTFProcessor< GoldenPatternType >::run ( unsigned int  iProcessor, l1t::tftype  mtfType, int  bx, OMTFinputMaker *  inputMaker, std::vector< std::unique_ptr< IOMTFEmulationObserver > > &  observers  )

overridevirtual

Implements IProcessorEmulator.

Definition at line 762 of file OMTFProcessor.cc.

References OMTFinputMaker::buildInputForProcessor(), nano_mu_digi_cff::bx, submitPVResolutionJobs::count, input, LogTrace, and RPCpg::pts.

                                                                  {
  //uncomment if you want to check execution time of each method
  //boost::timer::auto_cpu_timer t("%ws wall, %us user in getProcessorCandidates\n");

  for (auto& obs : observers)
    obs->observeProcesorBegin(iProcessor, mtfType);

  //input is shared_ptr because the observers may need them after the run() method execution is finished
  std::shared_ptr<OMTFinput> input = std::make_shared<OMTFinput>(this->myOmtfConfig);
  inputMaker->buildInputForProcessor(input->getMuonStubs(), iProcessor, mtfType, bx, bx, observers);

  if (this->myOmtfConfig->cleanStubs()) {
    //this has sense for the pattern generation from the tracks with the secondaries
    //if more than one stub is in a given layer, all stubs are removed from this layer
    for (unsigned int iLayer = 0; iLayer < input->getMuonStubs().size(); ++iLayer) {
      auto& layerStubs = input->getMuonStubs()[iLayer];
      int count = std::count_if(layerStubs.begin(), layerStubs.end(), [](auto& ptr) { return ptr != nullptr; });
      if (count > 1) {
        for (auto& ptr : layerStubs)
          ptr.reset();

        LogTrace("OMTFReconstruction") << __FUNCTION__ << ":" << __LINE__ << "cleaning stubs in the layer " << iLayer
                                       << " stubs count :" << count << std::endl;
      }
    }
  }

  //LogTrace("l1tOmtfEventPrint")<<"buildInputForProce "; t.report();
  processInput(iProcessor, mtfType, *(input.get()), observers);

  //LogTrace("l1tOmtfEventPrint")<<"processInput       "; t.report();
  AlgoMuons algoCandidates = sortResults(iProcessor, mtfType);

  if (ptAssignment) {
    for (auto& myCand : algoCandidates) {
      if (myCand->isValid()) {
        auto pts = ptAssignment->getPts(myCand, observers);
        /*for (unsigned int i = 0; i < pts.size(); i++) {
        trackAddr[10 + i] = this->myOmtfConfig->ptGevToHw(pts[i]);
      }*/
      }
    }
  }

  //LogTrace("l1tOmtfEventPrint")<<"sortResults        "; t.report();
  // perform GB
  //watch out: etaBits2HwEta is used in the ghostBust to convert the AlgoMuons eta, it affect algoCandidates as they are pointers
  AlgoMuons gbCandidates = ghostBust(algoCandidates);

  //LogTrace("l1tOmtfEventPrint")<<"ghostBust"; t.report();
  // fill RegionalMuonCand colleciton
  std::vector<l1t::RegionalMuonCand> candMuons = getFinalcandidates(iProcessor, mtfType, gbCandidates);

  //LogTrace("l1tOmtfEventPrint")<<"getFinalcandidates "; t.report();
  //fill outgoing collection
  for (auto& candMuon : candMuons) {
    candMuon.setHwQual(candMuon.hwQual());
  }

  for (auto& obs : observers) {
    obs->observeProcesorEmulation(iProcessor, mtfType, input, algoCandidates, gbCandidates, candMuons);
  }

  return candMuons;
}

saveExtrapolFactors()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::saveExtrapolFactors

(

)

Definition at line 841 of file OMTFProcessor.cc.

References relativeConstraints::empty, to_string(), compare::tree, simOmtfPhase2Digis_cfi::useEndcapStubsRInExtr, simOmtfPhase2Digis_cfi::useStubQualInExtr, and MCScenario_CRAFT1_22X::write_xml().

                                                           {
  //if(this->myOmtfConfig->nProcessors() == 3) //phase2
  extrapolMultiplier = 512;

  boost::property_tree::ptree tree;
  auto& extrFactorsTree = tree.add("ExtrapolationFactors", "");
  extrFactorsTree.add("<xmlattr>.multiplier", extrapolMultiplier);

  edm::LogVerbatim("OMTFReconstruction") << "saving extrapolFactors to ExtrapolationFactors.xml" << std::endl;
  for (unsigned int iRefLayer = 0; iRefLayer < extrapolFactors.size(); iRefLayer++) {
    for (unsigned int iLayer = 0; iLayer < extrapolFactors[iRefLayer].size(); iLayer++) {
      edm::LogVerbatim("OMTFReconstruction") << " iRefLayer " << iRefLayer << " iLayer " << iLayer << std::endl;

      auto& layerTree = extrFactorsTree.add_child("Lut", boost::property_tree::ptree());
      layerTree.add("<xmlattr>.RefLayer", std::to_string(iRefLayer));
      layerTree.add("<xmlattr>.Layer", iLayer);

      if (useStubQualInExtr && (iLayer == 0 || iLayer == 2 || iLayer == 4))
        layerTree.add("<xmlattr>.KeyType", "quality");
      else if (useEndcapStubsRInExtr && ((iLayer >= 6 && iLayer <= 9) || (iLayer >= 15 && iLayer <= 17)))
        layerTree.add("<xmlattr>.KeyType", "eta");
      else
        layerTree.add("<xmlattr>.KeyType", "none");

      for (auto& extrFactors : extrapolFactors[iRefLayer][iLayer]) {
        int norm = 1;
        if (!extrapolFactorsNorm[iRefLayer][iLayer].empty())
          norm = extrapolFactorsNorm[iRefLayer][iLayer][extrFactors.first];
        auto& lutVal = layerTree.add_child("LutVal", boost::property_tree::ptree());
        if (useEndcapStubsRInExtr && ((iLayer >= 6 && iLayer <= 9) || (iLayer >= 15 && iLayer <= 17)))
          lutVal.add("<xmlattr>.key", extrFactors.first);
        else
          lutVal.add("<xmlattr>.key", extrFactors.first);

        double value = round(extrapolMultiplier * extrFactors.second / norm);
        lutVal.add("<xmlattr>.value", value);

        edm::LogVerbatim("OMTFReconstruction")
            << std::setw(4) << " key = " << extrFactors.first << " extrFactors.second " << std::setw(10)
            << extrFactors.second << " norm " << std::setw(6) << norm << " value/norm " << std::setw(10)
            << extrFactors.second / norm << " value " << value << std::endl;
      }
    }
  }

  boost::property_tree::write_xml("ExtrapolationFactors.xml",
                                  tree,
                                  std::locale(),
                                  boost::property_tree::xml_parser::xml_writer_make_settings<std::string>(' ', 2));
}

setGhostBuster()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::setGhostBuster ( IGhostBuster *  ghostBuster )

inlineoverridevirtual

allows to use other IGhostBuster implementation than the default one

Implements IProcessorEmulator.

Definition at line 109 of file OMTFProcessor.h.

References OMTFProcessor< GoldenPatternType >::ghostBuster.

{ this->ghostBuster.reset(ghostBuster); }

setPtAssignment()

template<class GoldenPatternType >

virtual void OMTFProcessor< GoldenPatternType >::setPtAssignment ( PtAssignmentBase *  ptAssignment )

inlinevirtual

Definition at line 111 of file OMTFProcessor.h.

References OMTFProcessor< GoldenPatternType >::ptAssignment.

Referenced by OmtfEmulation::addObservers().

{ this->ptAssignment = ptAssignment; }

setSorter()

template<class GoldenPatternType >

virtual void OMTFProcessor< GoldenPatternType >::setSorter ( SorterBase< GoldenPatternType > *  sorter )

inlinevirtual

allows to use other sorter implementation than the default one

Definition at line 106 of file OMTFProcessor.h.

{ this->sorter.reset(sorter); }

shiftGP()

template<class GoldenPatternType >

void OMTFProcessor< GoldenPatternType >::shiftGP ( GoldenPattern *  aGP, const GoldenPattern::vector2D &  meanDistPhiNew, const GoldenPattern::vector2D &  meanDistPhiOld  )

private

Shift pdf indexes by differecne between averaged and original meanDistPhi

Shift pdfs by differecne between original menaDistPhi, and the averaged value

Definition at line 182 of file OMTFProcessor.cc.

References GoldenPattern::getPdf(), OMTFProcessor< GoldenPatternType >::myOmtfConfig, OMTFConfiguration::nLayers(), OMTFConfiguration::nPdfAddrBits(), OMTFConfiguration::nRefLayers(), GoldenPattern::pdfValue(), and GoldenPattern::setPdf().

Referenced by OMTFProcessor< GoldenPatternType >::averagePatterns().

                                                                         {
  unsigned int nPdfBins = exp2(myOmtfConfig->nPdfAddrBits());
  GoldenPattern::vector3D pdfAllRef = aGP->getPdf();

  int indexShift = 0;
  for (unsigned int iLayer = 0; iLayer < myOmtfConfig->nLayers(); ++iLayer) {
    for (unsigned int iRefLayer = 0; iRefLayer < myOmtfConfig->nRefLayers(); ++iRefLayer) {
      indexShift = meanDistPhiOld[iLayer][iRefLayer] - meanDistPhiNew[iLayer][iRefLayer];
      for (unsigned int iPdfBin = 0; iPdfBin < nPdfBins; ++iPdfBin)
        pdfAllRef[iLayer][iRefLayer][iPdfBin] = 0;
      for (unsigned int iPdfBin = 0; iPdfBin < nPdfBins; ++iPdfBin) {
        if ((int)(iPdfBin) + indexShift >= 0 && iPdfBin + indexShift < nPdfBins)
          pdfAllRef[iLayer][iRefLayer][iPdfBin + indexShift] = aGP->pdfValue(iLayer, iRefLayer, iPdfBin);
      }
    }
  }
  aGP->setPdf(pdfAllRef);
}

sortResults()

template<class GoldenPatternType >

AlgoMuons OMTFProcessor< GoldenPatternType >::sortResults ( unsigned int  iProcessor, l1t::tftype  mtfType, int  charge = 0  )

overridevirtual

Implements IProcessorEmulator.

Definition at line 343 of file OMTFProcessor.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge.

                                                                                                              {
  unsigned int procIndx = this->myOmtfConfig->getProcIndx(iProcessor, mtfType);
  return sorter->sortResults(procIndx, this->getPatterns(), charge);
}

Member Data Documentation

extrapolFactors

template<class GoldenPatternType >

std::vector<std::vector<std::map<int, double> > > OMTFProcessor< GoldenPatternType >::extrapolFactors

private

Definition at line 151 of file OMTFProcessor.h.

extrapolFactorsNorm

template<class GoldenPatternType >

std::vector<std::vector<std::map<int, int> > > OMTFProcessor< GoldenPatternType >::extrapolFactorsNorm

private

Definition at line 152 of file OMTFProcessor.h.

extrapolMultiplier

template<class GoldenPatternType >

int OMTFProcessor< GoldenPatternType >::extrapolMultiplier = 128

private

Definition at line 149 of file OMTFProcessor.h.

ghostBuster

template<class GoldenPatternType >

std::unique_ptr<IGhostBuster> OMTFProcessor< GoldenPatternType >::ghostBuster

private

Definition at line 137 of file OMTFProcessor.h.

Referenced by OMTFProcessor< GoldenPatternType >::ghostBust(), and OMTFProcessor< GoldenPatternType >::setGhostBuster().

myOmtfConfig

template<class GoldenPatternType >

const OMTFConfiguration* OMTFProcessor< GoldenPatternType >::myOmtfConfig

private

Configuration of the algorithm. This object does not contain the patterns data.

Definition at line 86 of file OMTFProcessor.h.

Referenced by OMTFProcessor< GoldenPatternType >::addGP(), OMTFProcessor< GoldenPatternType >::averagePatterns(), OMTFProcessor< GoldenPatternType >::configure(), OMTFProcessor< GoldenPatternType >::fillCounts(), OMTFProcessor< GoldenPatternType >::processInput(), OMTFProcessor< GoldenPatternType >::restrictInput(), and OMTFProcessor< GoldenPatternType >::shiftGP().

myResults

template<class GoldenPatternType >

std::vector<OMTFProcessor::resultsMap> OMTFProcessor< GoldenPatternType >::myResults

private

Map holding results on current event data for each GP. Reference hit number is isued as a vector index.

Definition at line 82 of file OMTFProcessor.h.

Referenced by OMTFProcessor< GoldenPatternType >::addGP(), OMTFProcessor< GoldenPatternType >::configure(), OMTFProcessor< GoldenPatternType >::processInput(), and OMTFProcessor< GoldenPatternType >::resetConfiguration().

ptAssignment

template<class GoldenPatternType >

PtAssignmentBase* OMTFProcessor< GoldenPatternType >::ptAssignment = nullptr

private

Definition at line 140 of file OMTFProcessor.h.

Referenced by OMTFProcessor< GoldenPatternType >::setPtAssignment().

sorter

template<class GoldenPatternType >

std::unique_ptr<SorterBase<GoldenPatternType> > OMTFProcessor< GoldenPatternType >::sorter

private

Definition at line 135 of file OMTFProcessor.h.

theGPs

template<class GoldenPatternType >

std::map<Key, GoldenPattern *> OMTFProcessor< GoldenPatternType >::theGPs

private

Map holding Golden Patterns.

Definition at line 77 of file OMTFProcessor.h.

Referenced by OMTFProcessor< GoldenPatternType >::addGP(), OMTFProcessor< GoldenPatternType >::averagePatterns(), OMTFProcessor< GoldenPatternType >::fillCounts(), OMTFProcessor< GoldenPatternType >::getPatterns(), OMTFProcessor< GoldenPatternType >::processInput(), OMTFProcessor< GoldenPatternType >::resetConfiguration(), and OMTFProcessor< GoldenPatternType >::~OMTFProcessor().

useEndcapStubsRInExtr

template<class GoldenPatternType >

bool OMTFProcessor< GoldenPatternType >::useEndcapStubsRInExtr = false

private

Definition at line 143 of file OMTFProcessor.h.

useFloatingPointExtrapolation

template<class GoldenPatternType >

bool OMTFProcessor< GoldenPatternType >::useFloatingPointExtrapolation = false

private

Definition at line 147 of file OMTFProcessor.h.

useStubQualInExtr

template<class GoldenPatternType >

bool OMTFProcessor< GoldenPatternType >::useStubQualInExtr = false

private

Definition at line 142 of file OMTFProcessor.h.