CMS 3D CMS Logo
Public Member Functions | Private Member Functions | Private Attributes

L1MuGMTMerger Class Reference

#include <L1MuGMTMerger.h>

List of all members.

Public Member Functions

const std::vector
< L1MuGMTExtendedCand * > & 		Cands () const
 return std::vector with all muon candidates
int id () const
 return identifier
 L1MuGMTMerger (const L1MuGlobalMuonTrigger &gmt, int id)
 constructor
void print () const
 print results after selection
void reset ()
 clear Merger
void run ()
 run GMT Merger
virtual ~L1MuGMTMerger ()
 destructor

Private Member Functions

unsigned convertedEta (const L1MuRegionalCand *mu) const
void createDTCSCCand (int idx_dtcsc)
void createMergedCand (int idx_dtcsc, int idx_rpc)
void createRPCCand (int idx_rpc)
int doANDMerge (unsigned MMconfig) const
int doSpecialMerge (unsigned MMconfig) const
void load ()
void merge ()
int merge_rank (const L1MuRegionalCand *) const
 Merge Rank Table.
unsigned projectedPhi (const L1MuRegionalCand *mu) const
int selectDTCSC (unsigned MMconfig, int by_rank, int by_pt, int by_combi) const
unsigned sysign (const L1MuRegionalCand *mu) const

Private Attributes

std::vector< const
L1MuRegionalCand * > 		dtcsc_mu
const L1MuGlobalMuonTrigger & m_gmt
int m_id
std::vector
< L1MuGMTExtendedCand * > 		m_MuonCands
std::vector< const
L1MuRegionalCand * > 		rpc_mu
std::vector< int > singleRank

Detailed Description

L1 Global Muon Trigger Merger.

There are two GMT Mergers. One for the barrel (id=0) and one for the forward (id=1) part.

The merger receives four DT(CSC) muons and four RPC muons. Each DT(CSC) muon is either merged with an RPC one or passed through.

To simplify the C++ model, the merger conatians also the conversion units and sort rank units as well as the MergeMethodSelection unit which are separate units in the hardware.

Definition at line 63 of file L1MuGMTMerger.h.

Constructor & Destructor Documentation

L1MuGMTMerger::L1MuGMTMerger ( const L1MuGlobalMuonTrigger &  gmt,
int  id 
)

constructor

Definition at line 73 of file L1MuGMTMerger.cc.

References dtcsc_mu, and rpc_mu.

                                                                     : 
  m_gmt(gmt), m_id(id), dtcsc_mu(4), rpc_mu(4) {

  dtcsc_mu.reserve(4);
  rpc_mu.reserve(4);
   
}
L1MuGMTMerger::~L1MuGMTMerger ( ) [virtual]

destructor

Definition at line 84 of file L1MuGMTMerger.cc.

References reset().

                              { 
  reset();  
}

Member Function Documentation

const std::vector<L1MuGMTExtendedCand*>& L1MuGMTMerger::Cands ( ) const [inline]

return std::vector with all muon candidates

Definition at line 86 of file L1MuGMTMerger.h.

References m_MuonCands.

Referenced by L1MuGMTSorter::run().

{ return m_MuonCands; }
unsigned L1MuGMTMerger::convertedEta ( const L1MuRegionalCand *  mu) const [private]

Definition at line 321 of file L1MuGMTMerger.cc.

References L1MuRegionalCand::eta_packed(), L1MuGMTConfig::getLFEtaConvLUT(), L1MuGMTLFEtaConvLUT::SpecificLookup_eta_gmt(), and L1MuRegionalCand::type_idx().

Referenced by createDTCSCCand(), createMergedCand(), and createRPCCand().

                                                                     {
  L1MuGMTLFEtaConvLUT* etaconv_lut = L1MuGMTConfig::getLFEtaConvLUT();
  return etaconv_lut->SpecificLookup_eta_gmt (mu->type_idx(), mu->eta_packed() );
}
void L1MuGMTMerger::createDTCSCCand ( int  idx_dtcsc) [private]

Definition at line 242 of file L1MuGMTMerger.cc.

References convertedEta(), dtcsc_mu, L1MuGMTSortRankUnit::getVeryLowQualityLevel(), L1MuGMTMipIsoAU::ISO(), m_gmt, m_id, m_MuonCands, L1MuGMTMipIsoAU::MIP(), L1MuGlobalMuonTrigger::MipIsoAU(), projectedPhi(), L1MuGMTCand::setBx(), L1MuGMTCand::setChargePacked(), L1MuGMTExtendedCand::setDTCSCIndex(), L1MuGMTCand::setEtaPacked(), L1MuGMTExtendedCand::setFwdBit(), L1MuGMTCand::setIsolation(), L1MuGMTCand::setMIP(), L1MuGMTCand::setPhiPacked(), L1MuGMTCand::setPtPacked(), L1MuGMTCand::setQuality(), L1MuGMTExtendedCand::setRank(), L1MuGMTExtendedCand::setRPCBit(), L1MuGMTExtendedCand::setRPCIndex(), L1MuGMTSortRankUnit::sort_rank(), and sysign().

Referenced by merge().

                                                 {
  L1MuGMTExtendedCand* tmpmuon = new L1MuGMTExtendedCand();
   
  tmpmuon->setBx(dtcsc_mu[idx_dtcsc]->bx());
  tmpmuon->setPhiPacked( projectedPhi(dtcsc_mu[idx_dtcsc]) ); 
  tmpmuon->setEtaPacked( convertedEta(dtcsc_mu[idx_dtcsc]) ); 
  tmpmuon->setPtPacked( dtcsc_mu[idx_dtcsc]->pt_packed() ); 
  tmpmuon->setChargePacked( sysign(dtcsc_mu[idx_dtcsc]) ); 
  tmpmuon->setMIP( m_gmt.MipIsoAU(m_id)->MIP(idx_dtcsc) );    
  tmpmuon->setIsolation ( m_gmt.MipIsoAU(m_id)->ISO(idx_dtcsc) );    
  tmpmuon->setRank( L1MuGMTSortRankUnit::sort_rank(dtcsc_mu[idx_dtcsc]) );

  unsigned quality = 0;
  switch ( L1MuGMTSortRankUnit::getVeryLowQualityLevel(dtcsc_mu[idx_dtcsc]) ) {
  case 0: quality = 6; break; //DT/CSC
  case 1: quality = 2; break; //VERY LOW QUALITY LEVEL 1
  case 2: quality = 3; break; //VERY LOW QUALITY LEVEL 2
  case 3: quality = 4; break; //VERY LOW QUALITY LEVEL 3
  }

  if ( (m_id==1) && (dtcsc_mu[idx_dtcsc]->finehalo_packed() == 1) ) 
    quality = 1; // HALO quality

  tmpmuon->setQuality( quality ); // RPC
  tmpmuon->setDTCSCIndex(idx_dtcsc);
  tmpmuon->setRPCIndex(0);
  tmpmuon->setFwdBit( m_id );
  tmpmuon->setRPCBit( 0 );

  m_MuonCands.push_back(tmpmuon);
}
void L1MuGMTMerger::createMergedCand ( int  idx_dtcsc,
int  idx_rpc 
) [private]

Definition at line 359 of file L1MuGMTMerger.cc.

References abs, convertedEta(), doANDMerge(), doSpecialMerge(), dtcsc_mu, eta(), L1MuGMTConfig::getLFPtMixLUT(), L1MuGMTConfig::getRegMMConfigCharge(), L1MuGMTConfig::getRegMMConfigEta(), L1MuGMTConfig::getRegMMConfigISO(), L1MuGMTConfig::getRegMMConfigMIP(), L1MuGMTConfig::getRegMMConfigPhi(), L1MuGMTConfig::getRegMMConfigPt(), L1MuGMTConfig::getRegMMConfigSRK(), L1MuGMTConfig::getRegSortRankOffset(), L1MuGMTReg::getValue(), L1MuGMTMipIsoAU::ISO(), m_gmt, m_id, m_MuonCands, merge_rank(), L1MuGMTMipIsoAU::MIP(), L1MuGlobalMuonTrigger::MipIsoAU(), phi, projectedPhi(), rpc_mu, selectDTCSC(), L1MuGMTCand::setBx(), L1MuGMTCand::setChargePacked(), L1MuGMTExtendedCand::setDTCSCIndex(), L1MuGMTCand::setEtaPacked(), L1MuGMTExtendedCand::setFwdBit(), L1MuGMTCand::setIsolation(), L1MuGMTCand::setMIP(), L1MuGMTCand::setPhiPacked(), L1MuGMTCand::setPtPacked(), L1MuGMTCand::setQuality(), L1MuGMTExtendedCand::setRank(), L1MuGMTExtendedCand::setRPCBit(), L1MuGMTExtendedCand::setRPCIndex(), L1MuGMTSortRankUnit::sort_rank(), L1MuGMTLFPtMixLUT::SpecificLookup_pt_mixed(), and sysign().

Referenced by merge().

                                                               {

  // In the hardware matrices of select_bits are calculated for all
  // possible pairings. 
  // In ORCA we only calculate selec-bits for the actual 
  // pairs to save time.

  // look up merge ranks
  int merge_rank_dtcsc = merge_rank(dtcsc_mu[idx_dtcsc]);
  int merge_rank_rpc   = merge_rank(rpc_mu[idx_rpc]);

  // calculate select-bits (1: take DT/CSC, 0: take RPC)
  int selected_by_rank = abs(merge_rank_dtcsc) > abs(merge_rank_rpc);
  int selected_by_pt   = dtcsc_mu[idx_dtcsc]->pt_packed() <= rpc_mu[idx_rpc]->pt_packed();

  // Selection by combination of min pt and higher rank
  // select by rank if both flags are set, otherwise by min pt
  // in other words: select by minpt if one flag is not set
  int selected_by_combi = (merge_rank_dtcsc<0 && merge_rank_rpc<0) ? selected_by_rank : selected_by_pt;

  
  L1MuGMTExtendedCand* tmpmuon = new L1MuGMTExtendedCand();
   
  tmpmuon->setBx(dtcsc_mu[idx_dtcsc]->bx());


  // merge phi
  //  unsigned MMConfig_phi = 32; // take DT
  unsigned MMConfig_phi = L1MuGMTConfig::getRegMMConfigPhi()->getValue(m_id);

  unsigned phi = 0;
  
  if (selectDTCSC(MMConfig_phi, selected_by_rank, selected_by_pt, selected_by_combi) )
    phi = projectedPhi( dtcsc_mu[idx_dtcsc] );
  else
    phi = projectedPhi( rpc_mu[idx_rpc] );

  tmpmuon->setPhiPacked( phi );
   
  // merge eta
  unsigned MMConfig_eta = L1MuGMTConfig::getRegMMConfigEta()->getValue(m_id);

  unsigned eta = 0;
  
  if (doSpecialMerge( MMConfig_eta ) ) {
    if ( (m_id == 1) || dtcsc_mu[idx_dtcsc]->finehalo_packed() ) 
      eta = convertedEta( dtcsc_mu[idx_dtcsc] );
    else
      eta = convertedEta( rpc_mu[idx_rpc] );
  } 
  else {
    if ( selectDTCSC(MMConfig_eta, selected_by_rank, selected_by_pt, selected_by_combi) )
      eta = convertedEta( dtcsc_mu[idx_dtcsc] );
    else
      eta = convertedEta( rpc_mu[idx_rpc] );    
  }
  tmpmuon->setEtaPacked( eta );

  // merge pt
  unsigned MMConfig_pt = L1MuGMTConfig::getRegMMConfigPt()->getValue(m_id);

  unsigned pt = 0;

  if (doSpecialMerge( MMConfig_pt ) ) { // mix pt
    L1MuGMTLFPtMixLUT* ptmix_lut = L1MuGMTConfig::getLFPtMixLUT();
    pt = ptmix_lut->SpecificLookup_pt_mixed (m_id, dtcsc_mu[idx_dtcsc]->pt_packed(), 
                                             rpc_mu[idx_rpc]->pt_packed());    
  } 
  else {
    if ( selectDTCSC(MMConfig_pt, selected_by_rank, selected_by_pt, selected_by_combi) )
      pt = dtcsc_mu[idx_dtcsc]->pt_packed();
    else
      pt = rpc_mu[idx_rpc]->pt_packed();
  }
  tmpmuon->setPtPacked( pt );

  // merge charge
  unsigned MMConfig_charge = L1MuGMTConfig::getRegMMConfigCharge()->getValue(m_id);

  unsigned sy_sign = 0;

  if (doSpecialMerge( MMConfig_charge ) ) {
    // based on charge valid bits
    if ( rpc_mu[idx_rpc]->charge_valid_packed()==1 && dtcsc_mu[idx_dtcsc]->charge_valid_packed()==0 )
      sy_sign = sysign(rpc_mu[idx_rpc]);
    else 
      sy_sign = sysign(dtcsc_mu[idx_dtcsc]);
  } 
  else {
    if ( selectDTCSC(MMConfig_charge, selected_by_rank, selected_by_pt, selected_by_combi) )
      sy_sign = sysign(dtcsc_mu[idx_dtcsc]);
    else
      sy_sign = sysign(rpc_mu[idx_rpc]);
  }
  tmpmuon->setChargePacked( sy_sign );
   
  // merge quality  

  
  // merge MIP
  unsigned MMConfig_MIP = L1MuGMTConfig::getRegMMConfigMIP()->getValue(m_id);

  bool mip_bit = false;

  bool mip_bit_dtcsc = m_gmt.MipIsoAU(m_id)->MIP(idx_dtcsc);
  bool mip_bit_rpc   = m_gmt.MipIsoAU(m_id)->MIP(idx_rpc+4);
 
  if ( doSpecialMerge(MMConfig_MIP) ) {
    if ( doANDMerge(MMConfig_MIP) )
      mip_bit =  mip_bit_dtcsc && mip_bit_rpc;
    else
      mip_bit =  mip_bit_dtcsc || mip_bit_rpc;
  }
  else {
    if ( selectDTCSC(MMConfig_MIP, selected_by_rank, selected_by_pt, selected_by_combi) )
      mip_bit = mip_bit_dtcsc;
    else
      mip_bit = mip_bit_rpc;
  }

  tmpmuon->setMIP(mip_bit);
 
  // merge ISO
  unsigned MMConfig_ISO = L1MuGMTConfig::getRegMMConfigISO()->getValue(m_id);

  bool iso_bit = false;

  bool iso_bit_dtcsc = m_gmt.MipIsoAU(m_id)->ISO(idx_dtcsc);
  bool iso_bit_rpc   = m_gmt.MipIsoAU(m_id)->ISO(idx_rpc+4);
 
  if ( doSpecialMerge(MMConfig_ISO) ) {
    if ( doANDMerge(MMConfig_ISO) )
      iso_bit =  iso_bit_dtcsc && iso_bit_rpc;
    else
      iso_bit =  iso_bit_dtcsc || iso_bit_rpc;
  }
  else {
    if ( selectDTCSC(MMConfig_ISO, selected_by_rank, selected_by_pt, selected_by_combi) )
      iso_bit = iso_bit_dtcsc;
    else
      iso_bit = iso_bit_rpc;
  }

  tmpmuon->setIsolation(iso_bit);

  // merge sort rank
  unsigned MMConfig_SRK = L1MuGMTConfig::getRegMMConfigSRK()->getValue(m_id);

  unsigned rank_offset = L1MuGMTConfig::getRegSortRankOffset()->getValue(m_id);

  unsigned rank = 0;
  if ( selectDTCSC(MMConfig_SRK, selected_by_rank, selected_by_pt, selected_by_combi) )
    rank = L1MuGMTSortRankUnit::sort_rank(dtcsc_mu[idx_dtcsc]) + rank_offset; 
  else
    rank = L1MuGMTSortRankUnit::sort_rank(rpc_mu[idx_rpc])     + rank_offset;

  tmpmuon->setRank(rank); 

  // quality of merged candidate
  tmpmuon->setQuality( 7 ); // code for matched muons

  tmpmuon->setDTCSCIndex(idx_dtcsc);
  tmpmuon->setRPCIndex(idx_rpc);
  tmpmuon->setFwdBit( m_id );
  tmpmuon->setRPCBit( 0 );
   
  m_MuonCands.push_back(tmpmuon);
}
void L1MuGMTMerger::createRPCCand ( int  idx_rpc) [private]

Definition at line 275 of file L1MuGMTMerger.cc.

References convertedEta(), L1MuGMTSortRankUnit::getVeryLowQualityLevel(), L1MuGMTMipIsoAU::ISO(), m_gmt, m_id, m_MuonCands, L1MuGMTMipIsoAU::MIP(), L1MuGlobalMuonTrigger::MipIsoAU(), projectedPhi(), rpc_mu, L1MuGMTCand::setBx(), L1MuGMTCand::setChargePacked(), L1MuGMTExtendedCand::setDTCSCIndex(), L1MuGMTCand::setEtaPacked(), L1MuGMTExtendedCand::setFwdBit(), L1MuGMTCand::setIsolation(), L1MuGMTCand::setMIP(), L1MuGMTCand::setPhiPacked(), L1MuGMTCand::setPtPacked(), L1MuGMTCand::setQuality(), L1MuGMTExtendedCand::setRank(), L1MuGMTExtendedCand::setRPCBit(), L1MuGMTExtendedCand::setRPCIndex(), L1MuGMTSortRankUnit::sort_rank(), and sysign().

Referenced by merge().

                                             {
  L1MuGMTExtendedCand* tmpmuon = new L1MuGMTExtendedCand();
   
  tmpmuon->setBx(rpc_mu[idx_rpc]->bx()); 
  tmpmuon->setPhiPacked( projectedPhi(rpc_mu[idx_rpc]) ); 
  tmpmuon->setEtaPacked( convertedEta(rpc_mu[idx_rpc]) ); 
  tmpmuon->setPtPacked( rpc_mu[idx_rpc]->pt_packed() ); 
  tmpmuon->setChargePacked( sysign( rpc_mu[idx_rpc] ) ); 
  tmpmuon->setMIP( m_gmt.MipIsoAU(m_id)->MIP(idx_rpc+4) );    
  tmpmuon->setIsolation ( m_gmt.MipIsoAU(m_id)->ISO(idx_rpc+4) );    
  tmpmuon->setRank( L1MuGMTSortRankUnit::sort_rank(rpc_mu[idx_rpc]) );

  unsigned quality = 0;
  switch ( L1MuGMTSortRankUnit::getVeryLowQualityLevel(rpc_mu[idx_rpc]) ) {
  case 0: quality = 5; break; //RPC
  case 1: quality = 2; break; //VERY LOW QUALITY LEVEL1
  case 2: quality = 3; break; //VERY LOW QUALITY LEVEL2
  case 3: quality = 4; break; //VERY LOW QUALITY LEVEL3
  }

  tmpmuon->setQuality( quality ); // RPC
  tmpmuon->setDTCSCIndex(0);
  tmpmuon->setRPCIndex(idx_rpc);
  tmpmuon->setFwdBit( m_id );
  tmpmuon->setRPCBit( 1 );
    
  m_MuonCands.push_back(tmpmuon);
}
int L1MuGMTMerger::doANDMerge ( unsigned  MMconfig) const [private]

Definition at line 317 of file L1MuGMTMerger.cc.

Referenced by createMergedCand().

                                                     {
  return (MMconfig & 64) == 64;
}
int L1MuGMTMerger::doSpecialMerge ( unsigned  MMconfig) const [private]

Definition at line 313 of file L1MuGMTMerger.cc.

Referenced by createMergedCand().

                                                         {
  return (MMconfig & 1) == 1;
}
int L1MuGMTMerger::id ( void  ) const [inline]

return identifier

Definition at line 83 of file L1MuGMTMerger.h.

References m_id.

{ return m_id; }
void L1MuGMTMerger::load ( ) [private]

Definition at line 142 of file L1MuGMTMerger.cc.

References L1MuGMTPSB::CSCMuon(), L1MuGlobalMuonTrigger::Data(), L1MuGMTPSB::DTBXMuon(), dtcsc_mu, m_gmt, m_id, MAXCSC, MAXDTBX, L1MuGMTConfig::MAXRPCbarrel, L1MuGMTConfig::MAXRPCendcap, rpc_mu, and L1MuGMTPSB::RPCMuon().

Referenced by run().

                         {

  // barrel Merger gets DTBX and barrel RPC muons
  if ( m_id == 0 ) {
    for ( unsigned idt = 0; idt < L1MuGMTConfig::MAXDTBX; idt++ ) {
      dtcsc_mu[idt] = m_gmt.Data()->DTBXMuon(idt);
    }
    for ( unsigned irpc = 0; irpc < L1MuGMTConfig::MAXRPCbarrel; irpc++ ) {
      rpc_mu[irpc] = m_gmt.Data()->RPCMuon(irpc);
    }
  }
  
  // endcap Merger gets CSC and endcap RPC muons
  if ( m_id == 1 ) {  
    for ( unsigned icsc = 0; icsc < L1MuGMTConfig::MAXCSC; icsc++ ) {
      dtcsc_mu[icsc] = m_gmt.Data()->CSCMuon(icsc);
    }
    for ( unsigned irpc = 0; irpc < L1MuGMTConfig::MAXRPCendcap; irpc++ ) {
      rpc_mu[irpc] = m_gmt.Data()->RPCMuon(irpc+4);
    }
  }
 
}
void L1MuGMTMerger::merge ( ) [private]

Definition at line 167 of file L1MuGMTMerger.cc.

