ana::SBNOsc::OpDetWaveformMaker Class Reference

Electron neutrino event selection. More...

Inheritance diagram for ana::SBNOsc::OpDetWaveformMaker:

Public Member Functions
	OpDetWaveformMaker ()
void	Initialize (fhicl::ParameterSet *config=NULL)
void	Finalize ()
bool	ProcessEvent (const gallery::Event &ev, const std::vector< event::Interaction > &truth, std::vector< event::RecoInteraction > &reco)
Public Member Functions inherited from core::SelectionBase
	SelectionBase ()
virtual	~SelectionBase ()
Public Member Functions inherited from core::ProcessorBase
	ProcessorBase ()
virtual	~ProcessorBase ()
virtual void	FillTree ()
virtual void	FillRecoTree ()
virtual void	EventCleanup ()
template<class T >
TBranch *	AddBranch (std::string name, T *obj)
template<class T >
TBranch *	AddRecoBranch (std::string name, T *obj)

Protected Attributes
std::string	fOpDetWaveformTag
float	clock_freq
unsigned	event_ind
Protected Attributes inherited from core::ProcessorBase
unsigned long	fEventIndex
	An incrementing index. More...
Experiment	fExperimentID
	Experiment identifier. More...
ProviderManager *	fProviderManager
	Interface for provider access. More...
std::string	fOutputFilename
	The output filename. More...
std::string	fProviderConfig
	A custom provider config fcl file. More...
std::vector< geo::BoxBoundedGeo >	fActiveVolumes
	List of active volumes in configured detector. More...
bool	fWriteTree
	Enable writing of the main tree. More...
TFile *	fOutputFile
	The output ROOT file. More...
TTree *	fTree
	The output ROOT tree. More...
event::Event *	fEvent
	The standard output event data structure. More...
bool	fWriteRecoTree
	Enable writing of the reco tree. More...
TTree *	fRecoTree
	The output reco ROOT tree. More...
event::RecoEvent *	fRecoEvent
	The standard output reco event data structure. More...
TTree *	fSubRunTree
	Subrun output tree. More...
SubRun *	fSubRun
	Standard output subrun structure. More...
TTree *	fFileMetaTree
	File metadata output tree. More...
FileMeta *	fFileMeta
	standard output file metadata structure More...
TParameter< int > *	fExperimentParameter
	Saves value of experiment enum. More...
std::set< std::pair< int, int > >	fSubRunCache
	Cache stored subruns. More...
art::InputTag	fTruthTag
	art tag for MCTruth information More...
art::InputTag	fFluxTag
	art tag for MCFlux information More...
std::vector< art::InputTag >	fWeightTags
	art tag(s) for MCEventWeight information More...
art::InputTag	fMCTrackTag
	art tag for MCTrack More...
art::InputTag	fMCShowerTag
	art tag for MCShower More...
art::InputTag	fMCParticleTag
	art tag for MCParticle More...
std::string	fGeneratorProcess
	process_name of process used to run genie. Used to extract subrun/POT information More...

Additional Inherited Members
Public Attributes inherited from core::ProcessorBase
std::vector < event::RecoInteraction > *	fReco
	Reco interaction list. More...
Protected Member Functions inherited from core::ProcessorBase
virtual void	Initialize (char *config=NULL)
virtual void	Setup (char *config=NULL)
virtual void	Setup (fhicl::ParameterSet *config=NULL)
virtual void	Teardown ()
void	BuildEventTree (gallery::Event &ev)
void	SetupServices (gallery::Event &ev)
void	UpdateSubRuns (gallery::Event &ev)
void	UpdateFileMeta (gallery::Event &ev)

Detailed Description

Electron neutrino event selection.

Definition at line 31 of file OpDetWaveformMaker.cc.

Constructor & Destructor Documentation

ana::SBNOsc::OpDetWaveformMaker::OpDetWaveformMaker ( )

inline

Constructor.

Definition at line 34 of file OpDetWaveformMaker.cc.

{}

Member Function Documentation

void ana::SBNOsc::OpDetWaveformMaker::Finalize ( )

inlinevirtual

Finalize and write objects to the output file.

Implements core::ProcessorBase.

Definition at line 49 of file OpDetWaveformMaker.cc.

                  {
  }

void ana::SBNOsc::OpDetWaveformMaker::Initialize ( fhicl::ParameterSet *  config = NULL )

inlinevirtual

Initialization.

Parameters

config

A configuration, as a FHiCL ParameterSet object

Implements core::ProcessorBase.

