#include <RawData311InputDriver.h>

Public Member Functions
	RawData311InputDriver (fhicl::ParameterSet const &pset, art::ProductRegistryHelper &helper, art::SourceHelper const &pm)

void	closeCurrentFile ()

void	readFile (std::string const &name, art::FileBlock *&fb)

bool	readNext (art::RunPrincipal const &inR, art::SubRunPrincipal const &inSR, art::RunPrincipal &outR, art::SubRunPrincipal &outSR, art::EventPrincipal *&outE)

Private Member Functions
void	process_Event311 (std::vector< raw::RawDigit > &digitList, dlardaq::evheader_t &event_head, uint16_t evt_num)

double	GetPedMean (size_t LAr_chan, std::vector< std::pair< double, double > > *fPedMap)

double	GetPedRMS (size_t LAr_chan, std::vector< std::pair< double, double > > *fPedMap)

Private Attributes
art::SourceHelper const &	fSourceHelper

art::SubRunID	fCurrentSubRunID

uint16_t	fEventCounter

uint16_t	fNEvents

std::string	filename

short	nchannels = 1280

short	nsamples = 1667

dlardaq::EventDecoder	DataDecode

dlardaq::runheader_t	file_head

dlardaq::footer_t	file_foot

std::vector< std::pair< double, double > >	fPedMap

std::string	fPedestalFile

Detailed Description

Definition at line 32 of file RawData311InputDriver.h.

Constructor & Destructor Documentation

lris::RawData311InputDriver::RawData311InputDriver	(	fhicl::ParameterSet const &	pset,
		art::ProductRegistryHelper &	helper,
		art::SourceHelper const &	pm
	)

Definition at line 154 of file RawData311InputDriver_service.cc.

     :
     fSourceHelper(pm),
     fCurrentSubRunID(),
     fEventCounter(0),
     DataDecode(nchannels, nsamples)
   {
     fPedestalFile = p.get<std::string>("PedestalFile");
     helper.reconstitutes<std::vector<raw::RawDigit>, art::InEvent>("daq");
   }

Member Function Documentation

void lris::RawData311InputDriver::closeCurrentFile ( )

Definition at line 169 of file RawData311InputDriver_service.cc.

   {
     mf::LogInfo(__FUNCTION__)<<"File boundary: processed " <<fEventCounter <<" events out of " <<fNEvents <<"\n";
     DataDecode.m_file.close();
   }

double lris::RawData311InputDriver::GetPedMean	(	size_t	LAr_chan,
		std::vector< std::pair< double, double > > *	fPedMap
	)

inlineprivate

Definition at line 73 of file RawData311InputDriver.h.

73 { return fPedMap->at(LAr_chan).first; }

double lris::RawData311InputDriver::GetPedRMS	(	size_t	LAr_chan,
		std::vector< std::pair< double, double > > *	fPedMap
	)

inlineprivate

Definition at line 75 of file RawData311InputDriver.h.

75 { return fPedMap->at(LAr_chan).second; }

void lris::RawData311InputDriver::process_Event311	(	std::vector< raw::RawDigit > &	digitList,
		dlardaq::evheader_t &	event_head,
		uint16_t	evt_num
	)

private

Definition at line 118 of file RawData311InputDriver_service.cc.

   {
     // one digit for every wire on each plane
     digitList.clear();
     digitList.resize(nchannels);
     // Get the data.
     std::vector<dlardaq::adc16_t> ADCvec311;
     DataDecode.GetEvent(evt_num, event_head, ADCvec311);
     std::vector<dlardaq::adc16_t> *ADCvec311Pointer = &ADCvec311;
     // fill the wires
     std::vector<short> adclist;
 
     for(size_t LAr_chan = 0; LAr_chan < (size_t)nchannels; LAr_chan++)
     {
       adclist.clear();
       size_t Chan311 = Get311Chan(LAr_chan);
       SplitAdc(ADCvec311Pointer, Chan311, nsamples, adclist);
       short unsigned int nTickReadout = nsamples;
       raw::ChannelID_t channel = LAr_chan;
       raw::Compress_t comp = raw::kNone;
       raw::RawDigit rd(channel, nTickReadout, adclist, comp);
 
       double pedval = RawData311InputDriver::GetPedMean(Chan311, &fPedMap);
       //std::cout << "Pedval: " << pedval << "\n";
       double pedrms = RawData311InputDriver::GetPedRMS(Chan311, &fPedMap);
       rd.SetPedestal(pedval, pedrms);
 
       digitList[LAr_chan] = rd;
     }
   }// process_Event311

