L1CaloRegionDetId.cc

Go to the documentation of this file.

#include "DataFormats/L1CaloTrigger/interface/L1CaloRegionDetId.h"

unsigned const L1CaloRegionDetId::N_PHI=18;
unsigned const L1CaloRegionDetId::N_ETA=22;


// default constructor - null DetId
L1CaloRegionDetId::L1CaloRegionDetId() : DetId() { }


// construct from raw id
L1CaloRegionDetId::L1CaloRegionDetId(uint32_t rawid) : DetId(rawid) { }


// construct from ieta, iphi indices
// ieta runs from 0 (at -z) to 21 (at +z)
L1CaloRegionDetId::L1CaloRegionDetId(unsigned ieta, unsigned iphi) :
  DetId(Calo, 2) 
{ 
  id_ |= (ieta & 0x1f) | ((iphi & 0x1f)<<5);
}


// construct from RCT crate, card, region IDs
L1CaloRegionDetId::L1CaloRegionDetId(unsigned icrate, unsigned icard, unsigned irgn) :
  DetId(Calo, 2)
{

  int ieta=0;
  int iphi=0;

  // Calculate iphi
  int phi_index = icrate % 9;
  if ((icard == 0) || (icard == 2) || (icard == 4))
    phi_index = phi_index * 2;
  else if ((icard == 1) || (icard == 3) || (icard == 5))
    phi_index = phi_index * 2 + 1;
  else if (icard == 6)
    phi_index = phi_index * 2 + irgn;
  // for HF
  else if (icard == 999)
    phi_index = phi_index * 2 + (irgn/4);
  iphi = (22 - phi_index) % 18;

  // Calculate ieta
  int eta_index = 0;
  if (icard < 6)
    eta_index = (icard/2) * 2 + irgn;
  else if (icard == 6)
    eta_index = 6;
  // HF
  else if (icard == 999)
    eta_index = (irgn % 4) + 7;
  
  if (icrate < 9)
    ieta = 10 - eta_index;
  else if (icrate >= 9)
    ieta = 11 + eta_index;

  id_ |= (ieta & 0x1f) | ((iphi & 0x1f)<<5);
}

// return RCT crate ID
unsigned L1CaloRegionDetId::rctCrate() const { // TODO - check this is correct!
  unsigned phiCrate = ((N_PHI + 4 - iphi()) % N_PHI) / 2;
  return (ieta()<(N_ETA/2) ? phiCrate : phiCrate + N_PHI/2) ;
}

// return RCT card number
unsigned L1CaloRegionDetId::rctCard() const {
  unsigned card = 999;
  unsigned rct_phi_index = (22 - iphi()) % 18;
  if ((ieta() == 4) || (ieta() == 17)){
    card = 6;
  }
  else if ((ieta() > 4) && (ieta() <= 10)){
    unsigned index = (ieta() - 5)/2;
    card = ((2 - index) * 2) + (rct_phi_index % 2);
  }
  else if ((ieta() >= 11) && (ieta() < 17)){
    unsigned index = (ieta() - 11)/2;
    card = (index * 2) + (rct_phi_index % 2);
  }
  return card;
}

// return RCT region number
unsigned L1CaloRegionDetId::rctRegion() const {
  unsigned rgn = 999;
  unsigned rct_phi_index = (22 - iphi()) % 18;
  if (ieta() < 4){
    rgn = (3 - ieta()) + 4 * (rct_phi_index % 2);
  }
  else if (ieta() > 17){
    rgn = (ieta() - 18) + 4 * (rct_phi_index % 2);
  }
  else if ((ieta() == 4) || (ieta() == 17)){
    rgn = (rct_phi_index % 2);
  }
  else if ((ieta() > 4) && (ieta() <= 10)){
    rgn = (ieta() % 2);
  }
  else if ((ieta() >= 11) && (ieta() < 17)){
    rgn = ((ieta() - 1) % 2);
  }
  return rgn;
}