#include <Stage1TauIsolationLUT.h>

Public Member Functions
unsigned	lutAddress (unsigned int, unsigned int) const

int	lutPayload (unsigned int) const

	Stage1TauIsolationLUT (CaloParamsHelper *params)

virtual	~Stage1TauIsolationLUT ()

Static Public Attributes
static const unsigned int	lut_version = 1

static const unsigned int	nbits_data = 1

static const unsigned int	nbitsJet = 8

static const unsigned int	nbitsTau = 8

Private Attributes
CaloParamsHelper *const	params_

Detailed Description

Definition at line 19 of file Stage1TauIsolationLUT.h.

Constructor & Destructor Documentation

Stage1TauIsolationLUT::Stage1TauIsolationLUT ( CaloParamsHelper * params )

Definition at line 15 of file Stage1TauIsolationLUT.cc.

15 : params_(params)

16 {};

l1t::Stage1TauIsolationLUT::params_

CaloParamsHelper *const params_

Definition: Stage1TauIsolationLUT.h:33

Stage1TauIsolationLUT::~Stage1TauIsolationLUT ( )

virtual

Definition at line 84 of file Stage1TauIsolationLUT.cc.

84 {};

Member Function Documentation

unsigned Stage1TauIsolationLUT::lutAddress	(	unsigned int	tauPt,
		unsigned int	jetPt
	)		const

Definition at line 18 of file Stage1TauIsolationLUT.cc.

References funct::abs(), gather_cfg::cout, l1t::CaloParamsHelper::jetLsb(), nbitsJet, nbitsTau, params_, and funct::pow().

 {
   const unsigned int maxJet = pow(2,nbitsJet)-1;
   const unsigned int maxTau = pow(2,nbitsTau)-1;
   
   // lut only defined for 8 bit Jet ET
   if ( (nbitsJet != 8) ||  (nbitsTau != 8) ) return 0;
 
   double jetLsb=params_->jetLsb();
   if (std::abs(jetLsb-0.5) > 0.0001){
     std::cout << "%Stage1TauIsolationLUT-E-Unexpected jetLsb " << jetLsb << " IsoTau calculation will be broken"<< std::endl;
     return 0;
   }
 
   tauPt=tauPt>>1;
   jetPt=jetPt>>1;
 
   if (jetPt>maxJet) jetPt=maxJet;
   if (tauPt>maxTau) tauPt=maxTau;
 
   unsigned int address= (jetPt << nbitsTau) + tauPt;
   return address;
 }

int Stage1TauIsolationLUT::lutPayload ( unsigned int address ) const

Definition at line 42 of file Stage1TauIsolationLUT.cc.

References gather_cfg::cout, nbitsJet, nbitsTau, params_, funct::pow(), l1t::CaloParamsHelper::tauMaxJetIsolationA(), l1t::CaloParamsHelper::tauMaxJetIsolationB(), and l1t::CaloParamsHelper::tauMinPtJetIsolationB().

Referenced by l1t::L1TCaloStage1LutWriter::writeIsoTauLut().

 {
 
   const unsigned int maxJet = pow(2,nbitsJet)-1;
   const unsigned int maxTau = pow(2,nbitsTau)-1;
 
   const double tauMaxJetIsolationA = params_->tauMaxJetIsolationA();
   const double tauMaxJetIsolationB = params_->tauMaxJetIsolationB();
   const double tauMinPtJetIsolationB = params_->tauMinPtJetIsolationB();
   
   unsigned int maxAddress = pow(2,nbitsJet+nbitsTau)-1;    
   if (address > maxAddress){  // check that address is right length
     std::cout << "%Stage1TauIsolationLUT-E-Address: " << address 
           << " exceeds maximum value allowed. Setting value to maximum (" << maxAddress << ")" << std::endl;
     address = maxAddress;
   }
   // extract the jet and tau et from the address
   int ijetet = address >> nbitsTau;
   int itauet = address & 0xff;
 
   
   double jet_pt = static_cast <float> (ijetet);  // no need convert to physical eT, as right shift (>>1) operation
   double tau_pt = static_cast <float> (itauet);  // in lutAddress automatically converts to physical eT, assuming lsb=0.5
 
   //std::cout << "ijetet: " << ijetet << "\titauet: " << itauet << std::endl;
   //std::cout << "jetet: " << jet_pt << "\ttauet: " << tau_pt << std::endl;
 
   int isol=0;
   if (maxTau == (unsigned) itauet){
     isol=1;
   } else if (maxJet == (unsigned) ijetet){
     isol=1;
   } else {
     double relativeJetIsolationTau = (jet_pt / tau_pt) -1;
         
     double isolCut=tauMaxJetIsolationA;
     if (tau_pt >= tauMinPtJetIsolationB)isolCut=tauMaxJetIsolationB;
     if (relativeJetIsolationTau < isolCut) isol=1;
   }
   return isol;
 }