L1RCTTPGProvider Class Reference

#include <L1RCTTPGProvider.h>

Inheritance diagram for L1RCTTPGProvider:

Public Member Functions
	L1RCTTPGProvider (const edm::ParameterSet &)
	~L1RCTTPGProvider () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () 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
std::vector< bool > const &	recordProvenanceList () 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)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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 Member Functions
void	endJob () override
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::InputTag	ecalTPG_
int	hbShift
edm::InputTag	hcalTPG_
int	hfShift
int	postSamples
int	preSamples
bool	useEcalCosmicTiming
bool	useHcalCosmicTiming

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

L1Trigger/L1RCTTPGProvider/src/L1RCTTPGProvider.cc

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 35 of file L1RCTTPGProvider.h.

Constructor & Destructor Documentation

L1RCTTPGProvider()

L1RCTTPGProvider::L1RCTTPGProvider ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 5 of file L1RCTTPGProvider.cc.

    : ecalTPG_(iConfig.getParameter<edm::InputTag>("ecalTPGs")),
      hcalTPG_(iConfig.getParameter<edm::InputTag>("hcalTPGs")),
      useHcalCosmicTiming(iConfig.getParameter<bool>("useHCALCosmicTiming")),
      useEcalCosmicTiming(iConfig.getParameter<bool>("useECALCosmicTiming")),
      preSamples(iConfig.getParameter<int>("preSamples")),
      postSamples(iConfig.getParameter<int>("postSamples")),
      hfShift(iConfig.getParameter<int>("HFShift")),
      hbShift(iConfig.getParameter<int>("HBShift")) {
  // Output :The new manipulated TPGs
  // make it smart - to name the collections
  // correctly
  char ecal_label[200];
  char hcal_label[200];
 
  for (int i = preSamples; i > 0; --i) {
    sprintf(ecal_label, "ECALBxminus%d", i);
    sprintf(hcal_label, "HCALBxminus%d", i);
    produces<EcalTrigPrimDigiCollection>(ecal_label);
    produces<HcalTrigPrimDigiCollection>(hcal_label);
  }
 
  produces<EcalTrigPrimDigiCollection>("ECALBx0");
  produces<HcalTrigPrimDigiCollection>("HCALBx0");
 
  for (int i = 0; i < postSamples; ++i) {
    sprintf(ecal_label, "ECALBxplus%d", i + 1);
    sprintf(hcal_label, "HCALBxplus%d", i + 1);
    produces<EcalTrigPrimDigiCollection>(ecal_label);
    produces<HcalTrigPrimDigiCollection>(hcal_label);
  }
}

References mps_fire::i, postSamples, and preSamples.

~L1RCTTPGProvider()

L1RCTTPGProvider::~L1RCTTPGProvider ( )

override

Definition at line 38 of file L1RCTTPGProvider.cc.

                                    {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

endJob()

void L1RCTTPGProvider::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 310 of file L1RCTTPGProvider.cc.

{}

produce()

void L1RCTTPGProvider::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 48 of file L1RCTTPGProvider.cc.

                                                                            {
  using namespace edm;
  // Declare handles
  Handle<EcalTrigPrimDigiCollection> ecal;
  Handle<HcalTrigPrimDigiCollection> hcal;
 
  // Declare vector of collection to send for output !
 
  std::vector<EcalTrigPrimDigiCollection> ecalColl(preSamples + 1 + postSamples);
  std::vector<HcalTrigPrimDigiCollection> hcalColl(preSamples + 1 + postSamples);
 
  unsigned nSamples = preSamples + postSamples + 1;
 
  if (iEvent.getByLabel(ecalTPG_, ecal))
    if (ecal.isValid()) {
      // loop through all ecal digis
      for (EcalTrigPrimDigiCollection::const_iterator ecal_it = ecal->begin(); ecal_it != ecal->end(); ecal_it++) {
        short zside = ecal_it->id().zside();
        unsigned short ietaAbs = ecal_it->id().ietaAbs();
        short iphi = ecal_it->id().iphi();
        unsigned short digiSize = ecal_it->size();
        unsigned short nSOI = (unsigned short)(ecal_it->sampleOfInterest());
        if (digiSize < nSamples || nSOI < preSamples || ((int)(digiSize - nSOI) < (int)(nSamples - preSamples))) {
          unsigned short preLoopsZero = (unsigned short)(preSamples)-nSOI;
 
          // fill extra bx's at beginning with zeros
          for (int sample = 0; sample < preLoopsZero; sample++) {
            // fill first few with zeros
            EcalTriggerPrimitiveDigi ecalDigi(
                EcalTrigTowerDetId((int)zside, EcalTriggerTower, (int)ietaAbs, (int)iphi));
            ecalDigi.setSize(1);
            ecalDigi.setSample(0, EcalTriggerPrimitiveSample(0, false, 0));
            ecalColl[sample].push_back(ecalDigi);
          }
 
          // loop through existing data
          for (int sample = preLoopsZero; sample < (preLoopsZero + digiSize); sample++) {
            // go through data
            EcalTriggerPrimitiveDigi ecalDigi(
                EcalTrigTowerDetId((int)zside, EcalTriggerTower, (int)ietaAbs, (int)iphi));
            ecalDigi.setSize(1);
            if (useEcalCosmicTiming && iphi >= 1 && iphi <= 36) {
              if (nSOI < (preSamples + 1)) {
                edm::LogWarning("TooLittleData") << "ECAL data needs at least one presample "
                                                 << "more than the number requested "
                                                 << "to use ecal cosmic timing mod!  "
                                                 << "reverting to useEcalCosmicTiming = false "
                                                 << "for rest of job.";
                useEcalCosmicTiming = false;
              } else {
                ecalDigi.setSample(0,
                                   EcalTriggerPrimitiveSample(ecal_it->sample(nSOI + sample - preSamples - 1).raw()));
              }
            }
            // else
            if ((!useEcalCosmicTiming) || (iphi >= 37 && iphi <= 72)) {
              ecalDigi.setSample(0, EcalTriggerPrimitiveSample(ecal_it->sample(nSOI + sample - preSamples).raw()));
            }
            ecalColl[sample].push_back(ecalDigi);
          }
 
          // fill extra bx's at end with zeros
          for (unsigned int sample = (preLoopsZero + digiSize); sample < nSamples; sample++) {
            // fill zeros!
            EcalTriggerPrimitiveDigi ecalDigi(
                EcalTrigTowerDetId((int)zside, EcalTriggerTower, (int)ietaAbs, (int)iphi));
            ecalDigi.setSize(1);
            ecalDigi.setSample(0, EcalTriggerPrimitiveSample(0, false, 0));
            ecalColl[sample].push_back(ecalDigi);
          }
        } else {
          for (unsigned short sample = 0; sample < nSamples; sample++) {
            // put each time sample into its own digi
            short zside = ecal_it->id().zside();
            unsigned short ietaAbs = ecal_it->id().ietaAbs();
            short iphi = ecal_it->id().iphi();
            EcalTriggerPrimitiveDigi ecalDigi(
                EcalTrigTowerDetId((int)zside, EcalTriggerTower, (int)ietaAbs, (int)iphi));
            ecalDigi.setSize(1);
 
            if (useEcalCosmicTiming && iphi >= 1 && iphi <= 36) {
              if (nSOI < (preSamples + 1)) {
                edm::LogWarning("TooLittleData") << "ECAL data needs at least one presample "
                                                 << "more than the number requested "
                                                 << "to use ecal cosmic timing mod!  "
                                                 << "reverting to useEcalCosmicTiming = false "
                                                 << "for rest of job.";
                useEcalCosmicTiming = false;
              } else {
                ecalDigi.setSample(0,
                                   EcalTriggerPrimitiveSample(
                                       ecal_it->sample(ecal_it->sampleOfInterest() + sample - preSamples - 1).raw()));
              }
            }
            // else
            if ((!useEcalCosmicTiming) || (iphi >= 37 && iphi <= 72)) {
              ecalDigi.setSample(
                  0,
                  EcalTriggerPrimitiveSample(ecal_it->sample(ecal_it->sampleOfInterest() + sample - preSamples).raw()));
            }
            // push back each digi into correct "time sample" of coll
            ecalColl[sample].push_back(ecalDigi);
          }
        }
      }
    }
 
  if (iEvent.getByLabel(hcalTPG_, hcal))
    if (hcal.isValid()) {
      // loop through all hcal digis
      for (HcalTrigPrimDigiCollection::const_iterator hcal_it = hcal->begin(); hcal_it != hcal->end(); hcal_it++) {
        short ieta = hcal_it->id().ieta();
        short iphi = hcal_it->id().iphi();
        // loop through time samples for each digi
        unsigned short digiSize = hcal_it->size();
        // (size of each digi must be no less than nSamples)
        unsigned short nSOI = (unsigned short)(hcal_it->presamples());
        if (digiSize < nSamples || nSOI < preSamples || ((int)(digiSize - nSOI) < (int)(nSamples - preSamples))) {
          unsigned short preLoopsZero = (unsigned short)(preSamples)-nSOI;
          // fill extra bx's at beginning with zeros
          for (int sample = 0; sample < preLoopsZero; sample++) {
            // fill first few with zeros
            HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int)ieta, (int)iphi));
            hcalDigi.setSize(1);
            hcalDigi.setPresamples(0);
            hcalDigi.setSample(0, HcalTriggerPrimitiveSample(0, false, 0, 0));
            hcalColl[sample].push_back(hcalDigi);
          }
 
          // loop through existing data
          for (int sample = preLoopsZero; sample < (preLoopsZero + digiSize); sample++) {
            // go through data
            HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int)ieta, (int)iphi));
            hcalDigi.setSize(1);
            hcalDigi.setPresamples(0);
 
            if (useHcalCosmicTiming && iphi >= 1 && iphi <= 36) {
              if (nSOI < (preSamples + 1)) {
                edm::LogWarning("TooLittleData") << "HCAL data needs at least one presample "
                                                 << "more than the number requested "
                                                 << "to use hcal cosmic timing mod!  "
                                                 << "reverting to useHcalCosmicTiming = false "
                                                 << "for rest of job.";
                useHcalCosmicTiming = false;
              } else {
                hcalDigi.setSample(
                    0,
                    HcalTriggerPrimitiveSample(hcal_it->sample(hcal_it->presamples() + sample - preSamples - 1).raw()));
              }
            }
            // else
            if ((!useHcalCosmicTiming) || (iphi >= 37 && iphi <= 72)) {
              hcalDigi.setSample(
                  0, HcalTriggerPrimitiveSample(hcal_it->sample(hcal_it->presamples() + sample - preSamples).raw()));
            }
            hcalColl[sample].push_back(hcalDigi);
          }
 
          // fill extra bx's at end with zeros
          for (unsigned int sample = (preLoopsZero + digiSize); sample < nSamples; sample++) {
            // fill zeros!
            HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int)ieta, (int)iphi));
            hcalDigi.setSize(1);
            hcalDigi.setPresamples(0);
            hcalDigi.setSample(0, HcalTriggerPrimitiveSample(0, false, 0, 0));
            hcalColl[sample].push_back(hcalDigi);
          }
        } else {
          for (unsigned short sample = 0; sample < nSamples; sample++) {
            // put each (relevant) time sample into its own digi
            HcalTriggerPrimitiveDigi hcalDigi(HcalTrigTowerDetId((int)ieta, (int)iphi));
            hcalDigi.setSize(1);
            hcalDigi.setPresamples(0);
 
            if (useHcalCosmicTiming && iphi >= 1 && iphi <= 36) {
              if (nSOI < (preSamples + 1)) {
                edm::LogWarning("TooLittleData") << "HCAL data needs at least one presample "
                                                 << "more than the number requested "
                                                 << "to use hcal cosmic timing mod!  "
                                                 << "reverting to useHcalCosmicTiming = false "
                                                 << "for rest of job.";
                useHcalCosmicTiming = false;
              } else {
                hcalDigi.setSample(
                    0,
                    HcalTriggerPrimitiveSample(hcal_it->sample(hcal_it->presamples() + sample - preSamples - 1).raw()));
              }
            }
            // else
            if ((!useHcalCosmicTiming) || (iphi >= 37 && iphi <= 72)) {
              if (ieta > -29 && ieta < 29)
                hcalDigi.setSample(0,
                                   HcalTriggerPrimitiveSample(
                                       hcal_it->sample(hcal_it->presamples() + sample - preSamples + hbShift).raw()));
              if (ieta <= -29 || ieta >= 29)
                hcalDigi.setSample(0,
                                   HcalTriggerPrimitiveSample(
                                       hcal_it->sample(hcal_it->presamples() + sample - preSamples + hfShift).raw()));
            }
            hcalColl[sample].push_back(hcalDigi);
          }
        }
      }
    }
 
  // Now put the events back to file
 
  for (int i = 0; i < preSamples; ++i) {
    char ecal_label[200];
    char hcal_label[200];
 
    sprintf(ecal_label, "ECALBxminus%d", preSamples - i);
    sprintf(hcal_label, "HCALBxminus%d", preSamples - i);
 
    std::unique_ptr<EcalTrigPrimDigiCollection> ecalIn(new EcalTrigPrimDigiCollection);
    std::unique_ptr<HcalTrigPrimDigiCollection> hcalIn(new HcalTrigPrimDigiCollection);
    for (unsigned int j = 0; j < ecalColl[i].size(); ++j) {
      ecalIn->push_back((ecalColl[i])[j]);
    }
    for (unsigned int j = 0; j < hcalColl[i].size(); ++j)
      hcalIn->push_back((hcalColl[i])[j]);
 
    iEvent.put(std::move(ecalIn), ecal_label);
    iEvent.put(std::move(hcalIn), hcal_label);
  }
 
  std::unique_ptr<EcalTrigPrimDigiCollection> ecal0(new EcalTrigPrimDigiCollection);
  std::unique_ptr<HcalTrigPrimDigiCollection> hcal0(new HcalTrigPrimDigiCollection);
  for (unsigned int j = 0; j < ecalColl[preSamples].size(); ++j)
    ecal0->push_back((ecalColl[preSamples])[j]);
  for (unsigned int j = 0; j < hcalColl[preSamples].size(); ++j)
    hcal0->push_back((hcalColl[preSamples])[j]);
 
  iEvent.put(std::move(ecal0), "ECALBx0");
  iEvent.put(std::move(hcal0), "HCALBx0");
 
  for (int i = preSamples + 1; i < preSamples + postSamples + 1; ++i) {
    char ecal_label[200];
    char hcal_label[200];
 
    sprintf(ecal_label, "ECALBxplus%d", i - preSamples);
    sprintf(hcal_label, "HCALBxplus%d", i - preSamples);
 
    std::unique_ptr<EcalTrigPrimDigiCollection> ecalIn2(new EcalTrigPrimDigiCollection);
    std::unique_ptr<HcalTrigPrimDigiCollection> hcalIn2(new HcalTrigPrimDigiCollection);
 
    for (unsigned int j = 0; j < ecalColl[i].size(); ++j)
      ecalIn2->push_back((ecalColl[i])[j]);
 
    for (unsigned int j = 0; j < hcalColl[i].size(); ++j)
      hcalIn2->push_back((hcalColl[i])[j]);
 
    iEvent.put(std::move(ecalIn2), ecal_label);
    iEvent.put(std::move(hcalIn2), hcal_label);
  }
}

References bsc_activity_cfg::ecal, ecalTPG_, EcalTriggerTower, hbShift, patCandidatesForDimuonsSequences_cff::hcal, hcalTPG_, hfShift, mps_fire::i, LEDCalibrationChannels::ieta, iEvent, LEDCalibrationChannels::iphi, dqmiolumiharvest::j, eostools::move(), PresampleTask_cfi::nSamples, postSamples, preSamples, simplePhotonAnalyzer_cfi::sample, HcalTriggerPrimitiveDigi::setPresamples(), EcalTriggerPrimitiveDigi::setSample(), HcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), HcalTriggerPrimitiveDigi::setSize(), useEcalCosmicTiming, useHcalCosmicTiming, and ecaldqm::zside().

Member Data Documentation

ecalTPG_

edm::InputTag L1RCTTPGProvider::ecalTPG_

private

Definition at line 45 of file L1RCTTPGProvider.h.

Referenced by produce().

hbShift

int L1RCTTPGProvider::hbShift

private

Definition at line 52 of file L1RCTTPGProvider.h.

Referenced by produce().

hcalTPG_

edm::InputTag L1RCTTPGProvider::hcalTPG_

private

Definition at line 46 of file L1RCTTPGProvider.h.

Referenced by produce().

hfShift

int L1RCTTPGProvider::hfShift

private

Definition at line 51 of file L1RCTTPGProvider.h.

Referenced by produce().

postSamples

int L1RCTTPGProvider::postSamples

private

Definition at line 50 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

preSamples

int L1RCTTPGProvider::preSamples

private

Definition at line 49 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

useEcalCosmicTiming

bool L1RCTTPGProvider::useEcalCosmicTiming

private

Definition at line 48 of file L1RCTTPGProvider.h.

Referenced by produce().

useHcalCosmicTiming

bool L1RCTTPGProvider::useHcalCosmicTiming

private

Definition at line 47 of file L1RCTTPGProvider.h.

Referenced by produce().