#include <Stage2Layer2JetSumAlgorithmFirmware.h>

Inheritance diagram for l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1:

Public Member Functions
void	processEvent (const std::vector< l1t::Jet > &alljets, std::vector< l1t::EtSum > &htsums) override

	Stage2Layer2JetSumAlgorithmFirmwareImp1 (CaloParamsHelper const *params)

	~Stage2Layer2JetSumAlgorithmFirmwareImp1 () override=default

Public Member Functions inherited from l1t::Stage2Layer2JetSumAlgorithm
virtual	~Stage2Layer2JetSumAlgorithm ()

Private Attributes
int32_t	httEtaMax_

int32_t	httEtaMaxHF_

int32_t	httJetThresholdHw_

int32_t	mhtEtaMax_

int32_t	mhtEtaMaxHF_

int32_t	mhtJetThresholdHw_

Detailed Description

Author: : Jim Brooke

Description: first iteration of stage 2 jet algo

Definition at line 22 of file Stage2Layer2JetSumAlgorithmFirmware.h.

Constructor & Destructor Documentation

l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::Stage2Layer2JetSumAlgorithmFirmwareImp1 ( CaloParamsHelper const * params )

Definition at line 12 of file Stage2Layer2JetSumAlgorithmFirmwareImp1.cc.

References l1t::CaloParamsHelper::etSumEtaMax(), l1t::CaloParamsHelper::etSumEtThreshold(), httEtaMax_, httEtaMaxHF_, httJetThresholdHw_, l1t::CaloParamsHelper::jetLsb(), l1t::CaloTools::kHFEnd, mhtEtaMax_, mhtEtaMaxHF_, and mhtJetThresholdHw_.

 {
   httJetThresholdHw_ = floor(params->etSumEtThreshold(1)/params->jetLsb());
   mhtJetThresholdHw_ = floor(params->etSumEtThreshold(3)/params->jetLsb());
 
   httEtaMax_  = params->etSumEtaMax(1);
   httEtaMaxHF_ = CaloTools::kHFEnd;
   mhtEtaMax_  = params->etSumEtaMax(3);
   mhtEtaMaxHF_ = CaloTools::kHFEnd;
 }

l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::~Stage2Layer2JetSumAlgorithmFirmwareImp1 ( )

overridedefault

Member Function Documentation

void l1t::Stage2Layer2JetSumAlgorithmFirmwareImp1::processEvent	(	const std::vector< l1t::Jet > &	alljets,
		std::vector< l1t::EtSum > &	htsums
	)

overridevirtual

Implements l1t::Stage2Layer2JetSumAlgorithm.

Definition at line 24 of file Stage2Layer2JetSumAlgorithmFirmwareImp1.cc.

References funct::abs(), l1t::CaloTools::cos_coeff, httEtaMax_, httEtaMaxHF_, httJetThresholdHw_, l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPt(), createfilelist::int, l1t::CaloTools::kHBHENrPhi, l1t::CaloTools::kHFEnd, l1t::CaloTools::kSatJet, L1Analysis::kTotalHt, L1Analysis::kTotalHtHF, L1Analysis::kTotalHtx, L1Analysis::kTotalHtxHF, L1Analysis::kTotalHty, L1Analysis::kTotalHtyHF, mhtEtaMax_, mhtEtaMaxHF_, mhtJetThresholdHw_, l1t::CaloTools::mpEta(), p4, and l1t::CaloTools::sin_coeff.

 {
 
   // etaSide=1 is positive eta, etaSide=-1 is negative eta
   for (int etaSide=1; etaSide>=-1; etaSide-=2) {
 
     int hx(0), hy(0), ht(0);
     int hxHF(0), hyHF(0), htHF(0);
 
     bool satMht(false), satMhtHF(false), satHt(false), satHtHF(false);
   
     // loop over rings    
     for (unsigned absieta=1; absieta<=(unsigned int)CaloTools::mpEta(CaloTools::kHFEnd); absieta++) {
 
       int ieta = etaSide * absieta;
 
       int ringHx(0), ringHy(0), ringHt(0); 
       int ringHxHF(0), ringHyHF(0), ringHtHF(0); 
       
       // loop over phi
       for (int iphi=1; iphi<=CaloTools::kHBHENrPhi; iphi++) {
     
         // find the jet at this (eta,phi)
         l1t::Jet thisJet;
         bool foundJet = false;
         for (unsigned jetIt=0; jetIt<alljets.size(); jetIt++) {
           if (CaloTools::mpEta(alljets.at(jetIt).hwEta())==ieta && alljets.at(jetIt).hwPhi()==iphi) {
             thisJet = alljets.at(jetIt);
             foundJet = true;
           }
         }
         if (!foundJet) continue;
     
           // x- and -y coefficients are truncated by after multiplication of Et by trig coefficient.
           // The trig coefficients themselves take values [-1023,1023] and so were scaled by
           // 2^10 = 1024, which requires bitwise shift to the right of the final value by 10 bits.
           // The 4 below account for part of that and the rest is accounted for at ouput of demux
           // (see Stage2Layer2DemuxSumsAlgoFirmwareImp1.cc)
 
         if (thisJet.hwPt()>mhtJetThresholdHw_ && CaloTools::mpEta(abs(thisJet.hwEta()))<=CaloTools::mpEta(mhtEtaMax_)) {
           if(thisJet.hwPt()==CaloTools::kSatJet){
             satMht=true;
             satMhtHF=true;
           }else{
             ringHx += (int) (( thisJet.hwPt() * CaloTools::cos_coeff[iphi - 1] ) >> 4 );
             ringHy += (int) (( thisJet.hwPt() * CaloTools::sin_coeff[iphi - 1] ) >> 4 );
           }
         }
         if (thisJet.hwPt()>mhtJetThresholdHw_ && CaloTools::mpEta(abs(thisJet.hwEta()))<=CaloTools::mpEta(mhtEtaMaxHF_)) {
           if(thisJet.hwPt()==CaloTools::kSatJet) satMhtHF=true;
           else{
             ringHxHF += (int) (( thisJet.hwPt() * CaloTools::cos_coeff[iphi - 1] ) >> 4 );
             ringHyHF += (int) (( thisJet.hwPt() * CaloTools::sin_coeff[iphi - 1] ) >> 4 );
           }
         }
         
         if (thisJet.hwPt()>httJetThresholdHw_ && CaloTools::mpEta(abs(thisJet.hwEta()))<=CaloTools::mpEta(httEtaMax_)) {
           if(thisJet.hwPt()==CaloTools::kSatJet){
             satHt=true;
             satHtHF=true;
           }
           else ringHt += thisJet.hwPt();
         }
         if (thisJet.hwPt()>httJetThresholdHw_ && CaloTools::mpEta(abs(thisJet.hwEta()))<=CaloTools::mpEta(httEtaMaxHF_)) {
           if(thisJet.hwPt()==CaloTools::kSatJet) satHtHF=true;
           else ringHtHF += thisJet.hwPt();
         }
       }
 
       hx += ringHx;
       hy += ringHy;
       ht += ringHt;
       
       hxHF += ringHxHF;
       hyHF += ringHyHF;
       htHF += ringHtHF;
 
     }
 
     if(satHt) ht = 0xffff;
     if(satHtHF) htHF = 0xffff;
 
     if(satMht){ 
       hx = 0x7fffffff;
       hy = 0x7fffffff;
     }
     if(satMhtHF){
       hxHF = 0x7fffffff;
       hyHF = 0x7fffffff;
     }
     
     math::XYZTLorentzVector p4;
     
     l1t::EtSum htSumHt(p4,l1t::EtSum::EtSumType::kTotalHt,ht,0,0,0);
     l1t::EtSum htSumHx(p4,l1t::EtSum::EtSumType::kTotalHtx,hx,0,0,0);
     l1t::EtSum htSumHy(p4,l1t::EtSum::EtSumType::kTotalHty,hy,0,0,0);
 
     l1t::EtSum htSumHtHF(p4,l1t::EtSum::EtSumType::kTotalHtHF,htHF,0,0,0);
     l1t::EtSum htSumHxHF(p4,l1t::EtSum::EtSumType::kTotalHtxHF,hxHF,0,0,0);
     l1t::EtSum htSumHyHF(p4,l1t::EtSum::EtSumType::kTotalHtyHF,hyHF,0,0,0);
     
     htsums.push_back(htSumHt);
     htsums.push_back(htSumHx);
     htsums.push_back(htSumHy);
      
     htsums.push_back(htSumHtHF);
     htsums.push_back(htSumHxHF);
     htsums.push_back(htSumHyHF);
 
   }
 }