#include <EcalDCC07UnpackingModule.h>

Inheritance diagram for EcalDCCTB07UnpackingModule:

Public Member Functions
void	beginJob ()

	EcalDCCTB07UnpackingModule (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob (void)

void	produce (edm::Event &e, const edm::EventSetup &c)
	Produce digis out of raw data. More...

virtual	~EcalDCCTB07UnpackingModule ()
	Destructor. More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

	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 &	itemsToGetFromEvent () 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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
CamacTBDataFormatter *	camacTBformatter_

EcalSupervisorTBDataFormatter *	ecalSupervisorFormatter_

edm::InputTag	fedRawDataCollectionTag_

EcalTB07DaqFormatter *	formatter_

MatacqTBDataFormatter *	matacqFormatter_

bool	ProduceEBDigis_

bool	ProduceEEDigis_

TableDataFormatter *	tableFormatter_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
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)

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

Definition at line 26 of file EcalDCC07UnpackingModule.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 41 of file EcalDCC07UnpackingModule.cc.

References camacTBformatter_, h2_mapping_cfi::ccuIDs, h2_mapping_cfi::channelIDs, ecalSupervisorFormatter_, edm::false, formatter_, edm::ParameterSet::getUntrackedParameter(), i, h2_mapping_cfi::ics, matacqFormatter_, h2_mapping_cfi::positionIDs, ProduceEBDigis_, ProduceEEDigis_, h2_mapping_cfi::statusIDs, AlCaHLTBitMon_QueryRunRegistry::string, h2_mapping_cfi::stripIDs, tableFormatter_, h2_mapping_cfi::towerIDs, and funct::true.

                                                                                   :
   fedRawDataCollectionTag_(pset.getParameter<edm::InputTag>("fedRawDataCollectionTag")) {
 
   std::string tbName = pset.getUntrackedParameter<std::string >("tbName", std::string("h2") );
 
   ProduceEEDigis_ = pset.getUntrackedParameter<bool >("produceEEdigi", true );
   ProduceEBDigis_ = pset.getUntrackedParameter<bool >("produceEBdigi", false );
 
   // index of crystal <-> tower ID (DQM plots) position <-> stripIDs <-> channelIDs for the test beam (2007)
   std::vector<int> ics        = pset.getUntrackedParameter<std::vector<int> >("ics",          std::vector<int>());
   std::vector<int> towerIDs   = pset.getUntrackedParameter<std::vector<int> >("towerIDs",     std::vector<int>());
   std::vector<int> stripIDs   = pset.getUntrackedParameter<std::vector<int> >("stripIDs",     std::vector<int>());
   std::vector<int> channelIDs = pset.getUntrackedParameter<std::vector<int> >("channelIDs",   std::vector<int>());
 
   // status id <-> tower CCU ID <-> DQM plots position mapping for the test beam (2007)
   std::vector<int> statusIDs   = pset.getUntrackedParameter<std::vector<int> >("statusIDs",   std::vector<int>());
   std::vector<int> ccuIDs      = pset.getUntrackedParameter<std::vector<int> >("ccuIDs",      std::vector<int>());
   std::vector<int> positionIDs = pset.getUntrackedParameter<std::vector<int> >("positionIDs", std::vector<int>());
 
   // check if vectors are filled
   if ( ics.size() == 0 || towerIDs.size() == 0 || stripIDs.size() == 0 || channelIDs.size() == 0 ){
     edm::LogError("EcalDCCTB07UnpackingModule") << "Some of the mapping info is missing! Check config files! " <<
       " Size of IC vector is " << ics.size() <<
       " Size of Tower ID vector is " << towerIDs.size() <<
       " Size of Strip ID vector is " << stripIDs.size() <<
       " Size of Channel ID vector is " << channelIDs.size();
   }
   if ( statusIDs.size() == 0 || ccuIDs.size() == 0 || positionIDs.size() == 0 ) {
     edm::LogError("EcalDCCTB07UnpackingModule") << "Some of the mapping info is missing! Check config files! " <<
       " Size of status ID vector is " << statusIDs.size() <<
       " Size of ccu ID vector is " << ccuIDs.size() <<
       " positionIDs size is " << positionIDs.size();
   }
   
   // check if vectors have the same size
   if ( ics.size() != towerIDs.size() || ics.size() != stripIDs.size() || ics.size() != channelIDs.size() ||
        towerIDs.size() != stripIDs.size() || towerIDs.size() != channelIDs.size() ||
        stripIDs.size() != channelIDs.size() )
     edm::LogError("EcalDCCTB07UnpackingModule") << "Mapping information is corrupted. " <<
       "Tower/DQM position/strip/channel vectors are of different size! Check cfi files! \n" <<
       " Size of IC vector is " << ics.size() <<
       " Size of Tower ID vector is " << towerIDs.size() <<
       " Size of Strip ID vector is " << stripIDs.size() <<
       " Size of Channel ID vector is " << channelIDs.size();
   
   if ( statusIDs.size() != ccuIDs.size() || statusIDs.size() != positionIDs.size() ||
        ccuIDs.size() != positionIDs.size() ) 
     edm::LogError("EcalDCCTB07UnpackingModule") << "Mapping information is corrupted. " <<
       "Status/CCU ID/DQM position vectors are of different size! Check cfi files! \n" <<
       " Size of status ID vector is " << statusIDs.size() <<
       " Size of ccu ID vector is " << ccuIDs.size() <<
       " positionIDs size is " << positionIDs.size();
   
   int cryIcMap[68][5][5];
   int tbStatusToLocation[71];
   int tbTowerIDToLocation[201];
   for (unsigned it=1; it <= 68; ++it )
     for (unsigned is=1; is <=5; ++is )
       for (unsigned ic=1; ic <=5; ++ic)
     cryIcMap[it-1][is-1][ic-1] = 1700;
 
   for (unsigned it=1; it <=71; ++it)
     tbStatusToLocation[it-1] = it - 1;
 
   for (unsigned it=1; it <= 201; ++it)
     tbTowerIDToLocation[it-1] = it - 1;
 
   // Fill the cry IC map
   for(unsigned int i=0; i < ics.size(); ++i) {
     int tower = towerIDs[i];
     int strip = stripIDs[i];
     int channel = channelIDs[i];
     int ic = ics[i];
     cryIcMap[tower-1][strip-1][channel-1] = ic;
   }
   for(unsigned int i = 0; i < statusIDs.size(); ++i) {
     int is = statusIDs[i];
     int it = ccuIDs[i];
     int itEB = positionIDs[i];
 
     tbStatusToLocation[is] = itEB;
     tbTowerIDToLocation[it] = itEB;
   }
 
   formatter_ = new EcalTB07DaqFormatter(tbName, cryIcMap, tbStatusToLocation, tbTowerIDToLocation);
   ecalSupervisorFormatter_ = new EcalSupervisorTBDataFormatter();
   camacTBformatter_ = new CamacTBDataFormatter();
   tableFormatter_ = new TableDataFormatter();
   matacqFormatter_ = new MatacqTBDataFormatter();
 
 
   // digis
   produces<EBDigiCollection>("ebDigis");
   produces<EEDigiCollection>("eeDigis");
   produces<EcalMatacqDigiCollection>();
   produces<EcalPnDiodeDigiCollection>();
   produces<EcalRawDataCollection>();
   produces<EcalTrigPrimDigiCollection>("EBTT");
 
   //TB specifics data
   produces<EcalTBHodoscopeRawInfo>();
   produces<EcalTBTDCRawInfo>();
   produces<EcalTBEventHeader>();
 
   // crystals' integrity
   produces<EBDetIdCollection>("EcalIntegrityDCCSizeErrors");
   produces<EcalElectronicsIdCollection>("EcalIntegrityTTIdErrors");
   produces<EcalElectronicsIdCollection>("EcalIntegrityBlockSizeErrors");
   produces<EBDetIdCollection>("EcalIntegrityChIdErrors");
   produces<EBDetIdCollection>("EcalIntegrityGainErrors");
   produces<EBDetIdCollection>("EcalIntegrityGainSwitchErrors");
 
   // mem channels' integrity
   produces<EcalElectronicsIdCollection>("EcalIntegrityMemTtIdErrors");
   produces<EcalElectronicsIdCollection>("EcalIntegrityMemBlockSize");
   produces<EcalElectronicsIdCollection>("EcalIntegrityMemChIdErrors");
   produces<EcalElectronicsIdCollection>("EcalIntegrityMemGainErrors");
 }

EcalDCCTB07UnpackingModule::~EcalDCCTB07UnpackingModule ( )

virtual

Destructor.

Definition at line 161 of file EcalDCC07UnpackingModule.cc.

References formatter_.

                                                        {
 
   delete formatter_;
 
 }

Member Function Documentation

void EcalDCCTB07UnpackingModule::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 167 of file EcalDCC07UnpackingModule.cc.

                                          {
 
 }

void EcalDCCTB07UnpackingModule::endJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 171 of file EcalDCC07UnpackingModule.cc.

                                        {
 
 }

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

virtual

Produce digis out of raw data.

Implements edm::EDProducer.

Definition at line 175 of file EcalDCC07UnpackingModule.cc.

References EcalDCCHeaderBlock::BEAMH4, BEG_DCC_FED_ID, BEG_DCC_FED_ID_GLOBAL, camacTBformatter_, EcalDCCHeaderBlock::COSMIC, gather_cfg::cout, data, ECAL_SUPERVISOR_FED_ID, ecalSupervisorFormatter_, END_DCC_FED_ID, END_DCC_FED_ID_GLOBAL, fedRawDataCollectionTag_, formatter_, edm::Event::getByLabel(), EcalSupervisorTBDataFormatter::interpretRawData(), TableDataFormatter::interpretRawData(), MatacqTBDataFormatter::interpretRawData(), EcalTB07DaqFormatter::interpretRawData(), CamacTBDataFormatter::interpretRawData(), LogDebug, MATACQ_FED_ID, matacqFormatter_, FEDNumbering::MAXFEDID, eostools::move(), ProduceEBDigis_, ProduceEEDigis_, edm::Event::put(), lumiPlot::rawdata, runPedHist::runType, FEDRawData::size(), TABLE_FED_ID, tableFormatter_, TBCAMAC_FED_ID, ECALTBParserException::what(), ECALTBParserBlockException::what(), and cms::Exception::what().

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

                                                                             {
 
   edm::Handle<FEDRawDataCollection> rawdata;
   e.getByLabel(fedRawDataCollectionTag_, rawdata);
   
 
   // create the collection of Ecal Digis
   auto productEb = std::make_unique<EBDigiCollection>();
 
   // YM create the collection of Ecal Endcap Digis
   auto productEe = std::make_unique<EEDigiCollection>();
 
   // create the collection of Matacq Digi
   auto productMatacq = std::make_unique<EcalMatacqDigiCollection>();
 
   // create the collection of Ecal PN's
   auto productPN = std::make_unique<EcalPnDiodeDigiCollection>();
   
   //create the collection of Ecal DCC Header
   auto productDCCHeader = std::make_unique<EcalRawDataCollection>();
 
   // create the collection with trigger primitives, bits and flags
   auto productTriggerPrimitives = std::make_unique<EcalTrigPrimDigiCollection>();
 
   // create the collection of Ecal Integrity DCC Size
   auto productDCCSize = std::make_unique<EBDetIdCollection>();
 
   // create the collection of Ecal Integrity TT Id
   auto productTTId = std::make_unique<EcalElectronicsIdCollection>();
 
   // create the collection of Ecal Integrity TT Block Size
   auto productBlockSize = std::make_unique<EcalElectronicsIdCollection>();
 
   // create the collection of Ecal Integrity Ch Id
   auto productChId = std::make_unique<EBDetIdCollection>();
 
   // create the collection of Ecal Integrity Gain
   auto productGain = std::make_unique<EBDetIdCollection>();
 
   // create the collection of Ecal Integrity Gain Switch
   auto productGainSwitch = std::make_unique<EBDetIdCollection>();
 
   // create the collection of Ecal Integrity Mem towerBlock_id errors
   auto productMemTtId = std::make_unique<EcalElectronicsIdCollection>();
   
   // create the collection of Ecal Integrity Mem gain errors
   auto productMemBlockSize = std::make_unique< EcalElectronicsIdCollection>();
 
   // create the collection of Ecal Integrity Mem gain errors
   auto productMemGain = std::make_unique< EcalElectronicsIdCollection>();
   
   // create the collection of Ecal Integrity Mem ch_id errors
   auto productMemChIdErrors = std::make_unique<EcalElectronicsIdCollection>();
   
   // create the collection of TB specifics data
   auto productHodo = std::make_unique<EcalTBHodoscopeRawInfo>();         
   auto productTdc = std::make_unique<EcalTBTDCRawInfo>();                      
   auto productHeader = std::make_unique<EcalTBEventHeader>();                      
 
 
   try {
 
   for (int id= 0; id<=FEDNumbering::MAXFEDID; ++id){ 
 
     //    edm::LogInfo("EcalDCCTB07UnpackingModule") << "EcalDCCTB07UnpackingModule::Got FED ID "<< id <<" ";
     const FEDRawData& data = rawdata->FEDData(id);
     //    edm::LogInfo("EcalDCCTB07UnpackingModule") << " Fed data size " << data.size() ;
    
     //std::cout <<"1 Fed id: "<<dec<<id<< " Fed data size: " <<data.size() << std::endl;
 //    const unsigned char * pData = data.data();
 //    int length = data.size();
 //    if(length >0 ){
 //      if(length >= 40){length = 40;}
 //    std::cout<<"##############################################################"<<std::endl;
 //    for( int i=0; i<length; i++ ) {
 //      std::cout << std::hex << std::setw(8) << int(pData[i]) << " ";
 //      if( (i+1)%8 == 0 ) std::cout << std::endl;
 //     }
 //    std::cout<<"##############################################################"<<std::endl;
 //    } 
     if (data.size()>16){
 
       if (    (id >= BEG_DCC_FED_ID && id <= END_DCC_FED_ID) ||
       ( BEG_DCC_FED_ID_GLOBAL <= id &&  id <= END_DCC_FED_ID_GLOBAL )
       )
     {   // do the DCC data unpacking and fill the collections
       
       (*productHeader).setSmInBeam(id);
       // YM add productEe to the list of arguments of the formatter
       formatter_->interpretRawData(data,  *productEb, *productEe, *productPN, 
                        *productDCCHeader, 
                        *productDCCSize, 
                        *productTTId, *productBlockSize, 
                        *productChId, *productGain, *productGainSwitch, 
                        *productMemTtId,  *productMemBlockSize,
                        *productMemGain,  *productMemChIdErrors,
                        *productTriggerPrimitives);
       int runType = (*productDCCHeader)[0].getRunType();
       if ( runType == EcalDCCHeaderBlock::COSMIC || runType == EcalDCCHeaderBlock::BEAMH4 ) 
         (*productHeader).setTriggerMask(0x1);
       else if ( runType == 4 || runType == 5 || runType == 6 ) //laser runs
         (*productHeader).setTriggerMask(0x2000);
       else if ( runType == 9 || runType == 10 || runType == 11 ) //pedestal runs
         (*productHeader).setTriggerMask(0x800);
       LogDebug("EcalDCCTB07UnpackingModule") << "Event type is " << (*productHeader).eventType() << " dbEventType " << (*productHeader).dbEventType();
     } 
       else if ( id == ECAL_SUPERVISOR_FED_ID )
     ecalSupervisorFormatter_->interpretRawData(data, *productHeader);
       else if ( id == TBCAMAC_FED_ID )
     camacTBformatter_->interpretRawData(data, *productHeader,*productHodo, *productTdc );
       else if ( id == TABLE_FED_ID )
     tableFormatter_->interpretRawData(data, *productHeader);
       else if ( id == MATACQ_FED_ID )     
     matacqFormatter_->interpretRawData(data, *productMatacq);
     }// endif 
   }//endfor
   
 
   // commit to the event  
   e.put(std::move(productPN));
   if (ProduceEBDigis_)  e.put(std::move(productEb),"ebDigis");
   if (ProduceEEDigis_)  e.put(std::move(productEe),"eeDigis");
   e.put(std::move(productMatacq));
   e.put(std::move(productDCCHeader));
   e.put(std::move(productTriggerPrimitives), "EBTT");
   
   if (ProduceEBDigis_)  e.put(std::move(productDCCSize),"EcalIntegrityDCCSizeErrors");
   if (ProduceEBDigis_)  e.put(std::move(productTTId),"EcalIntegrityTTIdErrors");
   if (ProduceEBDigis_)  e.put(std::move(productBlockSize),"EcalIntegrityBlockSizeErrors");
   if (ProduceEBDigis_)  e.put(std::move(productChId),"EcalIntegrityChIdErrors");
   if (ProduceEBDigis_)  e.put(std::move(productGain),"EcalIntegrityGainErrors");
   if (ProduceEBDigis_)  e.put(std::move(productGainSwitch),"EcalIntegrityGainSwitchErrors");
   
   if (ProduceEBDigis_)  e.put(std::move(productMemTtId),"EcalIntegrityMemTtIdErrors");
   if (ProduceEBDigis_)  e.put(std::move(productMemBlockSize),"EcalIntegrityMemBlockSize");
   if (ProduceEBDigis_)  e.put(std::move(productMemChIdErrors),"EcalIntegrityMemChIdErrors");
   if (ProduceEBDigis_)  e.put(std::move(productMemGain),"EcalIntegrityMemGainErrors");
 
   e.put(std::move(productHodo));
   e.put(std::move(productTdc));
   e.put(std::move(productHeader));
 
   } catch (ECALTBParserException &e) {
     std::cout << "[EcalDCCTB07UnpackingModule] " << e.what() << std::endl;
   } catch (ECALTBParserBlockException &e) {
     std::cout << "[EcalDCCTB07UnpackingModule] " << e.what() << std::endl;
   } catch (cms::Exception &e) {
     std::cout << "[EcalDCCTB07UnpackingModule] " << e.what() << std::endl;
   } catch (...) {
     std::cout << "[EcalDCCTB07UnpackingModule] Unknown exception ..." << std::endl;
   }
 
 }