void lris::RawData311InputDriver::readFile	(	std::string const &	name,
		art::FileBlock *&	fb
	)

Definition at line 177 of file RawData311InputDriver_service.cc.

   {
     // Read in the pedestal file
     ReadPedestalFile(fPedestalFile, fPedMap);
 
     filename = name;
     // Fill and return a new Fileblock
     fb = new art::FileBlock(art::FileFormatVersion(1, "311 RawInput 2017"), name);
 
     // Prepare the EventDecoder
     //short nchannels_Vplane = 320;
     //short nchannels_Zplane = 960;
     //short nchannels = nchannels_Vplane + nchannels_Zplane;
     //uint32_t nsamples = 1667;
     //DataDecode(nchannels, nsamples);
 
     DataDecode.m_file.open(name.c_str(), std::ios_base::in | std::ios_base::binary);
     if( !DataDecode.m_file.is_open() )
     {
       throw art::Exception( art::errors::FileReadError )
         << "failed to open input file " << name << "\n";
     }
 
     // Read in the file header.
     std::vector<dlardaq::BYTE> head_buf;
     head_buf.resize(dlardaq::RunHeadSz);
 
     DataDecode.m_file.read(&head_buf[0], head_buf.size());
     dlardaq::decode_runhead(&head_buf[0], file_head);
 
     // Define start of event data.
     std::streampos data_start = DataDecode.m_file.tellg();
 
     // Read in the file footer.
     std::vector<dlardaq::BYTE> foot_buf;
     foot_buf.resize(dlardaq::FileFootSz);
 
     DataDecode.m_file.seekg(-dlardaq::FileFootSz, DataDecode.m_file.end);
     DataDecode.m_file.read(&foot_buf[0], foot_buf.size());
     dlardaq::decode_filefoot(&foot_buf[0], file_foot);
     DataDecode.m_file.seekg(data_start);
 
     fNEvents = file_foot.num_events;
     if(fNEvents > 0 && fNEvents < 1000) //There should be 335 events at most in one file.
     {
       mf::LogInfo("")<<"Opened file " <<name <<" with " << fNEvents <<" events." <<"\n";
     } else
     {
       throw art::Exception( art::errors::FileReadError )
         <<"File " <<name <<" seems to have too many events: " <<fNEvents <<"\n";
     }
   }

bool lris::RawData311InputDriver::readNext	(	art::RunPrincipal *const &	inR,
		art::SubRunPrincipal *const &	inSR,
		art::RunPrincipal *&	outR,
		art::SubRunPrincipal *&	outSR,
		art::EventPrincipal *&	outE
	)

