/** * CbmMQTsaMultiSampler.cpp * * @since 2017-11-17 * @author F. Uhlig */ #include "CbmMQTsaMultiSampler.h" #include "CbmMQDefs.h" #include "FairMQLogger.h" #include "FairMQProgOptions.h" // device->fConfig #include "TimesliceSubscriber.hpp" #include "TimesliceMultiSubscriber.hpp" #include "TimesliceInputArchive.hpp" #include "TimesliceMultiInputArchive.hpp" #include #include #include #include namespace filesys = boost::filesystem; #include #include #include // this_thread::sleep_for #include #include #include using namespace std; #include struct InitTaskError : std::runtime_error { using std::runtime_error::runtime_error; }; CbmMQTsaMultiSampler::CbmMQTsaMultiSampler() : FairMQDevice() , fMaxTimeslices(0) , fFileName("") , fDirName("") , fInputFileList() , fFileCounter(0) , fHost("") , fPort(0) , fHighWaterMark(1) , fTSNumber(0) , fTSCounter(0) , fMessageCounter(0) , fSource(nullptr) , fTime() { } void CbmMQTsaMultiSampler::InitTask() try { // Get the values from the command line options (via fConfig) fFileName = fConfig->GetValue("filename"); fDirName = fConfig->GetValue("dirname"); fHost = fConfig->GetValue("flib-host"); fPort = fConfig->GetValue("flib-port"); fHighWaterMark = fConfig->GetValue("high-water-mark"); fMaxTimeslices = fConfig->GetValue("max-timeslices"); if (0 == fMaxTimeslices) fMaxTimeslices = UINT_MAX; // Check which input is defined // Posibilities // filename && ! dirname : single file // filename with wildcards && diranme : all files with filename regex in the directory // host && port : connect to the flim server bool isGoodInputCombi{false}; if ( 0 != fFileName.size() && 0 == fDirName.size() && 0 == fHost.size() && 0 == fPort ) { isGoodInputCombi=true; fInputFileList.push_back(fFileName); } else if ( 0 != fFileName.size() && 0 != fDirName.size() && 0 == fHost.size() && 0 == fPort ) { isGoodInputCombi=true; fInputFileList.push_back(fFileName); } else if ( 0 == fFileName.size() && 0 == fDirName.size() && 0 != fHost.size() && 0 != fPort) { isGoodInputCombi=true; LOG(info) << "Host: " << fHost; LOG(info) << "Port: " << fPort; } else if ( 0 == fFileName.size() && 0 == fDirName.size() && 0 != fHost.size() && 0 == fPort) { isGoodInputCombi=true; LOG(info) << "Host string: " << fHost; } else { isGoodInputCombi=false; } if (!isGoodInputCombi) { throw InitTaskError("Wrong combination of inputs. Either file or wildcard file + directory or host + port are allowed combination."); } LOG(info) << "MaxTimeslices: " << fMaxTimeslices; // Get the information about created channels from the device // Check if the defined channels from the topology (by name) // are in the list of channels which are possible/allowed // for the device // The idea is to check at initilization if the devices are // properly connected. For the time beeing this is done with a // nameing convention. It is not avoided that someone sends other // data on this channel. int noChannel = fChannels.size(); LOG(info) << "Number of defined output channels: " << noChannel; for(auto const &entry : fChannels) { LOG(info) << "Channel name: " << entry.first; if (!IsChannelNameAllowed(entry.first)) throw InitTaskError("Channel name does not match."); } for(auto const& value: fComponentsToSend) { LOG(info) << "Value : " << value; if (value > 1) { throw InitTaskError("Sending same data to more than one output channel not implemented yet."); } } if ( 0 == fFileName.size() && 0 != fHost.size() && 0 != fPort ) { // Don't add the protocol since this is done now in the TimesliceMultiSubscriber //std::string connector = "tcp://" + fHost + ":" + std::to_string(fPort); std::string connector = fHost + ":" + std::to_string(fPort); LOG(info) << "Open TSPublisher at " << connector; fSource = new fles::TimesliceMultiSubscriber(connector); } else if ( 0 == fFileName.size() && 0 != fHost.size() ) { std::string connector = fHost; LOG(info) << "Open TSPublisher with host string: " << connector; fSource = new fles::TimesliceMultiSubscriber(connector, fHighWaterMark); } else { // Create a ";" separated string with all file names std::string fileList{""}; for(const auto obj: fInputFileList) { std::string fileName = obj; fileList += fileName; fileList += ";"; } fileList.pop_back(); // Remove the last ; LOG(info) << "Input File String: " << fileList; fSource = new fles::TimesliceMultiInputArchive(fileList, fDirName); if ( !fSource) { throw InitTaskError("Could open files from file list."); } } fTime = std::chrono::steady_clock::now(); } catch (InitTaskError& e) { LOG(error) << e.what(); // Wrapper defined in CbmMQDefs.h to support different FairMQ versions cbm::mq::ChangeState(this, cbm::mq::Transition::ErrorFound); } bool CbmMQTsaMultiSampler::IsChannelNameAllowed(std::string channelName) { for(auto const &entry : fAllowedChannels) { LOG(info) << "Looking for name " << channelName << " in " << entry; std::size_t pos1 = channelName.find(entry); if (pos1!=std::string::npos) { const vector::const_iterator pos = std::find(fAllowedChannels.begin(), fAllowedChannels.end(), entry); const vector::size_type idx = pos-fAllowedChannels.begin(); LOG(info) << "Found " << entry << " in " << channelName; LOG(info) << "Channel name " << channelName << " found in list of allowed channel names at position " << idx; fComponentsToSend[idx]++; fChannelsToSend[idx].push_back(channelName); return true; } } LOG(info) << "Channel name " << channelName << " not found in list of allowed channel names."; LOG(error) << "Stop device."; return false; } bool CbmMQTsaMultiSampler::ConditionalRun() { auto timeslice = fSource->get(); if (timeslice) { if (fTSCounter < fMaxTimeslices) { fTSCounter++; if (fTSCounter % 10000 == 0) { LOG(info) << "Analyse Event " << fTSCounter; } const fles::Timeslice& ts = *timeslice; // auto tsIndex = ts.index(); LOG(debug) << "Found " << ts.num_components() << " different components in timeslice"; // CheckTimeslice(ts); for (unsigned int nrComp = 0; nrComp < ts.num_components(); ++nrComp) { CreateAndSendComponent(ts, nrComp); } return true; } else { CalcRuntime(); return false; } } else { CalcRuntime(); return false; } } bool CbmMQTsaMultiSampler::CreateAndSendComponent(const fles::Timeslice& ts, int nrComp) { // Check if component has to be send. If the corresponding channel // is connected create the new timeslice and send it to the // correct channel LOG(debug) << "SysID: " << static_cast(ts.descriptor(nrComp,0).sys_id); const vector::const_iterator pos = std::find(fSysId.begin(), fSysId.end(), static_cast(ts.descriptor(nrComp,0).sys_id)); if (pos != fSysId.end() ) { const vector::size_type idx = pos-fSysId.begin(); if (fComponentsToSend[idx]>0) { LOG(debug) << "Create timeslice component for link " << nrComp; fles::StorableTimeslice component{static_cast(ts.num_core_microslices()), ts.index()}; component.append_component(ts.num_microslices(0)); for (size_t m = 0; m < ts.num_microslices(nrComp); ++m) { component.append_microslice( 0, m, ts.descriptor(nrComp, m), ts.content(nrComp, m) ); } /* LOG(info) << "Number of core microslices before: " << ts.num_core_microslices(); LOG(info) << "Number of core microslices after : " << component.num_core_microslices(); LOG(info) << "Number of microslices: " << component.num_microslices(0); */ if ( ! SendData(component, idx) ) return false; return true; } } return true; } bool CbmMQTsaMultiSampler::SendData(const fles::StorableTimeslice& component, int idx) { // serialize the timeslice and create the message std::stringstream oss; boost::archive::binary_oarchive oa(oss); oa << component; std::string* strMsg = new std::string(oss.str()); FairMQMessagePtr msg(NewMessage(const_cast(strMsg->c_str()), // data strMsg->length(), // size [](void* /*data*/, void* object){ delete static_cast(object); }, strMsg)); // object that manages the data // TODO: Implement sending same data to more than one channel // Need to create new message (copy message??) if (fComponentsToSend[idx]>1) { LOG(info) << "Need to copy FairMessage"; } // in case of error or transfer interruption, // return false to go to IDLE state // successfull transfer will return number of bytes // transfered (can be 0 if sending an empty message). LOG(debug) << "Send data to channel " << fChannelsToSend[idx][0]; if (Send(msg, fChannelsToSend[idx][0]) < 0) { LOG(error) << "Problem sending data"; return false; } fMessageCounter++; LOG(debug) << "Send message " << fMessageCounter << " with a size of " << msg->GetSize(); return true; } CbmMQTsaMultiSampler::~ CbmMQTsaMultiSampler() { } void CbmMQTsaMultiSampler::CalcRuntime() { std::chrono::duration run_time = std::chrono::steady_clock::now() - fTime; LOG(info) << "Runtime: " << run_time.count(); LOG(info) << "No more input data"; } void CbmMQTsaMultiSampler::PrintMicroSliceDescriptor(const fles::MicrosliceDescriptor& mdsc) { LOG(info) << "Header ID: Ox" << std::hex << static_cast(mdsc.hdr_id) << std::dec; LOG(info) << "Header version: Ox" << std::hex << static_cast(mdsc.hdr_ver) << std::dec; LOG(info) << "Equipement ID: " << mdsc.eq_id; LOG(info) << "Flags: " << mdsc.flags; LOG(info) << "Sys ID: Ox" << std::hex << static_cast(mdsc.sys_id) << std::dec; LOG(info) << "Sys version: Ox" << std::hex << static_cast(mdsc.sys_ver) << std::dec; LOG(info) << "Microslice Idx: " << mdsc.idx; LOG(info) << "Checksum: " << mdsc.crc; LOG(info) << "Size: " << mdsc.size; LOG(info) << "Offset: " << mdsc.offset; } bool CbmMQTsaMultiSampler::CheckTimeslice(const fles::Timeslice& ts) { if ( 0 == ts.num_components() ) { LOG(error) << "No Component in TS " << ts.index(); return 1; } LOG(info) << "Found " << ts.num_components() << " different components in timeslice"; for (size_t c = 0; c < ts.num_components(); ++c) { LOG(info) << "Found " << ts.num_microslices(c) << " microslices in component " << c; LOG(info) << "Component " << c << " has a size of " << ts.size_component(c) << " bytes"; LOG(info) << "Component " << c << " has the system id 0x" << std::hex << static_cast(ts.descriptor(c,0).sys_id) << std::dec; LOG(info) << "Component " << c << " has the system id 0x" << static_cast(ts.descriptor(c,0).sys_id); /* for (size_t m = 0; m < ts.num_microslices(c); ++m) { PrintMicroSliceDescriptor(ts.descriptor(c,m)); } */ } return true; }