#include <LArRawInputDriverJP250L.h>

Public Member Functions
	LArRawInputDriverJP250L (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 Attributes
art::SourceHelper const &	principalMaker_

TTree *	m_eventTree
	TTree containing information from each trigger. More...

unsigned int	m_nEvent
	number of triggers in the TTree More...

unsigned int	m_current
	current entry in the TTree More...

unsigned short	m_runID
	run ID, has to start from 1 More...

unsigned int	m_unixtime
	unix timestamp of the start of the run More...

unsigned short	m_nChannels
	number of channels in the detector More...

unsigned short	m_nSamples
	number of time samples per channel More...

unsigned short *	m_data
	the ADC of each time sample for each channel More...

Detailed Description

Definition at line 27 of file LArRawInputDriverJP250L.h.

Constructor & Destructor Documentation

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

Class to fill the constraints on a template argument to the class, FileReaderSource

Definition at line 35 of file LArRawInputDriverJP250L.cxx.

     : principalMaker_(pm)
     , m_current(0)
     , m_data(0)
   {
     helper.reconstitutes<raw::DAQHeader,              art::InEvent>("daq");
     helper.reconstitutes<std::vector<raw::RawDigit>,  art::InEvent>("daq");
     helper.reconstitutes<sumdata::RunData,            art::InRun>  ("daq");
   }

Member Function Documentation

void lris::LArRawInputDriverJP250L::closeCurrentFile ( )

Definition at line 48 of file LArRawInputDriverJP250L.cxx.

   {
     delete [] m_data;
   }

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

Definition at line 54 of file LArRawInputDriverJP250L.cxx.

   {
     TFile m_f(name.c_str());
     TTree* runTree = dynamic_cast<TTree*>(m_f.Get("runTree"));
 
     m_eventTree    = dynamic_cast<TTree*>(m_f.Get("eventTree"));
     m_eventTree->SetDirectory(0);
     m_nEvent = m_eventTree->GetEntries();
 
     // run information
     runTree->SetBranchAddress("runID",&m_runID);
     runTree->SetBranchAddress("unixtime",&m_unixtime);
     runTree->SetBranchAddress("nChannels",&m_nChannels);
     runTree->SetBranchAddress("nSamples",&m_nSamples);
     runTree->GetEntry(0);
 
     // have to add 1 to the runID because it can't be zero
     // adding 1 to every runID ensures that all runIDs are bumped in the same way
     m_runID += 1;
 
     const int nLength=(m_nSamples+4)*m_nChannels;
     m_data = new unsigned short[nLength];
     m_eventTree->SetBranchAddress("data",m_data);
 
     // Fill and return a new Fileblock.
     // The string tells you what the version of the LArRawInputDriver is
     fb = new art::FileBlock(art::FileFormatVersion(1, "LArRawInputJP250L 2013_01"),
                             name);
 
 
     return;
   }

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

Definition at line 89 of file LArRawInputDriverJP250L.cxx.

   {
 
     if(m_current > m_nEvent) return false;
 
     m_eventTree->GetEntry(m_current);
 
     raw::DAQHeader daqHeader;
     daqHeader.SetRun(m_runID);
     daqHeader.SetTimeStamp(m_unixtime);
     daqHeader.SetNChannels(m_nChannels);
     // set the event number to start at 1 - ART likes to start numbering
     // things from 1
     daqHeader.SetEvent(m_current+1);
 
     // the following DAQHeader fields are not used by JP250L at this time
     //daqHeader.SetStatus(1);
     //daqHeader.SetFixedWord(h1.fixed);
     //daqHeader.SetFileFormat(h1.format);
     //daqHeader.SetSoftwareVersion(h1.software);
     //daqHeader.SetSpareWord(h1.spare);
 
     // make unique_ptrs for the data products to store in the event
     std::unique_ptr<raw::DAQHeader>              daqcol( new raw::DAQHeader(daqHeader)  );
     std::unique_ptr<std::vector<raw::RawDigit> >  rdcol( new std::vector<raw::RawDigit> );
 
     // loop over the signals and break them into
     // one RawDigit for each channel
     std::vector<short> adcVec(m_nSamples,0);
     for(unsigned int n = 0; n < m_nChannels; ++n){
       adcVec.clear();
       for(unsigned int i = 0; i < m_nSamples; ++i){
         adcVec.push_back(m_data[(m_nSamples+4)*n+(4+i)]);
       }
       rdcol->push_back(raw::RawDigit(n,m_nSamples,adcVec));
     }
 
     art::RunNumber_t    rn     = daqHeader.GetRun();
     art::SubRunNumber_t sn     = 1;
     art::EventNumber_t  en     = daqHeader.GetEvent();
     art::Timestamp      tstamp = daqHeader.GetTimeStamp();
 
     // Make the Run and SubRun principals
     // this step is done once per run.
     if (m_current < 1){
       std::unique_ptr<sumdata::RunData> rundata(new sumdata::RunData("jpl250l") );
       outR = principalMaker_.makeRunPrincipal(rn, tstamp);
       outSR = principalMaker_.makeSubRunPrincipal(rn,
                                                   sn,
                                                   tstamp);
       art::put_product_in_principal(std::move(rundata), *outR, "daq");
     }
 
 
     // make the Event principal
     outE = principalMaker_.makeEventPrincipal(rn,
                                               sn,
                                               en,
                                               tstamp);
 
     // Put products in the event.
     // first argument places the desired data product in the file
     // second argument is event record to associate it to
     // third argument is the label of the module storing the
     // information.  "daq" is standard in LArSoft
     art::put_product_in_principal(std::move(rdcol),
                                   *outE,
                                   "daq"); // Module label
     art::put_product_in_principal(std::move(daqcol),
                                   *outE,
                                   "daq"); // Module label
 
     // increment the entry to get out of the TTree for the next time through the loop
     ++m_current;
 
     return true;
   }

Member Data Documentation

unsigned int lris::LArRawInputDriverJP250L::m_current

private

current entry in the TTree

Definition at line 54 of file LArRawInputDriverJP250L.h.

unsigned short* lris::LArRawInputDriverJP250L::m_data

private

the ADC of each time sample for each channel

Definition at line 59 of file LArRawInputDriverJP250L.h.

TTree* lris::LArRawInputDriverJP250L::m_eventTree

private

TTree containing information from each trigger.

Definition at line 52 of file LArRawInputDriverJP250L.h.

unsigned short lris::LArRawInputDriverJP250L::m_nChannels

private

number of channels in the detector

Definition at line 57 of file LArRawInputDriverJP250L.h.

unsigned int lris::LArRawInputDriverJP250L::m_nEvent

private

number of triggers in the TTree

Definition at line 53 of file LArRawInputDriverJP250L.h.

unsigned short lris::LArRawInputDriverJP250L::m_nSamples

private

number of time samples per channel

Definition at line 58 of file LArRawInputDriverJP250L.h.

unsigned short lris::LArRawInputDriverJP250L::m_runID

private

run ID, has to start from 1

Definition at line 55 of file LArRawInputDriverJP250L.h.

unsigned int lris::LArRawInputDriverJP250L::m_unixtime

private

unix timestamp of the start of the run

Definition at line 56 of file LArRawInputDriverJP250L.h.

art::SourceHelper const& lris::LArRawInputDriverJP250L::principalMaker_

private

Definition at line 49 of file LArRawInputDriverJP250L.h.

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation