#include <HcalTriggerPrimitiveAlgo.h>

Classes
struct	HFDetails

struct	HFUpgradeDetails

Public Member Functions
template<typename D >
void	addDigis (const HcalDataFrameContainer< D > &collection)

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

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

	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, double 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 IntegerCaloSamples &samples)

void	addSignal (const QIE10DataFrame &frame)

void	addSignal (const QIE11DataFrame &frame)

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< double, 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.

     : incoder_(nullptr),
       outcoder_(nullptr),
       theThreshold(0),
       peakfind_(pf),
       weights_(w),
       latency_(latency),
       FG_threshold_(FG_threshold),
       FG_HF_thresholds_(FG_HF_thresholds),
       ZS_threshold_(ZS_threshold),
       numberOfSamples_(numberOfSamples),
       numberOfPresamples_(numberOfPresamples),
       numberOfFilterPresamplesHBQIE11_(numberOfFilterPresamplesHBQIE11),
       numberOfFilterPresamplesHEQIE11_(numberOfFilterPresamplesHEQIE11),
       numberOfSamplesHF_(numberOfSamplesHF),
       numberOfPresamplesHF_(numberOfPresamplesHF),
       useTDCInMinBiasBits_(useTDCInMinBiasBits),
       minSignalThreshold_(minSignalThreshold),
       PMT_NoiseThreshold_(PMT_NoiseThreshold),
       NCTScaleShift(0),
       RCTScaleShift(0),
       peak_finder_algorithm_(2),
       override_parameters_() {
   //No peak finding setting (for Fastsim)
   if (!peakfind_) {
     numberOfSamples_ = 1;
     numberOfPresamples_ = 0;
     numberOfSamplesHF_ = 1;
     numberOfPresamplesHF_ = 0;
   }
   // Switch to integer for comparisons - remove compiler warning
   ZS_threshold_I_ = ZS_threshold_;
 }

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

◆ ~HcalTriggerPrimitiveAlgo()

HcalTriggerPrimitiveAlgo::~HcalTriggerPrimitiveAlgo ( )

Definition at line 71 of file HcalTriggerPrimitiveAlgo.cc.

71 {}

Member Function Documentation

◆ addDigis() [1/3]

template<typename D >

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

inline

Definition at line 70 of file HcalTriggerPrimitiveAlgo.h.

                                                              {
     for (auto i = collection.begin(); i != collection.end(); ++i) {
       D digi(*i);
       addSignal(digi);
     }
   };

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

◆ addDigis() [2/3]

template<typename T >

void HcalTriggerPrimitiveAlgo::addDigis ( const T & collection )

inline

Definition at line 63 of file HcalTriggerPrimitiveAlgo.h.

                                      {
     for (const auto& digi : collection) {
       addSignal(digi);
     }
   };

References addSignal(), and universalConfigTemplate::collection.

◆ addDigis() [3/3]

template<typename T , typename... Args>

void HcalTriggerPrimitiveAlgo::addDigis	(	const T &	collection,
		const Args &...	digis
	)

inline

Definition at line 57 of file HcalTriggerPrimitiveAlgo.h.

                                                            {
     addDigis(collection);
     addDigis(digis...);
   };

References universalConfigTemplate::collection.

Referenced by run().

◆ addFG()

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

private

Definition at line 833 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                        {
   FGbitMap::iterator itr = fgMap_.find(id);
   if (itr != fgMap_.end()) {
     std::vector<bool>& _msb = itr->second;
     for (size_t i = 0; i < msb.size(); ++i)
       _msb[i] = _msb[i] || msb[i];
   } else
     fgMap_[id] = msb;
 }

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

Referenced by addSignal(), and addUpgradeFG().

◆ addSignal() [1/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const HBHEDataFrame & frame )

private

adds the signal to the map

Definition at line 94 of file HcalTriggerPrimitiveAlgo.cc.

                                                                    {
   // TODO: Need to add support for seperate 28, 29 in HE
   //Hack for 300_pre10, should be removed.
   if (frame.id().depth() == 5)
     return;
  
   std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
   assert(ids.size() == 1 || ids.size() == 2);
   IntegerCaloSamples samples1(ids[0], int(frame.size()));
  
   samples1.setPresamples(frame.presamples());
   incoder_->adc2Linear(frame, samples1);
  
   std::vector<bool> msb;
   incoder_->lookupMSB(frame, msb);
  
   if (ids.size() == 2) {
     // make a second trigprim for the other one, and split the energy
     IntegerCaloSamples samples2(ids[1], samples1.size());
     for (int i = 0; i < samples1.size(); ++i) {
       samples1[i] = uint32_t(samples1[i] * 0.5);
       samples2[i] = samples1[i];
     }
     samples2.setPresamples(frame.presamples());
     addSignal(samples2);
     addFG(ids[1], msb);
   }
   addSignal(samples1);
   addFG(ids[0], msb);
 }

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

Referenced by addDigis(), and addSignal().

◆ addSignal() [2/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const HFDataFrame & frame )

private

