l1t::Stage1Layer2EGammaAlgorithmImpHW Class Reference

#include <Stage1Layer2EGammaAlgorithmImp.h>

Inheritance diagram for l1t::Stage1Layer2EGammaAlgorithmImpHW:

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

Private Member Functions
int	AssociatedJetPt (int ieta, int iphi, const std::vector< l1t::Jet > *jets) const
unsigned	isoLutIndex (unsigned int etaPt, unsigned int jetPt) const

Private Attributes
CaloParamsHelper *const	params_

Detailed Description

Definition at line 66 of file Stage1Layer2EGammaAlgorithmImp.h.

Constructor & Destructor Documentation

Stage1Layer2EGammaAlgorithmImpHW::Stage1Layer2EGammaAlgorithmImpHW

(

CaloParamsHelper *

params

)

Definition at line 25 of file Stage1Layer2EGammaAlgorithmImpHW.cc.

: params_(params) {};

Stage1Layer2EGammaAlgorithmImpHW::~Stage1Layer2EGammaAlgorithmImpHW ( )

virtual

Definition at line 27 of file Stage1Layer2EGammaAlgorithmImpHW.cc.

{};

Member Function Documentation

int l1t::Stage1Layer2EGammaAlgorithmImpHW::AssociatedJetPt ( int  ieta, int  iphi, const std::vector< l1t::Jet > *  jets  ) const

private

Definition at line 154 of file Stage1Layer2EGammaAlgorithmImpHW.cc.

References gather_cfg::cout, Debug, reco::btau::jetEta, reco::btau::jetPhi, and EnergyCorrector::pt.

                                                                         {

  bool Debug=false;

  if (Debug) cout << "Number of jets: " << jets->size() << endl;
  int pt = 0;


  for(JetBxCollection::const_iterator itJet = jets->begin();
      itJet != jets->end(); ++itJet){

    int jetEta = itJet->hwEta();
    int jetPhi = itJet->hwPhi();
    if (Debug) cout << "Matching ETA: " << ieta << " " << jetEta << endl;
    if (Debug) cout << "Matching PHI: " << iphi << " " << jetPhi << endl;
    if ((jetEta == ieta) && (jetPhi == iphi)){
      pt = itJet->hwPt();
      break;
    }
  }

  // set output
  return pt;
}

unsigned l1t::Stage1Layer2EGammaAlgorithmImpHW::isoLutIndex ( unsigned int  etaPt, unsigned int  jetPt  ) const

private

Definition at line 142 of file Stage1Layer2EGammaAlgorithmImpHW.cc.

{
  const unsigned int nbitsEG=6;  // number of bits used for EG bins in LUT file (needed for left shift operation)
  //  const unsigned int nbitsJet=9; // not used but here for info  number of bits used for Jet bins in LUT file

  //jetPt &= 511; // Take only the LSB 9 bits to match firmware.
  if(jetPt > 511) jetPt = 511;
  unsigned int address= (jetPt << nbitsEG) + egPt;
  // std::cout << address << "\t## " << egPt << " " << jetPt << std::endl;
  return address;
}

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

virtual

Implements l1t::Stage1Layer2EGammaAlgorithm.

Definition at line 31 of file Stage1Layer2EGammaAlgorithmImpHW.cc.

References l1t::EGammaToGtScales(), cmsHarvester::index, l1t::pack15bits(), l1t::RegionCorrection(), l1t::SortEGammas(), cond::to_string(), and l1t::TwelveByTwelveFinder().

                                                                                                                                                                                                        {


  std::vector<l1t::CaloRegion> *subRegions = new std::vector<l1t::CaloRegion>();
  std::vector<l1t::EGamma> *preSortEGammas = new std::vector<l1t::EGamma>();
  std::vector<l1t::EGamma> *preGtEGammas = new std::vector<l1t::EGamma>();


  //Region Correction will return uncorrected subregions if
  //regionPUSType is set to None in the config
  RegionCorrection(regions, subRegions, params_);

  // ----- need to cluster jets in order to compute jet isolation ----
  std::vector<l1t::Jet> *unCorrJets = new std::vector<l1t::Jet>();
  TwelveByTwelveFinder(0, subRegions, unCorrJets);


  //std::cout << "Begin egamma cand debug dump:" << std::endl;
  //std::cout << "pt eta phi index assocJetPt lutAddress lutResult" << std::endl;
  for(CaloEmCandBxCollection::const_iterator egCand = EMCands.begin();
      egCand != EMCands.end(); egCand++) {

    int eg_et = egCand->hwPt();
    int eg_eta = egCand->hwEta();
    int eg_phi = egCand->hwPhi();
    int index = (egCand->hwIso()*4 + egCand->hwQual()) ;

    //std::cout << "JetRankMax: " << params_->jetScale().rankScaleMax()<< " EmRankMax: " << params_->emScale().rankScaleMax()<< std::endl;
    //std::cout << "JetLinMax: " << params_->jetScale().linScaleMax()<< " EmLinMax: " << params_->emScale().linScaleMax()<< std::endl;

    ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double> > egLorentz(0,0,0,0);

    int isoFlag = 0;
    int ijet_pt=AssociatedJetPt(eg_eta,eg_phi,unCorrJets);
    bool isinBarrel = (eg_eta>=7 && eg_eta<=14);
    unsigned int lutAddress = isoLutIndex(eg_et,ijet_pt);

    // Combined Barrel/Endcap LUT uses upper bit to indicate Barrel / Endcap:
    enum {MAX_LUT_ADDRESS = 0x7fff};
    enum {LUT_BARREL_OFFSET = 0x0, LUT_ENDCAP_OFFSET = 0x8000};
    enum {LUT_RCT_OFFSET = 0x10000};

    unsigned int rct_offset=0;
    if (egCand->hwIso()) rct_offset=LUT_RCT_OFFSET;

    if (eg_et >0){
      if (lutAddress > MAX_LUT_ADDRESS) lutAddress = MAX_LUT_ADDRESS;

      if (isinBarrel){
        isoFlag= params_->egIsolationLUT()->data(LUT_BARREL_OFFSET + rct_offset + lutAddress);
      } else{
        isoFlag= params_->egIsolationLUT()->data(LUT_ENDCAP_OFFSET + rct_offset + lutAddress);
      }
    }

    // std::cout << eg_et << " " << eg_eta << " " << eg_phi << " " << index << " "
    //        << ijet_pt << " " << lutAddress << " " << isoFlag << std::endl;
    l1t::EGamma theEG(*&egLorentz, eg_et, eg_eta, eg_phi, index, isoFlag);
    preSortEGammas->push_back(theEG);
  }

  // printf("Pre-Sort\n");
  // for(std::vector<l1t::EGamma>::const_iterator itEGamma = preSortEGammas->begin();
  //     itEGamma != preSortEGammas->end(); ++itEGamma){
  //   std::cout << itEGamma->hwPt() << " " << itEGamma->hwEta() << " " << itEGamma->hwPhi() << std::endl;
  // }

  SortEGammas(preSortEGammas, preGtEGammas);

  EGammaToGtScales(params_, preGtEGammas, egammas);

  // printf("Post-Sort\n");
  // for(std::vector<l1t::EGamma>::const_iterator itEGamma = egammas->begin();
  //     itEGamma != egammas->end(); ++itEGamma){
  //   std::cout << itEGamma->hwPt() << " " << itEGamma->hwEta() << " " << itEGamma->hwPhi() << std::endl;
  // }

  const bool verbose = false;
  if(verbose)
  {
    int cEGammas = 0;
    int fEGammas = 0;
    printf("EGammas Isolated\n");
    for(std::vector<l1t::EGamma>::const_iterator itEGamma = egammas->begin();
    itEGamma != egammas->end(); ++itEGamma){
      if(itEGamma->hwIso() != 1) continue;
      cEGammas++;
      unsigned int packed = pack15bits(itEGamma->hwPt(), itEGamma->hwEta(), itEGamma->hwPhi());
      cout << bitset<15>(packed).to_string() << endl;
      if(cEGammas == 4) break;
    }

    printf("EGammas Non-isolated\n");
    for(std::vector<l1t::EGamma>::const_iterator itEGamma = egammas->begin();
    itEGamma != egammas->end(); ++itEGamma){
      if(itEGamma->hwIso() != 0) continue;
      fEGammas++;
      unsigned int packed = pack15bits(itEGamma->hwPt(), itEGamma->hwEta(), itEGamma->hwPhi());
      cout << bitset<15>(packed).to_string() << endl;
      if(fEGammas == 4) break;
    }
  }

  delete subRegions;
  delete unCorrJets;
  delete preSortEGammas;
  delete preGtEGammas;

}

Member Data Documentation

CaloParamsHelper* const l1t::Stage1Layer2EGammaAlgorithmImpHW::params_

private

Definition at line 75 of file Stage1Layer2EGammaAlgorithmImp.h.