l1t::Stage1Layer2TauAlgorithmImpHW Class Reference

#include <Stage1Layer2TauAlgorithmImp.h>

Inheritance diagram for l1t::Stage1Layer2TauAlgorithmImpHW:

Public Member Functions
virtual void	processEvent (const std::vector< l1t::CaloEmCand > &EMCands, const std::vector< l1t::CaloRegion > &regions, std::vector< l1t::Tau > *isoTaus, std::vector< l1t::Tau > *taus)
	Stage1Layer2TauAlgorithmImpHW (CaloParamsHelper *params)
virtual	~Stage1Layer2TauAlgorithmImpHW ()
Public Member Functions inherited from l1t::Stage1Layer2TauAlgorithm
virtual	~Stage1Layer2TauAlgorithm ()

Private Member Functions
int	AssociatedJetPt (int ieta, int iphi, const std::vector< l1t::Jet > *jets) const
string	findNESW (int ieta, int iphi, int neta, int nphi) const
unsigned	isoLutIndex (unsigned int tauPt, unsigned int jetPt) const
double	JetIsolation (int et, int ieta, int iphi, const std::vector< l1t::Jet > &jets) const

Private Attributes
Stage1TauIsolationLUT *	isoTauLut
CaloParamsHelper *const	params_

Additional Inherited Members
Public Attributes inherited from l1t::Stage1Layer2TauAlgorithm
std::vector< double >	regionPUSParams
std::string	regionPUSType

Detailed Description

Definition at line 54 of file Stage1Layer2TauAlgorithmImp.h.

Constructor & Destructor Documentation

Stage1Layer2TauAlgorithmImpHW::Stage1Layer2TauAlgorithmImpHW

(

CaloParamsHelper *

params

)

Definition at line 25 of file Stage1Layer2TauAlgorithmImpHW.cc.

References isoTauLut, and params_.

                                                                                     : params_(params)
{

  isoTauLut = new Stage1TauIsolationLUT(params_);
}

Stage1Layer2TauAlgorithmImpHW::~Stage1Layer2TauAlgorithmImpHW ( )

virtual

Definition at line 31 of file Stage1Layer2TauAlgorithmImpHW.cc.

References isoTauLut.

                                                             {
  delete isoTauLut;
};

Member Function Documentation

int l1t::Stage1Layer2TauAlgorithmImpHW::AssociatedJetPt ( int  ieta, int  iphi, const std::vector< l1t::Jet > *  jets  ) const

private

Definition at line 250 of file Stage1Layer2TauAlgorithmImpHW.cc.

References gather_cfg::cout, Debug, reco::btau::jetEta, reco::btau::jetPhi, and EnergyCorrector::pt.

                                                                         {

  bool Debug=false;

  if (Debug) cout << "Number of jets: " << jets->size() << endl;
  int pt = -1;


  for(JetBxCollection::const_iterator itJet = jets->begin();
      itJet != jets->end(); ++itJet){

    int jetEta = itJet->hwEta();
    int jetPhi = itJet->hwPhi();
    if (Debug) cout << "Matching ETA: " << ieta << " " << jetEta << endl;
    if (Debug) cout << "Matching PHI: " << iphi << " " << jetPhi << endl;
    if ((jetEta == ieta) && (jetPhi == iphi)){
      pt = itJet->hwPt();
      break;
    }
  }

  // set output
  return pt;
}

string l1t::Stage1Layer2TauAlgorithmImpHW::findNESW ( int  ieta, int  iphi, int  neta, int  nphi  ) const

private

Definition at line 221 of file Stage1Layer2TauAlgorithmImpHW.cc.

References funct::abs(), HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, and L1CaloRegionDetId::N_PHI.

                                                                                              {

  int deltaPhi = iphi - nphi;
  if (std::abs(deltaPhi) == L1CaloRegionDetId::N_PHI-1)
    deltaPhi = -deltaPhi/std::abs(deltaPhi); //18 regions in phi

  int deltaEta = ieta - neta;

  if ((std::abs(deltaPhi) +  std::abs(deltaEta)) < 2) {
    if (deltaEta==-1) {
      return "isEast";
    }
    else if (deltaEta==0) {
      if (deltaPhi==-1) {
    return "isNorth";
      }
      if (deltaPhi==1) {
    return "isSouth";
      }
    }
    else {
      return "isWest";
    }
  }

  return "999";

}

unsigned l1t::Stage1Layer2TauAlgorithmImpHW::isoLutIndex ( unsigned int  tauPt, unsigned int  jetPt  ) const

private

Definition at line 276 of file Stage1Layer2TauAlgorithmImpHW.cc.

References funct::pow().

{
  //const unsigned int nbitsTau=9;  // number of bits used for et in LUT file (needed for left shift operation)
  //const unsigned int nbitsJet=9;

  const unsigned int nbitsTau=8;  // number of bits used for et in LUT file (needed for left shift operation)
  const unsigned int nbitsJet=8;

  const unsigned int maxJet = pow(2,nbitsJet)-1;
  const unsigned int maxTau = pow(2,nbitsTau)-1;

  if (nbitsTau < 9)
  {
    if (nbitsTau == 6)
      {
      tauPt=tauPt>>3;
      }
    else if (nbitsTau == 7)
      {
      tauPt=tauPt>>2;
      }
    else if (nbitsTau == 8)
      {
    tauPt=tauPt>>1;
      }
  }

  if (nbitsJet < 9)// no need to do shift if nbits>=9
  {
    if (nbitsJet == 6)
      {
      jetPt=jetPt>>3;
      }
    else if (nbitsJet == 7)
      {
      jetPt=jetPt>>2;
      }
    else if (nbitsJet == 8)
      {
    jetPt=jetPt>>1;
      }
  }

  if (jetPt>maxJet) jetPt=maxJet;
  if (tauPt>maxTau) tauPt=maxTau;

  unsigned int address= (jetPt << nbitsTau) + tauPt;

  // std::cout << address << "\t## " << tauPt << " " << jetPt << std::endl;
  return address;
}

double l1t::Stage1Layer2TauAlgorithmImpHW::JetIsolation ( int  et, int  ieta, int  iphi, const std::vector< l1t::Jet > &  jets  ) const

private

Definition at line 202 of file Stage1Layer2TauAlgorithmImpHW.cc.

References metsig::jet.

                                                                  {

  for(JetBxCollection::const_iterator jet = jets.begin();
      jet != jets.end(); jet++) {

    if (ieta==jet->hwEta() && iphi==jet->hwPhi()){

      double isolation = (double) (jet->hwPt() - et);
      return isolation/et;
    }
  }

  // set output
  return 999.;
}

void l1t::Stage1Layer2TauAlgorithmImpHW::processEvent ( const std::vector< l1t::CaloEmCand > &  EMCands, const std::vector< l1t::CaloRegion > &  regions, std::vector< l1t::Tau > *  isoTaus, std::vector< l1t::Tau > *  taus  )

virtual

Implements l1t::Stage1Layer2TauAlgorithm.

Definition at line 38 of file Stage1Layer2TauAlgorithmImpHW.cc.

References funct::abs(), l1t::calibrateAndRankTaus(), gather_cfg::cout, HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, L1CaloRegionDetId::N_PHI, l1t::pack15bits(), HLT_25ns14e33_v1_cff::quality, HLT_25ns14e33_v3_cff::region, l1t::RegionCorrection(), l1t::SortTaus(), l1t::TauToGtEtaScales(), cond::to_string(), and l1t::TwelveByTwelveFinder().

                                                        {

  double towerLsb = params_->towerLsbSum();

  int tauSeedThreshold= floor( params_->tauSeedThreshold()/towerLsb + 0.5); // convert GeV to HW units
  int tauNeighbourThreshold= floor( params_->tauNeighbourThreshold()/towerLsb + 0.5); // convert GeV to HW units
  int tauMaxPtTauVeto = floor( params_->tauMaxPtTauVeto()/towerLsb + 0.5);
  int isoTauEtaMin = params_->isoTauEtaMin();
  int isoTauEtaMax = params_->isoTauEtaMax();

  std::vector<l1t::CaloRegion> *subRegions = new std::vector<l1t::CaloRegion>();


  //Region Correction will return uncorrected subregions if
  //regionPUSType is set to None in the config
  RegionCorrection(regions, subRegions, params_);



  // ----- need to cluster jets in order to compute jet isolation ----
  std::vector<l1t::Jet> *unCorrJets = new std::vector<l1t::Jet>();
  TwelveByTwelveFinder(0, subRegions, unCorrJets);

  std::vector<l1t::Tau> *preGtTaus = new std::vector<l1t::Tau>();
  std::vector<l1t::Tau> *preSortTaus = new std::vector<l1t::Tau>();
  std::vector<l1t::Tau> *sortedTaus = new std::vector<l1t::Tau>();
  std::vector<l1t::Tau> *preGtIsoTaus = new std::vector<l1t::Tau>();
  std::vector<l1t::Tau> *preSortIsoTaus = new std::vector<l1t::Tau>();
  std::vector<l1t::Tau> *sortedIsoTaus = new std::vector<l1t::Tau>();

  for(CaloRegionBxCollection::const_iterator region = subRegions->begin();
      region != subRegions->end(); region++) {
    if(region->hwEta() < 4 || region->hwEta() > 17)
      continue;

    int regionEt = region->hwPt();
    if(regionEt < tauSeedThreshold) continue;

    int regionEta = region->hwEta();
    int regionPhi = region->hwPhi();

    int tauEt = regionEt;
    int isoFlag = 0;  // is 1 if it passes the relative jet iso requirement
    int quality = 0;  //doesn't really mean anything and isn't used

    int highestNeighborEt=-1;
    int highestNeighborEta=-1;
    int highestNeighborPhi=-1;
    int highestNeighborTauVeto=-1;

    //Find neighbor with highest Et
    for(CaloRegionBxCollection::const_iterator neighbor = subRegions->begin();
    neighbor != subRegions->end(); neighbor++) {

      int neighborPhi = neighbor->hwPhi();
      int neighborEta = neighbor->hwEta();
      if(neighborEta < 4 || neighborEta > 17)
    continue;

      int deltaPhi = regionPhi - neighborPhi;
      if (std::abs(deltaPhi) == L1CaloRegionDetId::N_PHI-1)
    deltaPhi = -deltaPhi/std::abs(deltaPhi); //18 regions in phi

      deltaPhi = std::abs(deltaPhi);
      int deltaEta = std::abs(regionEta - neighborEta);

      if (deltaPhi + deltaEta > 0 && deltaPhi + deltaEta < 2) {  //nondiagonal neighbors
    if (neighbor->hwPt() > highestNeighborEt) {
      highestNeighborEt = neighbor->hwPt();
      highestNeighborEta = neighbor->hwEta();
      highestNeighborPhi = neighbor->hwPhi();
      highestNeighborTauVeto = neighbor->hwQual() & 0x1; // tauVeto should be the first bit of quality integer
    }
      }
    }

    string NESW = findNESW(regionEta, regionPhi, highestNeighborEta, highestNeighborPhi);

    if((tauEt > highestNeighborEt && (NESW=="isEast" || NESW=="isNorth"))
       || (tauEt >= highestNeighborEt && (NESW=="isSouth" || NESW=="isWest"))
       || highestNeighborEt == 0 ) {

      if (highestNeighborEt >= tauNeighbourThreshold) tauEt += highestNeighborEt;

      int regionTauVeto = region->hwQual() & 0x1;  // tauVeto should be the first bit of quality integer

      // compute relative jet isolation
      if (region->hwEta() >= isoTauEtaMin && region->hwEta() <= isoTauEtaMax ){
    if ((highestNeighborTauVeto == 0 && regionTauVeto == 0) || tauEt > tauMaxPtTauVeto) {
      int jetEt=AssociatedJetPt(region->hwEta(), region->hwPhi(),unCorrJets);
      if (jetEt>0){
        unsigned int MAX_LUT_ADDRESS = params_->tauIsolationLUT()->maxSize()-1;
        // unsigned int lutAddress = isoLutIndex(tauEt,jetEt);
        unsigned lutAddress = isoTauLut->lutAddress(tauEt,jetEt);
        if (tauEt >0){
          if (lutAddress > MAX_LUT_ADDRESS) lutAddress = MAX_LUT_ADDRESS;
          isoFlag = params_->tauIsolationLUT()->data(lutAddress);
          // isoFlag= isoTauLut->lutPayload(lutAddress);
          // if (isoFlag != params_->tauIsolationLUT()->data(lutAddress)) std::cout << "XXX -- isoFlag: " << isoFlag << "\tisoFlag2: " <<  params_->tauIsolationLUT()->data(lutAddress) << std::endl;
        }

      }else{ // no associated jet
        isoFlag=1;
      }
    }
      }

      ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > tauLorentz(0,0,0,0);

      l1t::Tau theTau(*&tauLorentz, tauEt, region->hwEta(), region->hwPhi(), quality, isoFlag);

      preGtTaus->push_back(theTau);
      if(isoFlag)
    preGtIsoTaus->push_back(theTau);
    }
  }

  // TauToGtPtScales(params_, preGtTaus, preRankTaus);
  // TauToGtPtScales(params_, preGtIsoTaus, preRankIsoTaus);
  calibrateAndRankTaus(params_,preGtTaus,preSortTaus);
  calibrateAndRankTaus(params_,preGtIsoTaus, preSortIsoTaus);

  SortTaus(preSortTaus, sortedTaus);
  SortTaus(preSortIsoTaus, sortedIsoTaus);

  TauToGtEtaScales(params_, sortedTaus, taus);
  TauToGtEtaScales(params_, sortedIsoTaus, isoTaus);

  // set all filler taus to have isolation bit set
  for(std::vector<l1t::Tau>::iterator iTau = isoTaus->begin(); iTau != isoTaus->end(); ++iTau)
    iTau->setHwIso(1);

  delete subRegions;
  delete unCorrJets;
  delete preGtTaus;
  delete preSortTaus;
  delete sortedTaus;
  delete preGtIsoTaus;
  delete preSortIsoTaus;
  delete sortedIsoTaus;

  const bool verbose = false;
  if(verbose)
  {
    std::cout << "Taus" << std::endl;
    for(std::vector<l1t::Tau>::const_iterator iTau = taus->begin(); iTau != taus->end(); ++iTau)
    {
      unsigned int packed = pack15bits(iTau->hwPt(), iTau->hwEta(), iTau->hwPhi());
      std::cout << bitset<15>(packed).to_string() << std::endl;
    }
    std::cout << "Isolated Taus" << std::endl;
    for(std::vector<l1t::Tau>::const_iterator iTau = isoTaus->begin(); iTau != isoTaus->end(); ++iTau)
    {
      unsigned int packed = pack15bits(iTau->hwPt(), iTau->hwEta(), iTau->hwPhi());
      std::cout << bitset<15>(packed).to_string() << std::endl;
    }
  }
}

Member Data Documentation

Stage1TauIsolationLUT* l1t::Stage1Layer2TauAlgorithmImpHW::isoTauLut

private

Definition at line 66 of file Stage1Layer2TauAlgorithmImp.h.

Referenced by Stage1Layer2TauAlgorithmImpHW(), and ~Stage1Layer2TauAlgorithmImpHW().

CaloParamsHelper* const l1t::Stage1Layer2TauAlgorithmImpHW::params_

private

Definition at line 65 of file Stage1Layer2TauAlgorithmImp.h.

Referenced by Stage1Layer2TauAlgorithmImpHW().