l1t::Stage1Layer2DiTauAlgorithm Class Reference

#include <Stage1Layer2HFRingSumAlgorithmImp.h>

Inheritance diagram for l1t::Stage1Layer2DiTauAlgorithm:

Public Member Functions
void	processEvent (const std::vector< l1t::CaloRegion > &regions, const std::vector< l1t::CaloEmCand > &EMCands, const std::vector< l1t::Tau > *taus, l1t::CaloSpare *spare) override
	Stage1Layer2DiTauAlgorithm (CaloParamsHelper const *params)
	~Stage1Layer2DiTauAlgorithm () override=default
Public Member Functions inherited from l1t::Stage1Layer2HFRingSumAlgorithm
virtual	~Stage1Layer2HFRingSumAlgorithm ()

Private Attributes
CaloParamsHelper const *	params_

Detailed Description

Authors

: R. Alex Barbieri

Description: DiTau Algorithm

Definition at line 55 of file Stage1Layer2HFRingSumAlgorithmImp.h.

Constructor & Destructor Documentation

l1t::Stage1Layer2DiTauAlgorithm::Stage1Layer2DiTauAlgorithm

(

CaloParamsHelper const *

params

)

Definition at line 15 of file Stage1Layer2DiTauAlgorithm.cc.

                                                                                        : params_(params)
{
}

l1t::Stage1Layer2DiTauAlgorithm::~Stage1Layer2DiTauAlgorithm ( )

overridedefault

Member Function Documentation

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

overridevirtual

Implements l1t::Stage1Layer2HFRingSumAlgorithm.

Definition at line 20 of file Stage1Layer2DiTauAlgorithm.cc.

References l1t::LUT::data(), L1CaloEtScale::et(), l1t::CaloParamsHelper::HfRingScale(), l1t::CaloParamsHelper::jetScale(), l1t::LUT::maxSize(), params_, L1CaloEtScale::rank(), L1GctHFRingEtSums::raw(), alignCSCRings::s, L1GctHFRingEtSums::setEtSum(), l1t::L1Candidate::setHwPt(), and l1t::CaloParamsHelper::tauEtToHFRingEtLUT().

                                                  {

  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);
}

Member Data Documentation

CaloParamsHelper const* l1t::Stage1Layer2DiTauAlgorithm::params_

private

Definition at line 64 of file Stage1Layer2HFRingSumAlgorithmImp.h.

Referenced by processEvent().