lris::LArRawInputDriverJP250L Class Reference

#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:

• lardata/lardata/RawData/utils/LArRawInputDriverJP250L.h
• lardata/lardata/RawData/utils/LArRawInputDriverJP250L.cxx