#include <EventFilter/GctRawToDigi/src/GctDigiToRaw.cc>

Inheritance diagram for GctDigiToRaw:

Public Member Functions
	GctDigiToRaw (const edm::ParameterSet &)

	~GctDigiToRaw () 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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override

void	endJob () override

void	print (FEDRawData &data)

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

Private Attributes
int	counter_

int	fedId_

GctFormatTranslateMCLegacy	formatTranslator_

bool	packRctCalo_

bool	packRctEm_

edm::EDGetTokenT< L1CaloEmCollection >	tokenCaloEm_

edm::EDGetTokenT< L1CaloRegionCollection >	tokenCaloRegion_

edm::EDGetTokenT< L1GctEtHadCollection >	tokenGctEtHad_

edm::EDGetTokenT< L1GctEtMissCollection >	tokenGctEtMiss_

edm::EDGetTokenT< L1GctEtTotalCollection >	tokenGctEtTotal_

edm::EDGetTokenT< L1GctHFBitCountsCollection >	tokenGctHFBitCounts_

edm::EDGetTokenT< L1GctHFRingEtSumsCollection >	tokenGctHFRingEtSums_

edm::EDGetTokenT< L1GctHtMissCollection >	tokenGctHtMiss_

edm::EDGetTokenT< L1GctJetCandCollection >	tokenGctJetCand_cenJets_

edm::EDGetTokenT< L1GctJetCandCollection >	tokenGctJetCand_forJets_

edm::EDGetTokenT< L1GctJetCandCollection >	tokenGctJetCand_tauJets_

edm::EDGetTokenT< L1GctJetCountsCollection >	tokenGctJetCounts_

edm::EDGetTokenT< L1GctEmCandCollection >	tokenL1GctEmCand_isoEm_

edm::EDGetTokenT< L1GctEmCandCollection >	tokenL1GctEmCand_nonIsoEm_

bool	verbose_

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::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 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

Description: Produce fake GCT raw data from digis

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

Definition at line 41 of file GctDigiToRaw.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 40 of file GctDigiToRaw.cc.

References fedId_, edm::ParameterSet::getParameter(), LogDebug, packRctCalo_, packRctEm_, AlCaHLTBitMon_QueryRunRegistry::string, tokenCaloEm_, tokenCaloRegion_, tokenGctEtHad_, tokenGctEtMiss_, tokenGctEtTotal_, tokenGctHFBitCounts_, tokenGctHFRingEtSums_, tokenGctHtMiss_, tokenGctJetCand_cenJets_, tokenGctJetCand_forJets_, tokenGctJetCand_tauJets_, tokenGctJetCounts_, tokenL1GctEmCand_isoEm_, and tokenL1GctEmCand_nonIsoEm_.

                                                          :
   packRctEm_(iConfig.getUntrackedParameter<bool>("packRctEm", true)),
   packRctCalo_(iConfig.getUntrackedParameter<bool>("packRctCalo", true)),
   fedId_(iConfig.getParameter<int>("gctFedId")),
   verbose_(iConfig.getUntrackedParameter<bool>("verbose",false)),
   counter_(0),
   formatTranslator_()
 {
   LogDebug("GCT") << "GctDigiToRaw will pack FED Id " << fedId_;
 
   //register the products
   produces<FEDRawDataCollection>();
   const edm::InputTag rctInputTag = iConfig.getParameter<edm::InputTag>("rctInputLabel");
   const edm::InputTag gctInputTag = iConfig.getParameter<edm::InputTag>("gctInputLabel");
   const std::string& gctInputLabelStr = gctInputTag.label();
   tokenL1GctEmCand_isoEm_ = consumes<L1GctEmCandCollection>(edm::InputTag(gctInputLabelStr, "isoEm"));
   tokenL1GctEmCand_nonIsoEm_ = consumes<L1GctEmCandCollection>(edm::InputTag(gctInputLabelStr, "nonIsoEm"));
   tokenGctJetCand_cenJets_ = consumes<L1GctJetCandCollection>(edm::InputTag(gctInputLabelStr, "cenJets"));
   tokenGctJetCand_forJets_ = consumes<L1GctJetCandCollection>(edm::InputTag(gctInputLabelStr, "forJets"));
   tokenGctJetCand_tauJets_ = consumes<L1GctJetCandCollection>(edm::InputTag(gctInputLabelStr, "tauJets"));
   tokenGctEtTotal_ = consumes<L1GctEtTotalCollection>(gctInputTag);
   tokenGctEtHad_ = consumes<L1GctEtHadCollection>(gctInputTag);
   tokenGctEtMiss_ = consumes<L1GctEtMissCollection>(gctInputTag);
   tokenGctHFRingEtSums_ = consumes<L1GctHFRingEtSumsCollection>(gctInputTag);
   tokenGctHFBitCounts_ = consumes<L1GctHFBitCountsCollection>(gctInputTag);
   tokenGctHtMiss_ = consumes<L1GctHtMissCollection>(gctInputTag);
   tokenGctJetCounts_ = consumes<L1GctJetCountsCollection>(gctInputTag);
   if(packRctEm_) {
     tokenCaloEm_ = consumes<L1CaloEmCollection>(rctInputTag);
   }
   if(packRctCalo_) {
     tokenCaloRegion_ = consumes<L1CaloRegionCollection>(rctInputTag);
   }
 }

GctDigiToRaw::~GctDigiToRaw ( )

override

Definition at line 76 of file GctDigiToRaw.cc.

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

Member Function Documentation

void GctDigiToRaw::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 232 of file GctDigiToRaw.cc.

233 {

234 }

void GctDigiToRaw::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 238 of file GctDigiToRaw.cc.

References DEFINE_FWK_MODULE.

238 {

239 }

void GctDigiToRaw::print ( FEDRawData & data )

private

Definition at line 218 of file GctDigiToRaw.cc.

References gather_cfg::cout, edmIntegrityCheck::d, FEDRawData::data(), mps_fire::i, FEDRawData::size(), and w.

Referenced by produce().

                                          {
 
   const unsigned char * d = data.data();
 
   for (unsigned int i=0; i<data.size(); i=i+4) {
     uint32_t w = (uint32_t)d[i] + (uint32_t)(d[i+1]<<8) + (uint32_t)(d[i+2]<<16) + (uint32_t)(d[i+3]<<24);
     cout << std::hex << std::setw(4) << i/4 << " " << std::setw(8) << w << endl;
   }
 
 }

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

overrideprivate

Definition at line 89 of file GctDigiToRaw.cc.

References evf::compute_crc(), counter_, FEDRawData::data(), edm::EventID::event(), edm::HandleBase::failedToGet(), fedId_, formatTranslator_, edm::Event::getByToken(), edm::EventBase::id(), LogDebug, eostools::move(), packRctCalo_, packRctEm_, print(), edm::Handle< T >::product(), edm::Event::put(), FEDRawData::resize(), FEDHeader::set(), FEDTrailer::set(), GctFormatTranslateBase::setPackingBxId(), GctFormatTranslateBase::setPackingEventId(), tokenCaloEm_, tokenCaloRegion_, tokenGctEtHad_, tokenGctEtMiss_, tokenGctEtTotal_, tokenGctHFBitCounts_, tokenGctHFRingEtSums_, tokenGctHtMiss_, tokenGctJetCand_cenJets_, tokenGctJetCand_forJets_, tokenGctJetCand_tauJets_, tokenGctJetCounts_, tokenL1GctEmCand_isoEm_, tokenL1GctEmCand_nonIsoEm_, verbose_, GctFormatTranslateMCLegacy::writeAllRctCaloRegionBlock(), GctFormatTranslateMCLegacy::writeGctOutEmAndEnergyBlock(), GctFormatTranslateMCLegacy::writeGctOutJetBlock(), and GctFormatTranslateMCLegacy::writeRctEmCandBlocks().

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

 {
   using namespace edm;
   
   counter_++; // To "simulate" bunch crossings for now...
   unsigned int bx = counter_ % 3564;  // What's the proper way of doing this?
   EventNumber_t eventNumber = iEvent.id().event();
   
   // Supply bx and EvID to the packer so it can make internal capture block headers.
   formatTranslator_.setPackingBxId(bx);
   formatTranslator_.setPackingEventId(eventNumber);
  
   // get GCT digis
   edm::Handle<L1GctEmCandCollection> isoEm;
   iEvent.getByToken(tokenL1GctEmCand_isoEm_, isoEm);
   edm::Handle<L1GctEmCandCollection> nonIsoEm;
   iEvent.getByToken(tokenL1GctEmCand_nonIsoEm_, nonIsoEm);
   edm::Handle<L1GctJetCandCollection> cenJets;
   iEvent.getByToken(tokenGctJetCand_cenJets_, cenJets);
   edm::Handle<L1GctJetCandCollection> forJets;
   iEvent.getByToken(tokenGctJetCand_forJets_, forJets);
   edm::Handle<L1GctJetCandCollection> tauJets;
   iEvent.getByToken(tokenGctJetCand_tauJets_, tauJets);
   edm::Handle<L1GctEtTotalCollection> etTotal;
   iEvent.getByToken(tokenGctEtTotal_,  etTotal);
   edm::Handle<L1GctEtHadCollection> etHad;
   iEvent.getByToken(tokenGctEtHad_,  etHad);
   edm::Handle<L1GctEtMissCollection> etMiss;
   iEvent.getByToken(tokenGctEtMiss_,  etMiss);
   edm::Handle<L1GctHFRingEtSumsCollection> hfRingSums;
   iEvent.getByToken(tokenGctHFRingEtSums_,  hfRingSums);
   edm::Handle<L1GctHFBitCountsCollection> hfBitCounts;
   iEvent.getByToken(tokenGctHFBitCounts_,  hfBitCounts);
   edm::Handle<L1GctHtMissCollection> htMiss;
   iEvent.getByToken(tokenGctHtMiss_,  htMiss);
   edm::Handle<L1GctJetCountsCollection> jetCounts;
   iEvent.getByToken(tokenGctJetCounts_, jetCounts);
 
   // get RCT EM Cand digi
   bool packRctEmThisEvent = packRctEm_;
   edm::Handle<L1CaloEmCollection> rctEm;
   if(packRctEmThisEvent)
   {
     iEvent.getByToken(tokenCaloEm_, rctEm);
     if(rctEm.failedToGet())
     {
       packRctEmThisEvent = false;
       LogDebug("GCT") << "RCT EM Candidate packing requested, but failed to get them from event!";
     }
   }
 
   // get RCT Calo region digi
   bool packRctCaloThisEvent = packRctCalo_;
   edm::Handle<L1CaloRegionCollection> rctCalo;
   if(packRctCaloThisEvent)
   {
     iEvent.getByToken(tokenCaloRegion_, rctCalo);
     if(rctCalo.failedToGet())
     {
       packRctCaloThisEvent = false;
       LogDebug("GCT") << "RCT Calo Region packing requested, but failed to get them from event!";
     }
   }
   
   // create the raw data collection
   std::unique_ptr<FEDRawDataCollection> rawColl(new FEDRawDataCollection()); 
  
   // get the GCT buffer
   FEDRawData& fedRawData=rawColl->FEDData(fedId_);
  
   // set the size & make pointers to the header, beginning of payload, and footer.
   unsigned int rawSize = 88;  // MUST BE MULTIPLE OF 8! (slink packets are 64 bit, but using 8-bit data struct).
   if(packRctEmThisEvent) { rawSize += 232; }  // Space for RCT EM Cands.
   if(packRctCaloThisEvent) { rawSize += 800; }  // Space for RCT Calo Regions (plus a 32-bit word of padding to make divisible by 8)
   fedRawData.resize(rawSize);
   unsigned char * pHeader = fedRawData.data();  
   unsigned char * pPayload = pHeader + 16;  //  16 = 8 for slink header + 8 for Greg's versioning header.
   unsigned char * pFooter = pHeader + rawSize - 8;
  
   // Write CDF header (exactly as told by Marco Zanetti)
   FEDHeader fedHeader(pHeader);
   fedHeader.set(pHeader, 1, eventNumber, bx, fedId_);  // what should the bx_ID be?
  
   // Pack GCT jet output digis
   formatTranslator_.writeGctOutJetBlock(pPayload, 
                                         cenJets.product(),
                                         forJets.product(),
                                         tauJets.product(),
                                         hfRingSums.product(), 
                                         hfBitCounts.product(),
                                         htMiss.product());
 
   pPayload += 36; //advance payload pointer
   
   // Pack GCT EM and energy sums digis.
   formatTranslator_.writeGctOutEmAndEnergyBlock(pPayload,
                                                 isoEm.product(), 
                                                 nonIsoEm.product(),
                                                 etTotal.product(), 
                                                 etHad.product(), 
                                                 etMiss.product());
 
   pPayload += 28; //advance payload pointer
 
   // Pack RCT EM Cands
   if(packRctEmThisEvent)
   {
     formatTranslator_.writeRctEmCandBlocks(pPayload, rctEm.product());
     pPayload+=232;  //advance payload pointer
   }
 
   // Pack RCT Calo Regions
   if(packRctCaloThisEvent)
   {
     formatTranslator_.writeAllRctCaloRegionBlock(pPayload, rctCalo.product());
   }
   
   // Write CDF footer (exactly as told by Marco Zanetti)
   FEDTrailer fedTrailer(pFooter);
   fedTrailer.set(pFooter, rawSize/8, evf::compute_crc(pHeader, rawSize), 0, 0);
  
   // Debug output.
   if (verbose_) { print(fedRawData); }
  
   // Put the collection in the event.
   iEvent.put(std::move(rawColl));
 }