Inheritance diagram for dune::RceRawDecoder:

Public Member Functions
	RceRawDecoder (fhicl::ParameterSet const &p)

	RceRawDecoder (RceRawDecoder const &)=delete

	RceRawDecoder (RceRawDecoder &&)=delete

RceRawDecoder &	operator= (RceRawDecoder const &)=delete

RceRawDecoder &	operator= (RceRawDecoder &&)=delete

void	produce (art::Event &e) override

void	reconfigure (const fhicl::ParameterSet &pset)

void	beginJob () override

void	setRootObjects ()

Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDProducer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept

	Producer (fhicl::ParameterSet const &)

	Producer (Producer const &)=delete

	Producer (Producer &&)=delete

Producer &	operator= (Producer const &)=delete

Producer &	operator= (Producer &&)=delete

void	doBeginJob (SharedResources const &resources)

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::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

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 Types
typedef std::vector< raw::RawDigit >	RawDigits

Private Member Functions
bool	_process (const artdaq::Fragment &frag, RawDigits &raw_digits)

float	rmsADC (raw::RawDigit::ADCvector_t &wav)

float	meanADC (raw::RawDigit::ADCvector_t &wav)

Private Attributes
std::string	_input_label

std::string	_output_label

bool	_expect_container_fragments

TH1D *	_h_nticks

TH1I *	_h_all_adc_values

TH1I *	_h_crate_numbers

TH1I *	_h_slot_IDs

TH1I *	_h_fiber_IDs

TH1I *	_h_online_channels

TH1I *	_h_offline_channels

std::vector< TH2F * >	_h_mean_slot_channels

std::vector< TH2F * >	_h_rms_slot_channels

std::vector< TProfile * >	_h_mean_slot_channels_pfx

std::vector< TProfile * >	_h_rms_slot_channels_pfx

std::vector< TH2I * >	_h_fiber_persistent_wav

std::vector< int16_t >	_buffer

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)

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 36 of file RceRawDecoder_module.cc.

Member Typedef Documentation

typedef std::vector<raw::RawDigit> dune::RceRawDecoder::RawDigits

private

Definition at line 56 of file RceRawDecoder_module.cc.

Constructor & Destructor Documentation

dune::RceRawDecoder::RceRawDecoder ( fhicl::ParameterSet const & p )

explicit

Definition at line 87 of file RceRawDecoder_module.cc.

   : EDProducer{pset}
 // Initialize member data here.
 {
   art::ServiceHandle<art::TFileService> fs;
   //fs->registerFileSwitchCallback(this, &RceRawDecoder::setRootObjects);
   setRootObjects();
 
   reconfigure(pset);
   // Call appropriate produces<>() functions here.
   //produces< std::vector<raw::RawDigit> > ( _output_label );  
   produces<RawDigits>( _output_label );  
 }

dune::RceRawDecoder::RceRawDecoder ( RceRawDecoder const & )

delete

dune::RceRawDecoder::RceRawDecoder ( RceRawDecoder && )

delete

Member Function Documentation

bool dune::RceRawDecoder::_process	(	const artdaq::Fragment &	frag,
		RawDigits &	raw_digits
	)

private

