L1MuBMEUX Class Reference

#include <L1MuBMEUX.h>

Inheritance diagram for L1MuBMEUX:

Classes
class	EUX_Comp
	helper class for finding the best and second best extrapolation More...

Public Member Functions
unsigned short int	address () const
	return extrapolation address (0-11) (address = 15 indicates negative ext. result) More...
int	id () const
	return Extrapolation identifier More...
	L1MuBMEUX (const L1MuBMSectorProcessor &sp, const L1MuBMSEU &seu, int id)
	constructor More...
void	load (const L1MuBMTrackSegPhi *start_ts, const L1MuBMTrackSegPhi *target_ts)
	load data into EUX More...
bool	operator== (const L1MuBMEUX &) const
	equal operator More...
unsigned int	quality () const
	return extrapolation quality More...
void	reset () override
	reset Extrapolator More...
bool	result () const
	return extrapolation result More...
void	run (const edm::EventSetup &c) override
	run Extrapolator More...
std::pair< const L1MuBMTrackSegPhi *, const L1MuBMTrackSegPhi * >	ts () const
	return pointer to start and target track segment More...
	~L1MuBMEUX () override
	destructor More...
Public Member Functions inherited from L1AbstractProcessor
virtual void	run ()
	run processor logic More...
virtual void	run ()
	run processor logic More...
virtual	~L1AbstractProcessor ()
	destructor More...
virtual	~L1AbstractProcessor ()
	destructor More...

Private Member Functions
int	sec_mod (int) const
	output sector numbers in the range -6 to +5 More...

Private Attributes
edm::ESHandle< L1TMuonBarrelParams >	bmtfParamsHandle
unsigned short int	m_address
int	m_id
unsigned short int	m_quality
bool	m_result
const L1MuBMSEU &	m_seu
const L1MuBMSectorProcessor &	m_sp
const L1MuBMTrackSegPhi *	m_start
const L1MuBMTrackSegPhi *	m_target
unsigned short const	nbit_phi = 12
unsigned short const	nbit_phib = 10
L1MuDTTFParameters	pars
int const	theExtFilter = 1
L1MuBMLUTHandler *	theExtLUTs

Detailed Description

Extrapolator:

The Extrapolator receives a start and a target track segment and performs the actual extrapolation

N. Neumeister CERN EP

Definition at line 52 of file L1MuBMEUX.h.

Constructor & Destructor Documentation

L1MuBMEUX::L1MuBMEUX	(	const L1MuBMSectorProcessor &	sp,
		const L1MuBMSEU &	seu,
		int	id
	)

constructor

Definition at line 50 of file L1MuBMEUX.cc.

                                                                                  :
    m_sp(sp), m_seu(seu), m_id(id),
    m_result(false), m_quality(0), m_address(15),
    m_start(nullptr), m_target(nullptr),
    theExtFilter(L1MuBMTFConfig::getExtTSFilter()),
    nbit_phi(L1MuBMTFConfig::getNbitsExtPhi()),
    nbit_phib(L1MuBMTFConfig::getNbitsExtPhib())
{
}

L1MuBMEUX::~L1MuBMEUX ( )

override

destructor

Definition at line 65 of file L1MuBMEUX.cc.

{}

Member Function Documentation

unsigned short int L1MuBMEUX::address ( ) const

inline

return extrapolation address (0-11) (address = 15 indicates negative ext. result)

Definition at line 101 of file L1MuBMEUX.h.

References m_address, and sec_mod().

Referenced by operator==().

{ return m_address; }

int L1MuBMEUX::id ( void  ) const

inline

return Extrapolation identifier

Definition at line 92 of file L1MuBMEUX.h.

References m_id.

Referenced by operator==().

{ return m_id; }

void L1MuBMEUX::load	(	const L1MuBMTrackSegPhi *	start_ts,
		const L1MuBMTrackSegPhi *	target_ts
	)

load data into EUX

Definition at line 207 of file L1MuBMEUX.cc.

References EX21, EX34, L1MuBMSEU::ext(), m_seu, m_start, and m_target.

                                                         {

  m_start  = start_ts;
  m_target = target_ts;

  // in case of EX34 and EX21 exchange start and target
  if ( ( m_seu.ext() == EX34 ) || ( m_seu.ext() == EX21 ) ) {
    m_start  = target_ts;
    m_target = start_ts;
  }

}

bool L1MuBMEUX::operator==

(

const L1MuBMEUX &

eux

)

const

equal operator

Definition at line 75 of file L1MuBMEUX.cc.

References address(), id(), m_address, m_id, m_quality, m_result, quality(), and result().

                                                     {

  if ( m_id      != eux.id() )      return false;
  if ( m_result  != eux.result() )  return false;
  if ( m_quality != eux.quality() ) return false;
  if ( m_address != eux.address() ) return false;
  return true;

}

unsigned int L1MuBMEUX::quality ( void  ) const

inline

return extrapolation quality

Definition at line 98 of file L1MuBMEUX.h.

References m_quality.

Referenced by L1MuBMEUX::EUX_Comp::operator()(), and operator==().

{ return m_quality; }

void L1MuBMEUX::reset ( void  )

overridevirtual

reset Extrapolator

Implements L1AbstractProcessor.

Definition at line 226 of file L1MuBMEUX.cc.

References m_address, m_quality, m_result, m_start, and m_target.

                      {

  m_result  = false;
  m_quality = 0;
  m_address = 15;

  m_start  = nullptr;
  m_target = nullptr;

}

bool L1MuBMEUX::result ( ) const

inline

return extrapolation result

Definition at line 95 of file L1MuBMEUX.h.

References m_result.

Referenced by L1MuBMEUX::EUX_Comp::operator()(), and operator==().

{ return m_result; }

void L1MuBMEUX::run ( const edm::EventSetup &  c )

overridevirtual

run Extrapolator

< ext look-up tables

Reimplemented from L1AbstractProcessor.

Definition at line 89 of file L1MuBMEUX.cc.

References funct::abs(), bmtfParamsHandle, gather_cfg::cout, L1MuBMTFConfig::Debug(), diffTreeTool::diff, EX12, EX13, EX14, EX15, EX21, EX23, EX24, EX25, EX34, L1MuBMSEU::ext(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), L1MuDTTFParameters::get_soc_openlut_extr(), L1MuDTTFParameters::get_soc_qcut_st1(), L1MuDTTFParameters::get_soc_qcut_st2(), L1MuDTTFParameters::get_soc_qcut_st4(), L1MuBMLUTHandler::getHigh(), L1MuBMLUTHandler::getLow(), L1MuBMTFConfig::getopenLUTs(), L1MuBMSectorProcessor::id(), L1TMuonBarrelParams::l1mudttfparams, m_address, m_id, m_quality, m_result, m_seu, m_sp, m_start, m_target, nbit_phi, nbit_phib, PFRecoTauDiscriminationByIsolation_cfi::offset, pars, L1MuBMTrackSegPhi::phi(), L1MuBMTrackSegPhi::phib(), edm::ESHandle< T >::product(), L1MuBMTrackSegPhi::quality(), sec_mod(), L1MuBMSecProcId::sector(), L1MuBMTrackSegPhi::sector(), theExtLUTs, L1MuBMSecProcId::wheel(), and L1MuBMTrackSegPhi::wheel().

                                          {

  //c.get< L1MuDTExtLutRcd >().get( theExtLUTs );
//  c.get< L1MuDTTFParametersRcd >().get( pars );

  const L1TMuonBarrelParamsRcd& bmtfParamsRcd = c.get<L1TMuonBarrelParamsRcd>();
  bmtfParamsRcd.get(bmtfParamsHandle);
  const L1TMuonBarrelParams& bmtfParams = *bmtfParamsHandle.product();

  pars =  bmtfParams.l1mudttfparams;
  theExtLUTs =  new L1MuBMLUTHandler(bmtfParams);   


  if ( L1MuBMTFConfig::Debug(4) ) cout << "Run EUX "  << m_id << endl;
  if ( L1MuBMTFConfig::Debug(4) ) cout << "start :  " << *m_start  << endl;
  if ( L1MuBMTFConfig::Debug(4) ) cout << "target : " << *m_target << endl;

  if ( m_start == nullptr || m_target == nullptr ) {
    if ( L1MuBMTFConfig::Debug(4) ) cout << "Error: EUX has no data loaded" << endl;
    delete theExtLUTs;
    return;
  }

  // start sector
  int sector_st = m_start->sector();

  // target sector
  int sector_ta = m_target->sector();

  // get index of look-up table
  int lut_idx = m_seu.ext();
  if ( abs(m_target->wheel()) == 3 ) {

    switch ( m_seu.ext() ) {
      case EX13 : { lut_idx = EX15; break; }
      case EX23 : { lut_idx = EX25; break; }
      default   : { lut_idx = m_seu.ext(); break; }
    }

  }

  if ( L1MuBMTFConfig::Debug(5) ) cout << "EUX : using look-up table : "
                                       << static_cast<Extrapolation>(lut_idx)
                                       << endl;

  // Extrapolation TS quality filter
  int qcut = 0;
  if ( m_seu.ext() == EX12 ) qcut = pars.get_soc_qcut_st1(m_sp.id().wheel(), m_sp.id().sector());
  if ( m_seu.ext() == EX13 ) qcut = pars.get_soc_qcut_st1(m_sp.id().wheel(), m_sp.id().sector());
  if ( m_seu.ext() == EX14 ) qcut = pars.get_soc_qcut_st1(m_sp.id().wheel(), m_sp.id().sector());
  if ( m_seu.ext() == EX21 ) qcut = pars.get_soc_qcut_st2(m_sp.id().wheel(), m_sp.id().sector());
  if ( m_seu.ext() == EX23 ) qcut = pars.get_soc_qcut_st2(m_sp.id().wheel(), m_sp.id().sector());
  if ( m_seu.ext() == EX24 ) qcut = pars.get_soc_qcut_st2(m_sp.id().wheel(), m_sp.id().sector());
  if ( m_seu.ext() == EX34 ) qcut = pars.get_soc_qcut_st4(m_sp.id().wheel(), m_sp.id().sector());

  if ( m_start->quality() < qcut ) {
    delete theExtLUTs;
    return;
  }
  // calculate bit shift
  int sh_phi  = 12 - nbit_phi;
  int sh_phib = 10 - nbit_phib;

  int phi_target = m_target->phi() >> sh_phi;
  int phi_start  = m_start->phi()  >> sh_phi;
  int phib_start = (m_start->phib() >> sh_phib) << sh_phib;
  if ( phib_start < 0 ) phib_start += (1 << sh_phib) -1;

  // compute difference in phi
  int diff = phi_target - phi_start;

  // get low and high values from look-up table
  // and add offset (30 degrees ) for extrapolation to adjacent sector
  int offset = -2144 >> sh_phi;
  offset  *= sec_mod(sector_ta - sector_st);
  int low  = theExtLUTs->getLow(lut_idx,phib_start );
  int high = theExtLUTs->getHigh(lut_idx,phib_start );
  if ( low  < 0 ) low  += (1 << sh_phi) - 1;
  if ( high < 0 ) high += (1 << sh_phi) - 1;
  low  = (low  >> sh_phi) + offset;
  high = (high >> sh_phi) + offset;

  int phi_offset = phi_target - offset;
  if ( ( lut_idx == EX34 ) || ( lut_idx == EX21 ) )
    phi_offset = phi_start + offset;
  if ( phi_offset >= (1 << (nbit_phi-1)) -1 ) {
        delete theExtLUTs;
    return;
  }
  if ( phi_offset < -(1 << (nbit_phi-1)) +1 ) {
    delete theExtLUTs;
    return;
  }
  // is phi-difference within the extrapolation window?
  bool openlut = pars.get_soc_openlut_extr(m_sp.id().wheel(), m_sp.id().sector());
  if (( diff >= low && diff <= high ) || L1MuBMTFConfig::getopenLUTs() || openlut ) {
    m_result = true;
    int qual_st = m_start->quality();
    int qual_ta = m_target->quality();
    if ( m_seu.ext() == EX34 || m_seu.ext() == EX21 ) {
      m_quality = ( qual_st == 7 ) ? 0 : qual_st + 1;
    }
    else {
      m_quality = ( qual_ta == 7 ) ? 0 : qual_ta + 1;
    }
    m_address = m_id;
  }
  delete theExtLUTs;
  if ( L1MuBMTFConfig::Debug(5) ) cout << "diff : "   << low  << " "
                                       << diff << " " << high << " : "
                                       << m_result << " " << endl;

}

int L1MuBMEUX::sec_mod ( int  sector ) const

private

output sector numbers in the range -6 to +5

Definition at line 252 of file L1MuBMEUX.cc.

Referenced by address(), and run().

                                       {

  int new_sector = sector%12;
  if ( new_sector >= 6 )
    new_sector = new_sector - 12;
  if ( new_sector < -6 )
    new_sector = new_sector + 12;

  return new_sector;

}

pair< const L1MuBMTrackSegPhi *,const L1MuBMTrackSegPhi * > L1MuBMEUX::ts

(

)

const

return pointer to start and target track segment

Definition at line 241 of file L1MuBMEUX.cc.

References m_start, and m_target.

                                                                             {

  return pair<const L1MuBMTrackSegPhi*, const L1MuBMTrackSegPhi*>(m_start,m_target);

}

Member Data Documentation

edm::ESHandle< L1TMuonBarrelParams > L1MuBMEUX::bmtfParamsHandle

private

Definition at line 121 of file L1MuBMEUX.h.

Referenced by run().

unsigned short int L1MuBMEUX::m_address

private

Definition at line 116 of file L1MuBMEUX.h.

Referenced by address(), operator==(), reset(), and run().

int L1MuBMEUX::m_id

private

Definition at line 112 of file L1MuBMEUX.h.

Referenced by id(), operator==(), and run().

unsigned short int L1MuBMEUX::m_quality

private

Definition at line 115 of file L1MuBMEUX.h.

Referenced by operator==(), quality(), reset(), and run().

bool L1MuBMEUX::m_result

private

Definition at line 114 of file L1MuBMEUX.h.

Referenced by operator==(), reset(), result(), and run().

const L1MuBMSEU& L1MuBMEUX::m_seu

private

Definition at line 111 of file L1MuBMEUX.h.

Referenced by load(), and run().

const L1MuBMSectorProcessor& L1MuBMEUX::m_sp

private

Definition at line 110 of file L1MuBMEUX.h.

Referenced by run().

const L1MuBMTrackSegPhi* L1MuBMEUX::m_start

private

Definition at line 118 of file L1MuBMEUX.h.

Referenced by load(), reset(), run(), and ts().

const L1MuBMTrackSegPhi* L1MuBMEUX::m_target

private

Definition at line 119 of file L1MuBMEUX.h.

Referenced by load(), reset(), run(), and ts().

unsigned short const L1MuBMEUX::nbit_phi = 12

private

Definition at line 126 of file L1MuBMEUX.h.

Referenced by run().

unsigned short const L1MuBMEUX::nbit_phib = 10

private

Definition at line 127 of file L1MuBMEUX.h.

Referenced by run().

L1MuDTTFParameters L1MuBMEUX::pars

private

Definition at line 130 of file L1MuBMEUX.h.

Referenced by run().

int const L1MuBMEUX::theExtFilter = 1

private

Definition at line 125 of file L1MuBMEUX.h.

L1MuBMLUTHandler* L1MuBMEUX::theExtLUTs

private

Definition at line 124 of file L1MuBMEUX.h.

Referenced by run().