#include <RawDigitAdderAna.hh>

Public Member Functions
	RawDigitAdderAna (size_t sample=100, std::vector< raw::ChannelID_t > special_channels=std::vector< raw::ChannelID_t >(0), bool print_bad=true, std::string in1label="in1", std::string in2label="in2", std::string sumlabel="sum")
	Default constructor. More...

virtual	~RawDigitAdderAna ()
	Default destructor. More...

size_t	CheckOverlay (std::vector< raw::RawDigit > const &in1, std::vector< raw::RawDigit > const &in2, std::vector< raw::RawDigit > const &sum)

void	CreateOutputHistograms (std::vector< TH1S * > const histoPtrVector, std::vector< raw::RawDigit > const &in1, std::vector< raw::RawDigit > const &in2, std::vector< raw::RawDigit > const &sum, unsigned int run, unsigned int event)

Private Member Functions
void	ResetOutput ()

void	CreateHistogram (TH1S *histo, raw::RawDigit const &waveform, unsigned int run, unsigned int event, unsigned int channel, std::string label)

Private Attributes
size_t	fChannelSampleInterval

std::vector< raw::ChannelID_t >	fChannelsSpecial

bool	fPrintBadOverlays

std::string	fInput1Label

std::string	fInput2Label

std::string	fSumLabel

std::vector< raw::ChannelID_t >	fChannelsToPrint

std::vector< raw::ChannelID_t >	fChannelsBadOverlay

Detailed Description

Definition at line 26 of file RawDigitAdderAna.hh.

Constructor & Destructor Documentation

mix::RawDigitAdderAna::RawDigitAdderAna	(	size_t	sample = `100`,
		std::vector< raw::ChannelID_t >	special_channels = `std::vector<raw::ChannelID_t>(0)`,
		bool	print_bad = `true`,
		std::string	in1label = `"in1"`,
		std::string	in2label = `"in2"`,
		std::string	sumlabel = `"sum"`
	)

Default constructor.

Definition at line 11 of file RawDigitAdderAna.cxx.

   : fChannelSampleInterval(sample),
     fChannelsSpecial(special_channels),
     fPrintBadOverlays(print_bad),
     fInput1Label(in1label),
     fInput2Label(in2label),
     fSumLabel(sumlabel)
 {
 }

virtual mix::RawDigitAdderAna::~RawDigitAdderAna ( )

inlinevirtual

Default destructor.

Definition at line 39 of file RawDigitAdderAna.hh.

39 {};

Member Function Documentation

size_t mix::RawDigitAdderAna::CheckOverlay	(	std::vector< raw::RawDigit > const &	in1,
		std::vector< raw::RawDigit > const &	in2,
		std::vector< raw::RawDigit > const &	sum
	)

Definition at line 31 of file RawDigitAdderAna.cxx.

 {
 
   ResetOutput();
   
   //check size of all the rawdigit vectors
   if(in1.size()!=in2.size() ||
      in1.size()!=sum.size() )
     throw std::runtime_error("Error in RawDigitAdderAna::CheckOverlay : Input vector lists not equal size.");
 
   for(size_t i_ch=0; i_ch<sum.size(); i_ch++){
 
     //check to make sure channels are same
     if(in1[i_ch].Channel()!=in2[i_ch].Channel() ||
        in1[i_ch].Channel()!=sum[i_ch].Channel())
       throw std::runtime_error("Error in RawDigitAdderAna::CheckOverlay : Input vector lists out of order.");
     
     //check to make sure channels have same ADC vector
     if(in1[i_ch].ADCs().size()!=in2[i_ch].ADCs().size() ||
        in1[i_ch].ADCs().size()!=sum[i_ch].ADCs().size())
       throw std::runtime_error("Error in RawDigitAdderAna::CheckOverlay : Input channels have different sizes.");
 
     bool print= ( (in1[i_ch].Channel()%fChannelSampleInterval==0) ||
                   (std::find(fChannelsSpecial.begin(),fChannelsSpecial.end(),in1[i_ch].Channel())!=fChannelsSpecial.end()) );
 
     for(size_t i_t=0; i_t<sum[i_ch].ADCs().size(); i_t++){
       if( (in1[i_ch].ADCs()[i_t]+in2[i_ch].ADCs()[i_t])==sum[i_ch].ADCs()[i_t] )
         continue;
 
       fChannelsBadOverlay.emplace_back(in1[i_ch].Channel());
       if(fPrintBadOverlays) print=true;
       break;
     }
 
     if(print) fChannelsToPrint.emplace_back(in1[i_ch].Channel());
     
   }
 
   return fChannelsToPrint.size()*3;
   
 }

void mix::RawDigitAdderAna::CreateHistogram	(	TH1S *	histo,
		raw::RawDigit const &	waveform,
		unsigned int	run,
		unsigned int	event,
		unsigned int	channel,
		std::string	label
	)

