#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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_

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 11 of file HcalTTPDigiProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 23 of file HcalTTPDigiProducer.cc.

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_, and tok_hf_.

 {
     tok_hf_ = consumes<HFDigiCollection>(ps.getParameter<edm::InputTag>("HFDigiCollection")); 
     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>();
 }

HcalTTPDigiProducer::~HcalTTPDigiProducer ( )

overridedefault

Member Function Documentation

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

private

Definition at line 89 of file HcalTTPDigiProducer.cc.

References calc_, funct::false, mReq_, nHFm_, nHFp_, nHits_, pmLogic_, and pReq_.

Referenced by produce().

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

bool HcalTTPDigiProducer::isMasked ( HcalDetId id )

private

Definition at line 82 of file HcalTTPDigiProducer.cc.

References funct::false, mps_fire::i, and maskedChannels_.

Referenced by produce().

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

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

override

Definition at line 109 of file HcalTTPDigiProducer.cc.

References HcalTPGCoder::adc2Linear(), edm::SortedCollection< T, SORT >::begin(), decision(), edm::SortedCollection< T, SORT >::end(), fwAlgo_, edm::EventSetup::get(), edm::Event::getByToken(), edmIntegrityCheck::group, mps_fire::i, id_, iEtaMax_, iEtaMin_, inputs_, isMasked(), eostools::move(), PFRecoTauDiscriminationByIsolation_cfi::offset, presamples_, edm::Event::put(), hcalTTPDigis_cfi::samples, samples_, HcalTTPDigi::setSample(), edm::shift, SoI_, threshold_, and tok_hf_.

                                                                               {
     
     // Step A: Get Inputs
     edm::Handle<HFDigiCollection> hfDigiCollection ; 
     e.getByToken(tok_hf_,hfDigiCollection) ;
     edm::ESHandle<HcalTPGCoder> inputCoder ;
     eventSetup.get<HcalTPGRecord>().get(inputCoder) ;
 
     // 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));
 }

Member Data Documentation

std::string HcalTTPDigiProducer::bit_[4]

private

Definition at line 27 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer().

int HcalTTPDigiProducer::calc_[4]

private

Definition at line 28 of file HcalTTPDigiProducer.h.

Referenced by decision(), and HcalTTPDigiProducer().

int HcalTTPDigiProducer::fwAlgo_

private

Definition at line 32 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

int HcalTTPDigiProducer::id_

private

Definition at line 31 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

int HcalTTPDigiProducer::iEtaMax_

private

Definition at line 33 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

int HcalTTPDigiProducer::iEtaMin_

private

Definition at line 33 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and produce().

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 38 of file HcalTTPDigiProducer.h.

Referenced by produce().

std::vector<unsigned int> HcalTTPDigiProducer::maskedChannels_

private

Definition at line 26 of file HcalTTPDigiProducer.h.

Referenced by HcalTTPDigiProducer(), and isMasked().

char HcalTTPDigiProducer::mReq_[4]

private

Definition at line 30 of file HcalTTPDigiProducer.h.

Referenced by decision(), and HcalTTPDigiProducer().

int HcalTTPDigiProducer::nHFm_[4]

private