Inheritance diagram for DAQToOffline::SSPDiagnosticAna:

Public Member Functions
	SSPDiagnosticAna (fhicl::ParameterSet const &pset)

void	analyze (art::Event const &evt) override

void	reconfigure (fhicl::ParameterSet const &pset)

void	printParameterSet ()

	SSPDiagnosticAna (SSPDiagnosticAna const &)=delete

	SSPDiagnosticAna (SSPDiagnosticAna &&)=delete

SSPDiagnosticAna &	operator= (SSPDiagnosticAna const &)=delete

SSPDiagnosticAna &	operator= (SSPDiagnosticAna &&)=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

void	beginEvent (art::EventNumber_t eventNumber)

void	endEvent (art::EventNumber_t eventNumber)

Private Attributes
std::string	fFragType

std::string	fRawDataLabel

std::string	fOutputDataLabel

std::string	fInputModule

std::string	fInputLabel

double	fSampleFreq

const int	m1_window = 10

const int	i1_window = 500

const int	i2_window = 500

std::map< int, TH2D * >	averageWaveforms

std::map< int, TH1D * >	waveformFFTs

TH1D *	fPulseAmplitude

std::map< int, TH1D * >	PulseAmplitudePerChannel

std::map< int, TH1D * >	IntegratedChargePerChannel

std::map< int, TH2D * >	PulseAmplitudeVsIntegratedChargePerChannel

std::map< int, TH2D * >	PulseAmpVsRun

TSpectrum *	specAnalyzer = new TSpectrum( 100 )

SSPReformatterAlgs	sspReform

unsigned long int	firstTime

unsigned long int	lastTime

std::map< int, long int >	triggerCount

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 59 of file SSPDiagnosticAna_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 118 of file SSPDiagnosticAna_module.cc.

                                          : art::EDAnalyzer(pset),
                                            sspReform(pset.get<fhicl::ParameterSet>("SSPReformatter"))
 {
 
   this->reconfigure(pset);
 
   //first_FirstSample = -1;
   //first_TimeStamp = -1;
 }

DAQToOffline::SSPDiagnosticAna::SSPDiagnosticAna ( SSPDiagnosticAna const & )

delete

DAQToOffline::SSPDiagnosticAna::SSPDiagnosticAna ( SSPDiagnosticAna && )

delete

Member Function Documentation

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

overridevirtual

< Derived Optical channel

< first sample time in ticks

Implements art::EDAnalyzer.

Definition at line 364 of file SSPDiagnosticAna_module.cc.

 {
 
   art::ServiceHandle<art::TFileService> tfs;
   art::TFileDirectory RunDir = tfs->mkdir( TString::Format("r%03i", evt.run()).Data(), "SSP Diagnostics by Run" );
   //art::ServiceHandle<geo::Geometry> geo;
 
   art::Handle<artdaq::Fragments> rawFragments;
   evt.getByLabel(fRawDataLabel, fFragType, rawFragments);
 
   mf::LogInfo("SSPDiagnosticAna") << "Starting event analysis";
 
   // 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("SSPDiagnosticAna") << "WARNING: Raw SSP data not found in event " << evt.event();
     return;
   }
 
   // Check that the data is valid
   if(!rawFragments.isValid()){
     mf::LogError("SSPDiagnosticAna") << "Run: " << evt.run()
                                  << ", SubRun: " << evt.subRun()
                                  << ", Event: " << evt.event()
                                  << " is NOT VALID";
     throw cet::exception("raw NOT VALID");
     return;
   }
 
   mf::LogInfo("SSPDiagnosticAna") << "Data is valid!";
 
   // Get OpDetWaveforms from the event
   art::Handle< std::vector< raw::OpDetWaveform > > waveformHandle;
   evt.getByLabel(fInputModule, fInputLabel, waveformHandle);
 
   for (size_t i = 0; i < waveformHandle->size(); i++)
     {
       
       // This was probably required to overcome the "const" problem 
       // with OpDetPulse::Waveform()
       art::Ptr< raw::OpDetWaveform > waveformPtr(waveformHandle, i);
       raw::OpDetWaveform pulse = *waveformPtr;
       int channel = pulse.ChannelNumber();
 
       // Create the TH1 if it doesn't exist
       auto waveform = averageWaveforms.find( channel );
       if ( waveform == averageWaveforms.end() ) {
         TString histName = TString::Format("avgwaveform_channel_%03i", channel);
         //averageWaveforms[channel] =  RunDir.make< TH1D >(histName, ";t (us);", pulse.size(), 0, double(pulse.size()) / fSampleFreq);
         TString histTitle = TString::Format("Average Waveform for OP Channel %03i;t (us);amplitude (ADC)", channel);
         averageWaveforms[channel] =  RunDir.make< TH2D >(histName, histTitle, pulse.size(), 0, double(pulse.size()) / fSampleFreq, 2000, 1200, 5200);
       }
 
       // Add this waveform to this histogram
       for (size_t tick = 0; tick < pulse.size(); tick++) {
         averageWaveforms[channel]->Fill(double(tick)/fSampleFreq, pulse[tick]);
 
       }
   //     // Count number of waveforms on each channel
   //     waveformCount[channel]++;
       
             
     }
 
 
   unsigned int numFragments = rawFragments->size();
   mf::LogInfo("SSPDiagnosticAna") << "Number of fragments = " << numFragments;
 
   for (size_t idx = 0; idx < numFragments; ++idx) {
 
      mf::LogInfo("SSPDiagnosticAna") << "Number of fragments = " << idx;
 
     const auto& frag((*rawFragments)[idx]);
     lbne::SSPFragment sspf(frag);
 
     unsigned int nTriggers = sspReform.CheckAndGetNTriggers(frag, sspf);
 
     mf::LogInfo("SSPDiagnosticAna") << "Triggers = " << nTriggers;
       
     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
       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=(daqHeader->length-sizeof(SSPDAQ::EventHeader)/sizeof(unsigned int))*2;
 
       //get the information from the header
       try {
         OpChannel = sspReform.GetOpChannel(daqHeader);
 
         FirstSample = sspReform.GetGlobalFirstSample(daqHeader);
         //TimeStamp = ((double)FirstSample)/fNOvAClockFrequency;
 
         double peakSum = sspReform.GetPeakSum(daqHeader);
         double prerise = sspReform.GetBaselineSum(daqHeader);
         double integratedSum = sspReform.GetIntegratedSum(daqHeader);
 
         double PulseAmplitude = -prerise/i2_window*1.0 + peakSum/m1_window*1.0;
 
         int channel = (int) OpChannel; 
 
         // if(channel==94){
         //   if(PulseAmplitude>50&&PulseAmplitude<55)//event with 5 p.e. peak in channel 94
         //     std::cout << "Channel 94 has 5 p.e. peak in event" << evt.event() << std::endl;
           
         //   if(PulseAmplitude>75&&PulseAmplitude<85)//event with 5 p.e. peak in channel 94
         //     std::cout << "Channel 94 has 7 p.e. peak in event" << evt.event() << std::endl;
         // }
 
         double IntegratedCharge = integratedSum - prerise/i2_window*i1_window*1.0;
 
         fPulseAmplitude->Fill(PulseAmplitude);
         mf::LogInfo("SSPDiagnosticAna") << "Pulse Amplitude: " << PulseAmplitude;
         
 
         auto pulse_amplitude_per_channel = PulseAmplitudePerChannel.find( channel);
         if( pulse_amplitude_per_channel == PulseAmplitudePerChannel.end() ) {
           TString histName = TString::Format("pulse_amplitude_channel_%03i", channel);
           TString histTitle = TString::Format("Pulse Amplitude for OP Channel %03i;leading-edge amplitude [ADC]", channel);
           PulseAmplitudePerChannel[channel] =  RunDir.make< TH1D >(histName, histTitle ,125,-20,230);
         }
         
         PulseAmplitudePerChannel[channel]->Fill(PulseAmplitude);
 
         auto integrated_charge_per_channel = IntegratedChargePerChannel.find( channel);
         if( integrated_charge_per_channel == IntegratedChargePerChannel.end() ) {
           TString histName = TString::Format("integrated_charge_channel_%03i", channel);
           TString histTitle = TString::Format("Integrated Charge on OP Channel %03i;integrated charge [ADC*tick]", channel);
           IntegratedChargePerChannel[channel] =  RunDir.make< TH1D >(histName, histTitle ,300,0,3e4);
         }
         
         IntegratedChargePerChannel[channel]->Fill(IntegratedCharge);
 
         auto pulse_amplitude_vs_integrated_charge_per_channel = PulseAmplitudeVsIntegratedChargePerChannel.find( channel);
         if( pulse_amplitude_vs_integrated_charge_per_channel == PulseAmplitudeVsIntegratedChargePerChannel.end() ) {
           TString histName = TString::Format("pulse_amplitude_vs_integrated_charge_channel_%03i", channel);
           TString histTitle = TString::Format("Pulse Amplitude vs. Integrated Charge on OP Channel %03i;integrated charge [ADC*tick];leading-edge amplitude [ADC]", channel);
           PulseAmplitudeVsIntegratedChargePerChannel[channel] =  RunDir.make< TH2D >(histName, histTitle ,300,0,3e4,125,-20,230);
         }
         
         PulseAmplitudeVsIntegratedChargePerChannel[channel]->Fill(IntegratedCharge,PulseAmplitude);
 
         
 
         if (FirstSample < 1e16) {
           sspReform.PrintHeaderInfo(daqHeader);
           mf::LogInfo("SSPDiagnosticAna") << "Problem timestamp at " << FirstSample << std::endl;
           continue;
         }
       } 
       catch (cet::exception const&) {
         continue;
       }
 
       firstTime = std::min(firstTime, FirstSample);
       lastTime  = std::max(lastTime,  FirstSample);
       triggerCount[OpChannel]++;
       
       // Advance the dataPointer to the next header
       dataPointer+=nADC/2;
       
     } // End of loop over triggers
   } // End of loop over fragments (rawFragments)
 
 }

