Class for the decoding of the SIS3320 sampling ADC data. More...

#include <QwSIS3320_Channel.h>

Inheritance diagram for QwSIS3320_Channel:

Collaboration diagram for QwSIS3320_Channel:

Public Member Functions
	QwSIS3320_Channel (UInt_t channel=0, TString name="auto")

	QwSIS3320_Channel (const QwSIS3320_Channel &source)

virtual	~QwSIS3320_Channel ()

void	AddLogicalAccumulator (const TString name, const std::vector< TString > accums, const std::vector< Double_t > weights)

void	InitializeChannel (UInt_t channel, TString name)

void	ClearEventData ()

void	RandomizeEventData (int helicity=0, double time=0.0)
	Internally generate random event data. More...

void	EncodeEventData (std::vector< UInt_t > &buffer)
	Encode the event data into a CODA buffer. More...

Int_t	ProcessEvBuffer (UInt_t *buffer, UInt_t num_words_left, UInt_t index=0)

void	ProcessEvent ()

const QwSIS3320_Channel	operator+ (const Double_t &value) const

const QwSIS3320_Channel	operator- (const Double_t &value) const

const QwSIS3320_Channel	operator+ (const QwSIS3320_Channel &value) const

const QwSIS3320_Channel	operator- (const QwSIS3320_Channel &value) const

QwSIS3320_Channel &	operator= (const QwSIS3320_Channel &value)

QwSIS3320_Channel &	operator+= (const Double_t &value)

QwSIS3320_Channel &	operator-= (const Double_t &value)

QwSIS3320_Channel &	operator+= (const QwSIS3320_Channel &value)

QwSIS3320_Channel &	operator-= (const QwSIS3320_Channel &value)

void	Sum (QwSIS3320_Channel &value1, QwSIS3320_Channel &value2)

void	Difference (QwSIS3320_Channel &value1, QwSIS3320_Channel &value2)

void	Ratio (QwSIS3320_Channel &numer, QwSIS3320_Channel &denom)

void	Offset (Double_t Offset)

void	Scale (Double_t Offset)

void	ConstructHistograms (TDirectory *folder, TString &prefix)
	Construct the histograms for this data element. More...

void	FillHistograms ()
	Fill the histograms for this data element. More...

void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)

void	FillTreeVector (std::vector< Double_t > &values) const

QwSIS3320_Samples &	GetSamples (size_t i)

QwSIS3320_Samples &	GetSamplesRaw (size_t i)

size_t	GetNumberOfEvents () const

void	SetNumberOfEvents (UInt_t nevents)

size_t	GetNumberOfAccumulators () const

void	SetNumberOfAccumulators (UInt_t naccumulators)

void	SetPedestal (const Double_t ped)

Double_t	GetPedestal () const

void	SetCalibrationFactor (const Double_t factor)

Double_t	GetCalibrationFactor () const

Bool_t	IsGoodEvent ()

void	PrintValue () const

void	PrintInfo () const

Public Member Functions inherited from VQwDataElement
	VQwDataElement ()
	Default constructor. More...

	VQwDataElement (const VQwDataElement &value)
	Copy constructor. More...

virtual	~VQwDataElement ()
	Virtual destructor. More...

Bool_t	IsNameEmpty () const
	Is the name of this element empty? More...

void	SetElementName (const TString &name)
	Set the name of this element. More...

virtual const TString &	GetElementName () const
	Get the name of this element. More...

virtual void	LoadChannelParameters (QwParameterFile &paramfile)

size_t	GetNumberOfDataWords ()
	Get the number of data words in this data element. More...

UInt_t	GetGoodEventCount () const

virtual void	AssignValueFrom (const VQwDataElement *valueptr)

virtual VQwDataElement &	operator+= (const VQwDataElement &value)
	Addition-assignment operator. More...

virtual VQwDataElement &	operator-= (const VQwDataElement &value)
	Subtraction-assignment operator. More...

virtual void	Sum (const VQwDataElement &value1, const VQwDataElement &value2)
	Sum operator. More...

virtual void	Difference (const VQwDataElement &value1, const VQwDataElement &value2)
	Difference operator. More...

virtual void	Ratio (const VQwDataElement &numer, const VQwDataElement &denom)
	Ratio operator. More...

virtual void	SetSingleEventCuts (UInt_t errorflag, Double_t min, Double_t max, Double_t stability)
	set the upper and lower limits (fULimit and fLLimit), stability % and the error flag on this channel More...

virtual void	PrintErrorCounters () const
	report number of events failed due to HW and event cut failure More...

virtual UInt_t	GetEventcutErrorFlag ()
	return the error flag on this channel/device More...

virtual UInt_t	UpdateErrorFlag ()
	Update the error flag based on the error flags of internally contained objects Return paramter is the "Eventcut Error Flag". More...

virtual void	SetNeedsExternalClock (Bool_t needed)

virtual Bool_t	NeedsExternalClock ()

virtual std::string	GetExternalClockName ()

virtual void	SetExternalClockPtr (const VQwHardwareChannel *clock)

virtual void	SetExternalClockName (const std::string name)

virtual Double_t	GetNormClockValue ()

TString	GetSubsystemName () const
	Return the name of the inheriting subsystem name. More...

void	SetSubsystemName (TString sysname)
	Set the name of the inheriting subsystem name. More...

TString	GetModuleType () const
	Return the type of the beam instrument. More...

void	SetModuleType (TString ModuleType)
	set the type of the beam instrument More...

Public Member Functions inherited from MQwHistograms
void	ShareHistograms (const MQwHistograms *source)
	Share histogram pointers between objects. More...

Public Member Functions inherited from MQwMockable
	MQwMockable ()

virtual	~MQwMockable ()

void	SetRandomEventDriftParameters (Double_t amplitude, Double_t phase, Double_t frequency)
	Set a single set of harmonic drift parameters. More...

void	AddRandomEventDriftParameters (Double_t amplitude, Double_t phase, Double_t frequency)
	Add drift parameters to the internal set. More...

void	SetRandomEventParameters (Double_t mean, Double_t sigma)
	Set the normal random event parameters. More...

void	SetRandomEventAsymmetry (Double_t asymmetry)
	Set the helicity asymmetry. More...

Double_t	GetRandomValue ()

void	UseExternalRandomVariable ()
	Set the flag to use an externally provided random variable. More...

void	SetExternalRandomVariable (Double_t random_variable)
	Set the externally provided random variable. More...

Private Attributes
UInt_t	fChannel

Bool_t	fHasSamplingData

Bool_t	fHasAccumulatorData

Double_t	fPedestal

Double_t	fCalibrationFactor

Int_t	fCurrentEvent

UInt_t	fNumberOfEvents
	Current triggered event (allow for negative sentinel) More...

UInt_t	fSampleFormat
	Number of triggered events. More...

std::vector< QwSIS3320_Samples >	fSamples

std::vector< QwSIS3320_Samples >	fSamplesRaw

QwSIS3320_Samples	fAverageSamples

QwSIS3320_Samples	fAverageSamplesRaw

std::vector< Double_t >	fTimeWindowAverages

std::vector< std::pair< UInt_t, UInt_t > >	fTimeWindows

std::vector< Double_t >	fSampleWindowAverages

std::vector< std::pair < Double_t, Double_t > >	fSampleWindows

Int_t	fNumberOfAccumulators

Int_t	fAccumulatorDAC

Int_t	fAccumulatorThreshold1

Int_t	fAccumulatorThreshold2

Int_t	fAccumulatorTimingBefore5

Int_t	fAccumulatorTimingAfter5

Int_t	fAccumulatorTimingBefore6

Int_t	fAccumulatorTimingAfter6

std::vector < QwSIS3320_Accumulator >	fAccumulators

std::vector < QwSIS3320_Accumulator >	fAccumulatorsRaw

std::vector < QwSIS3320_LogicalAccumulator >	fLogicalAccumulators

Static Private Attributes
static const Bool_t	kDEBUG = kFALSE

static const Double_t	kVoltsPerBit = 5.0 / pow(2.0, 12)

static const Double_t	kNanoSecondsPerSample = 4.0

static const unsigned int	MODE_ACCUM_EVENT = 0x1

static const unsigned int	MODE_MULTI_EVENT = 0x3

static const unsigned int	MODE_SINGLE_EVENT = 0x4

static const unsigned int	MODE_NOTREADY = 0xda00

static const unsigned int	FORMAT_ACCUMULATOR = 0x0

static const unsigned int	FORMAT_LONG_WORD_SAMPLING = 0x1

static const unsigned int	FORMAT_SHORT_WORD_SAMPLING = 0x2

Additional Inherited Members
Public Types inherited from VQwDataElement
enum	EDataToSave { kRaw = 0, kDerived }

Protected Member Functions inherited from VQwDataElement
void	SetNumberOfDataWords (const UInt_t &numwords)
	Set the number of data words in this data element. More...

virtual VQwDataElement &	operator= (const VQwDataElement &value)
	Arithmetic assignment operator: Should only copy event-based data. More...

virtual void	UpdateErrorFlag (const UInt_t &error)

Protected Member Functions inherited from MQwHistograms
	MQwHistograms ()
	Default constructor. More...

	MQwHistograms (const MQwHistograms &source)
	Copy constructor. More...

virtual	~MQwHistograms ()
	Virtual destructor. More...

virtual MQwHistograms &	operator= (const MQwHistograms &value)

void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)

void	AddHistogram (TH1 *h)
	Register a histogram. More...

Protected Attributes inherited from VQwDataElement
TString	fElementName
	Name of this data element. More...

UInt_t	fNumberOfDataWords
	Number of raw data words in this data element. More...

Int_t	fGoodEventCount
	Number of good events accumulated in this element. More...

TString	fSubsystemName

TString	fModuleType

UInt_t	fErrorFlag
	This the standard error code generated for the channel that contains the global/local/stability flags and the Device error code (Unique error code for HW failures) More...

UInt_t	fErrorConfigFlag
	contains the global/local/stability flags More...

Protected Attributes inherited from MQwHistograms
std::vector< TH1_ptr >	fHistograms
	Histograms associated with this data element. More...

Protected Attributes inherited from MQwMockable
bool	fUseExternalRandomVariable
	Flag to use an externally provided normal random variable. More...

double	fExternalRandomVariable
	Externally provided normal random variable. More...

Double_t	fMockAsymmetry
	Helicity asymmetry. More...

Double_t	fMockGaussianMean
	Mean of normal distribution. More...

Double_t	fMockGaussianSigma
	Sigma of normal distribution. More...

std::vector< Double_t >	fMockDriftAmplitude
	Harmonic drift amplitude. More...

std::vector< Double_t >	fMockDriftFrequency
	Harmonic drift frequency. More...

std::vector< Double_t >	fMockDriftPhase
	Harmonic drift phase. More...

Static Protected Attributes inherited from MQwMockable
static boost::mt19937	fRandomnessGenerator
	Internal randomness generator. More...

static boost::normal_distribution < double >	fNormalDistribution
	Internal normal probability distribution. More...

static boost::variate_generator < boost::mt19937, boost::normal_distribution < double > >	fNormalRandomVariable
	Internal normal random variable. More...

Detailed Description

Class for the decoding of the SIS3320 sampling ADC data.

Author: W. Deconinck

Date: 2009-09-04 18:06:23

The QwSIS3320_Channel class is defined to read the integrated and sampled data from the Compton photon detector. Because the scope of this module is similar the the VQWK ADC module (integration and asymmetries), parts of this class are very similar to QwVQWK_Channel.

The main data members of a QwSIS3320_Channel are the vector fSamples of QwSIS3320_Samples where each entry stores a sample event, and the vector of fAccumulators of QwSIS3320_Accumulator where each entry stores an accumulator block. Both of these data members are derived from their –Raw counterparts by subtraction of pedestals and multiplication with calibration constants.

Definition at line 44 of file QwSIS3320_Channel.h.

Constructor & Destructor Documentation

QwSIS3320_Channel::QwSIS3320_Channel	(	UInt_t	channel = `0`,
		TString	name = `"auto"`
	)

inline

Definition at line 48 of file QwSIS3320_Channel.h.

References InitializeChannel().

     : VQwDataElement(),MQwMockable() {
       InitializeChannel(channel, name);
     };

Here is the call graph for this function:

QwSIS3320_Channel::QwSIS3320_Channel ( const QwSIS3320_Channel & source )

inline

Definition at line 52 of file QwSIS3320_Channel.h.

     : VQwDataElement(source),MQwMockable(source),
       fChannel(source.fChannel),
       fHasSamplingData(source.fHasSamplingData),
       fHasAccumulatorData(source.fHasAccumulatorData),
       fPedestal(source.fPedestal),
       fCalibrationFactor(source.fCalibrationFactor),
       fCurrentEvent(source.fCurrentEvent),
       fNumberOfEvents(source.fNumberOfEvents),
       fSampleFormat(source.fSampleFormat),
       fSamples(source.fSamples),
       fSamplesRaw(source.fSamplesRaw),
       fAverageSamples(source.fAverageSamples),
       fTimeWindowAverages(source.fTimeWindowAverages),
       fTimeWindows(source.fTimeWindows),
       fSampleWindowAverages(source.fTimeWindowAverages),
       fSampleWindows(source.fSampleWindows),
       fAccumulatorDAC(source.fAccumulatorDAC),
       fAccumulatorThreshold1(source.fAccumulatorThreshold1),
       fAccumulatorThreshold2(source.fAccumulatorThreshold2),
       fAccumulatorTimingBefore5(source.fAccumulatorTimingBefore5),
       fAccumulatorTimingAfter5(source.fAccumulatorTimingAfter5),
       fAccumulatorTimingBefore6(source.fAccumulatorTimingBefore6),
       fAccumulatorTimingAfter6(source.fAccumulatorTimingAfter6),
       fAccumulators(source.fAccumulators),
       fAccumulatorsRaw(source.fAccumulatorsRaw)
     { }

virtual QwSIS3320_Channel::~QwSIS3320_Channel ( )

inlinevirtual

Definition at line 79 of file QwSIS3320_Channel.h.

79 { };

Member Function Documentation

void QwSIS3320_Channel::AddLogicalAccumulator	(	const TString	name,
		const std::vector< TString >	accums,
		const std::vector< Double_t >	weights
	)

Create a logical accumulator

Definition at line 756 of file QwSIS3320_Channel.cc.

References QwSIS3320_LogicalAccumulator::AddAccumulatorReference(), QwLog::endl(), fAccumulators, fLogicalAccumulators, VQwDataElement::GetElementName(), and QwWarning.

 {
   // Create new logical accumulator with requested name
   QwSIS3320_LogicalAccumulator localLogAccum(GetElementName() + "_accum" + name);
 
   // Loop over all entries in this logical accumulator
   for (size_t i = 0; i < accums.size(); i++) {
 
     // Find the accumulator with the correct name
     size_t j;
     for (j = 0; j < fAccumulators.size(); j++) {
       if (accums[i] == fAccumulators[j].GetElementName())
         break;
     }
     // Not found
     if (j == fAccumulators.size()) {
       QwWarning << "Logical accumulator " << name << " has undefined components: " << QwLog::endl;
       QwWarning << accums[i] << " could not be found." << QwLog::endl;
       return;
     }
 
     // Add reference and weight
     localLogAccum.AddAccumulatorReference(&(fAccumulators[j]),weights[i]);
   }
 
   // Add logical accumulator to the list
   fLogicalAccumulators.push_back(localLogAccum);
 }

Here is the call graph for this function:

void QwSIS3320_Channel::ClearEventData ( )

virtual

Clear the event data in sampling buffer and accumulators

Reimplemented from VQwDataElement.

Definition at line 118 of file QwSIS3320_Channel.cc.

References QwSIS3320_Samples::ClearEventData(), fAccumulators, fAccumulatorsRaw, fAverageSamples, fAverageSamplesRaw, fLogicalAccumulators, fSamples, and fSamplesRaw.

 {
   // Clear the sampled events
   for (size_t i = 0; i < fSamples.size(); i++)
     fSamples.at(i).ClearEventData();
   fSamples.clear(); // and back to zero events
   for (size_t i = 0; i < fSamplesRaw.size(); i++)
     fSamplesRaw.at(i).ClearEventData();
   fSamplesRaw.clear(); // and back to zero events
   // Clear the average samples
   fAverageSamples.ClearEventData();
   fAverageSamplesRaw.ClearEventData();
 
   // Clear the accumulators
   for (size_t i = 0; i < fAccumulators.size(); i++)
     fAccumulators.at(i).ClearEventData();
   for (size_t i = 0; i < fAccumulatorsRaw.size(); i++)
     fAccumulatorsRaw.at(i).ClearEventData();
   for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
     fLogicalAccumulators.at(i).ClearEventData();
   // (the number of accumulators is constant, don't clear them)
 }

Here is the call graph for this function:

void QwSIS3320_Channel::ConstructBranchAndVector	(	TTree *	tree,
		TString &	prefix,
		std::vector< Double_t > &	values
	)

Definition at line 673 of file QwSIS3320_Channel.cc.

References fAccumulators, fAccumulatorsRaw, fLogicalAccumulators, fSamples, and VQwDataElement::GetElementName().

 {
   // Accumulators
   for (size_t i = 0; i < fAccumulators.size(); i++) {
     fAccumulators[i].ConstructBranchAndVector(tree, prefix, values);
     fAccumulatorsRaw[i].ConstructBranchAndVector(tree, prefix, values);
   }
   // Logical Accumulators
   for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
     fLogicalAccumulators[i].ConstructBranchAndVector(tree, prefix, values);
 
   // Samples (only collected when running over data, so structure does not
   // actually exist yet at time of branch construction)
   TString basename = prefix + GetElementName() + "_samples";
   tree->Branch(basename, &fSamples);
   //tree->Branch(basename + "_avg", &fAverageSamples);
 }

Here is the call graph for this function:

void QwSIS3320_Channel::ConstructHistograms	(	TDirectory *	folder,
		TString &	prefix
	)

virtual

Construct the histograms for this data element.

Implements VQwDataElement.

Definition at line 615 of file QwSIS3320_Channel.cc.

References QwHistogramHelper::Construct1DHist(), fAccumulators, fAccumulatorsRaw, MQwHistograms::fHistograms, fLogicalAccumulators, VQwDataElement::GetElementName(), gQwHists, and VQwDataElement::IsNameEmpty().

 {
   //  If we have defined a subdirectory in the ROOT file, then change into it.
   if (folder != NULL) folder->cd();
 
   if (IsNameEmpty()) {
     //  This channel is not used, so skip filling the histograms.
   } else {
 
     // Accumulators
     for (size_t i = 0; i < fAccumulators.size(); i++) {
       fAccumulators[i].ConstructHistograms(folder,prefix);
       fAccumulatorsRaw[i].ConstructHistograms(folder,prefix);
     }
 
     // Logical Accumulators
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i].ConstructHistograms(folder,prefix);
 
     TString basename = prefix + GetElementName();
 
     fHistograms.resize(3, NULL);
     size_t index = 0;
     fHistograms[index++] = gQwHists.Construct1DHist(basename+Form("_sum"));
     fHistograms[index++] = gQwHists.Construct1DHist(basename+Form("_ped"));
     fHistograms[index++] = gQwHists.Construct1DHist(basename+Form("_ped_raw"));
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::Difference	(	QwSIS3320_Channel &	value1,
		QwSIS3320_Channel &	value2
	)

Difference of two channels

Parameters

value1
value2

Definition at line 568 of file QwSIS3320_Channel.cc.

 {
   *this  = value1;
   *this -= value2;
 }

void QwSIS3320_Channel::EncodeEventData ( std::vector< UInt_t > & buffer )

virtual

Encode the event data into a CODA buffer.

Implements MQwMockable.

Definition at line 336 of file QwSIS3320_Channel.cc.

References fChannel, fNumberOfEvents, fSampleFormat, and VQwDataElement::IsNameEmpty().

 {
   //  Long_t sample;
   std::vector<UInt_t> header;
   std::vector<UInt_t> samples;
 
   if (IsNameEmpty()) {
     //  This channel is not used, but is present in the data stream.
     //  Skip over this data.
   } else {
     header.push_back(fChannel); // Channel number
     header.push_back(fSampleFormat); // Sample format (e.g. 0x20003)
 // TODO
 //    header.push_back(fNumberOfSamples * fNumberOfEvents); // Total number of samples
 //    header.push_back(fNumberOfSamples); // Number of samples per event
     header.push_back(fNumberOfEvents); // Number of events
     // Data is stored with two short words per long word
 //     for (size_t i = 0; i < fSamplesRaw.size() / 2; i++) {
 //       sample  = fSamplesRaw.at(2*i) << 16;
 //       sample |= fSamplesRaw.at(2*i+1);
 //       samples.push_back(sample);
 //     }
 
     for (size_t i = 0; i < header.size(); i++)
       buffer.push_back(header.at(i));
     for (size_t i = 0; i < samples.size(); i++)
       buffer.push_back(samples.at(i));
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::FillHistograms ( )

virtual

Fill the histograms for this data element.

Implements VQwDataElement.

Definition at line 644 of file QwSIS3320_Channel.cc.

References fAccumulators, fAccumulatorsRaw, fAverageSamples, fAverageSamplesRaw, MQwHistograms::fHistograms, fLogicalAccumulators, fSamples, QwSIS3320_Samples::GetSample(), QwSIS3320_Samples::GetSum(), and VQwDataElement::IsNameEmpty().

 {
   size_t index = -1;
   if (fSamples.size() == 0) return;
 
   if (IsNameEmpty()) {
       //  This channel is not used, so skip creating the histograms.
   } else {
 
     // Accumulators
     for (size_t i = 0; i < fAccumulators.size(); i++) {
       fAccumulators[i].FillHistograms();
       fAccumulatorsRaw[i].FillHistograms();
     }
 
     // Logical Accumulators
     for (size_t i = 0; i< fLogicalAccumulators.size(); i++ ) {
       fLogicalAccumulators[i].FillHistograms();
     }
 
     if (fHistograms[++index] != NULL)
       fHistograms[index]->Fill(fAverageSamples.GetSum());
     if (fHistograms[++index] != NULL)
       fHistograms[index]->Fill(fAverageSamples.GetSample(0));
     if (fHistograms[++index] != NULL)
       fHistograms[index]->Fill(fAverageSamplesRaw.GetSample(0));
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::FillTreeVector ( std::vector< Double_t > & values ) const

Definition at line 691 of file QwSIS3320_Channel.cc.

References fAccumulators, fAccumulatorsRaw, and fLogicalAccumulators.

 {
   // Accumulators
   for (size_t i = 0; i < fAccumulators.size(); i++) {
     fAccumulators[i].FillTreeVector(values);
     fAccumulatorsRaw[i].FillTreeVector(values);
   }
 
   // Logical Accumulators
   for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
     fLogicalAccumulators[i].FillTreeVector(values);
 }

Double_t QwSIS3320_Channel::GetCalibrationFactor ( ) const

inline

Definition at line 136 of file QwSIS3320_Channel.h.

References fCalibrationFactor.

136 { return fCalibrationFactor; };

QwSIS3320_Channel::fCalibrationFactor

Double_t fCalibrationFactor

Definition: QwSIS3320_Channel.h:156

size_t QwSIS3320_Channel::GetNumberOfAccumulators ( ) const

inline

Definition at line 126 of file QwSIS3320_Channel.h.

References fNumberOfAccumulators.

126 { return (fNumberOfAccumulators); };

QwSIS3320_Channel::fNumberOfAccumulators

Int_t fNumberOfAccumulators

Definition: QwSIS3320_Channel.h:175

size_t QwSIS3320_Channel::GetNumberOfEvents ( ) const

inline

Definition at line 119 of file QwSIS3320_Channel.h.

References fNumberOfEvents.

Referenced by PrintInfo(), PrintValue(), and ProcessEvBuffer().

119 { return (fNumberOfEvents); };

QwSIS3320_Channel::fNumberOfEvents

UInt_t fNumberOfEvents

Current triggered event (allow for negative sentinel)

Definition: QwSIS3320_Channel.h:160

Here is the caller graph for this function:

Double_t QwSIS3320_Channel::GetPedestal ( ) const

inline

Definition at line 134 of file QwSIS3320_Channel.h.

References fPedestal.

134 { return fPedestal; };

QwSIS3320_Channel::fPedestal

Double_t fPedestal

Definition: QwSIS3320_Channel.h:155

QwSIS3320_Samples& QwSIS3320_Channel::GetSamples ( size_t i )

inline

Definition at line 116 of file QwSIS3320_Channel.h.

References fSamples.

116 { return fSamples.at(i); };

QwSIS3320_Channel::fSamples

std::vector< QwSIS3320_Samples > fSamples

Definition: QwSIS3320_Channel.h:164

QwSIS3320_Samples& QwSIS3320_Channel::GetSamplesRaw ( size_t i )

inline

Definition at line 117 of file QwSIS3320_Channel.h.

References fSamplesRaw.

117 { return fSamplesRaw.at(i); };

QwSIS3320_Channel::fSamplesRaw

std::vector< QwSIS3320_Samples > fSamplesRaw

Definition: QwSIS3320_Channel.h:165

void QwSIS3320_Channel::InitializeChannel	(	UInt_t	channel,
		TString	name
	)

Initialize the QwSIS3320_Channel by assigning it a name

Parameters

channel	Number of the channel
name	Name for the channel

Definition at line 59 of file QwSIS3320_Channel.cc.

References QwSIS3320_Samples::ClearEventData(), fAccumulators, fAccumulatorsRaw, fAverageSamples, fAverageSamplesRaw, fCalibrationFactor, fChannel, fCurrentEvent, fHasAccumulatorData, fHasSamplingData, MQwMockable::fMockAsymmetry, MQwMockable::fMockGaussianMean, MQwMockable::fMockGaussianSigma, fPedestal, fSamples, fSamplesRaw, VQwDataElement::GetElementName(), VQwDataElement::SetElementName(), and VQwDataElement::SetNumberOfDataWords().

Referenced by QwSIS3320_Channel().

 {
   // Inherited from VQwDataElement
   SetElementName(name);
   SetNumberOfDataWords(0);
 
   // Set channel number
   fChannel = channel;
 
   // Set accumulator names
   for (size_t i = 0; i < fAccumulators.size(); i++) {
     TString name = GetElementName() + TString("_accum") + Form("%ld",i);
     fAccumulators[i].SetElementName(name);
     fAccumulatorsRaw[i].SetElementName(name + "_raw");
   }
 
   // Start with zero samples
   fSamples.clear(); fSamplesRaw.clear();
   // Clear the average samples
   fAverageSamples.ClearEventData();
   fAverageSamplesRaw.ClearEventData();
 
   // Enabled by MMD
   //for (size_t i = 0; i < fSamples.size(); i++) {
   //  TString name = GetElementName() + TString("_samples") + Form("%ld",i);
   //  fSamples[i].SetElementName(name);
   //  fSamplesRaw[i].SetElementName(name + "_raw");
   //}
 
   // Default values when no event read yet
   fCurrentEvent = -1;
   fHasSamplingData = kFALSE;
   fHasAccumulatorData = kFALSE;
 
   // Pedestal calibration
   fPedestal = 0.0;
   fCalibrationFactor = 1.0;
 
   // Mock data parameters
   fMockAsymmetry = 0.0;
   fMockGaussianMean = 0.0;
   fMockGaussianSigma = 0.0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Bool_t QwSIS3320_Channel::IsGoodEvent ( )

Check whether the event is a good event from the number of read samples

Returns: kTRUE if samples have been read

Definition at line 107 of file QwSIS3320_Channel.cc.

References fSamples.

 {
   Bool_t fEventIsGood = kTRUE;
   for (size_t i = 0; i < fSamples.size(); i++)
     fEventIsGood &= (fSamples[i].GetNumberOfSamples() > 0);
   return fEventIsGood;
 }

void QwSIS3320_Channel::Offset ( Double_t offset )

Addition of a offset

Parameters

offset

Definition at line 589 of file QwSIS3320_Channel.cc.

References fAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] += offset;
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i] += offset;
   }
 }

Here is the call graph for this function:

const QwSIS3320_Channel QwSIS3320_Channel::operator+ ( const Double_t & value ) const

Addition of offset

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 418 of file QwSIS3320_Channel.cc.

 {
   QwSIS3320_Channel result = *this;
   result += value;
   return result;
 }

const QwSIS3320_Channel QwSIS3320_Channel::operator+ ( const QwSIS3320_Channel & value ) const

Addition

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 442 of file QwSIS3320_Channel.cc.

 {
   QwSIS3320_Channel result = *this;
   result += value;
   return result;
 }

QwSIS3320_Channel & QwSIS3320_Channel::operator+= ( const Double_t & value )

Addition assignment of offset

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 485 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] += value;
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i] += value;
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i] += value;
   }
   return *this;
 }

Here is the call graph for this function:

QwSIS3320_Channel & QwSIS3320_Channel::operator+= ( const QwSIS3320_Channel & value )

Addition assignment

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 521 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] += value.fSamples.at(i);
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i] += value.fAccumulators.at(i);
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i] += value.fLogicalAccumulators.at(i);
   }
   return *this;
 }

Here is the call graph for this function:

const QwSIS3320_Channel QwSIS3320_Channel::operator- ( const Double_t & value ) const

Subtraction of offset

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 430 of file QwSIS3320_Channel.cc.

 {
   QwSIS3320_Channel result = *this;
   result -= value;
   return result;
 }

const QwSIS3320_Channel QwSIS3320_Channel::operator- ( const QwSIS3320_Channel & value ) const

Subtraction

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 454 of file QwSIS3320_Channel.cc.

 {
   QwSIS3320_Channel result = *this;
   result -= value;
   return result;
 }

QwSIS3320_Channel & QwSIS3320_Channel::operator-= ( const Double_t & value )

Subtraction assignment of offset

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 503 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] -= value;
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i] -= value;
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i] -= value;
   }
   return *this;
 }

Here is the call graph for this function:

QwSIS3320_Channel & QwSIS3320_Channel::operator-= ( const QwSIS3320_Channel & value )

Subtraction assignment

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 539 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] -= value.fSamples.at(i);
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i] -= value.fAccumulators.at(i);
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i] -= value.fLogicalAccumulators.at(i);
   }
   return *this;
 }

Here is the call graph for this function:

QwSIS3320_Channel & QwSIS3320_Channel::operator= ( const QwSIS3320_Channel & value )

Assignment

Parameters

value Right-hand side

Returns: Left-hand side

Definition at line 466 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] = value.fSamples.at(i);
     for (size_t i = 0; i < fAccumulators.size(); i++ ) {
       fAccumulators[i] = value.fAccumulators.at(i);
     }
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++ )
       fLogicalAccumulators[i] = value.fLogicalAccumulators.at(i);
   }
   return *this;
 }

Here is the call graph for this function:

void QwSIS3320_Channel::PrintInfo ( ) const

virtual

Print some debugging information about the QwSIS3320_Channel

Reimplemented from VQwDataElement.

Definition at line 723 of file QwSIS3320_Channel.cc.

References QwLog::endl(), fAccumulators, fAccumulatorsRaw, fAverageSamples, fAverageSamplesRaw, fLogicalAccumulators, fSamples, fSamplesRaw, VQwDataElement::GetElementName(), GetNumberOfEvents(), QwSIS3320_Samples::GetSum(), QwSIS3320_Samples::GetSumInTimeWindow(), and QwOut.

 {
   QwOut << "QwSIS3320_Channel \"" << GetElementName() << "\":" << QwLog::endl;
   QwOut << "Number of events: " << GetNumberOfEvents() << QwLog::endl;
   QwOut << "fSamplesRaw:" << QwLog::endl;
   for (size_t i = 0; i < fSamplesRaw.size(); i++) {
     QwOut << "event "  << i << ": " << fSamplesRaw.at(i) << QwLog::endl;
   }
   QwOut << "avg: " << fAverageSamplesRaw << QwLog::endl;
   QwOut << "avg sum: " << fAverageSamplesRaw.GetSum() << QwLog::endl;
   QwOut << "fSamples:" << QwLog::endl;
   for (size_t i = 0; i < fSamples.size(); i++) {
     QwOut << "event " << i << ": " << fSamples.at(i) << QwLog::endl;
   }
   QwOut << "avg: " << fAverageSamples << QwLog::endl;
   QwOut << "avg sum: " << fAverageSamples.GetSum() << QwLog::endl;
   QwOut << "fAccumulators -> fAccumulatorsRaw:" << QwLog::endl;
   for (size_t i = 0; i < fAccumulators.size(); i++) {
     QwOut << i << ": " << fAccumulatorsRaw.at(i) << " -> " << fAccumulators.at(i) << QwLog::endl;
   }
   QwOut << "fLogicalAccumulators:" << QwLog::endl;
   for (size_t i = 0; i < fLogicalAccumulators.size(); i++) {
     QwOut << i << ": " << fLogicalAccumulators.at(i) << QwLog::endl;
   }
   UInt_t nped = 10;
   if (fSamplesRaw.size() > 0) {
     QwOut << "Average pedestal is: " << fAverageSamplesRaw.GetSumInTimeWindow(0, nped) / (Double_t) nped << QwLog::endl;
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::PrintValue ( ) const

virtual

Print value of the QwSIS3320_Channel

Reimplemented from VQwDataElement.

Definition at line 707 of file QwSIS3320_Channel.cc.

References QwLog::endl(), fAccumulators, fLogicalAccumulators, VQwDataElement::GetElementName(), GetNumberOfEvents(), and QwMessage.

 {
   QwMessage << std::setprecision(4)
             << std::setw(18) << std::left << GetElementName() << ", "
             << std::setw(15) << std::left << GetNumberOfEvents();
   for (size_t i = 0; i < fAccumulators.size(); i++)
     QwMessage << ", " << i << ":" << fAccumulators[i].GetAccumulatorSum();
   for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
     QwMessage << ", " << i << ":" << fLogicalAccumulators[i].GetAccumulatorSum();
   QwMessage << QwLog::endl;
 }

Here is the call graph for this function:

Int_t QwSIS3320_Channel::ProcessEvBuffer	(	UInt_t *	buffer,
		UInt_t	num_words_left,
		UInt_t	index = `0`
	)

virtual

Extract the sampling and accumulator data from the CODA buffer

Parameters

buffer	Input buffer
num_words_left	Number of words left in the input buffer
index	Starting position in the buffer

Returns: Number of words processed in this method

Implements VQwDataElement.

Definition at line 148 of file QwSIS3320_Channel.cc.

References QwLog::endl(), fAccumulatorDAC, fAccumulatorsRaw, fAccumulatorThreshold1, fAccumulatorThreshold2, fAccumulatorTimingAfter5, fAccumulatorTimingAfter6, fAccumulatorTimingBefore5, fAccumulatorTimingBefore6, fChannel, fHasAccumulatorData, fHasSamplingData, VQwDataElement::fNumberOfDataWords, FORMAT_ACCUMULATOR, FORMAT_LONG_WORD_SAMPLING, FORMAT_SHORT_WORD_SAMPLING, fSamplesRaw, GetNumberOfEvents(), VQwDataElement::IsNameEmpty(), kDEBUG, MODE_ACCUM_EVENT, MODE_MULTI_EVENT, MODE_NOTREADY, MODE_SINGLE_EVENT, QwError, QwOut, QwWarning, and SetNumberOfEvents().

 {
   UInt_t words_read = 0;
 
   if (IsNameEmpty()) {
     // This channel is not used, but is present in the data stream.
     // Skip over this data.
     words_read = fNumberOfDataWords;
   } else if (num_words_left >= fNumberOfDataWords) {
 
     // The first word of the buffer is in the format: 0xFADC/id/ch/, where id
     // is the module identification, and ch the channel number in that module.
 
     // Check whether the first word has 0xFADC as the upper half
     UInt_t local_tag = (buffer[0] >> 16) & 0xFFFF;
     if (local_tag != 0xFADC) return words_read;
 
     // Check whether the first word contains the correct channel number
     //UInt_t local_module = (buffer[0] >> 8) & 0xFF;
     UInt_t local_channel = buffer[0] & 0xFF;
  
     if (local_channel != fChannel) return words_read;
 
     // Determine whether we are receiving a sampling or accumulator buffer:
     // Most significant short word (upper half):
     // - 0x0 in accumulator mode
     // - 0x1 when an expanded sampling buffer follows (12-bit sample / 32-bit word)
     // - 0x2 when a packed sampling buffer follows (2 * 12-bit samples / 32-bit word)
     // Least significant short word (lower half): 0x/stop mode/setup mode/ where
     // - stop mode is 1 for hardware trigger
     // - setup mode is 1,2,3,4...
     UInt_t local_format = (buffer[1] >> 16) & 0xFFFF;
     //UInt_t local_stopmode = (buffer[1] >> 8) & 0xFF;
 
     UInt_t local_setupmode = (buffer[1]) & 0xFF;
     words_read = 2;
     UInt_t packedtiming; // locals declared outside of switch/case
     switch (local_format) {
 
       // This is a sampling buffer using long words
       case FORMAT_LONG_WORD_SAMPLING:
         QwWarning << "QwSIS3320_Channel: Expanded word format not implemented"
                   << QwLog::endl;
         break; // end of FORMAT_LONG_WORD_SAMPLING
 
       // This is a sampling buffer using short words
       case FORMAT_SHORT_WORD_SAMPLING:
         UInt_t numberofsamples, numberofevents_expected, numberofevents_actual;
         UInt_t samplepointer;
         switch (local_setupmode) {
 
           // This is a sampling buffer in multi event mode:
           // - many events are saved in one buffer for a complete helicity event
           case MODE_ACCUM_EVENT:
           case MODE_MULTI_EVENT:
             samplepointer = buffer[2];
             numberofsamples = buffer[3];
             numberofevents_expected = buffer[4];
             words_read = 5;
 
             // For all events in this buffer
             SetNumberOfEvents(numberofevents_expected);
             numberofevents_actual = 0; // double check while reading the events
             for (size_t event = 0; event < GetNumberOfEvents(); event++) {
               // create a new raw sampled event
               fSamplesRaw[event].SetNumberOfSamples(numberofsamples);
               fSamplesRaw[event].SetSamplePointer(samplepointer);
               // pass the buffer to read the samples
               UInt_t samples_read = fSamplesRaw[event].ProcessEvBuffer(&(buffer[words_read]), num_words_left-words_read);
               // check whether we actually read any data
               if (samples_read == 0) break;
               // an actual sampled event was read
               words_read += samples_read;
               numberofevents_actual++;
               if (kDEBUG) QwOut << "Samples " << event << ": " << fSamplesRaw[event] << QwLog::endl;
             }
             if (numberofevents_expected != numberofevents_actual) {
               QwWarning << "QwSIS3320_Channel: Expected " << numberofevents_expected << " events, "
                         << "but only read " << numberofevents_actual << "." << QwLog::endl;
               SetNumberOfEvents(numberofevents_actual);
             }
 
             break;
 
           // This is a sampling buffer in single event mode:
           // - one event only is saved in the buffer, and each event triggers
           //   a read-out so multiple event can occur in one helicity period
           // - because a circular buffer is used the pointer to the last sample
           //   is stored
           case MODE_SINGLE_EVENT:
             samplepointer = buffer[2];
             numberofsamples = buffer[3];
             numberofevents_expected = 1;
 
             // Create a new raw sampled event
             SetNumberOfEvents(numberofevents_expected);
             // Pass the buffer to read the samples (only one event)
             fSamplesRaw.at(0).SetNumberOfSamples(numberofsamples);
             fSamplesRaw.at(0).SetSamplePointer(samplepointer);
             fSamplesRaw.at(0).ProcessEvBuffer(buffer, num_words_left, samplepointer);
 
             break;
 
           // Default
           default:
             QwError << "QwSIS3320_Channel: Received unknown sampling format: "
                     << std::hex << local_setupmode << std::dec << "." << QwLog::endl;
             words_read = 0;
             return words_read;
 
         } // end of switch (local_setupmode)
 
         fHasSamplingData = kTRUE;
         break; // end of FORMAT_SHORT_WORD_SAMPLING
 
       // This is an accumulator buffer
       case FORMAT_ACCUMULATOR:
 
         // Definition of the accumulators: (numbered from 1)
         // 1: digital sum of all samples in the event
         //
         // 2: threshold 1 < sample value
         // 3: threshold 2 < sample value <= threshold 1
         // 4: sample value <= threshold 2
         // Test: 1 == 2 + 3 + 4
         //
         // 6: sample value <= threshold 2
         //    and add samples before and after
         // 5: threshold 2 < sample value <= threshold 1
         //    and add samples before and after
         //    but samples added in accumulator 6 are not added to accumulator 5
 
         // Read the accumulator blocks
         for (size_t i = 0; i < fAccumulatorsRaw.size(); i++) {
           words_read += fAccumulatorsRaw[i].ProcessEvBuffer(&(buffer[words_read]), num_words_left-words_read);
           if (kDEBUG) QwOut << "Accum " << i+1 << ": " << fAccumulatorsRaw[i] << QwLog::endl;
         }
         // Read the threshold information
         fAccumulatorDAC = buffer[words_read++];
         // Thresholds are maximum 12 bits, so they seem to have wasted a word here :-)
         fAccumulatorThreshold1 = buffer[words_read++];
         fAccumulatorThreshold2 = buffer[words_read++];
         // Read the timing information, packed in one word
         packedtiming = buffer[words_read++];
         fAccumulatorTimingAfter6  = (packedtiming >>= 8) & 0xFF;
         fAccumulatorTimingBefore6 = (packedtiming >>= 8) & 0xFF;
         fAccumulatorTimingAfter5  = (packedtiming >>= 8) & 0xFF;
         fAccumulatorTimingBefore5 = (packedtiming >>= 8) & 0xFF;
         if (kDEBUG) {
           QwOut << "DAC: " << fAccumulatorDAC << QwLog::endl;
           QwOut << "Thresholds: " << fAccumulatorThreshold1 << ", "
                                   << fAccumulatorThreshold2 << QwLog::endl;
           QwOut << "Timings on accum 5: " << fAccumulatorTimingAfter5 << ", "
                                           << fAccumulatorTimingBefore5 << QwLog::endl;
           QwOut << "Timings on accum 6: " << fAccumulatorTimingAfter6 << ", "
                                           << fAccumulatorTimingBefore6 << QwLog::endl;
         }
         if (kDEBUG) {
           QwSIS3320_Accumulator sum("sum");
           sum = fAccumulatorsRaw[1] + fAccumulatorsRaw[2] + fAccumulatorsRaw[3];
           QwOut << "Accum 1: " << fAccumulatorsRaw[0] << QwLog::endl;
           QwOut << "Accum 2 + 3 + 4: " << sum << QwLog::endl;
         }
 
         fHasAccumulatorData = kTRUE;
         break; // end of FORMAT_ACCUMULATOR
 
       // This is an incomplete buffer
       case MODE_NOTREADY:
         QwError << "QwSIS3320_Channel: SIS3320 was not ready yet!" << QwLog::endl;
         break; // end of FORMAT_NOTREADY
 
       // Default
       default:
         QwError << "QwSIS3320_Channel: Received unknown mode: "
                 << std::hex << local_format << std::dec << "." << QwLog::endl;
         words_read = 0;
         return words_read;
 
     } // end of switch (local_format)
 
   } else {
     QwError << "QwSIS3320_Channel: Not enough words while processing buffer!" << QwLog::endl;
   }
 
   return words_read;
 }

Here is the call graph for this function:

void QwSIS3320_Channel::ProcessEvent ( )

Process the event by removing pedestals and applying calibration

Definition at line 369 of file QwSIS3320_Channel.cc.

References fAccumulators, fAccumulatorsRaw, fAverageSamples, fAverageSamplesRaw, fCalibrationFactor, fLogicalAccumulators, fPedestal, fSamples, fSamplesRaw, fTimeWindowAverages, fTimeWindows, and QwSIS3320_Samples::SetNumberOfSamples().

 {
   // Correct for pedestal and calibration factor
   if (fSamplesRaw.size() > 0) fSamples.resize(fSamplesRaw.size(), fSamplesRaw[0]);
   for (size_t i = 0; i < fSamplesRaw.size(); i++) {
     fSamples[i] = fSamplesRaw[i];
     fSamples[i] -= fPedestal;
     fSamples[i] *= fCalibrationFactor;
     fSamples[i].UpdateGraph();
   }
   for (size_t i = 0; i < fAccumulatorsRaw.size(); i++) {
     fAccumulators[i] = fAccumulatorsRaw[i];
     fAccumulators[i] -= fPedestal * fAccumulatorsRaw[i].GetNumberOfSamples();
     fAccumulators[i] *= fCalibrationFactor;
   }
 
   for (size_t i = 0; i < fLogicalAccumulators.size(); i++) {
     fLogicalAccumulators[i].ProcessEvent();
   }
 
   // Calculate the average sample snapshot
   if (fSamples.size() > 0) {
     fAverageSamples.SetNumberOfSamples(fSamples[0].GetNumberOfSamples());
     fAverageSamplesRaw.SetNumberOfSamples(fSamples[0].GetNumberOfSamples());
   }
   for (size_t i = 0; i < fSamples.size(); i++) {
     fAverageSamples += fSamples[i];
     fAverageSamplesRaw += fSamplesRaw[i];
   }
   if (fSamples.size() > 0) {
     fAverageSamples /= fSamples.size();
     fAverageSamplesRaw /= fSamplesRaw.size();
   }
 
   // Calculate the windowed sample sums
   for (size_t i = 0; i < fSamples.size(); i++) {
     for (size_t timewindow = 0; timewindow < fTimeWindows.size(); timewindow++) {
       UInt_t start = fTimeWindows[timewindow].first;
       UInt_t stop = fTimeWindows[timewindow].second;
       fTimeWindowAverages[timewindow] += fSamples[i].GetSumInTimeWindow(start, stop);
     }
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::RandomizeEventData	(	int	helicity = `0`,
		double	time = `0.0`
	)

inlinevirtual

Internally generate random event data.

Implements MQwMockable.

Definition at line 89 of file QwSIS3320_Channel.h.

89 { };

void QwSIS3320_Channel::Ratio	(	QwSIS3320_Channel &	numer,
		QwSIS3320_Channel &	denom
	)

Definition at line 574 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i].Ratio(numer.fAccumulators[i],denom.fAccumulators[i]);
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i].Ratio(numer.fLogicalAccumulators[i],
           denom.fLogicalAccumulators[i]);
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::Scale ( Double_t scale )

Scaling by a scale factor

Parameters

scale

Definition at line 603 of file QwSIS3320_Channel.cc.

References fAccumulators, fLogicalAccumulators, fSamples, and VQwDataElement::IsNameEmpty().

 {
   if (!IsNameEmpty()) {
     for (size_t i = 0; i < fSamples.size(); i++)
       fSamples[i] *= scale;
     for (size_t i = 0; i < fAccumulators.size(); i++)
       fAccumulators[i] *= scale;
     for (size_t i = 0; i < fLogicalAccumulators.size(); i++)
       fLogicalAccumulators[i] *= scale;
   }
 }

Here is the call graph for this function:

void QwSIS3320_Channel::SetCalibrationFactor ( const Double_t factor )

inline

Definition at line 135 of file QwSIS3320_Channel.h.

References fCalibrationFactor.

135 { fCalibrationFactor = factor; };

QwSIS3320_Channel::fCalibrationFactor

Double_t fCalibrationFactor

Definition: QwSIS3320_Channel.h:156

void QwSIS3320_Channel::SetNumberOfAccumulators ( UInt_t naccumulators )

inline

Definition at line 127 of file QwSIS3320_Channel.h.

References fAccumulators, fAccumulatorsRaw, and fNumberOfAccumulators.

                                                        {
       fNumberOfAccumulators = naccumulators;
       fAccumulatorsRaw.resize(naccumulators);
       fAccumulators.resize(naccumulators);
     };

void QwSIS3320_Channel::SetNumberOfEvents ( UInt_t nevents )

inline

Definition at line 120 of file QwSIS3320_Channel.h.

References fNumberOfEvents, fSamples, and fSamplesRaw.

Referenced by ProcessEvBuffer().

                                            {
       fNumberOfEvents = nevents;
       fSamplesRaw.resize(nevents);
       fSamples.resize(nevents);
     };

Here is the caller graph for this function:

void QwSIS3320_Channel::SetPedestal ( const Double_t ped )

inline

Definition at line 133 of file QwSIS3320_Channel.h.

References fPedestal.

133 { fPedestal = ped; };

QwSIS3320_Channel::fPedestal

Double_t fPedestal

Definition: QwSIS3320_Channel.h:155

void QwSIS3320_Channel::Sum	(	QwSIS3320_Channel &	value1,
		QwSIS3320_Channel &	value2
	)

Sum of two channels

Parameters

value1
value2

Definition at line 557 of file QwSIS3320_Channel.cc.

 {
   *this  = value1;
   *this += value2;
 }

Field Documentation

Int_t QwSIS3320_Channel::fAccumulatorDAC

private

Definition at line 176 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

std::vector<QwSIS3320_Accumulator> QwSIS3320_Channel::fAccumulators

private

Definition at line 180 of file QwSIS3320_Channel.h.

Referenced by AddLogicalAccumulator(), ClearEventData(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), InitializeChannel(), Offset(), operator+=(), operator-=(), operator=(), PrintInfo(), PrintValue(), ProcessEvent(), Ratio(), Scale(), and SetNumberOfAccumulators().

std::vector<QwSIS3320_Accumulator> QwSIS3320_Channel::fAccumulatorsRaw

private

Definition at line 181 of file QwSIS3320_Channel.h.

Referenced by ClearEventData(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), InitializeChannel(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), and SetNumberOfAccumulators().

Int_t QwSIS3320_Channel::fAccumulatorThreshold1

private

Definition at line 177 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

Int_t QwSIS3320_Channel::fAccumulatorThreshold2

private

Definition at line 177 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

Int_t QwSIS3320_Channel::fAccumulatorTimingAfter5

private

Definition at line 178 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

Int_t QwSIS3320_Channel::fAccumulatorTimingAfter6

private

Definition at line 179 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

Int_t QwSIS3320_Channel::fAccumulatorTimingBefore5

private

Definition at line 178 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

Int_t QwSIS3320_Channel::fAccumulatorTimingBefore6

private

Definition at line 179 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

QwSIS3320_Samples QwSIS3320_Channel::fAverageSamples

private

Definition at line 167 of file QwSIS3320_Channel.h.

Referenced by ClearEventData(), FillHistograms(), InitializeChannel(), PrintInfo(), and ProcessEvent().

QwSIS3320_Samples QwSIS3320_Channel::fAverageSamplesRaw

private

Definition at line 168 of file QwSIS3320_Channel.h.

Referenced by ClearEventData(), FillHistograms(), InitializeChannel(), PrintInfo(), and ProcessEvent().

Double_t QwSIS3320_Channel::fCalibrationFactor

private

Definition at line 156 of file QwSIS3320_Channel.h.

Referenced by GetCalibrationFactor(), InitializeChannel(), ProcessEvent(), and SetCalibrationFactor().

UInt_t QwSIS3320_Channel::fChannel

private

Definition at line 148 of file QwSIS3320_Channel.h.

Referenced by EncodeEventData(), InitializeChannel(), and ProcessEvBuffer().

Int_t QwSIS3320_Channel::fCurrentEvent

private

Definition at line 159 of file QwSIS3320_Channel.h.

Referenced by InitializeChannel().

Bool_t QwSIS3320_Channel::fHasAccumulatorData

private

Definition at line 150 of file QwSIS3320_Channel.h.

Referenced by InitializeChannel(), and ProcessEvBuffer().

Bool_t QwSIS3320_Channel::fHasSamplingData

private

Definition at line 149 of file QwSIS3320_Channel.h.

Referenced by InitializeChannel(), and ProcessEvBuffer().

std::vector<QwSIS3320_LogicalAccumulator> QwSIS3320_Channel::fLogicalAccumulators

private

Definition at line 182 of file QwSIS3320_Channel.h.

Referenced by AddLogicalAccumulator(), ClearEventData(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), operator+=(), operator-=(), operator=(), PrintInfo(), PrintValue(), ProcessEvent(), Ratio(), and Scale().

Int_t QwSIS3320_Channel::fNumberOfAccumulators

private

Definition at line 175 of file QwSIS3320_Channel.h.

Referenced by GetNumberOfAccumulators(), and SetNumberOfAccumulators().

UInt_t QwSIS3320_Channel::fNumberOfEvents

private

Current triggered event (allow for negative sentinel)

Definition at line 160 of file QwSIS3320_Channel.h.

Referenced by EncodeEventData(), GetNumberOfEvents(), and SetNumberOfEvents().

const unsigned int QwSIS3320_Channel::FORMAT_ACCUMULATOR = 0x0

staticprivate

Definition at line 191 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

const unsigned int QwSIS3320_Channel::FORMAT_LONG_WORD_SAMPLING = 0x1

staticprivate

Definition at line 192 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

const unsigned int QwSIS3320_Channel::FORMAT_SHORT_WORD_SAMPLING = 0x2

staticprivate

Definition at line 193 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

Double_t QwSIS3320_Channel::fPedestal

private

Definition at line 155 of file QwSIS3320_Channel.h.

Referenced by GetPedestal(), InitializeChannel(), ProcessEvent(), and SetPedestal().

UInt_t QwSIS3320_Channel::fSampleFormat

private

Number of triggered events.

Definition at line 163 of file QwSIS3320_Channel.h.

Referenced by EncodeEventData().

std::vector<QwSIS3320_Samples> QwSIS3320_Channel::fSamples

private

Definition at line 164 of file QwSIS3320_Channel.h.

Referenced by ClearEventData(), ConstructBranchAndVector(), FillHistograms(), GetSamples(), InitializeChannel(), IsGoodEvent(), Offset(), operator+=(), operator-=(), operator=(), PrintInfo(), ProcessEvent(), Scale(), and SetNumberOfEvents().

std::vector<QwSIS3320_Samples> QwSIS3320_Channel::fSamplesRaw

private

Definition at line 165 of file QwSIS3320_Channel.h.

Referenced by ClearEventData(), GetSamplesRaw(), InitializeChannel(), PrintInfo(), ProcessEvBuffer(), ProcessEvent(), and SetNumberOfEvents().

std::vector<Double_t> QwSIS3320_Channel::fSampleWindowAverages

private

Definition at line 171 of file QwSIS3320_Channel.h.

std::vector<std::pair<Double_t, Double_t> > QwSIS3320_Channel::fSampleWindows

private

Definition at line 172 of file QwSIS3320_Channel.h.

std::vector<Double_t> QwSIS3320_Channel::fTimeWindowAverages

private

Definition at line 169 of file QwSIS3320_Channel.h.

Referenced by ProcessEvent().

std::vector<std::pair<UInt_t, UInt_t> > QwSIS3320_Channel::fTimeWindows

private

Definition at line 170 of file QwSIS3320_Channel.h.

Referenced by ProcessEvent().

const Bool_t QwSIS3320_Channel::kDEBUG = kFALSE

staticprivate

Definition at line 145 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

const Double_t QwSIS3320_Channel::kNanoSecondsPerSample = 4.0

staticprivate

Conversion factor to translate the single sampling period to time. The ADC will sample at 250 MHz, corresponding with 4 ns per sample.

Definition at line 154 of file QwSIS3320_Channel.h.

const Double_t QwSIS3320_Channel::kVoltsPerBit = 5.0 / pow(2.0, 12)

staticprivate

Conversion factor to translate the average bit count in an ADC channel of the SIS3320 into average voltage. There are 2^12 possible states over the full 5 V range.

Definition at line 153 of file QwSIS3320_Channel.h.

const unsigned int QwSIS3320_Channel::MODE_ACCUM_EVENT = 0x1

staticprivate

Definition at line 185 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

const unsigned int QwSIS3320_Channel::MODE_MULTI_EVENT = 0x3

staticprivate

Definition at line 186 of file QwSIS3320_Channel.h.

Referenced by ProcessEvBuffer().

const unsigned int QwSIS3320_Channel::MODE_NOTREADY = 0xda00

staticprivate

Definition at line 188 of file QwSIS3320_Channel.h.