HcalFeatureHFEMBit.cc

Go to the documentation of this file.

 
//HcalFeatureHFEMBit
//version 2.0
 
#include "SimCalorimetry/HcalTrigPrimAlgos/interface/HcalFeatureHFEMBit.h"
#include "CondFormats/HcalObjects/interface/HcalQIECoder.h"
 
HcalFeatureHFEMBit::HcalFeatureHFEMBit(double ShortMinE,
                                       double LongMinE,
                                       double ShortLongCutSlope,
                                       double ShortLongCutOffset,
                                       const HcalDbService& conditions)
    : conditions_(conditions) {
  ShortMinE_ = ShortMinE;  //minimum energy deposited
  LongMinE_ = LongMinE;
  ShortLongCutSlope_ =
      ShortLongCutSlope;  // this is a the slope of the cut line related to energy deposited in short fibers vrs long fibers
  ShortLongCutOffset_ = ShortLongCutOffset;  // this is the offset of said line.
}
 
HcalFeatureHFEMBit::~HcalFeatureHFEMBit() {}
 
bool HcalFeatureHFEMBit::fineGrainbit(const HFDataFrame& shortDigi, const HFDataFrame& longDigi, int idx) const {
  float shortE = getE(shortDigi, idx);
  float longE = getE(longDigi, idx);
 
  if (shortE < ShortMinE_)
    return false;
  if (longE < LongMinE_)
    return false;
 
  return (shortE < (longE - ShortLongCutOffset_) * ShortLongCutSlope_);
}
 
bool HcalFeatureHFEMBit::fineGrainbit(const QIE10DataFrame& short1,
                                      const QIE10DataFrame& short2,
                                      const QIE10DataFrame& long1,
                                      const QIE10DataFrame& long2,
                                      bool validShort1,
                                      bool validShort2,
                                      bool validLong1,
                                      bool validLong2,
                                      int idx) const {
  float shortE = 0;
  if (validShort1)
    shortE += getE(short1, idx);
  if (validShort2)
    shortE += getE(short2, idx);
  if (validShort1 and validShort2)
    shortE *= .5;
 
  float longE = 0;
  if (validLong1)
    longE += getE(long1, idx);
  if (validLong2)
    longE += getE(long2, idx);
  if (validLong1 and validLong2)
    longE *= .5;
 
  if (shortE < ShortMinE_)
    return false;
  if (longE < LongMinE_)
    return false;
 
  return (shortE < (longE - ShortLongCutOffset_) * ShortLongCutSlope_);
}