Inheritance diagram for opdet::AverageWaveform:

Public Member Functions
	AverageWaveform (fhicl::ParameterSet const &)

virtual	~AverageWaveform ()

void	analyze (art::Event const &) override

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	beginJob () override

void	endJob () override

Private Attributes
std::string	fInputModule

std::string	fInstanceName

double	fSampleFreq

double	fTimeBegin

int	fBaselineSubtract

int	fNticks

std::map< int, TH1D * >	averageWaveforms

std::map< int, int >	waveformCount

TH1D *	averageWaveformAll

int	eventCount

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 43 of file AverageWaveform_module.cc.

Constructor & Destructor Documentation

opdet::AverageWaveform::AverageWaveform ( fhicl::ParameterSet const & pset )

Definition at line 89 of file AverageWaveform_module.cc.

         : EDAnalyzer(pset)
         , fInputModule      { pset.get< std::string >("InputModule") }
         , fInstanceName     { pset.get< std::string >("InstanceName") }
         , fTimeBegin        { 0 }
         , fBaselineSubtract { pset.get< int >("BaselineSubtract", 0) }
         , fNticks           { pset.get< int >("Nticks", 0) }
         , averageWaveformAll{ NULL }
         , eventCount        { 0 }
     {
 
         // Obtain parameters from TimeService
         auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
         fSampleFreq = clockData.OpticalClock().Frequency();
 
         // Assume starting at 0
         fTimeBegin  = 0;
 
     }

opdet::AverageWaveform::~AverageWaveform ( )

virtual

Definition at line 111 of file AverageWaveform_module.cc.

112 {

113 }

Member Function Documentation

void opdet::AverageWaveform::analyze ( art::Event const & evt )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 137 of file AverageWaveform_module.cc.

     {
 
         // Get OpDetWaveforms from the event
         art::InputTag itag1(fInputModule, fInstanceName);
         auto waveformHandle = evt.getHandle< std::vector< raw::OpDetWaveform > >(itag1);
 
         // Access ART's TFileService, which will handle creating and writing
         // histograms for us
         art::ServiceHandle< art::TFileService > tfs;
 
         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();
 
             if (fNticks == 0) fNticks = pulse.size();
 
 
             // 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] =  tfs->make< TH1D >(histName, ";t (us);", fNticks, 0, float(fNticks+1) / fSampleFreq);
             }
             if (!averageWaveformAll) {
                 averageWaveformAll =  tfs->make< TH1D >("avgwaveform_channel_all", ";t (us);", fNticks, 0, float(fNticks+1) / fSampleFreq);
             }
 
             // Add this waveform to this histogram
             for (size_t tick = 0; tick < pulse.size(); tick++) {
                 averageWaveforms[channel]->Fill(double(tick)/fSampleFreq, pulse[tick] - fBaselineSubtract);
                 averageWaveformAll->Fill(double(tick)/fSampleFreq, pulse[tick] - fBaselineSubtract);
             }
 
             // Count number of waveforms on each channel
             waveformCount[channel]++;
         }
         eventCount++;
 
     }

void opdet::AverageWaveform::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 117 of file AverageWaveform_module.cc.

118 {

119 }

void opdet::AverageWaveform::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 122 of file AverageWaveform_module.cc.

     {
         double waveformCountAll = 0;
         for (auto iter = averageWaveforms.begin(); iter != averageWaveforms.end(); iter++)
         {
             waveformCountAll += waveformCount[iter->first];
             mf::LogInfo("AverageWaveform") << "Scaling down channel " << iter->first << " by 1./" << waveformCount[iter->first] << std::endl;
             iter->second->Scale(1./waveformCount[iter->first]);
         }
         mf::LogInfo("AverageWaveform") << "Scaling down all channels by 1./" << eventCount << std::endl;
         averageWaveformAll->Scale(1./eventCount);
     }

Member Data Documentation

TH1D* opdet::AverageWaveform::averageWaveformAll

private

Definition at line 69 of file AverageWaveform_module.cc.

std::map< int, TH1D* > opdet::AverageWaveform::averageWaveforms

private

Definition at line 67 of file AverageWaveform_module.cc.

int opdet::AverageWaveform::eventCount

private

Definition at line 70 of file AverageWaveform_module.cc.

int opdet::AverageWaveform::fBaselineSubtract

private

Definition at line 63 of file AverageWaveform_module.cc.

std::string opdet::AverageWaveform::fInputModule

private

Definition at line 59 of file AverageWaveform_module.cc.

std::string opdet::AverageWaveform::fInstanceName

private

Definition at line 60 of file AverageWaveform_module.cc.

int opdet::AverageWaveform::fNticks

private

Definition at line 64 of file AverageWaveform_module.cc.

double opdet::AverageWaveform::fSampleFreq

private

Definition at line 61 of file AverageWaveform_module.cc.

double opdet::AverageWaveform::fTimeBegin

private

Definition at line 62 of file AverageWaveform_module.cc.

std::map< int, int > opdet::AverageWaveform::waveformCount

private

Definition at line 68 of file AverageWaveform_module.cc.

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

duneopdet/duneopdet/OpticalDetector/AverageWaveform_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation