#include <RawDigitMixer.h>

Classes
struct	RD_Info

Public Member Functions
	RawDigitMixer (bool p=false)
	Default constructor. More...

void	DeclareData (std::vector< raw::RawDigit > const &dataVector)

void	Mix (std::vector< raw::RawDigit > const &mcVector, std::unordered_map< raw::ChannelID_t, float > const &map)

void	FillRawDigitOutput (std::vector< raw::RawDigit > &output)

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

void	SetStuckBitRetentionMethod (bool b)

void	SetDataMixTicks (size_t start, size_t end)

void	SetMCMixTicks (size_t start, size_t end)

Private Attributes
bool	_stuckRetention

bool	_printWarnings

size_t	_datastart

size_t	_dataend

size_t	_mcstart

size_t	_mcend

std::vector< RD_Info >	fOutputWaveforms

std::unordered_map< raw::ChannelID_t, size_t >	fChannelIndexMap

RawDigitAdder_35t	fRDAdderAlg

Detailed Description

Definition at line 42 of file RawDigitMixer.h.

Constructor & Destructor Documentation

mix::RawDigitMixer::RawDigitMixer ( bool p = false )

inline

Default constructor.

Definition at line 47 of file RawDigitMixer.h.

47 :

48 _printWarnings(p){};

test.p

p

Definition: test.py:223

mix::RawDigitMixer::_printWarnings

bool _printWarnings

Definition: RawDigitMixer.h:66

virtual mix::RawDigitMixer::~RawDigitMixer ( )

inlinevirtual

Default destructor.

Definition at line 57 of file RawDigitMixer.h.

57 {};

Member Function Documentation

void RawDigitMixer::DeclareData ( std::vector< raw::RawDigit > const & dataVector )

Definition at line 10 of file RawDigitMixer.cxx.

                                                                           {
 
   fChannelIndexMap.clear();
   fOutputWaveforms.clear();
   fOutputWaveforms.resize(dataVector.size());
 
   for(size_t i_rd=0; i_rd<dataVector.size(); i_rd++){
     
     std::vector<short> waveform;
     if(dataVector[i_rd].Compression()!=raw::kNone)
       {
         raw::Uncompress(dataVector[i_rd].ADCs(),
                         waveform,
                         dataVector[i_rd].GetPedestal(),
                         dataVector[i_rd].Compression());
       }
     else 
       {
         waveform = dataVector[i_rd].ADCs();
       }    
 
     if (_datastart > waveform.size() || _dataend > waveform.size())
       throw std::runtime_error("RawDigitMixer::DeclareData : Input data RawDigit shape does not match what is expected.");
 
     fChannelIndexMap[dataVector[i_rd].Channel()] = i_rd;
     
     //initialize adc vector for output
     fOutputWaveforms[i_rd].waveform.resize(dataVector[i_rd].Samples(),0);
     fOutputWaveforms[i_rd].channel = dataVector[i_rd].Channel();
     fOutputWaveforms[i_rd].ped     = dataVector[i_rd].GetPedestal();
     fOutputWaveforms[i_rd].sigma   = dataVector[i_rd].GetSigma();
     
     
     //do the steps for filling that output vector
     //This is data, so set scales to one.
     fRDAdderAlg.SetScaleInputs(1.0,1.0);
     fRDAdderAlg.SetPedestalInputs(0.0,0.0);
     fRDAdderAlg.AddRawDigits(waveform,fOutputWaveforms[i_rd].waveform);
   }
   
 }

void RawDigitMixer::FillRawDigitOutput ( std::vector< raw::RawDigit > & output )

Definition at line 84 of file RawDigitMixer.cxx.

                                                                         {
 
   for(auto const& rd_i : fOutputWaveforms){
     
     //now emplace back onto output collection...
     output.emplace_back(rd_i.channel,
                         rd_i.waveform.size(),
                         rd_i.waveform);
     
     //set pedestal and rms to be same as data
     output.back().SetPedestal(rd_i.ped,rd_i.sigma);
   }
   
 }

void RawDigitMixer::Mix	(	std::vector< raw::RawDigit > const &	mcVector,
		std::unordered_map< raw::ChannelID_t, float > const &	map
	)

Definition at line 52 of file RawDigitMixer.cxx.

                                                                                      {
 
   for( auto const& rd : mcVector){
 
     if(rd.Compression()!=raw::kNone)
       throw std::runtime_error("Error in RawDigitMixer::Mix : Compressed MC waveforms not supported. Turn it off and try again."); 
 
     if (_mcstart > rd.Samples() || _mcend > rd.Samples())
       throw std::runtime_error("RawDigitMixer::Mix : Input MC RawDigit shape does not match what is expected.");
 
     auto it_ch = fChannelIndexMap.find(rd.Channel());
 
     //if this channel is not in the data, skip this channel!
     if(it_ch==fChannelIndexMap.end())
       continue;
 
     size_t i_output = it_ch->second;
 
     fRDAdderAlg.SetPedestalInputs(rd.GetPedestal(),0.0);
     fRDAdderAlg.SetScaleInputs(scale_map.at(rd.Channel()),1.0);
     fRDAdderAlg.SetStuckBitRetentionMethod(_stuckRetention);
     
     std::vector<short> const data_trimmed(fOutputWaveforms[i_output].waveform.begin()+_datastart,
                                           fOutputWaveforms[i_output].waveform.begin()+_dataend);
     std::vector<short> const mc_trimmed(rd.ADCs().begin()+_mcstart,
                                         rd.ADCs().begin()+_mcend);
     fRDAdderAlg.AddRawDigits(mc_trimmed,data_trimmed,fOutputWaveforms[i_output].waveform);
   }
 
 }