private

Definition at line 75 of file RawDigitAdderAna.cxx.

 {
     std::stringstream hname,htitle;
     hname << "h_" << label
           << "_run" << run
           << "_ev" << event
           << "_ch" << channel;
     htitle << "Waveform, Channel " << channel
            << ", Event " << event
            << ", Run " << run
            << ", Input " << label;
     histo->SetName(hname.str().c_str());
     histo->SetTitle(htitle.str().c_str());
     histo->SetBins(waveform.ADCs().size(),0,waveform.ADCs().size());
 
     for(size_t i_t=0; i_t<waveform.ADCs().size(); i_t++)
       histo->SetBinContent(i_t,waveform.ADCs()[i_t]);
     
 }

void mix::RawDigitAdderAna::CreateOutputHistograms	(	std::vector< TH1S * > const	histoPtrVector,
		std::vector< raw::RawDigit > const &	in1,
		std::vector< raw::RawDigit > const &	in2,
		std::vector< raw::RawDigit > const &	sum,
		unsigned int	run,
		unsigned int	event
	)

Definition at line 101 of file RawDigitAdderAna.cxx.

 {
 
   if(histoPtrVector.size()!=fChannelsToPrint.size())
     throw std::runtime_error("Error in RawDigitAdderAna::CreateOutputHistograms : Histogram output vector not equal to channels to print.");
   
 
   //check size of all the rawdigit vectors
   if(in1.size()!=in2.size() ||
      in1.size()!=sum.size() )
     throw std::runtime_error("Error in RawDigitAdderAna::CheckOverlay : Input vector lists not equal size.");
 
   size_t histo_count=0;
   
   for(size_t i_ch=0; i_ch<sum.size(); i_ch++){
     
     //check to make sure channels are same
     if(in1[i_ch].Channel()!=in2[i_ch].Channel() ||
        in1[i_ch].Channel()!=sum[i_ch].Channel())
       throw std::runtime_error("Error in RawDigitAdderAna::CheckOverlay : Input vector lists out of order.");
 
     if(std::find(fChannelsToPrint.begin(),fChannelsToPrint.end(),in1[i_ch].Channel())==fChannelsToPrint.end())
       continue;
     
     //check to make sure channels have same ADC vector
     if(in1[i_ch].ADCs().size()!=in2[i_ch].ADCs().size() ||
        in1[i_ch].ADCs().size()!=sum[i_ch].ADCs().size())
       throw std::runtime_error("Error in RawDigitAdderAna::CheckOverlay : Input channels have different sizes.");
 
     CreateHistogram(histoPtrVector.at(histo_count),
                     in1[i_ch],
                     run,event,in1[i_ch].Channel(),
                     fInput1Label);
     histo_count++;
     CreateHistogram(histoPtrVector.at(histo_count),
                     in2[i_ch],
                     run,event,in2[i_ch].Channel(),
                     fInput2Label);
     histo_count++;
     CreateHistogram(histoPtrVector.at(histo_count),
                     sum[i_ch],
                     run,event,sum[i_ch].Channel(),
                     fSumLabel);
     histo_count++;
     
   }
 
 }

void mix::RawDigitAdderAna::ResetOutput ( )

private

Definition at line 26 of file RawDigitAdderAna.cxx.

                                      {
   fChannelsToPrint.clear();
   fChannelsBadOverlay.clear();
 }

Member Data Documentation

size_t mix::RawDigitAdderAna::fChannelSampleInterval

private

Definition at line 53 of file RawDigitAdderAna.hh.

std::vector<raw::ChannelID_t> mix::RawDigitAdderAna::fChannelsBadOverlay

private

Definition at line 60 of file RawDigitAdderAna.hh.

std::vector<raw::ChannelID_t> mix::RawDigitAdderAna::fChannelsSpecial

private

Definition at line 54 of file RawDigitAdderAna.hh.

std::vector<raw::ChannelID_t> mix::RawDigitAdderAna::fChannelsToPrint

private

Definition at line 59 of file RawDigitAdderAna.hh.

std::string mix::RawDigitAdderAna::fInput1Label

private

Definition at line 56 of file RawDigitAdderAna.hh.

std::string mix::RawDigitAdderAna::fInput2Label

private

Definition at line 57 of file RawDigitAdderAna.hh.

bool mix::RawDigitAdderAna::fPrintBadOverlays

private

Definition at line 55 of file RawDigitAdderAna.hh.

std::string mix::RawDigitAdderAna::fSumLabel

private

Definition at line 58 of file RawDigitAdderAna.hh.

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

dunetpc/dune/DataOverlay/DataOverlay/RawDigitAdderAna.hh
dunetpc/dune/DataOverlay/DataOverlay/RawDigitAdderAna.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation