l1t::Stage1Layer2JetAlgorithmImpHI Class Reference

#include <Stage1Layer2JetAlgorithmImp.h>

Inheritance diagram for l1t::Stage1Layer2JetAlgorithmImpHI:

Public Member Functions
virtual void	processEvent (const std::vector< l1t::CaloRegion > &regions, const std::vector< l1t::CaloEmCand > &EMCands, std::vector< l1t::Jet > *jets, std::vector< l1t::Jet > *preGtJets)
	Stage1Layer2JetAlgorithmImpHI (CaloParamsHelper *params)
virtual	~Stage1Layer2JetAlgorithmImpHI ()
Public Member Functions inherited from l1t::Stage1Layer2JetAlgorithm
virtual	~Stage1Layer2JetAlgorithm ()

Private Attributes
CaloParamsHelper *const	params_

Detailed Description

Author

: R. Alex Barbieri MIT

Definition at line 27 of file Stage1Layer2JetAlgorithmImp.h.

Constructor & Destructor Documentation

Stage1Layer2JetAlgorithmImpHI::Stage1Layer2JetAlgorithmImpHI

(

CaloParamsHelper *

params

)

Definition at line 19 of file Stage1Layer2JetAlgorithmImpHI.cc.

: params_(params) { };

Stage1Layer2JetAlgorithmImpHI::~Stage1Layer2JetAlgorithmImpHI ( )

virtual

Definition at line 21 of file Stage1Layer2JetAlgorithmImpHI.cc.

{};

Member Function Documentation

void Stage1Layer2JetAlgorithmImpHI::processEvent ( const std::vector< l1t::CaloRegion > &  regions, const std::vector< l1t::CaloEmCand > &  EMCands, std::vector< l1t::Jet > *  jets, std::vector< l1t::Jet > *  preGtJets  )

virtual

Implements l1t::Stage1Layer2JetAlgorithm.

Definition at line 25 of file Stage1Layer2JetAlgorithmImpHI.cc.

References gather_cfg::cout, l1t::HICaloRingSubtraction(), l1t::CaloParamsHelper::jetRegionMask(), l1t::CaloParamsHelper::jetSeedThreshold(), l1t::JetToGtEtaScales(), l1t::JetToGtPtScales(), l1t::pack15bits(), l1t::pack16bits(), params_, l1t::SortJets(), cond::to_string(), and l1t::TwoByTwoFinder().

                                                        {

  //std::vector<double> regionPUSParams = params_->regionPUSParams();
  int jetThreshold = params_->jetSeedThreshold();

  unsigned int etaMask = params_->jetRegionMask();

  std::vector<l1t::CaloRegion> *subRegions = new std::vector<l1t::CaloRegion>();
  std::vector<l1t::Jet> *unSortedJets = new std::vector<l1t::Jet>();
  std::vector<l1t::Jet> *preGtEtaJets = new std::vector<l1t::Jet>();
  std::vector<l1t::Jet> *preRankJets = new std::vector<l1t::Jet>();

  HICaloRingSubtraction(regions, subRegions, params_);
  TwoByTwoFinder(jetThreshold, etaMask, subRegions, preRankJets);
  //slidingWindowJetFinder(0, subRegions, unSortedJets);
  JetToGtPtScales(params_, preRankJets, unSortedJets);
  //verboseDumpJets(*unSortedJets);
  SortJets(unSortedJets, preGtEtaJets);
  JetToGtEtaScales(params_, preGtEtaJets, preGtJets);
  JetToGtEtaScales(params_, preGtEtaJets, jets);
  //JetToGtPtScales(params_, preGtJets, jets);

  delete subRegions;
  delete unSortedJets;
  delete preGtEtaJets;

  const bool verbose = false;
  const bool hex = true;
  if(verbose)
  {
    if(!hex)
    {
      int cJets = 0;
      int fJets = 0;
      printf("Jets Central\n");
      //printf("pt\teta\tphi\n");
      for(std::vector<l1t::Jet>::const_iterator itJet = jets->begin();
      itJet != jets->end(); ++itJet){
    if((itJet->hwQual() & 2) == 2) continue;
    cJets++;
    unsigned int packed = pack15bits(itJet->hwPt(), itJet->hwEta(), itJet->hwPhi());
    cout << bitset<15>(packed).to_string() << endl;
    if(cJets == 4) break;
      }

      printf("Jets Forward\n");
      //printf("pt\teta\tphi\n");
      for(std::vector<l1t::Jet>::const_iterator itJet = jets->begin();
      itJet != jets->end(); ++itJet){
    if((itJet->hwQual() & 2) != 2) continue;
    fJets++;
    unsigned int packed = pack15bits(itJet->hwPt(), itJet->hwEta(), itJet->hwPhi());
    cout << bitset<15>(packed).to_string() << endl;

    if(fJets == 4) break;
      }
    } else {
      std::cout << "Jets" << std::endl;
      l1t::Jet ajets[8];
      for(std::vector<l1t::Jet>::const_iterator itJet = jets->begin();
      itJet != jets->end(); ++itJet){
    ajets[itJet - jets->begin()] = *itJet;
      }
      //std::cout << "Jets (hex)" << std::endl;
      std::cout << std::hex << pack16bits(ajets[0].hwPt(), ajets[0].hwEta(), ajets[0].hwPhi());
      std::cout << " ";
      std::cout << std::hex << pack16bits(ajets[1].hwPt(), ajets[1].hwEta(), ajets[1].hwPhi());
      std::cout << " ";
      std::cout << std::hex << pack16bits(ajets[4].hwPt(), ajets[4].hwEta(), ajets[4].hwPhi());
      std::cout << " ";
      std::cout << std::hex << pack16bits(ajets[5].hwPt(), ajets[5].hwEta(), ajets[5].hwPhi());
      std::cout << std::endl;
      std::cout << std::hex << pack16bits(ajets[2].hwPt(), ajets[2].hwEta(), ajets[2].hwPhi());
      std::cout << " ";
      std::cout << std::hex << pack16bits(ajets[3].hwPt(), ajets[3].hwEta(), ajets[3].hwPhi());
      std::cout << " ";
      std::cout << std::hex << pack16bits(ajets[6].hwPt(), ajets[6].hwEta(), ajets[6].hwPhi());
      std::cout << " ";
      std::cout << std::hex << pack16bits(ajets[7].hwPt(), ajets[7].hwEta(), ajets[7].hwPhi());
      std::cout << std::endl;
    }
  }
}

Member Data Documentation

CaloParamsHelper* const l1t::Stage1Layer2JetAlgorithmImpHI::params_

private

Definition at line 36 of file Stage1Layer2JetAlgorithmImp.h.

Referenced by processEvent().