CMS 3D CMS Logo
List of all members | Classes | Public Member Functions | Private Types | Private Member Functions | Private Attributes | Static Private Attributes
HcalTriggerPrimitiveAlgo Class Reference

#include <HcalTriggerPrimitiveAlgo.h>

Classes

struct  HFDetails
 
struct  HFUpgradeDetails
 

Public Member Functions

template<typename T , typename... Args>
void addDigis (const T &collection, const Args &... digis)
 
template<typename T >
void addDigis (const T &collection)
 
template<typename D >
void addDigis (const HcalDataFrameContainer< D > &collection)
 
 HcalTriggerPrimitiveAlgo (bool pf, const std::vector< double > &w, int latency, uint32_t FG_threshold, const std::vector< uint32_t > &FG_HF_thresholds, uint32_t ZS_threshold, int numberOfSamples, int numberOfPresamples, int numberOfFilterPresamplesHBQIE11, int numberOfFilterPresamplesHEQIE11, int numberOfSamplesHF, int numberOfPresamplesHF, bool useTDCInMinBiasBits, uint32_t minSignalThreshold=0, uint32_t PMT_NoiseThreshold=0)
 
void overrideParameters (const edm::ParameterSet &ps)
 
template<typename... Digis>
void run (const HcalTPGCoder *incoder, const HcalTPGCompressor *outcoder, const HcalDbService *conditions, HcalTrigPrimDigiCollection &result, const HcalTrigTowerGeometry *trigTowerGeometry, float rctlsb, const HcalFeatureBit *LongvrsShortCut, const Digis &... digis)
 
void runFEFormatError (const FEDRawDataCollection *rawraw, const HcalElectronicsMap *emap, HcalTrigPrimDigiCollection &result)
 
void runZS (HcalTrigPrimDigiCollection &tp)
 
void setFixSaturationFlag (bool fix_saturation)
 
void setNCTScaleShift (int)
 
void setNumFilterPresamplesHBQIE11 (int presamples)
 
void setNumFilterPresamplesHEQIE11 (int presamples)
 
void setPeakFinderAlgorithm (int algo)
 
void setRCTScaleShift (int)
 
void setUpgradeFlags (bool hb, bool he, bool hf)
 
void setWeightQIE11 (int aieta, int weight)
 
void setWeightsQIE11 (const edm::ParameterSet &weightsQIE11)
 
 ~HcalTriggerPrimitiveAlgo ()
 

Private Types

typedef std::map< HcalTrigTowerDetId, std::vector< bool > > FGbitMap
 
typedef std::vector< HcalFinegrainBit::Tower > FGUpgradeContainer
 
typedef std::map< HcalTrigTowerDetId, FGUpgradeContainer > FGUpgradeMap
 
typedef std::vector< HcalFinegrainBit::TowerTDC > FGUpgradeTDCContainer
 
typedef std::map< HcalTrigTowerDetId, FGUpgradeTDCContainer > FGUpgradeTDCMap
 
typedef std::map< HcalTrigTowerDetId, std::map< uint32_t, HFDetails > > HFDetailMap
 
typedef std::map< HcalTrigTowerDetId, std::map< uint32_t, std::array< HFUpgradeDetails, 4 > > > HFUpgradeDetailMap
 
typedef std::map< HcalTrigTowerDetId, std::vector< bool > > SatMap
 
typedef std::vector< IntegerCaloSamples > SumFGContainer
 
typedef std::map< HcalTrigTowerDetId, IntegerCaloSamples > SumMap
 
typedef std::map< HcalTrigTowerDetId, SumFGContainer > TowerMapFGSum
 
typedef std::map< uint32_t, std::vector< bool > > TowerMapVeto
 

Private Member Functions

void addFG (const HcalTrigTowerDetId &id, std::vector< bool > &msb)
 
void addSignal (const HBHEDataFrame &frame)
 adds the signal to the map More...
 
void addSignal (const HFDataFrame &frame)
 
void addSignal (const QIE10DataFrame &frame)
 
void addSignal (const QIE11DataFrame &frame)
 
void addSignal (const IntegerCaloSamples &samples)
 
void addUpgradeFG (const HcalTrigTowerDetId &id, int depth, const std::vector< std::bitset< 2 >> &bits)
 
void addUpgradeTDCFG (const HcalTrigTowerDetId &id, const QIE11DataFrame &frame)
 
void analyze (IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result)
 adds the actual digis More...
 
void analyzeHF (IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result, const int hf_lumi_shift)
 
void analyzeHF2016 (const IntegerCaloSamples &SAMPLES, HcalTriggerPrimitiveDigi &result, const int HF_LUMI_SHIFT, const HcalFeatureBit *HCALFEM)
 
void analyzeHFQIE10 (const IntegerCaloSamples &SAMPLES, HcalTriggerPrimitiveDigi &result, const int HF_LUMI_SHIFT, const HcalFeatureBit *HCALFEM)
 
void analyzeQIE11 (IntegerCaloSamples &samples, std::vector< bool > sample_saturation, HcalTriggerPrimitiveDigi &result, const HcalFinegrainBit &fg_algo)
 
bool needLegacyFG (const HcalTrigTowerDetId &id) const
 
bool needUpgradeID (const HcalTrigTowerDetId &id, int depth) const
 
bool passTDC (const QIE10DataFrame &digi, int ts) const
 
bool validChannel (const QIE10DataFrame &digi, int ts) const
 
bool validUpgradeFG (const HcalTrigTowerDetId &id, int depth) const
 

Private Attributes

const HcalDbService * conditions_
 
std::vector< uint32_t > FG_HF_thresholds_
 
uint32_t FG_threshold_
 
FGbitMap fgMap_
 
FGUpgradeMap fgUpgradeMap_
 
FGUpgradeTDCMap fgUpgradeTDCMap_
 
bool fix_saturation_ = false
 
TowerMapVeto HF_Veto
 
const HcaluLUTTPGCoder * incoder_
 
int latency_
 
HcalFeatureBit * LongvrsShortCut
 
uint32_t minSignalThreshold_
 
int NCTScaleShift
 
int numberOfFilterPresamplesHBQIE11_
 
int numberOfFilterPresamplesHEQIE11_
 
int numberOfPresamples_
 
int numberOfPresamplesHF_
 
int numberOfSamples_
 
int numberOfSamplesHF_
 
const HcalTPGCompressor * outcoder_
 
bool override_adc_hf_ = false
 
uint32_t override_adc_hf_value_
 
edm::ParameterSet override_parameters_
 
bool override_tdc_hf_ = false
 
unsigned long long override_tdc_hf_value_
 
int peak_finder_algorithm_
 
bool peakfind_
 
uint32_t PMT_NoiseThreshold_
 
int RCTScaleShift
 
HFDetailMap theHFDetailMap
 
HFUpgradeDetailMap theHFUpgradeDetailMap
 
SatMap theSatMap
 
SumMap theSumMap
 
double theThreshold
 
TowerMapFGSum theTowerMapFGSum
 
const HcalTrigTowerGeometry * theTrigTowerGeometry
 
bool upgrade_hb_ = false
 
bool upgrade_he_ = false
 
bool upgrade_hf_ = false
 
bool useTDCInMinBiasBits_
 
std::vector< double > weights_
 
std::array< std::array< int, 2 >, 29 > weightsQIE11_
 
uint32_t ZS_threshold_
 
int ZS_threshold_I_
 

Static Private Attributes

static const int FIRST_DEPTH7_TOWER = 26
 
static const int FIRST_FINEGRAIN_TOWER = 30
 
static const int HBHE_OVERLAP_TOWER = 16
 
static const int LAST_FINEGRAIN_DEPTH = 6
 
static const int LAST_FINEGRAIN_TOWER = 28
 
static const int QIE10_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE10_LUT_BITMASK
 
static const int QIE10_MAX_LINEARIZATION_ET = 0x7FF
 
static const int QIE11_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE11_LUT_BITMASK
 
static const int QIE11_MAX_LINEARIZATION_ET = 0x7FF
 
static const int QIE8_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE8_LUT_BITMASK
 

Detailed Description

Definition at line 27 of file HcalTriggerPrimitiveAlgo.h.

Member Typedef Documentation

◆ FGbitMap

typedef std::map<HcalTrigTowerDetId, std::vector<bool> > HcalTriggerPrimitiveAlgo::FGbitMap
private

Definition at line 210 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeContainer

typedef std::vector<HcalFinegrainBit::Tower> HcalTriggerPrimitiveAlgo::FGUpgradeContainer
private

Definition at line 213 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeMap

typedef std::map<HcalTrigTowerDetId, FGUpgradeContainer> HcalTriggerPrimitiveAlgo::FGUpgradeMap
private

Definition at line 214 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeTDCContainer

typedef std::vector<HcalFinegrainBit::TowerTDC> HcalTriggerPrimitiveAlgo::FGUpgradeTDCContainer
private

Definition at line 217 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeTDCMap

typedef std::map<HcalTrigTowerDetId, FGUpgradeTDCContainer> HcalTriggerPrimitiveAlgo::FGUpgradeTDCMap
private

Definition at line 218 of file HcalTriggerPrimitiveAlgo.h.

◆ HFDetailMap

typedef std::map<HcalTrigTowerDetId, std::map<uint32_t, HFDetails> > HcalTriggerPrimitiveAlgo::HFDetailMap
private

Definition at line 180 of file HcalTriggerPrimitiveAlgo.h.

◆ HFUpgradeDetailMap

typedef std::map<HcalTrigTowerDetId, std::map<uint32_t, std::array<HFUpgradeDetails, 4> > > HcalTriggerPrimitiveAlgo::HFUpgradeDetailMap
private

Definition at line 190 of file HcalTriggerPrimitiveAlgo.h.

◆ SatMap

typedef std::map<HcalTrigTowerDetId, std::vector<bool> > HcalTriggerPrimitiveAlgo::SatMap
private

Definition at line 171 of file HcalTriggerPrimitiveAlgo.h.

◆ SumFGContainer

typedef std::vector<IntegerCaloSamples> HcalTriggerPrimitiveAlgo::SumFGContainer
private

Definition at line 193 of file HcalTriggerPrimitiveAlgo.h.

◆ SumMap

typedef std::map<HcalTrigTowerDetId, IntegerCaloSamples> HcalTriggerPrimitiveAlgo::SumMap
private

Definition at line 168 of file HcalTriggerPrimitiveAlgo.h.

◆ TowerMapFGSum

typedef std::map<HcalTrigTowerDetId, SumFGContainer> HcalTriggerPrimitiveAlgo::TowerMapFGSum
private

Definition at line 194 of file HcalTriggerPrimitiveAlgo.h.

◆ TowerMapVeto

typedef std::map<uint32_t, std::vector<bool> > HcalTriggerPrimitiveAlgo::TowerMapVeto
private

Definition at line 207 of file HcalTriggerPrimitiveAlgo.h.

Constructor & Destructor Documentation

◆ HcalTriggerPrimitiveAlgo()

HcalTriggerPrimitiveAlgo::HcalTriggerPrimitiveAlgo ( bool  pf,
const std::vector< double > &  w,
int  latency,
uint32_t  FG_threshold,
const std::vector< uint32_t > &  FG_HF_thresholds,
uint32_t  ZS_threshold,
int  numberOfSamples,
int  numberOfPresamples,
int  numberOfFilterPresamplesHBQIE11,
int  numberOfFilterPresamplesHEQIE11,
int  numberOfSamplesHF,
int  numberOfPresamplesHF,
bool  useTDCInMinBiasBits,
uint32_t  minSignalThreshold = 0,
uint32_t  PMT_NoiseThreshold = 0 
)

Definition at line 23 of file HcalTriggerPrimitiveAlgo.cc.

References numberOfPresamples_, numberOfPresamplesHF_, numberOfSamples_, numberOfSamplesHF_, peakfind_, ZS_threshold_, and ZS_threshold_I_.

