#include <CastorTTRecord.h>

Inheritance diagram for CastorTTRecord:

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

void	getEnergy_fC (double energy[16][14], edm::Handle< CastorDigiCollection > &CastorDigiColl, edm::Event &e, const edm::EventSetup &c)

void	getTriggerDecisions (std::vector< bool > &decision, double energy[16][14]) const

void	getTriggerDecisionsPerOctant (std::vector< bool > tdps[16], double energy[16][14]) const

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

	~CastorTTRecord () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::EDGetTokenT< CastorDigiCollection >	CastorDigiColl_

unsigned int	CastorSignalTS_

double	reweighted_gain

std::vector< std::string >	TrigNames_

std::vector< double >	TrigThresholds_

std::vector< unsigned int >	ttpBits_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 11 of file CastorTTRecord.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file CastorTTRecord.cc.

References CastorDigiColl_, CastorSignalTS_, edm::ParameterSet::getParameter(), reweighted_gain, TrigNames_, TrigThresholds_, and ttpBits_.

                                                         {
   CastorDigiColl_ = consumes<CastorDigiCollection>(ps.getParameter<edm::InputTag>("CastorDigiCollection"));
   CastorSignalTS_ = ps.getParameter<unsigned int>("CastorSignalTS");
 
   ttpBits_ = ps.getParameter<std::vector<unsigned int>>("ttpBits");
   TrigNames_ = ps.getParameter<std::vector<std::string>>("TriggerBitNames");
   TrigThresholds_ = ps.getParameter<std::vector<double>>("TriggerThresholds");
 
   reweighted_gain = 1.0;
 
   produces<L1GtTechnicalTriggerRecord>();
 }

CastorTTRecord::~CastorTTRecord ( )

override

Definition at line 27 of file CastorTTRecord.cc.

27 {}

Member Function Documentation

void CastorTTRecord::getEnergy_fC	(	double	energy[16][14],
		edm::Handle< CastorDigiCollection > &	CastorDigiColl,
		edm::Event &	e,
		const edm::EventSetup &	c
	)

Definition at line 64 of file CastorTTRecord.cc.

References CastorCoderDb::adc2fC(), edm::SortedCollection< T, SORT >::begin(), AlignmentProducer_cff::calibrations, CastorSignalTS_, edm::SortedCollection< T, SORT >::end(), CastorCalibrations::gain(), edm::EventSetup::get(), CastorDbService::getCastorCalibrations(), CastorDbService::getCastorCoder(), CastorDbService::getCastorShape(), CastorDataFrame::id(), HcalCastorDetId::module(), CastorCalibrations::pedestal(), reweighted_gain, and HcalCastorDetId::sector().

Referenced by produce().

                                                                    {
   // std::cerr << "**** RUNNING THROUGH CastorTTRecord::getEnergy_fC" <<
   // std::endl;
 
   // Get Conditions
   edm::ESHandle<CastorDbService> conditions;
   eventSetup.get<CastorDbRecord>().get(conditions);
   const CastorQIEShape *shape = conditions->getCastorShape();  // this one is generic
 
   for (int isec = 0; isec < 16; isec++)
     for (int imod = 0; imod < 14; imod++)
       energy[isec][imod] = 0;
 
   // Loop over digis
   CastorDigiCollection::const_iterator idigi;
   for (idigi = CastorDigiColl->begin(); idigi != CastorDigiColl->end(); idigi++) {
     const CastorDataFrame &digi = (*idigi);
     HcalCastorDetId cell = digi.id();
 
     // Get Castor Coder
     const CastorQIECoder *channelCoder = conditions->getCastorCoder(cell);
     CastorCoderDb coder(*channelCoder, *shape);
 
     // Get Castor Calibration
     const CastorCalibrations &calibrations = conditions->getCastorCalibrations(cell);
 
     // convert adc to fC
     CaloSamples tool;
     coder.adc2fC(digi, tool);
 
     // pedestal substraction
     int capid = digi[CastorSignalTS_].capid();
     double fC = tool[CastorSignalTS_] - calibrations.pedestal(capid);
 
     // to correct threshold levels in fC for different gains
     reweighted_gain = calibrations.gain(capid) / 0.015;
 
     energy[digi.id().sector() - 1][digi.id().module() - 1] = fC;
   }
 }

