nlana::SSPMonitor Class Reference

Inheritance diagram for nlana::SSPMonitor:

Public Member Functions
	SSPMonitor (const fhicl::ParameterSet &)
virtual	~SSPMonitor ()
void	beginJob ()
void	analyze (const art::Event &)
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 (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::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 &)
fhicl::ParameterSetID	selectorConfig () 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	calculateFFT (TH1D *hist_waveform, TH1D *graph_frequency)

Private Attributes
std::string	fOpDetWaveformModuleLabel
std::string	fOpHitModuleLabel
double	ADC_max
double	ADC_min
int	min_time
int	max_time
int	timesamples
double	startTime
bool	haveStartTime
TH1D *	adc_values_
TH1D *	peaks_
TH1D *	areas_
TH1D *	hit_times_
std::map< size_t, TH1D * >	chan_peaks_
std::map< size_t, TH1D * >	chan_areas_
std::map< size_t, TH2D * >	persistent_waveform_
std::map< size_t, TH1D * >	fft_
std::set< size_t >	has_waveform
TH1I *	fHEventNumber

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
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 45 of file SSPMonitor_module.cc.

Constructor & Destructor Documentation

nlana::SSPMonitor::SSPMonitor

(

const fhicl::ParameterSet &

)

Definition at line 95 of file SSPMonitor_module.cc.

  : EDAnalyzer(pset)
{

  // Get channel map
  art::ServiceHandle<dune::PdspChannelMapService> channelMap;

  fOpDetWaveformModuleLabel = pset.get<std::string>("OpDetWaveformLabel");
  fOpHitModuleLabel = pset.get<std::string>("OpHitLabel");
  ADC_min=pset.get<int>("SSP_ADC_min");
  ADC_max=pset.get<int>("SSP_ADC_max");
  timesamples=pset.get<int>("SSP_TIMESAMPLES");
  min_time=pset.get<int>("SSP_min_time");
  max_time=pset.get<int>("SSP_max_time");

  haveStartTime = false;

  // summary histogram creation

  art::ServiceHandle<art::TFileService> tFileService;
  fHEventNumber = tFileService->make<TH1I>("EventNumber","SSP: EventNumber;Event Number",  100, 0, 10000);
  adc_values_ = tFileService->make<TH1D>("ssp_adc_values","SSP: ADC_Values;ADC Value",4096,-0.5,4095.5);  
  peaks_ = tFileService->make<TH1D>("peaks","Peak Amplitudes;Peak Amplitude",100,-30,50);
  areas_ = tFileService->make<TH1D>("areas","Hit Areas;Hit Area",100,-10000,15000);

  hit_times_ = tFileService->make<TH1D>("ssp_hit_times","Hit Times",1000,min_time,max_time);  
  hit_times_->GetYaxis()->SetTitle("Number of hits");
  hit_times_->GetXaxis()->SetTitle("Time [s]");
}

nlana::SSPMonitor::~SSPMonitor ( )

virtual

Definition at line 127 of file SSPMonitor_module.cc.

{}

Member Function Documentation

void nlana::SSPMonitor::analyze

(

const art::Event &

evt

)

Save the waveforms for all traces like an oscilloscope

Definition at line 135 of file SSPMonitor_module.cc.

{

  art::ServiceHandle<art::TFileService> tFileService;

  auto OpHitHandle = evt.getHandle< std::vector< recob::OpHit > >(fOpHitModuleLabel);

  auto OpDetWaveformHandle = evt.getHandle< std::vector< raw::OpDetWaveform > >(fOpDetWaveformModuleLabel);

  fHEventNumber->Fill(evt.event());

  for(unsigned int i = 0; i < OpHitHandle->size(); ++i)
    {
      art::Ptr< recob::OpHit > ohp(OpHitHandle, i);
      //recob::OpHit TheOpHit = *TheOpHitPtr;  -- see if we can use the ptr straight up.
      auto channel = ohp->OpChannel();
      if (chan_peaks_.find(channel) == chan_peaks_.end())
        {
          chan_peaks_[channel] = tFileService->make<TH1D>(Form("peaks_%d",channel),Form("Peak Amplitudes_%d;Peak Amplitude",channel),100,-30,50);
        }
      if (chan_areas_.find(channel) == chan_areas_.end())
        {
          chan_areas_[channel] = tFileService->make<TH1D>(Form("areas_%d",channel),Form("Peak Areas_%d;Peak Area",channel),100,-10000,15000);
        }
      peaks_->Fill(ohp->Amplitude());
      chan_peaks_[channel]->Fill(ohp->Amplitude());
      areas_->Fill(ohp->Area());
      chan_areas_[channel]->Fill(ohp->Area());

      double time = ohp->PeakTimeAbs()*1E-6;
      if (!haveStartTime) 
        {
          haveStartTime = true;
          startTime = time;
        }
      hit_times_->Fill(time - startTime);

    } // End loop over OpHits

  for(size_t iwaveform = 0; iwaveform < OpDetWaveformHandle->size(); ++iwaveform)
    {
      art::Ptr< raw::OpDetWaveform > odp(OpDetWaveformHandle, iwaveform);

      size_t nADC = odp->size();

      TH1D *hist = new TH1D("hist","hist",nADC,0,nADC);
      auto channel = odp->ChannelNumber();
      if (persistent_waveform_.find(channel) == persistent_waveform_.end())
        {
          TH2D* pwave = tFileService->make<TH2D>(Form("persistent_waveform_%d",channel),Form("persistent_waveform_%d",channel), 500,0,timesamples, (int)(ADC_max-ADC_min),ADC_min,ADC_max);
          pwave->SetTitle(Form("Persistent waveform - Channel %d",channel));
          pwave->GetYaxis()->SetTitle("ADC value");
          pwave->GetXaxis()->SetTitle("Time sample");
          persistent_waveform_[channel] = pwave;
        }
      if (fft_.find(channel) == fft_.end())
        {
          TH1D *fftp = tFileService->make<TH1D>(Form("fft_channel_%d",channel),Form("fft_channel_%d",channel), 100,0,4);
          fftp->SetTitle(Form("FFT - Channel %d",channel));
          fftp->GetXaxis()->SetTitle("Frequency [MHz]");          
          fft_[channel] = fftp;
        }

      TH2D *pwavep = persistent_waveform_[channel];
      for (size_t iadc=0; iadc < nADC; ++iadc)
        {
          auto adcval = odp->at(iadc);
          adc_values_->Fill(adcval);

          ///> Save the waveforms for all traces like an oscilloscope
          pwavep->Fill(iadc,adcval,1);  // (x,y,weight=1)
          hist->SetBinContent(iadc+1,adcval);
        }

      // save one waveform per channel

      if ( has_waveform.find(channel) == has_waveform.end() )
        {
          has_waveform.insert(channel);
          char histname[100];
          sprintf(histname,"evt%i_channel%d",evt.event(), channel);
          TH1D *htf = tFileService->make<TH1D>(histname,histname,nADC,0,nADC);
          htf->Reset();
          htf->Add(hist);
        }
          
      // FFT on the single waveform, output divided by channel
      calculateFFT(hist, fft_[channel]);
      
      hist->Delete();

    } // End loop over OpDetWaveforms

}

void nlana::SSPMonitor::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 131 of file SSPMonitor_module.cc.

{}

void nlana::SSPMonitor::calculateFFT ( TH1D *  hist_waveform, TH1D *  graph_frequency  )

private

Definition at line 230 of file SSPMonitor_module.cc.

                                                                            {
  
  int n_bins = hist_waveform->GetNbinsX();
  TH1* hist_transform = 0;

  // Create hist_transform from the input hist_waveform
  hist_transform = hist_waveform->FFT(hist_transform, "MAG");
  hist_transform -> Scale (1.0 / float(n_bins));
  int nFFT=hist_transform->GetNbinsX();
  
  Double_t frequency;
  Double_t amplitude;
  
  // Loop on the hist_transform to fill the hist_transform_frequency                                                                                        
  for (int k = 2; k <= nFFT/40; ++k){

    frequency =  (k-1)/(n_bins/150.); // MHz
    amplitude = hist_transform->GetBinContent(k);

    hist_frequency->Fill(frequency, amplitude);
  }

  hist_transform->Delete();
  
}

Member Data Documentation

double nlana::SSPMonitor::ADC_max

private

Definition at line 61 of file SSPMonitor_module.cc.

double nlana::SSPMonitor::ADC_min

private

Definition at line 62 of file SSPMonitor_module.cc.

TH1D* nlana::SSPMonitor::adc_values_

private

Definition at line 73 of file SSPMonitor_module.cc.

TH1D* nlana::SSPMonitor::areas_

private

Definition at line 75 of file SSPMonitor_module.cc.

std::map<size_t,TH1D*> nlana::SSPMonitor::chan_areas_

private

Definition at line 82 of file SSPMonitor_module.cc.

std::map<size_t,TH1D*> nlana::SSPMonitor::chan_peaks_

private

Definition at line 81 of file SSPMonitor_module.cc.

std::map<size_t,TH1D*> nlana::SSPMonitor::fft_

private

Definition at line 84 of file SSPMonitor_module.cc.

TH1I* nlana::SSPMonitor::fHEventNumber

private

Definition at line 87 of file SSPMonitor_module.cc.

std::string nlana::SSPMonitor::fOpDetWaveformModuleLabel

private

Definition at line 59 of file SSPMonitor_module.cc.

std::string nlana::SSPMonitor::fOpHitModuleLabel

private

Definition at line 60 of file SSPMonitor_module.cc.

std::set<size_t> nlana::SSPMonitor::has_waveform

private

Definition at line 85 of file SSPMonitor_module.cc.

bool nlana::SSPMonitor::haveStartTime

private

Definition at line 69 of file SSPMonitor_module.cc.

TH1D* nlana::SSPMonitor::hit_times_

private

Definition at line 76 of file SSPMonitor_module.cc.

int nlana::SSPMonitor::max_time

private

Definition at line 64 of file SSPMonitor_module.cc.

int nlana::SSPMonitor::min_time

private

Definition at line 63 of file SSPMonitor_module.cc.

TH1D* nlana::SSPMonitor::peaks_

private

Definition at line 74 of file SSPMonitor_module.cc.

std::map<size_t,TH2D*> nlana::SSPMonitor::persistent_waveform_

private

Definition at line 83 of file SSPMonitor_module.cc.

double nlana::SSPMonitor::startTime

private

Definition at line 68 of file SSPMonitor_module.cc.

int nlana::SSPMonitor::timesamples

private

Definition at line 66 of file SSPMonitor_module.cc.

---

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

• duneprototypes/duneprototypes/Protodune/singlephase/NearlineMonitor/SSPMonitor_module.cc