#include <HcalTTPDigiProducer.h>

Inheritance diagram for HcalTTPDigiProducer:

Public Member Functions
	HcalTTPDigiProducer (const edm::ParameterSet &ps)

void	produce (edm::Event &e, const edm::EventSetup &c) override

	~HcalTTPDigiProducer () override=default

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
bool	decision (int nP, int nM, int bit)

bool	isMasked (HcalDetId id)

Private Attributes
std::string	bit_ [4]

int	calc_ [4]

int	fwAlgo_

int	id_

int	iEtaMax_

int	iEtaMin_

std::vector< unsigned int >	maskedChannels_

char	mReq_ [4]

int	nHFm_ [4]

int	nHFp_ [4]

int	nHits_ [4]

char	pmLogic_ [4]

char	pReq_ [4]

int	presamples_

int	samples_

int	SoI_

unsigned int	threshold_

edm::EDGetTokenT< HFDigiCollection >	tok_hf_

edm::ESGetToken< HcalTPGCoder, HcalTPGRecord >	tok_tpgCoder_

Static Private Attributes
static const int	inputs_ []

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 14 of file HcalTTPDigiProducer.h.

Constructor & Destructor Documentation

◆ HcalTTPDigiProducer()

HcalTTPDigiProducer::HcalTTPDigiProducer ( const edm::ParameterSet & ps )

explicit

Definition at line 16 of file HcalTTPDigiProducer.cc.

                                                                   {
   tok_hf_ = consumes<HFDigiCollection>(ps.getParameter<edm::InputTag>("HFDigiCollection"));
   tok_tpgCoder_ = esConsumes<HcalTPGCoder, HcalTPGRecord>();
   maskedChannels_ = ps.getParameter<std::vector<unsigned int> >("maskedChannels");
   bit_[0] = ps.getParameter<std::string>("defTT8");
   bit_[1] = ps.getParameter<std::string>("defTT9");
   bit_[2] = ps.getParameter<std::string>("defTT10");
   bit_[3] = ps.getParameter<std::string>("defTTLocal");
  
   for (int i = 0; i < 4; i++) {
     nHits_[i] = -1;
     nHFp_[i] = -1;
     nHFm_[i] = -1;
     pReq_[i] = ' ';
     mReq_[i] = ' ';
     pmLogic_[i] = ' ';
     calc_[i] = sscanf(bit_[i].c_str(),
                       "hits>=%d:hfp%c=%d%chfm%c=%d",
                       &(nHits_[i]),
                       &(pReq_[i]),
                       &(nHFp_[i]),
                       &(pmLogic_[i]),
                       &(mReq_[i]),
                       &(nHFm_[i]));
     if (calc_[i] == 1) {
       if (nHits_[i] < 0)
         throw cms::Exception("HcalTTPDigiProducer") << "Unable to read logic for technical trigger";
     } else if (calc_[i] == 6) {
       if (nHits_[i] < 0 || nHFp_[i] < 0 || nHFm_[i] < 0)
         throw cms::Exception("HcalTTPDigiProducer") << "Unable to read logic for technical trigger";
       if ((pReq_[i] != '>' && pReq_[i] != '<') || (mReq_[i] != '>' && mReq_[i] != '<') ||
           (pmLogic_[i] != ':' && pmLogic_[i] != '|'))
         throw cms::Exception("HcalTTPDigiProducer") << "Technical Trigger logic must obey the following format:\n"
                                                        "\"hits>=[A1]:hfp[B1]=[A2][C]hfm[B2]=[A3]\",\n"
                                                        "or \"hits>=[A1]\",\n"
                                                        "with A# >= 0, B# = (</>) and C = (:/|)";
     } else {
       throw cms::Exception("HcalTTPDigiProducer") << "Unable to read logic for technical trigger";
     }
   }
  
   id_ = ps.getUntrackedParameter<int>("id", -1);
   samples_ = ps.getParameter<int>("samples");
   presamples_ = ps.getParameter<int>("presamples");
   iEtaMin_ = ps.getParameter<int>("iEtaMin");
   iEtaMax_ = ps.getParameter<int>("iEtaMax");
   threshold_ = ps.getParameter<unsigned int>("threshold");
   fwAlgo_ = ps.getParameter<int>("fwAlgorithm");
  
   SoI_ = ps.getParameter<int>("HFSoI");
  
   if (samples_ > 8) {
     samples_ = 8;
     edm::LogWarning("HcalTTPDigiProducer") << "Samples forced to maximum value of 8";
   }
   if (presamples_ - SoI_ > 0) {  // Too many presamples
     presamples_ = SoI_;
     edm::LogWarning("HcalTTPDigiProducer") << "Presamples reset to HF SoI value";
   }
  
   produces<HcalTTPDigiCollection>();
 }

References bit_, calc_, Exception, fwAlgo_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, id_, iEtaMax_, iEtaMin_, maskedChannels_, mReq_, nHFm_, nHFp_, nHits_, pmLogic_, pReq_, presamples_, samples_, SoI_, AlCaHLTBitMon_QueryRunRegistry::string, threshold_, tok_hf_, and tok_tpgCoder_.

◆ ~HcalTTPDigiProducer()

HcalTTPDigiProducer::~HcalTTPDigiProducer ( )

overridedefault

Member Function Documentation

◆ decision()

bool HcalTTPDigiProducer::decision	(	int	nP,
		int	nM,
		int	bit
	)

private

Definition at line 86 of file HcalTTPDigiProducer.cc.

                                                           {
   bool pOK = false;
   bool mOK = false;
   if ((nP + nM) < nHits_[bit])
     return false;
   if (calc_[bit] == 1)
     return ((nP + nM) >= nHits_[bit]);
  
   if (pReq_[bit] == '>')
     pOK = (nP >= nHFp_[bit]);
   else if (pReq_[bit] == '<')
     pOK = (nP <= nHFp_[bit]);
  
   if (mReq_[bit] == '>')
     mOK = (nM >= nHFm_[bit]);
   else if (mReq_[bit] == '<')
     mOK = (nM <= nHFm_[bit]);
  
   if (pmLogic_[bit] == ':')
     return (pOK && mOK);
   else if (pmLogic_[bit] == '|')
     return (pOK || mOK);
  
   // Should not ever get here...need to create a warning message
   edm::LogWarning("HcalTTPDigiProducer") << "Trigger logic exhausted.  Returning false";
   return false;
 }

References calc_, mReq_, nHFm_, nHFp_, nHits_, pmLogic_, and pReq_.

Referenced by produce().

◆ isMasked()

bool HcalTTPDigiProducer::isMasked ( HcalDetId id )

private

Definition at line 79 of file HcalTTPDigiProducer.cc.

                                                {
   for (unsigned int i = 0; i < maskedChannels_.size(); i++)
     if (id.rawId() == maskedChannels_.at(i))
       return true;
   return false;
 }

References mps_fire::i, and maskedChannels_.

Referenced by produce().

◆ produce()

void HcalTTPDigiProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

override

Definition at line 114 of file HcalTTPDigiProducer.cc.

                                                                               {
   // Step A: Get Inputs
   edm::Handle<HFDigiCollection> hfDigiCollection;
   e.getByToken(tok_hf_, hfDigiCollection);
   edm::ESHandle<HcalTPGCoder> inputCoder = eventSetup.getHandle(tok_tpgCoder_);
  
   // Step B: Create empty output
   std::unique_ptr<HcalTTPDigiCollection> ttpResult(new HcalTTPDigiCollection());
  
   // Step C: Compute TTP inputs
   uint16_t trigInputs[40];
   int nP[8];
   int nM[8];
   for (int i = 0; i < 8; i++) {
     nP[i] = 0;
     nM[i] = 0;
     for (int j = 0; j < 5; j++)
       trigInputs[j * 8 + i] = 0;
   }
   for (HFDigiCollection::const_iterator theDigi = hfDigiCollection->begin(); theDigi != hfDigiCollection->end();
        theDigi++) {
     HcalDetId id = HcalDetId(theDigi->id());
     if (isMasked(id))
       continue;
     if (id.ietaAbs() < iEtaMin_ || id.ietaAbs() > iEtaMax_)
       continue;
  
     IntegerCaloSamples samples(id, theDigi->size());
     inputCoder->adc2Linear(*theDigi, samples);
  
     for (int relSample = -presamples_; relSample < (samples_ - presamples_); relSample++) {
       if (samples[SoI_ + relSample] >= threshold_) {
         int linSample = presamples_ + relSample;
         int offset = (-1 + id.zside()) / 2;
         int shift = inputs_[id.iphi() + offset];
         int group = 0;
         while (shift >= 16) {
           shift -= 16;
           group++;
         }
         if (!(trigInputs[(linSample * 8) + group] & (1 << shift)))
           (id.ieta() > 0) ? (nP[linSample]++) : (nM[linSample]++);
         trigInputs[(linSample * 8) + group] |= (1 << shift);
       }
     }
   }
  
   // Step D: Compute trigger decision and fill TTP digi
   uint8_t trigOutput[8];
   uint32_t algoDepBits[8];
   HcalTTPDigi ttpDigi(id_, samples_, presamples_, 0, fwAlgo_, 0);
   for (int linSample = 0; linSample < 8; linSample++) {
     trigOutput[linSample] = 0;
     algoDepBits[linSample] = 0;
     if (linSample < samples_) {
       for (int j = 0; j < 4; j++)
         trigOutput[linSample] |= (decision(nP[linSample], nM[linSample], j) << j);
       int nT = nP[linSample] + nM[linSample];
  
       // Algorithm Dependent bits for FW flavor = 1
       // NOTE: this disagrees with the fw var. names that implies (LSB) T,M,P (MSB)
       if (fwAlgo_ == 1)
         algoDepBits[linSample] = (nT & 0x7F) | ((nP[linSample] & 0x3F) << 7) | ((nM[linSample] & 0x3F) << 13);
       ttpDigi.setSample(
           (linSample - presamples_), &trigInputs[linSample * 8], algoDepBits[linSample], trigOutput[linSample]);
     }
   }
   ttpResult->push_back(ttpDigi);
  
   // Step E: Put outputs into event
   e.put(std::move(ttpResult));
 }

References HcalTPGCoder::adc2Linear(), edm::SortedCollection< T, SORT >::begin(), decision(), MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), fwAlgo_, edm::EventSetup::getHandle(), watchdog::group, mps_fire::i, id_, iEtaMax_, iEtaMin_, inputs_, isMasked(), dqmiolumiharvest::j, eostools::move(), hltrates_dqm_sourceclient-live_cfg::offset, presamples_, EgammaValidation_cff::samples, samples_, HcalTTPDigi::setSample(), edm::shift, SoI_, threshold_, tok_hf_, and tok_tpgCoder_.

Member Data Documentation

◆ bit_

std::string HcalTTPDigiProducer::bit_[4]

private

Definition at line 28 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer().

◆ calc_

int HcalTTPDigiProducer::calc_[4]

private

Definition at line 29 of file HcalTTPDigiProducer.h.

Referenced by decision(), and HcalTTPDigiProducer().

◆ fwAlgo_

int HcalTTPDigiProducer::fwAlgo_

private

Definition at line 33 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

◆ id_

int HcalTTPDigiProducer::id_

private

Definition at line 32 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

◆ iEtaMax_

int HcalTTPDigiProducer::iEtaMax_

private

Definition at line 34 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

◆ iEtaMin_

int HcalTTPDigiProducer::iEtaMin_

private

Definition at line 34 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

◆ inputs_

const int HcalTTPDigiProducer::inputs_

staticprivate

Initial value:

= {30, 66, 4,  44, 4,  44, 0,  68, 0,  68, 16, 48, 16, 48, 6,  46, 6,  46,
                                            2,  70, 2,  70, 18, 50, 18, 50, 12, 40, 12, 40, 8,  52, 8,  52, 20, 36,
                                            20, 36, 14, 42, 14, 42, 10, 54, 10, 54, 22, 38, 22, 38, 24, 56, 24, 56,
                                            32, 60, 32, 60, 28, 64, 28, 64, 26, 58, 26, 58, 34, 62, 34, 62, 30, 66}

Definition at line 39 of file HcalTTPDigiProducer.h.

Referenced by produce().