References L1MuGMTCancelOutUnit::cancelMyChipMuon(), L1MuGMTCancelOutUnit::cancelOtherChipMuon(), L1MuGlobalMuonTrigger::CancelOutUnit(), L1MuGMTMatrix< T >::colAny(), createDTCSCCand(), createMergedCand(), createRPCCand(), dtcsc_mu, L1MuGMTCand::etaIndex(), L1MuScale::getCenter(), L1MuTriggerScales::getGMTEtaScale(), L1MuScale::getLowEdge(), L1MuTriggerScales::getPhiScale(), L1MuTriggerPtScale::getPtScale(), L1MuGMTConfig::getTriggerPtScale(), L1MuGMTConfig::getTriggerScales(), i, L1MuGMTSortRankUnit::isDisabled(), j, m_gmt, m_id, m_MuonCands, L1MuGlobalMuonTrigger::Matcher(), L1MuGMTMatcher::pairM(), L1MuGMTCand::phiIndex(), L1MuGMTCand::ptIndex(), L1MuGMTMatrix< T >::rowAny(), rpc_mu, L1MuGMTCand::setEtaValue(), L1MuGMTCand::setPhiValue(), and L1MuGMTCand::setPtValue().

Referenced by run().

                          {

  const L1MuGMTMatrix<bool>& pairM = m_gmt.Matcher(m_id)->pairM();

  // Handling of cancel-out and empty muons is different in software and hardware
  //
  // - in hardware, a candidate may be empty (pt_code == 0) and an empty bit
  //   is passed to the sorter in order to suppress the candidate. In the
  //   software no candidate is created in this case.
  //
  // - in hardware, RPC candidates are passed to the sorter even if they are 
  //   also used in a matched pair. They are then suppressed in the sorter. In
  //   software no RPC candidate is created if the muon is used in a pair.
  //
  // - in hardware, cancel-out signals from the cancel-out units in the own and   
  //   other Logic FPGA are passed to the sorter in order to cancel out muons.
  //   In software the cancel-out signals are alrady checked here in the merger 
  //   and no candidates are created for cancelled muons.
  //
  // There may therefore be less muons passed to the sorter in software than 
  // in Hardware. At the output of the first sorter stage the results should 
  // be comparable, again.
  
  unsigned HaloOverwritesMatched = 1;
  
  // loop over DT/CSC muons
  for (int i=0; i<4; i++) {
    if (dtcsc_mu[i] != 0) {

      int match_idx = pairM.rowAny(i);

      int csc_is_halo = (m_id==1) && (dtcsc_mu[i]->finehalo_packed() == 1);

      if ( (match_idx != -1) && // is it matched?
           (! (csc_is_halo && HaloOverwritesMatched) ) )
        createMergedCand(i, match_idx);         
      else {
        // check my first and the other chip's second cancel-out units
        if ((! m_gmt.CancelOutUnit(m_id)->cancelMyChipMuon(i) ) &&
            (! m_gmt.CancelOutUnit(3-m_id)->cancelOtherChipMuon(i) ) &&
            (! L1MuGMTSortRankUnit::isDisabled(dtcsc_mu[i]) ) )
          createDTCSCCand(i);
      }
    }
  }

  // additionally loop over RPC muons
  for (int j=0; j<4; j++) {
    if (rpc_mu[j] != 0) {      
      
      int match_idx = pairM.colAny(j);
      
      if (match_idx == -1) { // is it unmatched?
        if ((! m_gmt.CancelOutUnit(m_id+2)->cancelMyChipMuon(j) ) &&
            (! L1MuGMTSortRankUnit::isDisabled(rpc_mu[j]) ) )
          createRPCCand(j);
      }
    }
  }

  // set physical values in the GMT candidates for use in the analysis
  const L1MuTriggerScales* theTriggerScales = L1MuGMTConfig::getTriggerScales();
  const L1MuTriggerPtScale* theTriggerPtScale = L1MuGMTConfig::getTriggerPtScale();
  
  std::vector<L1MuGMTExtendedCand*>::const_iterator icand;
  for(icand=m_MuonCands.begin();icand!=m_MuonCands.end();icand++) {
    L1MuGMTExtendedCand* cand = (*icand);
    cand->setPhiValue( theTriggerScales->getPhiScale()->getLowEdge( cand->phiIndex() ));
    cand->setEtaValue( theTriggerScales->getGMTEtaScale()->getCenter( cand->etaIndex() ));
    cand->setPtValue( theTriggerPtScale->getPtScale()->getLowEdge( cand->ptIndex() ));
    // cand->setPtValue( theTriggerScales->getPtScale()->getLowEdge( cand->ptIndex() ));
  }

}
int L1MuGMTMerger::merge_rank ( const L1MuRegionalCand *  muon) const [private]

Merge Rank Table.

Definition at line 530 of file L1MuGMTMerger.cc.

References L1MuRegionalCand::empty(), eta(), L1MuRegionalCand::eta_packed(), L1MuGMTConfig::getLFMergeRankCombineLUT(), L1MuGMTConfig::getLFMergeRankEtaPhiLUT(), L1MuGMTConfig::getLFMergeRankEtaQLUT(), L1MuGMTConfig::getLFMergeRankPtQLUT(), phi, L1MuRegionalCand::phi_packed(), L1MuRegionalCand::pt_packed(), lumiQueryAPI::q, L1MuRegionalCand::quality_packed(), L1MuGMTLFMergeRankEtaQLUT::SpecificLookup_flag(), L1MuGMTLFMergeRankCombineLUT::SpecificLookup_merge_rank(), L1MuGMTLFMergeRankEtaPhiLUT::SpecificLookup_rank_etaphi(), L1MuGMTLFMergeRankEtaQLUT::SpecificLookup_rank_etaq(), L1MuGMTLFMergeRankPtQLUT::SpecificLookup_rank_ptq(), and L1MuRegionalCand::type_idx().

Referenced by createMergedCand().

                                                                {

  if ( muon == 0 || muon->empty() ) return 0;

  unsigned lut_idx= muon->type_idx();
  
  // obtain inputs as coded in HW
  unsigned eta = muon->eta_packed();
  unsigned q = muon->quality_packed();
  unsigned pt  = muon->pt_packed();
  unsigned phi = muon->phi_packed();

  // lookup eta-q
  L1MuGMTLFMergeRankEtaQLUT* etaq_lut = L1MuGMTConfig::getLFMergeRankEtaQLUT();
  unsigned rank_etaq = etaq_lut->SpecificLookup_rank_etaq (lut_idx, eta, q);
  unsigned flag = etaq_lut->SpecificLookup_flag (lut_idx, eta, q);
   
  // lookup pt-q
  L1MuGMTLFMergeRankPtQLUT* ptq_lut = L1MuGMTConfig::getLFMergeRankPtQLUT();
  unsigned rank_ptq = ptq_lut->SpecificLookup_rank_ptq (lut_idx, q, pt);
  
  // lookup etaphi
  L1MuGMTLFMergeRankEtaPhiLUT* etaphi_lut = L1MuGMTConfig::getLFMergeRankEtaPhiLUT();
  unsigned rank_etaphi = etaphi_lut->SpecificLookup_rank_etaphi (lut_idx, eta, phi);
  
  // combine
  L1MuGMTLFMergeRankCombineLUT* combine_lut = L1MuGMTConfig::getLFMergeRankCombineLUT();
  unsigned rank = combine_lut->SpecificLookup_merge_rank (lut_idx, rank_etaq, rank_ptq, rank_etaphi);

  int rank_signed = rank;

  if (flag == 1) rank_signed *= -1;
   
  return rank_signed;
}
void L1MuGMTMerger::print ( void  ) const

print results after selection

Definition at line 125 of file L1MuGMTMerger.cc.

References m_MuonCands.

Referenced by L1MuGlobalMuonTrigger::produce().

                                {

  edm::LogVerbatim("GMT_Merger_info") << " ";

  std::vector<L1MuGMTExtendedCand*>::const_iterator iter;
  for ( iter = m_MuonCands.begin(); iter != m_MuonCands.end(); iter++ ) {
    if ( *iter && !(*iter)->empty() ) (*iter)->print();
  }
  
  edm::LogVerbatim("GMT_Merger_info") << " ";

}
unsigned L1MuGMTMerger::projectedPhi ( const L1MuRegionalCand *  mu) const [private]

Definition at line 326 of file L1MuGMTMerger.cc.

References L1MuRegionalCand::charge_packed(), L1MuRegionalCand::eta_packed(), L1MuGMTConfig::getLFPhiProEtaConvLUT(), L1MuGMTConfig::getLFPhiProLUT(), L1MuSignedPacking< Bits >::idxFromPacked(), L1MuRegionalCand::phi_packed(), L1MuRegionalCand::pt_packed(), L1MuGMTLFPhiProLUT::SpecificLookup_dphi(), L1MuGMTLFPhiProEtaConvLUT::SpecificLookup_eta_out(), and L1MuRegionalCand::type_idx().

Referenced by createDTCSCCand(), createMergedCand(), and createRPCCand().

                                                                     {

  // convert eta
  L1MuGMTLFPhiProEtaConvLUT* phiproetaconv_lut = L1MuGMTConfig::getLFPhiProEtaConvLUT();
  unsigned eta4 = phiproetaconv_lut->SpecificLookup_eta_out (mu->type_idx(), mu->eta_packed() );
  
  // look up delta-phi 9 bit signed
  L1MuGMTLFPhiProLUT* phipro_lut = L1MuGMTConfig::getLFPhiProLUT();
  unsigned dphi9 = phipro_lut->SpecificLookup_dphi (mu->type_idx(), eta4, mu->pt_packed(), mu->charge_packed());
  
  // sign extend
  L1MuSignedPacking<9> DPhiPacking;
  int dphi = DPhiPacking.idxFromPacked( dphi9 );
  
  // add modulo 144
  int newphi = mu->phi_packed() + dphi; 
  if (newphi < 0) newphi += 144;
  if (newphi >= 144) newphi -= 144;

  return (unsigned) newphi;  
}
void L1MuGMTMerger::reset ( void  )

clear Merger

Definition at line 104 of file L1MuGMTMerger.cc.

References dtcsc_mu, i, m_MuonCands, and rpc_mu.

Referenced by L1MuGlobalMuonTrigger::reset(), and ~L1MuGMTMerger().

                          {

  for ( int i = 0; i < 4; i++ ) {  
    dtcsc_mu[i] = 0;
    rpc_mu[i] = 0;
  }
  
  std::vector<L1MuGMTExtendedCand*>::iterator iter;
  for ( iter = m_MuonCands.begin(); iter != m_MuonCands.end(); iter++ ) {
    if ( *iter ) delete (*iter); 
    *iter = 0;
  }

  m_MuonCands.clear();

}
void L1MuGMTMerger::run ( void  )

run GMT Merger

Definition at line 95 of file L1MuGMTMerger.cc.

References load(), and merge().

Referenced by L1MuGlobalMuonTrigger::produce().

                        {
  load();
  merge();  
}
int L1MuGMTMerger::selectDTCSC ( unsigned  MMconfig,
int  by_rank,
int  by_pt,
int  by_combi 
) const [private]

Definition at line 304 of file L1MuGMTMerger.cc.

Referenced by createMergedCand().

                                                                                            {
  return 
    (   (MMconfig & 32) == 32 ) ||
    ( ( (MMconfig &  8) ==  8 ) && by_rank ) ||
    ( ( (MMconfig &  4) ==  4 ) && by_pt   ) ||
    ( ( (MMconfig &  2) ==  2 ) && by_combi);

}
unsigned L1MuGMTMerger::sysign ( const L1MuRegionalCand *  mu) const [private]

Definition at line 349 of file L1MuGMTMerger.cc.

References L1MuRegionalCand::charge_packed(), and L1MuRegionalCand::charge_valid_packed().

Referenced by createDTCSCCand(), createMergedCand(), and createRPCCand().

                                                               {

  unsigned sysign = mu->charge_packed();

  if ( mu->charge_valid_packed() == 0 )
    sysign = 2; // undefined charge

  return sysign;
}

Member Data Documentation

std::vector<const L1MuRegionalCand*> L1MuGMTMerger::dtcsc_mu [private]

Definition at line 113 of file L1MuGMTMerger.h.

Referenced by createDTCSCCand(), createMergedCand(), L1MuGMTMerger(), load(), merge(), and reset().

const L1MuGlobalMuonTrigger& L1MuGMTMerger::m_gmt [private]

Definition at line 110 of file L1MuGMTMerger.h.

Referenced by createDTCSCCand(), createMergedCand(), createRPCCand(), load(), and merge().

int L1MuGMTMerger::m_id [private]

Definition at line 111 of file L1MuGMTMerger.h.

Referenced by createDTCSCCand(), createMergedCand(), createRPCCand(), id(), load(), and merge().

std::vector<L1MuGMTExtendedCand*> L1MuGMTMerger::m_MuonCands [private]

Definition at line 116 of file L1MuGMTMerger.h.

Referenced by Cands(), createDTCSCCand(), createMergedCand(), createRPCCand(), merge(), print(), and reset().

std::vector<const L1MuRegionalCand*> L1MuGMTMerger::rpc_mu [private]

Definition at line 114 of file L1MuGMTMerger.h.

Referenced by createMergedCand(), createRPCCand(), L1MuGMTMerger(), load(), merge(), and reset().

std::vector<int> L1MuGMTMerger::singleRank [private]

Definition at line 118 of file L1MuGMTMerger.h.

Generated for CMSSW Reference Manual by doxygen 1.7.3