38  : incoder_(nullptr),
39  outcoder_(nullptr),
40  theThreshold(0),
41  peakfind_(pf),
42  weights_(w),
43  latency_(latency),
44  FG_threshold_(FG_threshold),
45  FG_HF_thresholds_(FG_HF_thresholds),
46  ZS_threshold_(ZS_threshold),
47  numberOfSamples_(numberOfSamples),
48  numberOfPresamples_(numberOfPresamples),
49  numberOfFilterPresamplesHBQIE11_(numberOfFilterPresamplesHBQIE11),
50  numberOfFilterPresamplesHEQIE11_(numberOfFilterPresamplesHEQIE11),
51  numberOfSamplesHF_(numberOfSamplesHF),
52  numberOfPresamplesHF_(numberOfPresamplesHF),
53  useTDCInMinBiasBits_(useTDCInMinBiasBits),
54  minSignalThreshold_(minSignalThreshold),
55  PMT_NoiseThreshold_(PMT_NoiseThreshold),
56  NCTScaleShift(0),
57  RCTScaleShift(0),
58  peak_finder_algorithm_(2),
59  override_parameters_() {
60  //No peak finding setting (for Fastsim)
61  if (!peakfind_) {
62  numberOfSamples_ = 1;
63  numberOfPresamples_ = 0;
64  numberOfSamplesHF_ = 1;
65  numberOfPresamplesHF_ = 0;
66  }
67  // Switch to integer for comparisons - remove compiler warning
68  ZS_threshold_I_ = ZS_threshold_;
69 }
HcalTriggerPrimitiveAlgo::FG_HF_thresholds_
std::vector< uint32_t > FG_HF_thresholds_
Definition: HcalTriggerPrimitiveAlgo.h:142
w
T w() const
Definition: extBasic3DVector.h:225
config_102169_raw_cff.latency
latency
Definition: config_102169_raw_cff.py:29
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
hcaltpdigi_cfi.numberOfFilterPresamplesHEQIE11
numberOfFilterPresamplesHEQIE11
Definition: hcaltpdigi_cfi.py:56
packedPFCandidateRefMixer_cfi.pf
pf
Definition: packedPFCandidateRefMixer_cfi.py:4
hcaltpdigi_cfi.numberOfPresamples
numberOfPresamples
Definition: hcaltpdigi_cfi.py:52
hcaltpdigi_cfi.numberOfPresamplesHF
numberOfPresamplesHF
Definition: hcaltpdigi_cfi.py:54
hcaltpdigi_cfi.numberOfFilterPresamplesHBQIE11
numberOfFilterPresamplesHBQIE11
Definition: hcaltpdigi_cfi.py:55
hcal_dqm_sourceclient-live_cfg.ZS_threshold
ZS_threshold
Definition: hcal_dqm_sourceclient-live_cfg.py:107
hcaltpdigi_cff.FG_HF_thresholds
FG_HF_thresholds
Definition: hcaltpdigi_cff.py:27
hcaltpdigi_cfi.useTDCInMinBiasBits
useTDCInMinBiasBits
Definition: hcaltpdigi_cfi.py:57
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133
hcal_dqm_sourceclient-live_cfg.FG_threshold
FG_threshold
Definition: hcal_dqm_sourceclient-live_cfg.py:99
hcaltpdigi_cfi.numberOfSamplesHF
numberOfSamplesHF
Definition: hcaltpdigi_cfi.py:53

◆ ~HcalTriggerPrimitiveAlgo()

HcalTriggerPrimitiveAlgo::~HcalTriggerPrimitiveAlgo ( )

Definition at line 71 of file HcalTriggerPrimitiveAlgo.cc.

71 {}

Member Function Documentation

◆ addDigis() [1/3]

template<typename T , typename... Args>
void HcalTriggerPrimitiveAlgo::addDigis ( const T &  collection,
const Args &...  digis 
)
inline

Definition at line 57 of file HcalTriggerPrimitiveAlgo.h.

References universalConfigTemplate::collection.

Referenced by run().

57  {
58  addDigis(collection);
59  addDigis(digis...);
60  };
HcalTriggerPrimitiveAlgo::addDigis
void addDigis(const T &collection, const Args &... digis)
Definition: HcalTriggerPrimitiveAlgo.h:57

◆ addDigis() [2/3]

template<typename T >
void HcalTriggerPrimitiveAlgo::addDigis ( const T &  collection)
inline

Definition at line 63 of file HcalTriggerPrimitiveAlgo.h.

References addSignal(), and universalConfigTemplate::collection.

63  {
64  for (const auto& digi : collection) {
65  addSignal(digi);
66  }
67  };
HcalTriggerPrimitiveAlgo::addSignal
void addSignal(const HBHEDataFrame &frame)
adds the signal to the map
Definition: HcalTriggerPrimitiveAlgo.cc:94

◆ addDigis() [3/3]

template<typename D >
void HcalTriggerPrimitiveAlgo::addDigis ( const HcalDataFrameContainer< D > &  collection)
inline

Definition at line 70 of file HcalTriggerPrimitiveAlgo.h.

References addSignal(), universalConfigTemplate::collection, and mps_fire::i.

70  {
71  for (auto i = collection.begin(); i != collection.end(); ++i) {
72  D digi(*i);
73  addSignal(digi);
74  }
75  };
mps_fire.i
i
Definition: mps_fire.py:429
funct::D
DecomposeProduct< arg, typename Div::arg > D
Definition: Factorize.h:141
HcalTriggerPrimitiveAlgo::addSignal
void addSignal(const HBHEDataFrame &frame)
adds the signal to the map
Definition: HcalTriggerPrimitiveAlgo.cc:94

◆ addFG()

void HcalTriggerPrimitiveAlgo::addFG ( const HcalTrigTowerDetId &  id,
std::vector< bool > &  msb 
)
private

Definition at line 890 of file HcalTriggerPrimitiveAlgo.cc.

References fgMap_, mps_fire::i, and l1ctLayer2EG_cff::id.

Referenced by addSignal(), and addUpgradeFG().

890  {
891  FGbitMap::iterator itr = fgMap_.find(id);
892  if (itr != fgMap_.end()) {
893  std::vector<bool>& _msb = itr->second;
894  for (size_t i = 0; i < msb.size(); ++i)
895  _msb[i] = _msb[i] || msb[i];
896  } else
897  fgMap_[id] = msb;
898 }
mps_fire.i
i
Definition: mps_fire.py:429
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:85

◆ addSignal() [1/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const HBHEDataFrame &  frame)
private

adds the signal to the map

Definition at line 94 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addFG(), cms::cuda::assert(), mps_fire::i, incoder_, HcaluLUTTPGCoder::lookupMSB(), theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

Referenced by addDigis(), and addSignal().

94  {
95  // TODO: Need to add support for seperate 28, 29 in HE
96  //Hack for 300_pre10, should be removed.
97  if (frame.id().depth() == 5)
98  return;
99 
100  std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
101  assert(ids.size() == 1 || ids.size() == 2);
102  IntegerCaloSamples samples1(ids[0], int(frame.size()));
103 
104  samples1.setPresamples(frame.presamples());
105  incoder_->adc2Linear(frame, samples1);
106 
107  std::vector<bool> msb;
108  incoder_->lookupMSB(frame, msb);
109 
110  if (ids.size() == 2) {
111  // make a second trigprim for the other one, and split the energy
112  IntegerCaloSamples samples2(ids[1], samples1.size());
113  for (int i = 0; i < samples1.size(); ++i) {
114  samples1[i] = uint32_t(samples1[i] * 0.5);
115  samples2[i] = samples1[i];
116  }
117  samples2.setPresamples(frame.presamples());
118  addSignal(samples2);
119  addFG(ids[1], msb);
120  }
121  addSignal(samples1);
122  addFG(ids[0], msb);
123 }
mps_fire.i
i
Definition: mps_fire.py:429
HcaluLUTTPGCoder::lookupMSB
void lookupMSB(const HBHEDataFrame &df, std::vector< bool > &msb) const
Definition: HcaluLUTTPGCoder.cc:629
HcalTriggerPrimitiveAlgo::addFG
void addFG(const HcalTrigTowerDetId &id, std::vector< bool > &msb)
Definition: HcalTriggerPrimitiveAlgo.cc:890
cms::cuda::assert
assert(be >=bs)
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
HcaluLUTTPGCoder::adc2Linear
void adc2Linear(const HBHEDataFrame &df, IntegerCaloSamples &ics) const override
Definition: HcaluLUTTPGCoder.cc:565
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133
HcalTriggerPrimitiveAlgo::addSignal
void addSignal(const HBHEDataFrame &frame)
adds the signal to the map
Definition: HcalTriggerPrimitiveAlgo.cc:94
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15

◆ addSignal() [2/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const HFDataFrame &  frame)
private

Definition at line 125 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), amptDefault_cfi::frame, HF_Veto, mps_fire::i, incoder_, minSignalThreshold_, EgammaValidation_cff::samples, IntegerCaloSamples::setPresamples(), theHFDetailMap, theTowerMapFGSum, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), and HcalTrigTowerDetId::version().

125  {
126  if (frame.id().depth() == 1 || frame.id().depth() == 2) {
127  std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
128  std::vector<HcalTrigTowerDetId>::const_iterator it;
129  for (it = ids.begin(); it != ids.end(); ++it) {
130  HcalTrigTowerDetId trig_tower_id = *it;
131  IntegerCaloSamples samples(trig_tower_id, frame.size());
132  samples.setPresamples(frame.presamples());
133  incoder_->adc2Linear(frame, samples);
134 
135  // Don't add to final collection yet
136  // HF PMT veto sum is calculated in analyzerHF()
137  IntegerCaloSamples zero_samples(trig_tower_id, frame.size());
138  zero_samples.setPresamples(frame.presamples());
139  addSignal(zero_samples);
140 
141  // Pre-LS1 Configuration
142  if (trig_tower_id.version() == 0) {
143  // Mask off depths: fgid is the same for both depths
144  uint32_t fgid = (frame.id().maskDepth());
145 
146  if (theTowerMapFGSum.find(trig_tower_id) == theTowerMapFGSum.end()) {
147  SumFGContainer sumFG;
148  theTowerMapFGSum.insert(std::pair<HcalTrigTowerDetId, SumFGContainer>(trig_tower_id, sumFG));
149  }
150 
151  SumFGContainer& sumFG = theTowerMapFGSum[trig_tower_id];
152  SumFGContainer::iterator sumFGItr;
153  for (sumFGItr = sumFG.begin(); sumFGItr != sumFG.end(); ++sumFGItr) {
154  if (sumFGItr->id() == fgid) {
155  break;
156  }
157  }
158  // If find
159  if (sumFGItr != sumFG.end()) {
160  for (int i = 0; i < samples.size(); ++i) {
161  (*sumFGItr)[i] += samples[i];
162  }
163  } else {
164  //Copy samples (change to fgid)
165  IntegerCaloSamples sumFGSamples(DetId(fgid), samples.size());
166  sumFGSamples.setPresamples(samples.presamples());
167  for (int i = 0; i < samples.size(); ++i) {
168  sumFGSamples[i] = samples[i];
169  }
170  sumFG.push_back(sumFGSamples);
171  }
172 
173  // set veto to true if Long or Short less than threshold
174  if (HF_Veto.find(fgid) == HF_Veto.end()) {
175  vector<bool> vetoBits(samples.size(), false);
176  HF_Veto[fgid] = vetoBits;
177  }
178  for (int i = 0; i < samples.size(); ++i) {
179  if (samples[i] < minSignalThreshold_) {
180  HF_Veto[fgid][i] = true;
181  }
182  }
183  }
184  // HF 1x1
185  else if (trig_tower_id.version() == 1) {
186  uint32_t fgid = (frame.id().maskDepth());
187  HFDetails& details = theHFDetailMap[trig_tower_id][fgid];
188  // Check the frame type to determine long vs short
189  if (frame.id().depth() == 1) { // Long
190  details.long_fiber = samples;
191  details.LongDigi = frame;
192  } else if (frame.id().depth() == 2) { // Short
193  details.short_fiber = samples;
194  details.ShortDigi = frame;
195  } else {
196  // Neither long nor short... So we have no idea what to do
197  edm::LogWarning("HcalTPAlgo") << "Unable to figure out what to do with data frame for " << frame.id();
198  return;
199  }
200  }
201  // Uh oh, we are in a bad/unknown state! Things will start crashing.
202  else {
203  return;
204  }
205  }
206  }
207 }
amptDefault_cfi.frame
frame
Definition: amptDefault_cfi.py:12
mps_fire.i
i
Definition: mps_fire.py:429
IntegerCaloSamples::setPresamples
void setPresamples(int pre)
set presample information
Definition: IntegerCaloSamples.cc:13
HcalTriggerPrimitiveAlgo::SumFGContainer
std::vector< IntegerCaloSamples > SumFGContainer
Definition: HcalTriggerPrimitiveAlgo.h:193
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
details
Definition: helper.h:56
DetId
Definition: DetId.h:17
HcaluLUTTPGCoder::adc2Linear
void adc2Linear(const HBHEDataFrame &df, IntegerCaloSamples &ics) const override
Definition: HcaluLUTTPGCoder.cc:565
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:128
HcalTrigTowerDetId::version
int version() const
get the version code for the trigger tower
Definition: HcalTrigTowerDetId.h:60
HcalTriggerPrimitiveAlgo::addSignal
void addSignal(const HBHEDataFrame &frame)
adds the signal to the map
Definition: HcalTriggerPrimitiveAlgo.cc:94
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15

◆ addSignal() [3/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE10DataFrame &  frame)
private

Definition at line 209 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), hcalRecHitTable_cff::detId, amptDefault_cfi::frame, HcalForward, l1ctLayer2EG_cff::id, heavyIonCSV_trainingSettings::idx, incoder_, HcaluLUTTPGCoder::lookupMSB(), passTDC(), EgammaValidation_cff::samples, IntegerCaloSamples::setPresamples(), theHFUpgradeDetailMap, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), validChannel(), and BeamSplash_cfg::version.

209  {
210  HcalDetId detId = frame.detid();
211  // prevent QIE10 calibration channels from entering TP emulation
212  if (detId.subdet() != HcalForward)
213  return;
214 
215  auto ids = theTrigTowerGeometry->towerIds(frame.id());
216  for (const auto& id : ids) {
217  if (id.version() == 0) {
218  edm::LogError("HcalTPAlgo") << "Encountered QIE10 data frame mapped to TP version 0:" << id;
219  continue;
220  }
221 
222  int nsamples = frame.samples();
223 
224  IntegerCaloSamples samples(id, nsamples);
225  samples.setPresamples(frame.presamples());
226  incoder_->adc2Linear(frame, samples);
227 
228  // Don't add to final collection yet
229  // HF PMT veto sum is calculated in analyzerHF()
230  IntegerCaloSamples zero_samples(id, nsamples);
231  zero_samples.setPresamples(frame.presamples());
232  addSignal(zero_samples);
233 
234  auto fid = HcalDetId(frame.id());
235  auto& details = theHFUpgradeDetailMap[id][fid.maskDepth()];
236  auto& detail = details[fid.depth() - 1];
237  detail.samples = samples;
238  detail.digi = frame;
239  detail.validity.resize(nsamples);
240  detail.passTDC.resize(nsamples);
241  incoder_->lookupMSB(frame, detail.fgbits);
242  for (int idx = 0; idx < nsamples; ++idx) {
243  detail.validity[idx] = validChannel(frame, idx);
244  detail.passTDC[idx] = passTDC(frame, idx);
245  }
246  }
247 }
amptDefault_cfi.frame
frame
Definition: amptDefault_cfi.py:12
HcalTriggerPrimitiveAlgo::theHFUpgradeDetailMap
HFUpgradeDetailMap theHFUpgradeDetailMap
Definition: HcalTriggerPrimitiveAlgo.h:191
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
HcaluLUTTPGCoder::lookupMSB
void lookupMSB(const HBHEDataFrame &df, std::vector< bool > &msb) const
Definition: HcaluLUTTPGCoder.cc:629
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:85
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
details
Definition: helper.h:56
BeamSplash_cfg.version
version
Definition: BeamSplash_cfg.py:45
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
HcaluLUTTPGCoder::adc2Linear
void adc2Linear(const HBHEDataFrame &df, IntegerCaloSamples &ics) const override
Definition: HcaluLUTTPGCoder.cc:565
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18
HcalTriggerPrimitiveAlgo::passTDC
bool passTDC(const QIE10DataFrame &digi, int ts) const
Definition: HcalTriggerPrimitiveAlgo.cc:680
HcalTriggerPrimitiveAlgo::validChannel
bool validChannel(const QIE10DataFrame &digi, int ts) const
Definition: HcalTriggerPrimitiveAlgo.cc:696
HcalTriggerPrimitiveAlgo::addSignal
void addSignal(const HBHEDataFrame &frame)
adds the signal to the map
Definition: HcalTriggerPrimitiveAlgo.cc:94
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15

◆ addSignal() [4/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE11DataFrame &  frame)
private

Definition at line 249 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), addUpgradeFG(), addUpgradeTDCFG(), cms::cuda::assert(), hcalRecHitTable_cff::detId, HcalBarrel, HcalEndcap, mps_fire::i, incoder_, HcaluLUTTPGCoder::lookupMSB(), theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

249  {
250  HcalDetId detId(frame.id());
251  // prevent QIE11 calibration channels from entering TP emulation
252  if (detId.subdet() != HcalEndcap && detId.subdet() != HcalBarrel)
253  return;
254 
255  std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
256  assert(ids.size() == 1 || ids.size() == 2);
257  IntegerCaloSamples samples1(ids[0], int(frame.samples()));
258 
259  samples1.setPresamples(frame.presamples());
260  incoder_->adc2Linear(frame, samples1);
261 
262  std::vector<std::bitset<2>> msb(frame.samples(), 0);
263  incoder_->lookupMSB(frame, msb);
264 
265  if (ids.size() == 2) {
266  // make a second trigprim for the other one, and share the energy
267  IntegerCaloSamples samples2(ids[1], samples1.size());
268  for (int i = 0; i < samples1.size(); ++i) {
269  samples1[i] = uint32_t(samples1[i]);
270  samples2[i] = samples1[i];
271  }
272  samples2.setPresamples(frame.presamples());
273  addSignal(samples2);
274  addUpgradeFG(ids[1], detId.depth(), msb);
275  addUpgradeTDCFG(ids[1], frame);
276  }
277  addSignal(samples1);
278  addUpgradeFG(ids[0], detId.depth(), msb);
279  addUpgradeTDCFG(ids[0], frame);
280 }
HcalTriggerPrimitiveAlgo::addUpgradeTDCFG
void addUpgradeTDCFG(const HcalTrigTowerDetId &id, const QIE11DataFrame &frame)
Definition: HcalTriggerPrimitiveAlgo.cc:963
mps_fire.i
i
Definition: mps_fire.py:429
HcaluLUTTPGCoder::lookupMSB
void lookupMSB(const HBHEDataFrame &df, std::vector< bool > &msb) const
Definition: HcaluLUTTPGCoder.cc:629
cms::cuda::assert
assert(be >=bs)
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
HcaluLUTTPGCoder::adc2Linear
void adc2Linear(const HBHEDataFrame &df, IntegerCaloSamples &ics) const override
Definition: HcaluLUTTPGCoder.cc:565
HcalTriggerPrimitiveAlgo::addUpgradeFG
void addUpgradeFG(const HcalTrigTowerDetId &id, int depth, const std::vector< std::bitset< 2 >> &bits)
Definition: HcalTriggerPrimitiveAlgo.cc:929
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133
HcalTriggerPrimitiveAlgo::addSignal
void addSignal(const HBHEDataFrame &frame)
adds the signal to the map
Definition: HcalTriggerPrimitiveAlgo.cc:94
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15

◆ addSignal() [5/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const IntegerCaloSamples &  samples)
private

Definition at line 282 of file HcalTriggerPrimitiveAlgo.cc.

References cms::cuda::assert(), fix_saturation_, mps_fire::i, l1ctLayer2EG_cff::id, QIE11_LINEARIZATION_ET, EgammaValidation_cff::samples, theSatMap, and theSumMap.

282  {
283  HcalTrigTowerDetId id(samples.id());
284  SumMap::iterator itr = theSumMap.find(id);
285 
286  if (itr == theSumMap.end()) {
287  theSumMap.insert(std::make_pair(id, samples));
288  } else {
289  // wish CaloSamples had a +=
290  for (int i = 0; i < samples.size(); ++i) {
291  (itr->second)[i] += samples[i];
292  }
293  }
294 
295  // if fix_saturation == true, keep track of tower with saturated input LUT
296  if (fix_saturation_) {
297  SatMap::iterator itr_sat = theSatMap.find(id);
298 
299  assert((itr == theSumMap.end()) == (itr_sat == theSatMap.end()));
300 
301  if (itr_sat == theSatMap.end()) {
302  vector<bool> check_sat;
303  for (int i = 0; i < samples.size(); ++i) {
304  if (!(samples[i] < QIE11_LINEARIZATION_ET)) {
305  check_sat.push_back(true);
306  } else
307  check_sat.push_back(false);
308  }
309  theSatMap.insert(std::make_pair(id, check_sat));
310  } else {
311  for (int i = 0; i < samples.size(); ++i) {
312  if (!(samples[i] < QIE11_LINEARIZATION_ET))
313  (itr_sat->second)[i] = true;
314  }
315  }
316  }
317 }
mps_fire.i
i
Definition: mps_fire.py:429
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:85
cms::cuda::assert
assert(be >=bs)
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18

◆ addUpgradeFG()

void HcalTriggerPrimitiveAlgo::addUpgradeFG ( const HcalTrigTowerDetId &  id,
int  depth,
const std::vector< std::bitset< 2 >> &  bits 
)
private

Definition at line 929 of file HcalTriggerPrimitiveAlgo.cc.

References addFG(), gpuClustering::pixelStatus::bits, hcalRecHitTable_cff::depth, fgUpgradeMap_, mps_fire::i, needLegacyFG(), needUpgradeID(), and validUpgradeFG().

Referenced by addSignal().

931  {
932  if (not validUpgradeFG(id, depth)) {
933  if (needLegacyFG(id)) {
934  std::vector<bool> pseudo(bits.size(), false);
935  addFG(id, pseudo);
936  } else if (needUpgradeID(id, depth)) {
937  // If the tower id is not in the map yet
938  // then for safety's sake add it, otherwise, no need
939  // Likewise, we're here with non-fg depth 7 so the bits are not to be added
940  auto it = fgUpgradeMap_.find(id);
941  if (it == fgUpgradeMap_.end()) {
942  FGUpgradeContainer element;
943  element.resize(bits.size());
944  fgUpgradeMap_.insert(std::make_pair(id, element));
945  }
946  }
947 
948  return;
949  }
950 
951  auto it = fgUpgradeMap_.find(id);
952  if (it == fgUpgradeMap_.end()) {
953  FGUpgradeContainer element;
954  element.resize(bits.size());
955  it = fgUpgradeMap_.insert(std::make_pair(id, element)).first;
956  }
957  for (unsigned int i = 0; i < bits.size(); ++i) {
958  it->second[i][0][depth - 1] = bits[i][0];
959  it->second[i][1][depth - 1] = bits[i][1];
960  }
961 }
HcalTriggerPrimitiveAlgo::validUpgradeFG
bool validUpgradeFG(const HcalTrigTowerDetId &id, int depth) const
Definition: HcalTriggerPrimitiveAlgo.cc:900
mps_fire.i
i
Definition: mps_fire.py:429
gpuClustering::pixelStatus::bits
constexpr uint32_t bits
Definition: gpuClustering.h:25
HcalTriggerPrimitiveAlgo::needLegacyFG
bool needLegacyFG(const HcalTrigTowerDetId &id) const
Definition: HcalTriggerPrimitiveAlgo.cc:910
HcalTriggerPrimitiveAlgo::addFG
void addFG(const HcalTrigTowerDetId &id, std::vector< bool > &msb)
Definition: HcalTriggerPrimitiveAlgo.cc:890
HcalTriggerPrimitiveAlgo::needUpgradeID
bool needUpgradeID(const HcalTrigTowerDetId &id, int depth) const
Definition: HcalTriggerPrimitiveAlgo.cc:919
HcalTriggerPrimitiveAlgo::FGUpgradeContainer
std::vector< HcalFinegrainBit::Tower > FGUpgradeContainer
Definition: HcalTriggerPrimitiveAlgo.h:213
hcalRecHitTable_cff.depth
depth
Definition: hcalRecHitTable_cff.py:17

◆ addUpgradeTDCFG()

void HcalTriggerPrimitiveAlgo::addUpgradeTDCFG ( const HcalTrigTowerDetId &  id,
const QIE11DataFrame &  frame 
)
private

Definition at line 963 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), cms::cuda::assert(), hcalRecHitTable_cff::detId, fgUpgradeTDCMap_, HcaluLUTTPGCoder::group0FGbits(), HcalBarrel, HcalEndcap, mps_fire::i, incoder_, theTrigTowerGeometry, and HcalTrigTowerGeometry::towerIds().

Referenced by addSignal().

963  {
964  HcalDetId detId(frame.id());
965  if (detId.subdet() != HcalEndcap && detId.subdet() != HcalBarrel)
966  return;
967 
968  std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
969  assert(ids.size() == 1 || ids.size() == 2);
970  IntegerCaloSamples samples1(ids[0], int(frame.samples()));
971  samples1.setPresamples(frame.presamples());
972  incoder_->adc2Linear(frame, samples1); // use linearization LUT
973  std::vector<unsigned short> bits12_15 = incoder_->group0FGbits(frame); // get 4 energy bits (12-15) from group 0 LUT
974 
975  bool is_compressed = false;
976  if (detId.subdet() == HcalBarrel) {
977  is_compressed = (frame.flavor() == 3);
978  // 0 if frame.flavor is 0 (uncompressed), 1 if frame.flavor is 3 (compressed)
979  }
980 
981  auto it = fgUpgradeTDCMap_.find(id);
982  if (it == fgUpgradeTDCMap_.end()) {
983  FGUpgradeTDCContainer element;
984  element.resize(frame.samples());
985  it = fgUpgradeTDCMap_.insert(std::make_pair(id, element)).first;
986  }
987  for (int i = 0; i < frame.samples(); i++) {
988  it->second[i][detId.depth() - 1] =
989  std::make_pair(std::make_pair(bits12_15[i], is_compressed), std::make_pair(frame[i].tdc(), samples1[i]));
990  }
991 }
mps_fire.i
i
Definition: mps_fire.py:429
HcaluLUTTPGCoder::group0FGbits
std::vector< unsigned short > group0FGbits(const QIE11DataFrame &df) const
Definition: HcaluLUTTPGCoder.cc:602
cms::cuda::assert
assert(be >=bs)
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
HcalTriggerPrimitiveAlgo::FGUpgradeTDCContainer
std::vector< HcalFinegrainBit::TowerTDC > FGUpgradeTDCContainer
Definition: HcalTriggerPrimitiveAlgo.h:217
HcaluLUTTPGCoder::adc2Linear
void adc2Linear(const HBHEDataFrame &df, IntegerCaloSamples &ics) const override
Definition: HcaluLUTTPGCoder.cc:565
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15

◆ analyze()

void HcalTriggerPrimitiveAlgo::analyze ( IntegerCaloSamples &  samples,
HcalTriggerPrimitiveDigi &  result 
)
private

adds the actual digis

Definition at line 319 of file HcalTriggerPrimitiveAlgo.cc.

References HcalTPGCompressor::compress(), fgMap_, mps_fire::i, heavyIonCSV_trainingSettings::idx, createfilelist::int, numberOfPresamples_, numberOfSamples_, outcoder_, convertSQLitetoXML_cfg::output, peak_finder_algorithm_, peakfind_, QIE8_LINEARIZATION_ET, mps_fire::result, EgammaValidation_cff::samples, edm::shift, theThreshold, and weights_.

