DAQToOffline::SSPTriggerAna Class Reference

Inheritance diagram for DAQToOffline::SSPTriggerAna:

Public Member Functions
	SSPTriggerAna (fhicl::ParameterSet const &pset)
void	analyze (art::Event const &evt) override
void	reconfigure (fhicl::ParameterSet const &pset)
void	printParameterSet ()
	SSPTriggerAna (SSPTriggerAna const &)=delete
	SSPTriggerAna (SSPTriggerAna &&)=delete
SSPTriggerAna &	operator= (SSPTriggerAna const &)=delete
SSPTriggerAna &	operator= (SSPTriggerAna &&)=delete
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob ()
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	beginJob () override
void	endJob () override
void	beginRun (art::Run const &run) override
void	endRun (art::Run const &run) override

Private Attributes
std::string	fFragType
std::string	fRawDataLabel
std::string	fOutputDataLabel
SSPReformatterAlgs	sspReform
unsigned long int	firstTime
unsigned long int	lastTime
std::map< int, long int >	triggerCount
std::map< int, long int >	triggerCountHeaders
std::map< int, long int >	triggerCountWaveforms
TH1D *	fTriggerRateHist
TH1D *	fTriggerRateHistHeaders
TH1D *	fTriggerRateHistWaveforms

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &paths, fhicl::ParameterSet const &config)
detail::ProcessAndEventSelectors &	processAndEventSelectors ()
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 54 of file SSPTriggerAna_module.cc.

Constructor & Destructor Documentation

DAQToOffline::SSPTriggerAna::SSPTriggerAna ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 98 of file SSPTriggerAna_module.cc.

                                         : art::EDAnalyzer(pset),
                                           sspReform(pset.get<fhicl::ParameterSet>("SSPReformatter"))
{

  this->reconfigure(pset);

  //first_FirstSample = -1;
  //first_TimeStamp = -1;
}

DAQToOffline::SSPTriggerAna::SSPTriggerAna ( SSPTriggerAna const &  )

delete

DAQToOffline::SSPTriggerAna::SSPTriggerAna ( SSPTriggerAna &&  )

delete

Member Function Documentation

void DAQToOffline::SSPTriggerAna::analyze ( art::Event const &  evt )

overridevirtual

< Derived Optical channel

< first sample time in ticks

Implements art::EDAnalyzer.

Definition at line 219 of file SSPTriggerAna_module.cc.

{

  //  std::cout<<"Event: "<<evt.event()<<"Run: "<<evt.run()<<std::endl;

  art::Handle<artdaq::Fragments> rawFragments;
  evt.getByLabel(fRawDataLabel, fFragType, rawFragments);

  // Check if there is SSP data in this event
  // Don't crash code if not present, just don't save anything
  try { rawFragments->size(); }
  catch(std::exception const&) {
    mf::LogWarning("SSPTriggerAna") << "WARNING: Raw SSP data not found in event " << evt.event();
    return;
  }

  // Check that the data is valid
  if(!rawFragments.isValid()){
    mf::LogError("SSPTriggerAna") << "Run: " << evt.run()
                                 << ", SubRun: " << evt.subRun()
                                 << ", Event: " << evt.event()
                                 << " is NOT VALID";
    throw cet::exception("raw NOT VALID");
    return;
  }



  unsigned int numFragments = rawFragments->size();

  for (size_t idx = 0; idx < numFragments; ++idx) {
    const auto& frag((*rawFragments)[idx]);
    lbne::SSPFragment sspf(frag);

    unsigned int nTriggers = sspReform.CheckAndGetNTriggers(frag, sspf);
      
    const unsigned int* dataPointer = sspf.dataBegin();

        
    for (unsigned int triggersProcessed = 0;
         (nTriggers==0 || triggersProcessed < nTriggers) && dataPointer < sspf.dataEnd();
         ++triggersProcessed) {
      
      //
      // The elements of the OpDet Pulse
      //
      unsigned short     OpChannel = -1;       ///< Derived Optical channel
      unsigned long      FirstSample = 0;      ///< first sample time in ticks
        

      // Load the event header, advance the pointer
      auto daqHeader = sspReform.GetHeaderAndAdvance(dataPointer);
      //const SSPDAQ::EventHeader* daqHeader=reinterpret_cast<const SSPDAQ::EventHeader*>(dataPointer);
      //dataPointer += sizeof(SSPDAQ::EventHeader)/sizeof(unsigned int);

      // Get ADC Count, create pointer to adcs
      unsigned int nADC = sspReform.GetWaveformLength(daqHeader);

      //get the information from the header
      try {
        OpChannel = sspReform.GetOpChannel(daqHeader);

        FirstSample = sspReform.GetGlobalFirstSample(daqHeader);
        //TimeStamp = ((double)FirstSample)/fNOvAClockFrequency;

        if (FirstSample < 1e16) {
          sspReform.PrintHeaderInfo(daqHeader);
          mf::LogInfo("SSPTriggerAna") << "Problem timestamp at " << FirstSample << std::endl;
          continue;
        }
      }
      catch (cet::exception const&) {
        continue;
      }

      firstTime = std::min(firstTime, FirstSample);
      lastTime  = std::max(lastTime,  FirstSample);

      triggerCount[OpChannel]++;
      fTriggerRateHist->Fill(OpChannel);

      if (nADC == 0) {
        // Header only
        triggerCountHeaders[OpChannel]++;
        fTriggerRateHistHeaders->Fill(OpChannel);
      }
      else {
        // Waveform to skip
        triggerCountWaveforms[OpChannel]++;
        fTriggerRateHistWaveforms->Fill(OpChannel);
        
        // Advance the dataPointer to the next header
        dataPointer+=nADC/2;
      }

      
    } // End of loop over triggers
  } // End of loop over fragments (rawFragments)

    
}

void DAQToOffline::SSPTriggerAna::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 132 of file SSPTriggerAna_module.cc.

{
  //art::ServiceHandle<art::TFileService> tfs;
  //adc_values_ = tfs->make<TH1D>("adc_values","adc_values",4096,-0.5,4095.5);

  art::ServiceHandle<art::TFileService> tfs;
  fTriggerRateHist = tfs->make<TH1D>("triggerratehist",";Trigger Rate vs OpChannel;Rate (Hz)",116,-10,106);
  fTriggerRateHistHeaders = tfs->make<TH1D>("triggerratehist_headers",";Trigger Rate vs OpChannel;Rate (Hz)",116,-10,106);
  fTriggerRateHistWaveforms = tfs->make<TH1D>("triggerratehist_waveforms",";Trigger Rate vs OpChannel;Rate (Hz)",116,-10,106);

  firstTime = (((unsigned long int)1)<<63);
  lastTime = 0;
}

void DAQToOffline::SSPTriggerAna::beginRun ( art::Run const &  run )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 148 of file SSPTriggerAna_module.cc.

{
  fTriggerRateHist->Reset();
  fTriggerRateHistHeaders->Reset();
  fTriggerRateHistWaveforms->Reset();
}

void DAQToOffline::SSPTriggerAna::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 176 of file SSPTriggerAna_module.cc.

{
  //delete adc_values_;
  delete fTriggerRateHist;
  delete fTriggerRateHistHeaders;
  delete fTriggerRateHistWaveforms;
  
  long int deltaT = lastTime-firstTime;
  double deltaTus =  ((double)deltaT)/sspReform.ClockFrequency();


  mf::LogVerbatim("SSPTriggerAna") << "!! Trigger Rate Report." << std::endl;
  mf::LogVerbatim("SSPTriggerAna") << "!! Time: " << deltaTus / 60.e6 << " minutes." << std::endl;

  for (auto itr = triggerCount.begin(); itr != triggerCount.end(); itr++) {
    double freq = ((double)itr->second) / deltaTus * 1e6;
    mf::LogVerbatim("SSPTriggerAna") << "!!    Channel " << itr->first << ": " << std::setiosflags(std::ios::fixed) << std::setprecision(2) << std::setw(7) << freq << " Hz" << std::endl;
  }
  
  for (auto itr = triggerCountHeaders.begin(); itr != triggerCountHeaders.end(); itr++) {
    double freq = ((double)itr->second) / deltaTus * 1e6;
    mf::LogVerbatim("SSPTriggerAna") << "!!  H Channel " << itr->first << ": " << std::setiosflags(std::ios::fixed) << std::setprecision(2) << std::setw(7) << freq << " Hz" << std::endl;
  }
  
  for (auto itr = triggerCountWaveforms.begin(); itr != triggerCountWaveforms.end(); itr++) {
    double freq = ((double)itr->second) / deltaTus * 1e6;
    mf::LogVerbatim("SSPTriggerAna") << "!!  W Channel " << itr->first << ": " << std::setiosflags(std::ios::fixed) << std::setprecision(2) << std::setw(7) << freq << " Hz" << std::endl;
  }

  
  /*
  mf::LogInfo("SSPTriggerAna") << "firstSample:  " << firstTime << " samples\n"
                               << "lastSample:   " << lastTime  << " samples\n"
                               << "totalSamples: " << deltaT    << " samples\n"
                               << "totalTime:    " << deltaTus  << " us\n"
                               << "# Channels:   " << channels.size() << "\n"
                               << "# Triggers:   " << triggerCount << "\n"
                               << "Frequency:    " << freq << "kHz\n";
  */
}

void DAQToOffline::SSPTriggerAna::endRun ( art::Run const &  run )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 155 of file SSPTriggerAna_module.cc.

{
  long int deltaT = lastTime-firstTime;
  double deltaTus =  ((double)deltaT)/sspReform.ClockFrequency();

  fTriggerRateHist->Scale(1.e6/deltaTus);
  fTriggerRateHistHeaders->Scale(1.e6/deltaTus);
  fTriggerRateHistWaveforms->Scale(1.e6/deltaTus);
  
  fTriggerRateHist->SetTitle(Form("Run %i Trigger Rate",run.run()));
  fTriggerRateHistHeaders->SetTitle(Form("Run %i Header Trigger Rate",run.run()));
  fTriggerRateHistWaveforms->SetTitle(Form("Run %i Wavefrom Trigger Rate",run.run()));

  fTriggerRateHist->Write(Form("run_%i_trigger_rate",run.run()));
  fTriggerRateHistHeaders->Write(Form("run_%i_trigger_rate_headers",run.run()));
  fTriggerRateHistWaveforms->Write(Form("run_%i_trigger_rate_waveforms",run.run()));
  
}

SSPTriggerAna& DAQToOffline::SSPTriggerAna::operator= ( SSPTriggerAna const &  )

delete

SSPTriggerAna& DAQToOffline::SSPTriggerAna::operator= ( SSPTriggerAna &&  )

delete

void DAQToOffline::SSPTriggerAna::printParameterSet

(

)

Definition at line 122 of file SSPTriggerAna_module.cc.

                                                 {

  mf::LogDebug("SSPTriggerAna") << "===================================="   << "\n"
                               << "Parameter Set"                          << "\n"
                               << "===================================="   << "\n"
                               << "fFragType:        " << fFragType        << "\n"
                               << "fRawDataLabel:    " << fRawDataLabel    << "\n"
                               << "===================================="   << "\n";
}

void DAQToOffline::SSPTriggerAna::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 109 of file SSPTriggerAna_module.cc.

                                                                        {

  fFragType           = pset.get<std::string>("FragType");
  fRawDataLabel       = pset.get<std::string>("RawDataLabel");

  //fDebug = pset.get<bool>("Debug");
  //fZeroThreshold=0;
  //fCompression=raw::kNone;

  printParameterSet();
  
}

Member Data Documentation

std::string DAQToOffline::SSPTriggerAna::fFragType

private

Definition at line 79 of file SSPTriggerAna_module.cc.

unsigned long int DAQToOffline::SSPTriggerAna::firstTime

private

Definition at line 85 of file SSPTriggerAna_module.cc.

std::string DAQToOffline::SSPTriggerAna::fOutputDataLabel

private

Definition at line 81 of file SSPTriggerAna_module.cc.

std::string DAQToOffline::SSPTriggerAna::fRawDataLabel

private

Definition at line 80 of file SSPTriggerAna_module.cc.

TH1D* DAQToOffline::SSPTriggerAna::fTriggerRateHist

private

Definition at line 91 of file SSPTriggerAna_module.cc.

TH1D* DAQToOffline::SSPTriggerAna::fTriggerRateHistHeaders

private

Definition at line 92 of file SSPTriggerAna_module.cc.

TH1D* DAQToOffline::SSPTriggerAna::fTriggerRateHistWaveforms

private

Definition at line 93 of file SSPTriggerAna_module.cc.

unsigned long int DAQToOffline::SSPTriggerAna::lastTime

private

Definition at line 86 of file SSPTriggerAna_module.cc.

SSPReformatterAlgs DAQToOffline::SSPTriggerAna::sspReform

private

Definition at line 83 of file SSPTriggerAna_module.cc.

std::map<int, long int> DAQToOffline::SSPTriggerAna::triggerCount

private

Definition at line 87 of file SSPTriggerAna_module.cc.

std::map<int, long int> DAQToOffline::SSPTriggerAna::triggerCountHeaders

private

Definition at line 88 of file SSPTriggerAna_module.cc.

std::map<int, long int> DAQToOffline::SSPTriggerAna::triggerCountWaveforms

private

Definition at line 89 of file SSPTriggerAna_module.cc.

---

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

• dunetpc/dune/daqinput35t/SSPTriggerAna_module.cc