Definition at line 233 of file RawData311InputDriver_service.cc.

   {
     if(fEventCounter == fNEvents)
     {
       mf::LogInfo(__FUNCTION__)<<"All the files have been read in. Checking end of file..." <<"\n";
       std::streampos current_position = DataDecode.m_file.tellg();
       DataDecode.m_file.seekg(0, std::ios::end);
       std::streampos file_length = DataDecode.m_file.tellg();
       if( ((uint8_t)file_length - (uint8_t)current_position) > (uint8_t)100 )
       {
         throw art::Exception( art::errors::FileReadError )
           <<"Processed " <<fEventCounter <<" events out of " <<fNEvents <<" but there are still "
           <<(file_length - current_position) <<" bits left." <<"\n";
       }
       mf::LogInfo(__FUNCTION__)<<"Completed reading file and closing output file." <<"\n";
       return false; //Tells readNext that all events have been read.
     }
 
     mf::LogInfo(__FUNCTION__)<<"Reading event " << fEventCounter << " from " << fNEvents <<"\n";
 
     // Create empty result, then fill it from current file
     dlardaq::evheader_t event_head;
     std::unique_ptr< std::vector<raw::RawDigit> > tpc_raw_digits( new std::vector<raw::RawDigit> );
     process_Event311(*tpc_raw_digits, event_head, fEventCounter++);
 
 
 
     // Prepare some stuff
     art::RunNumber_t rn;
     rn = file_head.run_num;
 
     std::uint64_t tthi = event_head.trig_info.ts.tv_sec;
     std::uint64_t thilo = (tthi << 32) + event_head.trig_info.ts.tv_nsec;
     art::Timestamp tstamp(thilo);
 
     size_t pos = filename.find("-");
                 size_t end = filename.find(".");
     std::string str_subrn = filename.substr(pos+1, end-1);
     int int_subrn = std::stoi(str_subrn);
     art::SubRunNumber_t subrn = int_subrn;
     art::SubRunID newID(rn, subrn);
     if( fCurrentSubRunID.runID() != newID.runID() )
     {
       outR = fSourceHelper.makeRunPrincipal(rn, tstamp);
     }
     if( fCurrentSubRunID != newID )
     {
       outSR = fSourceHelper.makeSubRunPrincipal(rn, subrn, tstamp);
       fCurrentSubRunID = newID;
     }
 
     outE = fSourceHelper.makeEventPrincipal(fCurrentSubRunID.run(), fCurrentSubRunID.subRun(),
                                             fEventCounter - 1, tstamp);
 
     // Put products in the event.
     art::put_product_in_principal(std::move(tpc_raw_digits), *outE, "daq");
 
     return true;
    }

Member Data Documentation

dlardaq::EventDecoder lris::RawData311InputDriver::DataDecode

private

Definition at line 61 of file RawData311InputDriver.h.

art::SubRunID lris::RawData311InputDriver::fCurrentSubRunID

private

Definition at line 53 of file RawData311InputDriver.h.

uint16_t lris::RawData311InputDriver::fEventCounter

private

Definition at line 55 of file RawData311InputDriver.h.

dlardaq::footer_t lris::RawData311InputDriver::file_foot

private

Definition at line 64 of file RawData311InputDriver.h.

dlardaq::runheader_t lris::RawData311InputDriver::file_head

private

Definition at line 63 of file RawData311InputDriver.h.

std::string lris::RawData311InputDriver::filename

private

Definition at line 57 of file RawData311InputDriver.h.

uint16_t lris::RawData311InputDriver::fNEvents

private

Definition at line 56 of file RawData311InputDriver.h.

std::string lris::RawData311InputDriver::fPedestalFile

private

Definition at line 67 of file RawData311InputDriver.h.

std::vector< std::pair<double, double> > lris::RawData311InputDriver::fPedMap

private

Definition at line 66 of file RawData311InputDriver.h.

art::SourceHelper const& lris::RawData311InputDriver::fSourceHelper

private

Definition at line 52 of file RawData311InputDriver.h.

short lris::RawData311InputDriver::nchannels = 1280

private

Definition at line 59 of file RawData311InputDriver.h.

short lris::RawData311InputDriver::nsamples = 1667

private

Definition at line 60 of file RawData311InputDriver.h.

The documentation for this class was generated from the following files:

duneprototypes/duneprototypes/3x1x1dp/DataImport/Services/RawData311InputDriver.h
duneprototypes/duneprototypes/3x1x1dp/DataImport/Services/RawData311InputDriver_service.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation