HcalFinegrainBit Class Reference

#include <HcalFinegrainBit.h>

Public Types
typedef std::array< std::bitset< 6 >, 2 >	Tower
typedef std::array< std::pair< std::pair< int, bool >, std::pair< int, int > >, 7 >	TowerTDC

Public Member Functions
std::bitset< 2 >	compute (const Tower &) const
std::bitset< 6 >	compute (const TowerTDC &, const HcalTrigTowerDetId &) const
	HcalFinegrainBit (int version)

Private Attributes
const int	is_above_mip = 1
const int	is_mip = 0
const int	tdc_boundary [29][7]
int	version_

Detailed Description

Definition at line 8 of file HcalFinegrainBit.h.

Member Typedef Documentation

Tower

typedef std::array<std::bitset<6>, 2> HcalFinegrainBit::Tower

Definition at line 12 of file HcalFinegrainBit.h.

TowerTDC

typedef std::array<std::pair<std::pair<int, bool>, std::pair<int, int> >, 7> HcalFinegrainBit::TowerTDC

Definition at line 14 of file HcalFinegrainBit.h.

Constructor & Destructor Documentation

HcalFinegrainBit()

HcalFinegrainBit::HcalFinegrainBit ( int  version )

inline

Definition at line 16 of file HcalFinegrainBit.h.

: version_(version){};

Member Function Documentation

compute() [1/2]

std::bitset< 2 > HcalFinegrainBit::compute

(

const Tower &

tower

)

const

Definition at line 5 of file HcalFinegrainBit.cc.

References is_above_mip, is_mip, mps_fire::result, l1tHGCalTowerProducer_cfi::tower, and version_.

Referenced by HcalTriggerPrimitiveAlgo::analyzeQIE11().

                                                                               {
  if (version_ == 0) {
    std::bitset<2> result;

    // First layer consistent with a MIP
    result[0] = tower[is_mip][0];

    // First layer consistent with a MIP, at least one layer with more
    // than MIP energy deposition
    result[1] = result[0] && (tower[is_above_mip].count() > 0);

    // There layers consistent with a MIP
    //    result[2] = tower[is_mip].count() >= 3;

    // Unset
    //    result[3] = false;

    return result;
  }
  if (version_ == 1) {
    std::bitset<2> result;

    // All algorithms the same for testing purposes
    result[0] = result[1] = tower[is_mip][0];

    return result;
  }
  if (version_ == 2) {
    std::bitset<2> result;

    // All algorithms the same for testing purposes
    result[0] = result[1] = true;

    return result;
  }
  return 0;
}

compute() [2/2]

std::bitset< 6 > HcalFinegrainBit::compute	(	const TowerTDC &	tower,
		const HcalTrigTowerDetId &	id
	)		const

Definition at line 44 of file HcalFinegrainBit.cc.

References funct::abs(), mps_fire::i, gmt_cfi::Nprompt, mps_fire::result, tdc_boundary, and l1tHGCalTowerProducer_cfi::tower.

                                                                                                                {
  std::bitset<6> result;

  int tp_ieta = id.ieta();

  int Ndelayed = 0;
  int NveryDelayed = 0;
  int Nprompt = 0;
  int EarlyEnergy = 0;
  int DeepEnergy = 0;

  for (size_t i = 0; i < 7; i++) {
    int bit12_15set = tower[i].first.first;
    bool is_compressed = tower[i].first.second;
    int bit12 = (bit12_15set & 0b0001);       // low depth 1,2 energy
    int bit13 = (bit12_15set & 0b0010) >> 1;  // high depth 3+ energy
    int bit14 = (bit12_15set & 0b0100) >> 2;  // prompt energy passed
    int bit15 = (bit12_15set & 0b1000) >> 3;  // delayed energy passed
    int TDC = tower[i].second.first;
    int ADC = tower[i].second.second;

    // timing bits
    if (TDC < 50) {  // exclude error code for TDC in HE (unpacked)
      if ((abs(tp_ieta) <= 16) || (i >= 1)) {
        // count delayed / prompt hits either in HB, or in HE (excluding depth 1 due to backgrounds in HE)
        // Sim packing into Raw, has uncompressed TDC values (0-49) at the trigger primitive level. Packing (compressing HB TDC 6:2) happens in packer.
        // Hcal digis have compressed HB TDC (0-3)
        if (is_compressed == 1) {
          if (TDC == 1 && bit15 == 1)
            Ndelayed += 1;
          if (TDC == 2 && bit15 == 1)
            NveryDelayed += 1;
          if (TDC == 0 && bit14 == 1)
            Nprompt += 1;
        } else {
          if (TDC > tdc_boundary[abs(tp_ieta) - 1][i] && TDC <= tdc_boundary[abs(tp_ieta) - 1][i] + 2 && bit15 == 1)
            Ndelayed += 1;
          if (TDC > tdc_boundary[abs(tp_ieta) - 1][i] + 2 && bit15 == 1)
            NveryDelayed += 1;
          if (TDC <= tdc_boundary[abs(tp_ieta) - 1][i] && TDC >= 0 && bit14 == 1)
            Nprompt += 1;
        }
      }
    }

    // depth bit
    if (ADC > 0 && bit12 == 0)
      EarlyEnergy +=
          1;  // early layers, depth 1 and 2. If bit12 = 0, too much energy in early layers. Require ADC > 0 to ensure valid hit in cell
    if (ADC > 0 && bit13 == 1)
      DeepEnergy +=
          1;  // deep layers, 3+. If bit13 = 1, energy in deep layers. Require ADC > 0 to ensure valid hit in cell
  }

  // very delayed (001000), slightly delayed (000100), prompt (000010), depth flag (000001), 2 reserved bits (110000)
  if (DeepEnergy > 0 && EarlyEnergy == 0)
    result[0] = true;  // 000001
  else
    result[0] = false;
  if (Nprompt > 0)
    result[1] = true;  // 000010
  else
    result[1] = false;
  if (Ndelayed > 0)
    result[2] = true;  // 000100
  else
    result[2] = false;
  if (NveryDelayed > 0)
    result[3] = true;  // 001000
  else
    result[3] = false;
  result[4] = result[5] = false;  // 110000 in HcalTriggerPrimitiveAlgo.cc, set to MIP bits from above

  return result;
}

Member Data Documentation

is_above_mip

const int HcalFinegrainBit::is_above_mip = 1

private

Definition at line 24 of file HcalFinegrainBit.h.

Referenced by compute().

is_mip

const int HcalFinegrainBit::is_mip = 0

private

Definition at line 23 of file HcalFinegrainBit.h.

Referenced by compute().

tdc_boundary

const int HcalFinegrainBit::tdc_boundary[29][7]

private

Initial value:

= {
      {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0},
      {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0},
      {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0},
      {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0}, {12, 12, 12, 12, 0, 0, 0},
      {0, 12, 10, 0, 0, 0, 0},   {0, 9, 10, 9, 10, 0, 0},   {16, 9, 9, 9, 11, 10, 0},  {17, 8, 9, 8, 9, 10, 0},
      {9, 7, 7, 7, 9, 6, 0},     {8, 7, 7, 6, 6, 6, 0},     {8, 6, 6, 6, 7, 7, 0},     {7, 6, 6, 6, 7, 6, 0},
      {7, 6, 6, 6, 6, 6, 0},     {6, 6, 6, 6, 6, 6, 0},     {6, 5, 6, 6, 6, 7, 10},    {9, 9, 9, 5, 5, 6, 6},
      {0, 0, 0, 0, 0, 0, 0}}

Definition at line 29 of file HcalFinegrainBit.h.

Referenced by compute().

version_

int HcalFinegrainBit::version_

private

Definition at line 26 of file HcalFinegrainBit.h.

Referenced by compute().