#include <GEMRawToDigiModule.h>

Inheritance diagram for GEMRawToDigiModule:

Public Member Functions
	GEMRawToDigiModule (const edm::ParameterSet &pset)
	Constructor. More...

std::shared_ptr< GEMROmap >	globalBeginRun (edm::Run const &, edm::EventSetup const &) const override

void	globalEndRun (edm::Run const &, edm::EventSetup const &) const override

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

Public Member Functions inherited from edm::global::EDProducer< edm::RunCache< GEMROmap > >
	EDProducer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

ModuleDescription const &	moduleDescription () 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)

virtual	~ProducerBase () noexcept(false)

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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::EDGetTokenT< FEDRawDataCollection >	fed_token

bool	unPackStatusDigis_

bool	useDBEMap_

Additional Inherited Members
Public Types inherited from edm::global::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

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

on CSCDigiToRawModule

Author: J. Lee - UoS

Definition at line 32 of file GEMRawToDigiModule.h.

Constructor & Destructor Documentation

GEMRawToDigiModule::GEMRawToDigiModule ( const edm::ParameterSet & pset )

Constructor.

Definition at line 19 of file GEMRawToDigiModule.cc.

References unPackStatusDigis_.

                                                                    :
   fed_token(consumes<FEDRawDataCollection>( pset.getParameter<edm::InputTag>("InputLabel") )),
   useDBEMap_(pset.getParameter<bool>("useDBEMap")),
   unPackStatusDigis_(pset.getParameter<bool>("unPackStatusDigis"))
 {
   produces<GEMDigiCollection>(); 
   if (unPackStatusDigis_){
     produces<GEMVfatStatusDigiCollection>("vfatStatus");
     produces<GEMGEBStatusDigiCollection>("GEBStatus");
     produces<GEMAMCStatusDigiCollection>("AMCStatus"); 
   }
 }

Member Function Documentation

void GEMRawToDigiModule::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 32 of file GEMRawToDigiModule.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("InputLabel", edm::InputTag("rawDataCollector")); 
   desc.add<bool>("useDBEMap", false);
   desc.add<bool>("unPackStatusDigis", false);
   descriptions.add("muonGEMDigis", desc);  
 }

std::shared_ptr< GEMROmap > GEMRawToDigiModule::globalBeginRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)		const

override

Definition at line 41 of file GEMRawToDigiModule.cc.

References edm::EventSetup::get(), edm::ESHandle< T >::product(), and useDBEMap_.

 {
   auto gemORmap = std::make_shared<GEMROmap>();
   if (useDBEMap_){
     edm::ESHandle<GEMELMap> gemEMapRcd;
     iSetup.get<GEMELMapRcd>().get(gemEMapRcd);
     auto gemEMap = std::make_unique<GEMELMap>(*(gemEMapRcd.product()));
     gemEMap->convert(*gemORmap);
     gemEMap.reset();    
   }
   else {
     // no EMap in DB, using dummy
     auto gemEMap = std::make_unique<GEMELMap>();
     gemEMap->convertDummy(*gemORmap);
     gemEMap.reset();    
   }
   return gemORmap;
 }

void GEMRawToDigiModule::globalEndRun	(	edm::Run const &	,
		edm::EventSetup const &
	)		const

inlineoverride

Definition at line 40 of file GEMRawToDigiModule.h.

40 {};

void GEMRawToDigiModule::produce	(	edm::StreamID	iID,
		edm::Event &	iEvent,
		edm::EventSetup const &
	)		const

overridevirtual

Implements edm::global::EDProducerBase.

Definition at line 60 of file GEMRawToDigiModule.cc.

References GEMDigi::bx(), officialStyle::chan, GEMROmap::channelNum::chNum, FEDRawData::data(), data, RecoEcal_EventContent_cff::ec, fed_token, GEMROmap::dCoord::gemDetId, edm::Event::getByToken(), edm::Event::getRun(), mps_fire::i, triggerObjects_cff::id, edm::Run::index(), createfilelist::int, GEMROmap::dCoord::iPhi, gen::k, LogDebug, GEMELMap::maxChan_, GEMELMap::maxVFatGE11_, GEMELMap::maxVFatGE21_, FEDNumbering::MINGEMFEDID, eostools::move(), ecaldqm::nChannels, edm::Event::put(), FEDRawData::size(), GEMDetId::station(), GEMROmap::stripNum::stNum, unPackStatusDigis_, GEMROmap::eCoord::vfatId, GEMROmap::dCoord::vfatType, and globals_cff::x1.

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

 {
   auto outGEMDigis = std::make_unique<GEMDigiCollection>();
   auto outVfatStatus = std::make_unique<GEMVfatStatusDigiCollection>();
   auto outGEBStatus = std::make_unique<GEMGEBStatusDigiCollection>();
   auto outAMCStatus = std::make_unique<GEMAMCStatusDigiCollection>();
 
  // Take raw from the event
   edm::Handle<FEDRawDataCollection> fed_buffers;
   iEvent.getByToken( fed_token, fed_buffers );
   
   auto gemROMap = runCache(iEvent.getRun().index());
   
   for (unsigned int id=FEDNumbering::MINGEMFEDID; id<=FEDNumbering::MINGEMFEDID; ++id){ 
     const FEDRawData& fedData = fed_buffers->FEDData(id);
     
     int nWords = fedData.size()/sizeof(uint64_t);
     LogDebug("GEMRawToDigiModule") <<" words " << nWords;
     if (nWords<5) continue;
     const unsigned char * data = fedData.data();
     
     auto amc13Event = std::make_unique<AMC13Event>();
     
     const uint64_t* word = reinterpret_cast<const uint64_t* >(data);
     
     amc13Event->setCDFHeader(*word);
     amc13Event->setAMC13header(*(++word));
 
     // Readout out AMC headers
     for (unsigned short i = 0; i < amc13Event->nAMC(); ++i)
       amc13Event->addAMCheader(*(++word));
     
     // Readout out AMC payloads
     for (unsigned short i = 0; i < amc13Event->nAMC(); ++i){
       auto amcData = std::make_unique<AMCdata>();
       amcData->setAMCheader1(*(++word));      
       amcData->setAMCheader2(*(++word));
       amcData->setGEMeventHeader(*(++word));
       uint16_t amcId = amcData->boardId();
       uint16_t amcBx = amcData->bx();
 
       // Fill GEB
       for (unsigned short j = 0; j < amcData->gdCount(); ++j){
     auto gebData = std::make_unique<GEBdata>();
     gebData->setChamberHeader(*(++word));
     
     unsigned int m_nvb = gebData->vwh() / 3; // number of VFAT2 blocks
     uint16_t gebId = gebData->inputID();
     GEMDetId gemId(-1,1,1,1,1,0); // temp ID
     for (unsigned short k = 0; k < m_nvb; k++){
       auto vfatData = std::make_unique<VFATdata>();
       vfatData->read_fw(*(++word));
       vfatData->read_sw(*(++word));
       vfatData->read_tw(*(++word));
       gebData->addVFAT(*vfatData);
       
       uint16_t bc=vfatData->bc();
       uint8_t b1010=vfatData->b1010();
       uint8_t b1100=vfatData->b1100();
       uint8_t b1110=vfatData->b1110();
       uint16_t vfatId=vfatData->chipID();
       uint16_t crc = vfatData->crc();
       uint16_t crc_check = vfatData->checkCRC();
       bool Quality = (b1010==10) && (b1100==12) && (b1110==14) && (crc==crc_check);
 
       if (crc!=crc_check) edm::LogWarning("GEMRawToDigiModule") << "DIFFERENT CRC :"<<crc<<"   "<<crc_check;
       if (!Quality) edm::LogWarning("GEMRawToDigiModule") << "Quality "<< Quality
                                   << " b1010 "<< int(b1010)
                                   << " b1100 "<< int(b1100)
                                   << " b1110 "<< int(b1110);
       
       //check if ChipID exists.
       GEMROmap::eCoord ec = {amcId, gebId, vfatId};
 
       if (!gemROMap->isValidChipID(ec)){
         edm::LogWarning("GEMRawToDigiModule") << "InValid ChipID :"<<ec.vfatId;
         continue;
       }
       
       for (int chan = 0; chan < VFATdata::nChannels; ++chan) {
         uint8_t chan0xf = 0;
         if (chan < 64) chan0xf = ((vfatData->lsData() >> chan) & 0x1);
         else chan0xf = ((vfatData->msData() >> (chan-64)) & 0x1);
 
         // no hits
         if(chan0xf==0) continue;
         GEMROmap::dCoord dc = gemROMap->hitPosition(ec);
         // strip bx = vfat bx - amc bx
         int bx = bc-amcBx;
         gemId = dc.gemDetId;
          
             GEMROmap::channelNum chMap = {dc.vfatType, chan};
             GEMROmap::stripNum stMap = gemROMap->hitPosition(chMap);
 
             int maxVFat = GEMELMap::maxVFatGE11_;
             if (gemId.station() == 2) maxVFat = GEMELMap::maxVFatGE21_;
             int stripId = stMap.stNum + (dc.iPhi-1)%maxVFat*GEMELMap::maxChan_;    
 
         GEMDigi digi(stripId,bx);
         LogDebug("GEMRawToDigiModule") <<" vfatId "<<ec.vfatId
                        <<" gemDetId "<< gemId
                        <<" chan "<< chMap.chNum
                        <<" strip "<< stripId
                        <<" bx "<< digi.bx();
 
         outGEMDigis.get()->insertDigi(gemId,digi);      
       }
 
           if (unPackStatusDigis_){
         GEMVfatStatusDigi vfatStatus(vfatData->lsData(), vfatData->msData(),
                      crc, vfatData->crc_calc(),
                      b1010, b1100, b1110, vfatData->flag(),
                      vfatData->isBlockGood());
             outVfatStatus.get()->insertDigi(gebData->inputID(),vfatStatus);
       }
       
     }
     
     gebData->setChamberTrailer(*(++word));
         if (unPackStatusDigis_){
           GEMGEBStatusDigi gebStatus(gebData->zeroSup(),
                                      gebData->vwh(),
                                      gebData->errorC(),
                                      gebData->ohCRC(),
                                      gebData->vwt(),
                                      gebData->inFu(),
                                      gebData->inputID(),                     
                                      gebData->stuckd(),
                                      gebData->getGEBflag());
           outGEBStatus.get()->insertDigi(gebData->inputID(),gebStatus); 
         }
             
     amcData->addGEB(*gebData);
       }
       
       amcData->setGEMeventTrailer(*(++word));
       amcData->setAMCTrailer(*(++word));
 
       if (unPackStatusDigis_){
         GEMAMCStatusDigi amcStatus(amcData->amcNum(),
                                    amcData->l1A(),
                                    amcData->bx(),
                                    amcData->dlength(),
                                    amcData->fv(),
                                    amcData->runType(),
                                    amcData->param1(),
                                    amcData->param2(),
                                    amcData->param3(),
                                    amcData->orbitNum(),
                                    amcData->boardId(),
                                    amcData->gemDAV(),
                                    amcData->bStatus(),
                                    amcData->gdCount(),
                                    amcData->ttsState(),
                                    amcData->chamberTimeOut(),
                                    amcData->oosGLIB());
         outAMCStatus.get()->insertDigi(amcData->boardId(), amcStatus);
       }
 
       amc13Event->addAMCpayload(*amcData);
     }
     
     amc13Event->setAMC13trailer(*(++word));
     amc13Event->setCDFTrailer(*(++word));
   }
   
   iEvent.put(std::move(outGEMDigis));
   if (unPackStatusDigis_){
     iEvent.put(std::move(outVfatStatus), "vfatStatus");
     iEvent.put(std::move(outGEBStatus), "GEBStatus");
     iEvent.put(std::move(outAMCStatus), "AMCStatus");
   }
 }