GhostBuster Class Reference

#include <GhostBuster.h>

Inheritance diagram for GhostBuster:

Public Member Functions
	GhostBuster (const OMTFConfiguration *omtfConfig)
AlgoMuons	select (AlgoMuons refHitCands, int charge=0) override
std::vector< AlgoMuon >	select (std::vector< AlgoMuon > refHitCands, int charge=0) override
	~GhostBuster () override
	~GhostBuster () override
Public Member Functions inherited from IGhostBuster
virtual	~IGhostBuster ()
virtual	~IGhostBuster ()

Private Attributes
const OMTFConfiguration *	omtfConfig

Detailed Description

Definition at line 15 of file GhostBuster.h.

Constructor & Destructor Documentation

~GhostBuster() [1/2]

GhostBuster::~GhostBuster ( )

inlineoverride

Definition at line 17 of file GhostBuster.h.

{};

GhostBuster()

GhostBuster::GhostBuster ( const OMTFConfiguration *  omtfConfig )

inline

Definition at line 12 of file GhostBuster.h.

: omtfConfig(omtfConfig){};

~GhostBuster() [2/2]

GhostBuster::~GhostBuster ( )

inlineoverride

Definition at line 14 of file GhostBuster.h.

{};

Member Function Documentation

select() [1/2]

AlgoMuons GhostBuster::select ( AlgoMuons  refHitCands, int  charge = 0  )

overridevirtual

Implements IGhostBuster.

Definition at line 8 of file GhostBuster.cc.

References a, funct::abs(), b, ProcConfigurationBase::bits, ALCARECOTkAlJpsiMuMu_cff::charge, OMTFConfiguration::etaBits2HwEta(), ProcConfigurationBase::getStubEtaEncoding(), omtfConfig, OMTFConfiguration::procPhiToGmtPhi(), and jetUpdater_cfi::sort.

                                                               {
  //edm::LogImportant("OMTFReconstruction")<<"calling "<<__PRETTY_FUNCTION__ <<std::endl;

  AlgoMuons refHitCleanCands;
  // Sort candidates with decreased goodness,
  auto customLess = [&](const AlgoMuons::value_type& a, const AlgoMuons::value_type& b) -> bool {
    if (!a->isValid()) {
      return true;
    }
    if (!b->isValid()) {
      return false;
    }

    if (a->getQ() > b->getQ())
      return false;
    else if (a->getQ() == b->getQ() && a->getDisc() > b->getDisc())
      return false;
    else if (a->getQ() == b->getQ() && a->getDisc() == b->getDisc() &&
             a->getPatternNumConstr() > b->getPatternNumConstr())
      return false;
    else if (a->getQ() == b->getQ() && a->getDisc() == b->getDisc() &&
             a->getPatternNumConstr() == b->getPatternNumConstr() && a->getRefHitNumber() < b->getRefHitNumber())
      return false;
    else
      return true;
  };

  std::sort(refHitCands.rbegin(), refHitCands.rend(), customLess);

  for (AlgoMuons::iterator it1 = refHitCands.begin(); it1 != refHitCands.end(); ++it1) {
    bool isGhost = false;
    for (AlgoMuons::iterator it2 = refHitCleanCands.begin(); it2 != refHitCleanCands.end(); ++it2) {
      //do not accept candidates with similar phi (any charge combination)
      //veto window 5deg(=half of logic cone)=5/360*5760=80"logic strips"
      //veto window 5 degree in GMT scale is 5/360*576=8 units
      if (std::abs(omtfConfig->procPhiToGmtPhi((*it1)->getPhi()) - omtfConfig->procPhiToGmtPhi((*it2)->getPhi())) < 8) {
        //      if(std::abs(it1->getPhi() - it2->getPhi())<5/360.0*nPhiBins){
        isGhost = true;
        break;
        //which one candidate is killed depends only on the order in the refHitCands (the one with smaller index is taken), and this order is assured by the sort above
        //TODO here the candidate that is killed does not kill other candidates - check if the firmware does the same (KB)
      }
    }
    if ((*it1)->getQ() > 0 && !isGhost) {
      refHitCleanCands.emplace_back(new AlgoMuon(**it1));
      if (omtfConfig->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::bits)
        refHitCleanCands.back()->setEta(OMTFConfiguration::etaBits2HwEta((*it1)->getEtaHw()));
    }

    if (refHitCleanCands.size() >= 3)
      break;
  }

  while (refHitCleanCands.size() < 3)
    refHitCleanCands.emplace_back(new AlgoMuon());

  std::stringstream myStr;
  bool hasCandidates = false;
  for (unsigned int iRefHit = 0; iRefHit < refHitCands.size(); ++iRefHit) {
    if (refHitCands[iRefHit]->getQ()) {
      hasCandidates = true;
      break;
    }
  }
  for (unsigned int iRefHit = 0; iRefHit < refHitCands.size(); ++iRefHit) {
    if (refHitCands[iRefHit]->getQ())
      myStr << "Ref hit: " << iRefHit << " " << refHitCands[iRefHit] << std::endl;
  }
  myStr << "Selected Candidates with charge: " << charge << std::endl;
  for (unsigned int iCand = 0; iCand < refHitCleanCands.size(); ++iCand) {
    myStr << "Cand: " << iCand << " " << refHitCleanCands[iCand] << std::endl;
  }

  if (hasCandidates)
    edm::LogInfo("OMTF Sorter") << myStr.str();

  // update refHitCands with refHitCleanCands
  return refHitCleanCands;
}