Definition at line 232 of file RceRawDecoder_module.cc.

 {
   // FIXME: Remove hard-coded fragment type
   if((unsigned)frag.type() != 2) return false;
 
   MF_LOG_INFO("RceRawDecoder")
       << "   SequenceID = " << frag.sequenceID()
       << "   fragmentID = " << frag.fragmentID()
       << "   fragmentType = " << (unsigned)frag.type()
       << "   Timestamp =  " << frag.timestamp();
   art::ServiceHandle<dune::PdspChannelMapService> channelMap;
   dune::RceFragment rce(frag);
   
   uint32_t ch_counter = 0;
   for (int i = 0; i < rce.size(); ++i)
   {
      auto const * rce_stream = rce.get_stream(i);
      int n_ch = rce_stream->getNChannels();
      int n_ticks = rce_stream->getNTicks();
      auto const identifier = rce_stream->getIdentifier();
      uint32_t crateNumber = identifier.getCrate();
      uint32_t slotNumber = identifier.getSlot();
      uint32_t fiberNumber = identifier.getFiber();
      std::cout<<"crate, slot, fiber = "<<crateNumber<<", "<<slotNumber<<", "<<fiberNumber<<std::endl;
      uint32_t fiberID = (crateNumber*5+slotNumber)*4 + fiberNumber;
      _h_fiber_IDs->Fill(fiberID);
      MF_LOG_INFO("RceRawDecoder")
          << "RceFragment timestamp: " << rce_stream->getTimeStamp()
          << ", NChannels: " << n_ch
          << ", NTicks: " << n_ticks;
 
      _h_nticks->Fill(n_ticks);
 
      size_t buffer_size = n_ch * n_ticks;
      if (_buffer.capacity() < buffer_size)
      {
          MF_LOG_INFO("RceRawDecoder")
              << "Increase buffer size from " << _buffer.capacity()
              << " to " << buffer_size;
 
          _buffer.reserve(buffer_size);
      }
 
      int16_t* adcs = _buffer.data();
      rce_stream->getMultiChannelData(adcs);
 
      raw::RawDigit::ADCvector_t v_adc;
      for (int i_ch = 0; i_ch < n_ch; i_ch++)
      {
         if(i==0 && i_ch ==0) std::cout<<" ADCs for the the 100 ticks in the 1st channel of the 1st RCE "<<std::endl;
         v_adc.clear();
         for (int i_tick = 0; i_tick < n_ticks; i_tick++)
         {
             //print ADCs for the the 100 ticks in the 1st channel of the 1st RCE
             if(i==0 && i_ch==0 && i_tick<100) {
               std::cout<<adcs[i_tick]<<"\t";
               if(i_tick==99) std::cout<<std::endl;
             }
             v_adc.push_back(adcs[i_tick]);
             _h_all_adc_values->Fill(adcs[i_tick]);
             //save the 1st waveform in a fiber
             _h_fiber_persistent_wav.at(fiberID)->Fill(i_tick, adcs[i_tick]);
         }
         adcs += n_ticks;
 
         ch_counter++;
         int offlineChannel = -1;
         offlineChannel = channelMap->GetOfflineNumberFromDetectorElements(crateNumber, slotNumber, fiberNumber, i_ch, dune::PdspChannelMapService::kRCE);
         _h_offline_channels->Fill(offlineChannel);
         raw::RawDigit raw_digit(offlineChannel, n_ticks, v_adc);
         raw_digits.push_back(raw_digit);
         
         //=========== fill Mean and RMS with slotchannel number=====================
         uint32_t slotID = crateNumber*5 + slotNumber; //0 - 29
         _h_slot_IDs->Fill(slotID);
         uint32_t slotchannel = 0; // 0 - 127
         slotchannel = fiberNumber*128 + i_ch; //hard-coded. will be moved to analyzer module. Feb 7, 2018
         float mean = meanADC(v_adc);
         float rms = rmsADC(v_adc);
         _h_mean_slot_channels.at(slotID)->Fill(slotchannel, mean);
         _h_rms_slot_channels.at(slotID)->Fill(slotchannel, rms);
         _h_mean_slot_channels_pfx.at(slotID)->Fill(slotchannel, mean, 1);
         _h_rms_slot_channels_pfx.at(slotID)->Fill(slotchannel, rms, 1);
         //==========================================================================
         
         
      }
   }
 
   return true;
 }

void dune::RceRawDecoder::beginJob ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 145 of file RceRawDecoder_module.cc.

145 {

146 }

float dune::RceRawDecoder::meanADC ( raw::RawDigit::ADCvector_t & wav )

private

Definition at line 347 of file RceRawDecoder_module.cc.

   {
     int n = wav.size();
     float sum = 0.;
     for(int i = 0; i < n; i++)
       {
         if (wav[i]!=0) sum += abs(wav[i]);
       }
     return sum / n;
   }

RceRawDecoder& dune::RceRawDecoder::operator= ( RceRawDecoder const & )

delete

RceRawDecoder& dune::RceRawDecoder::operator= ( RceRawDecoder && )

delete

void dune::RceRawDecoder::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 148 of file RceRawDecoder_module.cc.

                                              {
     // TODO Use MF_LOG_DEBUG
     MF_LOG_INFO("RceRawDecoder")
       << "-------------------- RCE RawDecoder -------------------";
 
   RawDigits raw_digits;
   unsigned int n_rce_frags = 0;
 
   if (_expect_container_fragments) {
     art::InputTag itag1(_input_label, "ContainerTPC");
     auto cont_frags = evt.getHandle<artdaq::Fragments>(itag1);
     art::EventNumber_t eventNumber = evt.event();
     // Check if there is Timing data in this event
     // Don't crash code if not present, just don't save anything
     try { cont_frags->size(); }
     catch(std::exception const&) {
       std::cout << "WARNING: Container TPC/RCE data not found in event " << eventNumber << std::endl;
       std::vector<raw::RawDigit> digits;
       evt.put(std::make_unique<std::vector<raw::RawDigit>>(std::move(digits)), _output_label);
       return;
     }
     //Check that the data is valid
     if(!cont_frags){
       MF_LOG_ERROR("RceRawDecoder")
           << "Run: " << evt.run()
                   << ", SubRun: " << evt.subRun()
           << ", Event: " << evt.event()
                   << " Container Fragments is NOT VALID";
     }
     
     for (auto const& cont : *cont_frags)
     {
       artdaq::ContainerFragment cont_frag(cont);
       for (size_t ii = 0; ii < cont_frag.block_count(); ++ii)
         {
           //artdaq::Fragment frag;
           //size_t frag_size = cont_frag.fragSize(ii);
           //frag.resizeBytes(frag_size);
           //memcpy(frag.headerAddress(), cont_frag.at(ii), frag_size);
           if (_process(*cont_frag[ii], raw_digits)) ++n_rce_frags;
         }
     }
   }
   else
   {
     art::InputTag itag2(_input_label, "TPC");
     auto frags = evt.getHandle<artdaq::Fragments>(itag2);
 
     // Check if there is Timing data in this event
     // Don't crash code if not present, just don't save anything
     art::EventNumber_t eventNumber = evt.event();
     try { frags->size(); }
     catch(std::exception const&) {
       std::cout << "WARNING: Raw TPC/RCE data not found in event " << eventNumber << std::endl;
       std::vector<raw::RawDigit> digits;
       evt.put(std::make_unique<std::vector<raw::RawDigit>>(std::move(digits)), _output_label);
       return;
     }
 
     //Check that the data is valid
     if(!frags.isValid()){
       MF_LOG_ERROR("RceRawDecoder")
           << "Run: " << evt.run()
                   << ", SubRun: " << evt.subRun()
           << ", Event: " << evt.event()
                   << " Fragments is NOT VALID";
     }
 
     for(auto const& frag: *frags)
     {
       if (_process(frag, raw_digits)) ++n_rce_frags;
     }
   }
 
   MF_LOG_INFO("RceRawDecoder")
       << " Processed " << n_rce_frags
       << " RCE Fragments, "
       << raw_digits.size()
       << " RawDigits.";
 
   evt.put(std::make_unique<decltype(raw_digits)>(std::move(raw_digits)),
           _output_label);
 }

void dune::RceRawDecoder::reconfigure ( const fhicl::ParameterSet & pset )

Definition at line 101 of file RceRawDecoder_module.cc.

                                                                  {
 
   _input_label = pset.get<std::string>("RawDataLabel");
   _output_label = pset.get<std::string>("OutputDataLabel");
   _expect_container_fragments = pset.get<bool>("ExpectContainerFragments", true);
 }

float dune::RceRawDecoder::rmsADC ( raw::RawDigit::ADCvector_t & wav )

private

Definition at line 329 of file RceRawDecoder_module.cc.

   {
     int n = wav.size();
     float sum = 0.;
     for(int i = 0; i < n; i++){
       if(wav[i]!=0) sum += wav[i];
     }
     float mean = sum / n;
     sum = 0;
     for(int i = 0; i < n; i++)
       {
         if (wav[i]!=0)     sum += (wav[i]-mean)*(wav[i]-mean);
       }
     return sqrt(sum / n);
   }

void dune::RceRawDecoder::setRootObjects ( )

Definition at line 108 of file RceRawDecoder_module.cc.

                                       {
         //clear vectors
         _h_mean_slot_channels.clear();
         _h_rms_slot_channels.clear();
         _h_mean_slot_channels_pfx.clear();
         _h_rms_slot_channels_pfx.clear();
         _h_fiber_persistent_wav.clear();
         // place to define the histograms
   art::ServiceHandle<art::TFileService> file_srv;
   _h_nticks = file_srv->make<TH1D>("rce_NTicks","TPC: Number of ticks",  100, 0, 20000);
   _h_nticks->SetXTitle("NTicks");
   _h_all_adc_values = file_srv->make<TH1I>("rce_All_ADC_values","TPC: All_ADC_values",  5000, 0, 5000);
   _h_all_adc_values->SetXTitle("ADC");
   _h_crate_numbers = file_srv->make<TH1I>("rce_All_crate_numbers","TPC: Crate Numbers",  6, 0, 6);
   _h_crate_numbers->SetXTitle("Crate Number");
   _h_slot_IDs = file_srv->make<TH1I>("rce_All_slot_IDs","TPC: Slot IDs",  30, 0, 30);
   _h_slot_IDs->SetXTitle("Slot ID");
   _h_fiber_IDs = file_srv->make<TH1I>("rce_All_fiber_IDs","TPC: Fiber IDs",  120, 0, 120);
   _h_fiber_IDs->SetXTitle("Fiber ID");
   _h_online_channels = file_srv->make<TH1I>("rce_Online_Channels", "TPC: Online Channel ID", 5000, 0, 5000);
   _h_offline_channels = file_srv->make<TH1I>("rce_Offline_Channels", "TPC: Offline Channel ID", 5000, 0, 5000);
   for(int i=0;i<30;i++) {
         _h_mean_slot_channels.push_back(file_srv->make<TH2F>(Form("Slot%d_Mean", i), Form("Slot%d:Mean_vs_SlotChannel", i), 512, 0, 512, 5000, .0, 5000)); //hard-coded for now. Will be moved to the analyzer module. Feb 7, 2018
         _h_rms_slot_channels.push_back(file_srv->make<TH2F>(Form("Slot%d_RMS", i), Form("Slot%d:RMS_vs_SlotChannel", i), 512, 0, 512, 5000, .0, 5000)); //hard-coded for now. Will be moved to the analyzer module. Feb 7, 2018
         _h_mean_slot_channels_pfx.push_back(file_srv->make<TProfile>(Form("Slot%d_Mean_pfx", i), Form("Slot%d:Mean_vs_SlotChannel_pfx", i), 512, 0, 512)); //hard-coded for now. Will be moved to the analyzer module. Feb 7, 2018
         _h_rms_slot_channels_pfx.push_back(file_srv->make<TProfile>(Form("Slot%d_RMS_pfx", i), Form("Slot%d:RMS_vs_SlotChannel_pfx", i), 512, 0, 512)); //hard-coded for now. Will be moved to the analyzer module. Feb 7, 2018
         
         _h_mean_slot_channels[i]->GetXaxis()->SetTitle("Slot Channel"); _h_mean_slot_channels[i]->GetYaxis()->SetTitle("Raw Mean"); 
         _h_rms_slot_channels[i]->GetXaxis()->SetTitle("Slot Channel"); _h_rms_slot_channels[i]->GetYaxis()->SetTitle("Raw RMS"); 
         _h_mean_slot_channels_pfx[i]->GetXaxis()->SetTitle("Slot Channel"); _h_mean_slot_channels_pfx[i]->GetYaxis()->SetTitle("Profiled Mean"); 
         _h_rms_slot_channels_pfx[i]->GetXaxis()->SetTitle("Slot Channel"); _h_rms_slot_channels_pfx[i]->GetYaxis()->SetTitle("Profiled RMS"); 
   }
   for(int i=0;i<120;i++) {
     _h_fiber_persistent_wav.push_back(file_srv->make<TH2I>(Form("Persistent_Waveform_Fiber#%d", i), Form("Persistent_Waveform_Fiber#%d", i), 10000, 0, 10000, 500, 0, 5000));
   }
   
 }

Member Data Documentation

std::vector<int16_t> dune::RceRawDecoder::_buffer

private

Definition at line 79 of file RceRawDecoder_module.cc.

bool dune::RceRawDecoder::_expect_container_fragments

private

Definition at line 65 of file RceRawDecoder_module.cc.

TH1I* dune::RceRawDecoder::_h_all_adc_values

private

Definition at line 68 of file RceRawDecoder_module.cc.

TH1I* dune::RceRawDecoder::_h_crate_numbers

private

Definition at line 69 of file RceRawDecoder_module.cc.

TH1I* dune::RceRawDecoder::_h_fiber_IDs

private

Definition at line 71 of file RceRawDecoder_module.cc.

std::vector<TH2I*> dune::RceRawDecoder::_h_fiber_persistent_wav

private

Definition at line 78 of file RceRawDecoder_module.cc.

std::vector<TH2F*> dune::RceRawDecoder::_h_mean_slot_channels

private

Definition at line 74 of file RceRawDecoder_module.cc.

std::vector<TProfile*> dune::RceRawDecoder::_h_mean_slot_channels_pfx

private

Definition at line 76 of file RceRawDecoder_module.cc.

TH1D* dune::RceRawDecoder::_h_nticks

private

Definition at line 67 of file RceRawDecoder_module.cc.

TH1I* dune::RceRawDecoder::_h_offline_channels

private

Definition at line 73 of file RceRawDecoder_module.cc.

TH1I* dune::RceRawDecoder::_h_online_channels

private

Definition at line 72 of file RceRawDecoder_module.cc.

std::vector<TH2F*> dune::RceRawDecoder::_h_rms_slot_channels

private

Definition at line 75 of file RceRawDecoder_module.cc.

std::vector<TProfile*> dune::RceRawDecoder::_h_rms_slot_channels_pfx

private

Definition at line 77 of file RceRawDecoder_module.cc.

TH1I* dune::RceRawDecoder::_h_slot_IDs

private

Definition at line 70 of file RceRawDecoder_module.cc.

std::string dune::RceRawDecoder::_input_label

private

Definition at line 63 of file RceRawDecoder_module.cc.

std::string dune::RceRawDecoder::_output_label

private

Definition at line 64 of file RceRawDecoder_module.cc.

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

duneprototypes/duneprototypes/Protodune/singlephase/RawDecoding/RceRawDecoder_module.cc

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation