l1t::Stage1Layer2DiTauAlgorithm Class Reference

#include <Stage1Layer2HFRingSumAlgorithmImp.h>

Inheritance diagram for l1t::Stage1Layer2DiTauAlgorithm:

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

Private Attributes
CaloParamsHelper *	params_

Detailed Description

Authors

: R. Alex Barbieri

Description: DiTau Algorithm

Definition at line 56 of file Stage1Layer2HFRingSumAlgorithmImp.h.

Constructor & Destructor Documentation

l1t::Stage1Layer2DiTauAlgorithm::Stage1Layer2DiTauAlgorithm

(

CaloParamsHelper *

params

)

Definition at line 15 of file Stage1Layer2DiTauAlgorithm.cc.

                                                                                  : params_(params)
{
}

l1t::Stage1Layer2DiTauAlgorithm::~Stage1Layer2DiTauAlgorithm ( )

virtual

Definition at line 20 of file Stage1Layer2DiTauAlgorithm.cc.

{}

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  )

virtual

Implements l1t::Stage1Layer2HFRingSumAlgorithm.

Definition at line 23 of file Stage1Layer2DiTauAlgorithm.cc.

References gather_cfg::cout, l1t::L1Candidate::hwPt(), L1GctHFRingEtSums::raw(), alignCSCRings::s, L1GctHFRingEtSums::setEtSum(), l1t::L1Candidate::setHwPt(), and cond::to_string().

                                                  {

  std::vector<l1t::Tau> *isoTaus = new std::vector<l1t::Tau>();
  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->size()>0) {
    isoPtMax= (*isoTaus).at(0).hwPt();
    if (isoTaus->size()>1) diIsoPtMax  = (*isoTaus).at(1).hwPt();
    if (isoTaus->size()>2) triIsoPtMax = (*isoTaus).at(2).hwPt();
    if (isoTaus->size()>3) quadIsoPtMax= (*isoTaus).at(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);

  delete isoTaus;


  const bool verbose = false;
  if(verbose)
  {
    std::cout << "HF Ring Sums (Isolated Taus)" << std::endl;
    std::cout << bitset<12>(spares->hwPt()).to_string() << std::endl;
  }
}

Member Data Documentation

CaloParamsHelper* l1t::Stage1Layer2DiTauAlgorithm::params_

private

Definition at line 65 of file Stage1Layer2HFRingSumAlgorithmImp.h.