Definition at line 125 of file HcalTriggerPrimitiveAlgo.cc.

                                                                  {
   if (frame.id().depth() == 1 || frame.id().depth() == 2) {
     std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(frame.id());
     std::vector<HcalTrigTowerDetId>::const_iterator it;
     for (it = ids.begin(); it != ids.end(); ++it) {
       HcalTrigTowerDetId trig_tower_id = *it;
       IntegerCaloSamples samples(trig_tower_id, frame.size());
       samples.setPresamples(frame.presamples());
       incoder_->adc2Linear(frame, samples);
  
       // Don't add to final collection yet
       // HF PMT veto sum is calculated in analyzerHF()
       IntegerCaloSamples zero_samples(trig_tower_id, frame.size());
       zero_samples.setPresamples(frame.presamples());
       addSignal(zero_samples);
  
       // Pre-LS1 Configuration
       if (trig_tower_id.version() == 0) {
         // Mask off depths: fgid is the same for both depths
         uint32_t fgid = (frame.id().maskDepth());
  
         if (theTowerMapFGSum.find(trig_tower_id) == theTowerMapFGSum.end()) {
           SumFGContainer sumFG;
           theTowerMapFGSum.insert(std::pair<HcalTrigTowerDetId, SumFGContainer>(trig_tower_id, sumFG));
         }
  
         SumFGContainer& sumFG = theTowerMapFGSum[trig_tower_id];
         SumFGContainer::iterator sumFGItr;
         for (sumFGItr = sumFG.begin(); sumFGItr != sumFG.end(); ++sumFGItr) {
           if (sumFGItr->id() == fgid) {
             break;
           }
         }
         // If find
         if (sumFGItr != sumFG.end()) {
           for (int i = 0; i < samples.size(); ++i) {
             (*sumFGItr)[i] += samples[i];
           }
         } else {
           //Copy samples (change to fgid)
           IntegerCaloSamples sumFGSamples(DetId(fgid), samples.size());
           sumFGSamples.setPresamples(samples.presamples());
           for (int i = 0; i < samples.size(); ++i) {
             sumFGSamples[i] = samples[i];
           }
           sumFG.push_back(sumFGSamples);
         }
  
         // set veto to true if Long or Short less than threshold
         if (HF_Veto.find(fgid) == HF_Veto.end()) {
           vector<bool> vetoBits(samples.size(), false);
           HF_Veto[fgid] = vetoBits;
         }
         for (int i = 0; i < samples.size(); ++i) {
           if (samples[i] < minSignalThreshold_) {
             HF_Veto[fgid][i] = true;
           }
         }
       }
       // HF 1x1
       else if (trig_tower_id.version() == 1) {
         uint32_t fgid = (frame.id().maskDepth());
         HFDetails& details = theHFDetailMap[trig_tower_id][fgid];
         // Check the frame type to determine long vs short
         if (frame.id().depth() == 1) {  // Long
           details.long_fiber = samples;
           details.LongDigi = frame;
         } else if (frame.id().depth() == 2) {  // Short
           details.short_fiber = samples;
           details.ShortDigi = frame;
         } else {
           // Neither long nor short... So we have no idea what to do
           edm::LogWarning("HcalTPAlgo") << "Unable to figure out what to do with data frame for " << frame.id();
           return;
         }
       }
       // Uh oh, we are in a bad/unknown state! Things will start crashing.
       else {
         return;
       }
     }
   }
 }

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().

◆ addSignal() [3/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const IntegerCaloSamples & samples )

private

Definition at line 282 of file HcalTriggerPrimitiveAlgo.cc.

                                                                           {
   HcalTrigTowerDetId id(samples.id());
   SumMap::iterator itr = theSumMap.find(id);
  
   if (itr == theSumMap.end()) {
     theSumMap.insert(std::make_pair(id, samples));
   } else {
     // wish CaloSamples had a +=
     for (int i = 0; i < samples.size(); ++i) {
       (itr->second)[i] += samples[i];
     }
   }
  
   // if fix_saturation == true, keep track of tower with saturated input LUT
   if (fix_saturation_) {
     SatMap::iterator itr_sat = theSatMap.find(id);
  
     assert((itr == theSumMap.end()) == (itr_sat == theSatMap.end()));
  
     if (itr_sat == theSatMap.end()) {
       vector<bool> check_sat;
       for (int i = 0; i < samples.size(); ++i) {
         if (!(samples[i] < QIE11_LINEARIZATION_ET)) {
           check_sat.push_back(true);
         } else
           check_sat.push_back(false);
       }
       theSatMap.insert(std::make_pair(id, check_sat));
     } else {
       for (int i = 0; i < samples.size(); ++i) {
         if (!(samples[i] < QIE11_LINEARIZATION_ET))
           (itr_sat->second)[i] = true;
       }
     }
   }
 }

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

◆ addSignal() [4/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE10DataFrame & frame )

private

Definition at line 209 of file HcalTriggerPrimitiveAlgo.cc.

                                                                     {
   HcalDetId detId = frame.detid();
   // prevent QIE10 calibration channels from entering TP emulation
   if (detId.subdet() != HcalForward)
     return;
  
   auto ids = theTrigTowerGeometry->towerIds(frame.id());
   for (const auto& id : ids) {
     if (id.version() == 0) {
       edm::LogError("HcalTPAlgo") << "Encountered QIE10 data frame mapped to TP version 0:" << id;
       continue;
     }
  
     int nsamples = frame.samples();
  
     IntegerCaloSamples samples(id, nsamples);
     samples.setPresamples(frame.presamples());
     incoder_->adc2Linear(frame, samples);
  
     // Don't add to final collection yet
     // HF PMT veto sum is calculated in analyzerHF()
     IntegerCaloSamples zero_samples(id, nsamples);
     zero_samples.setPresamples(frame.presamples());
     addSignal(zero_samples);
  
     auto fid = HcalDetId(frame.id());
     auto& details = theHFUpgradeDetailMap[id][fid.maskDepth()];
     auto& detail = details[fid.depth() - 1];
     detail.samples = samples;
     detail.digi = frame;
     detail.validity.resize(nsamples);
     detail.passTDC.resize(nsamples);
     incoder_->lookupMSB(frame, detail.fgbits);
     for (int idx = 0; idx < nsamples; ++idx) {
       detail.validity[idx] = validChannel(frame, idx);
       detail.passTDC[idx] = passTDC(frame, idx);
     }
   }
 }

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

◆ addSignal() [5/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE11DataFrame & frame )

private

Definition at line 249 of file HcalTriggerPrimitiveAlgo.cc.

                                                                     {
   HcalDetId detId(frame.id());
   // prevent QIE11 calibration channels from entering TP emulation
   if (detId.subdet() != HcalEndcap && detId.subdet() != HcalBarrel)
     return;
  
   std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
   assert(ids.size() == 1 || ids.size() == 2);
   IntegerCaloSamples samples1(ids[0], int(frame.samples()));
  
   samples1.setPresamples(frame.presamples());
   incoder_->adc2Linear(frame, samples1);
  
   std::vector<std::bitset<2>> msb(frame.samples(), 0);
   incoder_->lookupMSB(frame, msb);
  
   if (ids.size() == 2) {
     // make a second trigprim for the other one, and share the energy
     IntegerCaloSamples samples2(ids[1], samples1.size());
     for (int i = 0; i < samples1.size(); ++i) {
       samples1[i] = uint32_t(samples1[i]);
       samples2[i] = samples1[i];
     }
     samples2.setPresamples(frame.presamples());
     addSignal(samples2);
     addUpgradeFG(ids[1], detId.depth(), msb);
     addUpgradeTDCFG(ids[1], frame);
   }
   addSignal(samples1);
   addUpgradeFG(ids[0], detId.depth(), msb);
   addUpgradeTDCFG(ids[0], frame);
 }

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

◆ addUpgradeFG()

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

private

Definition at line 872 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                  {
   if (not validUpgradeFG(id, depth)) {
     if (needLegacyFG(id)) {
       std::vector<bool> pseudo(bits.size(), false);
       addFG(id, pseudo);
     } else if (needUpgradeID(id, depth)) {
       // If the tower id is not in the map yet
       // then for safety's sake add it, otherwise, no need
       // Likewise, we're here with non-fg depth 7 so the bits are not to be added
       auto it = fgUpgradeMap_.find(id);
       if (it == fgUpgradeMap_.end()) {
         FGUpgradeContainer element;
         element.resize(bits.size());
         fgUpgradeMap_.insert(std::make_pair(id, element));
       }
     }
  
     return;
   }
  
   auto it = fgUpgradeMap_.find(id);
   if (it == fgUpgradeMap_.end()) {
     FGUpgradeContainer element;
     element.resize(bits.size());
     it = fgUpgradeMap_.insert(std::make_pair(id, element)).first;
   }
   for (unsigned int i = 0; i < bits.size(); ++i) {
     it->second[i][0][depth - 1] = bits[i][0];
     it->second[i][1][depth - 1] = bits[i][1];
   }
 }

References addFG(), LEDCalibrationChannels::depth, fgUpgradeMap_, mps_fire::i, needLegacyFG(), needUpgradeID(), and validUpgradeFG().

Referenced by addSignal().

◆ addUpgradeTDCFG()

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

private

Definition at line 906 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                                         {
   HcalDetId detId(frame.id());
   if (detId.subdet() != HcalEndcap && detId.subdet() != HcalBarrel)
     return;
  
   std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
   assert(ids.size() == 1 || ids.size() == 2);
   IntegerCaloSamples samples1(ids[0], int(frame.samples()));
   samples1.setPresamples(frame.presamples());
   incoder_->adc2Linear(frame, samples1);  // use linearization LUT
  
   auto it = fgUpgradeTDCMap_.find(id);
   if (it == fgUpgradeTDCMap_.end()) {
     FGUpgradeTDCContainer element;
     element.resize(frame.samples());
     it = fgUpgradeTDCMap_.insert(std::make_pair(id, element)).first;
   }
   for (int i = 0; i < frame.samples(); i++) {
     it->second[i][detId.depth() - 1] = std::make_pair(samples1[i], frame[i].tdc());
   }
 }

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

Referenced by addSignal().

◆ analyze()

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

private

adds the actual digis

Definition at line 319 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                                     {
   int shrink = weights_.size() - 1;
   std::vector<bool>& msb = fgMap_[samples.id()];
   IntegerCaloSamples sum(samples.id(), samples.size());
  
   //slide algo window
   for (int ibin = 0; ibin < int(samples.size()) - shrink; ++ibin) {
     int algosumvalue = 0;
     for (unsigned int i = 0; i < weights_.size(); i++) {
       //add up value * scale factor
       algosumvalue += int(samples[ibin + i] * weights_[i]);
     }
     if (algosumvalue < 0)
       sum[ibin] = 0;  // low-side
                       //high-side
     //else if (algosumvalue>QIE8_LINEARIZATION_ET) sum[ibin]=QIE8_LINEARIZATION_ET;
     else
       sum[ibin] = algosumvalue;  //assign value to sum[]
   }
  
   // Align digis and TP
   int dgPresamples = samples.presamples();
   int tpPresamples = numberOfPresamples_;
   int shift = dgPresamples - tpPresamples;
   int dgSamples = samples.size();
   int tpSamples = numberOfSamples_;
   if (peakfind_) {
     if ((shift < shrink) || (shift + tpSamples + shrink > dgSamples - (peak_finder_algorithm_ - 1))) {
       edm::LogInfo("HcalTriggerPrimitiveAlgo::analyze")
           << "TP presample or size from the configuration file is out of the accessible range. Using digi values from "
              "data instead...";
       shift = shrink;
       tpPresamples = dgPresamples - shrink;
       tpSamples = dgSamples - (peak_finder_algorithm_ - 1) - shrink - shift;
     }
   }
  
   std::vector<int> finegrain(tpSamples, false);
  
   IntegerCaloSamples output(samples.id(), tpSamples);
   output.setPresamples(tpPresamples);
  
   for (int ibin = 0; ibin < tpSamples; ++ibin) {
     // ibin - index for output TP
     // idx - index for samples + shift
     int idx = ibin + shift;
  
     //Peak finding
     if (peakfind_) {
       bool isPeak = false;
       switch (peak_finder_algorithm_) {
         case 1:
           isPeak = (samples[idx] > samples[idx - 1] && samples[idx] >= samples[idx + 1] && samples[idx] > theThreshold);
           break;
         case 2:
           isPeak = (sum[idx] > sum[idx - 1] && sum[idx] >= sum[idx + 1] && sum[idx] > theThreshold);
           break;
         default:
           break;
       }
  
       if (isPeak) {
         output[ibin] = std::min<unsigned int>(sum[idx], QIE8_LINEARIZATION_ET);
         finegrain[ibin] = msb[idx];
       }
       // Not a peak
       else
         output[ibin] = 0;
     } else {  // No peak finding, just output running sum
       output[ibin] = std::min<unsigned int>(sum[idx], QIE8_LINEARIZATION_ET);
       finegrain[ibin] = msb[idx];
     }
  
     // Only Pegged for 1-TS algo.
     if (peak_finder_algorithm_ == 1) {
       if (samples[idx] >= QIE8_LINEARIZATION_ET)
         output[ibin] = QIE8_LINEARIZATION_ET;
     }
   }
   outcoder_->compress(output, finegrain, result);
 }

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().

◆ analyzeHF()

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

private

Definition at line 508 of file HcalTriggerPrimitiveAlgo.cc.

                                                                   {
   HcalTrigTowerDetId detId(samples.id());
  
   // Align digis and TP
   int dgPresamples = samples.presamples();
   int tpPresamples = numberOfPresamplesHF_;
   int shift = dgPresamples - tpPresamples;
   int dgSamples = samples.size();
   int tpSamples = numberOfSamplesHF_;
   if (shift < 0 || shift + tpSamples > dgSamples) {
     edm::LogInfo("HcalTriggerPrimitiveAlgo::analyzeHF")
         << "TP presample or size from the configuration file is out of the accessible range. Using digi values from "
            "data instead...";
     tpPresamples = dgPresamples;
     shift = 0;
     tpSamples = dgSamples;
   }
  
   std::vector<int> finegrain(tpSamples, false);
  
   TowerMapFGSum::const_iterator tower2fg = theTowerMapFGSum.find(detId);
   assert(tower2fg != theTowerMapFGSum.end());
  
   const SumFGContainer& sumFG = tower2fg->second;
   // Loop over all L+S pairs that mapped from samples.id()
   // Note: 1 samples.id() = 6 x (L+S) without noZS
   for (SumFGContainer::const_iterator sumFGItr = sumFG.begin(); sumFGItr != sumFG.end(); ++sumFGItr) {
     const std::vector<bool>& veto = HF_Veto[sumFGItr->id().rawId()];
     for (int ibin = 0; ibin < tpSamples; ++ibin) {
       int idx = ibin + shift;
       // if not vetod, add L+S to total sum and calculate FG
       bool vetoed = idx < int(veto.size()) && veto[idx];
       if (!(vetoed && (*sumFGItr)[idx] > PMT_NoiseThreshold_)) {
         samples[idx] += (*sumFGItr)[idx];
         finegrain[ibin] = (finegrain[ibin] || (*sumFGItr)[idx] >= FG_threshold_);
       }
     }
   }
  
   IntegerCaloSamples output(samples.id(), tpSamples);
   output.setPresamples(tpPresamples);
  
   for (int ibin = 0; ibin < tpSamples; ++ibin) {
     int idx = ibin + shift;
     output[ibin] = samples[idx] >> hf_lumi_shift;
     static const int MAX_OUTPUT = QIE8_LINEARIZATION_ET;  // QIE8_LINEARIZATION_ET = 1023
     if (output[ibin] > MAX_OUTPUT)
       output[ibin] = MAX_OUTPUT;
   }
   outcoder_->compress(output, finegrain, result);
 }

References cms::cuda::assert(), HcalTPGCompressor::compress(), 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().

◆ analyzeHF2016()

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

private

Definition at line 562 of file HcalTriggerPrimitiveAlgo.cc.

                                                                           {
   // Align digis and TP
   const int SHIFT = samples.presamples() - numberOfPresamplesHF_;
   assert(SHIFT >= 0);
   assert((SHIFT + numberOfSamplesHF_) <= samples.size());
  
   // Try to find the HFDetails from the map corresponding to our samples
   const HcalTrigTowerDetId detId(samples.id());
   HFDetailMap::const_iterator it = theHFDetailMap.find(detId);
   // Missing values will give an empty digi
   if (it == theHFDetailMap.end()) {
     return;
   }
  
   std::vector<std::bitset<2>> finegrain(numberOfSamplesHF_, false);
  
   // Set up out output of IntergerCaloSamples
   IntegerCaloSamples output(samples.id(), numberOfSamplesHF_);
   output.setPresamples(numberOfPresamplesHF_);
  
   for (const auto& item : it->second) {
     auto& details = item.second;
     for (int ibin = 0; ibin < numberOfSamplesHF_; ++ibin) {
       const int IDX = ibin + SHIFT;
       int long_fiber_val = 0;
       if (IDX < details.long_fiber.size()) {
         long_fiber_val = details.long_fiber[IDX];
       }
       int short_fiber_val = 0;
       if (IDX < details.short_fiber.size()) {
         short_fiber_val = details.short_fiber[IDX];
       }
       output[ibin] += (long_fiber_val + short_fiber_val);
  
       uint32_t ADCLong = details.LongDigi[ibin].adc();
       uint32_t ADCShort = details.ShortDigi[ibin].adc();
  
       if (details.LongDigi.id().ietaAbs() >= FIRST_FINEGRAIN_TOWER) {
         finegrain[ibin][1] = (ADCLong > FG_HF_thresholds_[0] || ADCShort > FG_HF_thresholds_[0]);
  
         if (embit != nullptr)
           finegrain[ibin][0] = embit->fineGrainbit(details.ShortDigi, details.LongDigi, ibin);
       }
     }
   }
  
   for (int bin = 0; bin < numberOfSamplesHF_; ++bin) {
     static const unsigned int MAX_OUTPUT = QIE8_LINEARIZATION_ET;  // QIE8_LINEARIZATION_ET = 1023
     output[bin] = min({MAX_OUTPUT, output[bin] >> hf_lumi_shift});
   }
  
   std::vector<int> finegrain_converted;
   finegrain_converted.reserve(finegrain.size());
   for (const auto& fg : finegrain)
     finegrain_converted.push_back(fg.to_ulong());
   outcoder_->compress(output, finegrain_converted, result);
 }

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

Referenced by run().

◆ analyzeHFQIE10()

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

private

