Stage1Layer2DiTauAlgorithm.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "L1Trigger/L1TCalorimeter/interface/Stage1Layer2HFRingSumAlgorithmImp.h"
//#include "L1Trigger/L1TCalorimeter/interface/PUSubtractionMethods.h"
//#include "L1Trigger/L1TCalorimeter/interface/legacyGtHelper.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctHFRingEtSums.h"

l1t::Stage1Layer2DiTauAlgorithm::Stage1Layer2DiTauAlgorithm(CaloParamsHelper const* params) : params_(params)
{
}


void l1t::Stage1Layer2DiTauAlgorithm::processEvent(const std::vector<l1t::CaloRegion> & regions,
                           const std::vector<l1t::CaloEmCand> & EMCands,
                           const std::vector<l1t::Tau> * taus,
                           l1t::CaloSpare * spares) {

  std::vector<l1t::Tau> isoTaus;
  for(std::vector<l1t::Tau>::const_iterator itTau = taus->begin();
      itTau != taus->end(); ++itTau){
    if( !itTau->hwIso() ) continue;
    isoTaus.push_back( *itTau );
  }

  int isoPtMax=0;
  int diIsoPtMax=0;
  int triIsoPtMax=0;
  int quadIsoPtMax=0;
  if(!isoTaus.empty()) {
    isoPtMax= isoTaus[0].hwPt();
    if (isoTaus.size()>1) diIsoPtMax  = isoTaus[1].hwPt();
    if (isoTaus.size()>2) triIsoPtMax = isoTaus[2].hwPt();
    if (isoTaus.size()>3) quadIsoPtMax= isoTaus[3].hwPt();
  }

  int rankIso     = 0;
  int rankDiIso   = 0;
  int rankTriIso  = 0;
  int rankQuadIso = 0;

  bool useLut=true;
  // encode the highest pt Iso and DiIso in the HF ET rings
  if (useLut){
    unsigned int MAX_LUT_ADDRESS = params_->tauEtToHFRingEtLUT()->maxSize()-1;
    unsigned int lutAddress = isoPtMax;
    if (lutAddress > MAX_LUT_ADDRESS) lutAddress = MAX_LUT_ADDRESS;
    rankIso=params_->tauEtToHFRingEtLUT()->data(lutAddress);

    lutAddress = diIsoPtMax;
    if (lutAddress > MAX_LUT_ADDRESS) lutAddress = MAX_LUT_ADDRESS;
    rankDiIso=params_->tauEtToHFRingEtLUT()->data(lutAddress);

    lutAddress = triIsoPtMax;
    if (lutAddress > MAX_LUT_ADDRESS) lutAddress = MAX_LUT_ADDRESS;
    rankTriIso=params_->tauEtToHFRingEtLUT()->data(lutAddress);

    lutAddress = quadIsoPtMax;
    if (lutAddress > MAX_LUT_ADDRESS) lutAddress = MAX_LUT_ADDRESS;
    rankQuadIso=params_->tauEtToHFRingEtLUT()->data(lutAddress);

  }else{
    double etIso     = params_->jetScale().et( isoPtMax );  // convert from hwPt to Physical pT
    double etDiIso   = params_->jetScale().et( diIsoPtMax );
    double etTriIso  = params_->jetScale().et( triIsoPtMax );
    double etQuadIso = params_->jetScale().et( quadIsoPtMax );
    rankIso     = params_->HfRingScale().rank( etIso );  //convert to HfRingScale Rank
    rankDiIso   = params_->HfRingScale().rank( etDiIso );
    rankTriIso  = params_->HfRingScale().rank( etTriIso );
    rankQuadIso = params_->HfRingScale().rank( etQuadIso );
  }

  // std::cout << "Max Iso Tau pT: " << isoPtMax << "\t" << etIso << "\t" << rankIso << std::endl;

  L1GctHFRingEtSums s;
  s.setEtSum(0, rankIso);
  s.setEtSum(1, rankDiIso);
  s.setEtSum(2, rankTriIso);
  s.setEtSum(3, rankQuadIso);
  uint16_t raw = s.raw();
  spares->setHwPt(raw);
}