#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

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

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 ( const edm::ParameterSet & iConfig )

explicit

Definition at line 5 of file L1RCTTPGProvider.cc.

References mps_fire::i, postSamples, and preSamples.

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

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

void L1RCTTPGProvider::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 310 of file L1RCTTPGProvider.cc.

310 {}

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 48 of file L1RCTTPGProvider.cc.

References digitizers_cfi::ecal, ecalTPG_, EcalTriggerTower, edm::Event::getByLabel(), hbShift, digitizers_cfi::hcal, hcalTPG_, hfShift, mps_fire::i, dqmiolumiharvest::j, eostools::move(), hgc_digi::nSamples, postSamples, preSamples, edm::Event::put(), compare_using_db::sample, HcalTriggerPrimitiveDigi::setPresamples(), EcalTriggerPrimitiveDigi::setSample(), HcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), HcalTriggerPrimitiveDigi::setSize(), useEcalCosmicTiming, useHcalCosmicTiming, and ecaldqm::zside().

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

Member Data Documentation

edm::InputTag L1RCTTPGProvider::ecalTPG_

private

Definition at line 45 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::hbShift

private

Definition at line 52 of file L1RCTTPGProvider.h.

Referenced by produce().

edm::InputTag L1RCTTPGProvider::hcalTPG_

private

Definition at line 46 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::hfShift

private

Definition at line 51 of file L1RCTTPGProvider.h.

Referenced by produce().

int L1RCTTPGProvider::postSamples

private

Definition at line 50 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

int L1RCTTPGProvider::preSamples

private

Definition at line 49 of file L1RCTTPGProvider.h.

Referenced by L1RCTTPGProvider(), and produce().

bool L1RCTTPGProvider::useEcalCosmicTiming

private

Definition at line 48 of file L1RCTTPGProvider.h.

Referenced by produce().

bool L1RCTTPGProvider::useHcalCosmicTiming

private

Definition at line 47 of file L1RCTTPGProvider.h.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation