Inheritance diagram for WaveformPDSAndTruthDump:

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

	WaveformPDSAndTruthDump (WaveformPDSAndTruthDump const &)=delete

	WaveformPDSAndTruthDump (WaveformPDSAndTruthDump &&)=delete

WaveformPDSAndTruthDump &	operator= (WaveformPDSAndTruthDump const &)=delete

WaveformPDSAndTruthDump &	operator= (WaveformPDSAndTruthDump &&)=delete

void	analyze (art::Event const &e) 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 Attributes
std::string	m_inputTagGEANT

std::string	m_inputTagPDS

std::string	m_outputFilename_pds

std::string	m_outputFilename_true_pds

std::ofstream	m_outputFile_pds

std::ofstream	m_outputFile_true_pds

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 29 of file WaveformPDSAndTruthDump_module.cc.

Constructor & Destructor Documentation

WaveformPDSAndTruthDump::WaveformPDSAndTruthDump ( fhicl::ParameterSet const & p )

explicit

Definition at line 55 of file WaveformPDSAndTruthDump_module.cc.

   : EDAnalyzer(p),
     m_inputTagGEANT (p.get<std::string>("InputTagGEANT" , "largeant")),
     m_inputTagPDS   (p.get<std::string>("InputTagPDS"   , "opdigi"  )), 
     m_outputFilename_pds     (p.get<std::string>("OutputFilePDS"    ,"OutputFilePDS.txt"    )),
     m_outputFilename_true_pds(p.get<std::string>("OutputFileTruePDS","OutputFileTruePDS.txt")),
     m_outputFile_pds     (m_outputFilename_pds     ),
     m_outputFile_true_pds(m_outputFilename_true_pds)
 {
 }

WaveformPDSAndTruthDump::WaveformPDSAndTruthDump ( WaveformPDSAndTruthDump const & )

delete

WaveformPDSAndTruthDump::WaveformPDSAndTruthDump ( WaveformPDSAndTruthDump && )

delete

Member Function Documentation

void WaveformPDSAndTruthDump::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 66 of file WaveformPDSAndTruthDump_module.cc.

 {
   art::ServiceHandle<geo::Geometry> geo;
   auto const& digits_handle_pds=e.getValidHandle<std::vector<raw::OpDetWaveform>>(m_inputTagPDS);
   std::map<int, std::pair<int, int>> channel_timestamp;
   
   auto& digits_pds_in = *digits_handle_pds;
   for (auto&& digit: digits_pds_in) {
     m_outputFile_pds << e.event() << " "
                      << digit.ChannelNumber() << " ";
                      
     auto f = channel_timestamp.find(digit.ChannelNumber());
     if (f == channel_timestamp.end()) {
       channel_timestamp[digit.ChannelNumber()] = std::make_pair(digit.TimeStamp(), digit.size());
     } else {
       int begin_old = f->second.first;
       int begin_new = digit.TimeStamp();
       begin_new = (begin_old > begin_new) ? begin_new : begin_old;
       int end_old = f->second.first+f->second.second;
       int end_new = digit.TimeStamp()+digit.size();
       end_new = (end_old > end_new) ? end_old : end_new;
       f->second.first = begin_new;
       f->second.second = end_new - begin_new;
     }
     for(auto const& adc: digit) {
       m_outputFile_pds << adc << " ";
     }
     
     m_outputFile_pds << "\n";
   }
 
   // PDS truth
   auto const& truth_handle_pds=e.getValidHandle<std::vector<sim::OpDetBacktrackerRecord>>(m_inputTagGEANT);
   auto& truth_pds_in =*truth_handle_pds;
 
   for (auto const& interesting: channel_timestamp) {
     int channel = interesting.first;
     int timestamp = interesting.second.first;
     int nticks = interesting.second.second;
     bool found = false;
     for (auto&& truth: truth_pds_in) {
       if (truth.OpDetNum() != channel) continue;
       m_outputFile_true_pds << e.event() << " "
                             << truth.OpDetNum() << " ";
       for(int ichge=timestamp; ichge<nticks+timestamp; ++ichge){
         m_outputFile_true_pds << truth.Photons(ichge) << " ";
       }
       m_outputFile_true_pds << "\n";
       found=true;
       break;
     }
     if (!found) {
       std::cout << "didn't find this channel " << channel
                 << " timestamp " << timestamp
                 << " ticks " << nticks <<" in truth\n";
     }
   }
 }