Definition at line 651 of file HcalTriggerPrimitiveAlgo.cc.

                                                                            {
   // Align digis and TP
   const int shift = samples.presamples() - numberOfPresamplesHF_;
   assert(shift >= 0);
   assert((shift + numberOfSamplesHF_) <= samples.size());
   assert(hf_lumi_shift >= 2);
  
   // Try to find the HFDetails from the map corresponding to our samples
   const HcalTrigTowerDetId detId(samples.id());
   auto it = theHFUpgradeDetailMap.find(detId);
   // Missing values will give an empty digi
   if (it == theHFUpgradeDetailMap.end()) {
     return;
   }
  
   std::vector<std::bitset<2>> finegrain(numberOfSamplesHF_, false);
  
   // Set up out output of IntergerCaloSamples
   IntegerCaloSamples output(samples.id(), numberOfSamplesHF_);
   output.setPresamples(numberOfPresamplesHF_);
  
   for (const auto& item : it->second) {
     auto& details = item.second;
     for (int ibin = 0; ibin < numberOfSamplesHF_; ++ibin) {
       const int idx = ibin + shift;
  
       int long_fiber_val = 0;
       int long_fiber_count = 0;
       int short_fiber_val = 0;
       int short_fiber_count = 0;
  
       bool saturated = false;
  
       for (auto i : {0, 2}) {
         if (idx < details[i].samples.size() and details[i].validity[idx] and details[i].passTDC[idx]) {
           long_fiber_val += details[i].samples[idx];
           saturated = saturated || (details[i].samples[idx] == QIE10_LINEARIZATION_ET);
           ++long_fiber_count;
         }
       }
       for (auto i : {1, 3}) {
         if (idx < details[i].samples.size() and details[i].validity[idx] and details[i].passTDC[idx]) {
           short_fiber_val += details[i].samples[idx];
           saturated = saturated || (details[i].samples[idx] == QIE10_LINEARIZATION_ET);
           ++short_fiber_count;
         }
       }
  
       if (saturated) {
         output[ibin] = QIE10_MAX_LINEARIZATION_ET;
       } else {
         // For details of the energy handling, see:
         // https://cms-docdb.cern.ch/cgi-bin/DocDB/ShowDocument?docid=12306
         // If both readouts are valid, average of the two energies is taken
         // division by 2 is compensated by adjusting the total scale shift in the end
         if (long_fiber_count == 2)
           long_fiber_val >>= 1;
         if (short_fiber_count == 2)
           short_fiber_val >>= 1;
  
         auto sum = long_fiber_val + short_fiber_val;
         // Similar to above, if both channels are valid,
         // average of the two energies is calculated
         // division by 2 here is also compensated by adjusting the total scale shift in the end
         if (long_fiber_count > 0 and short_fiber_count > 0)
           sum >>= 1;
  
         output[ibin] += sum;
       }
  
       for (const auto& detail : details) {
         if (idx < int(detail.digi.size()) and detail.validity[idx] and
             HcalDetId(detail.digi.id()).ietaAbs() >= FIRST_FINEGRAIN_TOWER) {
           if (useTDCInMinBiasBits_ && !detail.passTDC[idx])
             continue;
           finegrain[ibin][1] = finegrain[ibin][1] or detail.fgbits[idx][0];
           // what is commonly called the "second" HF min-bias bit is
           // actually the 0-th bit, which can also be used instead for the EM bit
           // (called finegrain[ibin][0] below) in non-HI running
           finegrain[ibin][0] = finegrain[ibin][0] or detail.fgbits[idx][1];
         }
       }
       // the EM bit is only used if the "second" FG bit is disabled
       if (embit != nullptr and FG_HF_thresholds_.at(1) != 255) {
         finegrain[ibin][0] = embit->fineGrainbit(details[1].digi,
                                                  details[3].digi,
                                                  details[0].digi,
                                                  details[2].digi,
                                                  details[1].validity[idx],
                                                  details[3].validity[idx],
                                                  details[0].validity[idx],
                                                  details[2].validity[idx],
                                                  idx);
       }
     }
   }
  
   for (int bin = 0; bin < numberOfSamplesHF_; ++bin) {
     output[bin] = min({(unsigned int)QIE10_MAX_LINEARIZATION_ET, output[bin] >> (hf_lumi_shift - 2)});
   }
   std::vector<int> finegrain_converted;
   finegrain_converted.reserve(finegrain.size());
   for (const auto& fg : finegrain)
     finegrain_converted.push_back(fg.to_ulong());
   outcoder_->compress(output, finegrain_converted, result);
 }

References cms::cuda::assert(), newFWLiteAna::bin, HcalTPGCompressor::compress(), FG_HF_thresholds_, HcalFeatureBit::fineGrainbit(), FIRST_FINEGRAIN_TOWER, mps_fire::i, heavyIonCSV_trainingSettings::idx, createfilelist::int, B2GTnPMonitor_cfi::item, 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().

