#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, int numberOfSamplesZDC, int numberOfPresamplesZDC, 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	setCodedVetoThreshold (int aieta, int codedVetoThreshold)

void	setCodedVetoThresholds (const edm::ParameterSet &codedVetoThresholds)

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)

void	analyzeZDC (IntegerCaloSamples &samples, HcalTriggerPrimitiveDigi &result)

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
std::array< int, 29 >	codedVetoThresholds_

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	numberOfPresamplesZDC_

int	numberOfSamples_

int	numberOfSamplesHF_

int	numberOfSamplesZDC_

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	QIE10_ZDC_MAX_LINEARIZATION_ET = 0x3FF

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 220 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeContainer

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

private

Definition at line 223 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeMap

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

private

Definition at line 224 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeTDCContainer

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

private

Definition at line 227 of file HcalTriggerPrimitiveAlgo.h.

◆ FGUpgradeTDCMap

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

private

Definition at line 228 of file HcalTriggerPrimitiveAlgo.h.

◆ HFDetailMap

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

private

Definition at line 190 of file HcalTriggerPrimitiveAlgo.h.

◆ HFUpgradeDetailMap

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

private

Definition at line 200 of file HcalTriggerPrimitiveAlgo.h.

◆ SatMap

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

private

Definition at line 181 of file HcalTriggerPrimitiveAlgo.h.

◆ SumFGContainer

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

private

Definition at line 203 of file HcalTriggerPrimitiveAlgo.h.

◆ SumMap

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

private

Definition at line 178 of file HcalTriggerPrimitiveAlgo.h.

◆ TowerMapFGSum

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

private

Definition at line 204 of file HcalTriggerPrimitiveAlgo.h.

◆ TowerMapVeto

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

private

Definition at line 217 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,
		int	numberOfSamplesZDC,
		int	numberOfPresamplesZDC,
		bool	useTDCInMinBiasBits,
		uint32_t	minSignalThreshold = `0`,
		uint32_t	PMT_NoiseThreshold = `0`
	)

Definition at line 23 of file HcalTriggerPrimitiveAlgo.cc.

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

     : 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),
       numberOfSamplesZDC_(numberOfSamplesZDC),
       numberOfPresamplesZDC_(numberOfPresamplesZDC),
       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;
     numberOfSamplesZDC_ = 1;
     numberOfPresamplesZDC_ = 0;
   }
   // Switch to integer for comparisons - remove compiler warning
   ZS_threshold_I_ = ZS_threshold_;
 }

◆ ~HcalTriggerPrimitiveAlgo()

HcalTriggerPrimitiveAlgo::~HcalTriggerPrimitiveAlgo ( )

Definition at line 77 of file HcalTriggerPrimitiveAlgo.cc.

77 {}

Member Function Documentation

◆ addDigis() [1/3]

template<typename T , typename... Args>

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

inline

Definition at line 59 of file HcalTriggerPrimitiveAlgo.h.

References universalConfigTemplate::collection.

Referenced by run().

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

◆ addDigis() [2/3]

template<typename T >

void HcalTriggerPrimitiveAlgo::addDigis ( const T & collection )

inline

Definition at line 65 of file HcalTriggerPrimitiveAlgo.h.

References addSignal(), and universalConfigTemplate::collection.

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

◆ addDigis() [3/3]

template<typename D >

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

inline

Definition at line 72 of file HcalTriggerPrimitiveAlgo.h.

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

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

◆ addFG()

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

private

Definition at line 980 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by addSignal(), and addUpgradeFG().

                                                                                        {
   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;
 }

◆ addSignal() [1/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const HBHEDataFrame & frame )

private

adds the signal to the map

Definition at line 100 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().

                                                                    {
   // 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);
 }

◆ addSignal() [2/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const HFDataFrame & frame )

private

Definition at line 131 of file HcalTriggerPrimitiveAlgo.cc.

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

                                                                  {
   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;
       }
     }
   }
 }

◆ addSignal() [3/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE10DataFrame & frame )

private

Definition at line 215 of file HcalTriggerPrimitiveAlgo.cc.

