OMTFProcessor Class Reference

#include <OMTFProcessor.h>

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

Public Member Functions
void	averagePatterns (int charge)
bool	configure (XMLConfigReader *aReader)
	Fill GP map with patterns from XML file. More...
bool	configure (std::shared_ptr< L1TMuonOverlapParams > omtfParams)
	Fill GP map with patterns from CondFormats object. More...
void	fillCounts (unsigned int iProcessor, const OMTFinput &aInput, const SimTrack *aSimMuon)
const std::map< Key, GoldenPattern * > &	getPatterns () const
	Return map of GoldenPatterns. More...
	OMTFProcessor (const edm::ParameterSet &cfg)
const std::vector < OMTFProcessor::resultsMap > &	processInput (unsigned int iProcessor, const OMTFinput &aInput)
OMTFinput	shiftInput (unsigned int iProcessor, const OMTFinput &aInput)
	~OMTFProcessor ()

Private Member Functions
bool	addGP (GoldenPattern *aGP)
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)
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 < OMTFProcessor::resultsMap >	myResults
std::map< Key, GoldenPattern * >	theGPs
	Map holding Golden Patterns. More...

Detailed Description

Definition at line 20 of file OMTFProcessor.h.

Member Typedef Documentation

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

Definition at line 24 of file OMTFProcessor.h.

Constructor & Destructor Documentation

OMTFProcessor::OMTFProcessor

(

const edm::ParameterSet &

cfg

)

Definition at line 25 of file OMTFProcessor.cc.

References configure(), edm::ParameterSet::exists(), LaserDQM_cfg::fileNames, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), resetConfiguration(), and XMLConfigReader::setPatternsFile().

                                                             {

  if(theConfig.getParameter<bool>("configFromXML")){    
    if ( !theConfig.exists("patternsXMLFiles") ) return;
    std::vector<std::string> fileNames;
    for(auto it: theConfig.getParameter<std::vector<edm::ParameterSet> >("patternsXMLFiles")){
      fileNames.push_back(it.getParameter<edm::FileInPath>("patternsXMLFile").fullPath());
    }  

    resetConfiguration();
    XMLConfigReader myReader;
    for(auto it: fileNames){
      myReader.setPatternsFile(it);
      configure(&myReader);
    }
  }
}

OMTFProcessor::~OMTFProcessor

(

)

Definition at line 44 of file OMTFProcessor.cc.

References theGPs.

                             {

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

}

Member Function Documentation

bool OMTFProcessor::addGP ( GoldenPattern *  aGP )

private

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

Definition at line 130 of file OMTFProcessor.cc.

References Exception, GoldenPattern::key(), myResults, and theGPs.

Referenced by 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()] = new GoldenPattern(*aGP);

  for(auto & itRegion: myResults) itRegion[aGP->key()] = OMTFResult(); 

  return true;
}

void OMTFProcessor::averagePatterns

(

int

charge

)

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

Definition at line 146 of file OMTFProcessor.cc.

References GoldenPattern::getMeanDistPhi(), OMTFConfiguration::nLayers, OMTFConfiguration::nRefLayers, GoldenPattern::setMeanDistPhi(), shiftGP(), theGPs, and Key::thePtCode.

Referenced by L1TMuonOverlapTrackProducer::endJob().

                                              {

  Key aKey(1, 4, charge);

  while(theGPs.find(aKey)!=theGPs.end()){

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

    ++aKey.thePtCode;
    if(aKey.thePtCode<=31 && theGPs.find(aKey)!=theGPs.end()) aGP2 =  theGPs.find(aKey)->second;
    
    if(aKey.thePtCode>19){
      ++aKey.thePtCode;
      if(aKey.thePtCode<=31 && theGPs.find(aKey)!=theGPs.end()) aGP3 =  theGPs.find(aKey)->second;
      
      ++aKey.thePtCode;
      if(aKey.thePtCode<=31 && theGPs.find(aKey)!=theGPs.end()) aGP4 =  theGPs.find(aKey)->second;
    }
    else{
      aGP3 = aGP1;
      aGP4 = aGP2;
    }
    ++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();
    /*
     std::cout<<"Key: "
         <<aGP1->key()
         <<" "<<aGP2->key()
         <<" "<<aGP3->key()
         <<" "<<aGP4->key()<<std::endl;
    */
    for(unsigned int iLayer=0;iLayer<OMTFConfiguration::nLayers;++iLayer){
      for(unsigned int iRefLayer=0;iRefLayer<OMTFConfiguration::nRefLayers;++iRefLayer){
    meanDistPhi[iLayer][iRefLayer]+=meanDistPhi2[iLayer][iRefLayer];
    meanDistPhi[iLayer][iRefLayer]+=meanDistPhi3[iLayer][iRefLayer];
    meanDistPhi[iLayer][iRefLayer]+=meanDistPhi4[iLayer][iRefLayer];
    meanDistPhi[iLayer][iRefLayer]/=4;
    /*
    std::cout<<"Mean distPhi: "<<meanDistPhi[iLayer][iRefLayer]
         <<" "<<meanDistPhi1[iLayer][iRefLayer]
         <<" "<<meanDistPhi2[iLayer][iRefLayer]
         <<" "<<meanDistPhi3[iLayer][iRefLayer]
         <<" "<<meanDistPhi4[iLayer][iRefLayer]<<std::endl; 
    */
      }
    }
    
    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);   
    }
  }
  
}

bool OMTFProcessor::configure

(

XMLConfigReader *

aReader

)

Fill GP map with patterns from XML file.

Definition at line 58 of file OMTFProcessor.cc.

References addGP(), myResults, OMTFConfiguration::nTestRefHits, XMLConfigReader::readPatterns(), and resetConfiguration().

Referenced by L1TMuonOverlapTrackProducer::beginRun(), and OMTFProcessor().

                                                     {

  resetConfiguration();

  myResults.assign(OMTFConfiguration::nTestRefHits,OMTFProcessor::resultsMap());

  const std::vector<GoldenPattern *> & aGPs = aReader->readPatterns();
  for(auto it: aGPs){    
    if(!addGP(it)) return false;
  }
  
  return true;
}

bool OMTFProcessor::configure

(

std::shared_ptr< L1TMuonOverlapParams >

omtfParams

)

Fill GP map with patterns from CondFormats object.

Mean dist phi data

Pdf data

Definition at line 73 of file OMTFProcessor.cc.

References addGP(), l1t::LUT::data(), myResults, OMTFConfiguration::nGoldenPatterns, OMTFConfiguration::nLayers, OMTFConfiguration::nPdfAddrBits, l1t::LUT::nrBitsData(), OMTFConfiguration::nRefLayers, OMTFConfiguration::nTestRefHits, resetConfiguration(), GoldenPattern::setMeanDistPhi(), and GoldenPattern::setPdf().

                                                                            {

  resetConfiguration();
  
  myResults.assign(OMTFConfiguration::nTestRefHits,OMTFProcessor::resultsMap());
  
  l1t::LUT* chargeLUT =  omtfParams->chargeLUT();
  l1t::LUT* etaLUT =  omtfParams->etaLUT();
  l1t::LUT* ptLUT =  omtfParams->ptLUT();
  l1t::LUT* pdfLUT =  omtfParams->pdfLUT();
  l1t::LUT* meanDistPhiLUT =  omtfParams->meanDistPhiLUT();

  unsigned int nGPs = OMTFConfiguration::nGoldenPatterns;
  unsigned int address = 0;
  unsigned int iEta, iPt, iCharge;
  for(unsigned int iGP=0;iGP<nGPs;++iGP){
    address = iGP;
    iEta = etaLUT->data(address);
    iCharge = chargeLUT->data(address);
    if(iCharge==0) iCharge = -1;//FIXME: Standarise to 0,1
    iPt = ptLUT->data(address);

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

void OMTFProcessor::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.

Number of reference hits to be checked. Value read from XML configuration

Definition at line 346 of file OMTFProcessor.cc.

References funct::abs(), OMTFConfiguration::bendingLayers, KineDebug3::count(), OMTFinput::getLayerData(), OMTFinput::getRefHits(), RefHitDef::iInput, RPCConst::iptFromPt(), RefHitDef::iRefLayer, RefHitDef::iRegion, CoreSimTrack::momentum(), OMTFConfiguration::nLayers, OMTFConfiguration::nRefHits, OMTFConfiguration::refHitsDefs, OMTFConfiguration::refToLogicNumber, restrictInput(), theGPs, and CoreSimTrack::type().

                                            {

  int theCharge = (abs(aSimMuon->type()) == 13) ? aSimMuon->type()/-13 : 0;
  unsigned int  iPt =  RPCConst::iptFromPt(aSimMuon->momentum().pt());

  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<OMTFConfiguration::nLayers;++iLayer){
    const OMTFinput::vector1D & layerHits = aInput.getLayerData(iLayer);
    if(!layerHits.size()) continue;
    for(unsigned int iRefHit=0;iRefHit<OMTFConfiguration::nRefHits;++iRefHit){
      if(!refHitsBits[iRefHit]) continue;
      const RefHitDef & aRefHitDef = OMTFConfiguration::refHitsDefs[iProcessor][iRefHit];
      int phiRef = aInput.getLayerData(OMTFConfiguration::refToLogicNumber[aRefHitDef.iRefLayer])[aRefHitDef.iInput]; 
      unsigned int iRegion = aRefHitDef.iRegion;
      if(OMTFConfiguration::bendingLayers.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);
      }
    }
  }
}

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

private

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

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

private

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

(

)

const

Return map of GoldenPatterns.

Definition at line 245 of file OMTFProcessor.cc.

References theGPs.

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

{ return theGPs; }

const std::vector< OMTFProcessor::resultsMap > & OMTFProcessor::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. Value read from XML configuration

Definition at line 248 of file OMTFProcessor.cc.

References OMTFConfiguration::bendingLayers, KineDebug3::count(), OMTFinput::getLayerData(), OMTFinput::getRefHits(), RefHitDef::iInput, RefHitDef::iRefLayer, RefHitDef::iRegion, myResults, OMTFConfiguration::nLayers, OMTFConfiguration::nRefHits, OMTFConfiguration::nTestRefHits, OMTFConfiguration::refHitsDefs, OMTFConfiguration::refToLogicNumber, restrictInput(), and theGPs.

Referenced by L1TMuonOverlapTrackProducer::produce().

                                                                {

  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<OMTFConfiguration::nLayers;++iLayer){
    const OMTFinput::vector1D & layerHits = aInput.getLayerData(iLayer);
    if(!layerHits.size()) continue;
    unsigned int nTestedRefHits = OMTFConfiguration::nTestRefHits;
    for(unsigned int iRefHit=0;iRefHit<OMTFConfiguration::nRefHits;++iRefHit){
      if(!refHitsBits[iRefHit]) continue;
      if(nTestedRefHits--==0) break;
      const RefHitDef & aRefHitDef = OMTFConfiguration::refHitsDefs[iProcessor][iRefHit];
      int phiRef = aInput.getLayerData(OMTFConfiguration::refToLogicNumber[aRefHitDef.iRefLayer])[aRefHitDef.iInput];
      int etaRef = aInput.getLayerData(OMTFConfiguration::refToLogicNumber[aRefHitDef.iRefLayer],true)[aRefHitDef.iInput];
      unsigned int iRegion = aRefHitDef.iRegion;
      if(OMTFConfiguration::bendingLayers.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[OMTFConfiguration::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();

  //#ifndef NDEBUG
  std::ostringstream myStr;
  myStr<<"iProcessor: "<<iProcessor<<std::endl;
  myStr<<"Input: ------------"<<std::endl;
  myStr<<aInput<<std::endl;
  /*
  for(auto itRefHit: myResults){
    myStr<<"--- Reference hit ---"<<std::endl;
    for(auto itKey: itRefHit){      
      if(itKey.second.empty()) continue;
      myStr<<itKey.first<<std::endl;
      myStr<<itKey.second<<std::endl;
    }
    myStr<<"--------------------"<<std::endl;
  }
  */
  //LogDebug("OMTF processor")<<myStr.str();
  edm::LogInfo("OMTF processor")<<myStr.str();
  //#endif
  
  return myResults;
}   

void OMTFProcessor::resetConfiguration ( )

private

Reset all configuration parameters.

Definition at line 51 of file OMTFProcessor.cc.

References myResults, and theGPs.

Referenced by configure(), and OMTFProcessor().

                                      {

  myResults.clear();
  theGPs.clear();
}

OMTFinput::vector1D OMTFProcessor::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 328 of file OMTFProcessor.cc.

References OMTFConfiguration::connections, and OMTFConfiguration::nPhiBins.

Referenced by fillCounts(), and processInput().

                                                             {

  OMTFinput::vector1D myHits = layerHits;
 
  unsigned int iStart = OMTFConfiguration::connections[iProcessor][iRegion][iLayer].first;
  unsigned int iEnd = iStart + OMTFConfiguration::connections[iProcessor][iRegion][iLayer].second -1;

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

void OMTFProcessor::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 220 of file OMTFProcessor.cc.

References GoldenPattern::getPdf(), OMTFConfiguration::nLayers, OMTFConfiguration::nPdfAddrBits, OMTFConfiguration::nRefLayers, GoldenPattern::pdfValue(), and GoldenPattern::setPdf().

Referenced by averagePatterns().

                                                             {

  unsigned int nPdfBins =  exp2(OMTFConfiguration::nPdfAddrBits);

  GoldenPattern::vector3D pdfAllRef = aGP->getPdf();

  int indexShift = 0;
  for(unsigned int iLayer=0;iLayer<OMTFConfiguration::nLayers;++iLayer){
    for(unsigned int iRefLayer=0;iRefLayer<OMTFConfiguration::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);
}

OMTFinput OMTFProcessor::shiftInput	(	unsigned int	iProcessor,
		const OMTFinput &	aInput
	)

Shift phi values in input to fit the 11 bits range. For each processor the global phi beggining-511 is added, so it starts at -551

OMTFConfiguration::nPhiBins/2 to shift the minPhi to 0-nBins scale,

Definition at line 313 of file OMTFProcessor.cc.

References OMTFConfiguration::globalPhiStart(), OMTFConfiguration::nPhiBins, and OMTFinput::shiftMyPhi().

Referenced by L1TMuonOverlapTrackProducer::produce().

                                             {

  int minPhi =  OMTFConfiguration::globalPhiStart(iProcessor);

  if(minPhi<0) minPhi+=OMTFConfiguration::nPhiBins;

  OMTFinput myCopy = aInput;
  myCopy.shiftMyPhi(minPhi);
  
  return myCopy;
}

Member Data Documentation

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

private

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

Definition at line 103 of file OMTFProcessor.h.

Referenced by addGP(), configure(), processInput(), and resetConfiguration().

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

private

Map holding Golden Patterns.

Definition at line 98 of file OMTFProcessor.h.

Referenced by addGP(), averagePatterns(), fillCounts(), getPatterns(), processInput(), resetConfiguration(), and ~OMTFProcessor().