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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Attributes
CamacTBDataFormatter *	camacTBformatter_

EcalSupervisorTBDataFormatter *	ecalSupervisorFormatter_

EcalTB07DaqFormatter *	formatter_

MatacqTBDataFormatter *	matacqFormatter_

bool	ProduceEBDigis_

bool	ProduceEEDigis_

TableDataFormatter *	tableFormatter_

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

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 27 of file EcalDCC07UnpackingModule.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 43 of file EcalDCC07UnpackingModule.cc.

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

                                                                                  {
 
   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 162 of file EcalDCC07UnpackingModule.cc.

References formatter_.

                                                        {
 
   delete formatter_;
 
 }

Member Function Documentation

void EcalDCCTB07UnpackingModule::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 168 of file EcalDCC07UnpackingModule.cc.

                                          {
 
 }

void EcalDCCTB07UnpackingModule::endJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 172 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 176 of file EcalDCC07UnpackingModule.cc.

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

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

                                                                             {
 
   edm::Handle<FEDRawDataCollection> rawdata;
   e.getByType(rawdata);
   
 
   // create the collection of Ecal Digis
   std::auto_ptr<EBDigiCollection> productEb(new EBDigiCollection);
 
   // YM create the collection of Ecal Endcap Digis
   std::auto_ptr<EEDigiCollection> productEe(new EEDigiCollection);
 
   // create the collection of Matacq Digi
   std::auto_ptr<EcalMatacqDigiCollection> productMatacq(new EcalMatacqDigiCollection());
 
   // create the collection of Ecal PN's
   std::auto_ptr<EcalPnDiodeDigiCollection> productPN(new EcalPnDiodeDigiCollection);
   
   //create the collection of Ecal DCC Header
   std::auto_ptr<EcalRawDataCollection> productDCCHeader(new EcalRawDataCollection);
 
   // create the collection with trigger primitives, bits and flags
   std::auto_ptr<EcalTrigPrimDigiCollection> productTriggerPrimitives(new EcalTrigPrimDigiCollection);
 
   // create the collection of Ecal Integrity DCC Size
   std::auto_ptr<EBDetIdCollection> productDCCSize(new EBDetIdCollection);
 
   // create the collection of Ecal Integrity TT Id
   std::auto_ptr<EcalElectronicsIdCollection> productTTId(new EcalElectronicsIdCollection);
 
   // create the collection of Ecal Integrity TT Block Size
   std::auto_ptr<EcalElectronicsIdCollection> productBlockSize(new EcalElectronicsIdCollection);
 
   // create the collection of Ecal Integrity Ch Id
   std::auto_ptr<EBDetIdCollection> productChId(new EBDetIdCollection);
 
   // create the collection of Ecal Integrity Gain
   std::auto_ptr<EBDetIdCollection> productGain(new EBDetIdCollection);
 
   // create the collection of Ecal Integrity Gain Switch
   std::auto_ptr<EBDetIdCollection> productGainSwitch(new EBDetIdCollection);
 
   // create the collection of Ecal Integrity Mem towerBlock_id errors
   std::auto_ptr<EcalElectronicsIdCollection> productMemTtId(new EcalElectronicsIdCollection);
   
   // create the collection of Ecal Integrity Mem gain errors
   std::auto_ptr< EcalElectronicsIdCollection> productMemBlockSize(new EcalElectronicsIdCollection);
 
   // create the collection of Ecal Integrity Mem gain errors
   std::auto_ptr< EcalElectronicsIdCollection> productMemGain(new EcalElectronicsIdCollection);
   
   // create the collection of Ecal Integrity Mem ch_id errors
   std::auto_ptr<EcalElectronicsIdCollection> productMemChIdErrors(new EcalElectronicsIdCollection);
   
   // create the collection of TB specifics data
   std::auto_ptr<EcalTBHodoscopeRawInfo> productHodo(new EcalTBHodoscopeRawInfo());         
   std::auto_ptr<EcalTBTDCRawInfo> productTdc(new EcalTBTDCRawInfo());                      
   std::auto_ptr<EcalTBEventHeader> productHeader(new 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(productPN);
   if (ProduceEBDigis_)  e.put(productEb,"ebDigis");
   if (ProduceEEDigis_)  e.put(productEe,"eeDigis");
   e.put(productMatacq);
   e.put(productDCCHeader);
   e.put(productTriggerPrimitives, "EBTT");
   
   if (ProduceEBDigis_)  e.put(productDCCSize,"EcalIntegrityDCCSizeErrors");
   if (ProduceEBDigis_)  e.put(productTTId,"EcalIntegrityTTIdErrors");
   if (ProduceEBDigis_)  e.put(productBlockSize,"EcalIntegrityBlockSizeErrors");
   if (ProduceEBDigis_)  e.put(productChId,"EcalIntegrityChIdErrors");
   if (ProduceEBDigis_)  e.put(productGain,"EcalIntegrityGainErrors");
   if (ProduceEBDigis_)  e.put(productGainSwitch,"EcalIntegrityGainSwitchErrors");
   
   if (ProduceEBDigis_)  e.put(productMemTtId,"EcalIntegrityMemTtIdErrors");
   if (ProduceEBDigis_)  e.put(productMemBlockSize,"EcalIntegrityMemBlockSize");
   if (ProduceEBDigis_)  e.put(productMemChIdErrors,"EcalIntegrityMemChIdErrors");
   if (ProduceEBDigis_)  e.put(productMemGain,"EcalIntegrityMemGainErrors");
 
   e.put(productHodo);
   e.put(productTdc);
   e.put(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;
   }
 
 }