References HcaluLUTTPGCoder::adc2Linear(), addSignal(), DetId::Calo, hcalRecHitTable_cff::detId, HcalZDCDetId::EM, amptDefault_cfi::frame, HcalZDCDetId::HAD, DetId::Hcal, HcalForward, mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, heavyIonCSV_trainingSettings::idx, incoder_, HcaluLUTTPGCoder::lookupMSB(), passTDC(), EgammaValidation_cff::samples, IntegerCaloSamples::setPresamples(), HcalZDCDetId::SubdetectorId, theHFUpgradeDetailMap, theTrigTowerGeometry, HcalTrigTowerGeometry::towerIds(), HcalTrigTowerGeometry::towerIds_ZDC(), validChannel(), and relval_steps::version.

                                                                     {
   //HcalDetId detId = frame.detid();
   DetId detId = DetId(frame.detid());
   if (detId.det() == DetId::Hcal) {
     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, false);
 
       // 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);
       }
     }
   } else if (detId.det() == DetId::Calo && detId.subdetId() == HcalZDCDetId::SubdetectorId) {
     HcalZDCDetId detId = frame.detid();
     if (detId.section() != HcalZDCDetId::EM && detId.section() != HcalZDCDetId::HAD) {
       return;
     }
 
     auto ids = theTrigTowerGeometry->towerIds_ZDC(frame.id());
     for (const auto& id : ids) {
       int nsamples = frame.samples();
 
       IntegerCaloSamples samples(id, nsamples);
       IntegerCaloSamples samples_PUsub(id, nsamples);
 
       samples.setPresamples(frame.presamples());
       samples_PUsub.setPresamples(frame.presamples());
 
       incoder_->adc2Linear(frame, samples, false);
       incoder_->adc2Linear(frame, samples_PUsub, true);
 
       for (int i = 1; i < samples.size(); ++i) {
         if (samples_PUsub[i - 1] > samples[i])
           samples[i] = 0;
         else
           samples[i] -= samples_PUsub[i - 1];
       }
 
       addSignal(samples);
     }
   }
 }

◆ addSignal() [4/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const QIE11DataFrame & frame )

private

Definition at line 286 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().

                                                                     {
   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);
 }

◆ addSignal() [5/5]

void HcalTriggerPrimitiveAlgo::addSignal ( const IntegerCaloSamples & samples )

private

Definition at line 319 of file HcalTriggerPrimitiveAlgo.cc.

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

                                                                           {
   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;
       }
     }
   }
 }

◆ addUpgradeFG()

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

private

Definition at line 1019 of file HcalTriggerPrimitiveAlgo.cc.

References addFG(), ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits, hcalRecHitTable_cff::depth, fgUpgradeMap_, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, needLegacyFG(), needUpgradeID(), and validUpgradeFG().

Referenced by addSignal().

                                                                                  {
   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];
   }
 }

◆ addUpgradeTDCFG()

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

private

Definition at line 1053 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by addSignal().

                                                                                                         {
   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
   std::vector<unsigned short> bits12_15 = incoder_->group0FGbits(frame);  // get 4 energy bits (12-15) from group 0 LUT
 
   bool is_compressed = false;
   if (detId.subdet() == HcalBarrel) {
     is_compressed = (frame.flavor() == 3);
     // 0 if frame.flavor is 0 (uncompressed), 1 if frame.flavor is 3 (compressed)
   }
 
   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(std::make_pair(bits12_15[i], is_compressed), std::make_pair(frame[i].tdc(), samples1[i]));
   }
 }

◆ analyze()

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

private

adds the actual digis

Definition at line 356 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().

                                                                                                     {
   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);
 }

◆ analyzeHF()

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

private

Definition at line 655 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 relval_cleanedupgrade::veto.

Referenced by run().

                                                                   {
   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);
 }

◆ analyzeHF2016()

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

private

Definition at line 709 of file HcalTriggerPrimitiveAlgo.cc.

References cms::cuda::assert(), newFWLiteAna::bin, HcalTPGCompressor::compress(), hcalRecHitTable_cff::detId, FG_HF_thresholds_, HcalFeatureBit::fineGrainbit(), FIRST_FINEGRAIN_TOWER, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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().

                                                                           {
   // 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);
 }

◆ analyzeHFQIE10()

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

private

Definition at line 798 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, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), createfilelist::int, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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().

                                                                            {
   // 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);
 }

◆ analyzeQIE11()

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

private

Definition at line 438 of file HcalTriggerPrimitiveAlgo.cc.

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

                                                                              {
   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;
     //TP energy calculation for PFA2
     if (weightsQIE11_[theIeta][0] == 255) {
       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 = (check_sat | sample_saturation[ibin + i] | (sample > QIE11_MAX_LINEARIZATION_ET));
 
         if (sample > QIE11_MAX_LINEARIZATION_ET)
           sample = QIE11_MAX_LINEARIZATION_ET;
 
         // Usually use a segmentation factor of 1.0 but for ieta >= 21 use 2
         int segmentationFactor = 1;
         if (ids.size() == 2) {
           segmentationFactor = 2;
         }
 
         algosumvalue += int(sample / segmentationFactor);
       }
       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;
       //TP energy calculation for PFA1' and PFA1
     } else {
       //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
       int sampleTS = samples[ibin + 1];
       int sampleTSminus1 = samples[ibin];
 
       if (fix_saturation_ && (sample_saturation.size() > ibin + 1))
         check_sat |= sample_saturation[ibin + 1] | (sampleTS >= QIE11_MAX_LINEARIZATION_ET);
 
       if (sampleTS > QIE11_MAX_LINEARIZATION_ET)
         sampleTS = QIE11_MAX_LINEARIZATION_ET;
 
       if (sampleTSminus1 > QIE11_MAX_LINEARIZATION_ET || sample_saturation[ibin])
         sampleTSminus1 = QIE11_MAX_LINEARIZATION_ET;
 
       // Usually use a segmentation factor of 1.0 but for ieta >= 21 use 2
       int segmentationFactor = 1;
       if (ids.size() == 2) {
         segmentationFactor = 2;
       }
 
       // Based on the |ieta| of the sample, retrieve the correct region weight
       int theWeight = weightsQIE11_[theIeta][0];
 
       algosumvalue = ((sampleTS << 8) - (sampleTSminus1 * theWeight)) / 256 / segmentationFactor;
 
       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);
 
   // Based on the |ieta| of the sample, retrieve the correct region "coded" veto threshold
   // where two of the possible values have special meaning
   unsigned int codedVetoThreshold = codedVetoThresholds_[theIeta];
 
   // Anything in range (1, 2048) inclusive shall activate the veto
   unsigned int actualVetoThreshold = codedVetoThreshold;
   bool applyVetoThreshold = codedVetoThreshold > 0 && codedVetoThreshold <= 2048;
 
   // Special value to disable vetoing in the PFA1' algo is 0
   if (codedVetoThreshold > 0) {
     if (codedVetoThreshold <= 2048) {
       // Special value to run the veto in PFA1' with no threshold
       if (codedVetoThreshold == 2048)
         actualVetoThreshold = 0;
     } else {
       edm::LogWarning("HcalTPAlgo") << "Specified veto threshold value " << codedVetoThreshold
                                     << " is not in range (1, 2048) ! Vetoing in PFA1' will not be enabled !";
     }
   }
 
   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;
     }
 
     //Only run the peak-finder when the PFA2 FIR filter is running, which corresponds to weights = 1
     if (weightsQIE11_[theIeta][0] == 255) {
       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] && ids.size() == 2)
           output[ibin] = QIE11_MAX_LINEARIZATION_ET / 2;
         else if (fix_saturation_ && force_saturation[idx])
           output[ibin] = QIE11_MAX_LINEARIZATION_ET;
 
       } else {
         // Not a peak
         output[ibin] = 0;
       }
     } else {
       // Only if the sum for the future time sample is above the veto
       // threshold and the now sum is not a peak and the now sum is not
       // saturated does the current sum get zeroed
       if (applyVetoThreshold && sum[idx + 1] >= actualVetoThreshold &&
           (sum[idx] < sum[idx + 1] || force_saturation[idx + 1]) && !force_saturation[idx])
         output[ibin] = 0;
       else {
         // Here, either the "now" sum is a peak or the vetoing criteria are not satisfied
         // so assign the appropriate sum to the output
         output[ibin] = std::min<unsigned int>(sum[idx], QIE11_MAX_LINEARIZATION_ET);
         if (fix_saturation_ && force_saturation[idx]) {
           output[ibin] = QIE11_MAX_LINEARIZATION_ET;
           if (ids.size() == 2)
             output[ibin] /= 2;
         }
       }
     }
     // peak-finding is not applied for FG bits
     // compute(msb) returns two bits (MIP). compute(timingTDC,ids) returns 6 bits (1 depth, 1 prompt, 1 delayed 01, 1 delayed 10, 2 reserved)
     finegrain[ibin] = fg_algo.compute(timingTDC[idx + filterPresamples], ids[0]).to_ulong() |
                       fg_algo.compute(msb[idx + filterPresamples]).to_ulong() << 4;
     if (ibin == tpPresamples && (idx + filterPresamples) != dgPresamples)
       edm::LogError("HcalTriggerPritimveAlgo")
           << "TP SOI (tpPresamples = " << tpPresamples
           << ") is not aligned with digi SOI (dgPresamples = " << dgPresamples << ")";
   }
   outcoder_->compress(output, finegrain, result);
 }

◆ analyzeZDC()

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

private

Definition at line 631 of file HcalTriggerPrimitiveAlgo.cc.

References hcalRecHitTable_cff::detId, mps_fire::i, convertSQLitetoXML_cfg::output, QIE10_ZDC_MAX_LINEARIZATION_ET, mps_fire::result, and EgammaValidation_cff::samples.

Referenced by run().

                                                                                                        {
   HcalTrigTowerDetId detId(samples.id());
 
   unsigned int tpSamples;
   unsigned int tpPresamples;
 
   tpSamples = samples.size();
   tpPresamples = samples.presamples();
   result.setSize(tpSamples);
   result.setPresamples(tpPresamples);
 
   IntegerCaloSamples output(samples.id(), tpSamples);
   output.setPresamples(tpPresamples);
 
   for (int i = 0; i < samples.size(); i++) {
     if (samples[i] > QIE10_ZDC_MAX_LINEARIZATION_ET)
       output[i] = QIE10_ZDC_MAX_LINEARIZATION_ET;
     else
       output[i] = samples[i];
     HcalTriggerPrimitiveSample zdcSample(output[i]);
     result.setSample(i, zdcSample);
   }
 }

◆ needLegacyFG()

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

private

Definition at line 1000 of file HcalTriggerPrimitiveAlgo.cc.

References HBHE_OVERLAP_TOWER, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), and upgrade_hb_.

Referenced by addUpgradeFG().

                                                                               {
   // 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;
 }

◆ needUpgradeID()

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

private

Definition at line 1009 of file HcalTriggerPrimitiveAlgo.cc.

References FIRST_DEPTH7_TOWER, and LAST_FINEGRAIN_DEPTH.

Referenced by addUpgradeFG().

                                                                                           {
   // 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;
 }

◆ overrideParameters()

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

Definition at line 87 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().

                                                                            {
   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");
   }
 }

◆ passTDC()

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

private

Definition at line 770 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().

                                                                                {
   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;
 }

◆ 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 263 of file HcalTriggerPrimitiveAlgo.h.

References addDigis(), analyze(), analyzeHF(), analyzeHF2016(), analyzeHFQIE10(), analyzeQIE11(), analyzeZDC(), PixelBaryCentreAnalyzer_cfg::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 relval_steps::version.

Referenced by HcalTrigPrimDigiProducer::produce().

                                                           {
   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()) && detId.ietaAbs() < 42) {
       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 if (detId.ietaAbs() >= theTrigTowerGeometry->firstHFTower(detId.version()) && detId.ietaAbs() > 41) {
       analyzeZDC(item.second, result.back());
     } 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;
 }

◆ runFEFormatError()

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

Definition at line 924 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().

                                                                                     {
   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);
     }
   }
 }

◆ runZS()

void HcalTriggerPrimitiveAlgo::runZS ( HcalTrigPrimDigiCollection & tp )

Definition at line 908 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by HcalTrigPrimDigiProducer::produce().

                                                                        {
   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);
   }
 }

◆ setCodedVetoThreshold()

void HcalTriggerPrimitiveAlgo::setCodedVetoThreshold	(	int	aieta,
		int	codedVetoThreshold
	)

Definition at line 1109 of file HcalTriggerPrimitiveAlgo.cc.

References codedVetoThresholds_.

Referenced by HcalTrigPrimDigiProducer::beginRun().

                                                                                       {
   // Simple map of |ieta| in HBHE to veto threshold
   codedVetoThresholds_[aieta] = {codedVetoThreshold};
 }

◆ setCodedVetoThresholds()

void HcalTriggerPrimitiveAlgo::setCodedVetoThresholds ( const edm::ParameterSet & codedVetoThresholds )

Definition at line 1099 of file HcalTriggerPrimitiveAlgo.cc.

References hcaltpdigi_cfi::codedVetoThresholds, codedVetoThresholds_, and findQualityFiles::v.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

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

◆ setFixSaturationFlag()

void HcalTriggerPrimitiveAlgo::setFixSaturationFlag ( bool fix_saturation )

Definition at line 85 of file HcalTriggerPrimitiveAlgo.cc.

References fix_saturation_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

85 { fix_saturation_ = fix_saturation; }

HcalTriggerPrimitiveAlgo::fix_saturation_

bool fix_saturation_

Definition: HcalTriggerPrimitiveAlgo.h:235

◆ setNCTScaleShift()

void HcalTriggerPrimitiveAlgo::setNCTScaleShift ( int shift )

Definition at line 1120 of file HcalTriggerPrimitiveAlgo.cc.

References NCTScaleShift, and edm::shift.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

1120 { NCTScaleShift = shift; }

HcalTriggerPrimitiveAlgo::NCTScaleShift

int NCTScaleShift

Definition: HcalTriggerPrimitiveAlgo.h:164

edm::shift

static unsigned int const shift

Definition: LuminosityBlockID.cc:7

◆ setNumFilterPresamplesHBQIE11()

void HcalTriggerPrimitiveAlgo::setNumFilterPresamplesHBQIE11 ( int presamples )

inline

Definition at line 91 of file HcalTriggerPrimitiveAlgo.h.

References numberOfFilterPresamplesHBQIE11_, and hcalTTPDigis_cfi::presamples.

Referenced by HcalTrigPrimDigiProducer::beginRun().

91 { numberOfFilterPresamplesHBQIE11_ = presamples; }

hcalTTPDigis_cfi.presamples

presamples

Definition: hcalTTPDigis_cfi.py:9

HcalTriggerPrimitiveAlgo::numberOfFilterPresamplesHBQIE11_

int numberOfFilterPresamplesHBQIE11_

Definition: HcalTriggerPrimitiveAlgo.h:155

◆ setNumFilterPresamplesHEQIE11()

void HcalTriggerPrimitiveAlgo::setNumFilterPresamplesHEQIE11 ( int presamples )

inline

Definition at line 93 of file HcalTriggerPrimitiveAlgo.h.

References numberOfFilterPresamplesHEQIE11_, and hcalTTPDigis_cfi::presamples.

Referenced by HcalTrigPrimDigiProducer::beginRun().

93 { numberOfFilterPresamplesHEQIE11_ = presamples; }

hcalTTPDigis_cfi.presamples

presamples

Definition: hcalTTPDigis_cfi.py:9

HcalTriggerPrimitiveAlgo::numberOfFilterPresamplesHEQIE11_

int numberOfFilterPresamplesHEQIE11_

Definition: HcalTriggerPrimitiveAlgo.h:156

◆ setPeakFinderAlgorithm()

void HcalTriggerPrimitiveAlgo::setPeakFinderAlgorithm ( int algo )

Definition at line 1114 of file HcalTriggerPrimitiveAlgo.cc.

References Exception, and peak_finder_algorithm_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

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

◆ setRCTScaleShift()

void HcalTriggerPrimitiveAlgo::setRCTScaleShift ( int shift )

Definition at line 1122 of file HcalTriggerPrimitiveAlgo.cc.

References RCTScaleShift, and edm::shift.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

1122 { RCTScaleShift = shift; }

HcalTriggerPrimitiveAlgo::RCTScaleShift

int RCTScaleShift

Definition: HcalTriggerPrimitiveAlgo.h:165

edm::shift

static unsigned int const shift

Definition: LuminosityBlockID.cc:7

◆ setUpgradeFlags()

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

Definition at line 79 of file HcalTriggerPrimitiveAlgo.cc.

References Hcal_Conditions_forGlobalTag_cff::hb, Hcal_Conditions_forGlobalTag_cff::he, Hcal_Conditions_forGlobalTag_cff::hf, upgrade_hb_, upgrade_he_, and upgrade_hf_.

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

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

◆ setWeightQIE11()

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

Definition at line 1093 of file HcalTriggerPrimitiveAlgo.cc.

References mps_merge::weight, and weightsQIE11_.

Referenced by HcalTrigPrimDigiProducer::beginRun().

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

◆ setWeightsQIE11()

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

Definition at line 1083 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by HcalTrigPrimDigiProducer::HcalTrigPrimDigiProducer().

                                                                                   {
   // 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<int>>(ietaStr);
     weightsQIE11_[std::stoi(ietaStr.substr(4))] = {{v[0], v[1]}};
   }
 }

◆ validChannel()

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

private

Definition at line 786 of file HcalTriggerPrimitiveAlgo.cc.

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

Referenced by addSignal().

                                                                                     {
   // 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;
 }

◆ validUpgradeFG()

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

private

Definition at line 990 of file HcalTriggerPrimitiveAlgo.cc.

References hcalRecHitTable_cff::depth, HBHE_OVERLAP_TOWER, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), LAST_FINEGRAIN_DEPTH, LAST_FINEGRAIN_TOWER, and upgrade_hb_.

Referenced by addUpgradeFG().

                                                                                            {
   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;
 }

Member Data Documentation

◆ codedVetoThresholds_

std::array<int, 29> HcalTriggerPrimitiveAlgo::codedVetoThresholds_

private

Definition at line 147 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11(), setCodedVetoThreshold(), and setCodedVetoThresholds().

◆ conditions_

const HcalDbService* HcalTriggerPrimitiveAlgo::conditions_

private

Definition at line 142 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 150 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2016(), and analyzeHFQIE10().

◆ FG_threshold_

uint32_t HcalTriggerPrimitiveAlgo::FG_threshold_

private

Definition at line 149 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF().

◆ fgMap_

FGbitMap HcalTriggerPrimitiveAlgo::fgMap_

private

Definition at line 221 of file HcalTriggerPrimitiveAlgo.h.

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

◆ fgUpgradeMap_

FGUpgradeMap HcalTriggerPrimitiveAlgo::fgUpgradeMap_

private

Definition at line 225 of file HcalTriggerPrimitiveAlgo.h.

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

◆ fgUpgradeTDCMap_

FGUpgradeTDCMap HcalTriggerPrimitiveAlgo::fgUpgradeTDCMap_

private

Definition at line 229 of file HcalTriggerPrimitiveAlgo.h.

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

◆ FIRST_DEPTH7_TOWER

const int HcalTriggerPrimitiveAlgo::FIRST_DEPTH7_TOWER = 26

staticprivate

Definition at line 246 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needUpgradeID().

◆ FIRST_FINEGRAIN_TOWER

const int HcalTriggerPrimitiveAlgo::FIRST_FINEGRAIN_TOWER = 30

staticprivate

Definition at line 251 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF2016(), and analyzeHFQIE10().

◆ fix_saturation_

bool HcalTriggerPrimitiveAlgo::fix_saturation_ = false

private

Definition at line 235 of file HcalTriggerPrimitiveAlgo.h.

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

◆ HBHE_OVERLAP_TOWER

const int HcalTriggerPrimitiveAlgo::HBHE_OVERLAP_TOWER = 16

staticprivate

Definition at line 245 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needLegacyFG(), and validUpgradeFG().

◆ HF_Veto

TowerMapVeto HcalTriggerPrimitiveAlgo::HF_Veto

private

Definition at line 218 of file HcalTriggerPrimitiveAlgo.h.

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

◆ incoder_

const HcaluLUTTPGCoder* HcalTriggerPrimitiveAlgo::incoder_

private

Definition at line 140 of file HcalTriggerPrimitiveAlgo.h.

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

◆ LAST_FINEGRAIN_DEPTH

const int HcalTriggerPrimitiveAlgo::LAST_FINEGRAIN_DEPTH = 6

staticprivate

Definition at line 247 of file HcalTriggerPrimitiveAlgo.h.

Referenced by needUpgradeID(), and validUpgradeFG().

◆ LAST_FINEGRAIN_TOWER

const int HcalTriggerPrimitiveAlgo::LAST_FINEGRAIN_TOWER = 28

staticprivate

Definition at line 248 of file HcalTriggerPrimitiveAlgo.h.

Referenced by validUpgradeFG().

◆ latency_

int HcalTriggerPrimitiveAlgo::latency_

private

Definition at line 148 of file HcalTriggerPrimitiveAlgo.h.

◆ LongvrsShortCut

HcalFeatureBit* HcalTriggerPrimitiveAlgo::LongvrsShortCut

private

Definition at line 216 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run().

◆ minSignalThreshold_

uint32_t HcalTriggerPrimitiveAlgo::minSignalThreshold_

private

Definition at line 162 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal().

◆ NCTScaleShift

int HcalTriggerPrimitiveAlgo::NCTScaleShift

private

Definition at line 164 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setNCTScaleShift().

◆ numberOfFilterPresamplesHBQIE11_

int HcalTriggerPrimitiveAlgo::numberOfFilterPresamplesHBQIE11_

private

Definition at line 155 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11(), and setNumFilterPresamplesHBQIE11().

◆ numberOfFilterPresamplesHEQIE11_

int HcalTriggerPrimitiveAlgo::numberOfFilterPresamplesHEQIE11_

private

Definition at line 156 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11(), and setNumFilterPresamplesHEQIE11().

◆ numberOfPresamples_

int HcalTriggerPrimitiveAlgo::numberOfPresamples_

private

Definition at line 154 of file HcalTriggerPrimitiveAlgo.h.

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

◆ numberOfPresamplesHF_

int HcalTriggerPrimitiveAlgo::numberOfPresamplesHF_

private

Definition at line 158 of file HcalTriggerPrimitiveAlgo.h.

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

◆ numberOfPresamplesZDC_

int HcalTriggerPrimitiveAlgo::numberOfPresamplesZDC_

private

Definition at line 160 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo().

◆ numberOfSamples_

int HcalTriggerPrimitiveAlgo::numberOfSamples_

private

Definition at line 153 of file HcalTriggerPrimitiveAlgo.h.

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

◆ numberOfSamplesHF_

int HcalTriggerPrimitiveAlgo::numberOfSamplesHF_

private

Definition at line 157 of file HcalTriggerPrimitiveAlgo.h.

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

◆ numberOfSamplesZDC_

int HcalTriggerPrimitiveAlgo::numberOfSamplesZDC_

private

Definition at line 159 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo().

◆ outcoder_

const HcalTPGCompressor* HcalTriggerPrimitiveAlgo::outcoder_

private

Definition at line 141 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 239 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ override_adc_hf_value_

uint32_t HcalTriggerPrimitiveAlgo::override_adc_hf_value_

private

Definition at line 240 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ override_parameters_

edm::ParameterSet HcalTriggerPrimitiveAlgo::override_parameters_

private

Definition at line 237 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and run().

◆ override_tdc_hf_

bool HcalTriggerPrimitiveAlgo::override_tdc_hf_ = false

private

Definition at line 241 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 242 of file HcalTriggerPrimitiveAlgo.h.

Referenced by overrideParameters(), and passTDC().

◆ peak_finder_algorithm_

int HcalTriggerPrimitiveAlgo::peak_finder_algorithm_

private

Definition at line 171 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and setPeakFinderAlgorithm().

◆ peakfind_

bool HcalTriggerPrimitiveAlgo::peakfind_

private

Definition at line 144 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and HcalTriggerPrimitiveAlgo().

◆ PMT_NoiseThreshold_

uint32_t HcalTriggerPrimitiveAlgo::PMT_NoiseThreshold_

private

Definition at line 163 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHF().

◆ QIE10_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE10_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE10_LUT_BITMASK

staticprivate

Definition at line 254 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHFQIE10().

◆ QIE10_MAX_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE10_MAX_LINEARIZATION_ET = 0x7FF

staticprivate

Definition at line 257 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHFQIE10().

◆ QIE10_ZDC_MAX_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE10_ZDC_MAX_LINEARIZATION_ET = 0x3FF

staticprivate

Definition at line 259 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeZDC().

◆ QIE11_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE11_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE11_LUT_BITMASK

staticprivate

Definition at line 255 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal().

◆ QIE11_MAX_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE11_MAX_LINEARIZATION_ET = 0x7FF

staticprivate

Definition at line 258 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeQIE11().

◆ QIE8_LINEARIZATION_ET

const int HcalTriggerPrimitiveAlgo::QIE8_LINEARIZATION_ET = HcaluLUTTPGCoder::QIE8_LUT_BITMASK

staticprivate

Definition at line 253 of file HcalTriggerPrimitiveAlgo.h.

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

◆ RCTScaleShift

int HcalTriggerPrimitiveAlgo::RCTScaleShift

private

Definition at line 165 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setRCTScaleShift().

◆ theHFDetailMap

HFDetailMap HcalTriggerPrimitiveAlgo::theHFDetailMap

private

Definition at line 191 of file HcalTriggerPrimitiveAlgo.h.

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

◆ theHFUpgradeDetailMap

HFUpgradeDetailMap HcalTriggerPrimitiveAlgo::theHFUpgradeDetailMap

private

Definition at line 201 of file HcalTriggerPrimitiveAlgo.h.

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

◆ theSatMap

SatMap HcalTriggerPrimitiveAlgo::theSatMap

private

Definition at line 182 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), and run().

◆ theSumMap

SumMap HcalTriggerPrimitiveAlgo::theSumMap

private

Definition at line 179 of file HcalTriggerPrimitiveAlgo.h.

Referenced by addSignal(), and run().

◆ theThreshold

double HcalTriggerPrimitiveAlgo::theThreshold

private

Definition at line 143 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze(), and analyzeQIE11().

◆ theTowerMapFGSum

TowerMapFGSum HcalTriggerPrimitiveAlgo::theTowerMapFGSum

private

Definition at line 205 of file HcalTriggerPrimitiveAlgo.h.

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

◆ theTrigTowerGeometry

const HcalTrigTowerGeometry* HcalTriggerPrimitiveAlgo::theTrigTowerGeometry

private

Definition at line 176 of file HcalTriggerPrimitiveAlgo.h.

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

◆ upgrade_hb_

bool HcalTriggerPrimitiveAlgo::upgrade_hb_ = false

private

Definition at line 231 of file HcalTriggerPrimitiveAlgo.h.

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

◆ upgrade_he_

bool HcalTriggerPrimitiveAlgo::upgrade_he_ = false

private

Definition at line 232 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

◆ upgrade_hf_

bool HcalTriggerPrimitiveAlgo::upgrade_hf_ = false

private

Definition at line 233 of file HcalTriggerPrimitiveAlgo.h.

Referenced by run(), and setUpgradeFlags().

◆ useTDCInMinBiasBits_

bool HcalTriggerPrimitiveAlgo::useTDCInMinBiasBits_

private

Definition at line 161 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyzeHFQIE10().

◆ weights_

std::vector<double> HcalTriggerPrimitiveAlgo::weights_

private

Definition at line 145 of file HcalTriggerPrimitiveAlgo.h.

Referenced by analyze().

◆ weightsQIE11_

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

private

Definition at line 146 of file HcalTriggerPrimitiveAlgo.h.

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

◆ ZS_threshold_

uint32_t HcalTriggerPrimitiveAlgo::ZS_threshold_

private

Definition at line 151 of file HcalTriggerPrimitiveAlgo.h.

Referenced by HcalTriggerPrimitiveAlgo().

◆ ZS_threshold_I_

int HcalTriggerPrimitiveAlgo::ZS_threshold_I_

private

Definition at line 152 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()

◆ analyzeZDC()

◆ needLegacyFG()

◆ needUpgradeID()

◆ overrideParameters()

◆ passTDC()

◆ run()

◆ runFEFormatError()

◆ runZS()

◆ setCodedVetoThreshold()

◆ setCodedVetoThresholds()

◆ setFixSaturationFlag()

◆ setNCTScaleShift()

◆ setNumFilterPresamplesHBQIE11()

◆ setNumFilterPresamplesHEQIE11()

◆ setPeakFinderAlgorithm()

◆ setRCTScaleShift()

◆ setUpgradeFlags()

◆ setWeightQIE11()

◆ setWeightsQIE11()

◆ validChannel()

◆ validUpgradeFG()

Member Data Documentation

◆ codedVetoThresholds_

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