DSMAlgo_EM201_2013.cc

#include "DSM.hh"
#include "DSMAlgo_EM201_2013.hh"

int DSMAlgo_EM201_2013::ajpBarrel(DSM& dsm, int offset) const
{
  int jpBits[6];

  // BC101-106

  for (int ch = 0; ch < 6; ++ch)
    {
      jpBits[ch] = dsm.channels[ch] >> offset & 0x3;
      //      printf("The channel %d jp bit is %d\n", ch, jpBits[ch]); //Test by Z. Chang
    }
  const int R3 = dsm.registers[3];
  //  printf("R3 is %d out of %d\n", R3,dsm.registers[3]);

  return (((jpBits[0] > R3) && (jpBits[1] > R3)) ||
          ((jpBits[1] > R3) && (jpBits[2] > R3)) ||
          ((jpBits[2] > R3) && (jpBits[3] > R3)) ||
          ((jpBits[3] > R3) && (jpBits[4] > R3)) ||
          ((jpBits[4] > R3) && (jpBits[5] > R3)) ||
          ((jpBits[5] > R3) && (jpBits[0] > R3)));
}
/*
int DSMAlgo_EM201_2013::ajpEndcap(const DSM& dsm) const
{
  int jpBits[6];
  const int R3 = dsm.registers[3];

  // EE101

  jpBits[0] = dsm.channels[6]      & 0x3; // JPA (4 o'clock)
  jpBits[1] = dsm.channels[6] >> 2 & 0x3; // JPB (6 o'clock)
  jpBits[2] = dsm.channels[6] >> 4 & 0x3; // JPC (8 o'clock)

  // EE102

  jpBits[3] = dsm.channels[7]      & 0x3; // JPA (10 o'clock)
  jpBits[4] = dsm.channels[7] >> 2 & 0x3; // JPB (12 o'clock)
  jpBits[5] = dsm.channels[7] >> 4 & 0x3; // JPC (2  o'clock)

  return (((jpBits[0] > R3) && (jpBits[1] > R3)) ||
          ((jpBits[1] > R3) && (jpBits[2] > R3)) ||
          ((jpBits[2] > R3) && (jpBits[3] > R3)) ||
          ((jpBits[3] > R3) && (jpBits[4] > R3)) ||
          ((jpBits[4] > R3) && (jpBits[5] > R3)) ||
          ((jpBits[5] > R3) && (jpBits[0] > R3)));
}
*/
void DSMAlgo_EM201_2013::operator()(DSM& dsm)
{
  
  // INPUT:

  // EM201 - ch0 - BEMC BC101 - 10' - JP0 & JP6 (West & East)
  //         ch1 - BEMC BC102 - 12' - JP1 & JP7
  //         ch2 - BEMC BC103 -  2' - JP2 & JP8
  //         ch3 - BEMC BC104 -  4' - JP3 & JP9
  //         ch4 - BEMC BC105 -  6' - JP4 & JP10
  //         ch5 - BEMC BC106 -  8' - JP5 & JP11
  //         ch6 - EEMC EE101 -  4',  6' and 8' - JP3, JP4 & JP5
  //         ch7 - EEMC EE102 - 10', 12' and 2' - JP0, JP1 & JP2

  // From BC101-106 (16):

  // (0-1) JPX (east, -1 < eta < 0) threshold bits (2)
  // (2-3) JPY (middle, -0.6 < eta < 0.4) threshold bits (2)
  // (4-5) JPZ (west, 0 < eta < 1) threshold bits (2)
  // (6-11) JPpartial (0.4 < eta < 1) sum (6)
  // (12-15) HT bits (4)

  // From EE101 (16):

  // (0-1) JPA (4 o'clock) threshold bits (2)
  // (2-3) JPB (6 o'clock) threshold bits (2)
  // (4-5) JPC (8 o'clock) threshold bits (2)
  // (6-11) Selected partial jet patch sum (6)
  // (12-13) Partial jet patch ID (1=A, 2=B, 3=C) (2)
  // (14-15) HT bits (2)

  // From EE102 (16):

  // (0-1) JPA (10 o'clock) threshold bits (2)
  // (2-3) JPB (12 o'clock) threshold bits (2)
  // (4-5) JPC (2  o'clock) threshold bits (2)
  // (6-11) Selected partial jet patch sum (6)
  // (12-13) Partial jet patch ID (1=A, 2=B, 3=C) (2)
  // (14-15) HT bits (2)

  // REGISTERS:

  // R0: Hybrid jet patch threshold-0
  // R1: Hybrid jet patch threshold-1
  // R2: Hybrid jet patch threshold-2
  // R3: AJP-th-Sel (2)
  // ACTION:

  // Complete hybrid jet patches using partial jet patch ID from EEMC

  int jpSum1 = dsm.channels[6] >> 6 & 0x3f; // Partial sum from EE101
  int jpSum2 = dsm.channels[7] >> 6 & 0x3f; // Partial sum from EE102

  int jpId1 = dsm.channels[6] >> 12 & 0x3; // Partial jet patch ID from EE101
  int jpId2 = dsm.channels[7] >> 12 & 0x3; // Partial jet patch ID from EE102

  switch (jpId1) {
  case 1: jpSum1 += dsm.channels[3] >> 6 & 0x3f; break; // Add partial sum from BC104 (4')
  case 2: jpSum1 += dsm.channels[4] >> 6 & 0x3f; break; // Add partial sum from BC105 (6')
  case 3: jpSum1 += dsm.channels[5] >> 6 & 0x3f; break; // Add partial sum from BC106 (8')
  }

  switch (jpId2) {
  case 1: jpSum2 += dsm.channels[0] >> 6 & 0x3f; break; // Add partial sum from BC101 (10')
  case 2: jpSum2 += dsm.channels[1] >> 6 & 0x3f; break; // Add partial sum from BC102 (12')
  case 3: jpSum2 += dsm.channels[2] >> 6 & 0x3f; break; // Add partial sum from BC103 (2')
  }

  // Combine (OR) the HT bits from the six BEMC layer 1 DSM's

  int htBitsBarrel = 0;
  //0x3 for combined HT bits
  for (int ch = 0; ch < 6; ++ch){
    int packedHT  = dsm.channels[ch] >> 12 & 0x3;
    int unpackedHT = dsm.channels[ch] >> 14 & 0x1;
    htBitsBarrel |= ((1 << packedHT) - 1) | (unpackedHT << 3);
  }

  // Combine (OR) the HT bits from the two EEMC layer 1 DSM's

  int htBitsEndcap = 0;

  for (int ch = 6; ch < 8; ++ch)
    htBitsEndcap |= dsm.channels[ch] >> 14 & 0x3;

  // Combine (OR) the JP bits for the BEMC and EEMC separately

  int jpBitsBarrel = 0;

  for (int ch = 0; ch < 6; ++ch) {
    int jpx = dsm.channels[ch]      & 0x3;
    int jpy = dsm.channels[ch] >> 2 & 0x3;
    int jpz = dsm.channels[ch] >> 4 & 0x3;

    if (jpx > jpBitsBarrel) jpBitsBarrel = jpx;
    if (jpy > jpBitsBarrel) jpBitsBarrel = jpy;
    if (jpz > jpBitsBarrel) jpBitsBarrel = jpz;
  }

  int bjp1 = jpBitsBarrel > 1;
//  int bjp2 = jpBitsBarrel > 2;

  int jpBitsEndcap = 0;

  for (int ch = 6; ch < 8; ++ch) {
    int jpabc = dsm.channels[ch]      & 0x3;
    //    int jpa = dsm.channels[ch]      & 0x3;
    //    int jpb = dsm.channels[ch] >> 2 & 0x3;
    //    int jpc = dsm.channels[ch] >> 4 & 0x3;
    if (jpabc > jpBitsEndcap) jpBitsEndcap = jpabc;
    //    if (jpa > jpBitsEndcap) jpBitsEndcap = jpa;
    //    if (jpb > jpBitsEndcap) jpBitsEndcap = jpb;
    //    if (jpc > jpBitsEndcap) jpBitsEndcap = jpc;
  }

  int ejp1 = jpBitsEndcap > 1;
//  int ejp2 = jpBitsEndcap > 2;

  // Compare the two completed hybrid jet patches to three thresholds
  // and combine (OR) the results with the BEMC-only and EEMC-only bits

  int jpBits = 0;

  for (int reg = 0; reg < 3; ++reg)
    if (jpSum1 > dsm.registers[reg] || jpSum2 > dsm.registers[reg]) ++jpBits;
  //printf("em201: r0 = %d r1 = %d r2 = %d\n", dsm.registers[0], dsm.registers[1], dsm.registers[2]);
  if (jpBitsBarrel > jpBits) jpBits = jpBitsBarrel;
  if (jpBitsEndcap > jpBits) jpBits = jpBitsEndcap;

  int jp0 = jpBits > 0;
  int jp1 = jpBits > 1;
  int jp2 = jpBits > 2;

  // Adjacent jet patch logic

  //printf("R3 = %d\n", dsm.registers[3]);
  int ajpx = ajpBarrel(dsm, 0);
  int ajpy = ajpBarrel(dsm, 2);
  int ajpz = ajpBarrel(dsm, 4);
  int bajp = ajpx || ajpy || ajpz;

  int dijet0 = Dijet(dsm, 0);

  //printf("dijet0 = %d\n", dijet0);

  int dijet1 = Dijet(dsm, 1);

  //printf("dijet1 = %d\n", dijet1);

  int edijet = EndDijet(dsm, 0);
  //printf("edijet = %d\n", edijet);

  int daq10k = 0;
  //No endcap ajp in 2013
  //  int eajp = ajpEndcap(dsm);
  //  int  ajp = bajp || eajp;

  // OUTPUT (16):

  // (0:3) Barrel HT bits (4)
  // (4:5) Endcap HT bits (2)
  // (6) JP1, unified over the BEMC+EEMC (1)
  // (7) JP2, unified over the BEMC+EEMC (1)
  // (8) BJP1, for the 18 BEMC-only patches (1)
  // (9) EEMC-dijet, EEMC-only JP0-based dijet trigger bit (1)
  // (10) EJP1, for the 6 EEMC-only patches (1)
  // (11) JP1-dijet, JP1-based dijet trigger bit (1)
  // (12) JP0-dijet, JP0-based dijet trigger bit (1)
  // (13) BAJP, for the BEMC-only patches (1)
  // (14) DAQ10k, DAQ10k trigger bit (1)
  // (15) JP0, unified over the BEMC+EEMC (1)

  int out = 0;

  out |= htBitsBarrel;
  out |= htBitsEndcap << 4;
  out |= jp1  << 6;
  out |= jp2  << 7;
  out |= bjp1 << 8;
  out |= edijet << 9;
  out |= ejp1 << 10;
  out |= dijet1 << 11;
  out |= dijet0  << 12;
  out |= bajp << 13;
  out |= daq10k << 14;
  out |= jp0  << 15;

  dsm.output = out;

  // INFO

  dsm.info[0] = jpSum1;
  dsm.info[1] = jpSum2;

//  return 0;
}
int DSMAlgo_EM201_2013::Dijet(DSM &dsm, int jpTH = 0) const
{
  // Search for BEMC JP0 & JP1 di-jet signature

  int eemcDijetBits[6];

  eemcDijetBits[0] = dsm.channels[7] >> 2 & 1; // EE102 10 o'clock
  eemcDijetBits[1] = dsm.channels[7] >> 3 & 1; // EE102 12 o'clock
  eemcDijetBits[2] = dsm.channels[7] >> 4 & 1; // EE102  2 o'clock
  eemcDijetBits[3] = dsm.channels[6] >> 2 & 1; // EE101  4 o'clock
  eemcDijetBits[4] = dsm.channels[6] >> 3 & 1; // EE101  6 o'clock
  eemcDijetBits[5] = dsm.channels[6] >> 4 & 1; // EE101  8 o'clock

  int dijet = 0;
  //  int dijet1 = 0;

  for (int ch = 0; ch < 6; ++ch) {
    int jpx = dsm.channels[ch] >> 0 & 0x3;
    int jpy = dsm.channels[ch] >> 2 & 0x3;
    int jpz = dsm.channels[ch] >> 4 & 0x3;

    int bemcDijet2 = dsm.channels[(ch+2)%6] >> 15 & 1;
    int bemcDijet3 = dsm.channels[(ch+3)%6] >> 15 & 1;
    int bemcDijet4 = dsm.channels[(ch+4)%6] >> 15 & 1;

    int eemcDijet2 = eemcDijetBits[(ch+2)%6];
    int eemcDijet3 = eemcDijetBits[(ch+3)%6];
    int eemcDijet4 = eemcDijetBits[(ch+4)%6];

    int bjp = jpx > jpTH || jpy > jpTH || jpz > jpTH;
    //    int bjp = jpx > 1 || jpy > 1 || jpz > 1;

    int bemcDijet = bemcDijet2 || bemcDijet3 || bemcDijet4;
    int eemcDijet = eemcDijet2 || eemcDijet3 || eemcDijet4;

    dijet |= bjp && (bemcDijet || eemcDijet);
    //    dijet1 |= bjp1 && (bemcDijet || eemcDijet);
  }
  return dijet;
}
int DSMAlgo_EM201_2013::EndDijet(DSM &dsm, int jpTH = 0) const
{
  // Search for JP0 EEMC di-jet signature
  int eemcDijetBits[6];

  eemcDijetBits[0] = dsm.channels[7] >> 2 & 1; // EE102 10 o'clock
  eemcDijetBits[1] = dsm.channels[7] >> 3 & 1; // EE102 12 o'clock
  eemcDijetBits[2] = dsm.channels[7] >> 4 & 1; // EE102  2 o'clock
  eemcDijetBits[3] = dsm.channels[6] >> 2 & 1; // EE101  4 o'clock
  eemcDijetBits[4] = dsm.channels[6] >> 3 & 1; // EE101  6 o'clock
  eemcDijetBits[5] = dsm.channels[6] >> 4 & 1; // EE101  8 o'clock

  int ee101_jp0 = dsm.channels[6] & 3; // 4, 6 and 8 o'clock
  int ee102_jp0 = dsm.channels[7] & 3; // 10, 12 and 2 o'clock

  int ee101_dijet = ee101_jp0 && (eemcDijetBits[0] || eemcDijetBits[1] || eemcDijetBits[2]);
  int ee102_dijet = ee102_jp0 && (eemcDijetBits[3] || eemcDijetBits[4] || eemcDijetBits[5]);

  int eemc_dijet = ee101_dijet || ee102_dijet;

  return eemc_dijet;
}