|
|
|
#include <EcalDCC07UnpackingModule.h>
Public Member Functions | |
| void | beginJob () |
| EcalDCCTB07UnpackingModule (const edm::ParameterSet &pset) | |
| Constructor. | |
| void | endJob (void) |
| void | produce (edm::Event &e, const edm::EventSetup &c) |
| Produce digis out of raw data. | |
| virtual | ~EcalDCCTB07UnpackingModule () |
| Destructor. | |
Private Attributes | |
| CamacTBDataFormatter * | camacTBformatter_ |
| EcalSupervisorTBDataFormatter * | ecalSupervisorFormatter_ |
| EcalTB07DaqFormatter * | formatter_ |
| MatacqTBDataFormatter * | matacqFormatter_ |
| bool | ProduceEBDigis_ |
| bool | ProduceEEDigis_ |
| TableDataFormatter * | tableFormatter_ |
Definition at line 27 of file EcalDCC07UnpackingModule.h.
| 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_;
}
| void EcalDCCTB07UnpackingModule::beginJob | ( | void | ) | [virtual] |
| void EcalDCCTB07UnpackingModule::endJob | ( | void | ) | [virtual] |
| 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_, ExpressReco_HICollisions_FallBack::EcalTBEventHeader, ExpressReco_HICollisions_FallBack::EcalTBHodoscopeRawInfo, ExpressReco_HICollisions_FallBack::EcalTBTDCRawInfo, END_DCC_FED_ID, END_DCC_FED_ID_GLOBAL, formatter_, edm::Event::getByType(), ExpressReco_HICollisions_FallBack::id, MatacqTBDataFormatter::interpretRawData(), EcalTB07DaqFormatter::interpretRawData(), TableDataFormatter::interpretRawData(), EcalSupervisorTBDataFormatter::interpretRawData(), CamacTBDataFormatter::interpretRawData(), LogDebug, MATACQ_FED_ID, matacqFormatter_, FEDNumbering::MAXFEDID, ProduceEBDigis_, ProduceEEDigis_, edm::Event::put(), ExpressReco_HICollisions_FallBack::runType, FEDRawData::size(), TABLE_FED_ID, tableFormatter_, TBCAMAC_FED_ID, ECALTBParserException::what(), cms::Exception::what(), and ECALTBParserBlockException::what().
{
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;
}
}
Definition at line 48 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), and produce().
Definition at line 47 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), and produce().
Definition at line 46 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), produce(), and ~EcalDCCTB07UnpackingModule().
Definition at line 50 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), and produce().
bool EcalDCCTB07UnpackingModule::ProduceEBDigis_ [private] |
Definition at line 53 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), and produce().
bool EcalDCCTB07UnpackingModule::ProduceEEDigis_ [private] |
Definition at line 52 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), and produce().
Definition at line 49 of file EcalDCC07UnpackingModule.h.
Referenced by EcalDCCTB07UnpackingModule(), and produce().
1.7.3