Stage1Layer2SingleTrackHI.cc

Go to the documentation of this file.

// Stage1Layer2SingleTrackHI.cc
// Authors: Michael Northup
//          Alex Barbieri
//
// This is a special-purpose single-track seed trigger which uses the
// Tau channel to communicate with GT. Be wary of any naming scheme
// because we are masquerading as both a tau and track trigger when
// we are really just looking for the hottest region.


#include "L1Trigger/L1TCalorimeter/interface/Stage1Layer2TauAlgorithmImp.h"
#include "DataFormats/L1CaloTrigger/interface/L1CaloRegionDetId.h"

#include "L1Trigger/L1TCalorimeter/interface/PUSubtractionMethods.h"
#include "L1Trigger/L1TCalorimeter/interface/legacyGtHelper.h"
#include "L1Trigger/L1TCalorimeter/interface/HardwareSortingMethods.h"


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

void findRegions(const std::vector<l1t::CaloRegion> * sr, std::vector<l1t::Tau> * t, const int etaMask);

void l1t::Stage1Layer2SingleTrackHI::processEvent(const std::vector<l1t::CaloEmCand> & clusters,
                          const std::vector<l1t::CaloRegion> & regions,
                          std::vector<l1t::Tau> * isoTaus,
                          std::vector<l1t::Tau> * taus)
{
  int etaMask = params_->tauRegionMask();

  std::vector<l1t::CaloRegion> subRegions;
  std::vector<l1t::Tau> preGtEtaTaus;
  std::vector<l1t::Tau> preGtTaus;
  std::vector<l1t::Tau> unsortedTaus;


  HICaloRingSubtraction(regions, &subRegions, params_);
  findRegions(&subRegions, &preGtTaus, etaMask);
  TauToGtPtScales(params_, &preGtTaus, &unsortedTaus);
  SortTaus(&unsortedTaus, &preGtEtaTaus);
  //SortTaus(preGtTaus, unsortedTaus);
  //TauToGtPtScales(params_, unsortedTaus, preGtEtaTaus);
  TauToGtEtaScales(params_, &preGtEtaTaus, taus);

  isoTaus->resize(4);
  //taus->resize(4);

}

void findRegions(const std::vector<l1t::CaloRegion> * sr, std::vector<l1t::Tau> * t, const int etaMask)
{
  for(std::vector<l1t::CaloRegion>::const_iterator region = sr->begin(); region != sr->end(); region++)
  {
    int tauEta = region->hwEta();
    if(tauEta < 4 || tauEta > 17) continue; // taus CANNOT be in the forward region
    if((etaMask & (1<<tauEta))>>tauEta) continue;

    ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > TauLorentz(0,0,0,0);
    l1t::Tau taucand(*&TauLorentz,region->hwPt(),region->hwEta(),region->hwPhi());

    t->push_back(taucand);
  }
}