void CastorTTRecord::getTriggerDecisions	(	std::vector< bool > &	decision,
		double	energy[16][14]
	)		const

Definition at line 108 of file CastorTTRecord.cc.

References getTriggerDecisionsPerOctant().

Referenced by produce().

                                                                                                {
   // std::cerr << "**** RUNNING THROUGH CastorTTRecord::getTriggerDecisions" <<
   // std::endl;
 
   // check if number of bits is at least four
   if (decision.size() < 4)
     return;
 
   std::vector<bool> tdpo[8];  // TriggerDecisionsPerOctant
   getTriggerDecisionsPerOctant(tdpo, energy);
 
   // preset trigger decisions
   decision.at(0) = true;
   decision.at(1) = false;
   decision.at(2) = false;
   decision.at(3) = false;
 
   // canceld for low pt jet
   // bool EM_decision = false;
   // bool HAD_decision = false;
   // loop over castor octants
   for (int ioct = 0; ioct < 8; ioct++) {
     int next_oct = (ioct + 1) % 8;
     int prev_oct = (ioct + 8 - 1) % 8;
 
     // gap Trigger
     if (!tdpo[ioct].at(0))
       decision.at(0) = false;
     if (!tdpo[ioct].at(1))
       decision.at(0) = false;
 
     // jet Trigger
     if (tdpo[ioct].at(2))
       decision.at(1) = true;
 
     // electron
     // canceld for low pt jet
     // if( tdpo[ioct].at(3) ) EM_decision = true;
     // if( tdpo[ioct].at(4) ) HAD_decision = true;
 
     // iso muon
     if (tdpo[ioct].at(5)) {
       // was one of the other sectors
       // in the octant empty ?
       if (tdpo[ioct].at(0)) {
         if (tdpo[prev_oct].at(1) && tdpo[next_oct].at(0) && tdpo[next_oct].at(1))
           decision.at(3) = true;
       } else if (tdpo[ioct].at(1)) {
         if (tdpo[prev_oct].at(0) && tdpo[prev_oct].at(1) && tdpo[next_oct].at(0))
           decision.at(3) = true;
       }
       // when not no iso muon
     }
 
     // low pt jet Trigger
     if (tdpo[ioct].at(6))
       decision.at(2) = true;
   }
 
   // for EM Trigger whole castor not hadronic and somewhere EM
   // canceld for low pt jet
   // decision.at(2) = EM_decision && !HAD_decision;
 }

void CastorTTRecord::getTriggerDecisionsPerOctant	(	std::vector< bool >	tdps[16],
		double	energy[16][14]
	)		const

Definition at line 172 of file CastorTTRecord.cc.

References reweighted_gain, and TrigThresholds_.

Referenced by getTriggerDecisions().

                                                                                                       {
   // std::cerr << "**** RUNNING THROUGH
   // CastorTTRecord::getTriggerDecisionsPerOctant" << std::endl;
 
   // loop over octatants
   for (int ioct = 0; ioct < 8; ioct++) {
     // six bits from HTR card
     // 0. first sector empty
     // 1. second sector empty
     // 2. jet any sector
     // 3. EM any sector
     // 4. HAD any sector
     // 5. muon any sector
     // add instead of EM Trigger (not bit 6 in real)
     // 6. low pt jet any sector
     tdpo[ioct].resize(7);
 
     for (int ibit = 0; ibit < 7; ibit++)
       tdpo[ioct].at(ibit) = false;
 
     // loop over castor sectors in octant
     for (int ioctsec = 0; ioctsec < 2; ioctsec++) {
       // absolute sector number
       int isec = 2 * ioct + ioctsec;
 
       // init module sums for every sector
       double fCsum_mod = 0;
       double fCsum_em = 0, fCsum_ha = 0;
       double fCsum_jet_had = 0;
       double fCsum_col[3] = {0, 0, 0};
 
       // loop over modules
       for (int imod = 0; imod < 14; imod++) {
         // total sum
         fCsum_mod += energy[isec][imod];
 
         // EM & HAD sum
         if (imod < 2)
           fCsum_em += energy[isec][imod];
         if (imod > 2 && imod < 12)
           fCsum_ha += energy[isec][imod];
 
         // sum over three sector parts
         if (imod < 4)
           fCsum_col[0] += energy[isec][imod];
         else if (imod < 8)
           fCsum_col[1] += energy[isec][imod];
         else if (imod < 12)
           fCsum_col[2] += energy[isec][imod];
 
         // HAD sum for jet trigger v2
         if (imod > 1 && imod < 5)
           fCsum_jet_had += energy[isec][imod];
       }
 
       // gap Trigger
       if (fCsum_mod < TrigThresholds_.at(0)) {
         if (ioctsec == 0)
           tdpo[ioct].at(0) = true;
         else if (ioctsec == 1)
           tdpo[ioct].at(1) = true;
       }
 
       // jet Trigger
       // with gain correction
       /* old version of jet trigger ( deprecated because of saturation )
       if( fCsum_mod > TrigThresholds_.at(1) / reweighted_gain )
           tdpo[ioct].at(2) = true;
       */
       if (fCsum_jet_had > TrigThresholds_.at(1) / reweighted_gain)
         // additional high threshold near saturation for EM part
         if (energy[isec][0] > 26000 / reweighted_gain && energy[isec][1] > 26000 / reweighted_gain)
           tdpo[ioct].at(2) = true;
 
       // low pt jet Trigger
       if (fCsum_mod > TrigThresholds_.at(5) / reweighted_gain)
         tdpo[ioct].at(6) = true;
 
       // egamma Trigger
       // with gain correction only in the EM threshold
       if (fCsum_em > TrigThresholds_.at(2) / reweighted_gain)
         tdpo[ioct].at(3) = true;
       if (fCsum_ha > TrigThresholds_.at(3))
         tdpo[ioct].at(4) = true;
 
       // muon Trigger
       int countColumns = 0;
       for (int icol = 0; icol < 3; icol++)
         if (fCsum_col[icol] > TrigThresholds_.at(4))
           countColumns++;
       if (countColumns >= 2)
         tdpo[ioct].at(5) = true;
     }
   }
 }

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

override

Definition at line 29 of file CastorTTRecord.cc.

References CastorDigiColl_, edm::Event::getByToken(), getEnergy_fC(), getTriggerDecisions(), mps_fire::i, eostools::move(), convertSQLitetoXML_cfg::output, edm::Event::put(), TrigNames_, and ttpBits_.

                                                                          {
   // std::cerr << "**** RUNNING THROUGH CastorTTRecord::produce" << std::endl;
 
   std::vector<L1GtTechnicalTrigger> vecTT(ttpBits_.size());
 
   // Get Inputs
   edm::Handle<CastorDigiCollection> CastorDigiColl;
   e.getByToken(CastorDigiColl_, CastorDigiColl);
 
   if (!CastorDigiColl.failedToGet()) {
     double cas_efC[16][14];
     getEnergy_fC(cas_efC, CastorDigiColl, e, eventSetup);
 
     std::vector<bool> decision(ttpBits_.size());
 
     getTriggerDecisions(decision, cas_efC);
 
     for (unsigned int i = 0; i < ttpBits_.size(); i++) {
       // if( decision.at(i) ) std::cerr << "**** Something Triggered" <<
       // std::endl; std::cout << "Run CastorTTRecord::produce. TriggerBit = " <<
       // ttpBits_.at(i) << "; TriggerName = " << TrigNames_.at(i) << "; Decision
       // = " << decision[i] << std::endl;
       vecTT.at(i) = L1GtTechnicalTrigger(TrigNames_.at(i), ttpBits_.at(i), 0, decision.at(i));
     }
 
   } else {
     vecTT.clear();
   }
 
   // Put output into event
   std::unique_ptr<L1GtTechnicalTriggerRecord> output(new L1GtTechnicalTriggerRecord());
   output->setGtTechnicalTrigger(vecTT);
   e.put(std::move(output));
 }