Definition at line 41 of file OpDetWaveformMaker.cc.

                                                  {
    fOpDetWaveformTag = config ? config->get<std::string>("OpDetWaveformTag", "opdaq") : "opdaq";
    clock_freq = fProviderManager->GetDetectorClocksProvider()->OpticalClock().Frequency();
    std::cout << "CLOCK FREQ: " << clock_freq << std::endl;
    event_ind = 0;
  }

bool ana::SBNOsc::OpDetWaveformMaker::ProcessEvent ( const gallery::Event &  ev, const std::vector< event::Interaction > &  truth, std::vector< event::RecoInteraction > &  reco  )

inlinevirtual

Process one event.

Parameters

ev	A single event, as a gallery::Event
Reconstructed	interactions

Returns

True to keep event

Implements core::ProcessorBase.

Definition at line 59 of file OpDetWaveformMaker.cc.

                                                                                                                             {
    const std::vector<raw::OpDetWaveform> &waveforms = *ev.getValidHandle<std::vector<raw::OpDetWaveform>>(fOpDetWaveformTag);
    //std::map<raw::Channel_t, std::pair<std::vector<float>, std::vector<float>>> waveform_plots;
    std::map<raw::Channel_t, std::vector<std::pair<float, float>>> waveform_data;
    std::map<raw::Channel_t, std::vector<float>> wf_start_times;
    std::map<raw::Channel_t, std::vector<float>> wf_periods;
    for (const raw::OpDetWaveform &wfconst: waveforms) {
      raw::OpDetWaveform wf = wfconst;
      //if (waveform_plots.count(wf.ChannelNumber()) > 0) {
      //  std::cout << "Multi-waveform output! Channel: " << wf.ChannelNumber() << " Event: " << event_ind << std::endl;
      //}
      for (unsigned i = 0; i < wf.Waveform().size(); i++) {
        waveform_data[wf.ChannelNumber()].push_back({i / clock_freq + wf.TimeStamp(), (float)wf.Waveform()[i]});
      }
      wf_start_times[wf.ChannelNumber()].push_back(wf.TimeStamp());
      wf_periods[wf.ChannelNumber()].push_back( wf.Waveform().size() / clock_freq);
    }

    // now sort each one
    //for (auto &wf_data_pair: waveform_data) {
    //  std::sort(wf_data_pair.second.begin(), wf_data_pair.second.end(), [](auto lhs, auto rhs) {return lhs.first < rhs.first;});
    //}
    std::map<raw::Channel_t, std::pair<std::vector<float>, std::vector<float>>> waveform_plots;
    for (auto &wf_data_pair: waveform_data) {
      for (auto data_pair: wf_data_pair.second) {
        waveform_plots[wf_data_pair.first].first.push_back(data_pair.first);
        waveform_plots[wf_data_pair.first].second.push_back(data_pair.second);
      }
    }

    for (auto &start_time_pair: wf_start_times) {
      std::cout << "Event: " << event_ind << " Channel: " << start_time_pair.first << std::endl;
      const std::vector<float> &wf_period = wf_periods.at(start_time_pair.first); 
      for (unsigned i = 0; i < start_time_pair.second.size(); i++) {
        std::cout << start_time_pair.second[i] << " (" << wf_period[i] <<") "; 

      }
      std::cout << std::endl;
    }

    std::vector<TGraph *> plots;
    for (auto &plot_info: waveform_plots) {
      TGraph *plot = new TGraph(plot_info.second.first.size(), &plot_info.second.first[0], &plot_info.second.second[0]);
      plot->SetName((std::string("Waveform ch: ") + std::to_string(plot_info.first) + " ev: " + std::to_string(event_ind)).c_str());
      plot->SetTitle((std::string("Waveform ch: ") + std::to_string(plot_info.first) + " ev: " + std::to_string(event_ind)).c_str());
      plots.push_back(plot);
    }
    fOutputFile->cd();
    for (TGraph *plot: plots) {
      plot->Write();
      delete plot;
    }
    event_ind += 1;

    return false; 
  }

Member Data Documentation

float ana::SBNOsc::OpDetWaveformMaker::clock_freq

protected

Definition at line 118 of file OpDetWaveformMaker.cc.

unsigned ana::SBNOsc::OpDetWaveformMaker::event_ind

protected

Definition at line 119 of file OpDetWaveformMaker.cc.

std::string ana::SBNOsc::OpDetWaveformMaker::fOpDetWaveformTag

protected

Definition at line 117 of file OpDetWaveformMaker.cc.

---

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

• OpDetWaveformMaker.cc