void DAQToOffline::SSPDiagnosticAna::beginEvent ( art::EventNumber_t eventNumber )

private

Definition at line 179 of file SSPDiagnosticAna_module.cc.

 {
   //reset ADC histogram
   //adc_values_->Reset();
   //reset counters
   //n_adc_counter_  = 0;
   //adc_cumulative_ = 0;
 }

void DAQToOffline::SSPDiagnosticAna::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 161 of file SSPDiagnosticAna_module.cc.

 {
   art::ServiceHandle<art::TFileService> tfs;
   //adc_values_ = tfs->make<TH1D>("adc_values","adc_values",4096,-0.5,4095.5);
 
   fPulseAmplitude = tfs->make<TH1D>("pulseamplitude","Pulse Amplitude;leading-edge amplitude [ADC]",125,-50,200);
 
   firstTime = (((unsigned long int)1)<<63);
   lastTime = 0;
 }

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

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 172 of file SSPDiagnosticAna_module.cc.

 {
   PulseAmplitudePerChannel.clear();
   IntegratedChargePerChannel.clear();
   PulseAmplitudeVsIntegratedChargePerChannel.clear();
 }

void DAQToOffline::SSPDiagnosticAna::endEvent ( art::EventNumber_t eventNumber )

private

Definition at line 188 of file SSPDiagnosticAna_module.cc.

 {
   //write the ADC histogram for the given event
   //if(n_adc_counter_)
     //  adc_values_->Write(Form("adc_values:event_%d", eventNumber));
 }

void DAQToOffline::SSPDiagnosticAna::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 348 of file SSPDiagnosticAna_module.cc.

 {
   // Register all vs-run histograms.
   art::ServiceHandle<art::TFileService> tfs;
   for ( int i = 0; i < 100; ++i ) {
     auto histo = PulseAmpVsRun.find( i );
     if ( histo == PulseAmpVsRun.end() ) continue;
     TH2D *newHist = tfs->make< TH2D >( PulseAmpVsRun[i]->GetName(), PulseAmpVsRun[i]->GetTitle(),
                                        PulseAmpVsRun[i]->GetNbinsX(), PulseAmpVsRun[i]->GetXaxis()->GetBinLowEdge(1), PulseAmpVsRun[i]->GetXaxis()->GetBinLowEdge(PulseAmpVsRun[i]->GetNbinsX()+1),
                                        PulseAmpVsRun[i]->GetNbinsY(), PulseAmpVsRun[i]->GetYaxis()->GetBinLowEdge(1), PulseAmpVsRun[i]->GetYaxis()->GetBinLowEdge(PulseAmpVsRun[i]->GetNbinsY()+1) );
     newHist->Add( PulseAmpVsRun[i] );
   }
 }

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

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 195 of file SSPDiagnosticAna_module.cc.

 {
   //delete adc_values_;
   
   long int deltaT = lastTime-firstTime;
   double deltaTus =  ((double)deltaT)/sspReform.ClockFrequency();
 
 
   mf::LogVerbatim("SSPDiagnosticAna") << "!! Diagnostic Rate Report." << std::endl;
   mf::LogVerbatim("SSPDiagnosticAna") << "!! Time: " << deltaTus / 60.e6 << " minutes." << std::endl;
 
   for (auto itr = triggerCount.begin(); itr != triggerCount.end(); itr++) {
     double freq = ((double)itr->second) / deltaTus * 1000.;
     mf::LogVerbatim("SSPDiagnosticAna") << "!!    Channel " << std::setw(3) << itr->first << ": " << freq << " kHz" << std::endl;
   }
   
   
   /*
   mf::LogInfo("SSPDiagnosticAna") << "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";
   */
 
   // for (auto iter = averageWaveforms.begin(); iter != averageWaveforms.end(); iter++)
   //   {
   //     mf::LogInfo("Scaling down channel ") << iter->first << " by 1./" << waveformCount[iter->first] << std::endl;
   //     iter->second->Scale(1./waveformCount[iter->first]);
   //   }
 
   // Format output plots
   for ( int i = 0; i < 100; ++i ) {
     if ( PulseAmplitudePerChannel.find(i) == PulseAmplitudePerChannel.end() ) continue;
     PulseAmplitudePerChannel[i]->SetStats(false);
     PulseAmplitudePerChannel[i]->GetXaxis()->SetTitleSize(0.045);
     IntegratedChargePerChannel[i]->SetStats(false);
     IntegratedChargePerChannel[i]->GetXaxis()->SetTitleSize(0.045);
     PulseAmplitudeVsIntegratedChargePerChannel[i]->SetStats(false);
     PulseAmplitudeVsIntegratedChargePerChannel[i]->SetOption("colz");
     PulseAmplitudeVsIntegratedChargePerChannel[i]->GetXaxis()->SetTitleSize(0.045);
     PulseAmplitudeVsIntegratedChargePerChannel[i]->GetYaxis()->SetTitleSize(0.045);
     averageWaveforms[i]->SetStats(false);
     averageWaveforms[i]->SetOption("colz");
     averageWaveforms[i]->GetXaxis()->SetTitleSize(0.045);
     averageWaveforms[i]->GetYaxis()->SetTitleSize(0.045);
   }
 
   // Analyze plots
   for ( int i = 0; i < 100; ++i ) {
     if ( PulseAmplitudePerChannel.find(i) == PulseAmplitudePerChannel.end() ) continue;
     double ADCperPE(0);
     double INTperPE(0);
     { // Calculate leading-edge amplitude per photoelectron
       int nDiffs(0);
       int nPeaks = specAnalyzer->Search( PulseAmplitudePerChannel[i], 1.5/*sigma*/, "", 0.001 );
       auto *peaks = specAnalyzer->GetPositionX(); // ordered by height!
       std::vector<float> peakVec;
       for ( int p = 0; p < nPeaks; ++p ) peakVec.push_back(peaks[p]);
       std::sort(peakVec.begin(),peakVec.end());
       for ( int p = 1; p < nPeaks-1; ++p ) {
         double diff = peakVec.at(p+1) - peakVec.at(p);
         if ( diff > 20 || diff < 10 ) {
           //std::cout << "Caution! Outlier peak difference at " << p << "/" << p+1 << ". This pair will be skipped." << std::endl;
           continue;
         }
         ADCperPE += diff;
         nDiffs++;
       }
       ADCperPE /= double(nDiffs);
     }
     { // Calculate integrated charge per photoelectron
       int nDiffs(0);
       int nPeaks = specAnalyzer->Search( IntegratedChargePerChannel[i], 2.5/*sigma*/, "", 0.001 );
       auto *peaks = specAnalyzer->GetPositionX(); // ordered by height!
       std::vector<float> peakVec;
       for ( int p = 0; p < nPeaks; ++p ) peakVec.push_back(peaks[p]);
       std::sort(peakVec.begin(),peakVec.end());
       for ( int p = 1; p < nPeaks-1; ++p ) {
         double diff = peakVec.at(p+1) - peakVec.at(p);
         if ( diff > 1800 || diff < 1000 ) {
           //std::cout << "Caution! Outlier peak difference at " << p << "/" << p+1 << ". This pair will be skipped." << std::endl;
           continue;
         }
         INTperPE += diff;
         nDiffs++;
       }
       INTperPE /= double(nDiffs);
     }
     std::cout << "OpDet Channel " << i << " :\t LE " << ADCperPE << " ADC/PE" << "\t IC " << INTperPE << " charge/PE" << std::endl;
   }
 
   // FFT of average waveforms
   art::ServiceHandle<art::TFileService> tfs;
   art::TFileDirectory RunDir = tfs->mkdir( TString::Format("r%03i", run.run()).Data(), "SSP Diagnostics by Run" );
   for ( int i = 0; i < 100; ++i ) {
     if ( PulseAmplitudePerChannel.find(i) == PulseAmplitudePerChannel.end() ) continue;
     TProfile *profile = averageWaveforms[i]->ProfileX();
     auto fft = waveformFFTs.find( i );
     if ( fft == waveformFFTs.end() ) {
       TString histName = TString::Format("waveformFFT_channel_%03i", i);
       TString histTitle = TString::Format("Average Waveform FFT for OP Channel %03i;f (MHz);power", i);
       double dt = averageWaveforms[i]->GetXaxis()->GetBinLowEdge(2) - averageWaveforms[i]->GetXaxis()->GetBinLowEdge(1);
       double Fmax = 1. / (2*dt);
       waveformFFTs[i] = RunDir.make< TH1D >(histName, histTitle, averageWaveforms[i]->GetNbinsX()/2, 0., Fmax);
     }
     profile->FFT( waveformFFTs[i], "MAG" );
     waveformFFTs[i]->SetStats(false);
     waveformFFTs[i]->GetXaxis()->SetTitleSize(0.045);
     waveformFFTs[i]->GetYaxis()->SetTitleSize(0.045);
     delete profile;
   }
 
   // Add pulse amplitude histos to pulse amp vs run histograms
   for ( int i = 0; i < 100; ++i ) {
     if ( PulseAmplitudePerChannel.find(i) == PulseAmplitudePerChannel.end() ) continue;
     auto histo = PulseAmpVsRun.find( i );
     // Make a new histogram if it doesn't already exist.
     if ( histo == PulseAmpVsRun.end() ) {
       TString histName = TString::Format("PulseAmpDistVsRun_channel_%03i", i);
       TString histTitle = TString::Format("Pulse Amplitude Distribution vs Run Number for OP Channel %03i;run number;leading-edge amplitude [ADC]", i);
       PulseAmpVsRun[i] = new TH2D( histName, histTitle, 1, run.run(), run.run()+1, 125, -20, 230 );
     } else {
     // Remake the histogram with extended range
       TH2D* oldHist = PulseAmpVsRun[i];
       int firstRun = (int)std::min( double(run.run()), oldHist->GetXaxis()->GetBinLowEdge(1) );
       int lastRun = (int)std::max( double(run.run()), oldHist->GetXaxis()->GetBinLowEdge( oldHist->GetNbinsX() ) );
       PulseAmpVsRun[i] = new TH2D( oldHist->GetName(), oldHist->GetTitle(), lastRun-firstRun+1, firstRun, lastRun+1, 125, -20, 230 );
       for ( int binX = 0; binX <= oldHist->GetNbinsX(); ++binX ) {
         for ( int binY = 0; binY <= oldHist->GetNbinsY(); ++binY ) {
           double oldVal = oldHist->GetBinContent( binX, binY );
           int runNum = oldHist->GetXaxis()->GetBinLowEdge( binX );
           int targBin = PulseAmpVsRun[i]->GetXaxis()->FindBin( runNum );
           std::cout << "Transferring into new histogram: " << targBin << " " << binY << " " << oldVal << std::endl;
           PulseAmpVsRun[i]->SetBinContent( targBin, binY, oldVal );
         }
       }
       delete oldHist;
     }
     // Add the new data
     for ( int binY = 0; binY <= PulseAmplitudePerChannel[i]->GetNbinsX(); ++binY ) {
       double val = PulseAmplitudePerChannel[i]->GetBinContent(binY);
       double amp = PulseAmplitudePerChannel[i]->GetBinCenter(binY);
       std::cout << "Filling histogram with " << run.run() << " " << amp << " " << val << "( PulseAmplitudePerChannel bin " << binY << " )" << std::endl;
       PulseAmpVsRun[i]->Fill( run.run(), amp, val );
     }
     
   }
 
 }

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

delete

SSPDiagnosticAna& DAQToOffline::SSPDiagnosticAna::operator= ( SSPDiagnosticAna && )

delete

void DAQToOffline::SSPDiagnosticAna::printParameterSet ( )

Definition at line 151 of file SSPDiagnosticAna_module.cc.

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

void DAQToOffline::SSPDiagnosticAna::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 129 of file SSPDiagnosticAna_module.cc.

                                                                            {
 
   fFragType           = pset.get<std::string>("FragType");
   fRawDataLabel       = pset.get<std::string>("RawDataLabel");
   fInputModule        = pset.get<std::string>("InputModule");
   fInputLabel         = pset.get<std::string>("InputLabel");
   // m1_window           = pset.get<int>("SSPm1");
   // i1_window           = pset.get<int>("SSPi1");
   // i2_window           = pset.get<int>("SSPi2");
 
   //fDebug = pset.get<bool>("Debug");
   //fZeroThreshold=0;
   //fCompression=raw::kNone;
   // Obtain parameters from TimeService
   auto const* timeService = lar::providerFrom<detinfo::DetectorClocksService>();
   fSampleFreq = timeService->OpticalClock().Frequency();
 
 
   printParameterSet();
   
 }