Referenced by run().

319  {
320  int shrink = weights_.size() - 1;
321  std::vector<bool>& msb = fgMap_[samples.id()];
322  IntegerCaloSamples sum(samples.id(), samples.size());
323 
324  //slide algo window
325  for (int ibin = 0; ibin < int(samples.size()) - shrink; ++ibin) {
326  int algosumvalue = 0;
327  for (unsigned int i = 0; i < weights_.size(); i++) {
328  //add up value * scale factor
329  algosumvalue += int(samples[ibin + i] * weights_[i]);
330  }
331  if (algosumvalue < 0)
332  sum[ibin] = 0; // low-side
333  //high-side
334  //else if (algosumvalue>QIE8_LINEARIZATION_ET) sum[ibin]=QIE8_LINEARIZATION_ET;
335  else
336  sum[ibin] = algosumvalue; //assign value to sum[]
337  }
338 
339  // Align digis and TP
340  int dgPresamples = samples.presamples();
341  int tpPresamples = numberOfPresamples_;
342  int shift = dgPresamples - tpPresamples;
343  int dgSamples = samples.size();
344  int tpSamples = numberOfSamples_;
345  if (peakfind_) {
346  if ((shift < shrink) || (shift + tpSamples + shrink > dgSamples - (peak_finder_algorithm_ - 1))) {
347  edm::LogInfo("HcalTriggerPrimitiveAlgo::analyze")
348  << "TP presample or size from the configuration file is out of the accessible range. Using digi values from "
349  "data instead...";
350  shift = shrink;
351  tpPresamples = dgPresamples - shrink;
352  tpSamples = dgSamples - (peak_finder_algorithm_ - 1) - shrink - shift;
353  }
354  }
355 
356  std::vector<int> finegrain(tpSamples, false);
357 
358  IntegerCaloSamples output(samples.id(), tpSamples);
359  output.setPresamples(tpPresamples);
360 
361  for (int ibin = 0; ibin < tpSamples; ++ibin) {
362  // ibin - index for output TP
363  // idx - index for samples + shift
364  int idx = ibin + shift;
365 
366  //Peak finding
367  if (peakfind_) {
368  bool isPeak = false;
369  switch (peak_finder_algorithm_) {
370  case 1:
371  isPeak = (samples[idx] > samples[idx - 1] && samples[idx] >= samples[idx + 1] && samples[idx] > theThreshold);
372  break;
373  case 2:
374  isPeak = (sum[idx] > sum[idx - 1] && sum[idx] >= sum[idx + 1] && sum[idx] > theThreshold);
375  break;
376  default:
377  break;
378  }
379 
380  if (isPeak) {
381  output[ibin] = std::min<unsigned int>(sum[idx], QIE8_LINEARIZATION_ET);
382  finegrain[ibin] = msb[idx];
383  }
384  // Not a peak
385  else
386  output[ibin] = 0;
387  } else { // No peak finding, just output running sum
388  output[ibin] = std::min<unsigned int>(sum[idx], QIE8_LINEARIZATION_ET);
389  finegrain[ibin] = msb[idx];
390  }
391 
392  // Only Pegged for 1-TS algo.
393  if (peak_finder_algorithm_ == 1) {
394  if (samples[idx] >= QIE8_LINEARIZATION_ET)
395  output[ibin] = QIE8_LINEARIZATION_ET;
396  }
397  }
398  outcoder_->compress(output, finegrain, result);
399 }
mps_fire.i
i
Definition: mps_fire.py:429
mps_fire.result
result
Definition: mps_fire.py:311
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:72
HcalTPGCompressor::compress
void compress(const IntegerCaloSamples &ics, const std::vector< int > &fineGrain, HcalTriggerPrimitiveDigi &digi) const
Definition: HcalTPGCompressor.cc:6
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:7

◆ analyzeHF()

void HcalTriggerPrimitiveAlgo::analyzeHF ( IntegerCaloSamples &  samples,
HcalTriggerPrimitiveDigi &  result,
const int  hf_lumi_shift 
)
private

Definition at line 565 of file HcalTriggerPrimitiveAlgo.cc.

References cms::cuda::assert(), HcalTPGCompressor::compress(), hcalRecHitTable_cff::detId, FG_threshold_, HF_Veto, heavyIonCSV_trainingSettings::idx, createfilelist::int, numberOfPresamplesHF_, numberOfSamplesHF_, outcoder_, convertSQLitetoXML_cfg::output, PMT_NoiseThreshold_, QIE8_LINEARIZATION_ET, mps_fire::result, EgammaValidation_cff::samples, edm::shift, theTowerMapFGSum, and PbPb_ZMuSkimMuonDPG_cff::veto.

Referenced by run().

567  {
568  HcalTrigTowerDetId detId(samples.id());
569 
570  // Align digis and TP
571  int dgPresamples = samples.presamples();
572  int tpPresamples = numberOfPresamplesHF_;
573  int shift = dgPresamples - tpPresamples;
574  int dgSamples = samples.size();
575  int tpSamples = numberOfSamplesHF_;
576  if (shift < 0 || shift + tpSamples > dgSamples) {
577  edm::LogInfo("HcalTriggerPrimitiveAlgo::analyzeHF")
578  << "TP presample or size from the configuration file is out of the accessible range. Using digi values from "
579  "data instead...";
580  tpPresamples = dgPresamples;
581  shift = 0;
582  tpSamples = dgSamples;
583  }
584 
585  std::vector<int> finegrain(tpSamples, false);
586 
587  TowerMapFGSum::const_iterator tower2fg = theTowerMapFGSum.find(detId);
588  assert(tower2fg != theTowerMapFGSum.end());
589 
590  const SumFGContainer& sumFG = tower2fg->second;
591  // Loop over all L+S pairs that mapped from samples.id()
592  // Note: 1 samples.id() = 6 x (L+S) without noZS
593  for (SumFGContainer::const_iterator sumFGItr = sumFG.begin(); sumFGItr != sumFG.end(); ++sumFGItr) {
594  const std::vector<bool>& veto = HF_Veto[sumFGItr->id().rawId()];
595  for (int ibin = 0; ibin < tpSamples; ++ibin) {
596  int idx = ibin + shift;
597  // if not vetod, add L+S to total sum and calculate FG
598  bool vetoed = idx < int(veto.size()) && veto[idx];
599  if (!(vetoed && (*sumFGItr)[idx] > PMT_NoiseThreshold_)) {
600  samples[idx] += (*sumFGItr)[idx];
601  finegrain[ibin] = (finegrain[ibin] || (*sumFGItr)[idx] >= FG_threshold_);
602  }
603  }
604  }
605 
606  IntegerCaloSamples output(samples.id(), tpSamples);
607  output.setPresamples(tpPresamples);
608 
609  for (int ibin = 0; ibin < tpSamples; ++ibin) {
610  int idx = ibin + shift;
611  output[ibin] = samples[idx] >> hf_lumi_shift;
612  static const int MAX_OUTPUT = QIE8_LINEARIZATION_ET; // QIE8_LINEARIZATION_ET = 1023
613  if (output[ibin] > MAX_OUTPUT)
614  output[ibin] = MAX_OUTPUT;
615  }
616  outcoder_->compress(output, finegrain, result);
617 }
mps_fire.result
result
Definition: mps_fire.py:311
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
PbPb_ZMuSkimMuonDPG_cff.veto
veto
Definition: PbPb_ZMuSkimMuonDPG_cff.py:61
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:72
cms::cuda::assert
assert(be >=bs)
HcalTriggerPrimitiveAlgo::SumFGContainer
std::vector< IntegerCaloSamples > SumFGContainer
Definition: HcalTriggerPrimitiveAlgo.h:193
HcalTPGCompressor::compress
void compress(const IntegerCaloSamples &ics, const std::vector< int > &fineGrain, HcalTriggerPrimitiveDigi &digi) const
Definition: HcalTPGCompressor.cc:6
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:7

◆ analyzeHF2016()

void HcalTriggerPrimitiveAlgo::analyzeHF2016 ( const IntegerCaloSamples &  SAMPLES,
HcalTriggerPrimitiveDigi &  result,
const int  HF_LUMI_SHIFT,
const HcalFeatureBit *  HCALFEM 
)
private

Definition at line 619 of file HcalTriggerPrimitiveAlgo.cc.

References cms::cuda::assert(), newFWLiteAna::bin, HcalTPGCompressor::compress(), hcalRecHitTable_cff::detId, FG_HF_thresholds_, HcalFeatureBit::fineGrainbit(), FIRST_FINEGRAIN_TOWER, B2GTnPMonitor_cfi::item, SiStripPI::min, numberOfPresamplesHF_, numberOfSamplesHF_, outcoder_, convertSQLitetoXML_cfg::output, QIE8_LINEARIZATION_ET, mps_fire::result, EgammaValidation_cff::samples, and theHFDetailMap.

Referenced by run().

622  {
623  // Align digis and TP
624  const int SHIFT = samples.presamples() - numberOfPresamplesHF_;
625  assert(SHIFT >= 0);
626  assert((SHIFT + numberOfSamplesHF_) <= samples.size());
627 
628  // Try to find the HFDetails from the map corresponding to our samples
629  const HcalTrigTowerDetId detId(samples.id());
630  HFDetailMap::const_iterator it = theHFDetailMap.find(detId);
631  // Missing values will give an empty digi
632  if (it == theHFDetailMap.end()) {
633  return;
634  }
635 
636  std::vector<std::bitset<2>> finegrain(numberOfSamplesHF_, false);
637 
638  // Set up out output of IntergerCaloSamples
639  IntegerCaloSamples output(samples.id(), numberOfSamplesHF_);
640  output.setPresamples(numberOfPresamplesHF_);
641 
642  for (const auto& item : it->second) {
643  auto& details = item.second;
644  for (int ibin = 0; ibin < numberOfSamplesHF_; ++ibin) {
645  const int IDX = ibin + SHIFT;
646  int long_fiber_val = 0;
647  if (IDX < details.long_fiber.size()) {
648  long_fiber_val = details.long_fiber[IDX];
649  }
650  int short_fiber_val = 0;
651  if (IDX < details.short_fiber.size()) {
652  short_fiber_val = details.short_fiber[IDX];
653  }
654  output[ibin] += (long_fiber_val + short_fiber_val);
655 
656  uint32_t ADCLong = details.LongDigi[ibin].adc();
657  uint32_t ADCShort = details.ShortDigi[ibin].adc();
658 
659  if (details.LongDigi.id().ietaAbs() >= FIRST_FINEGRAIN_TOWER) {
660  finegrain[ibin][1] = (ADCLong > FG_HF_thresholds_[0] || ADCShort > FG_HF_thresholds_[0]);
661 
662  if (embit != nullptr)
663  finegrain[ibin][0] = embit->fineGrainbit(details.ShortDigi, details.LongDigi, ibin);
664  }
665  }
666  }
667 
668  for (int bin = 0; bin < numberOfSamplesHF_; ++bin) {
669  static const unsigned int MAX_OUTPUT = QIE8_LINEARIZATION_ET; // QIE8_LINEARIZATION_ET = 1023
670  output[bin] = min({MAX_OUTPUT, output[bin] >> hf_lumi_shift});
671  }
672 
673  std::vector<int> finegrain_converted;
674  finegrain_converted.reserve(finegrain.size());
675  for (const auto& fg : finegrain)
676  finegrain_converted.push_back(fg.to_ulong());
677  outcoder_->compress(output, finegrain_converted, result);
678 }
HcalTriggerPrimitiveAlgo::FG_HF_thresholds_
std::vector< uint32_t > FG_HF_thresholds_
Definition: HcalTriggerPrimitiveAlgo.h:142
mps_fire.result
result
Definition: mps_fire.py:311
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:72
cms::cuda::assert
assert(be >=bs)
B2GTnPMonitor_cfi.item
item
Definition: B2GTnPMonitor_cfi.py:148
HcalTPGCompressor::compress
void compress(const IntegerCaloSamples &ics, const std::vector< int > &fineGrain, HcalTriggerPrimitiveDigi &digi) const
Definition: HcalTPGCompressor.cc:6
details
Definition: helper.h:56
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18

◆ analyzeHFQIE10()

void HcalTriggerPrimitiveAlgo::analyzeHFQIE10 ( const IntegerCaloSamples &  SAMPLES,
HcalTriggerPrimitiveDigi &  result,
const int  HF_LUMI_SHIFT,
const HcalFeatureBit *  HCALFEM 
)
private

Definition at line 708 of file HcalTriggerPrimitiveAlgo.cc.

References cms::cuda::assert(), newFWLiteAna::bin, HcalTPGCompressor::compress(), hcalRecHitTable_cff::detId, FG_HF_thresholds_, HcalFeatureBit::fineGrainbit(), FIRST_FINEGRAIN_TOWER, mps_fire::i, heavyIonCSV_trainingSettings::idx, createfilelist::int, B2GTnPMonitor_cfi::item, SiStripPI::min, numberOfPresamplesHF_, numberOfSamplesHF_, or, outcoder_, convertSQLitetoXML_cfg::output, QIE10_LINEARIZATION_ET, QIE10_MAX_LINEARIZATION_ET, mps_fire::result, EgammaValidation_cff::samples, edm::shift, theHFUpgradeDetailMap, and useTDCInMinBiasBits_.

Referenced by run().

711  {
712  // Align digis and TP
713  const int shift = samples.presamples() - numberOfPresamplesHF_;
714  assert(shift >= 0);
715  assert((shift + numberOfSamplesHF_) <= samples.size());
716  assert(hf_lumi_shift >= 2);
717 
718  // Try to find the HFDetails from the map corresponding to our samples
719  const HcalTrigTowerDetId detId(samples.id());
720  auto it = theHFUpgradeDetailMap.find(detId);
721  // Missing values will give an empty digi
722  if (it == theHFUpgradeDetailMap.end()) {
723  return;
724  }
725 
726  std::vector<std::bitset<2>> finegrain(numberOfSamplesHF_, false);
727 
728  // Set up out output of IntergerCaloSamples
729  IntegerCaloSamples output(samples.id(), numberOfSamplesHF_);
730  output.setPresamples(numberOfPresamplesHF_);
731 
732  for (const auto& item : it->second) {
733  auto& details = item.second;
734  for (int ibin = 0; ibin < numberOfSamplesHF_; ++ibin) {
735  const int idx = ibin + shift;
736 
737  int long_fiber_val = 0;
738  int long_fiber_count = 0;
739  int short_fiber_val = 0;
740  int short_fiber_count = 0;
741 
742  bool saturated = false;
743 
744  for (auto i : {0, 2}) {
745  if (idx < details[i].samples.size() and details[i].validity[idx] and details[i].passTDC[idx]) {
746  long_fiber_val += details[i].samples[idx];
747  saturated = saturated || (details[i].samples[idx] == QIE10_LINEARIZATION_ET);
748  ++long_fiber_count;
749  }
750  }
751  for (auto i : {1, 3}) {
752  if (idx < details[i].samples.size() and details[i].validity[idx] and details[i].passTDC[idx]) {
753  short_fiber_val += details[i].samples[idx];
754  saturated = saturated || (details[i].samples[idx] == QIE10_LINEARIZATION_ET);
755  ++short_fiber_count;
756  }
757  }
758 
759  if (saturated) {
760  output[ibin] = QIE10_MAX_LINEARIZATION_ET;
761  } else {
762  // For details of the energy handling, see:
763  // https://cms-docdb.cern.ch/cgi-bin/DocDB/ShowDocument?docid=12306
764  // If both readouts are valid, average of the two energies is taken
765  // division by 2 is compensated by adjusting the total scale shift in the end
766  if (long_fiber_count == 2)
767  long_fiber_val >>= 1;
768  if (short_fiber_count == 2)
769  short_fiber_val >>= 1;
770 
771  auto sum = long_fiber_val + short_fiber_val;
772  // Similar to above, if both channels are valid,
773  // average of the two energies is calculated
774  // division by 2 here is also compensated by adjusting the total scale shift in the end
775  if (long_fiber_count > 0 and short_fiber_count > 0)
776  sum >>= 1;
777 
778  output[ibin] += sum;
779  }
780 
781  for (const auto& detail : details) {
782  if (idx < int(detail.digi.size()) and detail.validity[idx] and
783  HcalDetId(detail.digi.id()).ietaAbs() >= FIRST_FINEGRAIN_TOWER) {
784  if (useTDCInMinBiasBits_ && !detail.passTDC[idx])
785  continue;
786  finegrain[ibin][1] = finegrain[ibin][1] or detail.fgbits[idx][0];
787  // what is commonly called the "second" HF min-bias bit is
788  // actually the 0-th bit, which can also be used instead for the EM bit
789  // (called finegrain[ibin][0] below) in non-HI running
790  finegrain[ibin][0] = finegrain[ibin][0] or detail.fgbits[idx][1];
791  }
792  }
793  // the EM bit is only used if the "second" FG bit is disabled
794  if (embit != nullptr and FG_HF_thresholds_.at(1) != 255) {
795  finegrain[ibin][0] = embit->fineGrainbit(details[1].digi,
796  details[3].digi,
797  details[0].digi,
798  details[2].digi,
799  details[1].validity[idx],
800  details[3].validity[idx],
801  details[0].validity[idx],
802  details[2].validity[idx],
803  idx);
804  }
805  }
806  }
807 
808  for (int bin = 0; bin < numberOfSamplesHF_; ++bin) {
809  output[bin] = min({(unsigned int)QIE10_MAX_LINEARIZATION_ET, output[bin] >> (hf_lumi_shift - 2)});
810  }
811  std::vector<int> finegrain_converted;
812  finegrain_converted.reserve(finegrain.size());
813  for (const auto& fg : finegrain)
814  finegrain_converted.push_back(fg.to_ulong());
815  outcoder_->compress(output, finegrain_converted, result);
816 }
HcalTriggerPrimitiveAlgo::FG_HF_thresholds_
std::vector< uint32_t > FG_HF_thresholds_
Definition: HcalTriggerPrimitiveAlgo.h:142
HcalTriggerPrimitiveAlgo::theHFUpgradeDetailMap
HFUpgradeDetailMap theHFUpgradeDetailMap
Definition: HcalTriggerPrimitiveAlgo.h:191
mps_fire.i
i
Definition: mps_fire.py:429
mps_fire.result
result
Definition: mps_fire.py:311
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
HcalTriggerPrimitiveAlgo::QIE10_MAX_LINEARIZATION_ET
static const int QIE10_MAX_LINEARIZATION_ET
Definition: HcalTriggerPrimitiveAlgo.h:247
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:72
cms::cuda::assert
assert(be >=bs)
B2GTnPMonitor_cfi.item
item
Definition: B2GTnPMonitor_cfi.py:148
HcalTPGCompressor::compress
void compress(const IntegerCaloSamples &ics, const std::vector< int > &fineGrain, HcalTriggerPrimitiveDigi &digi) const
Definition: HcalTPGCompressor.cc:6
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
details
Definition: helper.h:56
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:7

◆ analyzeQIE11()

void HcalTriggerPrimitiveAlgo::analyzeQIE11 ( IntegerCaloSamples &  samples,
std::vector< bool >  sample_saturation,
HcalTriggerPrimitiveDigi &  result,
const HcalFinegrainBit &  fg_algo 
)
private

Definition at line 401 of file HcalTriggerPrimitiveAlgo.cc.

References HcalTPGCompressor::compress(), HcalFinegrainBit::compute(), hcalRecHitTable_cff::detId, fgUpgradeMap_, fgUpgradeTDCMap_, fix_saturation_, mps_fire::i, heavyIonCSV_trainingSettings::idx, createfilelist::int, HcalTopology::lastHBRing(), numberOfFilterPresamplesHBQIE11_, numberOfFilterPresamplesHEQIE11_, numberOfPresamples_, numberOfSamples_, outcoder_, convertSQLitetoXML_cfg::output, QIE11_MAX_LINEARIZATION_ET, mps_fire::result, ecalGpuTask_cfi::sample, EgammaValidation_cff::samples, edm::shift, theThreshold, theTrigTowerGeometry, HcalTrigTowerGeometry::topology(), HcalTrigTowerGeometry::towerIds(), and weightsQIE11_.

Referenced by run().

404  {
405  HcalDetId detId(samples.id());
406 
407  // Get the |ieta| for current sample
408  int theIeta = detId.ietaAbs();
409 
410  unsigned int dgSamples = samples.size();
411  unsigned int dgPresamples = samples.presamples();
412 
413  unsigned int tpSamples = numberOfSamples_;
414  unsigned int tpPresamples = numberOfPresamples_;
415 
416  unsigned int filterSamples = weightsQIE11_[theIeta].size();
417  unsigned int filterPresamples = theIeta > theTrigTowerGeometry->topology().lastHBRing()
418  ? numberOfFilterPresamplesHEQIE11_
419  : numberOfFilterPresamplesHBQIE11_;
420 
421  unsigned int shift = dgPresamples - tpPresamples;
422 
423  // shrink keeps the FIR filter from going off the end of the 8TS vector
424  unsigned int shrink = filterSamples - 1;
425 
426  auto& msb = fgUpgradeMap_[samples.id()];
427  auto& timingTDC = fgUpgradeTDCMap_[samples.id()];
428  IntegerCaloSamples sum(samples.id(), samples.size());
429 
430  std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
431 
432  // keep track of tower with saturated energy and force the total TP saturated
433  bool force_saturation[samples.size()];
434  for (int i = 0; i < samples.size(); i++) {
435  force_saturation[i] = false;
436  }
437 
438  //slide algo window
439  for (unsigned int ibin = 0; ibin < dgSamples - shrink; ++ibin) {
440  int algosumvalue = 0;
441  bool check_sat = false;
442  //TP energy calculation for PFA2
443  if (weightsQIE11_[theIeta][0] == 255) {
444  for (unsigned int i = 0; i < filterSamples; i++) {
445  //add up value * scale factor
446  // In addition, divide by two in the 10 degree phi segmentation region
447  // to mimic 5 degree segmentation for the trigger
448  unsigned int sample = samples[ibin + i];
449 
450  if (fix_saturation_ && (sample_saturation.size() > ibin + i))
451  check_sat = (check_sat | sample_saturation[ibin + i] | (sample > QIE11_MAX_LINEARIZATION_ET));
452 
453  if (sample > QIE11_MAX_LINEARIZATION_ET)
454  sample = QIE11_MAX_LINEARIZATION_ET;
455 
456  // Usually use a segmentation factor of 1.0 but for ieta >= 21 use 2
457  int segmentationFactor = 1;
458  if (ids.size() == 2) {
459  segmentationFactor = 2;
460  }
461 
462  algosumvalue += int(sample / segmentationFactor);
463  }
464  if (algosumvalue < 0)
465  sum[ibin] = 0; // low-side
466  //high-side
467  //else if (algosumvalue>QIE11_LINEARIZATION_ET) sum[ibin]=QIE11_LINEARIZATION_ET;
468  else
469  sum[ibin] = algosumvalue; //assign value to sum[]
470 
471  if (check_sat)
472  force_saturation[ibin] = true;
473  //TP energy calculation for PFA1' and PFA1
474  } else {
475  //add up value * scale factor
476  // In addition, divide by two in the 10 degree phi segmentation region
477  // to mimic 5 degree segmentation for the trigger
478  int sampleTS = samples[ibin + 1];
479  int sampleTSminus1 = samples[ibin];
480 
481  if (fix_saturation_ && (sample_saturation.size() > ibin + 1))
482  check_sat = (sample_saturation[ibin + 1] || (sampleTS >= QIE11_MAX_LINEARIZATION_ET) ||
483  sample_saturation[ibin] || (sampleTSminus1 >= QIE11_MAX_LINEARIZATION_ET));
484 
485  if (sampleTS > QIE11_MAX_LINEARIZATION_ET)
486  sampleTS = QIE11_MAX_LINEARIZATION_ET;
487 
488  if (sampleTSminus1 > QIE11_MAX_LINEARIZATION_ET)
489  sampleTSminus1 = QIE11_MAX_LINEARIZATION_ET;
490 
491  // Usually use a segmentation factor of 1.0 but for ieta >= 21 use 2
492  int segmentationFactor = 1;
493  if (ids.size() == 2) {
494  segmentationFactor = 2;
495  }
496 
497  // Based on the |ieta| of the sample, retrieve the correct region weight
498  int theWeight = weightsQIE11_[theIeta][0];
499 
500  algosumvalue = ((sampleTS << 8) - (sampleTSminus1 * theWeight)) / 256 / segmentationFactor;
501 
502  if (algosumvalue < 0)
503  sum[ibin] = 0; // low-side
504  //high-side
505  //else if (algosumvalue>QIE11_LINEARIZATION_ET) sum[ibin]=QIE11_LINEARIZATION_ET;
506  else
507  sum[ibin] = algosumvalue; //assign value to sum[]
508 
509  if (check_sat)
510  force_saturation[ibin] = true;
511  }
512  }
513 
514  std::vector<int> finegrain(tpSamples, false);
515 
516  IntegerCaloSamples output(samples.id(), tpSamples);
517  output.setPresamples(tpPresamples);
518 
519  for (unsigned int ibin = 0; ibin < tpSamples; ++ibin) {
520  // ibin - index for output TP
521  // idx - index for samples + shift - filterPresamples
522  int idx = ibin + shift - filterPresamples;
523 
524  // When idx is <= 0 peakfind would compare out-of-bounds of the vector. Avoid this ambiguity
525  if (idx <= 0) {
526  output[ibin] = 0;
527  continue;
528  }
529 
530  //Only run the peak-finder when the PFA2 FIR filter is running, which corresponds to weights = 1
531  if (weightsQIE11_[theIeta][0] == 255) {
532  bool isPeak = (sum[idx] > sum[idx - 1] && sum[idx] >= sum[idx + 1] && sum[idx] > theThreshold);
533  if (isPeak) {
534  output[ibin] = std::min<unsigned int>(sum[idx], QIE11_MAX_LINEARIZATION_ET);
535 
536  if (fix_saturation_ && force_saturation[idx] && ids.size() == 2)
537  output[ibin] = QIE11_MAX_LINEARIZATION_ET / 2;
538  else if (fix_saturation_ && force_saturation[idx])
539  output[ibin] = QIE11_MAX_LINEARIZATION_ET;
540 
541  } else {
542  // Not a peak
543  output[ibin] = 0;
544  }
545  } else {
546  output[ibin] = std::min<unsigned int>(sum[idx], QIE11_MAX_LINEARIZATION_ET);
547 
548  if (fix_saturation_ && force_saturation[idx] && ids.size() == 2)
549  output[ibin] = QIE11_MAX_LINEARIZATION_ET / 2;
550  else if (fix_saturation_ && force_saturation[idx])
551  output[ibin] = QIE11_MAX_LINEARIZATION_ET;
552  }
553  // peak-finding is not applied for FG bits
554  // compute(msb) returns two bits (MIP). compute(timingTDC,ids) returns 6 bits (1 depth, 1 prompt, 1 delayed 01, 1 delayed 10, 2 reserved)
555  finegrain[ibin] = fg_algo.compute(timingTDC[idx + filterPresamples], ids[0]).to_ulong() |
556  fg_algo.compute(msb[idx + filterPresamples]).to_ulong() << 4;
557  if (ibin == tpPresamples && (idx + filterPresamples) != dgPresamples)
558  edm::LogError("HcalTriggerPritimveAlgo")
559  << "TP SOI (tpPresamples = " << tpPresamples
560  << ") is not aligned with digi SOI (dgPresamples = " << dgPresamples << ")";
561  }
562  outcoder_->compress(output, finegrain, result);
563 }
mps_fire.i
i
Definition: mps_fire.py:429
mps_fire.result
result
Definition: mps_fire.py:311
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:72
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
HcalTriggerPrimitiveAlgo::weightsQIE11_
std::array< std::array< int, 2 >, 29 > weightsQIE11_
Definition: HcalTriggerPrimitiveAlgo.h:139
ecalGpuTask_cfi.sample
sample
Definition: ecalGpuTask_cfi.py:47
HcalTrigTowerGeometry::topology
const HcalTopology & topology() const
Definition: HcalTrigTowerGeometry.h:27
HcalTPGCompressor::compress
void compress(const IntegerCaloSamples &ics, const std::vector< int > &fineGrain, HcalTriggerPrimitiveDigi &digi) const
Definition: HcalTPGCompressor.cc:6
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
HcalTopology::lastHBRing
int lastHBRing() const
Definition: HcalTopology.h:92
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
HcalFinegrainBit::compute
std::bitset< 2 > compute(const Tower &) const
Definition: HcalFinegrainBit.cc:5
HcalTriggerPrimitiveAlgo::QIE11_MAX_LINEARIZATION_ET
static const int QIE11_MAX_LINEARIZATION_ET
Definition: HcalTriggerPrimitiveAlgo.h:248
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:18
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:7
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15

◆ needLegacyFG()

bool HcalTriggerPrimitiveAlgo::needLegacyFG ( const HcalTrigTowerDetId &  id) const
private

Definition at line 910 of file HcalTriggerPrimitiveAlgo.cc.

References HBHE_OVERLAP_TOWER, and upgrade_hb_.

Referenced by addUpgradeFG().

910  {
911  // This tower (ietaAbs == 16) does not accept upgraded FG bits,
912  // but needs pseudo legacy ones to ensure that the tower is processed
913  // even when the QIE8 depths in front of it do not have energy deposits.
914  if (id.ietaAbs() == HBHE_OVERLAP_TOWER and not upgrade_hb_)
915  return true;
916  return false;
917 }

◆ needUpgradeID()

bool HcalTriggerPrimitiveAlgo::needUpgradeID ( const HcalTrigTowerDetId &  id,
int  depth 
) const
private

Definition at line 919 of file HcalTriggerPrimitiveAlgo.cc.

References FIRST_DEPTH7_TOWER, and LAST_FINEGRAIN_DEPTH.

Referenced by addUpgradeFG().

919  {
920  // Depth 7 for TT 26, 27, and 28 is not considered a fine grain depth.
921  // However, the trigger tower for these ieta should still be added to the fgUpgradeMap_
922  // Otherwise, depth 7-only signal will not be analyzed.
923  unsigned int aieta = id.ietaAbs();
924  if (aieta >= FIRST_DEPTH7_TOWER and aieta <= LAST_FINEGRAIN_TOWER and depth > LAST_FINEGRAIN_DEPTH)
925  return true;
926  return false;
927 }

◆ overrideParameters()

void HcalTriggerPrimitiveAlgo::overrideParameters ( const edm::ParameterSet &  ps)

Definition at line 81 of file HcalTriggerPrimitiveAlgo.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), override_adc_hf_, override_adc_hf_value_, override_parameters_, override_tdc_hf_, and override_tdc_hf_value_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

81  {
82  override_parameters_ = ps;
83 
84  if (override_parameters_.exists("ADCThresholdHF")) {
85  override_adc_hf_ = true;
86  override_adc_hf_value_ = override_parameters_.getParameter<uint32_t>("ADCThresholdHF");
87  }
88  if (override_parameters_.exists("TDCMaskHF")) {
89  override_tdc_hf_ = true;
90  override_tdc_hf_value_ = override_parameters_.getParameter<unsigned long long>("TDCMaskHF");
91  }
92 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:681
HcalTriggerPrimitiveAlgo::override_tdc_hf_value_
unsigned long long override_tdc_hf_value_
Definition: HcalTriggerPrimitiveAlgo.h:232

◆ passTDC()

bool HcalTriggerPrimitiveAlgo::passTDC ( const QIE10DataFrame &  digi,
int  ts 
) const
private

Definition at line 680 of file HcalTriggerPrimitiveAlgo.cc.

References gpuClustering::adc, conditions_, HcalDbService::getHcalTPParameters(), override_adc_hf_, override_adc_hf_value_, override_tdc_hf_, and override_tdc_hf_value_.

Referenced by addSignal().

680  {
681  auto parameters = conditions_->getHcalTPParameters();
682  auto adc_threshold = parameters->getADCThresholdHF();
683  auto tdc_mask = parameters->getTDCMaskHF();
684 
685  if (override_adc_hf_)
686  adc_threshold = override_adc_hf_value_;
687  if (override_tdc_hf_)
688  tdc_mask = override_tdc_hf_value_;
689 
690  if (digi[ts].adc() < adc_threshold)
691  return true;
692 
693  return (1ul << digi[ts].le_tdc()) & tdc_mask;
694 }
HcalTriggerPrimitiveAlgo::conditions_
const HcalDbService * conditions_
Definition: HcalTriggerPrimitiveAlgo.h:135
HcalDbService::getHcalTPParameters
const HcalTPParameters * getHcalTPParameters() const
Definition: HcalDbService.cc:413
HcalTriggerPrimitiveAlgo::override_tdc_hf_value_
unsigned long long override_tdc_hf_value_
Definition: HcalTriggerPrimitiveAlgo.h:232
gpuClustering::adc
uint16_t *__restrict__ uint16_t const *__restrict__ adc
Definition: gpuClusterChargeCut.h:21

◆ run()

template<typename... Digis>
void HcalTriggerPrimitiveAlgo::run ( const HcalTPGCoder *  incoder,
const HcalTPGCompressor *  outcoder,
const HcalDbService *  conditions,
HcalTrigPrimDigiCollection &  result,
const HcalTrigTowerGeometry *  trigTowerGeometry,
float  rctlsb,
const HcalFeatureBit *  LongvrsShortCut,
const Digis &...  digis 
)

Definition at line 252 of file HcalTriggerPrimitiveAlgo.h.

References addDigis(), analyze(), analyzeHF(), analyzeHF2016(), analyzeHFQIE10(), analyzeQIE11(), submitPVValidationJobs::conditions, conditions_, hcalRecHitTable_cff::detId, edm::ParameterSet::exists(), fgMap_, fgUpgradeMap_, fgUpgradeTDCMap_, HcalTrigTowerGeometry::firstHFTower(), HcalTPParameters::getFGVersionHBHE(), HcalDbService::getHcalTPParameters(), edm::ParameterSet::getParameter(), HF_Veto, incoder_, B2GTnPMonitor_cfi::item, LongvrsShortCut, NCTScaleShift, or, outcoder_, override_parameters_, RCTScaleShift, mps_fire::result, theHFDetailMap, theHFUpgradeDetailMap, theSatMap, theSumMap, theTowerMapFGSum, theTrigTowerGeometry, upgrade_hb_, upgrade_he_, upgrade_hf_, and BeamSplash_cfg::version.

Referenced by HcalTrigPrimDigiProducer::produce().

259  {
260  theTrigTowerGeometry = trigTowerGeometry;
261 
262  incoder_ = dynamic_cast<const HcaluLUTTPGCoder*>(incoder);
263  outcoder_ = outcoder;
264  conditions_ = conditions;
265 
266  theSumMap.clear();
267  theSatMap.clear();
268  theTowerMapFGSum.clear();
269  HF_Veto.clear();
270  fgMap_.clear();
271  fgUpgradeMap_.clear();
272  fgUpgradeTDCMap_.clear();
273  theHFDetailMap.clear();
274  theHFUpgradeDetailMap.clear();
275 
276  // Add all digi collections
277  addDigis(digis...);
278 
279  // Prepare the fine-grain calculation algorithm for HB/HE
280  int version = 0;
281  if (upgrade_he_ or upgrade_hb_)
282  version = conditions_->getHcalTPParameters()->getFGVersionHBHE();
283  if (override_parameters_.exists("FGVersionHBHE"))
284  version = override_parameters_.getParameter<uint32_t>("FGVersionHBHE");
285  HcalFinegrainBit fg_algo(version);
286 
287  // VME produces additional bits on the front used by lumi but not the
288  // trigger, this shift corrects those out by right shifting over them.
289  for (auto& item : theSumMap) {
290  result.push_back(HcalTriggerPrimitiveDigi(item.first));
291  HcalTrigTowerDetId detId(item.second.id());
292  if (detId.ietaAbs() >= theTrigTowerGeometry->firstHFTower(detId.version())) {
293  if (detId.version() == 0) {
294  analyzeHF(item.second, result.back(), RCTScaleShift);
295  } else if (detId.version() == 1) {
296  if (upgrade_hf_)
297  analyzeHFQIE10(item.second, result.back(), NCTScaleShift, LongvrsShortCut);
298  else
299  analyzeHF2016(item.second, result.back(), NCTScaleShift, LongvrsShortCut);
300  } else {
301  // Things are going to go poorly
302  }
303  } else {
304  // Determine which energy reconstruction path to take based on the
305  // fine-grain availability:
306  // * QIE8 TP add entries into fgMap_
307  // * QIE11 TP add entries into fgUpgradeMap_
308  // (not for tower 16 unless HB is upgraded, too)
309  if (fgMap_.find(item.first) != fgMap_.end()) {
310  analyze(item.second, result.back());
311  } else if (fgUpgradeMap_.find(item.first) != fgUpgradeMap_.end()) {
312  SatMap::iterator item_sat = theSatMap.find(detId);
313  if (item_sat == theSatMap.end())
314  analyzeQIE11(item.second, std::vector<bool>(), result.back(), fg_algo);
315  else
316  analyzeQIE11(item.second, item_sat->second, result.back(), fg_algo);
317  }
318  }
319  }
320 
321  // Free up some memory
322  theSumMap.clear();
323  theSatMap.clear();
324  theTowerMapFGSum.clear();
325  HF_Veto.clear();
326  fgMap_.clear();
327  fgUpgradeMap_.clear();
328  fgUpgradeTDCMap_.clear();
329  theHFDetailMap.clear();
330  theHFUpgradeDetailMap.clear();
331 
332  return;
333 }
HcalTriggerPrimitiveAlgo::analyze
void analyze(IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result)
adds the actual digis
Definition: HcalTriggerPrimitiveAlgo.cc:319
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HcalTriggerPrimitiveAlgo::theHFUpgradeDetailMap
HFUpgradeDetailMap theHFUpgradeDetailMap
Definition: HcalTriggerPrimitiveAlgo.h:191
HcalTPParameters::getFGVersionHBHE
int getFGVersionHBHE() const
get FineGrain Algorithm Version for HBHE
Definition: HcalTPParameters.h:19
HcalTriggerPrimitiveAlgo::analyzeHF
void analyzeHF(IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result, const int hf_lumi_shift)
Definition: HcalTriggerPrimitiveAlgo.cc:565
mps_fire.result
result
Definition: mps_fire.py:311
HcalTriggerPrimitiveAlgo::analyzeHF2016
void analyzeHF2016(const IntegerCaloSamples &SAMPLES, HcalTriggerPrimitiveDigi &result, const int HF_LUMI_SHIFT, const HcalFeatureBit *HCALFEM)
Definition: HcalTriggerPrimitiveAlgo.cc:619
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:681
HcalTriggerPrimitiveAlgo::analyzeHFQIE10
void analyzeHFQIE10(const IntegerCaloSamples &SAMPLES, HcalTriggerPrimitiveDigi &result, const int HF_LUMI_SHIFT, const HcalFeatureBit *HCALFEM)
Definition: HcalTriggerPrimitiveAlgo.cc:708
HcalTriggerPrimitiveAlgo::outcoder_
const HcalTPGCompressor * outcoder_
Definition: HcalTriggerPrimitiveAlgo.h:134
B2GTnPMonitor_cfi.item
item
Definition: B2GTnPMonitor_cfi.py:148
HcalTrigTowerGeometry::firstHFTower
int firstHFTower(int version) const
Definition: HcalTrigTowerGeometry.h:18
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
HcalTriggerPrimitiveAlgo::conditions_
const HcalDbService * conditions_
Definition: HcalTriggerPrimitiveAlgo.h:135
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
BeamSplash_cfg.version
version
Definition: BeamSplash_cfg.py:45
HcalDbService::getHcalTPParameters
const HcalTPParameters * getHcalTPParameters() const
Definition: HcalDbService.cc:413
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
HcalTriggerPrimitiveAlgo::analyzeQIE11
void analyzeQIE11(IntegerCaloSamples &samples, std::vector< bool > sample_saturation, HcalTriggerPrimitiveDigi &result, const HcalFinegrainBit &fg_algo)
Definition: HcalTriggerPrimitiveAlgo.cc:401
HcalTriggerPrimitiveAlgo::addDigis
void addDigis(const T &collection, const Args &... digis)
Definition: HcalTriggerPrimitiveAlgo.h:57
HcalTriggerPrimitiveAlgo::incoder_
const HcaluLUTTPGCoder * incoder_
Definition: HcalTriggerPrimitiveAlgo.h:133

◆ runFEFormatError()

void HcalTriggerPrimitiveAlgo::runFEFormatError ( const FEDRawDataCollection *  rawraw,
const HcalElectronicsMap *  emap,
HcalTrigPrimDigiCollection &  result 
)

Definition at line 834 of file HcalTriggerPrimitiveAlgo.cc.

References HcalHTRData::check(), FEDRawData::data(), hcalRecHitTable_cff::detId, runTauDisplay::eid, FEDRawDataCollection::FEDData(), HcalHTRData::getErrorsWord(), HcalDCCHeader::getSourceId(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), HcalBarrel, HcalEndcap, HcalForward, HcalHTRData::htrSlot(), HcalHTRData::htrTopBottom(), mps_fire::i, HcalHTRData::isHistogramEvent(), HcalElectronicsMap::lookup(), FEDNumbering::MAXHCALFEDID, FEDNumbering::MINHCALFEDID, HcalHTRData::readoutVMECrateId(), mps_fire::result, FEDRawData::size(), HcalDCCHeader::SPIGOT_COUNT, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), cmsswSequenceInfo::tp, and validateGeometry_cfg::valid.

Referenced by HcalTrigPrimDigiProducer::produce().

836  {
837  std::set<uint32_t> FrontEndErrors;
838 
839  for (int i = FEDNumbering::MINHCALFEDID; i <= FEDNumbering::MAXHCALFEDID; ++i) {
840  const FEDRawData& raw = rawraw->FEDData(i);
841  if (raw.size() < 12)
842  continue;
843  const HcalDCCHeader* dccHeader = (const HcalDCCHeader*)(raw.data());
844  if (!dccHeader)
845  continue;
846  HcalHTRData htr;
847  for (int spigot = 0; spigot < HcalDCCHeader::SPIGOT_COUNT; spigot++) {
848  if (!dccHeader->getSpigotPresent(spigot))
849  continue;
850  dccHeader->getSpigotData(spigot, htr, raw.size());
851  int dccid = dccHeader->getSourceId();
852  int errWord = htr.getErrorsWord() & 0x1FFFF;
853  bool HTRError = (!htr.check() || htr.isHistogramEvent() || (errWord & 0x800) != 0);
854 
855  if (HTRError) {
856  bool valid = false;
857  for (int fchan = 0; fchan < 3 && !valid; fchan++) {
858  for (int fib = 0; fib < 9 && !valid; fib++) {
859  HcalElectronicsId eid(fchan, fib, spigot, dccid - FEDNumbering::MINHCALFEDID);
860  eid.setHTR(htr.readoutVMECrateId(), htr.htrSlot(), htr.htrTopBottom());
861  DetId detId = emap->lookup(eid);
862  if (detId.null())
863  continue;
864  HcalSubdetector subdet = (HcalSubdetector(detId.subdetId()));
865  if (detId.det() != 4 || (subdet != HcalBarrel && subdet != HcalEndcap && subdet != HcalForward))
866  continue;
867  std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
868  for (std::vector<HcalTrigTowerDetId>::const_iterator triggerId = ids.begin(); triggerId != ids.end();
869  ++triggerId) {
870  FrontEndErrors.insert(triggerId->rawId());
871  }
872  //valid = true;
873  }
874  }
875  }
876  }
877  }
878 
879  // Loop over TP collection
880  // Set TP to zero if there is FE Format Error
881  HcalTriggerPrimitiveSample zeroSample(0);
882  for (HcalTrigPrimDigiCollection::iterator tp = result.begin(); tp != result.end(); ++tp) {
883  if (FrontEndErrors.find(tp->id().rawId()) != FrontEndErrors.end()) {
884  for (int i = 0; i < tp->size(); ++i)
885  tp->setSample(i, zeroSample);
886  }
887  }
888 }
HcalHTRData::isHistogramEvent
bool isHistogramEvent() const
Is this event a histogram event? (do not call standard unpack in this case!!!!!)
Definition: HcalHTRData.cc:409
mps_fire.i
i
Definition: mps_fire.py:429
HcalDCCHeader::getSpigotData
int getSpigotData(int nspigot, HcalHTRData &decodeTool, int validSize) const
Definition: HcalDCCHeader.cc:27
mps_fire.result
result
Definition: mps_fire.py:311
HcalElectronicsMap::lookup
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
Definition: HcalElectronicsMap.cc:70
FEDRawData::size
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:45
HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
const HcalTrigTowerGeometry * theTrigTowerGeometry
Definition: HcalTriggerPrimitiveAlgo.h:166
HcalHTRData::htrSlot
unsigned int htrSlot() const
HcalElectronicsId-style HTR slot.
Definition: HcalHTRData.cc:365
HcalHTRData::htrTopBottom
unsigned int htrTopBottom() const
HcalElectronicsId-style HTR top/bottom (1=top/0=bottom)
Definition: HcalHTRData.cc:369
validateGeometry_cfg.valid
valid
Definition: validateGeometry_cfg.py:21
HcalDCCHeader::getSourceId
int getSourceId() const
Definition: HcalDCCHeader.h:32
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
FEDRawDataCollection::FEDData
const FEDRawData & FEDData(int fedid) const
retrieve data for fed
Definition: FEDRawDataCollection.cc:19
edm::SortedCollection::iterator
std::vector< T >::iterator iterator
Definition: SortedCollection.h:81
DetId
Definition: DetId.h:17
HcalDCCHeader::SPIGOT_COUNT
static const int SPIGOT_COUNT
Definition: HcalDCCHeader.h:19
HcalHTRData::getErrorsWord
unsigned int getErrorsWord() const
Get the errors word.
Definition: HcalHTRData.h:162
HcalDCCHeader::getSpigotPresent
bool getSpigotPresent(unsigned int nspigot) const
Read the "PRESENT" bit for this spigot.
Definition: HcalDCCHeader.h:106
FEDRawData::data
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:24
HcalHTRData::readoutVMECrateId
unsigned int readoutVMECrateId() const
HcalElectronicsId-style VME crate number.
Definition: HcalHTRData.cc:373
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
HcalTrigTowerGeometry::towerIds
std::vector< HcalTrigTowerDetId > towerIds(const HcalDetId &cellId) const
the mapping to and from DetIds
Definition: HcalTrigTowerGeometry.cc:15
HcalHTRData::check
bool check() const
Check for a good event Requires a minimum length, matching wordcount and length, not an empty event...
Definition: HcalHTRData.cc:63

◆ runZS()

void HcalTriggerPrimitiveAlgo::runZS ( HcalTrigPrimDigiCollection &  tp)

Definition at line 818 of file HcalTriggerPrimitiveAlgo.cc.

References mps_fire::i, mps_fire::result, cmsswSequenceInfo::tp, and ZS_threshold_I_.

Referenced by HcalTrigPrimDigiProducer::produce().

818  {
819  for (HcalTrigPrimDigiCollection::iterator tp = result.begin(); tp != result.end(); ++tp) {
820  bool ZS = true;
821  for (int i = 0; i < tp->size(); ++i) {
822  if (tp->sample(i).compressedEt() > ZS_threshold_I_) {
823  ZS = false;
824  break;
825  }
826  }
827  if (ZS)
828  tp->setZSInfo(false, true);
829  else
830  tp->setZSInfo(true, false);
831  }
832 }
mps_fire.i
i
Definition: mps_fire.py:429
mps_fire.result
result
Definition: mps_fire.py:311
edm::SortedCollection::iterator
std::vector< T >::iterator iterator
Definition: SortedCollection.h:81

◆ setFixSaturationFlag()

void HcalTriggerPrimitiveAlgo::setFixSaturationFlag ( bool  fix_saturation)

Definition at line 79 of file HcalTriggerPrimitiveAlgo.cc.

References fix_saturation_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

79 { fix_saturation_ = fix_saturation; }

◆ setNCTScaleShift()

void HcalTriggerPrimitiveAlgo::setNCTScaleShift ( int  shift)

Definition at line 1015 of file HcalTriggerPrimitiveAlgo.cc.

References NCTScaleShift, and edm::shift.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

1015 { NCTScaleShift = shift; }
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:7

◆ setNumFilterPresamplesHBQIE11()

void HcalTriggerPrimitiveAlgo::setNumFilterPresamplesHBQIE11 ( int  presamples)
inline

Definition at line 87 of file HcalTriggerPrimitiveAlgo.h.

References numberOfFilterPresamplesHBQIE11_, and hcalTTPDigis_cfi::presamples.

Referenced by HcalTrigPrimDigiProducer::beginRun().

87 { numberOfFilterPresamplesHBQIE11_ = presamples; }
hcalTTPDigis_cfi.presamples
presamples
Definition: hcalTTPDigis_cfi.py:9

◆ setNumFilterPresamplesHEQIE11()

void HcalTriggerPrimitiveAlgo::setNumFilterPresamplesHEQIE11 ( int  presamples)
inline

Definition at line 89 of file HcalTriggerPrimitiveAlgo.h.

References numberOfFilterPresamplesHEQIE11_, and hcalTTPDigis_cfi::presamples.

Referenced by HcalTrigPrimDigiProducer::beginRun().

89 { numberOfFilterPresamplesHEQIE11_ = presamples; }
hcalTTPDigis_cfi.presamples
presamples
Definition: hcalTTPDigis_cfi.py:9

◆ setPeakFinderAlgorithm()

void HcalTriggerPrimitiveAlgo::setPeakFinderAlgorithm ( int  algo)

Definition at line 1009 of file HcalTriggerPrimitiveAlgo.cc.

References Exception, and peak_finder_algorithm_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

1009  {
1010  if (algo <= 0 || algo > 2)
1011  throw cms::Exception("ERROR: Only algo 1 & 2 are supported.") << std::endl;
1012  peak_finder_algorithm_ = algo;
1013 }

◆ setRCTScaleShift()

void HcalTriggerPrimitiveAlgo::setRCTScaleShift ( int  shift)

Definition at line 1017 of file HcalTriggerPrimitiveAlgo.cc.

References RCTScaleShift, and edm::shift.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

1017 { RCTScaleShift = shift; }
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:7

◆ setUpgradeFlags()

void HcalTriggerPrimitiveAlgo::setUpgradeFlags ( bool  hb,
bool  he,
bool  hf 
)

Definition at line 73 of file HcalTriggerPrimitiveAlgo.cc.

References hcalSimParameters_cfi::hb, hcalSimParameters_cfi::he, photonIsolationHIProducer_cfi::hf, upgrade_hb_, upgrade_he_, and upgrade_hf_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

73  {
74  upgrade_hb_ = hb;
75  upgrade_he_ = he;
76  upgrade_hf_ = hf;
77 }
hcalSimParameters_cfi.hb
hb
Definition: hcalSimParameters_cfi.py:60
hcalSimParameters_cfi.he
he
Definition: hcalSimParameters_cfi.py:79
photonIsolationHIProducer_cfi.hf
hf
Definition: photonIsolationHIProducer_cfi.py:9

◆ setWeightQIE11()

void HcalTriggerPrimitiveAlgo::setWeightQIE11 ( int  aieta,
int  weight 
)

Definition at line 1003 of file HcalTriggerPrimitiveAlgo.cc.

References mps_merge::weight, and weightsQIE11_.

Referenced by HcalTrigPrimDigiProducer::beginRun().

1003  {
1004  // Simple map of |ieta| in HBHE to weight
1005  // Only one weight for SOI-1 TS
1006  weightsQIE11_[aieta] = {{weight, 255}};
1007 }
HcalTriggerPrimitiveAlgo::weightsQIE11_
std::array< std::array< int, 2 >, 29 > weightsQIE11_
Definition: HcalTriggerPrimitiveAlgo.h:139
mps_merge.weight
weight
Definition: mps_merge.py:88

◆ setWeightsQIE11()

void HcalTriggerPrimitiveAlgo::setWeightsQIE11 ( const edm::ParameterSet &  weightsQIE11)

Definition at line 993 of file HcalTriggerPrimitiveAlgo.cc.

References findQualityFiles::v, hcaltpdigi_cfi::weightsQIE11, and weightsQIE11_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

993  {
994  // Names are just abs(ieta) for HBHE
995  std::vector<std::string> ietaStrs = weightsQIE11.getParameterNames();
996  for (auto& ietaStr : ietaStrs) {
997  // Strip off "ieta" part of key and just use integer value in map
998  auto const& v = weightsQIE11.getParameter<std::vector<int>>(ietaStr);
999  weightsQIE11_[std::stoi(ietaStr.substr(4))] = {{v[0], v[1]}};
1000  }
1001 }
HcalTriggerPrimitiveAlgo::weightsQIE11_
std::array< std::array< int, 2 >, 29 > weightsQIE11_
Definition: HcalTriggerPrimitiveAlgo.h:139
hcaltpdigi_cfi.weightsQIE11
weightsQIE11
hardware algo
Definition: hcaltpdigi_cfi.py:16

◆ validChannel()

bool HcalTriggerPrimitiveAlgo::validChannel ( const QIE10DataFrame &  digi,
int  ts 
) const
private

Definition at line 696 of file HcalTriggerPrimitiveAlgo.cc.

References conditions_, HcalDbService::getHcalTPChannelParameter(), HcalTPChannelParameter::getMask(), QIE10DataFrame::id(), QIE10DataFrame::linkError(), gpuClustering::pixelStatus::mask, and QIE10DataFrame::samples().

Referenced by addSignal().

696  {
697  // channels with invalid data should not contribute to the sum
698  if (digi.linkError() || ts >= digi.samples() || !digi[ts].ok())
699  return false;
700 
701  auto mask = conditions_->getHcalTPChannelParameter(HcalDetId(digi.id()))->getMask();
702  if (mask)
703  return false;
704 
705  return true;
706 }
QIE10DataFrame::id
constexpr edm::DataFrame::id_type id() const
Definition: QIE10DataFrame.h:77
HcalTPChannelParameter::getMask
uint32_t getMask() const
get mask for channel validity and self trigger information
Definition: HcalTPChannelParameter.h:10
gpuClustering::pixelStatus::mask
constexpr uint32_t mask
Definition: gpuClustering.h:26
HcalTriggerPrimitiveAlgo::conditions_
const HcalDbService * conditions_
Definition: HcalTriggerPrimitiveAlgo.h:135
HcalDbService::getHcalTPChannelParameter
const HcalTPChannelParameter * getHcalTPChannelParameter(const HcalGenericDetId &fId, bool throwOnFail=true) const
Definition: HcalDbService.cc:391
QIE10DataFrame::linkError
constexpr bool linkError() const
Definition: QIE10DataFrame.h:101
QIE10DataFrame::samples
constexpr int samples() const
total number of samples in the digi
Definition: QIE10DataFrame.h:86

◆ validUpgradeFG()

bool HcalTriggerPrimitiveAlgo::validUpgradeFG ( const HcalTrigTowerDetId &  id,
int  depth 
) const
private

Definition at line 900 of file HcalTriggerPrimitiveAlgo.cc.

References hcalRecHitTable_cff::depth, HBHE_OVERLAP_TOWER, LAST_FINEGRAIN_DEPTH, LAST_FINEGRAIN_TOWER, and upgrade_hb_.

Referenced by addUpgradeFG().

900  {
901  if (depth > LAST_FINEGRAIN_DEPTH)
902  return false;
903  if (id.ietaAbs() > LAST_FINEGRAIN_TOWER)
904  return false;
905  if (id.ietaAbs() == HBHE_OVERLAP_TOWER and not upgrade_hb_)
906  return false;
907  return true;
908 }
hcalRecHitTable_cff.depth
depth
Definition: hcalRecHitTable_cff.py:17

Member Data Documentation

◆ conditions_

const HcalDbService* HcalTriggerPrimitiveAlgo::conditions_
private

Definition at line 135 of file HcalTriggerPrimitiveAlgo.h.

Referenced by passTDC(), run(), and validChannel().

◆ FG_HF_thresholds_

std::vector<uint32_t> HcalTriggerPrimitiveAlgo::FG_HF_thresholds_
private

Definition at line 142 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2016(), and analyzeHFQIE10().

◆ FG_threshold_

uint32_t HcalTriggerPrimitiveAlgo::FG_threshold_
private

Definition at line 141 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF().

◆ fgMap_

FGbitMap HcalTriggerPrimitiveAlgo::fgMap_
private

Definition at line 211 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addFG(), analyze(), and run().

◆ fgUpgradeMap_

FGUpgradeMap HcalTriggerPrimitiveAlgo::fgUpgradeMap_
private

Definition at line 215 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addUpgradeFG(), analyzeQIE11(), and run().

◆ fgUpgradeTDCMap_

FGUpgradeTDCMap HcalTriggerPrimitiveAlgo::fgUpgradeTDCMap_
private

Definition at line 219 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addUpgradeTDCFG(), analyzeQIE11(), and run().

◆ FIRST_DEPTH7_TOWER

const int HcalTriggerPrimitiveAlgo::FIRST_DEPTH7_TOWER = 26
staticprivate

Definition at line 236 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needUpgradeID().

◆ FIRST_FINEGRAIN_TOWER

const int HcalTriggerPrimitiveAlgo::FIRST_FINEGRAIN_TOWER = 30
staticprivate

Definition at line 241 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2016(), and analyzeHFQIE10().

◆ fix_saturation_

bool HcalTriggerPrimitiveAlgo::fix_saturation_ = false
private

Definition at line 225 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeQIE11(), and setFixSaturationFlag().

◆ HBHE_OVERLAP_TOWER

const int HcalTriggerPrimitiveAlgo::HBHE_OVERLAP_TOWER = 16
staticprivate

Definition at line 235 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needLegacyFG(), and validUpgradeFG().

◆ HF_Veto

TowerMapVeto HcalTriggerPrimitiveAlgo::HF_Veto
private

Definition at line 208 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF(), and run().

◆ incoder_

const HcaluLUTTPGCoder* HcalTriggerPrimitiveAlgo::incoder_
private

Definition at line 133 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), addUpgradeTDCFG(), and run().

◆ LAST_FINEGRAIN_DEPTH

const int HcalTriggerPrimitiveAlgo::LAST_FINEGRAIN_DEPTH = 6
staticprivate

Definition at line 237 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needUpgradeID(), and validUpgradeFG().

◆ LAST_FINEGRAIN_TOWER

const int HcalTriggerPrimitiveAlgo::LAST_FINEGRAIN_TOWER = 28
staticprivate

Definition at line 238 of file HcalTriggerPrimitiveAlgo.h.

Referenced by validUpgradeFG().

◆ latency_

int HcalTriggerPrimitiveAlgo::latency_
private

Definition at line 140 of file HcalTriggerPrimitiveAlgo.h.

◆ LongvrsShortCut

HcalFeatureBit* HcalTriggerPrimitiveAlgo::LongvrsShortCut
private

Definition at line 206 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run().

◆ minSignalThreshold_

uint32_t HcalTriggerPrimitiveAlgo::minSignalThreshold_
private

Definition at line 152 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal().

◆ NCTScaleShift

int HcalTriggerPrimitiveAlgo::NCTScaleShift
private

Definition at line 154 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setNCTScaleShift().

◆ numberOfFilterPresamplesHBQIE11_

int HcalTriggerPrimitiveAlgo::numberOfFilterPresamplesHBQIE11_
private

Definition at line 147 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11(), and setNumFilterPresamplesHBQIE11().

◆ numberOfFilterPresamplesHEQIE11_

int HcalTriggerPrimitiveAlgo::numberOfFilterPresamplesHEQIE11_
private

Definition at line 148 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11(), and setNumFilterPresamplesHEQIE11().

◆ numberOfPresamples_

int HcalTriggerPrimitiveAlgo::numberOfPresamples_
private

Definition at line 146 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyzeQIE11(), and HcalTriggerPrimitiveAlgo().

◆ numberOfPresamplesHF_

int HcalTriggerPrimitiveAlgo::numberOfPresamplesHF_
private

Definition at line 150 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF(), analyzeHF2016(), analyzeHFQIE10(), and HcalTriggerPrimitiveAlgo().

◆ numberOfSamples_

int HcalTriggerPrimitiveAlgo::numberOfSamples_
private

Definition at line 145 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyzeQIE11(), and HcalTriggerPrimitiveAlgo().

◆ numberOfSamplesHF_

int HcalTriggerPrimitiveAlgo::numberOfSamplesHF_
private

Definition at line 149 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF(), analyzeHF2016(), analyzeHFQIE10(), and HcalTriggerPrimitiveAlgo().

◆ outcoder_

const HcalTPGCompressor* HcalTriggerPrimitiveAlgo::outcoder_
private

Definition at line 134 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyzeHF(), analyzeHF2016(), analyzeHFQIE10(), analyzeQIE11(), and run().

◆ override_adc_hf_

bool HcalTriggerPrimitiveAlgo::override_adc_hf_ = false
private

Definition at line 229 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ override_adc_hf_value_

uint32_t HcalTriggerPrimitiveAlgo::override_adc_hf_value_
private

Definition at line 230 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ override_parameters_

edm::ParameterSet HcalTriggerPrimitiveAlgo::override_parameters_
private

Definition at line 227 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and run().

◆ override_tdc_hf_

bool HcalTriggerPrimitiveAlgo::override_tdc_hf_ = false
private

Definition at line 231 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ override_tdc_hf_value_

unsigned long long HcalTriggerPrimitiveAlgo::override_tdc_hf_value_
private

Definition at line 232 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ peak_finder_algorithm_

int HcalTriggerPrimitiveAlgo::peak_finder_algorithm_
private

Definition at line 161 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and setPeakFinderAlgorithm().

◆ peakfind_

bool HcalTriggerPrimitiveAlgo::peakfind_
private

Definition at line 137 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and HcalTriggerPrimitiveAlgo().

◆ PMT_NoiseThreshold_

uint32_t HcalTriggerPrimitiveAlgo::PMT_NoiseThreshold_
private

Definition at line 153 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF().

◆ QIE10_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE10_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE10_LUT_BITMASK
staticprivate

Definition at line 244 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHFQIE10().

◆ QIE10_MAX_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE10_MAX_LINEARIZATION_ET = 0x7FF
staticprivate

Definition at line 247 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHFQIE10().

◆ QIE11_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE11_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE11_LUT_BITMASK
staticprivate

Definition at line 245 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal().

◆ QIE11_MAX_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE11_MAX_LINEARIZATION_ET = 0x7FF
staticprivate

Definition at line 248 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11().

◆ QIE8_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE8_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE8_LUT_BITMASK
staticprivate

Definition at line 243 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), analyzeHF(), and analyzeHF2016().

◆ RCTScaleShift

int HcalTriggerPrimitiveAlgo::RCTScaleShift
private

Definition at line 155 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setRCTScaleShift().

◆ theHFDetailMap

HFDetailMap HcalTriggerPrimitiveAlgo::theHFDetailMap
private

Definition at line 181 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF2016(), and run().

◆ theHFUpgradeDetailMap

HFUpgradeDetailMap HcalTriggerPrimitiveAlgo::theHFUpgradeDetailMap
private

Definition at line 191 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHFQIE10(), and run().

◆ theSatMap

SatMap HcalTriggerPrimitiveAlgo::theSatMap
private

Definition at line 172 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), and run().

◆ theSumMap

SumMap HcalTriggerPrimitiveAlgo::theSumMap
private

Definition at line 169 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), and run().

◆ theThreshold

double HcalTriggerPrimitiveAlgo::theThreshold
private

Definition at line 136 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and analyzeQIE11().

◆ theTowerMapFGSum

TowerMapFGSum HcalTriggerPrimitiveAlgo::theTowerMapFGSum
private

Definition at line 195 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), analyzeHF(), and run().

◆ theTrigTowerGeometry

const HcalTrigTowerGeometry* HcalTriggerPrimitiveAlgo::theTrigTowerGeometry
private

Definition at line 166 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), addUpgradeTDCFG(), analyzeQIE11(), run(), and runFEFormatError().

◆ upgrade_hb_

bool HcalTriggerPrimitiveAlgo::upgrade_hb_ = false
private

Definition at line 221 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needLegacyFG(), run(), setUpgradeFlags(), and validUpgradeFG().

◆ upgrade_he_

bool HcalTriggerPrimitiveAlgo::upgrade_he_ = false
private

Definition at line 222 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

◆ upgrade_hf_

bool HcalTriggerPrimitiveAlgo::upgrade_hf_ = false
private

Definition at line 223 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

◆ useTDCInMinBiasBits_

bool HcalTriggerPrimitiveAlgo::useTDCInMinBiasBits_
private

Definition at line 151 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHFQIE10().

◆ weights_

std::vector<double> HcalTriggerPrimitiveAlgo::weights_
private

Definition at line 138 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze().

◆ weightsQIE11_

std::array<std::array<int, 2>, 29> HcalTriggerPrimitiveAlgo::weightsQIE11_
private

Definition at line 139 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11(), setWeightQIE11(), and setWeightsQIE11().

◆ ZS_threshold_

uint32_t HcalTriggerPrimitiveAlgo::ZS_threshold_
private

Definition at line 143 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo().

◆ ZS_threshold_I_

int HcalTriggerPrimitiveAlgo::ZS_threshold_I_
private

Definition at line 144 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo(), and runZS().

Generated for CMSSW Reference Manual by doxygen 1.8.14