◆ 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.

                                                                              {
   HcalDetId detId(samples.id());
  
   // Get the |ieta| for current sample
   int theIeta = detId.ietaAbs();
  
   unsigned int dgSamples = samples.size();
   unsigned int dgPresamples = samples.presamples();
  
   unsigned int tpSamples = numberOfSamples_;
   unsigned int tpPresamples = numberOfPresamples_;
  
   unsigned int filterSamples = weightsQIE11_[theIeta].size();
   unsigned int filterPresamples = theIeta > theTrigTowerGeometry->topology().lastHBRing()
                                       ? numberOfFilterPresamplesHEQIE11_
                                       : numberOfFilterPresamplesHBQIE11_;
  
   unsigned int shift = dgPresamples - tpPresamples;
  
   // shrink keeps the FIR filter from going off the end of the 8TS vector
   unsigned int shrink = filterSamples - 1;
  
   auto& msb = fgUpgradeMap_[samples.id()];
   auto& timingTDC = fgUpgradeTDCMap_[samples.id()];
   IntegerCaloSamples sum(samples.id(), samples.size());
  
   std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
  
   // keep track of tower with saturated energy and force the total TP saturated
   bool force_saturation[samples.size()];
   for (int i = 0; i < samples.size(); i++) {
     force_saturation[i] = false;
   }
  
   //slide algo window
   for (unsigned int ibin = 0; ibin < dgSamples - shrink; ++ibin) {
     int algosumvalue = 0;
     bool check_sat = false;
     for (unsigned int i = 0; i < filterSamples; i++) {
       //add up value * scale factor
       // In addition, divide by two in the 10 degree phi segmentation region
       // to mimic 5 degree segmentation for the trigger
       unsigned int sample = samples[ibin + i];
  
       if (fix_saturation_ && (sample_saturation.size() > ibin + i))
         check_sat = (sample_saturation[ibin + i] | (sample > QIE11_MAX_LINEARIZATION_ET));
       else if (sample > QIE11_MAX_LINEARIZATION_ET)
         sample = QIE11_MAX_LINEARIZATION_ET;
  
       // Usually use a segmentation factor of 1.0 but for ieta >= 21 use 0.5
       double segmentationFactor = 1.0;
       if (ids.size() == 2) {
         segmentationFactor = 0.5;
       }
  
       // Based on the |ieta| of the sample, retrieve the correct region weight
       double theWeight = weightsQIE11_[theIeta][i];
  
       algosumvalue += int(sample * segmentationFactor * theWeight);
     }
     if (algosumvalue < 0)
       sum[ibin] = 0;  // low-side
                       //high-side
     //else if (algosumvalue>QIE11_LINEARIZATION_ET) sum[ibin]=QIE11_LINEARIZATION_ET;
     else
       sum[ibin] = algosumvalue;  //assign value to sum[]
  
     if (check_sat)
       force_saturation[ibin] = true;
   }
  
   std::vector<int> finegrain(tpSamples, false);
  
   IntegerCaloSamples output(samples.id(), tpSamples);
   output.setPresamples(tpPresamples);
  
   for (unsigned int ibin = 0; ibin < tpSamples; ++ibin) {
     // ibin - index for output TP
     // idx - index for samples + shift - filterPresamples
     int idx = ibin + shift - filterPresamples;
  
     // When idx is <= 0 peakfind would compare out-of-bounds of the vector. Avoid this ambiguity
     if (idx <= 0) {
       output[ibin] = 0;
       continue;
     }
  
     bool isPeak = (sum[idx] > sum[idx - 1] && sum[idx] >= sum[idx + 1] && sum[idx] > theThreshold);
  
     if (isPeak) {
       output[ibin] = std::min<unsigned int>(sum[idx], QIE11_MAX_LINEARIZATION_ET);
       if (fix_saturation_ && force_saturation[idx])
         output[ibin] = QIE11_MAX_LINEARIZATION_ET;
     } else {
       // Not a peak
       output[ibin] = 0;
     }
     // peak-finding is not applied for FG bits
     // compute(msb) returns two bits (MIP). compute(timingTDC,ids) returns 6 bits (1 depth, 2 reserved, 1 prompt, 1 delayed 01, 1 delayed 10)
     finegrain[ibin] = fg_algo.compute(timingTDC[idx], ids[0]).to_ulong() | fg_algo.compute(msb[idx]).to_ulong() << 1;
   }
   outcoder_->compress(output, finegrain, result);
 }

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

Referenced by run().

◆ needLegacyFG()

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

private

Definition at line 853 of file HcalTriggerPrimitiveAlgo.cc.

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

References HBHE_OVERLAP_TOWER, and upgrade_hb_.

Referenced by addUpgradeFG().

◆ needUpgradeID()

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

private

Definition at line 862 of file HcalTriggerPrimitiveAlgo.cc.

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

References FIRST_DEPTH7_TOWER, and LAST_FINEGRAIN_DEPTH.

Referenced by addUpgradeFG().

◆ overrideParameters()

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

Definition at line 81 of file HcalTriggerPrimitiveAlgo.cc.

                                                                            {
   override_parameters_ = ps;
  
   if (override_parameters_.exists("ADCThresholdHF")) {
     override_adc_hf_ = true;
     override_adc_hf_value_ = override_parameters_.getParameter<uint32_t>("ADCThresholdHF");
   }
   if (override_parameters_.exists("TDCMaskHF")) {
     override_tdc_hf_ = true;
     override_tdc_hf_value_ = override_parameters_.getParameter<unsigned long long>("TDCMaskHF");
   }
 }

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().

◆ passTDC()

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

private

Definition at line 623 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                {
   auto parameters = conditions_->getHcalTPParameters();
   auto adc_threshold = parameters->getADCThresholdHF();
   auto tdc_mask = parameters->getTDCMaskHF();
  
   if (override_adc_hf_)
     adc_threshold = override_adc_hf_value_;
   if (override_tdc_hf_)
     tdc_mask = override_tdc_hf_value_;
  
   if (digi[ts].adc() < adc_threshold)
     return true;
  
   return (1ul << digi[ts].le_tdc()) & tdc_mask;
 }

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

Referenced by addSignal().

◆ 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.

                                                           {
   theTrigTowerGeometry = trigTowerGeometry;
  
   incoder_ = dynamic_cast<const HcaluLUTTPGCoder*>(incoder);
   outcoder_ = outcoder;
   conditions_ = conditions;
  
   theSumMap.clear();
   theSatMap.clear();
   theTowerMapFGSum.clear();
   HF_Veto.clear();
   fgMap_.clear();
   fgUpgradeMap_.clear();
   fgUpgradeTDCMap_.clear();
   theHFDetailMap.clear();
   theHFUpgradeDetailMap.clear();
  
   // Add all digi collections
   addDigis(digis...);
  
   // Prepare the fine-grain calculation algorithm for HB/HE
   int version = 0;
   if (upgrade_he_ or upgrade_hb_)
     version = conditions_->getHcalTPParameters()->getFGVersionHBHE();
   if (override_parameters_.exists("FGVersionHBHE"))
     version = override_parameters_.getParameter<uint32_t>("FGVersionHBHE");
   HcalFinegrainBit fg_algo(version);
  
   // VME produces additional bits on the front used by lumi but not the
   // trigger, this shift corrects those out by right shifting over them.
   for (auto& item : theSumMap) {
     result.push_back(HcalTriggerPrimitiveDigi(item.first));
     HcalTrigTowerDetId detId(item.second.id());
     if (detId.ietaAbs() >= theTrigTowerGeometry->firstHFTower(detId.version())) {
       if (detId.version() == 0) {
         analyzeHF(item.second, result.back(), RCTScaleShift);
       } else if (detId.version() == 1) {
         if (upgrade_hf_)
           analyzeHFQIE10(item.second, result.back(), NCTScaleShift, LongvrsShortCut);
         else
           analyzeHF2016(item.second, result.back(), NCTScaleShift, LongvrsShortCut);
       } else {
         // Things are going to go poorly
       }
     } else {
       // Determine which energy reconstruction path to take based on the
       // fine-grain availability:
       //  * QIE8 TP add entries into fgMap_
       //  * QIE11 TP add entries into fgUpgradeMap_
       //    (not for tower 16 unless HB is upgraded, too)
       if (fgMap_.find(item.first) != fgMap_.end()) {
         analyze(item.second, result.back());
       } else if (fgUpgradeMap_.find(item.first) != fgUpgradeMap_.end()) {
         SatMap::iterator item_sat = theSatMap.find(detId);
         if (item_sat == theSatMap.end())
           analyzeQIE11(item.second, std::vector<bool>(), result.back(), fg_algo);
         else
           analyzeQIE11(item.second, item_sat->second, result.back(), fg_algo);
       }
     }
   }
  
   // Free up some memory
   theSumMap.clear();
   theSatMap.clear();
   theTowerMapFGSum.clear();
   HF_Veto.clear();
   fgMap_.clear();
   fgUpgradeMap_.clear();
   fgUpgradeTDCMap_.clear();
   theHFDetailMap.clear();
   theHFUpgradeDetailMap.clear();
  
   return;
 }

References addDigis(), analyze(), analyzeHF(), analyzeHF2016(), analyzeHFQIE10(), analyzeQIE11(), submitPVValidationJobs::conditions, conditions_, 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().

◆ runFEFormatError()

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

Definition at line 777 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                     {
   std::set<uint32_t> FrontEndErrors;
  
   for (int i = FEDNumbering::MINHCALFEDID; i <= FEDNumbering::MAXHCALFEDID; ++i) {
     const FEDRawData& raw = rawraw->FEDData(i);
     if (raw.size() < 12)
       continue;
     const HcalDCCHeader* dccHeader = (const HcalDCCHeader*)(raw.data());
     if (!dccHeader)
       continue;
     HcalHTRData htr;
     for (int spigot = 0; spigot < HcalDCCHeader::SPIGOT_COUNT; spigot++) {
       if (!dccHeader->getSpigotPresent(spigot))
         continue;
       dccHeader->getSpigotData(spigot, htr, raw.size());
       int dccid = dccHeader->getSourceId();
       int errWord = htr.getErrorsWord() & 0x1FFFF;
       bool HTRError = (!htr.check() || htr.isHistogramEvent() || (errWord & 0x800) != 0);
  
       if (HTRError) {
         bool valid = false;
         for (int fchan = 0; fchan < 3 && !valid; fchan++) {
           for (int fib = 0; fib < 9 && !valid; fib++) {
             HcalElectronicsId eid(fchan, fib, spigot, dccid - FEDNumbering::MINHCALFEDID);
             eid.setHTR(htr.readoutVMECrateId(), htr.htrSlot(), htr.htrTopBottom());
             DetId detId = emap->lookup(eid);
             if (detId.null())
               continue;
             HcalSubdetector subdet = (HcalSubdetector(detId.subdetId()));
             if (detId.det() != 4 || (subdet != HcalBarrel && subdet != HcalEndcap && subdet != HcalForward))
               continue;
             std::vector<HcalTrigTowerDetId> ids = theTrigTowerGeometry->towerIds(detId);
             for (std::vector<HcalTrigTowerDetId>::const_iterator triggerId = ids.begin(); triggerId != ids.end();
                  ++triggerId) {
               FrontEndErrors.insert(triggerId->rawId());
             }
             //valid = true;
           }
         }
       }
     }
   }
  
   // Loop over TP collection
   // Set TP to zero if there is FE Format Error
   HcalTriggerPrimitiveSample zeroSample(0);
   for (HcalTrigPrimDigiCollection::iterator tp = result.begin(); tp != result.end(); ++tp) {
     if (FrontEndErrors.find(tp->id().rawId()) != FrontEndErrors.end()) {
       for (int i = 0; i < tp->size(); ++i)
         tp->setSample(i, zeroSample);
     }
   }
 }

References HcalHTRData::check(), FEDRawData::data(), DetId::det(), 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, DetId::null(), HcalHTRData::readoutVMECrateId(), mps_fire::result, FEDRawData::size(), HcalDCCHeader::SPIGOT_COUNT, DetId::subdetId(), theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), cmsswSequenceInfo::tp, and RunInfoPI::valid.

Referenced by HcalTrigPrimDigiProducer::produce().

◆ runZS()

void HcalTriggerPrimitiveAlgo::runZS ( HcalTrigPrimDigiCollection & tp )

Definition at line 761 of file HcalTriggerPrimitiveAlgo.cc.

                                                                        {
   for (HcalTrigPrimDigiCollection::iterator tp = result.begin(); tp != result.end(); ++tp) {
     bool ZS = true;
     for (int i = 0; i < tp->size(); ++i) {
       if (tp->sample(i).compressedEt() > ZS_threshold_I_) {
         ZS = false;
         break;
       }
     }
     if (ZS)
       tp->setZSInfo(false, true);
     else
       tp->setZSInfo(true, false);
   }
 }

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

Referenced by HcalTrigPrimDigiProducer::produce().

◆ setFixSaturationFlag()

void HcalTriggerPrimitiveAlgo::setFixSaturationFlag ( bool fix_saturation )

Definition at line 79 of file HcalTriggerPrimitiveAlgo.cc.

79 { fix_saturation_ = fix_saturation; }

References fix_saturation_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

◆ setNCTScaleShift()

void HcalTriggerPrimitiveAlgo::setNCTScaleShift ( int shift )

Definition at line 950 of file HcalTriggerPrimitiveAlgo.cc.

950 { NCTScaleShift = shift; }

References NCTScaleShift, and edm::shift.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

◆ setNumFilterPresamplesHBQIE11()

void HcalTriggerPrimitiveAlgo::setNumFilterPresamplesHBQIE11 ( int presamples )

inline

Definition at line 87 of file HcalTriggerPrimitiveAlgo.h.

87 { numberOfFilterPresamplesHBQIE11_ = presamples; }

References numberOfFilterPresamplesHBQIE11_, and hcalTTPDigis_cfi::presamples.

Referenced by HcalTrigPrimDigiProducer::beginRun().

◆ setNumFilterPresamplesHEQIE11()

void HcalTriggerPrimitiveAlgo::setNumFilterPresamplesHEQIE11 ( int presamples )

inline

Definition at line 89 of file HcalTriggerPrimitiveAlgo.h.

89 { numberOfFilterPresamplesHEQIE11_ = presamples; }

References numberOfFilterPresamplesHEQIE11_, and hcalTTPDigis_cfi::presamples.

Referenced by HcalTrigPrimDigiProducer::beginRun().

◆ setPeakFinderAlgorithm()

void HcalTriggerPrimitiveAlgo::setPeakFinderAlgorithm ( int algo )

Definition at line 944 of file HcalTriggerPrimitiveAlgo.cc.

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

References Exception, and peak_finder_algorithm_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

◆ setRCTScaleShift()

void HcalTriggerPrimitiveAlgo::setRCTScaleShift ( int shift )

Definition at line 952 of file HcalTriggerPrimitiveAlgo.cc.

952 { RCTScaleShift = shift; }

References RCTScaleShift, and edm::shift.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

◆ setUpgradeFlags()

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

Definition at line 73 of file HcalTriggerPrimitiveAlgo.cc.

                                                                         {
   upgrade_hb_ = hb;
   upgrade_he_ = he;
   upgrade_hf_ = hf;
 }

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

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

◆ setWeightQIE11()

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

Definition at line 938 of file HcalTriggerPrimitiveAlgo.cc.

                                                                       {
   // Simple map of |ieta| in HBHE to weight
   // Only one weight for SOI-1 TS
   weightsQIE11_[aieta] = {{weight, 1.0}};
 }

References mps_merge::weight, and weightsQIE11_.

Referenced by HcalTrigPrimDigiProducer::beginRun().

◆ setWeightsQIE11()

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

Definition at line 928 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                   {
   // Names are just abs(ieta) for HBHE
   std::vector<std::string> ietaStrs = weightsQIE11.getParameterNames();
   for (auto& ietaStr : ietaStrs) {
     // Strip off "ieta" part of key and just use integer value in map
     auto const& v = weightsQIE11.getParameter<std::vector<double>>(ietaStr);
     weightsQIE11_[std::stoi(ietaStr.substr(4))] = {{v[0], v[1]}};
   }
 }

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

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

◆ validChannel()

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

private

Definition at line 639 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                     {
   // channels with invalid data should not contribute to the sum
   if (digi.linkError() || ts >= digi.samples() || !digi[ts].ok())
     return false;
  
   auto mask = conditions_->getHcalTPChannelParameter(HcalDetId(digi.id()))->getMask();
   if (mask)
     return false;
  
   return true;
 }

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

Referenced by addSignal().

◆ validUpgradeFG()

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

private

Definition at line 843 of file HcalTriggerPrimitiveAlgo.cc.

                                                                                            {
   if (depth > LAST_FINEGRAIN_DEPTH)
     return false;
   if (id.ietaAbs() > LAST_FINEGRAIN_TOWER)
     return false;
   if (id.ietaAbs() == HBHE_OVERLAP_TOWER and not upgrade_hb_)
     return false;
   return true;
 }

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

Referenced by addUpgradeFG().

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<double, 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().

Classes

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ FGbitMap

◆ FGUpgradeContainer

◆ FGUpgradeMap

◆ FGUpgradeTDCContainer

◆ FGUpgradeTDCMap

◆ HFDetailMap

◆ HFUpgradeDetailMap

◆ SatMap

◆ SumFGContainer

◆ SumMap

◆ TowerMapFGSum

◆ TowerMapVeto

Constructor & Destructor Documentation

◆ HcalTriggerPrimitiveAlgo()

◆ ~HcalTriggerPrimitiveAlgo()

Member Function Documentation

◆ addDigis() [1/3]

◆ addDigis() [2/3]

◆ addDigis() [3/3]

◆ addFG()

◆ addSignal() [1/5]

◆ addSignal() [2/5]

◆ addSignal() [3/5]

◆ addSignal() [4/5]

◆ addSignal() [5/5]

◆ addUpgradeFG()

◆ addUpgradeTDCFG()

◆ analyze()

◆ analyzeHF()

◆ analyzeHF2016()

◆ analyzeHFQIE10()

◆ analyzeQIE11()

◆ needLegacyFG()

◆ needUpgradeID()

◆ overrideParameters()

◆ passTDC()

◆ run()

◆ runFEFormatError()

◆ runZS()

◆ setFixSaturationFlag()

◆ setNCTScaleShift()

◆ setNumFilterPresamplesHBQIE11()

◆ setNumFilterPresamplesHEQIE11()

◆ setPeakFinderAlgorithm()

◆ setRCTScaleShift()

◆ setUpgradeFlags()

◆ setWeightQIE11()

◆ setWeightsQIE11()

◆ validChannel()

◆ validUpgradeFG()

Member Data Documentation

◆ conditions_

◆ FG_HF_thresholds_

◆ FG_threshold_

◆ fgMap_

◆ fgUpgradeMap_

◆ fgUpgradeTDCMap_

◆ FIRST_DEPTH7_TOWER

◆ FIRST_FINEGRAIN_TOWER

◆ fix_saturation_

◆ HBHE_OVERLAP_TOWER

◆ HF_Veto

◆ incoder_

◆ LAST_FINEGRAIN_DEPTH

◆ LAST_FINEGRAIN_TOWER

◆ latency_

◆ LongvrsShortCut

◆ minSignalThreshold_

◆ NCTScaleShift

◆ numberOfFilterPresamplesHBQIE11_

◆ numberOfFilterPresamplesHEQIE11_

◆ numberOfPresamples_