select() [2/2]

std::vector< AlgoMuon > GhostBuster::select ( std::vector< AlgoMuon >  refHitCands, int  charge = 0  )

overridevirtual

Implements IGhostBuster.

Definition at line 13 of file GhostBuster.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, and jetUpdater_cfi::sort.

                                                                                   {
  std::vector<AlgoMuon> refHitCleanCands;
  // Sort candidates with decreased goodness,
  // where goodness definied in < operator of AlgoMuon
  std::sort(refHitCands.rbegin(), refHitCands.rend());

  for (std::vector<AlgoMuon>::iterator it1 = refHitCands.begin(); it1 != refHitCands.end(); ++it1) {
    bool isGhost = false;
    for (std::vector<AlgoMuon>::iterator it2 = refHitCleanCands.begin(); it2 != refHitCleanCands.end(); ++it2) {
      //do not accept candidates with similar phi (any charge combination)
      //veto window 5deg(=half of logic cone)=5/360*5760=80"logic strips"
      //veto window 5 degree in GMT scale is 5/360*576=8 units
      if (std::abs(phiGMT(it1->getPhi()) - phiGMT(it2->getPhi())) < 8) {
        //      if(std::abs(it1->getPhi() - it2->getPhi())<5/360.0*nPhiBins){
        isGhost = true;
        break;
        //which one candidate is killed depends only on the order in the refHitCands (the one with smaller index is taken), and this order is assured by the sort above
        //TODO here the candidate that is killed does not kill other candidates - check if the firmware does the same (KB)
      }
    }
    if (it1->getQ() > 0 && !isGhost)
      refHitCleanCands.push_back(*it1);
  }

  refHitCleanCands.resize(3, AlgoMuon(0, 999, 9999, 0, 0, 0, 0, 0));  //FIXME
  //refHitCleanCands.resize( 3, AlgoMuon() );

  std::stringstream myStr;
  bool hasCandidates = false;
  for (unsigned int iRefHit = 0; iRefHit < refHitCands.size(); ++iRefHit) {
    if (refHitCands[iRefHit].getQ()) {
      hasCandidates = true;
      break;
    }
  }
  for (unsigned int iRefHit = 0; iRefHit < refHitCands.size(); ++iRefHit) {
    if (refHitCands[iRefHit].getQ())
      myStr << "Ref hit: " << iRefHit << " " << refHitCands[iRefHit] << std::endl;
  }
  myStr << "Selected Candidates with charge: " << charge << std::endl;
  for (unsigned int iCand = 0; iCand < refHitCleanCands.size(); ++iCand) {
    myStr << "Cand: " << iCand << " " << refHitCleanCands[iCand] << std::endl;
  }

  if (hasCandidates)
    edm::LogInfo("OMTF Sorter") << myStr.str();

  // update refHitCands with refHitCleanCands
  return refHitCleanCands;
}

Member Data Documentation

omtfConfig

const OMTFConfiguration* GhostBuster::omtfConfig

private

Definition at line 9 of file GhostBuster.h.

Referenced by select().